MATLAB解常微分方程与其应用.pptVIP

格式一： y=dsolve(f1, f2, …, fm) 格式二： y=dsolve(f1, f2, …, fm ,x) %指明自变量 fi既可以描述微分方程，又可描述初始条件或边界条件。如： 描述微分方程时 描述条件时 例5 求方程 的解析解. 解： syms x,t; x=dsolve(D1x=x*(1-x^2),t) x = 1 0 -1 (-1/(exp(C3 - 2*t) - 1))^(1/2) 例题： 4． 几种方法的比较 1 欧拉方法 对于一般的常微分方程： 取前一点的斜率作为平均斜率，即： 整体截断误差： O(h) 2 向后的欧拉方法 取后一点的斜率作为平均斜率，即： 整体截断误差： O(h) 3 改进的欧拉方法 取两点斜率平均值，即： 整体截断误差： O(h2) 欧拉方法（前、后）和改进的欧拉方法优点是算法简单，但计算精度较低，不能满足实际问题求解对精度的要求，所以使用较少。 从以上三种方法的比较，可以推测：若取多点处斜率的加权平均作为平均斜率，精度会更高，这就是龙格—库塔法的基本思想。 四阶龙格—库塔法计算公式为： 整体截断误差： O(h4) 例题： time_Euler = 0.2500 time_EulerPro = 0.5000 time_RK4 = 1.0150 5． MATLAB的常微分方程函数 ode45 ode23 ode113 ode15s ode23s ode23t ode23tb 格式 [x,y]=ode45(′ fun′, [x0,xn], y0,option] 说明： 适用于求解一阶常微分方程组 fun定义微分方程组的函数文件名 [x0,xn]求解区域 y0初始条件向量 option可选参数，由ODESET函数设置，比较复杂 x输出自变量向量，y输出[y, y ′, y ″,..] 没有一种算法可以有效地解决所有的 ODE 问题，因此MATLAB 提供了多种ODE函数。 函数 ODE类型 特点 说明 ode45 非刚性 单步法；4，5 阶 R-K 方法；累计截断误差为 (△x)3 大部分场合的首选方法 ode23 非刚性 单步法；2，3 阶 R-K 方法；累计截断误差为 (△x)3 使用于精度较低的情形 ode113 非刚性 多步法；Adams算法；高低精度均可到 10-3～10-6 计算时间比 ode45 短 ode23t 适度刚性 采用梯形算法 适度刚性情形 ode15s 刚性 多步法；Gear’s 反向数值微分；精度中等 若 ode45 失效时，可尝试使用 ode23s 刚性 单步法；2 阶Rosebrock 算法；低精度 当精度较低时，计算时间比 ode15s 短 ode23tb 刚性 梯形算法；低精度 当精度较低时，计算时间比ode15s短 MATLAB解常微分方程及其应用 授课教师：唐跃龙 Telephone QQ: 794218136 Email:tangyuelonga@163.com 1． 微分方程的基本概念 定义1 含有未知函数的导数(或微分)，同时也可能含有未知函数与自变量的方程，称为微分方程，未知函数是一元函数的，称为常微分方程；未知函数是多元函数的，称为偏微分方程. 定义2 微分方程中所出现的未知函数的最高阶导(或微分)的阶数，叫做微分方程的阶. 定义3 若 是 的一次有理整式， 称为 阶线性常微分方程. 定义4 不是线性的微分方程称为非线性常微分方程. 定义8 若某函数代入微分方程能使该方程成为恒等 式，该函数就叫做该微分方程的解. 定义5 不含有未知函数及其导数(或微分)的项称为 微分方程自由项. 定义6 自由项为零的微分方程称为齐次方程. 定义7 自由项不为零的微分方程称为非齐次方程. 例1. 指出下列微分方程的阶数并判断是线性还是非线性、 齐次还是非齐次方程： 定义5 确定了通解中任意常数以后得到的解称为微分方程 的特解. 定义4 若微分方程的解中所含相互独立的任意常数的个数 与微分方程的阶数相同，这样的解叫做微分方程的通解. 例2. 验证函数 是微分方程 的解. 例3. 求方程 的通解.