matlab-常微分方程课件.pptx

Matlab—常微分方程求解;;一、常微分方程（组）符号解;;;;例1;例2;二、常微分方程（组）数值解;二、常微分方程（组）数值解;参数说明：;参数说明：;在区间tspan=[t0,tf]上，从t0到tf，用

初始条件y0求解显式微分方程y’=f(t,y)。

对于标量t与列向量y，函数f=odefun(t,y)

必须返回一f(t,y)列向量f。

解矩阵Y中每一行对应于返回时间列向量T中一个时间点。

要取得问题在其它指定时间点t0,t1,t2,…上解，则令tspan=[t0,t1,t2,…,tf]（要求是单调）。;用参数options（用命令odeset生成）设置属性（代替了缺省积分参数），再进行操作。惯用属性包含相对误差值RelTol（缺省值为1e-3）与绝对误差向量AbsTol（缺省值为每一元素为1e-6）。;将参数p1,p2,p3,..等传递给函数odefun，再进行计算。若没有参数设置，则令options=[]。;求解详细ODE基本过程：;求解详细ODE基本过程：;（3）依据（1）与（2）结果，编写能计算

导数M-函数文件odefile。

（4）将文件odefile与初始条件传递给求解

器Solver中一个，运行后就可得到ODE

、在指定时间区间上解列向量y（其中包含y及不一样阶导数）。;不一样求解器Solver特点;求解器Solver;不一样求解器Solver特点;参数设置;参数设置;参数设置;参数设置;参数设置;参数设置;例3;;例4;;;例5;创建函数function4，存在function4.m中

functionf=function4(t,x)

globalU

f=[x(2);U*(1-x(1)^2)*x(2)-x(1)];;