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 xdot=lorenz(t,x) xdot=[-8/3,0,x(2);0,-10,10;-x(2),28,-1]*x; 2、ltest.m x0=[0 0 0.1]; [t,x]=ode45(lorenz,[0,10],x0); plot(t,x(:,1),-,t,x(:,2),*,t,x(:,3),+) pause plot3(x(:,1),x(:,2),x(:,3)),grid on 计算结果如下图 范例 图中，x1的图形为实线(蓝），x2的图形为“*”线（绿）， x3的图形为“+”线（红）.取[t0，tf]=[0，10]。 若自变量区间取[0，20]、[0，40]，计算结果如下： 范例 曲线呈震荡发散状 三维图形的混沌状 ltest.m 观察结果： 1、该曲线包含两个“圆盘”，每一个都是由螺线形轨道构成。某些轨道几乎是垂直地离开圆盘中一个而进入另一个。 2、随着t的增加，x(t)先绕一个圆盘几圈，然后“跳”到另一个圆盘中。绕第二个圆盘几圈，又跳回原来的圆盘。 并以这样的方式继续下去，在每个圆盘上绕的圈数是随机的。 思考：该空间曲线与初始点x0的选择有关吗？ 1）x0=[0 0.1 0.1]‘;[t0,tf]=[0,30];解向量y 2）x00=[0.01 0.11 0.11]‘;[t0,tf]=[0,30];解向量x y – x = (y1-x1,y2-x2,y3-x3) 返 回 1、Apollo卫星的运动轨迹的绘制 实验内容 2、用向前欧拉公式和改进的欧拉公式求方程 y’=y-2x/y, y(0) = 1的数值解(0≤x≤1,h=0.1) 要求编写程序。 实验内容 3、Rossler微分方程组： 当固定参数b=2,c=4时，试讨论随参数a由小到大变化（如a∈(0,0.65])而方程解的变化情况，并且画出空间曲线图形，观察空间曲线是否形成混沌状？ 4、操练一。 返 回 实验内容 主 页 下一页 上一页 数学实验 定义：含有导数的方程称为微分方程。如 f(x, y(x), y’(x))=0 微分方程模型 1、微分方程的一般形式： F(x, y, y’,…,y(n) ) = 0 隐式 或 y(n) = f (x, y, y’,…,y (n-1) ) 显式 特殊情形： 2、一阶微分方程组的一般形式： 初始条件：y(x0) = y0 微分方程模型 ③　图形解 t y o ①简单的微分方程。 ②复杂、大型的微分方程。 返 回 ①　解析解 y = f(t) ②　数值解 (ti, yi) 欧拉方法 改进欧拉方法 梯形法 龙格-库塔法 微分方程求解方法简介 微分方程数值解 1、欧拉法 2、龙格—库塔法 数值求解思想：(变量离散化) 引入自变量点列{xn} → {yn}， 在x0?x1?x2?…?xn?…上求y(xn)的近似值yn.通常取等步长 h，即xn = x0+ n×h，或 xn = xn-1+ h，（n=1,2,…)。 1) 向前欧拉公式: （y’= f (x, y) ） y (xn+1) ? y(xn) + h f(xn, y(