东华大学matlab数学实验课件Chapter6资料.pptVIP

为了保持对准目标，导弹轨迹切线方向应为 由上面两个方程得下列微分方程 初始条件为x(0)=0,y(0)=0, 对于给定的a,b,d, ?进行计算。 当x(t)满足x(t)+?d+atcos?，则认为已击中目标。 这里?代表允许的误差，因为敌舰是有一定大小的。如果tT,则敌舰在打击范围内，可以发射。 例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。 (1)在线算法：对于测定的d和?，可用下式计算t。 如d=50,?=?/2，写出M函数eg6_8fun.m(P109) 由于在T小时内，横坐标没有突破x=50,所以敌舰不在有效打击范围，应等近一些再发射。 (2)离线算法：首先对于所有可能的d和?，计算击中所需时间，从而对不同?，得d的临界值。具体应用时直接查表判断。 x(t)+?d+atcos? 取?=0.1，运行得临界曲线。使用时查询即可。 (3)计算机模拟： 6.6 习题 P121 1（1），（3）、2、3 第六章 常微分方程 6.1 预备知识：常微分方程 6.2 解常微分方程的 MATLAB命令 6.3 计算实验：Euler法和刚性方程 6.4 建模实验：导弹系统的改进 6.5 习题 6.1 预备知识：常微分方程 1、微分方程的概念 2、初等积分法 3、常系数线性微分方程 4、初值问题数值解 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、常系数线性微分方程 线性常微分方程的解为一个特解和相应的齐次微分方程通解的叠加； 齐次微分方程的解可用特征根法求得； 特解用常数变易法求得； 例1、求x’’+ 0.2x’+3.92x = 0的通解 解：特征方程为?2 + 0.2? +3.92=0 roots([1 0.2 3.92]) 求得共轭复根 ?+?i=-0.1?1.9774i, 通解为 x(t) = Ae?tcos(?t) +Be?tsin(?t) =Ae-0.1tcos(1.9774t)+Be-0.1tsin(1.9774t) A、B为任意常数 4、初值问题数值解 数值解法：寻求解y(t)在一系列离散节点 t0t1…tntf上的近似值yk(k=0,1,…n)。 hk=tk+1－tk为步长，通常取为常量h。 其中 高阶常微分方程初值问题可以化为一阶常微分方程组，已给一个n阶方程y(n)=f(t,y,y’,…,y(n-1)) 设y1=y,y2=y’,…,yn=y(n-1),化为一阶方程组 6.2 解常微分方程的 MATLAB命令 1、初值问题求解 2、解边值问题 1、初值问题求解 [t,y]=ode45(odefun,tspan,y0)－常用格式 odefun：表示f(t,y)的函数句柄或内嵌函数 其中t是标量，y是标量或向量； tspan：若为[t0,tf],表自变量初值t0和终值tf 若为[t0,t1,?, tn],表输出节点列向量 y0：表示初值向量y0； t： 表示节点列向量(t0,t1,…,tn)T； y： 数值解矩阵,一列对应y的一个分量 若无输出参数，则作出图形 完整格式： [t,y]=ode45(odefun,tspan,y0,options, p1,p2,?) options－为计算参数(如精度要求)设置，默认可用空矩阵[]表示； p1,p2,?－为附加传递参数，这时odefun的表示为f(t,y,flag,p1,p2,?) 例2 解微分方程 y’ = y?2t/y, y(0)=1, 0t4 解： odefun=inline(y-2*t/y,t,y); [t,y]=ode45(odefun,[0,4],1); [t,y]