(2016mod student)常微分方程sxsy5.pptVIP

* * * * * * * * * * * * * * * * 海军方面要求改进现有的舰对舰导弹系统。目前的电子系统能迅速测出敌舰的种类、位置以及敌舰行驶速度和方向，且导弹自动制导系统能保证在发射后任一时刻都能对准目标。 根据情报，这种敌舰能在我军舰发射导弹后T分钟作出反应并摧毁导弹。 现在要求改进电子导弹系统，使我舰能自动计算出敌舰是否在有效打击范围之内。 一、引例 导弹系统的改进 设我舰发射导弹时位置在坐标原点，敌舰在x轴正向d（km）处, 其行驶速度为 a (km/h), 方向与x轴夹角为? , 导弹飞行线速度b (km/h) 。 ? d 易知 t 时刻敌舰位置为 (d+atcos?, atsin?)。 敌舰 (x(t)，y(t)) 设t 时刻时导弹位置为(x(t)，y(t)) , 则 为了保持对准目标，导弹轨迹切线方向的斜率应为 ? d 敌舰 (x(t)，y(t)) 将上面两个方程联立，得下列一阶常微分方程组 化简后得一阶常微分方程组： 初始条件为x(0)=0, y(0)=0, 对于给定的a,b,d,? 进行计算。 当x(t)满足 x(t) + ? d + a t cos? 则认为已击中目标。 ?代表允许的误差，因为敌舰是有一定大小的。如果t T,则敌舰在打击范围内，可以发射。 二、数学理论复习: 常微分方程（组） 1、常微分方程的概念 常微分方程: f (t,y,y’,y,…,y(n))=0 常微分方程组: 由若干个常微分方程联立而成的方程组 线性常微分方程: y(n) + a1 (t) y(n-1) +…+ an-1 (t) y’+an (t) y = b(t) 若ai (t) (i =1, …,n) 与t无关, 称为常系数的 若b(t)=0,称为齐次的 2、初等积分法 3、常系数线性常微分方程 导出方程（齐次） y(n) + a1 y(n-1) +…+ an-1 y+an y = 0 通解=特解+相应的齐次方程通解 常系数齐次微分方程的解可用特征根法求得 例1 求x+ 0.2 x+3.92x = 0的通解 解： 特征方程为?2 + 0.2? +3.92=0 ? roots([1, 0.2,3.92]） 求得共轭复根 ?+i ? =-0.1?1.9774i, 通解为 x(t) = Ae?t cos(?t) +Be?t sin(?t) 三、常微分方程（组）数值解法：Euler法 数值解法:寻求解y(t)在一系列离散节点 t0 t1 … tn 上的近似值yk (k=0,1,…n)。 hk = tk+1 - tk 为步长，通常取为常量h 。 Euler法：在节点处用差商近似代替导数 Euler格式 k=0,1,2…,n-1 教材上的m函数euler.m给出Euler法计算程序 使用格式为 [tOut,yOut] = euler(ypfun, tspan, y0,h) ypfun: 表示f(t, y)的m文件名 tspan=[t0, tn]: 表示自变量初值t0和终值tn y0: 表示初值列向量y0。 h是步长。 输出列向量tOut: 表示节点 (t0 , t1 , … , tn) 输出矩阵yOut: 表示数值解，每一列对应 y的一个分量 例2 解方程 y = y-2t/y, y(0)=1, 0t1 四、使用MATLAB命令求解 解: 先写m函数：eg5_2fun.m ?[t,y]=euler(eg5_2fun,[0,1],1,0.1) 1、数值解 [tOut,yOut] = ode45(yprime, tspan, y0) 四、五阶Rung-Kutta-Felhberg法，用法与euler.m相同。 若无输出参数，则作出图形。 ode23与ode45类似，只是精度低一些。 2、常微分方程（组）解析解（符号运算） 例3 (1)求y=ay+b的通解；(2)求解例2 (3)高阶方程 y=cos(2x)-y, y(0)=1, y(0)=0 (4)方程组 f =f+g, g=-f+g, f(0)=1, g(0)=2 s=dsolve(‘方程1’,‘方程2’,?,‘初始条 件 1’,‘初始条件2’,?,‘自变量’) 方程均用字符串方式表示,自变量缺省值为t 导数用D表示,2阶导数用D2表示,以此类推。 s返回解析解。方程组情形, s为符号型结构数组。 3、刚性(stiff)方程组解法 刚性方程组解法ode15s，使用格式同ode45 解 先将方程写为m函数eg5_4fun.m ? [t,y]=ode15s(eg5_4fun,[0,400],[2,1]’)