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...