matlab关于ode45解二阶微分方程的困惑
matlab关于ode45解二阶微分方程的困惑
一个二阶微分方程:
y''+y'+y=sin(t)
初始条件为y(0)=5,y'(0)=6。
过程:
先降阶为一阶微分方程组
y'=z
z'=-z-y+sin(t)
编制如下函数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'的图像
但是我用同样的思路,是不过在函数文件的编制时候将函数的代表意义换一下,结果出来的图形就完全不一样了?
即我的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解方程,根据函数文件只需将初始条件换一下位置
[t,y]=ode45(@weifen,[0 20],[6 5])
plot(t,y)就可以画出y和y'的图像,但是结果完全不一样了,怎么回事?
难道必须有一定的顺序么?