matlab微分方程的求解的方法.pptVIP

例2 Van der pol 方程： 令 y1=x (t), y2 = x’(t); 该方程无解析解！ 范例 （1）编写M文件 ( 文件名为 vdpol.m): function yp = vdpol(t,y); yp=[y(2);(1-y(1)^2)*y(2)-y(1)]; （2）编写程序如下:（vdj.m） [t,y]=ode23(vdpol,[0,20],[3,0]); y1=y(:,1); % 原方程的解 y2=y(:,2); plot(t,y1,t,y2,‘--’) % y1(t),y2(t) 曲线图 pause, plot(y1,y2),grid, % 相轨迹图，即y2(y1)曲线 范例 蓝色曲线 ——y（1）； （原方程解） 红色曲线 ——y（2）； 计算结果 范例 范例 例3 考虑Lorenz模型： 其中参数β=8/3，σ=10，ρ=28 解：1）编写M函数文件(lorenz.m); 2) 数值求解并画三维空间的相平面轨线； (ltest.m) 范例 1、 lorenz.m function