数学物理方法(数值方法)怎么用?

一、理论解,用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));

图像为: