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

Question

如题所述

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));

图像为：

Answer 2

方程组求解要降阶才可以求解