<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>Trajectory Modeling Calculation and Maneuverability Analysis of Proportional Navigation Method</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Xingwang Li</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Shuangjie Liu</string-name>
          <email>shuangjieliu@126.com</email>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Yongping Hao</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Shenyang Ligong University</institution>
          ,
          <addr-line>Shenyang 110159</addr-line>
          ,
          <country country="CN">China</country>
        </aff>
      </contrib-group>
      <fpage>72</fpage>
      <lpage>78</lpage>
      <abstract>
        <p>Aiming at the ballistic characteristics and related performance indexes of a certain type of antiaircraft missile, the differential equations of missile-target relative motion are established. Under different given initial conditions, the Runge-Kutta method is used to solve the system of micro-equations numerically. Thereby, the ballistic curve, ballistic parameters and performance index parameters of the corresponding missile interception are obtained. In the simulation, different missile/target maneuvering conditions, route shortcuts and interception methods are verified and calculated to measure the missile's ballistic characteristics and performance indicators comprehensively. The simulation results show that when the proportional guidance method is used, under different missile/target maneuvering conditions and interception methods, the greater the target maneuvering degree, the smaller the curvature radius of the missile's ballistic curve, and the greater the required normal overload and miss distance. The ballistic interception trajectory and the interception interception trajectory are relatively straight, and the required normal overload and misses are relatively small; when intercepting in the air, the radius of curvature of the ballistic trajectory is the smallest, and the required normal overload and misses are the largest.</p>
      </abstract>
      <kwd-group>
        <kwd>1 Ballistic characteristics</kwd>
        <kwd>Runge-Kutta method</kwd>
        <kwd>Proportional guidance，Normal overload， Miss distance</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>1. Introduction</title>
      <p>
        The power of the missile mainly depends on the guidance accuracy of the missile and the number of
warhead ammunition, and the guidance rules based on previous research are often referred to as
conventional guidance rules. There are three main induction rules: tracking method, line-of-sight
display induction method and proportional guidance method.The proportional guidance method is
simple, technically easy to implement, does not require too much information, its trajectory is relatively
flat, mechanical targets and low-altitude flying targets can be intercepted, the import accuracy is high,
and it is widely used[
        <xref ref-type="bibr" rid="ref1">1</xref>
        ].
      </p>
      <p>
        Generally speaking, in the guidance system of a missile, the role of the guidance law is to ensure
that the missile flies and hits the target's kinematic trajectory. The guidance law also determines the
flight trajectory of the missile center of mass in space.So, how to choose a suitable guidance method to
improve the guidance performance of the missile, improve the hit accuracy of the missile, and reduce
the amount of misses, has always been our special attention[
        <xref ref-type="bibr" rid="ref2 ref3">2-3</xref>
        ]. Some scholars briefly introduced the
commonly used missile guidance laws in their articles, took the proportional guidance method as an
example to conduct a two-dimensional plane simulation of the guidance ballistic, and analyzed the
proportional guidance coefficient, the ratio of the missile and the target speed, etc. The influence of the
change of the main ballistic parameters on the ballistics, some suggestions for improving the missile
guidance and control system are put forward[
        <xref ref-type="bibr" rid="ref4">4</xref>
        ].
      </p>
    </sec>
    <sec id="sec-2">
      <title>2. Vertical plane trajectory model</title>
      <p>This paper mainly studies the motion characteristics of the missile in the vertical plane. First, the
missile's full ballistic equation is decomposed on the vertical plane, and the dynamic mathematical
models of 10 vertical planes are obtained. Secondly, the dynamic mathematical model of the vertical
plane is replaced by the kinematic mathematical model of the vertical plane, that is, the change of the
tangential velocity and the normal maneuvering law of the missile are replaced by the known variation
law of the missile velocity movement. Then, the Runge-Kutta method is used to solve the two
differential equations about the sight angle of the target and the distance of the projectile, and the
simulation model of the ballistic trajectory is normalized.</p>
    </sec>
    <sec id="sec-3">
      <title>2.1. Dynamic mathematical model</title>
      <p>Usually the movement of the missile can be decomposed into the movement of the missile's vertical
plane and the movement of the missile's horizontal plane. If we require a missile to fly at high altitude
and high speed in the vertical plane, then a horizontal lateral force in the horizontal direction should be
equivalent to at zero. If we choose the ground coordinate system Ax axis is located in the vertical plane,
at this time, β, γ, etc. are correspondingly regarded as zero. When a particle in the vertical plane moves,
the missile only performs a parallel motion and a rotational motion along the axis and the OZ1 axis in
the rotational direction of the center of the inner particle on the vertical plane. At this time, z=0 ,   =
0,   = 0. When the missile moves along the vertical plane, the external forces of the missile are only:
the thrust P of the engine, the air resistance X, the lift Y, and the gravity G.Such as Formula 1.</p>
    </sec>
    <sec id="sec-4">
      <title>2.2. Kinematic mathematical model</title>
      <p>In the preliminary research stage of our anti-aircraft missile and guidance system technology, in
order to further simplify the research, we assume that the missile, the target and the guidance station
always move in the same fixed plane, that is, on the attack plane, which may be a vertical plane, or a
may be horizontal or inclined.</p>
      <p>We choose the attack plane as the plumb plane for research.According to the kinematic analysis
method of guided ballistics, it is assumed that a relative motion parameter equation between the missile
and the attack target can be defined by the motion parameters such as r, q in polar coordinates in a plane
called the missile attack force. The relative transformation motion law of , etc. can be used to describe
it accurately.
= − 
 =  − 
 1 = 0
 4 = 0
−  − 
+  −</p>
      <p>The relative motion equation of automatic aiming guidance refers to the equation describing the
relative distance r and the rate of change of the target line angle q. According to the relative motion
relationship between the missile and the target shown in Figure 1, the relative motion equation can be
established directly. The missile velocity vector VM and the target velocity vector VT are decomposed
along the target line direction and the normal line direction, respectively. The component   
points to the target along the target line, which reduces the relative distance r. The component   
deviates from the missile, which increases the relative distance r. The normal component along the
 
 
around the missile's position as the origin, thereby reducing the target line angle q.
target line,</p>
      <p>, rotates the target line counterclockwise around the origin where the target is
located to increase the target line angle q. The component   
  rotates the target line clockwise</p>
      <p>We also consider the relative geometric motion relationship between the two motion angles shown
in Figure 1. When guided proportionally, when the attack plane is a vertical plane, σ is the ballistic
inclination angle θ. Therefore, the mathematical model of the kinematics of the missile in the vertical
plane can be written as Formula 2.</p>
      <p>Some parameters in the ballistic dynamics mathematical model of the vertical plane are very
confidential, so it is difficult for us to obtain these parameters. For example, the first two equations of
Formula 1 can know the lift coefficient and drag coefficient in the lift Y and the drag X, so the ballistic
is used. It is difficult to directly study the motion law of the missile with the dynamic mathematical
model. Therefore, we will study the kinematic mathematical model of the vertical plane ballistic
trajectory, and use the change law of the missile velocity movement in the kinematic model to replace
the change law of the tangential velocity of the first two equations of the equation group in the vertical
plane ballistic dynamic mathematical model. and the law of motion in the normal direction.
+  
+</p>
      <p>In general, since the right side of the equation is a nonlinear function of the motion parameters, the
missile motion equations are nonlinear first-order ordinary differential equations. In this paper, the
Runge-Kutta method is used to carry out high-precision calculations. This algorithm is a high-precision
program that is easy to edit and calculate. It is also possible to quickly change the calculation step size.</p>
      <p>By normalizing the kinematic mathematical model of the vertical plane trajectory, we can design the
simulation algorithm according to the simulation model and initial conditions, and then carry out the
simulation. After obtaining the trajectory and the trajectory of the target, the final analysis of the results
is carried out, including the calculation of the missed target amount and the normal overload.
2.3.2. Simulation result</p>
    </sec>
    <sec id="sec-5">
      <title>2.3. Simulation analysis</title>
    </sec>
    <sec id="sec-6">
      <title>2.3.1. Simulation example</title>
      <p>The target moves in a straight line with horizontal uniform acceleration, and the acceleration is
  .The missile moves at a uniform acceleration of   .</p>
      <p>An example of missile head-on interception: the initial target speed is 300m/s, the initial missile
speed is 600m/s, the initial missile target distance is 7000m, the initial target line of sight angle is 50 °,
the initial lead angle is -3 °and the relevant parameters are shown in Table 1.</p>
      <p>An example of missile tail pursuit and interception: the initial target speed is 300m/s, the initial
missile speed is 600m/s, the initial missile target distance is 2000m, the initial target line of sight angle
is 50 °, and the initial lead angle is 3 °and the relevant parameters are shown in Table 2.In addition, we
also assume that in the process of missile flying to the target, the rotation of the earth, the curvature of
the earth and the influence of wind are not considered, and the missile also has an ideal control and
guidance system.</p>
      <p>When the scale factor is 3, the flight time and miss distance of the three calculation examples of the
missile intercepting the target head-on are shown in the Table 3 below.The last three are data with a
scale factor of 5.The data of tail pursuit interception are shown in Table 4.</p>
      <p>From the simulation results, we know that in the case of only increasing the flight speed of the
missile, the flight time of the missile decreases with the increase of the speed of the missile to a certain
extent. The proportional coefficient K also has a certain effect on the flight time, but the effect is not as
great as the effect of increasing the missile flight speed on the flight time. In the case of only increasing
the target flight speed, the flight time decreases with the increase of the target speed to a certain extent.
Table 3
Results of head-on interception of targets</p>
      <p>Miss distance(s)
3.16
2.70
4.82
4.65
2.78
6.34</p>
      <p>As the speed of the missile increases, the trajectory of the missile will be straighter, and the radius
of curvature will be larger; as the speed of the target increases, the trajectory of the missile will be more
curved, and the radius of curvature will be smaller; as the proportional coefficient increases, the
trajectory of the missile will be more Straight, the radius of curvature will be larger. In addition, a
relatively curved front section of the ballistic trajectory can be obtained. In the front section of the
missile trajectory, the proportional guidance method can reasonably and fully utilize the
maneuverability of the missile.</p>
    </sec>
    <sec id="sec-7">
      <title>3. Acknowledgements</title>
    </sec>
    <sec id="sec-8">
      <title>4. Conclusion</title>
      <p>When the missile speed is constant, the normal overload will decrease to a certain extent with the
increase of the proportional coefficient K. As shown in Figure 4.At the same time, the miss distance is
also an important index for us to study whether the missile can hit the target accurately. With the
increase of target maneuvering acceleration, the miss distance will also increase, and the missile
maneuvering acceleration and scale coefficient will also have different effects on the miss distance.The
simulation data of miss distance are shown in Table 3 and Table 4.</p>
      <p>The simulation diagram of the relevant tail chasing intercepting target is similar to the simulation
diagram of the head-on intercepting target, which will not be explained here.</p>
      <p>Thanks to the scholarly monographs cited in this article and to my teachers for their guidance.</p>
      <p>From the comprehensive point of view of the paper, the guidance law of the missile determines
whether the missile can accurately hit the target. In the process of attacking the target, the missile needs
to comprehensively consider the position state of the missile attack target, and select the appropriate
proportional coefficient K to determine whether the missile trajectory can be carried out. Reasonable
control, etc., and the process of missile attacking the target is a very complex process, and there are
many factors that affect it. In the process of modeling in this paper, we have simplified the model to
facilitate the analysis of various problems. Many factors We have not considered it, and this is also the
direction of our next research.</p>
    </sec>
    <sec id="sec-9">
      <title>5. References</title>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <given-names>Gao</given-names>
            <surname>Shang</surname>
          </string-name>
          .
          <article-title>Simulation of ideal ballistics with proportional guidance[J]</article-title>
          .
          <source>Computer Engineering and Design</source>
          ,
          <year>2003</year>
          (
          <volume>08</volume>
          ):
          <fpage>66</fpage>
          -
          <lpage>68</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2] BWicaksono,AHalim,Kurdianto.
          <article-title>Performanc e evaluation of proportional navigation methode during target pursuit with disturbance[J]</article-title>
          .
          <source>Journal of Physics: Conference Series</source>
          ,
          <year>2018</year>
          ,
          <volume>1130</volume>
          (
          <issue>1</issue>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <given-names>Feng</given-names>
            <surname>Tyan</surname>
          </string-name>
          .
          <source>Analysis of General Ideal Proportional Navigation Guidance Laws[J]</source>
          .
          <source>Asian Journal of Control</source>
          ,
          <year>2016</year>
          ,
          <volume>18</volume>
          (
          <issue>3</issue>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [4]
          <string-name>
            <given-names>Zhang</given-names>
            <surname>Zhongnan</surname>
          </string-name>
          , TongYoutang. Simulation research
          <article-title>on guided ballistics by proportional guidance method</article-title>
          [J].
          <source>Tactical Missile Technology</source>
          .
          <year>2005</year>
          (
          <volume>02</volume>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [5]
          <string-name>
            <given-names>Yunes</given-names>
            <surname>Sh. ALQUDSI</surname>
          </string-name>
          ,
          <string-name>
            <surname>Gamal</surname>
            <given-names>M.</given-names>
          </string-name>
          <article-title>EL- BAYOUMI. A Qualitative Comparison between the Proportional Navigation and Differential Geometry GuidanceAlgorithms</article-title>
          [J].
          <source>INCASBULLETIN</source>
          ,
          <year>2018</year>
          ,
          <volume>10</volume>
          (
          <issue>2</issue>
          ).
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>