matlab關(guān)于ode45解二階微分方程的困惑

matlab關(guān)于ode45解二階微分方程的困惑

一個(gè)二階微分方程：

y''+y'+y=sin(t)

初始條件為y(0)=5,y'(0)=6。

過程：

先降階為一階微分方程組

y'=z

z'=-z-y+sin(t)

編制如下函數(shù)m文件

function dy=weifen(t,x)

dy=zeros(2,1);

%y=x(1)

%z=x(2)

dy(1)=x(2);

dy(2)=sin(t)-x(2)-x(1);

*******************

然后用ode45解方程

[t,y]=ode45(@weifen,[0 20],[5 6])

plot(t,y)就可以畫出y和y'的圖像

但是我用同樣的思路，是不過在函數(shù)文件的編制時(shí)候?qū)⒑瘮?shù)的代表意義換一下，結(jié)果出來的圖形就完全不一樣了？

即我的m文件重新編制為

function dy=weifen(t,x)

dy=zeros(2,1);

%y=x(2)

%z=x(1)

dy(1)=x(1);

dy(2)=sin(t)-x(1)-x(2);

*******************

然后同樣用ode45解方程，根據(jù)函數(shù)文件只需將初始條件換一下位置

[t,y]=ode45(@weifen,[0 20],[6 5])

plot(t,y)就可以畫出y和y'的圖像，但是結(jié)果完全不一樣了，怎么回事？

難道必須有一定的順序么？

總結(jié)

以上是生活随笔為你收集整理的matlab ode45 二阶微分,matlab关于ode45解二阶微分方程的困惑的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。