matlab解常微分方程新的特别方法.pptVIP

应用matlab软件对常微分方程求解;前沿：;一、常微分方程（组）的符号解;(3)在微分方程（组）的表达式eq中，大写字母D表示对自变量(设为x)的微分算子：D=d/dx，D2=d2/dx2，…。微分算子D后面的字母则表示为因变量，即待求解的未知函数。初始和边界条件由字符串表示：y(a)=b，Dy(c)=d，D2y(e)=f，等等，分别表示;(4)若边界条件少于方程（组）的阶数，则返回的结果r中会出现任意常数C1，C2，…；;(6)若没有给定输出参量，则在命令窗口显示解列表。若该命令找不到解析解，则返回一警告信息，同时返回一空的sym对象。这时，用户可以用命令ode23或ode45求解方程组的数值解。;例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中 function f = function4(t,x) global U f = [x(2);U*(1-x(1)^2)*x(2)-x(1)];;