《数学实验》第五讲.pptVIP

初始条件为x(0)=0, y(0)=0, 对于给定的a,b,d, ?进行计算。当x(t)满足 x(t) ? d + a t cos? 出现交点，则认为已击中目标。 如果t T, 则敌舰在打击范围内，可以发射。 ? d 例8 在导弹系统中设 a=90km/h, b=450km/h, T=0.1h. 求d, ?的有效范围? 解 有两个极端情形容易算出。 若? =0, 即敌舰正好背向行驶, 即x轴正向。 那么导弹直线飞行, 击中时间 t=d/(b-a)T 得d=T(b-a)=36km。 若? =?, 即迎面驶来, 类似地有d=T(a+b)=54km 一般地, 有36d54。 在线算法：对于测定的d 和?，可用 计算出t。如d=50, ? =?/2，写出M函数eg5_8fun，再用ode45求解 在线算法：eg5_8a 离线算法： 首先对于所有可能的d和?，计算击中所需时间，从而对不同? ，得d的临界值。具体应用时直接查表判断。 计算机模拟： 对于任意选定的参数a,b,d,T和?，利用MATLAB动画制作指令getframe和动画播放指令movie制作一个导弹追击敌舰演示工具 eg5_8b.m eg5_8c(90,450,30,0.3*pi,0.1) * * MATLAB数学实验 第五章 常微分方程 第五章 常微分方程 5.1 预备知识：常微分方程 5.2 解常微分方程的MATLAB指令 5.3 计算实验：Euler法和刚性方程组 5.4 建模实验: 导弹系统的改进 本章主要MATLAB指令 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,称为齐次的 5.1 预备知识：常微分方程 2.初等积分法 分离变量法等 3.常系数线性微分方程 线性常微分方程的解为一个特解和相应 的齐次微分方程通解的叠加。 齐次微分方程的解可用特征根法求得 例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) 4.初值问题数值解 数值解法:寻求解y(t)在一系列离散节点 t0 t1 … tn tf 上的近似值yk (k=0,1,…n)。 hk = tk+1 －tk 为步长，通常取为常量h 。 其中 数值解 准确解y(xn)未知, 数值解yn作为y(xn)的近似 高阶常微分方程初值问题可以化为一阶常微分方程组，已给一个n阶方程 y(n) = f (t, y, y’,…,y(n-1)) 设y1 = y, y2 = y’, …, y n= y(n-1), 化为一阶方程组 y0： 表示初值向量y0； t： 表示节点列向量 (t0 , t1 , … , tn)T； y： 数值解矩阵,每一列对应y的一个分量 若无输出参数，则作出图形 初值问题求解? 常用格式[t,y] = ode45(odefun, tspan, y0) odefun：表示f (t, y)的函数句柄或Inline函数 t是标量，y是标量或向量； tspan： 若为[t0, tf], 表示自变量初值t0和终值tf 若为[t0,t1,?, tn], 表示输出节点列向量 5.2解常微分方程MATLAB指令 完整格式 [t,y] = ode45(odefun, tspan, y0,options, p1,p2,?) options——为计算参数(如精度要求)设 置，默认可用空矩阵[]表示； p1,p2,? ——为附加传递参数，这时 odefun的表示为f (t, y, p1,p2, ?) odefun= inline(y-2*t/y,t,y); [t,y]=ode45(odefun,[0,4],1); [t,y] plot(t,y,o-) [t,y]=ode45(odefun,0:1:4,1);[t,y] 例2 解微分方程 y’ =