用matlab解下列方程組(含積分)

時間 2021-08-14 05:14:58

1樓:匿名使用者

從原理上說,tx、ty應該可以是隨時間變化的函式,但那樣就成了一個泛函問題。如果tx、ty取常數值方程有解,而且沒有其它要求,當然應該使用常數解。

現在假設tx、ty為常數,可以使用fsolve求解:

m = 2400;

q = 7500/2490;

a = 1.63;

% tx = x1, ty = x2, t = x3

f = @(x) [quadl(@(t)x(1)./(m-q*t),0,x(3))-1692;

quadl(@(t)x(2)./(m-q*t)-a,0,x(3))-57;

sqrt(x(1)^2+x(2)^2)-7500];

x = fsolve(f,[5304,5304,700]);

fprintf('tx = %.2f\nty = %.2f\nt = %.2f\n',x)

得到的結果:

tx = 6876.27

ty = 2994.48

t = 417.07

2樓:匿名使用者

你給的問題還不是很清楚

初始條件,結束條件 也沒有給清楚

你的兩個微分方程應該是關於兩個速度 vx vy的,他們的導數就是加速度,其實就是兩個方向的牛頓方程dvx f vx---- = -------- ------ (1)

dt m-q*t v

dvy f vy---- = -------- ------ - a (2)

dt m-q*t v

其中 v = sqrt(vx^2+vy^2) 是速度總值然後就是要給出 邊界上的條件和約束

例如vx(0) vy(0) 的值(初始兩個方向的速度)最後時刻應該滿足的條件

這些條件,引數你都要給齊 ,而且注意引數的單位,最後都用標準國際單位這樣才能夠正確的得到解

matlab求解含積分的方程組

3樓:匿名使用者

用matlab求解含積分的方程組,可以這樣實現:

1、先用integral函式,求解定積分,即

integral(@(t)sqrt(a^2.*sin(t).^2+b^2.*cos(t).^2),0,pi/2)

2、再用vpasolve函式,求解β、θ值,即

[beta1,theta1]=vpasolve(eq1,eq2, [beta1,theta1])

這裡:eq1是指方程一,eq2是指方程二

3、詳細的執行**為

a=15;b=14;r=14.5;

syms theta1 beta1

eq1=integral(@(t)sqrt(a^2.*sin(t).^2+b^2.*cos(t).^2),0,pi/2)+beta1*pi*r/180;

eq2=(a*sin(theta1))/(b*cos(theta1));

[beta1,theta1]=vpasolve(eq1==100,eq2==tan(beta1), [beta1,theta1])

這裡:beta1是指β,theta1是指θ

執行結果

用matlab求一個含有積分的方程組 5

4樓:匿名使用者

用fsolve求一個含有積分的方程組,可以這樣來寫:

第一步:建立自定義函式 func()

第二步:用fsolve求解其近似值

[k,fval] =fsolve(@(k)func(k),k0)  %k為係數a、b,即a=k(1),b=k(2),fval是每個等式的值

第三步:程式設計執行

執行結果:a=9.01704582797954;b =25.897553850857

執行**及執行過程。

5樓:雲南美食吧

(1830658703514901060,替換網頁地址question後面的那串數字即可)。

這個問題應該很難求出解析解,應該使用fsolve求數值解。

參考**(題中所給條件u=1/6疑為u=g/6之誤):

m = 2400;q = 7500/2490;g = 9.8;a = g/6;% tx = x1, ty = x2, t = x3f = @(x) [quadl(@(t)x(1)./(m-q*t),0,x(3))-1692; quadl(@(t)x(2).

/(m-q*t)-a,0,x(3))-57; sqrt(x(1)^2+x(2)^2)-7500];x = fsolve(f,[5304,5304,700]);fprintf('tx = %.2f\nty = %.2f\nt = %.

2f\n',x)

求出的結果:

tx = 6873.97ty = 2999.75t = 417.17

6樓:獅子中國博物館

7(2x-1)-3(4x-1)=4(3x 2)-1; (5y 1) (1-y)= (9y 1) (1-3y); 20% (1-20%)(320-x)=320×40% 2(x-2) 2=x 1 2(x-2)-3(4x-1)=9(1-x) x/3 -5 = (5-x)/2 2(x 1) /3=5(x 1) /6 -1 (1/5)x 1 =(2x 1)/4 (5-2)/2 - (4 x)/3 =1 x/3 -1 = (1-x)/2 (x-2)/2 - (3x-2)/4 =-1 11x 64-2x=100-9x 15-(8-5x)=7x (4-3x) 3(x-7)-2[9-4(2-x)]=22 3/2[2/3(1/4x-1)-2]-x=2 2(x-2)-3(4x-1)=9(1-x) 11x 64-2x=100-9x 15-(8-5x)=7x (4-3x) 3(x-7)-2[9-4(2-x)]=22 3/2[2/3(1/4x-1)-2]-x=2 2(x-2) 2=x 1 7(2x-1)-3(4x-1)=4(3x 2)-1(5y 1) (1-y)= (9y 1) (1-3y)[ (- 2)-4 ]=x 220% (1-20%)(320-x)=320×40%2(x-2) 2=x 1 6。

2(x-2)-3(4x-1)=9(1-x) 7。11x 64-2x=100-9x 15-(8-5x)=7x (4-3x) 3(x-7)-2[9-4(2-x)]=22 3/2[2/3(1/4x-1)-2]-x=25x 1-2x=3x-23y-4=2y 187x*13=57z/93=41 15x 863-65x=54 58y*55=274892(x 2) 4=92(x 4)=103(x-5)=184x 8=2(x-1)3(x 3)=9 x6(x/2 1)=129(x 6)=632 x=2(x-1/2)8x 3(1-x)=-27 x-2(x-1)=1x/3 -5 = (5-x)/2 2(x 1) /3=5(x 1) /6 -1 (1/5)x 1 =(2x 1)/4 (5-2)/2 - (4 x)/3 =1 15x-8(5x 1。

怎麼用matlab解含有字母系數的方程組的解,舉個簡單例子

7樓:大野瘦子

用法以這個為例:

x+a*y=10

x-b*y=1

其中x,y為變數,a,b為字母系數.

只要在matlab中輸入

syms x,y,a,b

[x y]=solve('x+a*y=10','x-b*y=1','x','y')

即可求出解

x =(a + 10*b)/(a + b)

y =9/(a + b)

對於函式solve的具體用法,可以通過輸入help solve來學習。

matlab中方程求解的基本命令

1.roots(p)  %求多項式的根,其中p是多項式向量。

例求x3-x2+x-1=0的根

解:>>roots([1,-1,1,-1])

注: [1,-1,1,-1]在matlab中表示多項式 x3-x2+x-1

2.solve(fun)     %求方程fun=0的符號解,如果不能求得精確的符號解,可以計算可變精度的數值解

例:用solve求方程x9+x8+1=0的根

解:>>solve(‘x^9+x^8+1’)

給出了方程的數值解(32位有效數字的符號量)

3.solve(fun,var)     %對指定變數var求代數方程fun=0的符號解。

例:解方程 ax2+bx2+c=0

解:>>syms a b c x;

>>f=a*x^2+b*x+c;

>>solve(f)

如果不指明變數,系統預設為x,也可指定自變數,比如指定b為自變數

>>symsa b c x;

>> f=a*x^2+b*x+c;

>>solve(f,b)

4.fsolve(fun,x0)    %求非線性方程fun=0在估計值x0附近的近似解。

例:用fsolve求方程x=e-x在0附近的根

解:>>fsolve(‘x-exp(-x)’,0)

5.fzero(fun,x0)   %求函式fun在x0附近的零點

例:求方程x-10x+2=0在x0=0.5附近的根

解:>>fzero(‘x-10^x+2’,0.5)

8樓:匿名使用者

舉個簡單例子,解方程組

x+a*y=10

x-b*y=1

其中x,y為變數,a,b為字母系數。

只要在matlab中輸入

syms x,y,a,b

[x y]=solve('x+a*y=10','x-b*y=1','x','y')

即可求出解

x =(a + 10*b)/(a + b)y =9/(a + b)

對於函式solve的具體用法,可以通過輸入help solve來學習。

希望我的回答能夠解決您的疑問,謝謝。

9樓:匿名使用者

例:解方程組a*x+2*y=4,4*x+b*y=7matlab**:syms a b x y;

z=solve('a*x+2*y=4','4*x+b*y=7')即可。

解下列方程組 (1m,解下列方程組 (1) m

3 3 x 1 y 5 5 y 1 3 x 3 解 將原方程組整理得 3x y 8 3x 5y 14 得 4y 22,y 11 2 把y 11 2代入到 x 9 2 原方程組的解是 x 9 2 y 11 2 4 x 7 9 2 x 5 3y 3 解 由x 7 9得 x 2 把x 2代入到2 x 5 ...

解下列方程組m 2 2 2(m n 5m n)

m 96 17,n 45 17 m 5 n 2 2 兩邊同時乘以10得2m 5n 20 2 m n 5 m n 23 2m 2n 10 m n 23 3m n 10 23 3m n 13 5得 17m 45 m 45 17 把m 45 17代入 得 90 17 5n 20 5n 90 17 20 5...

初一數學 用代入法解下列方程組x y z 2,3x y 4z 5,2x 3y 2z

1.z 38 2x 3y 3x 4y 2 38 2x 3y 56 x 2y 20 x 2y 20 4x 5y 38 2x 3y 66 2x 2y 28 x y 14 y 14 x x 2 14 x 20 x 8 x 8y 6 z 42.z 9 2x 3y 3x 4 9 2x 3y 7 3x 36 8...