 Open Access
 Authors : B. Naveen Naidu , K. Aravind Kumar , K. Sai Priyanka , D. Govardhan
 Paper ID : IJERTV10IS060036
 Volume & Issue : Volume 10, Issue 06 (June 2021)
 Published (First Online): 08062021
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Design of Lateral Acceleration Control Autopilots for Missile Systems
B. Naveen Naidu1, K. Aravind Kumar2, K. Sai Priyanka3, D. Govardhan4
UG student1, UG student2, Assistant professor3, Head of the department4
Dept. of Aeronautical Engineering, Institute of Aeronautical Engineering, Telangana, India500055
Abstract: The equation defining the threeloop autopilot design and its airframe dynamics of a guided missile using a stable closed loop control system. The airframe dynamics of a guided missile are nonlinear, however the nonlinearities are in structured form (i.e., quadratic form), and thus novel approach has been considered in expressing the non linear dynamics in state space form. The derivation of airframe dynamics for a rear controlled autopilot missile has been derived using the aerodynamic transfer function and characteristic equations. The autopilot design is having been directly evaluated according to the gains acquired by the unit step responses. The proposed technique of integrated guidance and control system is implemented with the assistance of linearized state space model which can be further extracted using the linear techniques of modern control theory. Finally, the numerical simulation has been done using MATLAB and demonstrated the results for the effectiveness and feasibility of the proposed technique.
Keyword: Guided missiles , threeloop autopilots , step response , frequency responses , MATLAB

INTRODUCTION TO MISSILE SYSTEMS
Introduction
Any object thrown at a target with the aim of hitting it's a missile. Thus, a stone thrown at a bird may be a missile. The bird, by using its power of reasoning may evade the missile (the stone) by moving either to the Left, right, top or bottom with reference to the flight path (trajectory) of the missile. Thus, the missile, during this case, has been ineffective in its objective of hitting the bird (the target). Now, if the stone is imparted with some intelligence and quick response to maneuver with reference to the bird, to beat aiming errors and therefore the bird's evasive actions and hit it accurately, the stone now becomes a missile . Thus, guided missile are selfpropelled, unmanned space or air vehicle carrying an explosive warhead whose path are often adjusted during flight, either by automatic selfcontained controls or remote human control and are powered either by rocket engines or by reaction propulsion . Missiles could also be aerodynamic, i.e., controlled by aerodynamic surfaces and following a straightline trajectory to the target, or ballistic, i.e., powered during flight and following a parabolic trajectory. The trajectory of the missile could be predetermined or continuously varying, based on whether the target is stationary or moving. It is also probable that the missile deviates from its path due to rigging errors, disturbances in the atmosphere, or asymmetrical loading of control surfaces in supersonic flight. In both the situations, a guided missile should be able to maneuver back to its desired path in order
to hit the target. One of the tasks of the missile guidance system is to measure the deviations (in both vertical and in horizontal planes) from missiles desired path. This deviation or error is fed to the control system, which in response to these error signals maneuver the missile quickly and efficiently until these errors are reduced to zero, with the help of various controls available onboard.

Guidance and Control System
The guidance and control system are a key element that allows the missile to meet its system performance requirements. The objective of the guidance and control system is to force the missile to achieve the steering commands developed by the guidance system. The types of steering commands vary depending on the phase of flight and the type of interceptor. For example, in the boost phase the guidance and control system may be designed to force the missile to track a desired flightpath angle or attitude. In the midcourse and terminal phases, the system may be designed to track acceleration commands to affect an intercept of the target.
Figure 1 : Guidance and Control System Loop
The four basic elements of the flight control system are shown in the figure 1. The Inertial Measuring Unit senses the inertial motion of the missile. Its outputs and the inputs from the guidance law are combined in the autopilot to form a command input to the control effectors, such as the commanded deflection angle to an aerodynamic control surface. The actuator turns the autopilot command into the physical motion of the control effectors, which in turn influences the airframe dynamics to track the guidance command.

Autopilot
An autopilot is a closed loop system and it is a minor loop inside the main guidance loop; not all missile systems require an autopilot. A missile will maneuver updown or leftright in an apparently satisfactory manner if a control surface is moved or the direction of thrust altered. If the missile carries accelerometers and/or gyros to provide additional feedback into the missile servos to modify the missile motion then the missile control system consisting of servos, control surfaces or thrust vector elements, the airframe, and feedback instruments plus control electronics is usually called an autopilot. Broadly speaking autopilots either control the motion in the pitch and yaw planes, in which case they are called lateral autopilots, or they control the motion about the fore and aft axis in which case they are called roll autopilots. This contrasts with the usual definition of aircraft autopilots; those designed to control the motion in the pitch plane are called lateral autopilots. For instance, an aircraft autopilot designed to keep the heading constant would be called lateral autopilot. For symmetrical cruciform missile however pitch and yaw autopilots are often identical; one injects a g bias in the vertical plane to offset the effect of gravity but this does not affect the design of the autopilot. Roll autopilots serve quite a different purpose and will be considered separately.
Figure 2 : Types of Autopilots
The autopilots are mainly classified into three types as shown in the Fig: 2 the special autopilot is used in some of special design cases. But the roll autopilot and Lateral autopilot are mainly used. Roll autopilot is meant to maintain particular roll attitude or to keep the roll rate zero to prevent coupling effects between pitch, yaw and roll planes. Lateral Attitude Autopilot is mainly meant for the missile applications where specific attitude is required to be maintained. It is being used in ballistic missiles or strategic missiles. Lateral Acceleration autopilot is mainly meant for maintaining the trajectory. This is mostly used in tactical missiles.

Rigid Body : The body does not go any changes in size and shape.

Aerodynamic symmetry in Roll

Mass : A constant mass is assumed , dm/dt=0

Body axes coordinate frame is used to write the missile equations.

A spherical earth rotating at constant angular velocity is assumed .

Vehicle aerodynamics is nonlinear.

The undisturbed atmosphere rotates with the Earth.

The winds are defined with respect to the Earth.

An inversesquare gravitational law is used for the spherical Earth model.

The gradients of the lowfrequency winds are small enough to be neglected.

Eulers Equation of Motion for a Rigid Body
The rigid bodyequations are obtained from Newtons second law which states that the summation of all external forces acting on a body is equal to time rate of change of momentum of body, and the summation of external moments acting on the body is equal to time rate of change of moment of momentum (angular momentum). There are six equations of motion for a body with six degrees of freedom: three force equations and three moment equations. The equations are somewhat simpler if the mass is constant and radiiofgyration changes slowly. The translation and rotation of a rigid body may be expressed mathematically by the following equations:
Translation: F = ma (3.1) Rotational: M = (3.2)
These motions are shown in Figure, the translations being (u, v, w) and the rotations (P, Q, R).


METHODOLOGY
According to Newtons laws we will consider the missile and derive the rigid body equations and the following assumptions are considered :
Figure 3 : Representation of the missiles six degrees of freedom
Mathematically, we can write the missile vector velocity,
VM, in terms of the components as
VM =ui + vj + wk (3.3)
where (i, j, k) are the unit vectors along the respective missile body axes. The magnitude of the missile velocity is given by
VM=VM = (u2 +v2 +w2)1/2 (3.4)
In a similar manner, the missiles angular velocity vector can be broken up into the components P,Q, and R about the (Xb, Yb , Zb) axes, respectively, as follows:
=Pi + Qj + Rk (3.5)
where P is the roll rate, Q is the pitch rate, and R is the yaw rate.
Fx = m ( + wQ vR), (3.6)
Fy = m ( + uR wP ), (3.7)
Fz = m ( + vP uQ ). (3.8)
L= Ix + (Iz Iy)QR ( + PQ)Ixz (3.9)
M = Iy + (Ix Iz)PR + (P 2 R2)Ixz (3.10)
N = Iz + (Iy Ix)PQ ( QR)Ixz (3.11)
Thus the set of Equations (3.6), (3.7), (3.8), (3.9), (3.10), (3.11) represents the complete 6Degrees of Freedom missile equations of motion .
Aerodynamic Transfer functions
To design a missile control system, the linearity of hardware (i.e. the electronics, fin servos, instruments and equations of motion) is generally assumed. Aerodynamic transfer functions are obtained from aerodynamics derivatives defined above. With pitch, roll and yaw dynamics under consideration, aerodynamic derivatives are force derivative if they are used in force equation and moment derivatives if used in moment equation.
The reference axis system standardized in the guided weapons industry is centered on the center of gravity and fixed in the body as follows:
X axis, called the roll axis, forwards, along the axis of symmetry if one exists, but in any case, in the plane of symmetry.
Y axis, called the pitch axis, outwards and to the right if viewing the missile from behind.
Z axis, called the yaw axis, downwards in the plane of symmetry to form a righthanded orthogonal system with the
other two .
Table 1 defines the forces and moments acting on the linear and angular velocities, and the moments of inertia; these quantities are shown in fig 4 .The moments of inertia about O are defined as:
A= m (y2+z2) B= m (z2+x2) C= m (x2+y2)
The products of inertia are defined as: D= m yz
E= m xz F= m xy
The yaw plane is the Oxy plane and the pitch plane is the Oxz plane. The following angles are defined:
: Incidence in the pitch plane.
: Incidence in the yaw plane.
: Incidence plane angle.
: Total incidence, such that:
tan = tan cos and tan = tan sin
The reason why U, the missile velocity along the x axis is denoted by a capital letter is to emphasize that it is a large positive quantity changing at most only a few per cent per second. The angular rates and components of velocity along the pitch and yaw axes however, tend to be much smaller quantities which can be positive or negative and can have much larger rates of change.
Figure 4 : Earth fixed and Body axes coordinate systems

Notations and Conventions
The motion of a missile is mostly described using two sets of coordinate systems such as earthfixed coordinate system and body axes coordinate system. The Earthfixed coordinate system or inertial coordinate system is fixed to earth and may be used to define missiles position and attitude in a threedimensional space. While the body axes coordinate system is fixed to the body of missile and is centered at the centerofgravity (CG) of the missile. Figure 4 illustrates these two axes systems.
Where x, y and zaxes are called Roll axis, Pitch axis and Yaw axis respectively. The forces acting on the missile airframe such as aerodynamic, thrust and gravitational forces are resolved along these axes. The forces and moments acting on the missile, the linear and angular velocities and other motion variables are defined in Table 1. Figure 5 shows these quantities.
Table 1 : Notations
Figure 5 : Convention
These are six equations of motion for a body with six degrees of freedom: three force equations and three moment equations. The equations are somewhat simpler if the mass is constant and the radii of gyration change slowly, in which case the standard Euler equations can be used. During boost the mass rate of change can be relatively fast but the radii of gyration usually change slowly. In the great majority of design calculations, the standard equations are used. They are:
If an autopilot with acceleration on position due to the higher frequency. If an autopilot with acceleration in position due to the higher frequency. If an autopilot with acceleration feedback is used the effects will be reduced. In any case the guidance loop has a certain stiffness which will be reduced.
In any case the guidance loop has a certain stiffness which will tend to reduce this effect; the missile does not fly open loop.
Now consider the moment equations. Ideally these should read
A = L, B = M, C = N (4.7)

., moments about a given axis produce angular accelerations about that axis. All other terms in these equations are crosscoupling terms and are undesirable from the point of view of system accuracy. We note that three out of four of the crosscoupling terms in each equation disappear if the products of inertia are zero of there are two axes of symmetry, and two will be zero and the missile is reasonably symmetrical about another axis. With two planes of symmetry and a small roll rate therefore these equations reduce to.
m ( + qw – rv) = X
(4.8)
m ( + ru) = Y
(4.9)
m ( – qU) = Z
(4.10)
A = L
(4.11)
B = M
(4.12)
C = N
(4.13)

The justification for neglecting the terms pq, pr, pv, pw is that the terms q, r, v and w are not large terms and if p is small then their products can be neglected.

Equation 4.11 Shows that there is zero coupling between the pitch and roll and yaw and roll motions if there are two axes of symmetry (B=C), and unless the missile is very unsymmetrical the cross coupling should be weak.

Nevertheless, even with a missile with two axes of symmetry, the principal moments of inertia about the roll axis (A). Moderate roll rates will certainly affect the accuracy of the system.



Derivation of Aerodynamic Transfer Functions
We will now write down the control equations for missile, omitting the force equations along the xaxis as this neither affects the roll, pitch or yaw motion. We will also omit the gravity term and consider forces and moments which are purely of aerodynamics origin.
The forces equation in the direction of Yaxi is given by
(4.14)
But for small angles of tan is equal to ,
Thus, we get,
But Zw = Z / U, as U is the velocity along X axis which is very large, Zw tends to zero.
Similarly the body rate is equal to
(4.15)
(4.16)
The viscous effect of the air is assumed negligible; Thus, the body rate terms mq, Zq are negligible.
If we consider a Simple pendulum with initial deflection , then for the stable condition we get the frequency as shown below
Figure 6 : Signal Flow graph
By using the above signal flow graph and by using Masons gain Formula:
Transfer Function () by using the Masons gain formula
()
is given by
(4.17)
Taking the Inverse Laplace transform we get,
Where it vibrates with a frequency of
If we take the static stability condition then mw will be a negative term.
Where, wb2 = Umw
If we take the Center of gravity as the pinned joint and the missile as a simple pendulum as shown above, then for a small disturbance the missile will undergo oscillations with the frequency wb which is known as the weathercock frequency.
In the similar way we can derive the transfer function for the acceleration and the deflection of the control surface.
By doing some manipulations to the two transfer functions we get,
(4.18)
Lateral Autopilot. In contrast to missile autopilot, an aircraft autopilot designed to control the motion in the pitch plane and yaw plane are called longitudinal autopilot and lateral autopilot. Where the lateral autopilot is designed to keep the heading constant. However, in case of a symmetrical cruciform missile, pitch and yaw autopilot are often identical while one injects a g (acceleration due to gravity) bias in the vertical plane to counter the effect of gravity. But this does not affect the design of the autopilot.

Design of Lateral Acceleration Autopilot
We use the transfer functions in 4.18 , 4.19 for the design of the three loop Lateral Acceleration Autopilot using pole placement technique. The main objectives for designing the lateral acceleration (LATAX) autopilot:

The maintenance of nearconstant steady state aerodynamics gain.

To increase the weathercock frequency.

To increase the weathercock damping.

To reduce crosscoupling pitch and yaw motion .

To assist in gathering.
3.1.1 A Lateral Autopilot with an actuator:
Where,
In the similar manner we get,
Where,
(4.19)
Figure 7 : Signal flow graph for the threeloop autopilot
Equations (4.18), (4.19) are the required transfer functions. Thus, the transfer functions for the Acceleration, Body rates with respect to deflection of the control surface are derived.
3. AUTOPILOT DESIGN
Missile control system consisting of servos, control surfaces, the airframe and the feedback instruments plus the control electronics, all working together to automatically adjust the orientation of the missile in space can be termed as an autopilot. In general, an autopilot either control the motion about the fore and aft axis (called Roll Autopilot) or they control the motion in the pitch and yaw planes (called
Let the transfer function of the actuator be Ta(s) which is having a second order characteristic equation.
Characteristic equation for the figure 7 is as follows
(5.1)
As the above characteristic equation is of fifth order, we cannot solve the above equation, so we need a characteristic equation which is having a second order natural frequency and the remaining three variables are placed according to our design conditions.
(5.2)
We are taking the natural frequency of the actuator value is three times the value of the weathercock frequency. The desired frequency of the desired autopilot is again three times that of the natural frequency of the actuator natural frequency. Thus, we need to place the three variables in between the natural frequency of the autopilot and the frequency of the actuator.
By equating the equations (5.1) and (5.2), we get, Let a = wn.
\
Thus the above equations are used to solve for the unknown terms (kq, ka, wi, a, b, c) in the signal flow graph.
MATLAB :
A MATLAB code is developed in order to solve the aerodynamic transfer functions we derived and plot the step response of the three loop LATAX control autopilot and frequency response (BODE) plots of the following:

Analysis with closed loop

Stability Margins with innermost loop open

Stability Margins with inner/intermediate loop open

Stability Margins with outermost loop open
The program codes to maintain the phase angle between – 180Â° & +180Â° and the analysis of the design without errors are both inserted in the main code.
The flight parameters for some assumed missiles are taken and are solved using MATLAB.
The MATLAB code and the results are found in the Appendix.



. RESULT AND DISCUSSION
Analysis of LATAX control autopilot :
The relations obtained from our three loop LATAX autopilot design with the actuator Ta(s) are solved using MATLAB .
The required control parameters for Lateral Autopilot, with their first estimate values are as listed:
also, with estimated values (on the basis of Unitstep response) of various gains, natural frequencies and damping
ratios, which can be optimized are as follows:
The step response and frequency response plots obtained from MATLAB program code are :
Figure 8 : Step Response of LATAX Control Autopilot
This gives a unit step response with the following parameters:
Rise Time = 0.1803 sec
Settling Time = 0.0848 sec (within 2%) Frequency Response Plots :
Figure 9 : Frequency Response Plot of Closed Loop Autopilot
In the frequency response of the closed loop, the change in magnitude is observed at frequency value of 37.57 rad/s
Figure 10 : Frequency Response Plot with Inner Loop Open
For the bode plot with the inner loop open , in bandwidth of
0.026rad/s and 0.65 rad/s the gain margins are observed . The phase margin frequencies are 30.05 rad/s and 67.25 rad/s with the phase angles at 139.88Â° and 179.899Â° respectively.
Figure 11 : Frequency Response Plot with Inner most Loop Open
With the inner most loop open , the magnitude is 0.05012 and 7.688 at gain crossover frequencies at 61.07 rad/s and
119.1 rad/s. The phase angles are 152.7Â° and 179.99Â° with the frequencies of 61.33 rad/s and 119.749 rad/s respectively.
Figure 12 : Frequency Response Plot with Outer Loop Open
In the frequency response with outer loop open , the gain margin frequencies are 12.06 rad/s and 45.163rad/s along with the phase angles at 104.602Â° and 179.995Â° respectively.
P1 =791.8308
P2 = 1.0460e+06
Q1 =1.0460e+06
Q2 = 2.9053e+06
R1 = 1.3035e+06
R2 = 3.4645e+07
P1 =791.8308
P2 = 1.0460e+06
Q1 =1.0460e+06
Q2 = 2.9053e+06
R1 = 1.3035e+06
R2 = 3.4645e+07
The unknown parameters of the derived characteristic equations have been acquired by using MATLAB and are as follows :
MATLAB CODE :
clear all; close all;

APPENDIX

Table 2 : Characteristic Equation Solution Values
5. CONCLUSION
The accelerometer location relative to the CG greatly affects the measured acceleration and may be wont to shape the zeros within the acceleration transfer function. The missile autopilot design methodology should consider the unfavorable aspect of onboard hardware feature in high frequency reion and employs openloop crossover frequency to constrain high frequency dynamics. The pole placement technique used has produced with the proper values of the control parameters and the desired step response was achieved within short period of time. The actuator portion plays an important role in the autopilot by sending continuous signals to the control system loop and causing the necessary deviations of the missiles in its aerial path . It allows the missile to guide towards the target with no room for errors in trajectory , path or in any sort that may cause failures. These responses obtained in our design of LATAX says that these three loop autopilots are very much reliable in the missile systems in both ways of attacking and counterattack modes of the military.
FUTURE SCOPE
Further improvements of these three loop autopilots, in which the application of RSS technology on missile design can introduce much benefit to missile maneuverability and full flight performance. However, the missile with lower static stability requires a complicated control strategy. Both the twoloop and threeloop topologies can stabilize a static unstable missile and achieve desired performance requirement. The exhibited performance of the classical three loop topology is analogous to an inertial loop, and can be modified for the wants of next generation airtoair agile missile. The Iron Dome is an example of successful missile defense system with 90 percent accuracy , but it is for short range. It requires a lot of efforts to develop a longrange autopilot missile system which is accurate. Much more research and analysis will help in building a missile control system design of such accuracy and efficacy.
%%% plant parameters zq=0.0; mq=0.0;
mdelta =53.0;
zdelta = 20.06;
mw = 0.2948;
zw = 2.889;
U =500;%%%% Missile speed Ta=mdelta/(zdelta*mwzw*mdelta) wbsq = U*mw
kb=(zdelta*mwzw*mdelta)/wbsq sigmasq=zdelta/(U*(zdelta*mwzw*mdelta)) ls=0.0;%xcgxsense
%c=1;
%for ki=0.0:0.1:30
%wa=2*pi*wa;
%%%%% desired system % specified system
zn=0.9 wn=sqrt(abs(wbsq))*3.0%1.2*sqrt(abs(wbsq)); wnsq=wn*wn
znwn=zn*wn
sys_des = tf([wnsq],[1 2.0*znwn wnsq]);%%% Transfer function of desired system
%%%% actuator part
%wa=2*pi*10;
%wa=3.5*wn; zawa=3*znwn za=0.7 wa=zawa/0.7 wasq=wa*wa
%%%%%%%%%%%%%%%
alp=(1.0znwn/zawa); a=alp*(2*zawa)*0.5
%%%%%%%%%%%%%%%
apbpc=2*(zawaznwn) bpc=2*(zawaznwn)a
bc=wasq+wbsqwnsq2*znwn*apbpca*bpc abpbcpca=a*bpc+bc
abc=a*bc
%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
P1=kb*wasq*wbsq*(U*sigmasqls*Ta)/U
Q1=kb*wasq*wbsq*Ta
R1=wnsq*apbpc+2*znwn*abpbcpca+abc*(1+wnsq*(U*sig masqls*Ta)/U)2*zawa*wbsq
P2=kb*wasq*wbsq*(Tals/U) Q2=kb*wasq*wbsq
R2=wasq*wbsq+wnsq*abpbcpca+2*znwn*abc (abc*wnsq*ls)/U
%%%%%%%%%%%%%%%%%%%%%% kqf=(R2(P2/P1)*R1)/(Q2(P2/P1)*Q1)
kqfwi=(R1kqf*Q1)/P1 wi=kqfwi/kqf
kawi=(kqfwi/U)(abc*wnsq)/(U*kb*wasq*wbsq) ka=kawi/wi
%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%
CE=[1,…
2*zawa,…
wbsq+wasq,…
2*zawa*wbsq+kawi*U*kb*wasq*wbsq*sigmasq+kqf*kb* wasq*wbsq*Ta,…
wasq*wbsq+kqf*kb*wasq*wbsq+kqfwi*kb*wasq*wbsq*T a,…
kb*wasq*wbsq*(kqfwiU*kawi)] BCE=[1,…
apbpc+2*znwn,… wnsq+2*znwn*apbpc+abpbcpca,… wnsq*apbpc+2*znwn*abpbcpca+abc,… wnsq*abpbcpca+2*znwn*abc,… abc*wnsq]
%%%%%%%%%%%%%%%%%%%%
rootsCE=roots(CE) rootsBCE=roots(BCE)
analysis;%%%% to carry out time and frequency domain analysis
%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%
%lroots(c,:)=rootsBCE;
%c=c+1;
%end
%%
%
% for x = 1:10
%
% for x = 1:10
%
% for x = 1:10
%
% for x = 1:10
% disp(x)
% end
%
% disp(x)
% end
%
% disp(x)
% end
%
% disp(x)
% end
%
% %ki=0.025:0.2:30;
% %zn=0.5:0.01:1.5
% % wa=0.5:0.1:30;
% ki=0.0:0.1:30;
% % % %
% % % %
% figure;
%
% plot(ki,lroots);gzh;
REFERENCES

Optimization of three loop autopilot gain under cross over frequency constraint by LongJun Qian*and YuMing BO*, School of automation, University of Science and Technology , Ninjang 210094, China . (September 2015)

Design and Comparison of TwoLoop with PI and ThreeLoop Autopilot for Static Unstable Missile by Emad Mahrous Mohamed*, Lin Yan School of Automation Science and Electrical Engineering, Bei hang University BUAA), Beijing 100191, China. (2015)

Threeloop Autopilot Design and Simulation by CHU Hairong*
, ZHANG Yue* , Changchun* ; Institute of Optics, Fine Mechanics and Physics, Chinese Academy of Sciences Changchun Institute of Optics, Fine Mechanics and Physics, Chinese Academy of Sciences Changchun, Jilin, China Proceedings of 2012 IEEE International Conference on Mechatronics and Automation August 5 – 8, Chengdu, China.

Modified Design of Three Loop Lateral Missile Autopilot based on LQR and Reduced Order Observer by Parijat Bhowmik , Prof. Gourhari Das , Dept. of Electrical Eng. (Control Systems), Jadavpur University, Kolkata, India. International Journal of Engineering Research and Development eISSN ,Volume 2, Issue 6 (August 2012)

Analysis and improvement of missile threeloop autopilots, Lin Defu; Fan Junfang; Qi Zaikang; Mou Yu published in Journal of Systems Engineering and Electronics ( Volume: 20, Issue: 4, Aug. 2009)

Robust Stability Analysis of Adaptive Missile Autopilots by Kevin A. Wise in The Boeing Company P.O. Box 516, St. Louis, MO 63166 AIAA Guidance, Navigation and Control Conference and Exhibit 18 – 21 August 2008, Honolulu, Hawaii.

Integrated Design of Agile Missile Guidance and Control Systems by P. K. Menon* and E. J. Ohl Meyer ; Paper Presented at the 1999 Mediterranean Control Conference, June 2830, Haifa, Israel.