Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Calculation and Identification of The Aerodynamic Parameters

Download as pdf or txt
Download as pdf or txt
You are on page 1of 18

sensors

Article
Calculation and Identification of the Aerodynamic
Parameters for Small-Scaled Fixed-Wing UAVs
Jieliang Shen 1 , Yan Su 1, *, Qing Liang 2 and Xinhua Zhu 1, *
1 School of Mechanical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China;
perfect.sjlchg2008@163.com
2 School of Computer Technologies and Control, ITMO University, St. Petersburg 197101, Russia;
liangqing1688@gmail.com
* Correspondence: suyan@njust.edu.cn (Y.S.); zhuxinhua@njust.edu.cn (X.Z.)

Received: 9 October 2017; Accepted: 11 January 2018; Published: 13 January 2018

Abstract: The establishment of the Aircraft Dynamic Model (ADM) constitutes the prerequisite for
the design of the navigation and control system, but the aerodynamic parameters in the model could
not be readily obtained especially for small-scaled fixed-wing UAVs. In this paper, the procedure of
computing the aerodynamic parameters is developed. All the longitudinal and lateral aerodynamic
derivatives are firstly calculated through semi-empirical method based on the aerodynamics, rather
than the wind tunnel tests or fluid dynamics software analysis. Secondly, the residuals of each
derivative are proposed to be identified or estimated further via Extended Kalman Filter (EKF),
with the observations of the attitude and velocity from the airborne integrated navigation system.
Meanwhile, the observability of the targeted parameters is analyzed and strengthened through
multiple maneuvers. Based on a small-scaled fixed-wing aircraft driven by propeller, the airborne
sensors are chosen and the model of the actuators are constructed. Then, real flight tests are
implemented to verify the calculation and identification process. Test results tell the rationality
of the semi-empirical method and show the improvement of accuracy of ADM after the compensation
of the parameters.

Keywords: aerodynamic parameters; semi-empirical aerodynamic coefficient modeling; parameters


identification; EKF; real flight tests

1. Introduction
In contrast to the large-scaled high-altitude long-term Unmanned Aerial Vehicles (UAVs) for
remote surveillance and combat, the small-scaled fixed-wing UAVs, like “Raven” RQ-11 and Sand
Hawk of the US, or “Rainbow-802” of China, have unique superiority in civilian and military fields.
With the light weight, small volume and simple take-off, small UAVs could execute tasks like the
close-range monitoring in a covert and flexible way [1]. The establishment of Aircraft Dynamic Model
(ADM) should be accomplished first before any task designing. It contains the computation of aircrafts’
structural parameters, position of the center of gravity, moment of inertia, aerodynamic parameters
and the modeling of actuators.
For small-scaled fixed-wing UAVs, the structural parameters, like the wing span l, mean
aerodynamic chord c and the reference area S, could be measured precisely. The thrust, driven by the
propeller, can be determined by the diameter of the propeller D, rotation rate n P and airspeed Va [2].
However, the aerodynamic parameters, covering stability derivatives and the control derivatives, are
relatively hard to get, which is crucial in calculating the force and moment. Stability derivatives involve
partial derivatives with respect to states and control derivatives involve partial derivatives with respect
to control inputs. The stability derivatives could be further divided into static stability derivatives for
derivatives associated with air-relative velocity quantities and the dynamic stability derivatives for

Sensors 2018, 18, 206; doi:10.3390/s18010206 www.mdpi.com/journal/sensors


Sensors 2018, 18, 206 2 of 18

derivatives associated with angular rates and unsteady aerodynamics [3]. Generally, the aerodynamic
parameters are analyzed through wind tunnel tests or CFD software, which are not appropriate for
small ones. Wind tunnel test is high-cost and time-consuming, most of all, the results are not accurate
in the flight condition of low Reynolds number. And CFD software analysis usually has heavy work
on the 3D model meshing, and the gap effect as well as the frictional drag is often neglected [4,5].
DATCOM, introduced by the US Air Force, is another choice for the parameter calculation. On the
basis of the huge flight database, all the aerodynamic coefficients could be obtained when the required
parameters are imported. But DATCOM does have drawbacks. The input parameters are quite
complex and the analysis will somehow be affected by the computational interval of the angle-of-attack
α, resulting in limited accuracy for small UAVs [6,7].
Considering the special aerodynamic configuration and flight condition of the small UAVs, the
semi-empirical method based on the aerodynamic theory is proposed here, by means of which the
longitudinal and lateral derivatives could be calculated step by step according to the semi-empirical
formulas or diagrams. Semi-empirical method was also mentioned in Arifianto’s research [8].
The reason why it is called ”Semi-empirical” is that the formulas and diagrams used are actually
the combination of the fundamental aerodynamics and empirical summary from abundant flight
tests. For instance, the diagram for calculate the lift curve slope of wing CW α is fitted with the flight

data according to the theoretical relationship with the geometrical shape of the wing and the Mach
number. Monographs [9–14] introduce the basic aerodynamic theory of different aircrafts flying
under different circumstances, and describe the detailed procedure of calculating all the aerodynamic
derivatives. Specifically, the aerodynamic analysis of small-scaled aircrafts is often simplified for
quite small Mach number, plain aerodynamic structure and small angle-of-attack assumption. So,
this kind semi-empirical method could rationally calculate the initial parameters in a small amount
of computation with low cost and it will be adopted in the paper. Although the semi-empirical
analysis is an effectively feasible method, the accuracy of the parameters is sometimes suspicious
due to the structure simplification and the changing of the flight condition. After the theoretical
calculation, the aerodynamic parameters are nearly all constants, which is understandable in the small
Mach number condition and small angle-of-attack assumption. It means that the method owns poor
real-time performance. Therefore, the identification process is carried out to further compensate the
error of aerodynamic parameters with the real-time flight data.
Monographs [3,15] synthesize the basic issues of system identification for fixed-wing aircrafts or
rotorcrafts, including aircraft’s mathematical model, estimation theories, classification of identification
methods as well as engineering practices, like experiment design, data compatibility check and data
analysis. Generally, there exists many approaches to identify or estimate the error of the aerodynamic
parameters, and it could be classified into off-line and on-line way. Off-line way consists of the formula
error method and the output error method, the former is based on the Least Square (LS) while the
latter based on the Maximum Likelihood Method(MLE). And off-line way could be operated either in
the the time domain or the frequency domain [16,17]. Neural network is also chosen as the structure to
describe or identify the dynamic characteristics of aircraft [18,19]. All the methods above own a large
amount of computation and call for the whole accurate measurements corresponding to the system
states. Morelli [20] focuses on the aerodynamic model identification of the combat aircraft like F-16
using MLE in the frequency domain, and the measurements used are the 6 outputs of the IMU and the
aerodynamic states like Va , α and β, which requires high-precision IMUs and expensive sensors for
α, β detection. But it is quite an expectation for small-loaded aircrafts. Dorobantu’s research [21,22]
exhibit the system identification for small, low-cost, fixed-wing UAVs, which usually carry lightweight
devices, like the MEMS IMUs, but the accuracy of the measurement is usually hard to guarantee.
On-line identification, or the real-time estimation, could handle the problems above. Based on
the linear or non-linear Kalman Filter, errors of the aerodynamic parameters are expanded into the
motion states, and then estimated to some extent in a real time. Garcia [23] introduces that EKF
method breaks the limitation of the linearity of the system and the estimation results could be variant.
Sensors 2018, 18, 206 3 of 18

But, the inaccuracy of the statistical characteristics of noise and the observability issue affect the
estimating results. According to [24], maneuvering motion, like the changing of the acceleration
and angular velocity, will be able to increase the observability of system states. Hence, the specific
maneuvers excited by the changing of the inputs of the actuators definitely improve the observability
of each expanded states [25]. After the categorization of current identification methods for different
kinds of aircrafts equipped with different sensors, Hoffer [26] proposed that the Square-root UKF
would be a good choice for the on-line identification. Girish [27] employs both EKF and UKF
to estimate the aerodynamic parameters with the flight data, and it is pointed out that although
UKF owns a theoretically better precision, EKF shows the same performance with less computation.
Therefore, EKF is utilized in the paper to estimate and compensate the error of all the longitudinal and
lateral aerodynamic derivatives. And the observations are the relatively accurate airborne integrated
INS/GPS attitude and velocity. Meanwhile, maneuvering flight path is designed for the promotion
of observability. The identified parameters could be used to calculate the navigation results with the
dynamic and kinematic model of aircraft, or the motion model, then aiding the low-grade INS [28].
So the precision of identification results could be evaluated by comparing the navigation results
with that of INS/GPS. It is quite different for the evaluation of the identification of second-order
transfer function.
The outline of this paper is as follows. In Section 2, the aerodynamic characteristics of the
small-scaled fixed-wing UAVs are summarized and all the parameters are calculated semi-empirically.
In Section 3, we present the process of the identification of the aerodynamic parameters in details,
and also the observability analysis is discussed. Section 4 selects the airborne sensors and the data
acquisition board, then describes the modeling of the propeller and the control surfaces as a test
preparation. In Section 5, the results of the flight tests are demonstrated as a verification. At last, the
paper is concluded in Section 6.

2. Semi-Empirical Calculation
Aerodynamic parameters reside in the aerodynamic forces and moments acting on the aircraft,
containing the lift Z, drag X and pitch moment My in the longitudinal channel and the side force
Y, roll moment Mx and yaw moment Mz in the lateral channel, shown in Figure 1. It is a typical
fixed-wing aircraft with a conventional configuration. Table 1 lists the formulas for the calculation of
the thrust driven by the propeller and all the aerodynamic forces and moments, which are the product
of the dynamic pressure Q (Q = ρVa2 /2), reference area S, reference size l or c̄ and the non-dimensional
coefficients Ci or mi , i = X, Y, Z. As seen in the formulas, Ci or mi consists of several non-dimensional
and dimensional parameters, namely the aerodynamic derivatives, which are exactly the targeted
parameters in this paper.
Aerodynamic derivatives are determined by the flight condition, aerodynamic layout and
structural characteristic. Reynolds number of small UAVs, around 200,000, belongs to low Reynolds
number. Laminar separation occurs and the aerodynamic efficiency decreases. For instance, the
relationship between the lift coefficient and α, namely CZα , would be nonlinear as the α increases.
Hysteresis effect would happen around the stall angle resulting in the abnormal change of CZα [29].
Given the condition of small Mach number Ma, the real-time changing of the aerodynamic derivatives
caused by the velocity could be neglected, and the assumption of small α or β could make the
β
derivatives, for example CZα , CY , constants [12,13]. Therefore, the nonlinear aerodynamic derivatives,
2 αq
like CZα , CZ , CZ0 (α, β), are neglected for the ignoring of large-amplitude maneuvers and high
angle-of-attack, which could be a huge simplification. And dynamic derivatives referenced to the α̇,
β̇ could be small and also neglected, just as shown in Table 1. Moreover, small-scaled UAVs usually
own a conventional configuration, namely wing-body-horizontal tail-vertical tail structure, and the
sweepback angle χ or dihedral angle ψd doesn’t exist in the wing designing. So the complexity of
the aerodynamic analysis is reduced especially for the lateral channel. And the wing of small UAV
functions predominately in the aerodynamic performance, so substituting the wing for the whole
Sensors 2018, 18, 206 4 of 18

aircraft could be acceptable for the lift or pitch moment parameter calculation. The detailed process of
calculating all the aerodynamic derivatives will be analyzed then.

Figure 1. Aerodynamic Forces & Moments.

Table 1. Formulas of the Forces & Moments.

Forces/Moments Formulas
  2 
Va Ma Va Ma
Thrust FT FT = ρn2P D4 CFT1 + CFT2 Dπn P
+ CFT3 Dπn P

Lift Z Z = QSCZ CZ = CZ0 + CZα α + CZδe δe


Drag X X = QSCX 2
CX = CX0 + CXi = CX0 + ACZ
ω̄ ωy c̄
Pitch Moment My My = QSc̄my my = my0 + mαy α + mδye δe + my y Va
β
Side Force Y Y = QSCY CY = CY β + CYδr δr
β z ωz l ω̄x ωx l
Yaw Moment Mz Mz = QSlmz mz = mz β + mδzr δr + mω̄
z 2Va + mz 2Va
β z ωz l ω̄x ωx l
Roll Moment Mx Mx = QSlm x m x = m x β + mδxr δr + mδxa δa + mω̄
x 2Va + m x 2Va

2.1. Longitudinal Channel Analysis


Total lift of an aircraft is the collective effect of the wing, body and horizontal tail. With regard
to the small-scaled fixed-wing aircraft, wing constitutes the main source of the lift. And there often
exists interaction between the wing and body. As depicted in Table 1, lift coefficient consists of CZ0 , CZα
and CZδe . Zero lift coefficient, CZ0 , is determined by the airfoil camber, and the incidence angle of the
wing ϕW and horizontal tail ϕS , which means that lift still exists when the angle-of-attack is zero. CZα ,
known as the slope of lift curve, fuses the separate cantilever-wing slope CZCWs α , the separate body
α α
slope CZBs , the separate horizontal tail slope CZHTs and all the mutual effect between each of them.
α
In details, CZCWs α
, CZHTs could be computed with the same diagram in [13,14], both of which are the
function of the structural parameters of the wing and the Mach number Ma. While the lift slope of
body could be divided into the conical or ogive head, the cylindrical body and the shrinking tail part.
CZδe concerns with the aerodynamic efficiency ηe and the size of the elevator.
Lift of the wing-body part ZWB is the sum of the lift of separate body ZBs and the lift caused by
the existence of the cantilever-wing ZCW , while the latter contains the lift of separate cantilever-wing
ZCWs , the interference lift from the body to the cantilever-wing ZB−CW and the opposite lift ZCW − B .
Interference factors [13,14], Kαα , K ϕ0 , are introduced to describe the mutual effect numerically, both
of which concerns with the shape of the wing and the diameter-to-span ratio D̄. Generally, when we
have the condition that the aspect ratio λ is large enough (λ > 5) and D̄ is small enough ( D̄ 6 0.1),
neglecting the mutual effect hardly matters. In the linear range, the lift slope of the wing-body part are
shown in (1) and (2).
Sensors 2018, 18, 206 5 of 18

SB S
α
CZWB α
= CZBs α
+ CZCWs Kαα CW (1)
S S
ϕ α SCW
CZWB = CZCWs K ϕ0 (2)
S
The airflow block factor k HT and the down wash angle ε, caused by the wing part, should be
taken into account when analyzing the lift of the horizontal tail. The block factor k HT (k HT = 0.85)
should be multiplied to the dynamic pressure while the down wash angle ε should be eliminated to
obtain the real angle-of-attack of the horizontal tail α HT . The down wash angle after the wing-body
part ε WB are nearly linear to the α and the incidence angle of the wing, so the related derivatives are
calculated in (3) and (4), in which εαW , the down wash caused by the separate wing could be calculated
in [13,14]. Thus, the lift coefficients could be deduced in Formulas (5)–(7), among which (Kαα ) HT and
(K ϕ0 ) HT are the interference factor between the body and the horizontal tail.
α
SCW CZCWs
εαWB = Kαα εαW (3)
S CZWs
α

α
ϕ SCW CZCWs εαW
ε WB = K ϕ0 (4)
S CZWs 1 − D̄
α

ϕ ϕ S HT
α
CZ0 = CZWB ϕW + CZHTs (K ϕ0 ) HT ( ϕS − ε WB ϕW )k HT (5)
S
S HT
CZα = CZWB
α α
+ CZHTs (Kαα ) HT (1 − εαWB )k HT (6)
S
S HT
CZδe = CZHT
δe
k HT (7)
S
Generally, aircraft drag contains the parasitic drag and the induced drag, neglecting the wave
drag in supersonic flight. Parasitic drag usually consists of friction drag, form drag and interference
drag, among which the friction drag are the main source of the parasitic drag in the low speed flight.
Induced drag is caused by the lift of the wing. Under the small Mach number and small angle-of-attack
assumption, the drag coefficient CX is separated into zero lift coefficient CX0 and the lift-induced drag
coefficient CXi , depicted as the parabolic formula in Table 1. CX0 mainly contains the friction part and
the due-to-lift part in the parasitic drag is neglected. So CX0 is calculated with the surface friction
coefficient c f and correction factor η [13,14]. CXi , consistent with the formation ACZ2 , is calculated
considering the wing-body part according to [13,14], in which CXi of wing are closely related to the
wing’s shape and CXi of body is proportional to α2 .
Assuming that the vector of thrust points through the center of gravity, the pitch moments are
primarily caused by the lift of different parts. Prior to calculating the pitch moment coefficients,
the action point of the lift, namely the aerodynamic center, should be analyzed, which covers the
analysis of the wing, body, wing-body and the horizontal tail. Similar to the lift, the location of the
aerodynamic center of the separate wing and horizontal tail x pW , x pHT could be obtained by the same
diagram in [13,14] according to the size and shape. Referenced to the vertex in the front, the location
of aerodynamic center is usually positive. But for the body part, x pB could be negative as the tail
generates large minus pitch moment. The location of the wing-body part ( x pWB )αα , ( x pWB ) ϕ0 could
also be computed with the factor Kαα , K ϕ0 . For the elevator, x pHTδe locates at the leading edge of the
control surface.
Along with the aerodynamic center, the aerodynamic focus is also introduced to reflect the
longitudinal stability. The location of aerodynamic focus x F , only determined by the angle-of-attack,
are related to the longitudinal static stability derivative mαy and the location of the center of gravity xG ,
depicting in (8). It is explained that when the mαy is negative (xG < x F ) , the aircraft is statically stable.
While the aircraft deviates from its equilibrium position, the pitching moment makes itself back to the
equilibrium automatically. Three static aerodynamic derivatives m Z0 , mαZ , mδZe are deduced in (9)–(11).
Sensors 2018, 18, 206 6 of 18

ω̄y
For the dynamic derivative, the damping moment coefficient my caused by the non-dimensional ω̄y
should be taken into account to strengthen the stability which could be obtained in [14].

xG − x F
mαy = CZα (8)
c

ϕ
xG − ( x pWB ) ϕ0 S HT xG − x pHT
my0 = CZWB ϕW + CZ0HT k HT (9)
c S c
SB xG − x pB xG − ( x pWB )αα S xG − x pHT
mαy = CZBs
α α
+ CZWB α
+ CZHT k HT HT (10)
S c c S c
S x G − x pHTδe
δe
mδye = CZHT k HT HT (11)
S c

2.2. Lateral Channel Analysis


The aerodynamic analysis in the lateral channel is comparatively more complicated concerning
the coupling between the roll and yaw moments. But there are some easy ways to simplify the process.
Firstly, side force Y mainly stems from the body and the vertical tail. When the aircraft rotates around
the x-axis for 90◦ , Y is just similar to the lift Z. As a result, the angle of sideslip β resembles the
angle-of-attack α and the deflection angle of the rudder δr resembles that of elevator δe . Then, the
β
aerodynamic derivatives CY , CYδr could be deduced. Similarly, the calculation of the yaw moment
derivatives could be operated after the rotation of 90◦ and at this time, the side force Y is the main cause
β
of the yaw moment. The static derivatives mz , mδzr and the dynamic derivative mω̄ z
z are computed,
β β
among which mz (mz < 0) determines the stability of yawing channel [13,14].
β
Among the roll moment derivatives, m x is an important derivative that concerns with the α, β,
the aerodynamic configuration and structural size, such as the dihedral angle ψd , sweepback angle
χ0.5 , the wing-body interaction, the wing tip characteristic and so on, depicting in (12). The aileron
plays a vital role in rolling the aircraft and the relevant control derivative mδxa could be obtained in (13),
in which ηa is the aerodynamic efficiency of the aileron. The damping roll derivative mω̄ x
x could also
be obtained just like the damping one in the yaw moment. There does exist coupling effects between
the yaw and roll moment, described by the relevant derivatives mδza , mδxr , mω̄ x ω̄z
z , m x . They are usually
small but ought to be considered according to the specific situation [13,14].

β β β β ∂2 m x ∂2 m x ∂2 m x ∂2 m x
(m x )∗ = (m x )ψd + (m x )WB + (m x )VT + [( )χ0.5 + ( )Wtip + ( )ε ]α + ( )ε δe (12)
∂α∂β ∂α∂β ∂α∂β ∂δe ∂β

β
1 mx
mδxa = CZWs
α
ηa ( · ) (13)
CZα ψd

3. Parameter Identification
After the semi-empirical analysis, the initial aerodynamic parameters are obtained, but they
are not accurate enough for further research. An identification process is necessary to correct the
parameters. The nonlinear dynamic and kinematic differential equations of V b , ω b , and Euler angles
are to be used as the system equations in the identification process, shown in (14)–(16) , in which a
non-rotating and flat earth is assumed. Considering the possible fast dynamic change of the motion
states, the error states are chosen as system states, X M = [δu, δv, δw, δp, δq, δr, δφ, δθ, δψ] T . The error
of aerodynamic parameters XC are expanded into the system state, depicted in (17). Here, XC are
assumed as constants given the condition of small Mach number and angle-of-attack, and they are
exactly the parameters to be estimated and then compensated to the initial value . The linearization
of 1st-order Taylor Series approximation of the nonlinear model is operated and the related Jacobian
Sensors 2018, 18, 206 7 of 18

matrix F is formulated in (18), in which F could be partitioned into FM , FC and the zero matrix. Adding
up the related noise W9×1 , the linearized system equation is shown in (19).

V̇ b = Cnb gn + ( FT + Cab FA )/m − ω b × V b (14)

ω̇ b = I −1 ( M − ω b × Iω b ) (15)
    
φ̇ 1 tan θ sin φ tan θ cos φ p
 θ̇  =  0 cos φ − sin φ   q  (16)
    
sin φ cos φ
ψ̇ 0 cos θ cos θ
r

X = [δu, δv, δw, δp, δq, δr, δφ, δθ, δψ,


ω̄y
δCZ0 , δCZα , δCZδe , δCX0 , δCXA , δmy0 , δmαy , δmδye , δmy , (17)
β β β ω̄z T
δCY , δCYδr , δmz , δmδzr , δmω̄ z ω̄ x δa δr ω̄ x
z , δmz , δm x , δm x , δm x , δm x , δm x ]29×1
" #
FM 9×9 FC 9×20
F= (18)
020×9 020×20

Ẋ = FX + W (19)

There are 20 aerodynamic parameters in the aircraft dynamic model according to Table 1. All the
derivatives are directly fused with the airspeed Va residing in the dynamic pressure Q. It could be
excited through acceleration or deceleration by changing the rotation rate of the propeller. Some
β β
derivatives have close relations with the aerodynamic angle α, β, like CZα , mαy , CY , m x , and the damping
ω̄y
derivatives like my , are concerned with the rotation rate of the aircraft. All of them are supposed to
be stimulated by the uniform or non-uniform changing of the attitude. And the rest derivatives are
directly related to the deflection angle of the control surfaces δe , δr , δa .
As strong nonlinearity exists among the motion states, depicted in (14)–(16), and the number of
invariant targeted derivatives are too many, it is not easy to compensate all of them correctly at the
same time. Different parameters may have different impact on the model, and the observability is also
different. Lift coefficient CZ is taken as an example to analyze the observability of each aerodynamic
derivatives, depicted in (20). With the observations of non-gravitational acceleration and the thrust, the
left hand Z is potentially calculated according to (14). Right hand of the equation is a sum formation,
and it is easier to distinguish each element with different characteristics. So the aerodynamic derivative
CZ0 , CZα and CZδe could be observed just when the angle-of-attack and deflection angle δe vary differently,
namely they have distinguishing frequency spectrum. The analysis could be generalized into the
moment coefficients, like the roll moment in (21). Also it is not easy to tell apart ωx from ωz , so mω̄ x
x ,
ω̄z δa δr
m x may not be separated totally, same theory for m x , m x . Therefore, for the similar derivatives, they
could be observed when the corresponding actuators are excited at different time.

1
Z = QSCZ = ρS(Va2 CZ0 + Va2 CZα α + Va2 CZδe δe ) (20)
2
 
1 l l
ρSl Va2 m x β + Va2 mδxa δa + Va2 mδxr δr + Va mω̄
β x ω̄z
Mx = QSlm x = x ω x + Va m x ωz (21)
2 2 2
To sum up theoretically, if all the aerodynamic parameters are to be observed, the aerodynamic
forces and moments are calculated firstly and the observations of acceleration and angular velocity
are necessary, referring to the dynamic Equations (14) and (15). But, for small-scaled fixed-wing
aircraft, it is hard for the MEMS-IMU to provide precise observations of acceleration and angular
velocity on account of the unpredictable drift or bias. As MEMS-IMU integrated with GPS is the
most common airborne navigation devices, the integrated navigation results could be adopted as
measurements for further compensation of the aerodynamic coefficients. The acceleration and angular
velocity could be deduced by the first order derivative, instead of the outputs of MEMS-IMU. So,
Sensors 2018, 18, 206 8 of 18

it is reasonable to accomplish the identification process just with the observations of velocity and
attitude. Now, the differences between the relatively accurate velocity and attitude from the INS/GPS
fusion(subscript ”SF”) and that from the integral of the aircraft motion model of Equations (14)–(16)
(subscript ”AMM”) are used as measurements, shown in (22). uSF , vSF , wSF are deduced from the
INS/GPS ground velocity VN , VE , VD through a coordinate transformation by multiplying the rotation
matrix Cnb computed with Euler angles. Then, the identification process will be accomplished by EKF.
  
u AMM − uSF

δu

 δv  
  v AMM − vSF 

 δw   w AMM − wSF 
Z= = (22)
   
φ AMM − φSF

 δφ   
   
 δθ   θ AMM − θSF 
δψ ψ AMM − ψSF
The maneuvering flight is also indispensable to excite all the motion states and the inputs of the
actuators to increase the observability of all the aerodynamic derivatives. Relevant researches [3,8,15]
introduce the specific changes of the elevator, rudder and aileron to excite the expecting maneuvering,
like dutch roll, bank-to-bank and the short period motion of the aircraft, which is an effective approach
to identify the parameters in the time domain. Here in the paper, a maneuvering flight is designed,
including acceleration and deceleration, rolling, the rise and fall when circling around.

4. Experimental Preparation
A model UAV called Extra300, shown in Figure 2, is employed to carry out the real flight test. Two
goals are about to achieve in the practical tests: (1) Rationality of the empirically calculated aerodynamic
coefficients should be validated. (2) Parameters identification is operated and the improvement
should be proved compared with the semi-empirical results. For the empirical calculation of the
aerodynamic parameters, the 3D model of Extra300 is established and the structural size is measured.
The characteristic parameters of the aircraft like mass m, position of center of gravity ( xG , yG , zG ) and
moment of inertia I are tested successively on a Mass & Gravity Center Test Board and a Moment of
Inertia Turntable. Firstly, by setting each aircraft’s axis along the central axis of the Test Board, shown
in the left part of Figure 3, m and ( xG , yG , zG ) are obtained based on the moment balance principle.
It is tested that the thrust from the propeller FT points nearly through the center of gravity. Then,
the aircraft will be installed on the turntable. Make sure that the rotation axis of the turntable points
right through the obtained gravity center, shown in the right part of Figure 3. Moment of inertia is
calculated from the rotation period tested by an optoelectronic switch. The rotation is generated by the
turntable and slowed down by a spring device. Before the whole test, a standard cylinder with known
mass (about 3 kg) and moment of inertia (about 0.00058 kg · m2 ) is used to calibrate the parameters
of the two test tables, like the elastic coefficient of the spring device. Together with the above, the
empirical aerodynamic derivatives calculated according to Section 2 are listed in Table 2.
Sensors 2018, 18, 206 9 of 18

Extra300

XSENS IMU

Under the Cockpit


In the Cabin

Data Acquisition Board

Figure 2. Flight Test Devices.

zG yG I xx I yy

xG I zz

Figure 3. Characteristic Parameter Test.

The airborne sensors and devices are shown in Figure 2. A XSENS MTi-G module is chosen, which
owns a mature integrated navigation system of INS/GPS. It also contains an Attitude and Heading
Reference System (AHRS) processor. XSENS MTi-G deploys MEMS inertial sensors in 9 axises, in which
the outputs of gyroscopes own the 1◦ /s bias stability and 0.05◦ /s noise density, the accelerometers
have the 0.02 m/s2 bias stability and 0.002 m/s2 noise density, the magnetic sensors has 0.1 m Gauss
bias stability and 1.5 m Gauss noise density. MTi-G provides not only the outputs of different inertial
sensors but also the integrated navigation results of attitude, velocity and position. The accuracy of
the integrated position and velocity is less than 2.5 m and 0.1 m/s respectively. The dynamic accuracy
of the pitch/roll angle is less than 1◦ and that of the heading angle is less than 2◦ . Thus, the complex
initial alignment process for INS or ADM is left out. The data output rate of XSENS is usually set to be
100 Hz [30].
The data of MTi-G is collected through a data acquisition board. As shown in the picture,
the collecting board, placed in the cabin, has 2 layers, both of which has the core of DSP-BF506F.
The lower board is in charge of collecting the data of all the sensors and transmitting to the upper
board. The upper board will collect the data from the lower board and that of the signal receiver of
4 channels transmitted from the controller, then output all of them to a data recorder. The data recorder
is recommended, rather than the radio transmission, to avoid the loss of data or electromagnetic
interference. Furthermore, the upper board is responsible for the switching of the flight mode. In the
tests, the aircraft is controlled remotely by a controller, namely the RC mode [31,32].
Sensors 2018, 18, 206 10 of 18

Table 2. Extra300 Parameters.

Structural
S = 0.285 m2 , l = 1.2 m, c = 0.24 m
Parameters

Propeller ρ = 1.22 kg/m3 , D = 0.33 m


Parameters CFT1 = 0.0842, CFT2 = −0.136, CFT3 = −0.928

m = 1.72 kg, xG = 0.323 m, yG ≈ 0 m, zG ≈ 0 m


 
Inertia 0.001393 0 0
Parameters I= 0 0.003543 0  kg·m2
 

0 0 0.003631

CZ CX my

CZ0 = 0.1601 CX0 = 0.01602 my0 = 0.05267


CZα = 4.4851 rad−1 CXA = 0.1347 mαy = −0.1195 rad−1
CZδe = 0.5901 rad−1 mδye = −2.8648 rad−1
ω̄
Aerodynamic my y = −6.4099 rad−1
Parameters CY mz mx

= −0.7076 rad−1 mz = −0.007517 rad−1 m x = −0.05094 rad−1


β β β
CY
CYδr = −0.3893 rad−1 mδzr = −0.2152 rad−1 mδxa = −0.3022 rad−1
mω̄ −1 mδxr = −0.06492 rad−1
z = −0.5053 rad
z

mz = −0.05272 rad−1
ω̄x
mω̄
x = −0.3968 rad
x −1

m x = −0.007 rad−1
ω̄z

Aiming at the verification of the accuracy of the aerodynamic parameters, the aircraft dynamic
and kinematic model should be calculated with the input data of the actuators in the real flight. So the
correspondence between the width of PWM signal , d PW M , transmitted from the controller and the
actuators’ inputs, namely the rotation rate n and the deflection angles δ, has to be derived. Arifianto [8]
establishes the model by testing the force of the propeller corresponding to the rotation rate, but the
airspeed is ignored in the static experiment. Herein, the ”n-d PW M ” relationship is constructed directly,
and the thrust will be calculated with the equation in Table 1. An experiment of testing the rotation rate
of the propeller is firstly carried out, and the test and the relevant experimental apparatus is shown
in Figure 4. The PWM signal width of the channel of the throttle d PW M ranges from 2202 to 3879.
The rotation rate of the propeller is measured by a handheld digital tachometer corresponding to the
different position of the joystick, which is limited to 70% of the full range for safety. After several tests
and second order polynomial fit, the equation could be obtained in (23), and the maximal rotation rate
could reach 8966.3 RPM, namely 149.4 r/s. Then, the relationship between the deflection angle of the
control surfaces (measured by a digital goniometer) and the PWM width is fitted linearly in (24)–(26).

n P = (−0.001571d2n + 13.18dn − 21300)/60 (23)


(
49.2480 − 0.0162de 2202 ≤ de ≤ 3029
δe (◦ ) = (24)
60.1920 − 0.0198de 3029 < de ≤ 3879
(
◦ 60.7800 − 0.0200dr 2200 ≤ dr ≤ 3040
δr ( ) = (25)
60.7800 − 0.0200dr 3040 < dr ≤ 3879
(
◦ 55.5840 − 0.0180d a 2199 ≤ d a ≤ 3085
δa ( ) = (26)
58.6720 − 0.0190d a 3085 < d a ≤ 3878
Sensors 2018, 18, 206 11 of 18

Figure 4. Rotation Rate Test.

5. Flight Test and Analysis


Real flight tests were carried out in an open playground in Figure 5. The aircraft Extra300 was
basically set to circle around above the field, controlled by a pilot on the ground. Several flight tests
were completed and one of the experiment data was chosen to identify the aerodynamic derivatives,
and then to verify the improvement compared with that calculated by the semi-empirical method.
The 3D flight path is shown in Figure 5. It is truncated from the whole flight, lasting 3 min. To make it
more clear to recognize the changing of the path, it is decomposed into three parts, from point A to D.
In the flight route, there exists circling, slight climbing and diving, accelerating and decelerating and
other maneuvering, which contribute to the observation. Besides, the measured speed is estimated as
ground speed not the airspeed as the wind is neglected.

Path1 Path2 Path3

Figure 5. Flight Test.

The identification was operated according to the analysis in Section 3. The data fusion frequency
is 10 Hz and the standard deviation of white noise is set according to the precision of the INS/GPS
navigation results. The estimation results of the error of aerodynamic parameters are shown in
ω̄y
Figure 6a–h, and all of them could converge to a certain range. Some coefficients, like δmδye , δmy , δmδzr ,
δmω̄ x
x , converge piece-wisely. It may be caused by the changing of the real-time flight condition,
and it matters slightly in this estimation process. The average of the estimation results in the
converging section are calculated and approximated, shown in Table 3. Then, the estimated errors of
the aerodynamic derivatives are subtracted or compensated from the initial derivatives calculated by
semi-empirical method, listed in Table 2, thus further correcting the aircraft model.
Sensors 2018, 18, 206 12 of 18

δCZ0 δCX0
2 0.4

0.2
1
0
0
−0.2

−1 −0.4
0 20 40 60 80 100 120 140 160 180 0 20 40 60 80 100 120 140 160 180
t/s t/s
δCZα δCXA
6 0.4

4 0.2
rad−1

2 0

0 −0.2
−2
0 20 40 60 80 100 120 140 160 180 −0.4
0 20 40 60 80 100 120 140 160 180
t/s t/s

(a) Estimation Results of δCZ0 , δCZα (b) Estimation Results of δCX0 , δCXA
δmy0 δmδye
0.5 10

rad−1
0
0

−0.5 −5
0 20 40 60 80 100 120 140 160 180 0 20 40 60 80 100 120 140 160 180
t/s t/s
δmαy ω̄
δmy y
0.5 20

0
10
rad−1

rad−1

−0.5
0
−1

−1.5 −10
0 20 40 60 80 100 120 140 160 180 0 20 40 60 80 100 120 140 160 180
t/s t/s
ω̄y
(c) Estimation Results of δmy0 , δmαy (d) Estimation Results of δmδye , δmy
δCYβ δmβz
2 1.5

1 1
rad−1

rad−1

0 0.5

−1 0

−2 −0.5
0 20 40 60 80 100 120 140 160 180 0 20 40 60 80 100 120 140 160 180
t/s t/s
δCYδr δmδzr
3 1

2
0
rad−1

rad−1

1
−1
0

−1 −2
0 20 40 60 80 100 120 140 160 180 0 20 40 60 80 100 120 140 160 180
t/s t/s
β β
(e) Estimation Results of δCY , δCYδr (f) Estimation Results of δmz , δmδzr
δmω̄z z δmδxa
2 0.6

0.4
0
rad−1

rad−1

0.2
−2
0

−4 −0.2
0 20 40 60 80 100 120 140 160 180 0 20 40 60 80 100 120 140 160 180
t/s t/s
δmβx δmω̄xx
5 1

0.5
rad−1

rad−1

0 0

−0.5

−5 −1
0 20 40 60 80 100 120 140 160 180 0 20 40 60 80 100 120 140 160 180
t/s t/s
β
(g) Estimation Results of δmω̄ z
z , δm x (h) Estimation Results of δmδxa , δmω̄
x
x

Figure 6. Aerodynamic Parameters Error Estimation Results.


Sensors 2018, 18, 206 13 of 18

Table 3. Parameter Error Estimation Results.

Parameter Result Parameter Result


β
δCZ0 0.000 δCY 1.000
δCZα 3.000 δCYδr 1.600
β
δCZδe 0.000 δmz 0.000
δr
δCX0 −0.100 δmz −0.166
δCXA −0.200 δmω̄
z
z
−0.110
ω̄x
δmy0 0.150 δmz 0.147
β
δmαy 0.100 δm x 0.000
δmδye 1.000 δmδxa 0.200
ω̄
δmy y 5.000 δmδxr 0.000
δmω̄
x
x
0.200
ω̄z
δm x 0.593

With the initial motion states provided by the INS/GPS integration in XSENS MTi-G, the
navigation results from the integration of the aircraft motion equations of (14)–(16) is calculated.
By comparing with the INS/GPS navigation results, on one hand, the rationality of the aerodynamic
parameters computed by the semi-empirical method will be proved, on the other hand, the effect of
latter compensation process is also verified, which is a unique method to validate the correctness of
the calculation and identification results. Figures 7–9 show the comparison of the errors of navigation
results between the empirically calculated model and the compensated one by subtracting the INS/GPS
navigation results.
Firstly, it is demonstrated that the semi-empirically calculated parameters based on aerodynamics
are rational and useful, which is the premise for further analysis. After the compensation of the error
of the aerodynamic derivatives, there shows an improvement in the navigation results. Specifically, the
error of the latitude and longitude of the two processes are similar, and it is calculated that the average
error of the latitude and longitude are about −39.1824 m and 5.5703 m respectively. The average of
height error decreases to about 177.6206 m from 272.3769 m after the identification. The average errors
of the horizontal velocity are alike, namely δVn , δVe are about −0.9853 m/s, 0.8024 m/s respectively
while δVd changes from 3.1785 m/s to 1.8770 m/s. Although there exists just a slight improvement in
the horizontal velocity, it is obviously figured that the compensated model are more stable and owns a
better tendency compared with the INS/GPS results. The average error of roll angle decreases from
−44.4609◦ to−11.1175◦ . The average error of pitch angle, around 12.6113◦ , are similar but more stable.
Yaw angle, with the average error of −26.4751◦ , are more rational and precise after the compensation
as the wide swing are taken out.
Despite the improvement after the correction process, differences still exist between the
compensated aircraft model and the precise integration results. The model of the rotation rate
of propeller or the deflection angle of the actuators maybe inaccurate, but it has little effect on
the navigation results according to the experiment. It is largely caused by the inaccuracy of the
identification of aerodynamic parameters. The flight condition, like the large variation of the
aerodynamic angle, will influence the estimation precision. It could be seen from the changing of α, β in
Figure 10. As seen from the pictures, it doesn’t match the small angle-of-attack assumption sometimes,
so some aerodynamic derivatives are no longer constants. Compared with the navigation error in
Figure 8, the error just becomes larger when the angle-of-attack increases largely. The phenomenon
also take place in other flight tests, and the estimation results of the parameters’ error are a little
different in different flight sorties concerning the different flight condition. To address the problem,
the nonlinear aerodynamic derivatives should be taken into account as the simplified modeling of the
coefficients and derivatives leads to the uncertainty. For example, the slope of lift curve CZα is variant,
and some nonlinear derivatives associated with α, β should be augmented to the states to be estimated.
The variant coefficients should be estimated precisely and compensated piece-wisely in a real time.
Sensors 2018, 18, 206 14 of 18

Furthermore, more sensors measuring the aerodynamic angles α, β and the airspeed should be added
to take the changing circumstances into account.

Empirical
Latitude Error (δ L) Compensated
100

m
−100

0 20 40 60 80 100 120 140 160 180


Longitude Error (δ λ)
100
m

−100
0 20 40 60 80 100 120 140 160 180
Height Error (δh)
0
−200
m

−400
−600
0 20 40 60 80 100 120 140 160 180
t/s

Figure 7. Position Error.

Empirical
North Velocity Error (δ VN ) Compensated
20
m/s

0
−20

0 20 40 60 80 100 120 140 160 180


East Velocity Error (δ VE )
20
m/s

0
−20
0 20 40 60 80 100 120 140 160 180
Down Velocity Error (δVD )
20
10
m/s

0
−10
0 20 40 60 80 100 120 140 160 180
t/s

Figure 8. Velocity Error.

Empirical
Roll Error (δ φ) Compensated
200
deg

−200
0 20 40 60 80 100 120 140 160 180
Pitch Error (δθ)
100
deg

−100
0 20 40 60 80 100 120 140 160 180
Yaw Error (δ ψ)
400
200
deg

0
−200
−400
0 20 40 60 80 100 120 140 160 180
t/s

Figure 9. Attitude Error.


Sensors 2018, 18, 206 15 of 18

100
α

50

deg
0

−50

−100
0 20 40 60 80 100 120 140 160 180
t/s

40
β
20
0
deg

−20
−40
−60
0 20 40 60 80 100 120 140 160 180
t/s

Figure 10. Aerodynamic Angle in the Flight.

6. Conclusions
The paper mainly studied the method of calculating the aerodynamic parameters of the
small-scaled fixed-wing aircraft, and it is an important part in the establishment of the aircraft dynamic
and kinematic model. It contains 2 sequential processes:

(1) Prepared with the related structural parameters and the position of gravity center, the
aerodynamic derivatives are semi-empirically calculated based on fundamental aerodynamics.
The simplified structure and the specific flight condition reduce the complexity. So, the
aerodynamic analysis could be operated in a fast, effective and low-cost way to originally
provide a serial rational parameters. The paper introduces and analyzes the entire course of the
calculating of the longitudinal and lateral derivatives, and it is helpful for the modeling of all
the other similar small aircraft;
(2) The initially obtained parameters is further identified or compensated with the real flight data.
Given the airborne MEMS-based inertial sensors with finite precision, the relatively accurate
velocity and attitude from the integration of INS/GPS was employed as observations to estimate
the error of the aerodynamic parameters. It may have the problem of lack of observations, but
it could be solved with the abundant maneuvering flight. Then, the aerodynamic parameters
could be corrected with the estimated error of the derivatives and the aircraft model could be
more accurate and reliable.

The real flight test demonstrates a promising result of the above 2 methods. Abide by the
traditional aerodynamics and the empiricism from flight database, a rational aircraft model is
constructed. Although about 20% or more error exists in the aerodynamic coefficients, it still provides
a basic prototype or a baseline model. During the circling flight test, the controller adjusted the
model aircraft all the time, thus providing the required maneuvers to increase the observability of
all the derivatives in the identification process. And the accuracy of the compensated parameters
was verified by integrating the aircraft dynamic and kinematic model, then comparing with the
INS/GPS integration results, which is different from the traditional verification. The accuracy of the
method synthesizes the precision of the model itself including the structural and aerodynamic part
as well as the precision of the actuators’ inputs. But the correctness of the aerodynamic parameters
affect predominantly.
As discovered from the tests, the ideal aerodynamic parameters, namely the true values, are hard
to reach. Firstly, the inaccuracy of the structural parameters, the possible shifting of the position of
the center of gravity, and the related assumptions and simplifications could make an effect on the
Sensors 2018, 18, 206 16 of 18

semi-empirical analysis. Then, the variant flight condition makes it more difficult for the filter to
identify the variant parameters correctly and instantaneously, especially encountering the disturbance
of the wind or gust. In the subsequent study or experiment, more aerodynamic sensors, like the
low-cost or homemade vanes measuring the angle-of-attack and the side-slip angle, will be equipped
to take the wind into consideration. Then, the relatively accurate aircraft model are expected to be
used as an assistance for the fast divergent pure INS in the future research.

Acknowledgments: The study and the related experiment was supported by the MEMS Research Center of
Nanjing University of Science and Technology. The work was also supported by a Natural Science Foundation
Project (No.51375244), China.
Author Contributions: Jieliang Shen conducted all the research and wrote the paper. Yan Su guided and revised
the thoughts of the research. Qing Liang helped conceiving and designing the experiment. Xinhua Zhu contributed
to the improvement of the thoughts and the writing of the paper.
Conflicts of Interest: The authors declare no conflict of interest.

Nomenclature
Abbreviation
ADM Aircraft Dynamic Model
AMM Aircraft Motion Model
CFD Computational Fluid Dynamics
EKF Extended Kalman Filter
GPS Global Positioning System
I MU Inertial Measurement Unit
I NS Inertial Navigation System
MEMS Micro-Electro Mechanical Systems
MLE Maximum Likelihood Estimation
PMW Pulse Width Modulation
U AV Unmanned Aerial Vehicle
UKF Unscented Kalman Filter
Coordinate
g Earth-surface inertial reference frame
pointing North, East, Down
b Aircraft-body coordinate frame
a Wind coordinate frame
Cgb Coordinate rotation matrix from g to b
Cab Coordinate rotation matrix from a to b
Aircraft
m Total mass, kg
I Moment of inertia, kg·m2
xG , yG , zG Position of center of gravity, m
xF Aerodynamic focus, m
x pW , x pB , x pWB Aerodynamic center of wing, body and wing-body, m
x pHT , x pHTδe Aerodynamic center of horizontal tail and elevator, m
l Wing span, m
c̄ Mean aerodynamic chord, m
S Reference area, m2
SCW Cantilever wing area, m2
SB Fuselage maximum cross-sectional area, m2
S HT Horizontal tail area, m2
λ Aspect ratio
Sensors 2018, 18, 206 17 of 18

D̄ Fuselage diameter-to-wing span ratio


D Diameter of propeller, m
ϕW , ϕ S Incidence angle of wing, horizontal tail, degree
Kαα , Kφ0 Wing-body interference factor
(Kαα ) HT , (Kφ0 ) HT Horizontal tail-body interference factor
FT Thrust of propeller, N
FA = [ X, Y, Z ] T Aerodynamic forces in wind coordinates,
Drag, Side Force and Lift, N
Mb = [ M x , My , Mz ] T Aerodynamic moments in body coordinates,
Roll, Pitch and Yaw moment, N· m
CFT1 , CFT2 , CFT3 Dimensionless thrust coefficients
CX , CY , CZ Dimensionless aerodynamic-force coefficients
m x , my , mz Dimensionless aerodynamic-moment coefficients
Va Airspeed, m/s
α Angle-of-attack, rad
β Angle of sideslip, rad
ε Angle of downwash, rad
ρ Air density, kg/m3
Q Dynamic pressure, Pa
k HT Airflow block factor
δe , δr , δa Deflection angle of elevator, rudder and aileron, rad
nP Revolutions per second of propeller, r/s
d n , d e , dr , d a PWM width corresponding to the channel of
thrust, elevator, rudder and aileron
φ, θ, ψ Roll, pitch, yaw angle rad
V b = [u, v, w] T Body axis velocities, m/s
ω b = [ p, q, r ] T Body axis angular rates, rad/s

References
1. Austin, R. Unmanned Aircraft Systems: UAVS Design, Development and Deployment, 1st ed.; John Wiley & Sons:
Chichester, UK, 2011.
2. Ducard, G.J. Fault-Tolerant Flight Control and Guidance Systems: Practical Methods for Small Unmanned Aerial
Vehicles; Springer Science & Business Media: London, UK, 2009.
3. Klein, V.; Morelli, E.A. Aircraft System Identification: Theory and Practice; American Institute of Aeronautics
and Astronautics: Reston, VA, USA, 2006.
4. Marqués, P.; Da Ronch, A. Advanced Uav Aerodynamics, Flight Stability and Control: Novel Concepts, Theory and
Applications; John Wiley & Sons: Southport, UK, 2017.
5. Wisnoe, W.; Nasir, R.E.M.; Kuntjoro, W.; Mamat, A.M. Wind tunnel experiments and CFD analysis of
Blended Wing Body (BWB) Unmanned Aerial Vehicle (UAV) at mach 0.1 and mach 0.3. In Proceedings of
the 13th International Conference on Aerospace Sciences & Aviation Technology, ASAT-13, Cairo, Egypt,
26–28 May 2009.
6. Williams, J.E.; Vukelich, S.R. The USAF Stability and Control DATCOM, Users Manual; McDonnell Douglas
Astronautics Company: Berkeley, MO, USA, 1999; Volume 1.
7. STABILITY, USAF. Control Datcom. In Air Force Flight Dynamics Laboratory, Wright-Patterson Air Force
Base, Ohio; 1972.
8. Arifianto, O.; Farhood, M. Development and modeling of a low-cost unmanned aerial vehicle research
platform. J. Intell. Robot. Syst. 2015, 80, 139.
9. Etkin, B. Dynamics of Atmospheric Flight; Dover Publications: New York, NY, USA, 2005.
10. Nelson, R.C. Flight Stability and Automatic Control, 2nd ed.; WCB/McGraw Hill: New York, NY, USA, 1998.
11. Phillips, W.F. Mechanics of Flight; John Wiley & Sons: Chichester, UK, 2004.
Sensors 2018, 18, 206 18 of 18

12. Schmidt, L.V. Introduction to Aircraft Flight Dynamics; American Institute of Aeronautics and Astronautics
(AIAA): Reston, VA, USA, 1998.
13. Lebedev, A. Flight Dynamics of the Unmanned Aerial Vehicles; Zhang, B., Translator; National Defence Industry
Press: Beijing China, 1964.
14. Xu, M.; An, X. Analysis and Calculation of Aerodynamic Characteristics of Aircraft; Northwestern Polytechnic
University Press: Xi’an, China, 2012.
15. Tischler, M.B.; Remple, R.K. Aircraft and Rotorcraft System Identification; AIAA Education Series; American
Institute of Aeronautics and Astronautics: Reston, VA, USA, 2006.
16. Ljung, L. System identification. In Signal Analysis and Prediction; Springer: New York, NY, USA, 1998;
pp. 163–173.
17. Guo, M.; Su, Y.; Gu, D. System identification of the quadrotor with inner loop stabilisation system. Int. J.
Model. Identif. Control 2017, 28, 245–255.
18. Ghosh, A.; Raisinghani, S. Frequency-domain estimation of parameters from flight data using neural
networks. J. Guid. Control Dyn. 2001, 24, 525–530.
19. Chen, S.; Billings, S. Neural networks for nonlinear dynamic system modelling and identification.
Int. J. Control 1992, 56, 319–346.
20. Morelli, E.A. Real-time aerodynamic parameter estimation without air flow angle measurements. J. Aircr.
2012, 49, 1064–1074.
21. Dorobantu, A.; Ozdemir, A.A.; Turkoglu, K.; Freeman, P.; Murch, A.; Mettler, B.; Balas, G. Frequency
domain system identification for a small, low-cost, fixed-wing uav. In Proceedings of the AIAA Guidance,
Navigation, and Control Conference, Portland, OR, USA, 8–11 August 2011, p. 6719.
22. Dorobantu, A.; Murch, A.; Mettler, B.; Balas, G. System identification for small, low-cost, fixed-wing
unmanned aircraft. J. Aircr. 2013.
23. Garcia-Velo, J.; Walker, B.K. Aerodynamic parameter estimation for high-performance aircraft using extended
Kalman filtering. J. Guid. Control Dyn. 1997, 20, 1257–1260.
24. Tang, Y.; Wu, Y.; Wu, M.; Wu, W.; Hu, X.; Shen, L. INS/GPS integration: Global observability analysis.
IEEE Trans. Veh. Technol. 2009, 58, 1129–1142.
25. Jategaonkar, R. Flight Vehicle System Identification: A Time Domain Methodology; American Institute of
Aeronautics and Astronautics: Reston, VA, USA, 2006; Volume 216.
26. Von Hoffer, N. System Identification of a Small Low-Cost Unmanned Aerial Vehicle Using Flight Data from
Low-Cost Sensors; Utah State University: Logan, UT, USA, 2015.
27. Chowdhary, G.; Jategaonkar, R. Aerodynamic parameter estimation from flight data applying extended and
unscented Kalman filter. Aerosp. Sci. Technol. 2010, 14, 106–117.
28. Koifman, M.; Bar-Itzhack, I. Inertial navigation system aided by aircraft dynamics. IEEE Trans. Control
Syst. Technol. 1999, 7, 487–493.
29. Yang, Z.; Igarashi, H.; Martin, M.; Hu, H. An experimental investigation on aerodynamic hysteresis of
a low-Reynolds number airfoil. In Proceedings of the 46th AIAA Aerospace Sciences Meeting and Exhibit,
Reno, NV, USA, 7–10 January 2008.
30. Xsens, T. MTi-G User Manual and Technical Documentation; Xsens Technologies B.V.: Enschede,
The Netherlands, 2010.
31. Jung, D.; Levy, E.; Zhou, D.; Fink, R.; Moshe, J.; Earl, A.; Tsiotras, P. Design and development of a low-cost
test-bed for undergraduate education in UAVs. In Proceedings of the 44th IEEE Conference on Decision and
Control, and 2005 European Control Conference, Seville, Spain, 12–15 December 2005, pp. 2739–2744.
32. Owens, D.B.; Cox, D.E.; Morelli, E.A. Development of a low-cost sub-scale aircraft for flight research:
The FASER project. In Proceedings of the 25th AIAA Aerodynamic Measurement Technology and Ground
Testing Conference, Fluid Dynamics and Co-located Conferences, San Francisco, CA, USA, 5–8 June 2006.

c 2018 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access
article distributed under the terms and conditions of the Creative Commons Attribution
(CC BY) license (http://creativecommons.org/licenses/by/4.0/).

You might also like