微分方程的matlab求解..pptVIP

蓝色曲线 ——y（1）； （原方程解） 红色曲线 ——y（2）； 例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)曲线 范例 计算结果 范例 主 页 下一页 上一页 数学实验 定义：含有导数的方程称为微分方程。如 f(x, y(x), y’(x))=0 微分方程模型 1、微分方程的一般形式： F(x, y, y’,…,y(n) ) = 0 隐式 或