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

2021-08-14 05:14:31 字數 5091 閱讀 5842

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')即可。

初一數學 用代入法解下列方程組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...

用代入法解下列方程組 4 x 2 5y 1,2x 3 y

解答過程為 4 x 2 5y 1,2x 3 y 2 3方程組化簡可得 4x 5y 7 1 2x 3y 3 2 2 x2後得 4x 6y 6,所以4x 6 6y 3 將 3 帶入 1 得 6 6y 5y 7 解得 y 1 4 將 4 帶入 3 得 4x 12,所以x 3。所以方程組的解為 x 3,y ...

用克萊姆法則解下列線性方程組,克萊姆法則解線性方程組

d 2,2,1,1 4,3,1,2 8,3,3,4 3,3,2,2 28 d1 4,2,1,1 16,3,1,2 12,3 3,4 6,3,2,2 148 d2 2,4,1,1 4,16,1,2 8,12,3,4 3,6,2,2 84 d3 2,2,4,1 4,3,16,1 8,3,12,4 3,3...