MATLAB教程〉木6】微分方程.pptVIP

1.4.8 解常微分方程 (ordinary differential equation, ODE) ;Runge Kutta法 ;二、解题方法 1.编写odefile 文件格式：function ydot=odefile(t,y) ydot=[待求的函数];(1)建立ode函数文件 % m-function, g1.m function dy=g1(x,y) dy=3*x.^2;;例3 ：解二阶微分方程;为方便令x1=x，x2=x’分别对x1,x2求一 阶导数，整理后写成一阶微分方程组 形式 x1’=x2 x2’=x2(1-x12)-x1;例：研究有空气阻力时抛体运动的特征。比较下面三种情况下的抛体的轨道：没有空气阻力；空气阻力与速度一次方成正比；以及空气阻力与速度二次方成正比。;%%program ddqxn.m m=1; b=[0,0.1,0.1]; p=[0,0,1]; %设置参数 figure axis([0 9 0 4]); %设置坐标轴的范围 hold on for i=1:3 %解微分方程三次 [t,y]=ode45(ddqxnfun,[0:0.001:2],[0,5,0,8],[ ],b(i),p(i),m); comet(y(:,1),y(:,3)) %画轨道的彗星轨迹 end;3.确定求解的条件和要求 [T,Y]=ode23(‘F’,tspan,y0,options,p1,p2,…) Odeset 建立和修改优化选项 语句格式： options=odeset(‘name1’,value1, ‘name2’,value2,…) 指定各个参数的取值，对不指定取值的参数，取默认值 options=odeset(odeopts,‘name1’,value1,…) 使用原来的优化选项，但对其中部分参数指定新值 options=odeset(oldopts,newopts) 合并了两个优化选项oldopts和newopts，重复部分取newopts的指定值; odeset AbsTol: [ positive scalar or vector {1e-6} ] 绝对误差。默认1e－6 BDF: [ on | {off} ] 在使用ODE15s时是否使用反向差分公式 RelTol: [ positive scalar {1e-3} ] 相对误差。默认1e－3 Events: [ function ] 设置事件 InitialStep: [ positive scalar ] 解方程时使用的初始步长 Jacobian: [ matrix | function ] 设定是否计算雅各比矩阵 JPattern: [ sparse matrix ] 若返回稀疏雅各比矩阵，为‘on’ Mass: [ matrix | function ] 返回质量矩阵 MassSingular: [ yes | no | {maybe} ] 质量矩阵是否为奇异矩阵 MaxOrder: [ 1 | 2 | 3 | 4 | {5} ] ODE15s解刚性方程的最高阶数 MaxStep: [ positive scalar ] 步长上界 ;NormControl: [ on | {off} ] 解向量范数的误差控制 OutputSel: [ vector of integers ] 选择输出解向量哪个分量 Refine: [ positive integer ] 细化输出因子 Stats: [ on | {off} ] 显示计算成本统计结果 Vectorized: [ on | {off} ] 设定向量形式的ODE函数文件 OutputFcn: [ function ]; 气象学家Lorenz提出一篇论文，名叫「一只蝴蝶拍一下翅膀会不会在Taxas州引起龙卷风？」论述某系统如果初期条件差一点点，结果会很不稳定，他把这种现象戏称做「蝴蝶效应」。Lorenz为何要写这篇论文呢？ 这故事发生在1961年的某个冬天，他如往常一般在办公室操作气象电脑。平时，他只需要将温度、湿度、压力等气象数据输入，电脑就会依据三个内建的微分方程式，计算出下一刻可能的气象数据，因此模拟出气象变化图。 这一天，Lorenz想更进一步了解某段纪录的後续变化，他把某时刻的气象数据重新输入电脑，让电脑计算出更多的後续结果。当时，电脑处理数据资料的数度不快，在结果出来之前，足够