第3章节差分微分matlab幻灯片.pptVIP

（三）可以用Matlab软件求常微分方程的数值解 [t，x]=solver（’f’,ts,x0,options） ode45 ode23 ode113ode15sode23s 由待解方程写成的m-文件名 ts=[t0，tf]，t0、tf为自变量的初值和终值 函数的初值 ode23：组合的2/3阶龙格-库塔-芬尔格算法 ode45：运用组合的4/5阶龙格-库塔-芬尔格算法 自变量值 函数值 用于设定误差限(缺省时设定相对误差10-3, 绝对误差10-6), 命令为：options=odeset（’reltol’,rt,’abstol’,at）, rt，at：分别为设定的相对误差和绝对误差. 不同求解器Solver的特点 求解器 Solver ODE类型 特点 说明 ode45 非 刚 性 一步算法；4，5阶Runge-Kutta方程；累计截断误差达(△x)3 大部分场合的首选算法 ode23 非 刚 性 一步算法；2，3阶Runge-Kutta方程；累计截断误差达(△x)3 使用于精度较低的情形 求解器Solver ODE类型 特点 说明 ode113 非刚性 多步法；Adams算法；高低精度均可到10-3～10-6 计算时间比ode45短 ode23t 适度刚性 采用梯形算法 适度刚性情形 ode15s 刚性 多步法；Gear’s反向数值微分；精度中等 若ode45失效时，可尝试使用 不同求解器Solver的特点 求解器Solver ODE类型 特点 说明 ode23s 刚性 一步法；2阶Rosebrock算法；低精度 当精度较低时，计算时间比ode15s短 ode23tb 刚性 梯形算法；低精度 当精度较低时，计算时间比ode15s短 参数设置 属性名 取值 含义 AbsTol 有效值：正实数或向量 缺省值：1e-6 绝对误差对应于解向量中的所有元素；向量则分别对应于解向量中的每一分量 RelTol 有效值：正实数 缺省值：1e-3 相对误差对应于解向量中的所有元素。在每步(第k步)计算过程中，误差估计为： e(k)=max(RelTol*abs(y(k)),AbsTol(k)) 参数设置 属性名 取值 含义 NormControl 有效值：on、off 缺省值：off 为‘on’时，控制解向量范数的相对误差，使每步计算中，满足：norm(e)=max(RelTol*norm(y),AbsTol) Events 有效值：on、off 为‘on’时，返回相应的事件记录 参数设置 属性名 取值 含义 OutputFcn 有效值：odeplot、odephas2、odephas3、odeprint 缺省值：odeplot 若无输出参量，则solver将执行下面操作之一： 画出解向量中各元素随时间的变化； 画出解向量中前两个分量构成的相平面图； 画出解向量中前三个分量构成的三维相空间图； 随计算过程，显示解向量 参数设置 属性名 取值 含义 OutputSel 有效值：正整数向量 缺省值：[] 若不使用缺省设置，则OutputFcn所表现的是那些正整数指定的解向量中的分量的曲线或数据。若为缺省值时，则缺省地按上面情形进行操作 Refine 有效值：正整数k1 缺省值：k = 1 若k1，则增加每个积分步中的数据点记录，使解曲线更加的光滑 参数设置 属性名 取值 含义 Jacobian 有效值：on、off 缺省值：off 若为‘on’时，返回相应的ode函数的Jacobi矩阵 Jpattern 有效值：on、off 缺省值：off 为‘on’时，返回相应的ode函数的稀疏Jacobi矩阵 参数设置 属性名 取值 含义 Mass 有效值：none、M、 M(t)、M(t,y) 缺省值：none M：不随时间变化的常数矩阵 M(t)：随时间变化的矩阵 M(t,y)：随时间、地点变化的矩阵 MaxStep 有效值：正实数 缺省值：tspans/10 最大积分步长 创建函数function2, 保存在function2.m中function f=function2(t,x) f=-x.^2; 在命令窗口中执行 [t,x]=ode45(function2,[0,1],1); plot(t,x,-,t,x,o); xlabel(time t0=0,tt=1); ylabel(x values x(0)=1); 创建函数function3,保存在function3.m中： function f=function3(t,x) f=[x(1)-0.1*x(1)*x(2)+0.01*t;... -x(2)+0.02*x(1)*x(2)+0.04*t]; 运行命令文件runf3.m [t,x]=ode45(function3,[0,20