常微分方程模型的解.pptVIP

第五章微分方程模型 鲁东大学数学与信息学院 魏建新 wjx0426@163.com 第五章 微分方程模型 5.1 小序,引例 5.2 微分方程的解 应用已知物理定律 §5.2 微分方程的解 5.2.1微分方程的解析解 5.2.2微分方程的数值解 5.2.1微分方程的解析解 在Matlab中,符号运算工具箱提供了功能强大的求解常微分方程的符号运算命令dsolve. 格式： r = dsolve(eq1,eq2,..., cond1,cond2,..., v) r = dsolve(eq1,eq2,...,cond1,cond2,...,v) r=dsolve(‘方程1’, …‘方程n’, ‘初始条件’, ‘自变量’) 记号: 在表达微分方程时，用字母D表示求微分，D2、D3等表示求高阶微分.任何D后所跟的字母为因变量，自变量可以指定或由系统规则选定为确省. 几个求解例子 练习 equ1=D2f+3*g=sin(x);equ2=Dg+Df=cos(x); [general_f,general_g]=dsolve(equ1,equ2,x) [f,g]=dsolve(equ1,equ2,Df(2)=0,f(3)=3,g(5)=1,x) 5.2.2微分方程的数值解 建立微分方程只是解决问题的第一步,通常需要求出方程的解来说明实际现象,并加以检验.如果能得到解析形式的解固然是便于分析和应用的.但是我们知道,只有线性常系数微分方程,并且自由项是某些特殊类型的函数时,才可以肯定得到这样的解,而绝大多数变系数方程、非线性方程都是所谓“解不出来”的,即使看起来非常简单的方程,于是对于用微分方程解决实际问题来说,数值解法就是一个十分重要的手段. 1.3.预估-校正Euler 法 综合使用前两种方法：用1去求初步的近似 值，用之再代入2 格式 注：1.当一个数值公式的截断误差可表示为O(hk+1)时(k为正整数，h为步长),称它是一个k阶公式. 2.欧拉法是一阶公式,改进的欧拉法是二阶公式. 2 龙格—库塔（Runge—Kutta）方法 2.2 龙格-库塔方法一般形式 2.3 经典4级4阶龙格-库塔方法 （三）一阶微分方程组和高阶微分方程初值问题的数值解 设有一阶微分方程组的初值问题 高阶微分方程的初值问题可以通过变量代换化为一阶微分方程组初值问题。 设有m阶常微分方程初值问题 (四) 刚性常微分方程（组） 例题1 求初值问题 的数值解. 例3 ：考虑初值问题 (2)把一阶方程组写成接受两个参数 和 返回一个列向量的M文件djl.m： function dy=djl(t,y); dy=[y(2);y(3);3*y(3)+y(2)*y(1)]; (或者dy=zeros(3,1);dy(1)= y(2);dy(2)= y(3); dy(3)=3*y(3)+y(2)*y(1)];) (3)取t0=0，tf=1，输入命令： [T,Y]=ode45(rigid,[0 1],[0 1 -1]); plot(T,Y(:,1),-,T,Y(:,2),*,T,Y(:,3),+) 例4：求van der Pol 方程 的数值解，其中 . （2）书写M文件（对于 ）vdp1.m: function dy=vdp1(t,y); dy=[y(2);(1-y(1)^2)*y(2)-y(1)]; （3）调用Matlab函数。对于初值 ， 解为 [T,Y]=ode45(vdp1,[0 20],[2;0]); plot(T,Y(:,1),-,T,Y(:,2),--) To Matlab(ff7) * * §5.1 小序引例 在许多实际问题中,当直接导出变量之间的函数关系较为困难,但导出包含未知函数的导数或微分的关系式较为容易时,可用建立微分方程模型的方法来研究该问题. 常用建模方法 已知定律 模拟近似法 微元法 例1一个较热的物体置于室温为180c的房间内,该物体最初的温度是600c,3分钟以后降到500c .想知道它的温度降到300c 需要多少间?10分钟以后它的温度是多少? 牛顿冷却(加热)定律:将温度为T的物体放入处于常温m的介质中时,T的变化速率正比于T与周围介质的温度差. 分析:假设房间足够大,放入温度较低或较高的物体时,室内温度基本不受影响,即室温分布均衡,保持