1樓:
dy=zeros(1,3); 這說明返回值是個行向量,但是ode要求返回值是列向量。
改成 dy=zeros(3, 1);就對了。
matlab裡ode45解決二階常係數微分方程用法
求助,用 matlab 的 ode45 求解微分方程組
2樓:景志國敖央
例如:求解下列微分方程組。
求解步驟:①自定義函式。
functiondy=
rigid(t,y)dy=
zeros(3,1); a
column
vector
dy(1)=y(2)
*y(3);
dy(2)=-y(1)
*y(3);dy(3)
*y(1)*y(2);
②在執行視窗下,執行下列命令。
options
=odeset('reltol',1e-4,'abstol',[1e-4
1e-41e-5]);
[t,y]=ode45(@rigid,[0
11],options);
plot(t,y(:,1),'t,y(:,2),'t,y(:,3),'
③求解結果。
matlab中ode45方法求解微分方程組 20
3樓:匿名使用者
使用ode45解微分方程組,可以用下列格式來求解。
[t,x]=ode45(odefun,tspan,x0)式中:odefun——給定微分方程組的自定義函式;tspan——自變數t的範圍,即【t0,tf】;x0——x的初值向量,即x1(0),x1'(0),x2(0),x2'(0),x3(0),x3'(0)
由於沒有給出自變數t的範圍和x的初值,所以不能得到其x1,x2,x3的值。
4樓:來自天湖山嫋嫋娜娜的黑耳鳶
我想問一下**中的方程組,右邊若採用三角函式該怎麼程式設計序求解,能不能指導一下,謝謝。
求助,用 matlab 的 ode45 求解微分方程組
5樓:匿名使用者
例如:求解下列微分方程組。
求解步驟:①自定義函式
function dy = rigid(t,y)
dy = zeros(3,1); a column vector
dy(1) =y(2) *y(3);
dy(2) =y(1) *y(3);
dy(3) =y(1) *y(2);
②在執行視窗下,執行下列命令。
options = odeset('reltol',1e-4,'abstol',[1e-4 1e-4 1e-5]);
[t,y] =ode45(@rigid,[0 12],[0 1 1],options);
plot(t,y(:,1),'t,y(:,2),'t,y(:,3),'
③求解結果。
用matlab解二元二階微分方程(用ode45)
6樓:二中混過
function dy=dong22(t,y)m=1;
k=2;g= %注意輸入g l0
l0=dy=[y(2);
y(1)*y(4)^2+g*cos(y(3))-k/m)*(y(1)-l0);
y(4);-2/y(1))*y(2)*y(4)-(g/y(1))*sin(y(3))]
然後和一次。
微分方程的內一樣容。
matlab問題:使用ode45求如下微分方程的數值解,繪製y-t圖
7樓:匿名使用者
用ode45()求微分方程的格式為。
[t,y] =ode45(odefun,tspan,y0)odefun——微分方程自定義函式;tspan——t的區間;y0——y,dy的初始值。
用plot()繪出微分方程的解,即y(t)函式圖形。
對於本提問,可以按下列**實施。
下圖是精確解與數值解的圖形比較。
關於全微分方程,關於全微分方程的解
檀君博 不可能對,您的理解有問題,沒明白全微分方程的實質。全微分方程實際上是方程可以寫成d f x,y 0的形式,然後對兩邊同時取積分,解得f x,y c為原方程的解,例如2xdx 3y 2 方程可以化為d x 2 d y 3 0等價於d x 2 y 3 0直接積分得x 2 y 3 c,因此原方程也...
簡諧運動微分方程怎麼解,簡諧運動的微分方程是怎麼得到的?
本題的解題過程如下 1 無阻尼的簡諧自由運動的微分方程 mx kx 0 1 2 初始條件 x 0 x0 x 0 x 0 2 1 的特徵方程 ms 2 k 0 3 解出 s1 k m 0.5 s2 k m 0.5 4 3 1 的通x t c1e s1t c2e s2t 5 根據 2 c1 c2 x0 ...
什麼是解微分方程,什麼叫微分方程?如何理解?包含哪些形式?
微分方程指描述未知函式的導數與自變數之間的關係的方程。解微分方程就是解答微分方程的函式值,微分方程的解是一個符合方程的函式。而在初等數學的代數方程,其解是常數值。介紹含有未知函式的導數,符合定義式,一般的凡是表示未知函式 未知函式的導數與自變數之間的關係的方程,叫做微分方程。未知函式是一元函式的,叫...