三点法 比例导引法 课程设计.docx

沈阳理工大学装备工程学院综合课程设计说明书绪 论导弹制导规律即导引律是空战中实现战机追踪拦截导引的火控系统关键技术之一。导引律的选择对导弹能否精确打击目标至关重要,它根据双方的相对位置、速度和加速度等基本信息导引载机接近目标，实施攻击。针对机动目标的攻击导引技术是导引律研究的重点，这是因为实际空战中双方采取机动方式对抗，目标的机动往往难于预测。为此人们从不同角度采用不同的理论和方法研究针对机动目标的导引律，提高导引性能。本章对导弹导引律的研究状况进行了综述，以期为导弹制导与控制及相关问题研究提供参考。1反坦克导弹实际目标的运动特性是无法预先确定的。在导弹设计或研究问题时，往往对目标的运动规律进行假设。如假设目标平直等速飞行，或等速盘旋飞行等。导弹的飞行速度的变化，则由弹体结构、空气动力外形和发动机特性来确定。而决定理想弹道最重要的因素是导引法的选择。对于遥控导弹来说，一个好的导引法应具有以下特点：（1）由导引法确定的理想弹道必须通过目标；（2）理想弹道各点的法向加速度值在目标遭遇区附近应非常小；（3）目标机动飞行时，对遭遇区附近的弹道法向加速度的影响愈小愈好；（4）实现导引法的误差公式要简单，在技术上要易于实现，并具有一定的抗干扰性。目前，都是以这四项标准来衡量导引法的优劣。为此，需要深入研究导弹在各种导引法情况下所确定的理想弹道的运动特性。同时，在自寻的制导中，有三种经典导引方式，分别是追踪法、平行接近法、比例导引法。追踪法是指导弹在飞向目标的过程中，导弹的运动速度飞向始终指向目标。其优点在于制导系统工程实现容易，但缺点是导弹迎击目标或攻击近距离高速目标时，弹道弯曲严重，需要较大的法向过载。平行接近法是指导弹在飞向目标的过程中，目标视线在空间始终保持平行（即目标视线角保持不变）,采用平行将近导引律时，不需要太大的法向过载，导弹在空间飞行直至命中目标的飞行时间较短，这是它的优点，但这种导引规律实现起来很困难。比例导引法是在自寻的导弹上采用较多的一种导引规律，它是指在导弹飞向目标的过程中，导弹速度方向的变化率与目标视线的变化率成比例，这种导引规律易于工程实现，同时通过选择合适的导引比，就不会需要太太的法向过载，对不同机动特性的目标适应能力也较强，因此广泛应用于各类导弹上。当今和未来战场的大纵深、立体化、信息化、密集综合火力支援以及快速机动等突出特点。因此，未来战争对反坦克导弹的首发命中率、抗干扰能力、全天候作战能力等提出了更高的要求。反坦克导弹的发展趋势是“发射后不用管”、全天候作战能力、自动目标识别以及较强的抗干扰能力等。目前，激光制导反坦克导弹采取的制导方式主要有两类：寻的制导和指令制导。寻的制导有主动和半主动之分，迄今为止，主要是半主动式。激光半主动制导是用单独的激光目标指示器照射目标，弹上导引头接收目标发射的激光，经过信号处理形成控制指令控制导弹的飞行。激光半主动制导能实现间接瞄准，可采用准比例导引法，导弹弹道特性好，对目标机动有一定的适应性。比例导引法实现了打了不用管的作战要求，这就提高了武器系统的生存能力，同时增加了其攻击效率4。虽然，导弹在实际作战中是三维机动，但平面问题的比例导引关系是研究导弹空间运动的比例导引规律的基础，对于初学者来说，是研究比例导引关系的重要手段。本次课程设计就是通过采用matlab仿真软件，对主要弹道参数对反坦克导弹比例导引法导引弹道影响仿真。通过仿真结果揭示比例导引法的优越性，及其弹道特性。31 三点法1.1 三点法简介三点法导引是指导弹在攻击目标的导引过程中，导弹始终处于制导站与目标的连线上。如果观察者从制导站上看目标，则目标的影像正好被导弹的影像所覆盖。因此，三点法又称目标覆盖法或重合法。对导引弹道的研究是以经典力学为基础的。在导弹和制导系统初步设计阶段，为了简化研究，通常采用运动学分析方法。为此，我们通常作如下假设：导弹、目标和制导站的运动视为质点运动；制导系统的工作是理想的，无惯性，无延迟；导弹和目标始终在攻击平面内运动；导弹速度是时间的已知函数，目标和制导站的运动规律也是已知的。三点法属于遥控指令制导中的一种导引方法。遥控指令制导与自寻的制导的不同之处在于：导弹的运动受设在弹外的制导站的控制，其运动规律不仅取决于目标的运动，而且还与制导站的运动规律有关。研究遥控弹道时，既要考虑导弹相对目标的运动，还要考虑制导站的运动对导弹运动的影响。遥控导引时，导弹和目标的运动参数都由制导站来测量。2 1.2 三点法导引的导引关系方程和运动学方程组1.2.1 导引关系方程由于导弹始终处在目标和制导站的连线上，所以导弹与制导站连线的高低角 和方位角 要始终与目标与制导站连线的高低角 和方位角 相等，由此得出三点法导引关系方程为：图1.1 三点法1.2.2 运动学方程组为讨论方便，认为导弹在垂直平面内飞行，雷达坐标系 。原点取在制导站； 轴指向迎着制导站飞来的目标 ； 轴同轴垂直朝上。 轴与地平面的夹角，称为高低角。轴为地面上的某一参考线，它是度量雷达方位角的基准线。如下图1.2所示：图1.2 反坦克导弹三点法导引则三点法导引的相对运动方程组为： （1-1）其中： 目标到制导站距离； 导弹到制导站距离； 目标速度； 导弹速度； 目标高低角、目标前置角； 导弹高低角、导弹前置角； 、 与基准线之间的夹角。1.3 初始条件设坦克作水平等速直线运动，如图1.2所示，m/s ，反坦克导弹按三点法拦截目标，并作等速飞行，m/s 。攻击平面为一水平面，制导站静止。导弹开始导引瞬间的攻击条件为m ， ， 。选取地面坐标系oyz ，原点o与制导站重合，oz轴平行于目标的运动方向（如图1.2所示）。将方程组（1-1）改写成便于数值积分的形式，即（1-2）2 反坦克导弹三点法导引弹道仿真结果及分析2.1 给定条件下的弹道曲线图2.1水平面攻击弹道曲线图2.1是反坦克导弹三点法导引时的水平面弹道曲线，根据从工作空间的程序结果可以得到弹目交会距离约为=4417m。2.2 的变化曲线图2.2 变化曲线从图2.2可以看出视线与基准线间的夹角随着时间t是不断增大的。2.3 的变化曲线图2.3 变化曲线从图2.3可知，导弹距制导站的距离随时间t不断增大，弹目交会在 时，距离为4417m。2.4 的变化曲线 图2.4 变化曲线从图2.4可知，目标距制导站的距离随时间t不断减小。2.5 主要弹道参数对导引弹道的影响仿真由之前的导弹运动学方程组可知，某些弹道参数或初始条件不同会得到不同的弹道，如下就部分参数对弹道的影响进行了仿真分析。2.5.1 导弹与目标的速度比对弹道的影响图2.5不同速度比下的弹道曲线从图2.5中可知，在弹目速度比不同的情况下，弹道的弯曲程度有明显的差异。当速度比较小时，弹道会比较弯曲，从而导弹的法向过载会相应增大，而速度比较大时，弹道则更为平缓。此外，速度比不同对弹目的交会时间也会有影响，速度比越大，导弹击中目标所需的时间越短。2.5.2 对导引导弹的影响图2.6不同下的弹道曲线从图2.6中可知不同的初始瞄准角对导引弹道有不同程度的影响，在其他初始条件不变的前提下，越大，弹道越弯曲，反之弹道越平滑。2.6 结论分析导弹和目标运动轨迹，建立二者的运动学方程组，然后分析运动学方程组，简化导弹与目标的运动模型，使用Matlab软件对运动学方程组进行编程仿真。根据仿真结果求出弹目遭遇时间，弹目交会距离为4417m，视线与基准线间的夹角随着时间t是不断增大的；导弹距制导站的距离随时间t不断增大；目标距制导站的距离随时间t不断减小。根据主要弹道参数对导引弹道的影响仿真结果得出主要弹道参数对弹道的影响规律如下：（1）导弹与目标的速度比对弹道的影响：在弹目速度比不同的情况下，弹道的弯曲程度有着显著的差异。当速度比较小时，弹道比较弯曲，从而导弹的法向过载相应增大，而速度比较大时，弹道则更为平缓。此外，速度比不同对弹目的交会时间也会有影响，速度比越大，导弹击中目标所需的时间越短。（2）初始瞄准角对导引弹道的影响：在其他初始条件不变的前提下，越大，弹道越弯曲，遭遇时间越长，反之弹道越平滑。对导弹的弯曲程度影响不大，主要是对遭遇时间的影响比较明显。2.6.1 三点法的优点三点法导引的最显著的优点就是技术实施简单，抗干扰性能好。对涉及低速目标，射击从高空向低空滑行或俯冲的目标；被射击的目标释放干扰，导弹制导站不能测量到目标距离信息时；制导雷达波束宽度或扫描范围很窄时，在这些范围内应用三点法不仅简单易行，而且其性能往往优于其他一些制导规律。2.6.2 三点法的缺点（1）弹道较弯曲，迎击目标时，越是接近目标，弹道越弯曲，需用法向过载越大，命中点的需用法向过载最大。在接近目标过程中，可能出现导弹可用法向过载小于需用法向过载，导致导弹脱靶。（2）动态误差难以补偿。动态误差是指制导系统过渡过程中复现输入时的误差。由于目标机动所引起的动态误差难以补偿，往往会形成偏离波束中心线十几米的动态误差。（3）三点法导引时，发射角小，导弹离轨时飞行速度很小，此时操纵效率低，空气动力提供的法向力也比较小，导弹离轨时可能下沉，为了克服这一缺点，采用小高度三点法，提高初始段弹道高度。3比例导引法3.1 比例导引法简介比例导引法是指导弹飞行过程中速度向量V的转动角速度与目标视线的转动角速度成比例的一种导引方法。7反坦克导弹比例导引法导引数学模型是根据相对位置关系以及导弹与目标的运动状态建立的。所以，首先要画出目标、导弹与基坐标之间的关系图，而后才能写出微分方程组。3.2 比例导引法的数学模型3.2.1 目标、导弹与基准线的相对位置关系自动瞄准制导的相对运动方程实际上是描述导弹与目标之间相对运动的关系。如图3.1所示，假设在某一时刻，目标位于T点，导弹位于M点。链接M和T称为目标瞄准线。选取基准线Ax，它可以任意选择，它的位置的不同选择不会影响导弹与目标之间的相对运动特性，而只影响相对运动方程的繁简程度。一般选取目标的飞行方向为基准线方向最为简单。根据导引弹道的运动学分析方法，假设导弹与目标的相对运动方程可以用定义在攻击平面内的极坐标参数r、q的变化规律来描述。图3.1 导弹、目标与基准线的相对位置图其中，r弹道相对目标的距离。导弹命中目标是r=0.q目标与基准线之间的夹角称为目标线方位角（简称目标线角）。若从基准线逆时针转到目标线上时，则q为正。 分别为导弹、目标速度矢量与基准线之间的夹角，称之为导弹弹道角和目标航向角。分别以导弹、目标所在位置为原点，若由基准线逆时针旋转到各自的速度矢量上时，则为正，当攻击平面为铅垂面时，就是导弹倾角；当攻击平面为水平面时，就是导弹偏。 分别为导弹、目标速度矢量与目标之间的夹角，相应称之为导弹速度矢量前置角和目标速度矢量前置角（简称为前置角）。分别以导弹、目标为原点，若从各自的速度矢量逆时针旋转到目标线上时，则为正。53.2.2 建立相对的运动方程考虑到上图所示角度间的几何关系以及导引关系方程，就可以得到自动瞄准制导的相对运动方程组为： （1-3）上面方程中，1=0为描述导引方法的导引关系方程（或称理想控制关系方程）。在自动瞄准制导中常见的导引方法有：追踪法、平行接近法、比例导引法等，相应的导引关系方程为 追踪法： ； 平行接近法： ； 比例导引法： 。上述方程组中： （t）， （t）， （t）为已知，方程组中只含有5个未知参数；r(t),q(t), , (t), (t),因此方程是封闭的，可以求得确定解。根据r(t),q(t),可获得导弹相对目标的运动轨迹，称为导弹的相对弹道（即观察者在目标上所观察到的导弹运动轨迹）.若已知目标相对地面坐标系（惯性坐标系）的运动轨迹后，则通过换算可获得导弹相对地面坐标系的运动轨迹绝对弹道。4 弹道仿真弹道仿真是根据建立的数学模型编写程序，然后采用相应的软件进行仿真，得到各种条件下的弹道。在仿真过程中，做如下基本假设：（1）对于拦截弹在其射程范围内，不考虑地球曲率和自转的影响；（2）动能拦截器为刚体，忽略弹体弹性；（3）导弹具有理想的控制，制导系统；（4）大气静止不计风的影响。4.1 反坦克导弹比例导引弹道仿真按比例导引时，导弹与目标之间的相对运动方程组为：（1-4）设目标作水平等速直线运动，VM =27m/s，导弹等速飞行，VD =270m/s，比例系数k=4，攻击平面为一水平面如图4.1，设初始条件为：R0=4700m，q0=70， =45。选取基准线Az平行于目标的运动方向，根据上述已知条件得出导弹与目标的相对运动方程为： （1-5）确定绝对弹道，所选地面坐标系Oxyz的原点与导弹初始位置重合，弹道参数为（x，z），其表达式为 （1-6）根据上述方程组及任务书给出的参数，编写Matlab程序进行弹道仿真，根据不同的参数对导弹进行多次仿真，并观察不同的仿真结果。坦克目标水平等速直线运动， VM=27m/s ，反坦克导弹制导系统使用比例导引法制导拦截目标，速度V=270m/s。设攻击平面为一水平面，制导站静止，导弹开始引导瞬间的攻击条件为目标与基准线原点相距 ，导弹与基准线原点相距 ，导弹弹道角 ，视线角 导引比例系数k=4，仿真结果如图4.1图 4.1 比例k=4对应的弹道曲线4.2 主要弹道参数对导引弹道的影响比例导引法主要弹道参数是指对导弹弹道影响较大的参数，如导弹速度，目标速度，采用的目标线角，比例导引系数等，现主要分析目标导引系数k以及导弹与目标的速度比p这两个参数对弹道的影响。4.2.1 比例导引系数对导引弹道的影响 当k取不同值时的弹道如图2.2所示,从图中我们可以清楚的看到，当k=1时，弹道弯曲较大。当k=9时，弹道几乎成了一条直线，由此可以得出结论：当k较小时，导引弹道越弯曲，k较大，导引弹道几乎是平直的。6（1）当k=1，即 ，这就使导弹速度矢量与视线矢量相同的角速度转动，则弹道初始段的曲率半径较大；但随着弹目距离的缩小，q变大，而使弹道曲率半径变小。由于，这就意味着导弹速度矢量永远都与视线相重合，相当于直接追踪法。（2）当k=4，即 ，这就使导弹速度矢量以4倍于视线角速度转动，这就导致视线角速度q减小，从而使 减小，即保证导弹按q=0，q为常数的平行接近导引。（3）当k=6或9，这就使导弹速度矢量以较高倍于视线角速度转动，这将导致视线角速度q迅速减小，从而使 迅速减小，即保证导弹按q=0,q为常数的平行接近导引。图4.2 不同比例系数对应的导引弹道曲线由于 ，k1,表明这样的比例导引的初始段的导弹过载要大，及弹道曲率半径小；而到一定的时间，弹道过载为零，即弹道的曲率半径为无穷大，弹道为一直线。4.2.2 导弹与目标的速度比对导引弹道的影响假设目标的运动速度不变VM =27m/s，当导弹速度分别为270m/s、405m/s、675m/s时，取比例系数k=4，从图中可以看出，对于不同的速度比p，仿真得到的导引弹道曲线也是不同的，主要表现在以下几个方面：（1）p越小，弹道初始段曲率半径越大，弹道相对比较平直；p越大，弹道初始段曲率半径越小，弹道越弯曲。（2）p越大，弹道曲线与水平线夹角越大，导弹与目标相遇距离越近，相遇时间越短；p越小，弹道曲线与平面线的夹角越大，导弹与目标相遇距离越远，相遇时间越长。图4.3 不同速度比p的导引弹道曲线通过以上的分析，可以得出结论，仿真结果与实际情况基本相符。4.3 弹道参数对遭遇时间的影响的仿真4.3.1 目标线角对遭遇时间的影响在其他条件不变的情况下，改变初始目标线角 ，通过仿真可以得到遭遇时间随目标线角 变化的曲线，如图4.4所示图4.4不同目标线角对应的遭遇时间从图中可以看出目标线角q0越大，遭遇时间越长。4.3.2 弹目距离对遭遇时间的影响在其他条件不变的情况下，改变弹目距离，通过仿真可以得到遭遇时间随目标线角变化的曲线，如图4.5所示图 4.5不同弹目距离对应的遭遇时间4.4 比例导引弹道仿真结果4.4.1 仿真结果分析通过以上的分析我们可以看到，导弹与目标速度比p的变化对导引弹道有较大的影响。同时，比例导引系数k的大小，对整个弹道曲率半径的影响也是很大的。P越大，弹道初始段曲率半径越大，弹道相对比较平直，弹道曲线与水平线夹角越小，导弹与目标相遇距离越近，相遇时间越短；p越小，弹道初始段曲率半径越小，弹道越弯曲，弹道曲线与平面线的夹角越大，导弹与目标相遇距离越远，相遇时间越长。比例导引系数k越小，导引弹道越弯曲，导引弹道越弯曲；k越大，导引弹道越平直。由导引方程和以上仿真分析可以看出，k=1, 时为追踪法的弹道； 为有限量，则 ，则为平行接近法的导引弹道；时为比例导引法的弹道。换句话说，比例导引法是介于追踪法和平行接近法之间的一种导引方法。比例导引法的弹道特性也介于追踪法和平行接近法两者之间，随着比例系数k的增大，导引弹道越加平直，需用法向过载也就越小。从上面讨论可知，比例系数k的大小直接影响弹道能否直接命中目标。选择合适的k值除考虑这两个因素外，还需要考虑结构强度所允许的承受过载能力，以及制导系统能否稳定地工作等因素。因此k值选择应满足：（1）k值受可用法向过载的限制上式限制了比例系数k的下限值。但其上限值如果取得过大，由 可知，即使 值不太大，也可能使需用法向过载很大。导弹在飞行中的可用法向过载受到最大舵偏角的限制。若需用法向过载超过可用法向过载，则导弹不能沿比例导引弹道飞行。因此可用过载限制了k值上限。（2）k值下限应满足 收敛的条件 收敛使导弹在接近目标的过程中目标线的旋转角速度 不断减小，相应的需用法向过载也不断减小。 的收敛条件为：这就限制了k的下限值。综合考虑上述因素，才能选择出一个合适的k值。它可以是个常数，也可以是个变数。（3）制导系统的要求 如果比例系数K选得过大，那么外界干扰信号的作用会被放大，这将影响导弹的正常飞行。因此，从制导系统稳定工作的角度出发，K值的上限值也不能选得太大。综合考虑上述因素，才能选择出一个合适的K值。它可以是一个常数，也可以是一个变数。一般认为，K值通常在36范围内。 4.4.2 过载问题在弹体结构和控制系统设计中，常需要考虑导弹在飞行过程中能够承受的过载。根据战术技术要求的规定，飞行过程中过载不得超过某一数值。这个数值决定了弹体结构和弹上各部件能够承受的最大载荷。为保证导弹能正常飞行，飞行中的过载也必须小于这个数值。在导弹设计过程中，经常用到需用过载和可用过载的概念，下面分别加以叙述。（1）需用过载：所谓需用过载是指导弹按给定的弹道飞行时所需要的法向过载，用 nR 表示。导弹的需用过载是飞行弹道的一个重要特性。需用过载必须满足导弹的战术技术要求，例如，导弹要攻击机动性强的空中目标，则导弹按一定的导引规律飞行时必须具有较大的法向过载(即需用过载)；另一方面，从设计和制造的观点来看，希望需用过载在满足导弹战术技术要求的前提下越小越好。因为需用过载越小，导弹在飞行过程中所承受的载荷越小，这对防止弹体结构破坏、保证弹上仪器和设备的正常工作以及减小导引误差都是有利的。（2）可用过载：当操纵面的偏转角为最大时，导弹所能产生的法向过载称为可用过载。它表征着导弹产生法向控制力的实际能力。若要使导弹沿着导引规律所确定的弹道飞行，那么，在这条弹道的任一点上，导弹所能产生的可用过载都应大于需用过载。4.5 结论比例导引法具有平行接近法的优点，即导弹过载小，实现比例导引法的装置比较简单，因此这种导引方法得到了广泛的应用。比例导引法在满足 的条件，逐渐减小，弹道前段较弯曲，充分利用了导弹的机动能力。弹道后端较平直，使导弹具有较充裕的机动能力。只要导引比例系数k、导弹速度矢量前置角初始值、目标线方位角初始值 、弹目速度比p等参数组合适当，就可以是全弹道上所需法向过载均小于可用法向过载，因而能实现全向攻击。另外，与平行接近法相比，对瞄准发射时的初始条件要求不严。在技术上只需测量 ，实现比例导引比较容易。比例导引法的弹道也比较平直。因此自动瞄准制导的导弹都广泛采用比例导引法。比例导引律经过几十年的发展，在基础理论、实际应用中都取得了发展，尤其是其改进形式更是得到了广泛应用。为了满足高科技条件下的战争需求，应以现有的实用制导规律为基础，扬长避短，同时结合实际需要，利用最新的相关科学知识与技术，对其作进一步的改进研究和优化设计。所谓法向过载就是垂直于导弹速度的加速度，在仿真中是由弹目运动参数计算得到的，不过能否在工程上实现得看硬件。就是舵板和控制舵板的电动机能产生多大的力矩。所以法向过载越平稳越好、越小越好，这样电动机的能量消耗就小。4.5.1 比例导引法的优点（1）可以得到较为平直的弹道；（2）在满足的条件下,逐渐减小，弹道前段较弯曲，充分利用了导弹的机动能力；（3）弹道后段较为平直，导弹具有较充裕的机动能力；（4）只要等参数组合适当，就可以使全弹道上的需用过载均小于可用过载，从而实现全向攻击。（5）与平行接近法相比，它对发射瞄准时的初始条件要求不严，在技术实施上是可行的，因为只需测量 。因此，比例导引法得到了广泛的应用。4.5.2 比例导引法的缺点比例导引法还存在明显的缺点，即命中点导弹需用法向过载受导弹速度和攻击方向的影响。为了消除比例导引法的缺点，多年来人们一直致力于比例导引法的改进，研究出了很多形式的比例导引方法。例如，需用法向过载与目标视线旋转角速度成比例的广义比例导引法，其导引关系式为 或 ，式中，K1,K2为比例系数； 为导弹接近速度。参考文献1 潘云芝，潘传勇.导引律研究现状及其发展，科技信息.20092 任 波，胡晓阳，王 欣.弹箭制导原理与控制技术.3 钱杏芳，林瑞雄，赵亚男.导弹飞行力学M.北京理工大学出版社，2000.4 雷虎民.导弹制导与控制原理M.北京国防工业出版社，2006.5 高尚，比例导引理想弹道仿真J.计算机工程与设计，2003,24（4）.6 金永德，崔乃刚，关英姿.导弹与航天技术概论M.哈尔滨工业大学出版社，20027 百度百科，比例导引法/link?url=B8SxDvkPI25vVYl0O_BKwfdgQ0C2F3zLRz1ofxydFitEgZ1xMkMBJPerUDswemTcjlah6szX_SJHR92aIFzAf_附 录1 三点法相关程序1 运动微分方程组的自定义函数子程序：function f=dyfun(t,x) vt=27; v=270; faiv=x(3)+asin(vt/v)*(x(1)/x(2)*sin(x(3); f(1)=v*cos(x(3)-faiv); f(2)=-vt*cos(x(3);f(3)=-v/x(1)*sin(x(3)-faiv);f=f(:);主程序（方程组的求解程序）：clear;clc;x0=50,3500,70/180*pi; t,x=ode45(dyfun,0 26.10,x0); rm=x(:,1);rt=x(:,2);qm=x(:,3);ym=rm.*sin(qm);zm=rm.*cos(qm);yt=rt.*sin(qm);zt=rt.*cos(qm);plot(zm,ym,zt,yt,-); xlabel(Z/m);ylabel(Y/m);title(水平攻击面内的弹道曲线);legend(导弹弹道,目标航迹,Location,SouthEast);gtext(T0;T1);figure plot(t,qm); xlabel(t/s);ylabel(qm/rad);title(qmt曲线);figureplot(t,rm); xlabel(t/s);ylabel(rm/m);title(rmt曲线);figureplot(t,rt); xlabel(t/s);ylabel(rt/m);title(rtt曲线);2 分析不同速度比c的程序子程序：function f=csfz(t,x) vt=27; faiv=x(3)+asin(1/x(4)*(x(1)/x(2)*sin(x(3);f(1)=vt*x(4)*cos(x(3)-faiv);f(2)=-vt*cos(x(3);f(3)=-vt*x(4)/x(1)*sin(x(3)-faiv);f(4)=0;f=f(:); 主程序：clear;clc;c=3 5 7 9; tf=85.5 42.95 32.45 28.95; for i=1:4 x0=50,3500,70/180*pi,c(i); t,x=ode45(csfz,0 tf(i),x0); rm(:,i)=x(:,1); rt(:,i)=x(:,2);qm(:,i)=x(:,3);ym(:,i)=rm(:,i).*sin(qm(:,i); zm(:,i)=rm(:,i).*cos(qm(:,i);yt(:,i)=rt(:,i).*sin(qm(:,i);zt(:,i)=rt(:,i).*cos(qm(:,i);endfigure plot(zm(:,1),ym(:,1),zt(:,1),yt(:,1),-,zm(:,2),ym(:,2),zt(:,2),yt(:,2),-,zm(:,3),ym(:,3),zt(:,3),yt(:,3),-,zm(:,4),ym(:,4),zt(:,4),yt(:,4),-);xlabel(Z/m);ylabel(Y/m);title(不同速度比对应的弹道曲线);legend(导弹弹道,目标航迹,Location,SouthEast);gtext(c=3;c=5;c=7;c=9);gtext(T0;T1;T2;T3;T4);3 分析不同瞄准角qM0的程序clear;clc;qm0=45 55 60 80; tf=25.55 26.00 26.35 26.85; for i=1:4 x0=50,3500,qm0(i)./180*pi; t,x=ode45(dyfun,0 tf(i),x0); rm(:,i)=x(:,1); rt(:,i)=x(:,2);qm(:,i)=x(:,3);ym(:,i)=rm(:,i).*sin(qm(:,i);zm(:,i)=rm(:,i).*cos(qm(:,i);yt(:,i)=rt(:,i).*sin(qm(:,i);zt(:,i)=rt(:,i).*cos(qm(:,i);endfigure plot(zm(:,1),ym(:,1),zt(:,1),yt(:,1),-,zm(:,2),ym(:,2),zt(:,2),yt(:,2),-,zm(:,3),ym(:,3),zt(:,3),yt(:,3),-,zm(:,4),ym(:,4),zt(:,4),yt(:,4),-);xlabel(Z/m);ylabel(Y/m);title(不同初始视线角对应的弹道曲线);legend(导弹弹道,目标航迹);gtext(qm=45circ;qm=55circ;qm=60circ;qm=80circ);gtext(T1;T2;T3;T4);附 录2 比例导引法相关程序1 分析不同比例导引系数的程序主程序：clear;clc;k=1 4 6 9; tf=50 44 40 40; for i=1:4 x0=4700,1.202,k(i); t,x=ode45(kfun,0 tf(i),x0); zt(:,i)=4700*sin(70*pi/180);xt(:,i)=4700*cos(70*pi/180)+12.*t(1:60); xm(:,i)=xt(:,i)-x(1:60,1).*cos(x(1:60,2);zm(:,i)=zt(:,i)-x(1: