Question

如何用matlab求解一个二阶常系数微分方程组

Answer 1

一、理论解，用dsolve函数（在command window 中输入doc dsolve可以查看帮助）
示例：

代码：
[x,y]=dsolve('D2x+2*Dx=x+2*y-exp(-t)','Dy=4*x+3*y+4*exp(-t)')

二、数值解，用ode45，或ode23, ode15s其他函数。（在command window 中输入doc ode45可以查看帮助）
示例：

1. 建立M文件
function dy=vdp(t,y)
dy=[y(2);
1000*(1-y(1)^2)*y(2)-y(1);];

2. 调用MATLAB 函数ODE15S
[T,Y]=ode15s('vdp',[0 1000],[2 0]);
plot(T,Y(:,1));
图像为：