常微分方程数值解PPT.pptVIP

第四次：常微分方程数值解一： 引言： 1：微分方程在数模中有重要作用。 2：列出微分方程仅是第一步，求解微方程为第二步。 3：但仅有少数微分方程可解析解，大部分非线性方程，  变系数方程，均所谓“解不出来” x（t）甲 y（t）乙 其中 r为甲独立生存的增长率：a反映捕食者对食饵的捕食能力 。 d为乙无甲的死亡率； b反映食饵对捕食者的供养能力。 初值为 x（0）= x0 y（0）= y0 ……（2）2；试用数值解讨论以下问题：[（1）无解析解] ? 设r =1 ，d=0.5, a=0.1 b=0.02, x0 =25, y0=2 求模型 (1)在(2)下的数值解,画出函数 x(t),y(t)图形以及相图 (x,y),观察x(t),y(t)的周期变化,近似地确定争的周期和x,y的最大、小值,近似计算x,y在一周期内的平均值. ? 与(掠俘问题讨论过的理论值)比较. 3:导弹跟踪问题： 某军的一导弹基地（位于坐标原点（0，0））发现基地正北方向120km处海面（位于坐标原点（0，120））上有敌舰一艘正以90km/h的速度向正东方向行驶，该基地立即发射导弹跟踪追击敌舰,导弹速度为450km/h，自动导航系统使导弹在 任意时刻都能对准敌舰，试问导弹在何时何处击中敌舰？ （1）假设：在t时刻导弹位于P(x(t),y(t)),敌舰位于M(90t,H)(其中H=120) （2）：建模： 三:欧拉方法和龙格--库塔方法. 常微分方程初值问题的提法是:设有一阶方程  和初始条件: C:误差估计：假设 yn没有误差，则由（不考虑累积误差） B:(6)式的计算：因为（6）式为隐式，无法用它直接计算yn+1:(6)式通常用迭代法计算即先由向前欧拉公式（5）产生初值: C:误差估计：局部截断误差精度为1阶。D:几何意义： （3）：梯形公式： A: 方法： B:(7)式也需要象向后欧拉公式一样进行迭代求yn+1； C:误差估计：局部截断误差精度为2阶。 D:缺点：需要迭代，计算量大。 （4）：改进的欧拉公式： A:引言：向前欧拉公式（5）计算简单，但精度只有1阶；梯形公式精度提高，但迭代太繁；结合两者得改进的欧拉公式。 B:先由向前欧拉公式（5）计算yn+1的预测值 ； 再把它代入梯形公式（7）右端，作为校正，即： C:对于这四种欧拉方法：通常用向前欧拉公式（5）和改进的欧拉公式（8）。 （5）：实例： （6）：用欧拉方法解微分方程组（和高阶微分方程） A:例掠俘问题的方程组； B:例： C:解高阶微分方程，只需先把它化为方程组 2:龙格--库塔方法: (1):思想:欧拉方法的基本思想?用差商代替导数; 由微分中值定理: (2):2阶龙格--库塔公式: 可以证明2阶龙格--库塔公式就是改进的欧拉公式. (3):4阶龙格--库塔公式为: (4):龙格--库塔公式为也可以推广至解微分方程组和高阶方程. (5):龙格--库塔公式的MATLAB实现(MATLAB5.0) 2阶,3阶用[t,x]=ode23(‘f’,ts,x0,options); 4阶,5阶用[t,x]=ode45(‘f’,ts,x0,options); 参数意义:其中f 是由待解方程写成的m文件名; 其中ts=[t0,tf]; ts=[t0,t1,t2,……，tf]; ts=t0:k:tk; t0表示自变量的初值，tf为终值； x0表示函数的初值； options用于设定误差限，确省时误差限（相对误差为0。001，相对误差0.000001） 程序为：options=odeset(‘reltol’,rt,’abstol’,at) rt,at分别设定相对绝对误差； 用于解n个未知函数的方程组时，m文件中的待解方程组应以X的向量形式给出，x0亦然。 四：计算实际例： 把模型（1）改写成矩阵形式 1：对于给定数据r=1,d=0.5,a=0.1,b=0.02