Question

二阶微分方程组的MATLAB解法

Answer 1

一个二阶微分方程：

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'的图像