matlab画斜率场,微分方程的matlab求解教程.ppt

使用Matlab软件求数值解时，高阶微分方程必须等价地变换成一阶微分方程组. 注意: 选择一组状态变量 注意 1、建立M文件函数 function xdot = fun(t,x,y) xdot = [x2(t)；x3(t)；…；f(t, x1(t), x2(t),…xn(t))]； 2、数值计算(执行以下命令) [t,x1,x2,…,xn]=ode45(‘fun',[t0,tf], [x1(0),x2(0),…,xn(0)]) 例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) 欧拉方法 改进欧拉方法