1樓:花開勿敗的雨季
方程或者初始條件有問題,請認真核實。
原因:如果按照現在的初始條件,則第一個微分方程各項為0,方程不起作用。而且就演算法而言,二階(最高階)導數的係數為0,會導致數值積分演算法無法啟動。
參考**大致如下:12
3456
78910
1112
1314
1516
17m=1;
g=9.8;
k=100;
l0=1.1;
x0=[0.1 0 0 0];
% 定義x1=l, x2=l', x3=a, x4=a'
dx=#(t,x)[x(2); (m*x(1)*x(4)^2-m*g*cos(x(3))-k*(x(1)-l0))/m; ...
x(4); (-2*m*x(1)*x(2)*x(4)+m*g*x(1)*sin(x(3)))/(m*x(1)^2)];
[t,x]=ode15s(dx,[0 10],x0);
lstr = ''', '\alpha', '\alpha''' };
for i=1:length(lstr)
subplot(2,2,i)
plot(t, x(:,i));
xlabel('time')
ylabel( lstr )end
2樓:匿名使用者
function dz=yueh(t,z) %儲存為 yueh.m 檔案
dz=[z(2);
10-z(1)-0.001*z(3);
z(4);
-0.0001*z(1)-0.001*z(3)];
執行一下:
[t,z]=ode45('yueh',[0 1],[0 0 0 0]);
% 感覺初始值不完整啊,這裡假設x(0)=θ(0)=x'(0)=θ'(0)=0
plot(z(:,1),z(:,3)) % 這裡是 x-θ 影象
3樓:光之痕跡
微分方程初值問題,用ode45可以解
龍格庫塔求解二階微分方程組的matlab程式設計
4樓:匿名使用者
matlab求解x''+0.7x'+0.8x'|x'|+25.6x-25.6x³=0二階微分方程組的方法,可以按下列步驟進行:
1、建立自定義函式func()
function f = func(t,x)
%x''+0.7x'+0.8x'|x'|+25.6x-25.6x³=0
f(1)=x(2);
f(2)=25.6*x(1)^3-25.6*x(1)-0.8*x(2)*abs(x(2))-0.7*x(2);
f=f(:);
2、建立龍格庫塔演算法函式runge_kutta()
呼叫格式:[t,x] = runge_kutta(@(t,x)func(t,x),x0,h,a,b);
3、然後根據x和x'資料,繪製出x(t)、x′(t)的圖形。
plot(x(:,1),x(:,2))
怎麼用matlab求解二階微分方程組
如何用matlab求解二階微分方程,以及程式例項
5樓:用著追她
1、首先看一下matlab求解方程
的方法,首先指明所解方程的變數,然後指明方專程,未知數和限制條件,屬最後求解方程。
2、在matlab命令列視窗中輸入syms x [x,params,conds]=solve(sin(x)==1,'returnconditions', true) ,按回車鍵可以得到方程解。
3、轉換一下,可以看到sin(x)=1方程的解是如下圖所示 。
4、也可以求解下面的一個方程。
5、輸入syms a b c y x;[x,y]=solve([a*x^2+b*y+c==0,a*x+2*y==4],[x,y])。
6、按回車鍵可以得到方程解。
6樓:匿名使用者
^1、對
copy於求數值解的微分方程,你可以用ode45()函式求解。如求下列微分方程
func。m %自定義微分方程的函式
function z = func(t,y)z =[y(2);(1-y(1)^2)*y(2)-y(1)];
main。m %主程式
clear all;close all;clcy0 = [0.25;0];
h = 0.1;
a = 0;
b = 20;
[t1 y1] = ode45(@func,y0,h,a,b)2、對於求解析解的微分方程,你可以用dsolve()函式求解。如求微分方程x*y''+x﹡(y')^2-y'=0的解析解,可以下列步驟計算
>> syms y(x)
>>dy = diff(y);d2y = diff(y, 2);
>>dsolve(x*d2y+x*(dy)^2-dy==0,'x')
7樓:閻輝門曼雲
^dsolve('d2y=x','x')
ans=
x^3/6
+c2*x+c3
d2y的意思就是復y的二階微分項制
不明白你的問題什麼意思,要輸入的話直接定義符號變數輸入syms
d2xx
d2x=x;
如果是矩陣,那就用矩陣表示
求二階微分方程的通解,高等數學,二階微分方程,求通解,需要詳細步驟,謝謝
2y y y 3e x,先求齊次方程通解。令2t 2 t 1 0,解得t 1或1 2即齊次解為y a e x b e 1 2x 其中a,b r 再求1個特解即可。令y c e x,則2c c c 3,即c 3 2故問題的解為3 2 e x a e x b e x 2 其中a,b r 北極灬寒冰 可以...
關於微分方程已知某二階線性微分方程三解為X,X 1,eX X求其通解及原方程
通解寫成這樣的形式是線性方程解的結構決定的。高等數學 同濟大學第五版p296 若原方程為y p x y q x y f x 則對應的齊次方程y p x y q x y 0 y1 x,y2 x 1,y3 ex x為原方程的特解,則可知y1 y3,y2 y3是對應齊次方程的兩個線性無關的特解,於是 的通...
可降階的二階微分方程和二階常係數線性微分方程的區別
倔強的水蘿蔔 可降階的就是把y 換成y,算出y後再積分!實際上就是一階的! 不要說話 可降階的二階微分方程 1,y f x 型的微分方程 此類方程特點是 方程右端僅含有自變數x,只需積分兩次便可得到方程的通解。2,y f x,y 型的微分方程 此類方程特點是 方程右端不顯含未知函式y。作變數代換y ...