Computer-Implemented Method and System for Modelling Performance of a Fixed-Wing Aerial Vehicle with Six Degrees of Freedom

Abstract

Computer-implemented method and system for modelling performance of a fixed-wing aerial vehicle (AV) with six degrees of freedom. The system comprises a collecting unit (330) to collect data from a plurality of modelling measures and modelling maneuvers; a processing unit (340) to communicate with a collecting unit (330). The processing unit (340) further sequentially processes data sets from a plurality of modelling measures and to determine models to generate an accurate APM. Modelling measures and modelling maneuvers are designed to modify an influence on the of variables of a model of the AV (100).

Claims

1. A computer-implemented method for modeling performance of a fixed-wing aerial vehicle (AV) with six degrees of freedom, sequentially comprising the steps of: collecting a first data set from a plurality of fuel consumption modelling measures and determining a fuel consumption model based on the first data set; collecting a second data set from a plurality of thrust modelling maneuvers and determining a thrust model based on the second data set; collecting a third data set from a plurality of aerodynamic forces modelling maneuvers and determining aerodynamic forces model based on the third data set; collecting a fourth data set from a plurality of propulsive moments modelling maneuvers and inertia matrix modelling maneuvers and determining propulsive moments model and inertia matrix based on the fourth data set; and collecting a fifth data set from a plurality of aerodynamics moments modelling maneuvers and determining aerodynamics moments model based on the fifth data set; wherein modelling measures and modelling maneuvers are performed to modify an influence on the AV of one or more variables of a model to be determined.

2. The computer-implemented method of claim 1, further comprising a step of measuring variables of a model with a state estimator of the AV and air data system of the AV.

3. The computer-implemented method of claim 1, wherein determining a model further comprises applying a least square estimate to a collected data set.

4. The computer-implemented method of claim 1, further comprising a step of configuring a flight control system to automatically perform modelling maneuvers in flight.

5. The computer-implemented method of claim 1, wherein modelling maneuvers comprise at least one of the following control loops: AOS-on-rudder, altitude-on-bank speed-on-elevator, AOA-on-elevator, pitch rate-on-elevator and bank-on-ailerons.

6. The computer-implemented method of claim 1, wherein fuel consumption modelling measures are performed for a plurality of values of air density, outside air temperature (OAT) and throttle level.

7. The computer-implemented method of claim 1, wherein thrust modelling maneuvers are performed for a plurality of values of barometric altitude, mass variation and throttle level under a condition of coordinated flight.

8. The computer-implemented method of claim 7, wherein thrust modelling maneuvers comprise an altitude-on-bank control loop and a speed-on-elevator control loop.

9. The computer-implemented method of claim 1, wherein aerodynamic forces modelling maneuvers comprise AOA-on-elevator, AOS-on-rudder and speed-on-elevator control loops.

10. The computer-implemented method of claim 1, wherein propulsive moments modelling maneuvers comprise AOS-on-rudder and bank-on-ailerons control loops under a condition of coordinated flight.

11. The computer-implemented method of claim 1, wherein aerodynamic moments modelling maneuvers comprise AOS-on-rudder, bank-on-ailerons and q-on-elevator control loops.

12. A system for modeling performance of a fixed-wing aerial vehicle (AV) with six degrees of freedom, the system comprising: a collecting unit configured to collect data from a plurality of modelling measures and modelling maneuvers; a processing unit configured to communicate with a collecting unit, wherein the processing unit is further configured to sequentially process: a first data set from a plurality of fuel consumption modelling measures and to determine a fuel consumption model based on the first data set; a second data set from a plurality of thrust modelling maneuvers and to determine a thrust model based on the second data set; a third data set from a plurality of aerodynamic forces modelling maneuvers and to determine aerodynamic forces model based on the third data set; a fourth data set from a plurality of propulsive moments modelling maneuvers and inertia matrix modelling maneuvers and to determine propulsive moments model and inertia matrix based on the fourth data set; a fifth data set from a plurality of aerodynamics moments modelling maneuvers and to determine aerodynamics moments model based on the fifth data set; wherein modelling measures and modelling maneuvers are performed to modify an influence on the AV of one or more variables of a model to be determined.

13. The system of claim 12, wherein the collecting unit is further configured to instruct a state estimator and air data system to read out sensors values of the AV during modelling maneuvers.

14. The system of claim 12, wherein the processing unit is further configured to instruct a flight control system to automatically perform modelling maneuvers in flight.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0013] A series of drawings which aid in better understanding the disclosure and which are presented as non-limiting examples and are very briefly described below.

[0014] FIG. 1 illustrates a three-dimensional system used to define the orientation of an aircraft.

[0015] FIG. 2 illustrates a flow diagram depicting one embodiment of the disclosed method.

[0016] FIG. 3 illustrates a block system according to an embodiment of the disclosed system.

DETAILED DESCRIPTION

[0017] For the purposes of explanation, specific details are set forth in order to provide a thorough understanding of the present disclosure. However, it is apparent to one skilled in the art that the present disclosure may be practiced without these specific details or with equivalent arrangements.

[0018] FIG. 1 illustrates the body fixed system (BFS) of a fixed-wing aerial vehicle (AV) 100. A three-dimensional coordinate system is defined through the center of gravity (CoG) with each axis of this coordinate system perpendicular to the other two axes. The AV 100 typically includes a body or fuselage 102, fixed wings 104 that are coupled to the fuselage 102 and a tail 105 that includes a horizontal stabilizer 106 and a vertical stabilizer 108. The AV 100 can rotate about three axes, namely, a longitudinal x-axis 110, a lateral y-axis 112, and a vertical or directional z-axis 114. Roll, pitch and yaw refer to rotations of the AV 100 about the respective axes 110, 112, 114. Pitch refers to the rotation (nose up or down) of the AV 100 about the lateral y-axis 112, roll refers to the rotation of the AV 100 about the longitudinal x-axis 110, and yaw refers to the rotation of the AV 100 about the vertical or directional z-axis 114.

[0019] An AV 100 typically include four controls, namely throttle, elevator 116, ailerons 124, and vertical rudder 118. Throttle is used to operate the engines 134 and generate a forward force. Control surfaces: elevator 116, ailerons 124 and vertical rudder 118 allow to aerodynamically control the AV 100 for a stable flight.

[0020] Elevator 116 is used to raise or lower the nose of the AV 100. Ailerons 124 are located on the outside of the wings 104 and are used to rotate AV 100. The ailerons 124 move anti-symmetrically (i.e. one goes up, the one on the opposite side low). Vertical rudder 118 is used so that the aircraft does not skid.

[0021] It is assumed that the AV 100 is equipped with an aircraft flight control system. The AV 100 may be commanded with the flight control system to automatically perform certain maneuvers and collect data (via telemetry) in real time that serve to generate an APM. At least some of these tasks may be coded in a computer program. Some of these computer program instructions, when running, may steer the AV 100 to perform modelling maneuvers. Said maneuvers advantageously provide data that can be collected and used to generate APM features in a sequence of different steps to arrive at an accurate APM that reproduces the behavior of the AV 100.

[0022] FIG. 2 illustrates several steps of an embodiment according to the present disclosure to generate an APM for a particular fixed-wings AV. Generally, the model determined in a step serves for determining another model in a subsequent step.

[0023] A first step of fuel consumption model determination 210 is based on fuel consumption measures 212 performed under different conditions of air density, outside air temperature (OAT) and throttle level (ET).

[0024] A second step of thrust model determination 220 is based on thrust modeling maneuvers 222 performed under conditions of coordinated flight (a flight without sideslip), constant barometric altitude for a plurality of AV mass variations, a plurality throttle levels and a plurality of angle of attack (AOA) values for different air densities, OAT and throttle levels.

[0025] A third step of aerodynamic forces model determination 230 is based on data collected from aerodynamic forces modelling maneuvers 232 like AOA-on-elevator, AOS-on-rudder (angle of sideslip) and speed-on-elevator.

[0026] A fourth step of propulsive moments and inertia matrix determination 240 is based on data collected from propulsive moments modelling maneuvers 242 from control loops of AOS-on-rudder and bank-on-ailerons. Optionally, maneuvers 242 may be performed under conditions of coordinated flight or actuating the elevator.

[0027] Inertia matrix moments modelling maneuvers 243 may be performed under conditions of lateral perturbation via rudder actuation for several AV mass variations.

[0028] A fifth step of aerodynamic moments determination 250 is based on data collected from performing maneuvers under certain conditions 252 like AOS-on-rudder, bank-on-ailerons and q-on-elevator.

[0029] Once this sequence of steps is performed for the considered fixed-wings AV, its corresponding APM can be obtained.

[0030] FIG. 3 shows an embodiment of a system to generate an APM according to the present disclosure. As to instrumentation and data acquisition, the AV 100 is assumed to be equipped with a state estimator 320 an air data system 310. The state estimator 320 may fusion GPS 322, IMU 325 (e.g. solid-state accelerometers 329 and gyros 324), magnetometer 328, wind vanes 326 (providing direct observation of angle-of-attack and angle-of-sideslip). The air data system 310 provides observations of true airspeed, pressure and temperature to produce an estimate (e.g. in the sense of the Extended Kalman Filter).

[0031] The aircraft flight control system of the AV 100 includes a processing unit 340 that may instruct operating mechanisms that actuate control surfaces to modify direction in flight. The processing unit 340 also may manage throttle that controls engines to modify speed.

[0032] Several actions can be performed as outlined in FIG. 2, including flight modelling maneuvers 222, 232, 242, 243 and 252 and associated measures for collecting data in order to sequentially build an APM for the AV 100. At least some of them may be automated and implemented in flight via processing unit 340 to reduce complexity and facilitate computations to generate a realistic APM. Main maneuvers of assistance to generate the APM are: [0033] AOS-on-rudder [0034] Altitude-on-bank [0035] Speed-on-elevator [0036] AOA-on-elevator [0037] Pitch rate-on-elevator

[0038] These will be evidenced in the following example expressed with a more detailed mathematical formulation.

Formal Example of APM Identification

[0039] The present example mathematically demonstrates how a series of modelling maneuvers and measures under particular conditions are used to generate data to be collected in order to determine parameters of each particular model by means of least squares (LS).

Position (WGS84):

[0040] {, , h} Geodetic coordinates
Linear speed (LLS):
u Absolute horizontal speed
True bearing
w Absolute vertical speed
v Module of the absolute speed
Absolute (geometric) path angle
Absolute bank angle
u.sub.TAS Horizontal component of the true airspeed
.sub.TAS Bearing of the true airspeed
w.sub.TAS Vertical component of the true airspeed
v.sub.TAS Module of the true airspeed v.sub.TAS={square root over (u.sub.TAS.sup.2+w.sub.TAS.sup.2)}
.sub.TAS Aerodynamic path angle
.sub.TAS Aerodynamic bank angle
u.sub.WIND Horizontal component of the wind speed
.sub.WIND Bearing of the wind speed
w.sub.WIND Vertical component of the wind speed
v.sub.WIND Module of the wind speed v.sub.WIND={square root over (u.sub.WIND.sup.2+w.sub.WIND.sup.2)}
Linear acceleration (BFS):

[00001] $[\begin{matrix} a_{1}^{BFS} \\ a_{2}^{BFS} \\ a_{3}^{BFS} \end{matrix}]$

BFS components of the non-gravitational accelerations (likewise sensed by a 3-axis accelerometer)

Attitude (BFS):

[0041] {,,} Euler angles
Angle-of-attach (AOA)
Angle-of-sideslip (AOS)
Angular speed (BFS):

[00002] $[\begin{matrix} p \\ q \\ r \end{matrix}]$

BFS components of the angular speed of the AV with respect to LLS

[00003] $[\begin{matrix} \dot{} \\ \dot{} \end{matrix}]$

Derivatives of AOA and AOS

[0042] Angular acceleration (BFS): [angular acceleration not typically a native observable in low-cost IMUs, which implies numerical derivation]

[00004] $[\begin{matrix} \dot{p} \\ \dot{q} \\ \dot{r} \end{matrix}]$

BFS components of the angular acceleration of the AV with respect to LLS

[0043] The accelerometer measures {a.sub.1.sup.BFS, a.sub.2.sup.BFS, a.sub.3.sup.BFS} provide the means to split the 3 equations that govern the linear dynamics (i.e. the motion of the AV's gravity center) in two systems where either non-gravitational only forces (i.e. aerodynamic and propulsive) or gravitational-only forces (i.e. weight and Coriolis inertia force) appear. The first system, represented by expression [1] (with further details in expression [2]) is the one off interest to our purpose:

[00005] $\begin{matrix} [\begin{matrix} {Tt}_{1} - D \\ {Tt}_{2} - Q \\ {Tt}_{3} - L \end{matrix}] = m [\begin{matrix} a_{1}^{BFS} .Math. \cos .Math. .Math. .Math. .Math. \cos .Math. .Math. + a_{2}^{BFS} .Math. \sin .Math. .Math. + a_{3}^{BFS} .Math. .Math. \sin .Math. .Math. .Math. .Math. \cos .Math. .Math. \\ - a_{1}^{BFS} .Math. \cos .Math. .Math. .Math. .Math. \sin .Math. .Math. + a_{2}^{BFS} .Math. \cos .Math. .Math. - a_{3}^{BFS} .Math. .Math. \sin .Math. .Math. .Math. .Math. \sin .Math. .Math. \\ - a_{1}^{BFS} .Math. \sin .Math. .Math. + a_{3}^{BFS} .Math. \cos .Math. .Math. \end{matrix}] & [1] \\ t_{1} (,,, .Math.) = \cos .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. - \sin .Math. .Math. .Math. .Math. \sin .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. - \sin .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. .Math. \sin .Math. .Math. .Math. .Math. .Math. t_{2} (,,, .Math.) = - (\cos .Math. .Math. .Math. .Math. \sin .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. + \cos .Math. .Math. .Math. .Math. \sin .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. - \sin .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. .Math. \sin .Math. .Math. .Math.) .Math. .Math. t_{3} (,, .Math.) = - (\sin .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. + \cos .Math. .Math. .Math. .Math. \sin .Math. .Math. .Math.) & [2] \end{matrix}$

[0044] The variables that take part of expression [1] can be classified as follows: [0045] a.sub.1.sup.BFS, a.sub.2.sup.BFS, a.sub.3.sup.BFS, , are observable aspects of the AV state at any time, provided by the estate estimator [0046] , define the direction of the thrust force with respect to BFS; they are unknowns, in principle, but small constant angles (i.e., v<<1 and E<<1) [0047] m (the actual mass) is a slowly varying variable [0048] L, Q, D, T respectively represent the aerodynamic lift, side and drag forces and the thrust force

[0049] In the quasi-steady state assumption ({dot over ()}{dot over ()}0), aerodynamic and propulsive forces can be expressed as:

[00006] $\begin{matrix} L = \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. {SC}_{L} (,, M, \hat{p}, \hat{q}, \hat{r}, {.Math.}_{h}, {.Math.}_{a}, {.Math.}_{r}, {.Math.}_{e}) & [3] \\ Q = \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. {SC}_{Q} (,, M, \hat{p}, \hat{q}, \hat{r}, {.Math.}_{h}, {.Math.}_{a}, {.Math.}_{r}, {.Math.}_{e}) & [4] \\ D = \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. {SC}_{Q} (,, M, \hat{p}, \hat{q}, \hat{r}, {.Math.}_{h}, {.Math.}_{a}, {.Math.}_{r}, {.Math.}_{e}) & [5] \\ T = W_{MTOW} .Math. C_{T} (,, M, {.Math.}_{T}) & [6] \\ \hat{p} = \frac{{pb}_{w}}{2 .Math. v_{TAS}} .Math. .Math. b_{w} .Math. .Math. is .Math. .Math. the .Math. .Math. wingspan & [7] \\ \hat{q} = \frac{{qc}_{w}}{2 .Math. v_{TAS}} .Math. .Math. c_{w} .Math. .Math. is .Math. .Math. the .Math. .Math. wing .Math. .Math. mean .Math. .Math. aerodynamic .Math. .Math. chord & [8] \\ \hat{r} = \frac{{rb}_{w}}{2 .Math. v_{TAS}} & [9] \end{matrix}$

[0050] With typical symmetry assumptions:

L=p.sub.0M.sup.2SC.sub.L(,,M,{circumflex over (q)},.sub.h,.sub.e)[10]

Q=p.sub.0M.sup.2SC.sub.Q(,,M,{circumflex over (p)},{circumflex over (r)},.sub.a,.sub.r)[11]

D=p.sub.0M.sup.2SC.sub.D(,,M,{circumflex over (q)},.sub.h,.sub.e)[12]

Step 210 for Determination of a Fuel Consumption Model

[0051] The APM identification process starts with fuel consumption modelling measures 212 as a bench test intended to identify instantaneous fuel consumption as a function of air density, temperature and throttle level:

[00007] $\begin{matrix} \dot{m} = - F & [13] \\ m = m_{0} -_{t_{0}}^{t} .Math. Fdt & [14] \\ F = \frac{w_{MTOW} .Math. a_{0} .Math. .Math. \sqrt{}}{L_{HV}} .Math. C_{F} (,, {.Math.}_{T}) & [15] \end{matrix}$

[0052] In effect, C.sub.F(,,.sub.T) can be identified for the air pressure and temperature conditions existing at the moment of conducting the bench test, which shall measure mass variation at different levels of .sub.T within its range [0,1], fixed during a time interval wide enough to accurately measure fuel consumption. Alternatively (or complementarily) a direct measure of F can be obtained through a precise caudalimeter. Ideally, the bench tests should be repeated for significantly different temperature and pressure conditions in order to capture its dependence with such variables.

[0053] Once the fuel consumption model C.sub.F(,,.sub.T) is identified, expression [14] allows the calculation of the instantaneous mass at any time, as a function of the variation of , and .sub.T recorded over the time elapsed from the initial time t.sub.0 where the initial mass was m.sub.0.

Step 220 for Determination of a Thrust Model

[0054] Lateral aerodynamic force Q is symmetric with respect to its dependency with AOS or, in other words, Q is null for coordinated flight, i.e.:

0.Math.Q0[16]

[0055] Thus, if a basic control loop AOS-on-rudder is implemented to maintain coordinated flight for any value of .sub.a, expression [1] turns into the much simpler form:

[00008] $\begin{matrix} [\begin{matrix} {Tt}_{1} - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. {SC}_{D} \\ {Tt}_{2} \\ {Tt}_{3} - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. {SC}_{L} \end{matrix}] = m [\begin{matrix} a_{1}^{BFS} .Math. \cos .Math. .Math. + a_{3}^{BFS} .Math. \sin .Math. .Math. \\ a_{2}^{BFS} \\ - a_{1}^{BFS} .Math. \sin .Math. .Math. + a_{3}^{BFS} .Math. \cos .Math. .Math. \end{matrix}] & [17] \end{matrix}$

[0056] For AVs that fly at moderate speeds, i.e. at Mach numbers that allow neglecting compressibility effects, the dependency of C.sub.L, C.sub.Q and C.sub.D with the Mach number M can be neglected. Furthermore, for coordinated flights (null AOS), linear aerodynamics theory shows that at moderate AOA, C.sub.D depends quadratically with the AOA while C.sub.L depends linearly with the AOA, which enables relating both coefficients through the classical parabolic drag polar (first order approach):

C.sub.D=C.sub.D,0+C.sub.D,2C.sub.L.sup.2[18]

[0057] Such model can be extended to cope with compressibility effects to certain extent by tailoring the coefficients to a given Mach. Thus, for coordinated level flight at given Mach M and trim condition .sub.h:

C.sub.D(,0,M,0,.sub.h,0)=C.sub.D,0+K.sup.2[19]

C.sub.L(,0,M,0,.sub.h,0)=C.sub.L,1[20]

K(,0,M,0,.sub.h,0)=C.sub.D,2C.sub.L,1.sup.2[21]

with C.sub.D,0, K, C.sub.L,1 and C.sub.D,2 depending on the selected speed M and trim condition .sub.h

[00009] $\begin{matrix} [\begin{matrix} T (\cos .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. - \sin .Math. .Math. .Math. .Math. \sin .Math. .Math. .Math.) - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S (C_{D, 0} + K .Math. .Math.^{2}) \\ - T .Math. .Math. \sin .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. \\ - T (\sin .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. + \cos .Math. .Math. .Math. .Math. \sin .Math. .Math. .Math.) - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. {SC}_{L, 1} .Math. \end{matrix}] = m [\begin{matrix} a_{1}^{BFS} .Math. \cos .Math. .Math. + a_{3}^{BFS} .Math. \sin .Math. .Math. \\ a_{2}^{BFS} \\ - a_{1}^{BFS} .Math. \sin .Math. .Math. + a_{3}^{BFS} .Math. \cos .Math. .Math. \end{matrix}] & [22] \\ [\begin{matrix} T (\cos .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. .Math. .Math. \cos .Math. .Math. - T .Math. .Math. \sin .Math. .Math. .Math. .Math. .Math. \sin .Math. .Math.) - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S (C_{D, 0} + K .Math. .Math.^{2}) \\ - T .Math. .Math. \sin .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. \\ - T .Math. .Math. \cos .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. .Math. .Math. \sin .Math. .Math. - T .Math. .Math. \sin .Math. .Math. .Math. .Math. .Math. \cos .Math. .Math. - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. {SC}_{L, 1} .Math. \end{matrix}] = m [\begin{matrix} a_{1}^{BFS} .Math. \cos .Math. .Math. + a_{3}^{BFS} .Math. \sin .Math. .Math. \\ a_{2}^{BFS} \\ - a_{1}^{BFS} .Math. \sin .Math. .Math. + a_{3}^{BFS} .Math. \cos .Math. .Math. \end{matrix}] & [23] \\ [\begin{matrix} W_{MTOW} .Math. .Math. .Math. \cos .Math. .Math. & - W_{MTOW} .Math. .Math. .Math. \sin .Math. .Math. & 0 & - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S & 0 & - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S .Math. .Math.^{2} \\ 0 & 0 & - W_{MTOW} .Math. & 0 & 0 & 0 \\ - W_{MTOW} .Math. .Math. .Math. \sin .Math. .Math. & - W_{MTOW} .Math. .Math. .Math. \cos .Math. .Math. & 0 & 0 & - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S .Math. .Math. & 0 \end{matrix}] .Math. [\begin{matrix} a \\ b \\ c \\ C_{D, 0} \\ C_{L, 1} \\ K \end{matrix}] = m [\begin{matrix} a_{1}^{BFS} .Math. \cos .Math. .Math. + a_{3}^{BFS} .Math. \sin .Math. .Math. \\ a_{2}^{BFS} \\ - a_{1}^{BFS} .Math. \sin .Math. .Math. + a_{3}^{BFS} .Math. \cos .Math. .Math. \end{matrix}] & [24] \end{matrix}$

where:

a=C.sub.T cos cos

b=C.sub.T sin

c=C.sub.T sin cos [25]

[0058] The linear observation problem represented by expression [24] is underdetermined, as there are only 3 equations for 6 unknowns. Furthermore, while C.sub.D,0, C.sub.L,1 and K depend on the Mach number M and trim condition .sub.h, the unknowns a, b and c depend on the same variables as C.sub.T, i.e. {,,M,.sub.T}. The thrust model determination step 220 may be performed through flight testing when the flight is controlled so the mentioned variables are held constant, which can be achieved by performing thrust modelling maneuvers 222: [0059] and can be held almost constant by holding barometric altitude H through an altitude-on-bank control loop (as pressure varies relatively slowly with altitude, it is reserved the elevator to control speed through AOA); both positive and negative banks should be used to prevent asymmetry issues. Holding barometric altitude makes {dot over ()}{dot over ()}q0. [0060] The throttle level .sub.T is held constant [0061] For the selected throttle level, M is held constant through a speed-on-elevator control loop; the value of M is chosen so that level flight requires adopting a bank angle that leaves the altitude-on-bank loop room for altitude control through slightly increasing or decreasing bank angle as needed

[0062] Given the actual mass m there is a unique pair of values of a and trim condition .sub.h that make the flight condition described possible. Thus, in order to render the linear observation system of [24] determined or, preferably, overdetermined, k>2 different samples of AOA (ranging from low to high, but avoiding stall) and corresponding actual mass need to be considered, which add more observation equations without adding more unknowns. The idea is to perform a sequence of cases (ideally within the same flight test to avoid variations in the relationship between and ), for different values of {, , M, .sub.T}, which shall be repeated several times once the mass has experienced a significant change due to fuel consumption so the dataset collected ends up containing enough samples (i=1, . . . , k) of the same cases (j=1, . . . , l) for significantly different actual masses (and, thus, AOA values). In order to capture the variation of C.sub.T with and independently, the tasks should be repeated under different OAT (Outside Air Temperature) conditions with sufficient variation among each other covering the desired operational range (e.g. early in the morning vs. noon or selecting cold vs. hot days).

[00010] $\begin{matrix} H_{ij} = [\begin{matrix} W_{MTOW} .Math._{j} .Math. .Math. \cos .Math. .Math._{i} & - W_{MTOW} .Math._{j} .Math. \cos .Math. .Math._{i} & 0 & - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math._{j} .Math. M_{j}^{2} .Math. S & 0 & - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math._{j} .Math. M_{j}^{2} .Math. S .Math. .Math._{i}^{2} \\ 0 & 0 & - W_{MTOW} .Math._{j} & 0 & 0 & 0 \\ - W_{MTOW} .Math._{j} .Math. .Math. \sin .Math. .Math._{i} & - W_{MTOW} .Math. {.Math.}_{j} .Math. .Math. \cos .Math. .Math._{i} & 0 & 0 & - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math._{j} .Math. M_{j}^{2} .Math. S .Math. .Math._{i} & 0 \end{matrix}] & [26] \end{matrix}$

[0063] Observation Matrix for Sample i of Case j

[00011] $z_{j} = [\begin{matrix} a_{j} \\ b_{j} \\ c_{j} \\ C_{D, 0, j} \\ C_{L, 1, j} \\ K_{j} \end{matrix}]$

[0064] Vector of Unknowns for Case j

[00012] $O_{ij} = m_{i} [\begin{matrix} a_{1, ij}^{BFS} .Math. \cos .Math. .Math._{i} + a_{3, ij}^{BFS} .Math. \sin .Math. .Math._{i} \\ a_{1, ij}^{BFS} \\ - a_{1, ij}^{BFS} .Math. \sin .Math. .Math._{i} + a_{3, ij}^{BFS} .Math. \cos .Math. .Math._{i} \end{matrix}]$

[0065] Vector of Observations for Sample i of Case j

[0066] By combining all samples (for i=1, . . . , k) of case j, the overdetermined linear observation system below is obtained:

[00013] $\begin{matrix} H_{ij} = f (_{j},_{j}, M_{j},_{i}) & [29] \\ O_{ij} = f (m_{i},_{i}, a_{1, ij}^{BFS}, a_{1, ij}^{BFS}, a_{3, ij}^{BFS}) & [30] \\ H_{j} = [\begin{matrix} H_{1 .Math. j} \\ H_{2 .Math. j} \\ .Math. \\ H_{kj} \end{matrix}] & [31] \\ O_{j} = [\begin{matrix} O_{1 .Math. j} \\ O_{2 .Math. j} \\ .Math. \\ O_{kj} \end{matrix}] & [32] \\ H_{j} .Math. z_{j} = O_{j} & [33] \end{matrix}$

[0067] The best estimation of z.sub.j in the least-squares sense is given by the expression:

z.sub.1=(H.sub.j.sup.TH.sub.j).sup.1H.sub.j.sup.TO.sub.j for j=1, . . . ,l[34]

[0068] The LS solution obtained also provides a model of the aerodynamic drag and lift coefficients as a function of AOA and Mach valid for null AOS conditions.

[0069] To identify and (constant for all the cases) and C.sub.T(, , M, .sub.T) for each particular case further processing is required. In effect, expression [25] can then be used to obtain a linear observation system from [25], assuming that and are very small angles, the following approximation may be made:

a=C.sub.T

b=C.sub.T

C=C.sub.T[35]

and take logarithms in both sides of the equations:

ln a=ln C.sub.T

ln b=ln C.sub.T+ln

ln c=ln C.sub.T+ln [36]

[0070] Considering all the cases identified above, an overall linear observation system for and and C.sub.T,j for j=1, . . . , l can be composed as follows:

[00014] $\begin{matrix} H = [\begin{matrix} 0 & 0 & 1 & 0 & .Math. & 0 \\ 1 & 0 & 1 & 0 & .Math. & 0 \\ 0 & 1 & 1 & 0 & .Math. & 0 \\ 0 & 0 & 0 & 1 & .Math. & 0 \\ 1 & 0 & 0 & 1 & .Math. & 0 \\ 0 & 1 & 0 & 1 & .Math. & 0 \\ .Math. & .Math. & .Math. & .Math. & .Math. \\ 0 & 0 & 0 & 0 & .Math. & 1 \\ 1 & 0 & 0 & 0 & .Math. & 1 \\ 0 & 1 & 0 & 0 & .Math. & 1 \end{matrix}] .Math. .Math. \dim (H) = 3 .Math. l (2 + l) & [37] \\ O = [\begin{matrix} \ln .Math. .Math. a_{1} \\ \ln .Math. .Math. b_{1} \\ \ln .Math. .Math. c_{1} \\ \ln .Math. .Math. a_{2} \\ \ln .Math. .Math. b_{2} \\ \ln .Math. .Math. c_{2} \\ .Math. \\ \ln .Math. .Math. a_{l} \\ \ln .Math. .Math. b_{l} \\ \ln .Math. .Math. c_{l} \end{matrix}] .Math. .Math. \dim .Math. .Math. (O) = 3 .Math. l (2 + l) & [38] \\ z = [\begin{matrix} \ln .Math. .Math. .Math. \\ \ln .Math. .Math. \\ \ln .Math. .Math. C_{T, 1} \\ \ln .Math. .Math. C_{T, 2} \\ .Math. \\ \ln .Math. .Math. C_{T, l} \end{matrix}] .Math. .Math. \dim (z) = (2 + l) 1 & [39] \end{matrix}$

[0071] The LS solution is then obtained as:

z=(H.sup.TH).sup.1H.sup.TO[40]

which finally renders:

[00015] $.Math. = \exp (z [1])$ $= \exp (z [2])$ $Thrust .Math. .Math. orientation .Math. .Math. parameters$ $\begin{matrix} C_{T, 1} = C_{T} (_{1},_{1}, M_{1}, {.Math.}_{T, 1}) = \exp (z [2 + 1]) \\ C_{T, 2} = C_{T} (_{2},_{2}, M_{2}, {.Math.}_{T, 2}) = \exp (z [2 + 2]) \\ .Math. \\ C_{T, l} = C_{T} (_{l},_{l}, M_{l}, {.Math.}_{T, l}) = \exp (z [2 + l]) \end{matrix}$ $Thrust .Math. .Math. coefficients .Math. .Math. identified .Math. .Math. for .Math. .Math. all .Math. .Math. the .Math. .Math. test .Math. .Math. cases .Math. .Math. performed$

[0072] The thrust model determination step 220 described above has a limitation associated to the fact that all these test cases are performed at constant altitudes. In the experimental conditions defined, the thrust model derived cannot capture thrust levels below the minimum required to hold altitude in level flight with minimum actual mass at the speed that renders maximum aerodynamic efficiency (i.e. the minimum thrust to sustain level flight with minimum mass). Thus, no thrust information can be obtained for throttle levels that deliver a thrust lower than that one, for which it is needed an alternative approach that involves lower throttle levels.

[0073] A possible approach to identify the low throttle part of the thrust model may consist on performing decelerations at constant altitude to avoid variations of and . Starting from a flight condition at high speed (M) and high bank angle should give enough time to record speed variations in between the time when the engine transient (after having set and held a low throttle level) is over and the moment where level flight can no longer be sustained. Another approach may consist on performing descents holding low throttle level and, e.g speed (M), and record data at the altitudes of interest. In both cases, the thrust can be estimated through the following expression:

[00016] $W_{MTOW} .Math. [\begin{matrix} \cos .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. - \sin .Math. .Math. .Math. .Math. \sin .Math. .Math. .Math. \\ - \sin .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. \\ - (\sin .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. + \cos .Math. .Math. .Math. .Math. \sin .Math. .Math. .Math.) \end{matrix}] .Math. C_{T} = m [\begin{matrix} a_{1}^{BFS} .Math. \cos .Math. .Math. + a_{3}^{BFS} .Math. \sin .Math. .Math. \\ a_{2}^{BFS} \\ - a_{1}^{BFS} .Math. \sin .Math. .Math. + a_{3}^{BFS} .Math. \cos .Math. .Math. \end{matrix}] + \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S [\begin{matrix} (C_{D, 0} + K_{^{2}}) \\ 0 \\ C_{L, 1} .Math. \end{matrix}]$

with , , C.sub.D,0, K and C.sub.L,1 already known, which already makes expression [43] an overdetermined linear observation system for C.sub.T. In principle, no special provisions and or control for , , m or M are required as all these variables are observables, however, for practical reasons, it might be of interest to set up/control then and/or record data to fit the same cases identified for higher throttle levels.

[0074] In effect, to identify C.sub.T(.sub.j,.sub.j,M.sub.j,.sub.T,j) following the first approach, .sub.j and be set and held for the same cases as in high throttle and C.sub.T estimated when M matches the Mach M.sub.j associated to those cases. If the second approach is followed, M.sub.j can be set and held and then estimate C.sub.T when .sub.j matches the cases identified for high throttle levels. The dependency of thrust with is less than an issue for low level throttle. Notice that regardless the approach, this identification problem does only require controlling 2 control DOFs of the AV motion, namely, the throttle level and either altitude (in the first approach) or Mach (in the second approach), which allows using the remaining 3.sup.rd DOF to change AOA or bank angle to increase redundancy in the observation and, thus, achieve more robust estimations of C.sub.T. Alternatively, redundancy of observations in this problem can be, again, increased by repeating the test cases for different masses and corresponding AOAs, analogously to what has been done in the high throttle cases.

[0075] Expression [43] can be employed in general to further identify thrust model cases or refine them regardless the level of throttle, once , , C.sub.D,0, K and C.sub.L,1 are well known, as long as AOS is held null.

Step 230 for Determination of Aerodynamic Forces Model

[0076] Once the thrust model parameters have been identified, expression [1] can be used again, this time as a direct observer for the aerodynamic force coefficients as follows:

[00017] $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S [\begin{matrix} C_{D} \\ C_{Q} \\ C_{L} \end{matrix}] = [\begin{matrix} {Tt}_{1} \\ {Tt}_{2} \\ {Tt}_{3} \end{matrix}] - m [\begin{matrix} a_{1}^{BFS} .Math. \cos .Math. .Math. .Math. .Math. \cos .Math. .Math. + a_{2}^{BFS} .Math. \sin .Math. .Math. + a_{3}^{BFS} .Math. \sin .Math. .Math. \cos .Math. .Math. \\ - a_{1}^{BFS} .Math. \cos .Math. .Math. .Math. .Math. \sin .Math. .Math. + a_{2}^{BFS} .Math. \cos .Math. .Math. - a_{3}^{BFS} .Math. \sin .Math. .Math. \sin .Math. .Math. \\ - a_{1}^{BFS} .Math. \sin .Math. .Math. + a_{3}^{BFS} .Math. \cos .Math. .Math. \end{matrix}]$

with T given by expression [6], where the thrust model C.sub.T(,,M,.sub.T) is already known from previous step.

[0077] Expression [44] allows the direct estimation of the aerodynamic force coefficients C.sub.D, C.sub.Q and C.sub.L in terms of their dependency variables through flight testing, for which the respective domains have to be swept, which can be done manually and or with the help of control loops such as AOA-on-elevator and AOS-on-rudder and speed-on-elevator. If the mentioned control loops are available, an alternative approach to the identification of the aerodynamic force coefficients could take advantage of the two remaining control degrees of freedom (i.e. ailerons and throttle) and variables such as , M and to produce an overdetermined linear observation system, which would, expectedly, bring a best estimate where errors are further compensated.

[0078] Another validation that can be done is to check how well C.sub.D and C.sub.L fit the parabolic drag polar coefficients as a function of Mach identified in Step 220 for the case of null AOS.

[0079] Once the thrust model is known, expression [44] can still be used, if necessary, to estimate aerodynamic coefficients in the most general form:

C.sub.L=C.sub.L(,M,{dot over ({circumflex over ()})},{dot over ({circumflex over ()})},{circumflex over (p)},{circumflex over (q)},{circumflex over (r)},.sub.h,.sub.a,.sub.r,.sub.e)[45]

C.sub.Q=C.sub.Q(,M,{dot over ({circumflex over ()})},{dot over ({circumflex over ()})},{circumflex over (p)},{circumflex over (q)},{circumflex over (r)},.sub.h,.sub.a,.sub.r,.sub.e)[46]

C.sub.D=C.sub.D(,M,{dot over ({circumflex over ()})},{dot over ({circumflex over ()})},{circumflex over (p)},{circumflex over (q)},{circumflex over (r)},.sub.h,.sub.a,.sub.r,.sub.e)[47]

where:

[00018] $\hat{\dot{}} = \frac{\dot{} .Math. c_{w}}{2 .Math. .Math. v_{TAS}}$ $\hat{\dot{}} = \frac{\dot{} .Math. b_{w}}{2 .Math. .Math. v_{TAS}}$

[0080] However, for uncompressible aerodynamics and small AOS the main dependencies of the aerodynamics coefficients under the quasi-steady state and typical symmetry assumptions get simpler:

C.sub.L=C.sub.L(,{circumflex over (q)},.sub.h,.sub.e)[50]

C.sub.Q=C.sub.Q(,{circumflex over (p)},{circumflex over (r)},.sub.a,.sub.r)[51]

C.sub.D=C.sub.D(,{circumflex over (q)},h,.sub.e)[52]

Step 240 for Determination of Inertia Properties and Propulsive Moments

[0081] The 2.sup.nd Newton law applied to angular dynamics renders the following expression for the balance of moments:

[00019] $[\begin{matrix} R_{A} \\ P_{A} \\ Y_{A} \end{matrix}] = [\begin{matrix} I_{xx} & 0 & - I_{xz} \\ 0 & I_{yy} & 0 \\ - I_{xz} & 0 & I_{zz} \end{matrix}] [\begin{matrix} \dot{p} - D_{9} .Math. pq + D_{10} .Math. qr \\ \dot{q} + D_{11} (p^{2} - r^{2}) + D_{12} .Math. pr \\ \dot{r} - D_{13} .Math. pq + D_{14} .Math. qr \end{matrix}] - T [\begin{matrix} E_{1} \\ E_{2} \\ E_{3} \end{matrix}] - M_{T} [\begin{matrix} \cos .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. \\ - \sin .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. \\ - \sin .Math. .Math. .Math. \end{matrix}]$

Where:

[0082] [00020] $D_{9} = D_{7} .Math. I_{xz} - D_{3} .Math. D_{8}$ $D_{10} = D_{1} .Math. D_{7} + D_{8} .Math. I_{xz}$ $D_{11} = D_{6} .Math. I_{xz}$ $D_{12} = D_{2} .Math. D_{6}$ $D_{13} = D_{8} .Math. I_{xz} - D_{3} .Math. D_{5}$ $D_{14} = D_{1} .Math. D_{8} + D_{5} .Math. I_{xz}$ $D_{1} = I_{zz} - I_{yy}$ $D_{2} = I_{xx} - I_{zz}$ $D_{3} = I_{yy} - I_{xx}$ $D_{4} = \frac{I_{xx} .Math. I_{zz} - I_{xz}^{2}}{I_{yy}}$ $D_{5} = \frac{I_{xx}}{I_{yy} .Math. D_{4}} = \frac{I_{xx}}{I_{xx} .Math. I_{zz} - I_{xz}^{2}}$ $D_{6} = \frac{1}{I_{yy}}$ $D_{7} = \frac{I_{zz}}{I_{yy} .Math. D_{4}} = \frac{I_{zz}}{I_{xx} .Math. I_{zz} - I_{xz}^{2}}$ $D_{8} = \frac{I_{xz}}{I_{yy} .Math. D_{4}} = \frac{I_{xz}}{I_{xx} .Math. I_{zz} - I_{xz}^{2}}$ $I_{xx}, I_{yy}, I_{zz}, I_{xz} = f (m)$ $E_{1} = - y_{O_{T}}^{BFS} .Math. \sin .Math. .Math. .Math. + z_{O_{T}}^{BFS} .Math. \sin .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math.$ $E_{2} = x_{O_{T}}^{BFS} .Math. \sin .Math. .Math. .Math. + z_{O_{T}}^{BFS} .Math. \cos .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math.$ $E_{3} = - x_{O_{T}}^{BFS} .Math. \sin .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. - y_{O_{T}}^{BFS} .Math. \cos .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. .Math. [\begin{matrix} R_{T} \\ P_{T} \\ Y_{T} \end{matrix}] = T [\begin{matrix} E_{1} \\ E_{2} \\ E_{3} \end{matrix}] + M_{T} [\begin{matrix} \cos .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. \\ - \sin .Math. .Math. .Math. .Math. \cos .Math. .Math. .Math. \\ - \sin .Math. .Math. .Math. \end{matrix}]$

Propulsive or thrust-related moments

[0083] In the general case, aerodynamic moments can be expressed as:

R.sub.A=p.sub.0M.sup.2Sb.sub.wC.sub.R(,,M,{dot over ({circumflex over ()})},{dot over ({circumflex over ()})},.sub.h,.sub.a,.sub.r,.sub.e) Aerodynamic roll moment[61]

P.sub.A=p.sub.0M.sup.2Sb.sub.wC.sub.P(,,M,{dot over ({circumflex over ()})},{dot over ({circumflex over ()})},.sub.h,.sub.a,.sub.r,.sub.e) Aerodynamic pitch moment[62]

Y.sub.A=p.sub.0M.sup.2Sb.sub.wC.sub.Y(,,M,{dot over ({circumflex over ()})},{dot over ({circumflex over ()})},.sub.h,.sub.a,.sub.r,.sub.e) Aerodynamic yaw moment[61]

M.sub.T=m.sub.MTOWb.sub.w.sup.2N.sub.s.sup.2C.sub.M.sub.T Reactive torque due to rotating engine parts & propeller, where the rotating speed N.sub.s is assumed to be observable[64]

[0084] With the assumption that and are small angles:

[00021] $.Math. M_{T} [\begin{matrix} 1 \\ - \\ - .Math. \end{matrix}]$ $.Math. E_{1} = - y_{O_{T}}^{BFS} .Math. .Math. + z_{O_{T}}^{BFS} .Math.$ $.Math. E_{2} = x_{O_{T}}^{BFS} .Math. .Math. + z_{O_{T}}^{BFS}$ $.Math. E_{3} = - x_{O_{T}}^{BFS} .Math. - y_{O_{T}}^{BFS} .Math. [\begin{matrix} R_{A} \\ P_{A} \\ Y_{A} \end{matrix}] = [\begin{matrix} I_{xx} & 0 & - I_{xz} \\ 0 & I_{yy} & 0 \\ - I_{xz} & 0 & I_{zz} \end{matrix}] [\begin{matrix} \dot{p} - D_{9} .Math. pq + D_{10} .Math. qr \\ \dot{q} + D_{11} (p^{2} - r^{2}) + D_{12} .Math. pr \\ \dot{r} - D_{13} .Math. pq + D_{14} .Math. qr \end{matrix}] - T [\begin{matrix} E_{1} \\ E_{2} \\ E_{3} \end{matrix}] - M_{T} [\begin{matrix} 1 \\ - \\ - .Math. \end{matrix}]$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. {Sb}_{w} .Math. C_{R} (,, M, \hat{\dot{}}, \hat{\dot{}}, \hat{p}, \hat{q}, \hat{r}, {.Math.}_{h}, {.Math.}_{a}, {.Math.}_{r}, {.Math.}_{e}) = I_{xx} (m) .Math. (\dot{p} - D_{9} .Math. pq + D_{10} .Math. qr) - I_{xz} (m) .Math. (\dot{r} - D_{13} .Math. pq + D_{14} .Math. qr) - {TE}_{1} - M_{T}$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. {Sc}_{w} .Math. C_{P} (,, M, \hat{\dot{}}, \hat{\dot{}}, \hat{p}, \hat{q}, \hat{r}, {.Math.}_{h}, {.Math.}_{a}, {.Math.}_{r}, {.Math.}_{e}) = I_{yy} (m) .Math. (\dot{q} + D_{11} (p^{2} - r^{2}) + D_{12} .Math. pr) - {TE}_{2} - M_{T} .Math.$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. {Sb}_{w} .Math. C_{Y} (,, M, \hat{\dot{}}, \hat{\dot{}}, \hat{p}, \hat{q}, \hat{r}, {.Math.}_{h}, {.Math.}_{a}, {.Math.}_{r}, {.Math.}_{e}) = I_{xz} (m) .Math. (\dot{p} - D_{9} .Math. pq + D_{10} .Math. qr) + I_{zz} (m) .Math. (\dot{r} - D_{13} .Math. pq + D_{14} .Math. qr) - {TE}_{3} + M_{T} .Math. .Math.$

[0085] To isolate longitudinal dynamics, it is assumed that AOS-on-rudder and bank-on-ailerons control loops are active to maintain the condition:

{dot over ()}0[71]

pr0[72]

[0086] If a steady coordinated level flight is held at given altitude , Mach M.sub.0 and actual mass m.sub.0, which requires certain known AOA .sub.0 (q{dot over ()}0) and trim level .sub.h,0, then expression [69] further reduces to:

p.sub.0M.sub.0.sup.2Sc.sub.wC.sub.P,0=T.sub.0E.sub.2+M.sub.T,0=W.sub.MTOWC.sub.T,0E.sub.2+m.sub.MTOWb.sub.w.sup.2N.sub.s,0.sup.2C.sub.M.sub.T[73]

where:

C.sub.p,0=C.sub.P(.sub.0,0,M.sub.0,0,0,0,0,0,.sub.h,0,0,0,0)[74]

C.sub.T,0=C.sub.T,0(,,M.sub.0,.sub.T,0)[75]

[0087] If the elevator is actuated, an aerodynamic pitch moment appears, which breaks the balance and results in pitch acceleration; the expression that governs the longitudinal motion (after returning the elevator to the null position) reduces to:

p.sub.0M.sup.2Sc.sub.wC.sub.P(,0,M,{dot over ({circumflex over ()})},0,0,{circumflex over (q)},0,.sub.h,0,.sub.r,0)=I.sub.yy(m){dot over (q)}TE.sub.2+M.sub.T=I.sub.yy(m){dot over (q)}W.sub.MTOWC.sub.T(,,M,.sub.T)E.sub.2+m.sub.MTOWb.sub.w.sup.2N.sub.s.sup.2C.sub.M.sub.T[76]

[0088] Close enough to the balanced level flight condition and, as long as the flight is kept coordinated, .sub.e is held null and .sub.h,0 is held constant, the aerodynamic pitch moment and thrust coefficients can be approximated by a Taylor expansion in the form (first order approach):

C.sub.P(,0,M,{dot over ({circumflex over ()})},0,0,{circumflex over (q)},0,.sub.h,0,0,.sub.r,0)=C.sub.P,0+C.sub.P,+C.sub.P.sub.A.sub.,MM+C+C.sub.P,{dot over ({circumflex over ()})}{dot over ({circumflex over ()})}+C.sub.P,{circumflex over (q)}{circumflex over (q)}+C.sub.P,.sub.r.sub.r[77]

C.sub.T(,,M,.sub.T,0)=C.sub.T,0(,,M.sub.0,.sub.T,0)+C.sub.T,M(,,M.sub.0,.sub.T,0)M[78]

with:

[00022] $= -_{0}$ $.Math. .Math. M = M - M_{0}$ $C_{P, \hat{\dot{}}} .Math. \hat{\dot{}} = C_{P, \hat{\dot{}}} .Math. \frac{\dot{} .Math. c_{w}}{2 .Math. .Math. v_{TAS}} = C_{P, \hat{\dot{}}} .Math. \frac{c_{w}}{2 .Math. .Math. a_{0}} .Math. \frac{\dot{}}{M} = C_{P, \dot{}} .Math. \frac{\dot{}}{M}$ $C_{P, \dot{}} = C_{P, \hat{\dot{}}} .Math. \frac{c_{w}}{2 .Math. .Math. a_{0}}$ $C_{P, \hat{q}} .Math. \hat{q} = C_{P, \hat{q}} .Math. \frac{{qc}_{w}}{2 .Math. .Math. v_{TAS}} = C_{P, \hat{q}} .Math. \frac{c_{w}}{2 .Math. .Math. a_{0}} .Math. \frac{q}{M} = C_{P, q} .Math. \frac{q}{M}$ $C_{P, q} = C_{P, \hat{q}} .Math. \frac{c_{w}}{2 .Math. .Math. a_{0}}$

[0089] Assuming that the steady flight condition is perturbed with a longitudinal command (typically a so-called doublet) by means of the elevator and then return the elevator back to and hold it at the null position, without changing the throttle level (.sub.T,0 is held unchanged, thus, N.sub.s,0 does not change either), so the amplitude of the resulting oscillatory motion (the combination of the so-called short-term and the phugoid response modes) is small enough to allow approximating C.sub.P and C.sub.T by, respectively, expressions [77] and [78]:

[00023] $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. {Sc}_{w} (C_{P, 0} + C_{P,} .Math. + C_{P_{A}, M} .Math. .Math. .Math. M + C_{P, \dot{}} .Math. \frac{\dot{}}{M} + C_{P, q} .Math. \frac{q}{M} + C_{P, {.Math.}_{r}} .Math. {.Math.}_{r}) = I_{yy} (m) .Math. \dot{q} - W_{MTOW} .Math. (C_{T, 0} + C_{T, M} .Math. .Math. .Math. M) .Math. E_{2} + m_{MTOW} .Math. b_{w}^{2} .Math. N_{s, 0}^{2} .Math. C_{M_{T}} .Math.$ $M^{2} = {(M_{0} + .Math. .Math. M)}^{2} = M_{0}^{2} + 2 .Math. .Math. M_{0} .Math. .Math. .Math. M + .Math. .Math. M^{2} = M_{0}^{2} [1 + 2 .Math. \frac{.Math. .Math. M}{M_{0}} + {(\frac{.Math. .Math. M}{M_{0}})}^{2}]$ $.Math. \frac{.Math. .Math. M}{M_{0}} 1$ $.Math. M^{2} M_{0}^{2} (1 + 2 .Math. \frac{.Math. .Math. M}{M_{0}})$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. (M_{0}^{2} + 2 .Math. .Math. M_{0} .Math. .Math. .Math. M) .Math. {Sc}_{w} .Math. C_{P, 0} + \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. {Sc}_{w} (C_{P,} .Math. + C_{P_{A}, M} .Math. .Math. .Math. M + C_{P, \dot{}} .Math. \frac{\dot{}}{M} + C_{P, q} .Math. \frac{q}{M} + C_{P, {.Math.}_{r}} .Math. {.Math.}_{r}) = I_{yy} (m) .Math. \dot{q} - W_{MTOW} .Math. (C_{T, 0} + C_{T, M} .Math. .Math. .Math. M) .Math. E_{2} + m_{MTOW} .Math. b_{w}^{2} .Math. N_{s, 0}^{2} .Math. C_{M_{T}} .Math.$

[0090] Bearing in mind expression [73], the longitudinal moments balance yields:

[00024] $.Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w} .Math. M_{0} .Math. .Math. .Math. {MC}_{P, 0} + \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. {Sc}_{w} (C_{P,} .Math. + C_{P_{A}, M} .Math. .Math. .Math. M + C_{P, \dot{}} .Math. \frac{\dot{}}{M} + C_{P, q} .Math. \frac{q}{M} + C_{P, {.Math.}_{r}} .Math. {.Math.}_{r}) = I_{yy} (m) .Math. \dot{q} - W_{MTOW} .Math. .Math. .Math. C_{T, M} .Math. .Math. .Math. {ME}_{2}$ $.Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w} .Math. M_{0} .Math. C_{P, 0} .Math. .Math. .Math. M + \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w} (M^{2} .Math. C_{P,} .Math. + M^{2} .Math. C_{P_{A}, M} .Math. .Math. .Math. M + {MC}_{P, \dot{}} .Math. \dot{} + {MC}_{P, q} .Math. q + M^{2} .Math. C_{P, {.Math.}_{r}} .Math. {.Math.}_{r}) = I_{yy} (m) .Math. \dot{q} - W_{MTOW} .Math. .Math. .Math. C_{T, M} .Math. .Math. .Math. {ME}_{2}$ $M^{2} .Math. .Math. .Math. M = {(M_{0} + .Math. .Math. M)}^{2} .Math. .Math. .Math. M = (M_{0}^{2} + 2 .Math. .Math. M_{0} .Math. .Math. .Math. M + .Math. .Math. M^{2}) .Math. .Math. .Math. M = M_{0}^{2} [1 + 2 .Math. \frac{.Math. .Math. M}{M_{0}} + {(\frac{.Math. .Math. M}{M_{0}})}^{2}] .Math. .Math. .Math. M M_{0}^{2} (1 + 2 .Math. \frac{.Math. .Math. M}{M_{0}}) .Math. .Math. .Math. M$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w} (2 .Math. .Math. M_{0} .Math. C_{P, 0} + \frac{W_{MTOW} .Math. C_{T, M} .Math. E_{2}}{\frac{1}{2} .Math. .Math. .Math. P_{0} .Math. {Sc}_{w}}) .Math. .Math. .Math. M + \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w} [M^{2} .Math. C_{P,} .Math. + M_{0}^{2} (1 + 2 .Math. \frac{.Math. .Math. M}{M_{0}}) .Math. C_{P_{A}, M} .Math. .Math. .Math. M + {MC}_{P, \dot{}} .Math. \dot{} + {MC}_{P, q} .Math. q + M^{2} .Math. C_{P, {.Math.}_{r}} .Math. {.Math.}_{r}] = I_{yy} (m) .Math. \dot{q}$ $.Math. C_{P, M} = 2 .Math. \frac{C_{P, 0}}{M_{0}} + C_{P_{A}, M} + \frac{W_{MTOW} .Math. C_{T, M} .Math. E_{2}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. {Sc}_{w} .Math. M_{0}}$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w} (M_{0}^{2} .Math. C_{P, M} .Math. .Math. .Math. M + M^{2} .Math. C_{P,} .Math. + 2 .Math. .Math. M_{0} .Math. C_{P_{A}, M} .Math. .Math. .Math. M^{2} + {MC}_{P, \dot{}} .Math. \dot{} + {MC}_{P, q} .Math. q + M^{2} .Math. C_{P, {.Math.}_{r}} .Math. {.Math.}_{r}) = I_{yy} (m) .Math. \dot{q}$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w} [M_{0}^{2} .Math. C_{P, M} .Math. .Math. .Math. M + M_{0}^{2} (1 + 2 .Math. \frac{.Math. .Math. M}{M_{0}}) .Math. C_{P,} .Math. + 2 .Math. .Math. M_{0} .Math. C_{P_{A}, M} .Math. .Math. .Math. M^{2} + (M_{0} + .Math. .Math. M) .Math. C_{P, \dot{}} .Math. \dot{} + (M_{0} + .Math. .Math. M) .Math. C_{P, q} .Math. q + M_{0}^{2} (1 + 2 .Math. \frac{.Math. .Math. M}{M_{0}}) .Math. C_{P, {.Math.}_{r}} .Math. {.Math.}_{r}] = I_{yy} (m) .Math. \dot{q}$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w} [M_{0}^{2} .Math. C_{P, M} .Math. .Math. .Math. M + M_{0}^{2} .Math. C_{P,} .Math. + 2 .Math. .Math. M_{0} .Math. C_{P_{A}, M} .Math. .Math. .Math. M + 2 .Math. .Math. M_{0} .Math. C_{P_{A}, M} .Math. .Math. .Math. M^{2} + M_{0} .Math. C_{P, \dot{}} .Math. \dot{} + C_{P, \dot{}} .Math. \dot{} .Math. .Math. .Math. M + M_{0} .Math. C_{P, q} .Math. q + C_{P, q} .Math. q .Math. .Math. .Math. .Math. M + M_{0}^{2} .Math. C_{P, {.Math.}_{r}} .Math. {.Math.}_{r} + 2 .Math. .Math. M_{0} .Math. C_{P, {.Math.}_{r}} .Math. .Math. .Math. M] = I_{yy} (m) .Math. \dot{q}$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w} .Math. M_{0}^{2} [C_{P,} .Math. + C_{P, M} .Math. .Math. .Math. M + \frac{C_{P, \dot{}}}{M_{0}} .Math. \dot{} + \frac{C_{P, q}}{M_{0}} .Math. q + C_{P, {.Math.}_{r}} .Math. {.Math.}_{r} + (2 .Math. \frac{C_{P,}}{M_{0}} .Math. + 2 .Math. \frac{C_{P_{A}, M}}{M_{0}} .Math. .Math. .Math. M + \frac{C_{P, \dot{}}}{M_{0}^{2}} .Math. \dot{} + \frac{C_{P, q}}{M_{0}^{2}} .Math. q + 2 .Math. \frac{C_{P, {.Math.}_{r}}}{M_{0}} .Math. {.Math.}_{r}) .Math. .Math. .Math. M] = I_{yy} (m) .Math. \dot{q}$

[0091] Expression [93] is a linear observer for the unknowns (C.sub.P,M, C.sub.P,, C.sub.P.sub.A.sub.,M, C.sub.P,{dot over ()}, C.sub.P,q, C.sub.P,.sub.r) and I.sub.yy(m) valid for any Mach and actual mass. However, pitch acceleration {dot over (q)} and the derivative of the AOA, a are not typically native observables, which means numerical derivation should be required. To avoid the deriving noisy signals, a better idea is to integrate the known signals in [93] over time.

[0092] One possibility consists of neglecting the influence of a, which yields:

p.sub.0Sc.sub.w(M.sub.0.sup.2C.sub.P,MM+M.sup.2C.sub.P,+2M.sub.0C.sub.P.sub.A.sub.,MM.sup.2+MC.sub.P,qq+M.sup.2C.sub.P,.sub.r.sub.r)=I.sub.yy(m){dot over (q)}[97]

f.sup.1=Mdt

f.sup.2=M.sup.2dt

f.sup.3=M.sup.2dt

f.sup.4=Mqdt

f.sup.r=M.sup.2.sub.rdt[98]

p.sub.0Sc.sub.w(M.sub.0.sup.2C.sub.P,Mf.sup.1+C.sub.P,f.sup.2+2M.sub.0C.sub.P.sub.A.sub.,Mf.sup.3+C.sub.P,qf.sup.4+C.sub.P,.sub.rf.sup.r)=I.sub.yy(m)q+B[99]

[0093] Expression [99] is a linear observer for the unknowns {C.sub.P,M, C.sub.P,, C.sub.P.sub.A.sub.,M, C.sub.P,q, C.sub.P,.sub.r, B} and I.sub.yy(m) valid for any Mach and actual mass, where B is an integration constant.

[0094] A second possibility consists on computing M{dot over ()} in [93] as:

M{dot over ()}=M{dot over (M)}[Expression for {dot over (M)} TBD][100]

[0095] The third possibility consists of eliminating the dependency with {dot over ()} in [93] through, substituting expression [286] obtained from the linearization of the equations of linear motion, i.e.:

[00025] $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w} [M_{0}^{2} .Math. C_{P, M} .Math. .Math. .Math. M + M^{2} .Math. C_{P,} .Math. + 2 .Math. .Math. M_{0} .Math. C_{P_{A}, M} .Math. .Math. .Math. M^{2} + \frac{C_{P, \dot{}}}{F_{4}} .Math. (M_{0}^{2} .Math. F_{1} .Math. .Math. .Math. M + M^{2} .Math. F_{2} .Math. + {MF}_{3} .Math. q) + {MC}_{P, q} .Math. q + M^{2} .Math. C_{P, {.Math.}_{r}} .Math. {.Math.}_{r}] = I_{yy} (m) .Math. \dot{q} + m_{0} .Math. c_{w} .Math. \frac{C_{P, \dot{}}}{F_{4}} .Math. (F_{5} .Math. .Math. .Math. a_{1}^{BFS} + F_{6} .Math. .Math. .Math. a_{3}^{BFS})$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w} [M_{0}^{2} (C_{P, M} + C_{P, \dot{}} .Math. \frac{F_{1}}{F_{4}}) .Math. .Math. .Math. M + (C_{P,} + C_{P, \dot{}} .Math. \frac{F_{2}}{F_{4}}) .Math. M^{2} .Math. + 2 .Math. .Math. M_{0} .Math. C_{P_{A}, M} .Math. .Math. .Math. M^{2} + (C_{P, q} + C_{P, \dot{}} .Math. \frac{F_{3}}{F_{4}}) .Math. Mq + C_{P, {.Math.}_{r}} .Math. M^{2} .Math. {.Math.}_{r}] = I_{yy} (m) .Math. \dot{q} + m_{0} (C_{P, \dot{}} .Math. \frac{F_{5}}{F_{4}} .Math. c_{w} .Math. .Math. .Math. a_{1}^{BFS} + C_{P, \dot{}} .Math. \frac{F_{6}}{F_{4}} .Math. c_{w} .Math. .Math. .Math. a_{3}^{BFS})$ $.Math. G_{1} = C_{P, M} + C_{P, \dot{}} .Math. \frac{F_{1}}{F_{4}}$ $.Math. G_{2} = C_{P,} + C_{P, \dot{}} .Math. \frac{F_{2}}{F_{4}}$ $.Math. G_{3} = 2 .Math. .Math. C_{P_{A}, M}$ $.Math. G_{4} = C_{P, q} + C_{P, \dot{}} .Math. \frac{F_{3}}{F_{4}}$ $.Math. G_{5} = C_{P, \dot{}} .Math. \frac{F_{5}}{F_{4}} .Math. c_{w}$ $.Math. G_{6} = C_{P, \dot{}} .Math. \frac{F_{6}}{F_{4}} .Math. c_{w}$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w} [M_{0}^{2} .Math. G_{1} .Math. .Math. .Math. M + G_{2} .Math. M^{2} .Math. + M_{0} .Math. G_{3} .Math. .Math. .Math. M^{2} + G_{4} .Math. Mq + C_{P, {.Math.}_{r}} .Math. M^{2} .Math. {.Math.}_{r}] = I_{yy} (m) .Math. \dot{q} + m_{0} (G_{5} .Math. .Math. .Math. a_{1}^{BFS} + G_{6} .Math. .Math. .Math. a_{3}^{BFS})$ $.Math. f^{5} = .Math. .Math. a_{1}^{BFS} .Math. dt$ $.Math. f^{6} = .Math. .Math. a_{3}^{BFS} .Math. dt$

which allows rewriting expression [93] as:

p.sub.0Sc.sub.w[M.sub.0.sup.2G.sub.1f.sup.1+G.sub.2f.sup.2+M.sub.0G.sub.3f.sup.3+G.sub.4f.sup.4+C.sub.P,.sub.rf.sup.r]=I.sub.yy(m)q+m.sub.0(G.sub.5f.sup.5+G.sub.6f.sup.6)+B[106]

where B is an integration constant

[0096] Expression [106] is a linear observer for the unknowns {G.sub.1, G.sub.2, G.sub.3, G.sub.4, C.sub.P,.sub.r, G.sub.5, G.sub.6, B} and I.sub.yy(m) valid for any Mach and actual mass. Thus by performing maneuvers for different cases of Mach at the trim condition and actual mass, an overdetermined linear observation system is obtained that allows estimating the mentioned unknowns. In effect, a case may be considered for the identification of I.sub.yy(m) for the combination {m.sub.i, M.sub.j}, i.e., at the trim condition at Mach M.sub.j (j=1, . . . , l.sub.i) with actual mass m.sub.i (i=1, . . . , k). For each case of M.sub.j, n.sub.j samples of the response to a lateral perturbation are recorded. Thus:

[00026] $H_{ij} = [\begin{matrix} \begin{matrix} - \frac{q_{j .Math. .Math. 1}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} & M_{j}^{2} .Math. f_{j .Math. .Math. 1}^{1} & f_{j .Math. .Math. 1}^{2} & M_{j} .Math. f_{j .Math. .Math. 1}^{3} & f_{j .Math. .Math. 1}^{4} & f_{j .Math. .Math. 1}^{r} & - \frac{m_{i} .Math. f_{j .Math. .Math. 1}^{5}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} & - \frac{m_{i} .Math. f_{j .Math. .Math. 1}^{6}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} & - \frac{1}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} \end{matrix} \\ \begin{matrix} - \frac{q_{j .Math. .Math. 2}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} & M_{j}^{2} .Math. f_{j .Math. .Math. 2}^{1} & f_{j .Math. .Math. 2}^{2} & M_{j} .Math. f_{j .Math. .Math. 2}^{3} & f_{j .Math. .Math. 2}^{4} & f_{j .Math. .Math. 2}^{r} & - \frac{m_{i} .Math. f_{j .Math. .Math. 2}^{5}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} & - \frac{m_{i} .Math. f_{j .Math. .Math. 2}^{6}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} & - \frac{1}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} \\ .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. \end{matrix} \\ \begin{matrix} \begin{matrix} - \frac{q_{{jn}_{j}}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} & M_{j}^{2} .Math. f_{{jn}_{j}}^{1} & f_{{jn}_{j}}^{2} & M_{j} .Math. f_{{jn}_{j}}^{3} & f_{{jn}_{j}}^{4} & f_{{jn}_{j}}^{r} & - \frac{m_{i} .Math. f_{{jn}_{j}}^{5}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} & - \frac{m_{i} .Math. f_{{jn}_{j}}^{6}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} & - \frac{1}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} \end{matrix} \end{matrix} \end{matrix}] .Math. .Math. z_{ij} = [\begin{matrix} I_{yy} (m_{i}) \\ G_{1} (_{ij}, M_{j}, {.Math.}_{h, ij}) \\ G_{2} (_{ij}, M_{j}, {.Math.}_{h, ij}) \\ G_{3} (_{ij}, M_{j}, {.Math.}_{h, ij}) \\ G_{4} (_{ij}, M_{j}, {.Math.}_{h, ij}) \\ C_{P, {.Math.}_{r}} (_{ij}, M_{j}, {.Math.}_{h, ij}) \\ G_{5} (_{ij}, M_{j}, {.Math.}_{h, ij}) \\ G_{6} (_{ij}, M_{j}, {.Math.}_{h, ij}) \\ B_{ij} \end{matrix}] .Math. .Math. O_{ij} = [\begin{matrix} 0 \\ 0 \\ .Math. \\ 0 \end{matrix}]$

[0097] For the sake of robustness, I.sub.yy(m.sub.i) is obtained considering all the l.sub.i samples as follows:

[00027] $F_{ij} = [\begin{matrix} \begin{matrix} M_{j}^{2} .Math. f_{j .Math. .Math. 1}^{1} & f_{j .Math. .Math. 1}^{2} & M_{j} .Math. f_{j .Math. .Math. 1}^{3} & f_{j .Math. .Math. 1}^{4} & f_{j .Math. .Math. 1}^{r} & - \frac{m_{i} .Math. f_{j .Math. .Math. 1}^{5}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} & - \frac{m_{i} .Math. f_{j .Math. .Math. 1}^{6}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} & - \frac{1}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} \end{matrix} \\ \begin{matrix} M_{j}^{2} .Math. f_{j .Math. .Math. 2}^{1} & f_{j .Math. .Math. 2}^{2} & M_{j} .Math. f_{j .Math. .Math. 2}^{3} & f_{j .Math. .Math. 2}^{4} & f_{j .Math. .Math. 2}^{r} & - \frac{m_{i} .Math. f_{j .Math. .Math. 2}^{5}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} & - \frac{m_{i} .Math. f_{j .Math. .Math. 2}^{6}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} & - \frac{1}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} \\ .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. \end{matrix} \\ \begin{matrix} M_{j}^{2} .Math. f_{{jn}_{j}}^{1} & f_{{jn}_{j}}^{2} & M_{j} .Math. f_{{jn}_{j}}^{3} & f_{{jn}_{j}}^{4} & f_{{jn}_{j}}^{r} & - \frac{m_{i} .Math. f_{{jn}_{j}}^{5}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} & - \frac{m_{i} .Math. f_{{jn}_{j}}^{6}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} & - \frac{1}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} \end{matrix} \end{matrix}] .Math. .Math. G_{ij} = - \frac{1}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} [\begin{matrix} q_{j .Math. .Math. 1} \\ q_{j .Math. .Math. 2} \\ .Math. \\ q_{{jn}_{j}} \end{matrix}] .Math. .Math. C_{P, ij} = [\begin{matrix} G_{1} (_{ij}, M_{j}, {.Math.}_{h, ij}) \\ G_{2} (_{ij}, M_{j}, {.Math.}_{h, ij}) \\ G_{3} (_{ij}, M_{j}, {.Math.}_{h, ij}) \\ G_{4} (_{ij}, M_{j}, {.Math.}_{h, ij}) \\ C_{P, {.Math.}_{r}} (_{ij}, M_{j}, {.Math.}_{h, ij}) \\ G_{5} (_{ij}, M_{j}, {.Math.}_{h, ij}) \\ G_{6} (_{ij}, M_{j}, {.Math.}_{h, ij}) \\ B_{ij} \end{matrix}] .Math. .Math. H_{i} = [\begin{matrix} G_{i .Math. .Math. 1} & F_{i .Math. .Math. 1} & 0 & .Math. & 0 \\ G_{i .Math. .Math. 2} & 0 & F_{i .Math. .Math. 2} & .Math. & 0 \\ .Math. & .Math. & .Math. & .Math. \\ G_{{il}_{i}} & 0 & 0 & .Math. & F_{{il}_{i}} \end{matrix}] .Math. .Math. \dim (H_{i}) = ({.Math.}_{j = 1}^{l_{i}} .Math. .Math. n_{j}) (1 + 8 .Math. .Math. l_{i}) .Math. .Math. z_{i} = [\begin{matrix} I_{yy} (m_{i}) \\ C_{P, i .Math. .Math. 1} \\ C_{P, i .Math. .Math. 2} \\ .Math. \\ C_{P, {il}_{i}} \end{matrix}] .Math. .Math. \dim (z_{i}) = (1 + 8 .Math. .Math. l_{i}) 1 .Math. .Math. O_{i} = [\begin{matrix} 0 \\ 0 \\ .Math. \\ 0 \end{matrix}] .Math. .Math. \dim (O_{i}) = ({.Math.}_{j = 1}^{l_{i}} .Math. .Math. n_{j}) 1 .Math. .Math. z_{i} = {(H_{i}^{T} .Math. H_{i})}^{- 1} .Math. H_{i}^{T} .Math. O_{i} .Math. .Math. I_{yy} (m_{i}) = z_{i} [1]$

[0098] Repeating the maneuvers for different values of the actual mass would render the variation of I.sub.yy with m as I.sub.yy=f(m.sub.i) for i=1, . . . , k

[0099] Now, to isolate lateral dynamics, it is assumed that a new pitch rate-on-elevator control loop is active to ensure q0. If a steady level flight is held at given Mach M and actual mass m, which requires certain known AOA .sub.0 (q{dot over ()}0) and trim level .sub.h,0, then expressions [68] and [69] further reduce to:

p.sub.0M.sup.2Sb.sub.wC.sub.R,0=T.sub.0E.sub.1M.sub.T,0=W.sub.MTOWC.sub.T,0E.sub.1m.sub.MTOWb.sub.wN.sub.s,0.sup.2C.sub.M.sub.T[118]

p.sub.0M.sup.2Sb.sub.wC.sub.Y,0=T.sub.0E.sub.3+M.sub.T,0=W.sub.MTOWC.sub.T,0E.sub.3+m.sub.MTOWb.sub.w.sup.2N.sub.s,0.sup.2C.sub.M.sub.T [119]

where:

C.sub.R,0=C.sub.R(.sub.0,0,M.sub.0,0,0,0,0,0,.sub.h,0,0,0,0)[120]

C.sub.Y,0=C.sub.Y(.sub.0,0,M.sub.0,0,0,0,0,0,.sub.h,0,0,0,0)[121]

[0100] If the ailerons are actuated, an aerodynamic roll moment appears, which breaks the balance and results in roll acceleration; expression [68] reduces in that case to:

[00028] $\begin{matrix} \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. {Sb}_{w} .Math. C_{R} (,, M, 0, \hat{\dot{}}, \hat{p}, 0, \hat{r}, {.Math.}_{h, 0}, 0, 0, 0) = I_{xx} (m) .Math. \dot{p} - I_{xz} (m) .Math. \dot{r} - {TE}_{1} - M_{T} == I_{xx} (m) .Math. \dot{p} - I_{xz} (m) .Math. \dot{r} - W_{MTOW} .Math. .Math. .Math. C_{T} (,, M, {.Math.}_{T}) .Math. E_{1} - m_{MTOW} .Math. b_{w}^{2} .Math. N_{s}^{2} .Math. C_{M_{T}} & [122] \end{matrix}$

[0101] Close enough to the balanced flight condition and, as long as .sub.a and .sub.r are held null and, .sub.h,0 is held constant, the aerodynamic roll moment coefficient can be approximated by a Taylor expansion in the form:

C.sub.R(,,M,0,{dot over ({circumflex over ()})},{circumflex over (p)},0,{circumflex over (r)},.sub.h,0,0,0,0)=C.sub.R,0+C.sub.R,+C.sub.R,+C.sub.R.sub.A.sub.,+C.sub.R.sub.A.sub.,MM+C.sub.R,{dot over ({circumflex over ()})}{dot over ({circumflex over ()})}+C.sub.R,{circumflex over (p)}{circumflex over (p)}+C.sub.R,{circumflex over (f)}{circumflex over (r)}[123]

with:

[00029] $C_{R, \hat{\dot{}}} .Math. \hat{\dot{}} = C_{R, \hat{\dot{}}} .Math. \frac{\dot{} .Math. b_{w}}{2 .Math. .Math. v_{TAS}} = C_{R, \hat{\dot{}}} .Math. \frac{b_{w}}{2 .Math. .Math. a_{0}} .Math. \frac{\dot{}}{M} = C_{R, \dot{}} .Math. \frac{\dot{}}{M}$ $C_{R, \dot{}} .Math. C_{R, \hat{\dot{}}} .Math. \frac{b_{w}}{2 .Math. .Math. a_{0}}$ $C_{R, \hat{p}} .Math. \hat{p} = C_{R, \hat{p}} .Math. \frac{{pb}_{w}}{2 .Math. .Math. v_{TAS}} = C_{R, \hat{p}} .Math. \frac{b_{w}}{2 .Math. .Math. a_{0}} .Math. \frac{p}{M} = C_{R, p} .Math. \frac{p}{M}$ $C_{P, p} = C_{R, \hat{p}} .Math. \frac{b_{w}}{2 .Math. .Math. a_{0}}$ $C_{R, \hat{r}} .Math. \hat{r} = C_{R, \hat{r}} .Math. \frac{{rb}_{w}}{2 .Math. .Math. v_{TAS}} = C_{R, \hat{r}} .Math. \frac{b_{w}}{2 .Math. .Math. a_{0}} .Math. \frac{r}{M} = C_{R, r} .Math. \frac{r}{M}$ $C_{P, r} = C_{R, \hat{r}} .Math. \frac{b_{w}}{2 .Math. .Math. a_{0}}$

[0102] Analogously, if the rudder is actuated from the balanced condition, an aerodynamic yaw moment appears, which breaks the balance and results in yaw acceleration; expression [70] reduces in that case to:

[00030] $\begin{matrix} \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. {Sb}_{w} .Math. C_{Y} (,, M, 0, \hat{\dot{}}, \hat{p}, 0, \hat{r}, {.Math.}_{h, 0}, 0, 0, 0) = - I_{xz} (m) .Math. \dot{p} + I_{zz} (m) .Math. \dot{r} - {TE}_{3} + M_{T} .Math. .Math. = - I_{xz} (m) .Math. \dot{p} + I_{zz} (m) .Math. \dot{r} - W_{MTOW} .Math. .Math. .Math. C_{T} (,, M, {.Math.}_{T}) .Math. E_{3} + m_{MTOW} .Math. b_{w}^{2} .Math. N_{s}^{2} .Math. C_{M_{T}} .Math. .Math. & [127] \end{matrix}$

[0103] Again, close enough to the balanced flight condition and, as long as .sub.a and .sub.r are held null and, .sub.h,0 is held constant, the aerodynamic yaw moment coefficient can be approximated by a Taylor expansion in the form:

C.sub.Y(,{circumflex over (p)},{circumflex over (r)},.sub.a,.sub.r)=C.sub.Y,0+C.sub.Y,+C.sub.Y,+C.sub.Y.sub.A.sub.,MM+C.sub.Y,{dot over ({circumflex over ()})}{dot over ({circumflex over ()})}+C.sub.Y,{circumflex over (p)}{circumflex over (p)}+C.sub.Y,{circumflex over (r)}{circumflex over (r)}[128]

with:

[00031] $C_{Y, \hat{\dot{}}} .Math. \hat{\dot{}} = C_{Y, \hat{\dot{}}} .Math. \frac{\dot{} .Math. b_{w}}{2 .Math. .Math. v_{TAS}} = C_{Y, \hat{\dot{}}} .Math. \frac{b_{w}}{2 .Math. .Math. a_{0}} .Math. \frac{\dot{}}{M} = C_{Y, \dot{}} .Math. \frac{\dot{}}{M}$ $C_{Y, \dot{}} .Math. C_{Y, \hat{\dot{}}} .Math. \frac{b_{w}}{2 .Math. .Math. a_{0}}$ $C_{Y, \hat{p}} .Math. \hat{p} = C_{Y, \hat{p}} .Math. \frac{{pb}_{w}}{2 .Math. .Math. v_{TAS}} = C_{Y, \hat{p}} .Math. \frac{b_{w}}{2 .Math. .Math. a_{0}} .Math. \frac{p}{M} = C_{Y, p} .Math. \frac{p}{M}$ $C_{Y, p} = C_{Y, \hat{p}} .Math. \frac{b_{w}}{2 .Math. .Math. a_{0}}$ $C_{Y, \hat{r}} .Math. \hat{r} = C_{Y, \hat{r}} .Math. \frac{{rb}_{w}}{2 .Math. .Math. v_{TAS}} = C_{Y, \hat{r}} .Math. \frac{b_{w}}{2 .Math. .Math. a_{0}} .Math. \frac{r}{M} = C_{Y, r} .Math. \frac{r}{M}$ $C_{Y, r} = C_{Y, \hat{r}} .Math. \frac{b_{w}}{2 .Math. .Math. a_{0}}$

[0104] In both cases:

C.sub.T(,,M,.sub.T,0)=C.sub.T,0(,,M.sub.0,.sub.T,0)+C.sub.T,M(,,M.sub.0,.sub.T,0)M[132]

[0105] Expressions [122] and [127] along with the respective coefficients from [123], [128] and [132] govern the lateral-directional motion in the cases described, regardless that the perturbation that triggers such motion comes from ailerons or rudder actuations as rolling and yawing motions are coupled.

[0106] It is assumed that perturbation of the steady flight condition with a lateral-directional command (typically a so-called doublet) by means of either ailerons or rudder and then return the given control surface back to and hold it at the null position so the amplitude of the resulting oscillatory motion (in general, a combination of the so-called roll, Dutch roll and spiral response modes) is small enough to allow approximating C.sub.R and C.sub.y by, respectively expressions [123] and [128]:

[00032] $\begin{matrix} \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. {Sb}_{w} (C_{R, 0} + C_{R,} .Math. .Math. .Math. + C_{R,} .Math. + C_{R_{A}, M} .Math. .Math. .Math. M + C_{R, \dot{}} .Math. \frac{\dot{}}{M} + C_{R, p} .Math. \frac{p}{M} + C_{R, r} .Math. \frac{r}{M}) = I_{xx} (m) .Math. \dot{p} - I_{xz} (m) .Math. \dot{r} - W_{MTOW} .Math. (C_{T, 0} + C_{T, M} .Math. .Math. .Math. M) .Math. E_{1} - m_{MTOW} .Math. b_{w}^{2} .Math. N_{s, 0}^{2} .Math. C_{M_{T}} & [133] \\ \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. {Sb}_{w} (C_{Y, 0} + C_{Y,} .Math. .Math. .Math. + C_{Y,} .Math. + C_{Y_{A}, M} .Math. .Math. .Math. M + C_{Y, \dot{}} .Math. \frac{\dot{}}{M} + C_{Y, p} .Math. \frac{p}{M} + C_{Y, r} .Math. \frac{r}{M}) = - I_{xz} (m) .Math. \dot{p} + I_{zz} (m) .Math. \dot{r} - W_{MTOW} .Math. (C_{T, 0} + C_{T, M} .Math. .Math. .Math. M) .Math. E_{3} + m_{MTOW} .Math. b_{w}^{2} .Math. N_{s, 0}^{2} .Math. C_{M_{T}} .Math. .Math. & [134] \end{matrix}$

[0107] If expressions [84], [85] and [86] are taken into account:

[00033] $\begin{matrix} \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. (M_{0}^{2} + 2 .Math. M_{0} .Math. .Math. .Math. M) .Math. {Sb}_{w} .Math. C_{R, 0} .Math. + \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. {Sb}_{w} (C_{R,} .Math. .Math. .Math. + C_{R,} .Math. + C_{R_{A}, M} .Math. .Math. .Math. M + C_{R, \dot{}} .Math. \frac{\dot{}}{M} + C_{R, p} .Math. \frac{p}{M} + C_{R, r} .Math. \frac{r}{M}) = I_{xx} (m) .Math. \dot{p} - I_{xz} (m) .Math. \dot{r} - W_{MTOW} .Math. .Math. .Math. (C_{T, 0} + C_{T, M} .Math. .Math. .Math. M) .Math. E_{1} - m_{MTOW} .Math. b_{w}^{2} .Math. N_{s, 0}^{2} .Math. C_{M_{T}} & [135] \\ \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. (M_{0}^{2} + 2 .Math. M_{0} .Math. .Math. .Math. M) .Math. {Sb}_{w} .Math. C_{Y, 0} + \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. {Sb}_{w} (C_{Y,} .Math. .Math. .Math. + C_{Y,} .Math. + C_{Y_{A}, M} .Math. .Math. .Math. M + C_{Y, \dot{}} .Math. \frac{\dot{}}{M} + C_{Y, p} .Math. \frac{p}{M} + C_{Y, r} .Math. \frac{r}{M}) = - I_{xz} (m) .Math. \dot{p} + I_{zz} (m) .Math. \dot{r} - W_{MTOW} .Math. (C_{T, 0} + C_{T, M} .Math. .Math. .Math. M) .Math. E_{3} + m_{MTOW} .Math. b_{w}^{2} .Math. N_{s, 0}^{2} .Math. C_{M_{T}} .Math. .Math. & [136] \end{matrix}$

[0108] Bearing in mind expressions [118] and [119], the lateral-directional moments balance yields:

[00034] $\begin{matrix} .Math. .Math. p_{0} .Math. .Math. .Math. M_{0} .Math. .Math. .Math. {MSb}_{w} .Math. C_{R, 0} .Math. + \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. {Sb}_{w} (C_{R,} .Math. .Math. .Math. + C_{R,} .Math. + C_{R_{A}, M} .Math. .Math. .Math. M + C_{R, \dot{}} .Math. \frac{\dot{}}{M} + C_{R, p} .Math. \frac{p}{M} + C_{R, r} .Math. \frac{r}{M}) = I_{xx} (m) .Math. \dot{p} - I_{xz} (m) .Math. \dot{r} - W_{MTOW} .Math. .Math. .Math. C_{T, M} .Math. .Math. .Math. {ME}_{1} & [137] \\ .Math. .Math. p_{0} .Math. .Math. .Math. M_{0} .Math. .Math. .Math. {MSb}_{w} .Math. C_{Y, 0} + \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. {Sb}_{w} (C_{Y,} .Math. .Math. .Math. + C_{Y,} .Math. + C_{Y_{A}, M} .Math. .Math. .Math. M + C_{Y, \dot{}} .Math. \frac{\dot{}}{M} + C_{Y, p} .Math. \frac{p}{M} + C_{Y, r} .Math. \frac{r}{M}) = - I_{xz} (m) .Math. \dot{p} + I_{zz} (m) .Math. \dot{r} - W_{MTOW} .Math. .Math. .Math. C_{T, M} .Math. .Math. .Math. {ME}_{3} & [138] \\ .Math. .Math. p_{0} .Math. .Math. .Math. M_{0} .Math. .Math. {Sb}_{w} .Math. C_{R, 0} .Math. .Math. .Math. M + \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} (M^{2} .Math. C_{R,} .Math. .Math. .Math. + M^{2} .Math. C_{R,} .Math. + M^{2} .Math. C_{R_{A}, M} .Math. .Math. .Math. M + {MC}_{R, \dot{}} .Math. \dot{} + {MC}_{R, p} .Math. p + {MC}_{R, r} .Math. r) = I_{xx} (m) .Math. \dot{p} - I_{xz} (m) .Math. \dot{r} - W_{MTOW} .Math. .Math. .Math. C_{T, M} .Math. .Math. .Math. {ME}_{1} & [139] \\ .Math. .Math. p_{0} .Math. .Math. .Math. M_{0} .Math. {Sb}_{w} .Math. C_{Y, 0} .Math. .Math. .Math. M + \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} (M^{2} .Math. C_{Y,} .Math. .Math. .Math. + M^{2} .Math. C_{Y,} .Math. + M^{2} .Math. C_{Y_{A}, M} .Math. .Math. .Math. M + {MC}_{Y, \dot{}} .Math. \dot{} + {MC}_{Y, p} .Math. p + {MC}_{Y, r} .Math. r) = - I_{xz} (m) .Math. \dot{p} + I_{zz} (m) .Math. \dot{r} - W_{MTOW} .Math. .Math. .Math. C_{T, M} .Math. .Math. .Math. {ME}_{3} & [140] \\ .Math. .Math. p_{0} .Math. .Math. .Math. M_{0} .Math. {Sb}_{w} .Math. C_{R, 0} .Math. .Math. .Math. M + \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} (M^{2} .Math. C_{R,} .Math. .Math. .Math. + M^{2} .Math. C_{R,} .Math. + M^{2} .Math. C_{R_{A}, M} .Math. .Math. .Math. M + {MC}_{R, \dot{}} .Math. \dot{} + {MC}_{R, p} .Math. p + {MC}_{R, r} .Math. r) = I_{xx} (m) .Math. \dot{p} - I_{xz} (m) .Math. \dot{r} - W_{MTOW} .Math. .Math. .Math. C_{T, M} .Math. .Math. .Math. {ME}_{1} & [141] \\ .Math. .Math. p_{0} .Math. .Math. .Math. M_{0} .Math. {Sb}_{w} .Math. C_{Y, 0} .Math. .Math. .Math. M + \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} (M^{2} .Math. C_{Y,} .Math. .Math. .Math. + M^{2} .Math. C_{Y,} .Math. + M^{2} .Math. C_{Y_{A}, M} .Math. .Math. .Math. M + {MC}_{Y, \dot{}} .Math. \dot{} + {MC}_{Y, p} .Math. p + {MC}_{Y, r} .Math. r) = - I_{xz} (m) .Math. \dot{p} + I_{zz} (m) .Math. \dot{r} - W_{MTOW} .Math. .Math. .Math. C_{T, M} .Math. .Math. .Math. {ME}_{3} & [142] \end{matrix}$

[0109] Taking into account approximation [90]:

[00035] $\begin{matrix} \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} (2 .Math. M_{0} .Math. C_{R, 0} + \frac{W_{MTOW} .Math. C_{T, M} .Math. E_{1}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. {Sc}_{w}}) .Math. .Math. .Math. M + \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} (M^{2} .Math. C_{R,} .Math. .Math. .Math. + M^{2} .Math. C_{R,} .Math. + M_{0}^{2} (1 + 2 .Math. \frac{.Math. .Math. M}{M_{0}}) .Math. C_{R_{A}, M} .Math. .Math. .Math. M + {MC}_{R, \dot{}} .Math. \dot{} + {MC}_{R, p} .Math. p + {MC}_{R, r} .Math. r) = I_{xx} (m) .Math. \dot{p} - I_{xz} (m) .Math. \dot{r} & [143] \\ \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} (2 .Math. M_{0} .Math. C_{Y, 0} + \frac{W_{MTOW} .Math. C_{T, M} .Math. E_{3}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. {Sc}_{w}}) .Math. .Math. .Math. M + \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} (M^{2} .Math. C_{Y,} .Math. .Math. .Math. + M^{2} .Math. C_{Y,} .Math. + M_{0}^{2} (1 + 2 .Math. \frac{.Math. .Math. M}{M_{0}}) .Math. C_{Y_{A}, M} .Math. .Math. .Math. M + {MC}_{Y, \dot{}} .Math. \dot{} + {MC}_{Y, p} .Math. p + {MC}_{Y, r} .Math. r) = - I_{xz} (m) .Math. \dot{p} + I_{zz} (m) .Math. \dot{r} & [144] \\ .Math. C_{R, M} = 2 .Math. \frac{C_{R, 0}}{M_{0}} + C_{R_{A}, M} + \frac{W_{MTOW} .Math. C_{T, M} .Math. E_{1}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. {Sc}_{w} .Math. M_{0}^{2}} & [145] \\ .Math. C_{Y, M} = 2 .Math. \frac{C_{Y, 0}}{M_{0}} + C_{Y_{A}, M} + \frac{W_{MTOW} .Math. C_{T, M} .Math. E_{3}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. {Sc}_{w} .Math. M_{0}^{2}} & [146] \\ \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} (M_{0}^{2} .Math. C_{R, M} .Math. .Math. .Math. M + M^{2} .Math. C_{R,} .Math. .Math. .Math. + M^{2} .Math. C_{R,} .Math. + 2 .Math. M_{0} .Math. C_{R_{A}, M} .Math. .Math. .Math. M^{2} + {MC}_{R, \dot{}} .Math. \dot{} + {MC}_{R, p} .Math. p + {MC}_{R, r} .Math. r) = I_{xx} (m) .Math. \dot{p} - I_{xz} (m) .Math. \dot{r} & [147] \\ \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} (M_{0}^{2} .Math. C_{Y, M} .Math. .Math. .Math. M + M^{2} .Math. C_{Y,} .Math. .Math. .Math. + M^{2} .Math. C_{Y,} .Math. + 2 .Math. M_{0} .Math. C_{R_{A}, M} .Math. .Math. .Math. M^{2} + {MC}_{Y, \dot{}} .Math. \dot{} + {MC}_{Y, p} .Math. p + {MC}_{Y, r} .Math. r) = - I_{xz} (m) .Math. \dot{p} + I_{zz} (m) .Math. \dot{r} & [148] \\ \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} [M_{0}^{2} .Math. C_{R, M} .Math. .Math. .Math. M + M_{0}^{2} (1 + 2 .Math. \frac{.Math. .Math. M}{M_{0}}) .Math. C_{R,} .Math. + M_{0}^{2} (1 + 2 .Math. \frac{.Math. .Math. M}{M_{0}}) .Math. C_{R,} .Math. + 2 .Math. M_{0} .Math. C_{R_{A}, M} .Math. .Math. .Math. M^{2} + (M_{0} + .Math. .Math. M) .Math. C_{R, \dot{}} .Math. \dot{} + (M_{0} + .Math. .Math. M) .Math. C_{R, p} .Math. p + (M_{0} + .Math. .Math. M) .Math. C_{R, r} .Math. r] = I_{xx} (m) .Math. \dot{p} - I_{xz} (m) .Math. \dot{r} & [149] \\ \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} [M_{0}^{2} .Math. C_{Y, M} .Math. .Math. .Math. M + M_{0}^{2} (1 + 2 .Math. \frac{.Math. .Math. M}{M_{0}}) .Math. C_{Y,} .Math. + M_{0}^{2} (1 + 2 .Math. \frac{.Math. .Math. M}{M_{0}}) .Math. C_{Y,} .Math. + 2 .Math. M_{0} .Math. C_{R_{A}, M} .Math. .Math. .Math. M^{2} + (M_{0} + .Math. .Math. M) .Math. C_{Y, \dot{}} .Math. \dot{} + (M_{0} + .Math. .Math. M) .Math. C_{Y, p} .Math. p + (M_{0} + .Math. .Math. M) .Math. C_{Y, r} .Math. r] = - I_{xz} (m) .Math. \dot{p} + I_{zz} (m) .Math. \dot{r} & [150] \\ \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} .Math. M_{0}^{2} [C_{R,} .Math. .Math. .Math. + C_{R,} .Math. + C_{R, M} .Math. .Math. .Math. M + \frac{C_{R, \dot{}}}{M_{0}} .Math. \dot{} + \frac{C_{R, p}}{M_{0}} .Math. p + \frac{C_{R, r}}{M_{0}} .Math. r + (2 .Math. \frac{C_{R,}}{M_{0}} .Math. .Math. .Math. + 2 .Math. \frac{C_{R,} .Math.}{M_{0}} .Math. + 2 .Math. \frac{C_{R_{A}, M}}{M_{0}} .Math. .Math. .Math. M + \frac{C_{R, \dot{}}}{M_{0}^{2}} .Math. \dot{} + \frac{C_{R, p}}{M_{0}^{2}} .Math. p + \frac{C_{R, r}}{M_{0}^{2}} .Math. r) .Math. .Math. .Math. M] = I_{xx} (m) .Math. \dot{p} - I_{xz} (m) .Math. \dot{r} & [151] \\ \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} .Math. M_{0}^{2} [C_{Y,} .Math. .Math. .Math. + C_{Y,} .Math. + C_{Y, M} .Math. .Math. .Math. M + \frac{C_{Y, \dot{}}}{M_{0}} .Math. \dot{} + \frac{C_{Y, p}}{M_{0}} .Math. p + \frac{C_{Y, r}}{M_{0}} .Math. r + (2 .Math. \frac{C_{Y,}}{M_{0}} .Math. .Math. .Math. + 2 .Math. \frac{C_{Y,} .Math.}{M_{0}} .Math. + 2 .Math. \frac{C_{Y_{A}, M}}{M_{0}} .Math. .Math. .Math. M + \frac{C_{Y, \dot{}}}{M_{0}^{2}} .Math. \dot{} + \frac{C_{Y, p}}{M_{0}^{2}} .Math. p + \frac{C_{Y, r}}{M_{0}^{2}} .Math. r) .Math. .Math. .Math. M] = I_{xz} (m) .Math. \dot{p} + I_{zz} (m) .Math. \dot{r} & [152] \end{matrix}$

[0110] Expressions [147] and [148] constitute a linear observer for the unknowns {C.sub.R,M, C.sub.R,, C.sub.R,, C.sub.R.sub.A.sub.,M, C.sub.R,{dot over ()}, C.sub.R,p, C.sub.R,r, C.sub.Y,M, C.sub.Y,, C.sub.Y,, C.sub.Y.sub.A.sub.,M, C.sub.Y,{dot over ()}, C.sub.Y,p, C.sub.Y,r} and {I.sub.xx(m),I.sub.xz(m),I.sub.zz(m).sub.} valid for any Mach and actual mass. However, roll and yaw accelerations {dot over (p)} and {dot over (r)} and the derivative of the AOS, {dot over ()} are not typically native observables, which means numerical derivation should be required. To avoid the deriving noisy signals, a better idea is to integrate the known signals in [147] and [148] over time. Again, one option consists of neglecting the influence of {dot over ()}, which yields:

p.sub.0Sb.sub.w(M.sub.0.sup.2C.sub.R,MM+M.sup.2C.sub.R,+M.sup.2C.sub.R,+2M.sub.0C.sub.R.sub.A,MM.sup.2+MC.sub.R,pp+MC.sub.R,rr)=I.sub.xx(m){dot over (p)}I.sub.xz(m){dot over (r)}[153]

p.sub.0Sb.sub.w(M.sub.0.sup.2C.sub.Y,MM+M.sup.2C.sub.Y,+M.sup.2C.sub.Y,+2M.sub.0C.sub.R.sub.A,MM.sup.2+MC.sub.Y,pp+MC.sub.Y,rr)=I.sub.xx(m){dot over (p)}I.sub.xz(m){dot over (r)}[154]

f.sup.1=Mdt

f.sup.2=M.sup.2dt

f.sup.3=M.sup.2dt

f.sup.4=M.sup.2dt

f.sup.5=Mpdt

f.sup.6=Mrdt[155]

p.sub.0Sb.sub.w(M.sub.0.sup.2C.sub.R,Mf.sup.1+C.sub.R,f.sup.2+C.sub.R,f.sup.3+2M.sub.0C.sub.R.sub.A.sub.,Mf.sup.4+C.sub.R,pf.sup.5+C.sub.R,rf.sup.6)=I.sub.xx(m){dot over (p)}I.sub.xz(m){dot over (r)}+A[156]

p.sub.0Sb.sub.w(M.sub.0.sup.2C.sub.Y,Mf.sup.1+C.sub.Y,f.sup.2+C.sub.Y,f.sup.3+2M.sub.0C.sub.R.sub.A.sub.,Mf.sup.4+C.sub.Y,pf.sup.5+C.sub.Y,rf.sup.6)=I.sub.xx(m){dot over (p)}I.sub.xz(m){dot over (r)}+A[157]

[0111] Expressions [156] and [157] constitute a linear observer for the unknowns {C.sub.R,M, C.sub.R,, C.sub.R,, C.sub.R.sub.A.sub.,M, C.sub.R,p, C.sub.R,r, A, C.sub.Y,M, C.sub.Y,, C.sub.Y,, C.sub.Y.sub.A.sub.,M, C.sub.Y,p, C.sub.Y,r, C} and {I.sub.xx(m),I.sub.xz(m),I.sub.zz(m)} valid for any Mach and actual mass, where A and C are integration constants.

[0112] A second possibility consists on computing M{dot over ()} in expressions [147] and [148] as:

M{dot over ()}=M{dot over (M)} [Expression for {dot over (M)} TBD][158]

[0113] The third possibility consists of eliminating the dependency with {dot over ()} in [147] and [148][93] through, substituting expression [315] obtained from the linearization of the equations of linear motion, i.e.:

[00036] $\begin{matrix} .Math. F_{1} = C_{R, M} - C_{R, \dot{}} .Math. \frac{C_{y, M}}{C_{Q, \dot{}}} .Math. .Math. .Math. F_{2} = C_{R,} - C_{R, \dot{}} .Math. \frac{C_{Q,}}{C_{Q, \dot{}}} .Math. .Math. .Math. F_{3} = C_{R,} - C_{R, \dot{}} .Math. \frac{C_{Q,}}{C_{Q, \dot{}}} .Math. .Math. .Math. F_{4} = 2 .Math. (C_{R_{A}, M} - C_{R, \dot{}} .Math. \frac{C_{Q, M}}{C_{Q, \dot{}}}) .Math. .Math. .Math. F_{5} = C_{R, p} - C_{R, \dot{}} .Math. \frac{C_{Q, p}}{C_{Q, \dot{}}} .Math. .Math. .Math. F_{6} = C_{R, r} - C_{R, \dot{}} .Math. \frac{C_{Q, r}}{C_{Q, \dot{}}} .Math. .Math. .Math. F_{7} = - C_{R, \dot{}} .Math. \frac{C_{D,}}{C_{Q, \dot{}}} .Math. .Math. .Math. F_{8} = - C_{R, \dot{}} .Math. \frac{C_{D,}}{C_{Q, \dot{}}} .Math. .Math. .Math. F_{9} = - C_{R, \dot{}} .Math. \frac{C_{D, M}}{C_{Q, \dot{}}} .Math. .Math. .Math. F_{10} = - C_{R, \dot{}} .Math. \frac{C_{D, p}}{C_{Q, \dot{}}} .Math. .Math. .Math. F_{11} = - C_{R, \dot{}} .Math. \frac{C_{D, r}}{C_{Q, \dot{}}} .Math. .Math. .Math. F_{12} = C_{R, \dot{}} .Math. \frac{1}{C_{Q, \dot{}}} .Math. b_{w} & [159] \\ .Math. G_{1} = C_{Y, M} - C_{Y, \dot{}} .Math. \frac{C_{y, M}}{C_{Q, \dot{}}} .Math. .Math. .Math. G_{2} = C_{Y,} - C_{Y, \dot{}} .Math. \frac{C_{Q,}}{C_{Q, \dot{}}} .Math. .Math. .Math. G_{3} = C_{Y,} - C_{Y, \dot{}} .Math. \frac{C_{Q,}}{C_{Q, \dot{}}} .Math. .Math. .Math. G_{4} = 2 .Math. (C_{Y_{A}, M} - C_{Y, \dot{}} .Math. \frac{C_{Q, M}}{C_{Q, \dot{}}}) .Math. .Math. .Math. G_{5} = C_{Y, p} - C_{Y, \dot{}} .Math. \frac{C_{Q, p}}{C_{Q, \dot{}}} .Math. .Math. .Math. G_{6} = C_{Y, r} - C_{Y, \dot{}} .Math. \frac{C_{Q, r}}{C_{Q, \dot{}}} .Math. .Math. .Math. G_{7} = - C_{Y, \dot{}} .Math. \frac{C_{D,}}{C_{Q, \dot{}}} .Math. .Math. .Math. G_{8} = - C_{Y, \dot{}} .Math. \frac{C_{D,}}{C_{Q, \dot{}}} .Math. .Math. .Math. G_{9} = - C_{Y, \dot{}} .Math. \frac{C_{D, M}}{C_{Q, \dot{}}} .Math. .Math. .Math. G_{10} = - C_{Y, \dot{}} .Math. \frac{C_{D, p}}{C_{Q, \dot{}}} .Math. .Math. .Math. G_{11} = - C_{Y, \dot{}} .Math. \frac{C_{D, r}}{C_{Q, \dot{}}} .Math. .Math. .Math. G_{12} = C_{Y, \dot{}} .Math. \frac{1}{C_{Q, \dot{}}} .Math. b_{w} & [160] \\ \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} (M_{0}^{2} .Math. F_{1} .Math. .Math. .Math. M + F_{2} .Math. M^{2} .Math. .Math. .Math. + F_{3} .Math. M^{2} .Math. + M_{0} .Math. F_{4} .Math. .Math. .Math. M^{2} + F_{5} .Math. Mp + F_{6} .Math. Mr + F_{7} .Math. M^{2} .Math. .Math. .Math. .Math. .Math. + F_{8} .Math. M^{2} .Math.^{2} + F_{9} .Math. M^{2} .Math. .Math. .Math. .Math. .Math. M + F_{10} .Math. Mp .Math. .Math. + F_{11} .Math. Mr .Math. .Math.) = I_{xx} (m) .Math. \dot{p} - I_{xz} (m) .Math. \dot{r} + m_{0} .Math. F_{12} .Math. .Math. .Math. a_{2}^{BFS} & [161] \\ \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} (M_{0}^{2} .Math. G_{1} .Math. .Math. .Math. M + G_{2} .Math. M^{2} .Math. .Math. .Math. + G_{3} .Math. M^{2} .Math. + M_{0} .Math. G_{4} .Math. .Math. .Math. M^{2} + G_{5} .Math. Mp + G_{6} .Math. Mr + G_{7} .Math. M^{2} .Math. .Math. .Math. .Math. .Math. + G_{8} .Math. M^{2} .Math.^{2} + G_{9} .Math. M^{2} .Math. .Math. .Math. .Math. .Math. M + G_{10} .Math. Mp .Math. .Math. + G_{11} .Math. Mr .Math. .Math.) = - I_{xz} (m) .Math. \dot{p} + I_{zz} (m) .Math. \dot{r} + m_{0} .Math. G_{12} .Math. .Math. .Math. a_{2}^{BFS} & [162] \\ .Math. f^{7} = M^{2} .Math. .Math. .Math. .Math. .Math. dt .Math. .Math. .Math. f^{8} = M^{2} .Math.^{2} .Math. dt .Math. .Math. .Math. f^{9} = M^{2} .Math. .Math. .Math. Mdt .Math. .Math. .Math. f^{10} = Mp .Math. .Math. .Math. .Math. dt .Math. .Math. .Math. f^{11} = Mr .Math. .Math. .Math. .Math. dt .Math. .Math. .Math. f^{12} = .Math. .Math. a_{2}^{BFS} .Math. dt & [163] \end{matrix}$

Which allows rewriting expression [161] and [162] as:

[00037] $\begin{matrix} \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} (M_{0}^{2} .Math. F_{1} .Math. f^{1} + F_{2} .Math. f^{2} + F_{3} .Math. f^{3} + .Math. M_{0} .Math. F_{4} .Math. f^{4} + F_{5} .Math. f^{5} + F_{6} .Math. f^{6} + F_{7} .Math. f^{7} + F_{8} .Math. f^{8} + F_{9} .Math. f^{9} + F_{10} .Math. f^{10} + F_{11} .Math. f^{11}) = I_{xx} (m) .Math. p - I_{xz} (m) .Math. r + m_{0} .Math. F_{12} .Math. f^{12} + A & [164] \\ \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} (M_{0}^{2} .Math. G_{1} .Math. f^{1} + G_{2} .Math. f^{2} + G_{3} .Math. f^{3} + .Math. M_{0} .Math. G_{4} .Math. f^{4} + G_{5} .Math. f^{5} + G_{6} .Math. f^{6} + G_{7} .Math. f^{7} + G_{8} .Math. f^{8} + G_{9} .Math. f^{9} + G_{10} .Math. f^{10} + G_{11} .Math. f^{11}) = - I_{xz} (m) .Math. p + I_{zz} (m) .Math. r + m_{0} .Math. G_{12} .Math. f^{12} + C & [165] \end{matrix}$

where {A,C} are integration constants.

[0114] Expressions [164] and [165] constitute a linear observer for the unknowns {F.sub.1, . . . , F.sub.12, A, G.sub.1, . . . , G.sub.12, C} and {I.sub.xx(m), I.sub.xz (m), I.sub.zz(m)} valid for any Mach and actual mass. Thus by performing both experiments (i.e. the ailerons and rudder doublets) for different cases of Mach at the trim condition and actual mass, an overdetermined linear observation system is obtained that allows estimating the mentioned unknowns. In effect, if it is considered a case for the identification of {I.sub.xx(m),I.sub.xz(m),I.sub.zz(m)} for the combination {m.sub.i,M.sub.j}, i.e., at the trim condition at Mach M.sub.j (j=1, . . . , l.sub.i) with actual mass m.sub.i (i=1, . . . , k) through aileron perturbation (expression [164]). For each case of M.sub.j, n.sub.j samples of the response to a lateral perturbation are recorded. Thus:

[00038] $\begin{matrix} H_{ij} = [\begin{matrix} - \frac{p_{j .Math. .Math. 1}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & \frac{r_{j .Math. .Math. 1}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & 0 & M_{j}^{2} .Math. f_{j .Math. .Math. 1}^{1} & f_{j .Math. .Math. 1}^{2} & f_{j .Math. .Math. 1}^{3} & M_{j} .Math. f_{j .Math. .Math. 1}^{4} & f_{j .Math. .Math. 1}^{5} & f_{j .Math. .Math. 1}^{6} & f_{j .Math. .Math. 1}^{7} & f_{j .Math. .Math. 1}^{8} & f_{j .Math. .Math. 1}^{9} & f_{j .Math. .Math. 1}^{10} & f_{j .Math. .Math. 1}^{11} & - \frac{m_{i} .Math. f_{j .Math. .Math. 1}^{12}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & - \frac{1}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & 0 \\ - \frac{p_{j .Math. .Math. 2}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & \frac{r_{j .Math. .Math. 2}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & 0 & M_{j}^{2} .Math. f_{j .Math. .Math. 2}^{1} & f_{j .Math. .Math. 2}^{2} & f_{j .Math. .Math. 2}^{3} & M_{j} .Math. f_{j .Math. .Math. 2}^{4} & f_{j .Math. .Math. 2}^{5} & f_{j .Math. .Math. 2}^{6} & f_{j .Math. .Math. 2}^{7} & f_{j .Math. .Math. 2}^{8} & f_{j .Math. .Math. 2}^{9} & f_{j .Math. .Math. 2}^{10} & f_{j .Math. .Math. 2}^{11} & - \frac{m_{i} .Math. f_{j .Math. .Math. 2}^{12}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & - \frac{1}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & 0 \\ .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. \\ - \frac{p_{j .Math. .Math. n_{j}}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & \frac{r_{j .Math. .Math. n_{j}}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & 0 & M_{j}^{2} .Math. f_{j .Math. .Math. n_{j}}^{1} & f_{j .Math. .Math. n_{j}}^{2} & f_{j .Math. .Math. n_{j}}^{3} & M_{j} .Math. f_{j .Math. .Math. n_{j}}^{4} & f_{j .Math. .Math. n_{j}}^{5} & f_{j .Math. .Math. n_{j}}^{6} & f_{j .Math. .Math. n_{j}}^{7} & f_{j .Math. .Math. n_{j}}^{8} & f_{j .Math. .Math. n_{j}}^{9} & f_{j .Math. .Math. n_{j}}^{10} & f_{j .Math. .Math. n_{j}}^{11} & - \frac{m_{i} .Math. f_{j .Math. .Math. n_{j}}^{12}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & - \frac{1}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & 0 \end{matrix}] & [166] \\ z_{ij} = [\begin{matrix} I_{xx} (m_{i}) \\ I_{xz} (m_{i}) \\ I_{zz} (m_{i}) \\ F_{1} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ F_{2} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ F_{3} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ F_{4} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ F_{5} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ F_{7} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ F_{8} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ F_{9} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ F_{10} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ F_{11} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ F_{12} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ A_{ij} \\ B_{ij} \end{matrix}] & [167] \\ O_{ij} = [\begin{matrix} 0 \\ 0 \\ .Math. \\ 0 \end{matrix}] & [168] \\ F_{ij} = [\begin{matrix} M_{j}^{2} .Math. f_{j .Math. .Math. 1}^{1} & f_{j .Math. .Math. 1}^{2} & f_{j .Math. .Math. 1}^{3} & M_{j} .Math. f_{j .Math. .Math. 1}^{4} & f_{j .Math. .Math. 1}^{5} & f_{j .Math. .Math. 1}^{6} & f_{j .Math. .Math. 1}^{7} & f_{j .Math. .Math. 1}^{8} & f_{j .Math. .Math. 1}^{9} & f_{j .Math. .Math. 1}^{10} & f_{j .Math. .Math. 1}^{11} & - \frac{m_{i} .Math. f_{j .Math. .Math. 1}^{12}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & - \frac{1}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & 0 \\ M_{j}^{2} .Math. f_{j .Math. .Math. 2}^{1} & f_{j .Math. .Math. 2}^{2} & f_{j .Math. .Math. 2}^{3} & M_{j} .Math. f_{j .Math. .Math. 2}^{4} & f_{j .Math. .Math. 2}^{5} & f_{j .Math. .Math. 2}^{6} & f_{j .Math. .Math. 2}^{7} & f_{j .Math. .Math. 2}^{8} & f_{j .Math. .Math. 2}^{9} & f_{j .Math. .Math. 2}^{10} & f_{j .Math. .Math. 2}^{11} & - \frac{m_{i} .Math. f_{j .Math. .Math. 2}^{12}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & - \frac{1}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & 0 \\ .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. \\ M_{j}^{2} .Math. f_{j .Math. .Math. n_{j}}^{1} & f_{j .Math. .Math. n_{j}}^{2} & f_{j .Math. .Math. n_{j}}^{3} & M_{j} .Math. f_{j .Math. .Math. n_{j}}^{4} & f_{j .Math. .Math. n_{j}}^{5} & f_{j .Math. .Math. n_{j}}^{6} & f_{j .Math. .Math. n_{j}}^{7} & f_{j .Math. .Math. n_{j}}^{8} & f_{j .Math. .Math. n_{j}}^{9} & f_{j .Math. .Math. n_{j}}^{10} & f_{j .Math. .Math. n_{j}}^{11} & - \frac{m_{i} .Math. f_{j .Math. .Math. n_{j}}^{12}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & - \frac{1}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & 0 \end{matrix}] & [169] \\ G_{ij} = \frac{1}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} [\begin{matrix} - p_{j .Math. .Math. 1} & r_{j .Math. .Math. 1} & 0 \\ - p_{j .Math. .Math. 2} & r_{j .Math. .Math. 2} & 0 \\ .Math. & .Math. & .Math. \\ - p_{{jn}_{j}} & r_{{jn}_{j}} & 0 \end{matrix}] & [170] \\ C_{R, ij} = [\begin{matrix} F_{1} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ F_{2} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ F_{3} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ F_{4} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ F_{5} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ F_{7} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ F_{8} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ F_{9} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ F_{10} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ F_{11} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ F_{12} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ A_{ij} \\ B_{ij} \end{matrix}] & [171] \end{matrix}$

[0115] It is considered the same case for the identification of {I.sub.xx(m),I.sub.xz(m),I.sub.zz(m)} for the combination {m.sub.i,M.sub.j}, i.e., at the trim condition at Mach M.sub.j (j=1, . . . , l.sub.i) with actual mass m.sub.i (i=1, . . . , k) through ruder perturbation (expression [165]). For each case of samples of the response to a lateral perturbation are recorded. Thus:

[00039] $\begin{matrix} H_{ij} = [\begin{matrix} 0 & \frac{p_{j .Math. .Math. 1}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & - \frac{r_{j .Math. .Math. 1}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & M_{j}^{2} .Math. f_{j .Math. .Math. 1}^{1} & f_{j .Math. .Math. 1}^{2} & f_{j .Math. .Math. 1}^{3} & M_{j} .Math. f_{j .Math. .Math. 1}^{4} & f_{j .Math. .Math. 1}^{5} & f_{j .Math. .Math. 1}^{6} & f_{j .Math. .Math. 1}^{7} & f_{j .Math. .Math. 1}^{8} & f_{j .Math. .Math. 1}^{9} & f_{j .Math. .Math. 1}^{10} & f_{j .Math. .Math. 1}^{11} & - \frac{m_{i} .Math. f_{j .Math. .Math. 1}^{12}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & 0 & - \frac{1}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} \\ 0 & \frac{p_{j .Math. .Math. 2}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & - \frac{r_{j .Math. .Math. 2}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & M_{j}^{2} .Math. f_{j .Math. .Math. 2}^{1} & f_{j .Math. .Math. 2}^{2} & f_{j .Math. .Math. 2}^{3} & M_{j} .Math. f_{j .Math. .Math. 2}^{4} & f_{j .Math. .Math. 2}^{5} & f_{j .Math. .Math. 2}^{6} & f_{j .Math. .Math. 2}^{7} & f_{j .Math. .Math. 2}^{8} & f_{j .Math. .Math. 2}^{9} & f_{j .Math. .Math. 2}^{10} & f_{j .Math. .Math. 2}^{11} & - \frac{m_{i} .Math. f_{j .Math. .Math. 2}^{12}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & 0 & - \frac{1}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} \\ .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. \\ 0 & \frac{p_{j .Math. .Math. h_{j}}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & - \frac{r_{j .Math. .Math. h_{j}}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & M_{j}^{2} .Math. f_{j .Math. .Math. h_{j}}^{1} & f_{j .Math. .Math. h_{j}}^{2} & f_{j .Math. .Math. h_{j}}^{3} & M_{j} .Math. f_{j .Math. .Math. h_{j}}^{4} & f_{j .Math. .Math. h_{j}}^{5} & f_{j .Math. .Math. h_{j}}^{6} & f_{j .Math. .Math. h_{j}}^{7} & f_{j .Math. .Math. h_{j}}^{8} & f_{j .Math. .Math. h_{j}}^{9} & f_{j .Math. .Math. h_{j}}^{10} & f_{j .Math. .Math. h_{j}}^{11} & - \frac{m_{i} .Math. f_{j .Math. .Math. h_{j}}^{12}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & 0 & - \frac{1}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} \end{matrix}] & [172] \\ z_{ij} = [\begin{matrix} I_{xx} (m_{i}) \\ I_{xz} (m_{i}) \\ I_{zz} (m_{i}) \\ G_{1} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ G_{2} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ G_{3} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ G_{4} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ G_{5} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ G_{6} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ G_{7} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ G_{8} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ G_{9} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ G_{10} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ G_{11} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ G_{12} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ A_{ij} \\ B_{ij} \end{matrix}] & [173] \\ O_{ij} = [\begin{matrix} 0 \\ 0 \\ .Math. \\ 0 \end{matrix}] & [174] \\ J_{ij} = [\begin{matrix} M_{j}^{2} .Math. f_{j .Math. .Math. 1}^{1} & f_{j .Math. .Math. 1}^{2} & f_{j .Math. .Math. 1}^{3} & M_{j} .Math. f_{j .Math. .Math. 1}^{4} & f_{j .Math. .Math. 1}^{5} & f_{j .Math. .Math. 1}^{6} & f_{j .Math. .Math. 1}^{7} & f_{j .Math. .Math. 1}^{8} & f_{j .Math. .Math. 1}^{9} & f_{j .Math. .Math. 1}^{10} & f_{j .Math. .Math. 1}^{11} & - \frac{m_{i} .Math. f_{j .Math. .Math. 1}^{12}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & 0 & - \frac{1}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} \\ M_{j}^{2} .Math. f_{j .Math. .Math. 2}^{1} & f_{j .Math. .Math. 2}^{2} & f_{j .Math. .Math. 2}^{3} & M_{j} .Math. f_{j .Math. .Math. 2}^{4} & f_{j .Math. .Math. 2}^{5} & f_{j .Math. .Math. 2}^{6} & f_{j .Math. .Math. 2}^{7} & f_{j .Math. .Math. 2}^{8} & f_{j .Math. .Math. 2}^{9} & f_{j .Math. .Math. 2}^{10} & f_{j .Math. .Math. 2}^{11} & - \frac{m_{i} .Math. f_{j .Math. .Math. 2}^{12}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & 0 & - \frac{1}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} \\ .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. \\ M_{j}^{2} .Math. f_{j .Math. .Math. h_{j}}^{1} & f_{j .Math. .Math. h_{j}}^{2} & f_{j .Math. .Math. h_{j}}^{3} & M_{j} .Math. f_{j .Math. .Math. h_{j}}^{4} & f_{j .Math. .Math. h_{j}}^{5} & f_{j .Math. .Math. h_{j}}^{6} & f_{j .Math. .Math. h_{j}}^{7} & f_{j .Math. .Math. h_{j}}^{8} & f_{j .Math. .Math. h_{j}}^{9} & f_{j .Math. .Math. h_{j}}^{10} & f_{j .Math. .Math. h_{j}}^{11} & - \frac{m_{i} .Math. f_{j .Math. .Math. h_{j}}^{12}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} & 0 & - \frac{1}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} \end{matrix}] & [175] \\ K_{ij} = \frac{1}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w}} [\begin{matrix} 0 & p_{j .Math. .Math. 1} & - r_{j .Math. .Math. 1} \\ 0 & p_{j .Math. .Math. 2} & - r_{j .Math. .Math. 2} \\ .Math. & .Math. & .Math. \\ 0 & p_{j .Math. .Math. h_{j}} & - r_{j .Math. .Math. h_{j}} \end{matrix}] & [176] \\ C_{Y, ij} = [\begin{matrix} G_{1} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ G_{2} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ G_{3} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ G_{4} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ G_{5} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ G_{6} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ G_{7} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ G_{8} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ G_{9} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ G_{10} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ G_{11} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ G_{12} (_{ij}, M_{ij}, {.Math.}_{h, ij}) \\ A_{ij} \\ B_{ij} \end{matrix}] & [177] \\ H_{i} = [\begin{matrix} G_{i .Math. .Math. 1} & F_{i .Math. .Math. 1} & 0 & .Math. & 0 & 0 & 0 & .Math. & 0 \\ G_{i .Math. .Math. 2} & 0 & F_{i .Math. .Math. 2} & .Math. & 0 & 0 & 0 & .Math. & 0 \\ .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. \\ G_{{il}_{i}} & 0 & 0 & .Math. & F_{{il}_{i}} & 0 & 0 & .Math. & 0 \\ K_{i .Math. .Math. 1} & 0 & 0 & .Math. & 0 & J_{i .Math. .Math. 1} & 0 & .Math. & 0 \\ K_{i .Math. .Math. 2} & 0 & 0 & .Math. & 0 & 0 & J_{i .Math. .Math. 2} & .Math. & 0 \\ .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. \\ K_{{il}_{i}} & 0 & 0 & .Math. & 0 & 0 & 0 & .Math. & J_{{il}_{i}} \end{matrix}] .Math. .Math. \dim (H_{i}) = ({.Math.}_{j = 1}^{l_{i}} .Math. n_{j} + {.Math.}_{j = 1}^{l_{i}} .Math. h_{j}) (3 + 14 .Math. l_{i}) & [178] \\ z_{i} = [\begin{matrix} I_{xx} (m_{i}) \\ I_{xz} (m_{i}) \\ I_{zz} (m_{i}) \\ C_{R, i .Math. .Math. 1} \\ C_{R, i .Math. .Math. 2} \\ .Math. \\ C_{R, i .Math. .Math. l_{i}} \\ C_{Y, i .Math. .Math. 1} \\ C_{Y, i .Math. .Math. 2} \\ .Math. \\ C_{Y, i .Math. .Math. l_{i}} \end{matrix}] .Math. .Math. \dim (z_{i}) = (3 + 14 .Math. l_{i}) 1 & [179] \\ O_{i} = [\begin{matrix} 0 \\ 0 \\ .Math. \\ 0 \end{matrix}] .Math. .Math. \dim (O_{i}) = ({.Math.}_{j = 1}^{l_{i}} .Math. n_{j} + {.Math.}_{j = 1}^{l_{i}} .Math. h_{j}) 1 & [180] \\ z_{i} = {(H_{i}^{T} .Math. H_{i})}^{- 1} .Math. H_{i}^{T} .Math. O_{i} & [181] \\ I_{xx} (m_{i}) = z_{i} [1] & [182] \\ I_{xz} (m_{i}) = z_{i} [2] & [183] \\ I_{zz} (m_{i}) = z_{i} [3] & [184] \end{matrix}$

[0116] Repeating the process for different values of the actual mass would render the variation of {I.sub.xx(m),I.sub.xz(m),I.sub.zz(m)} with m for i=1, . . . , k

[0117] Now, referring back to the balance condition governed by expressions [73], [118] and [119] for the identification cases considered:

[00040] $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M_{ij}^{2} .Math. {Sb}_{w} .Math. C_{R, 0} (_{ij}, M_{ij}, {.Math.}_{h, ij}) = - T_{ij} .Math. E_{1} - M_{T, ij} = - T_{ij} .Math. E_{1} - m_{MTOW} .Math. b_{w}^{2} .Math. N_{s, ij}^{2} .Math. C_{M_{T}}$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M_{ij}^{2} .Math. {Sc}_{w} .Math. C_{P, 0} (_{ij}, M_{ij}, {.Math.}_{h, ij}) = - T_{ij} .Math. E_{2} + M_{T, ij} .Math. = - T_{ij} .Math. E_{2} + m_{MTOW} .Math. b_{w}^{2} .Math. N_{s, ij}^{2} .Math. C_{M_{T}} .Math.$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M_{ij}^{2} .Math. {Sb}_{w} .Math. C_{Y, 0} (_{ij}, M_{ij}, {.Math.}_{h, ij}) = - T_{ij} .Math. E_{3} + M_{T, ij} .Math. .Math. = - T_{ij} .Math. E_{3} + m_{MTOW} .Math. b_{w}^{2} .Math. N_{s, ij}^{2} .Math. C_{M_{T}} .Math. .Math. .Math.$

with T.sub.ij for the different flight conditions given by the thrust model already known and M.sub.T given by expression [64].

[0118] Expressions [185], [186] and [187] contain 4 more unknowns than equations. To reduce the system encompassed by [185], [186] and [187] to an overdetermined system, it is introduced a Taylor expansion for C.sub.R,0, C.sub.p,0 and C.sub.Y,0 in the form (approximations of higher orders could be employed, if necessary):

C.sub.R,0(.sub.ij,M.sub.ij,.sub.h,ij)=r.sub.0+r.sub.1.sub.ij+r.sub.M1M.sub.ij+r.sub..sub.h.sub.1.sub.h,ij[188]

C.sub.P,0(.sub.ij,M.sub.ij,.sub.h,ij)=p.sub.0+p.sub.1.sub.ij+p.sub.M1M.sub.ij+p.sub..sub.h.sub.1.sub.h,ij[189]

C.sub.Y,0(.sub.ij,M.sub.ij,.sub.h,ij)=y.sub.0+y.sub.1.sub.ij+y.sub.M1M.sub.ij+y.sub..sub.h.sub.1.sub.h,ij[190]

which renders the following linear observation system for {r.sub.0, r.sub.1, r.sub.M1, r.sub..sub.h.sub.1, p.sub.0, p.sub.1, p.sub.M1, p.sub..sub.h.sub.1, y.sub.0, y.sub.1, y.sub.M1, y.sub..sub.h.sub.1)} and {E.sub.1, E.sub.2, E.sub.3, C.sub.M.sub.T}:

[00041] $H_{i} = [.Math. \begin{matrix} \frac{T_{i .Math. .Math. 1}}{\begin{matrix} \frac{1}{2} \\ .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} \end{matrix}} & 0 & 0 & \frac{\begin{matrix} m_{MTOW} \\ b_{w}^{2} .Math. N_{s, i .Math. .Math. 1}^{2} \end{matrix}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} & M_{i .Math. .Math. 1}^{2} & M_{i .Math. .Math. 1}^{2} .Math._{i .Math. .Math. 1} & M_{i .Math. .Math. 1}^{3} & M_{i .Math. .Math. 1}^{2} .Math. {.Math.}_{h, i .Math. .Math. 1} & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 \\ \frac{T_{i .Math. .Math. 2}}{\begin{matrix} \frac{1}{2} \\ .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} \end{matrix}} & 0 & 0 & \frac{\begin{matrix} m_{MTOW} \\ b_{w}^{2} .Math. N_{s, i .Math. .Math. 2}^{2} \end{matrix}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} & M_{i .Math. .Math. 2}^{2} & M_{i .Math. .Math. 2}^{2} .Math._{i .Math. .Math. 2} & M_{i .Math. .Math. 2}^{3} & M_{i .Math. .Math. 2}^{2} .Math. {.Math.}_{h, i .Math. .Math. 2} & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 \\ .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. \\ \frac{T_{{il}_{i}}}{\begin{matrix} \frac{1}{2} \\ .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} \end{matrix}} & 0 & 0 & \frac{\begin{matrix} m_{MTOW} \\ b_{w}^{2} .Math. N_{s, {il}_{i}}^{2} \end{matrix}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} & M_{{il}_{i}}^{2} & M_{{il}_{i}}^{2} .Math._{{il}_{i}} & M_{{il}_{i}}^{3} & M_{{il}_{i}}^{2} .Math. {.Math.}_{h, {il}_{i}} & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 \\ 0 & \frac{T_{i .Math. .Math. 1}}{\begin{matrix} \frac{1}{2} \\ .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} \end{matrix}} & 0 & - \frac{\begin{matrix} m_{MTOW} \\ b_{w}^{2} .Math. N_{s, i .Math. .Math. 1}^{2} .Math. \end{matrix}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} & 0 & 0 & 0 & 0 & M_{i .Math. .Math. 1}^{2} & M_{i .Math. .Math. 1}^{2} .Math._{i .Math. .Math. 1} & M_{i .Math. .Math. 1}^{3} & M_{i .Math. .Math. 1}^{2} .Math. {.Math.}_{h, i .Math. .Math. 1} & 0 & 0 & 0 & 0 \\ 0 & \frac{T_{i .Math. .Math. 2}}{\begin{matrix} \frac{1}{2} \\ .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} \end{matrix}} & 0 & - \frac{\begin{matrix} m_{MTOW} \\ b_{w}^{2} .Math. N_{s, i .Math. .Math. 2}^{2} .Math. \end{matrix}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} & 0 & 0 & 0 & 0 & M_{i .Math. .Math. 2}^{2} & M_{i .Math. .Math. 2}^{2} .Math._{i .Math. .Math. 2} & M_{i .Math. .Math. 2}^{3} & M_{i .Math. .Math. 2}^{2} .Math. {.Math.}_{h, i .Math. .Math. 2} & 0 & 0 & 0 & 0 \\ .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. \\ 0 & \frac{T_{{il}_{i}}}{\begin{matrix} \frac{1}{2} \\ .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} \end{matrix}} & 0 & - \frac{\begin{matrix} m_{MTOW} \\ b_{w}^{2} .Math. N_{s, {il}_{i}}^{2} .Math. \end{matrix}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} & 0 & 0 & 0 & 0 & M_{{il}_{i}}^{2} & M_{{il}_{i}}^{2} .Math._{{il}_{i}} & M_{{il}_{i}}^{3} & M_{{il}_{i}}^{2} .Math. {.Math.}_{h, {il}_{i}} & 0 & 0 & 0 & 0 \\ 0 & 0 & \frac{T_{i .Math. .Math. 1}}{\begin{matrix} \frac{1}{2} \\ .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} \end{matrix}} & - \frac{\begin{matrix} m_{MTOW} \\ b_{w}^{2} .Math. N_{s, i .Math. .Math. 1}^{2} .Math. .Math. \end{matrix}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 & M_{i .Math. .Math. 1}^{2} & M_{i .Math. .Math. 1}^{2} .Math._{i .Math. .Math. 1} & M_{i .Math. .Math. 1}^{3} & M_{i .Math. .Math. 1}^{2} .Math. {.Math.}_{h, i .Math. .Math. 1} \\ 0 & 0 & \frac{T_{i .Math. .Math. 1}}{\begin{matrix} \frac{1}{2} \\ .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} \end{matrix}} & - \frac{\begin{matrix} m_{MTOW} \\ b_{w}^{2} .Math. N_{s, i .Math. .Math. 2}^{2} .Math. .Math. \end{matrix}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 & M_{i .Math. .Math. 2}^{2} & M_{i .Math. .Math. 2}^{2} .Math._{i .Math. .Math. 2} & M_{i .Math. .Math. 2}^{3} & M_{i .Math. .Math. 2}^{2} .Math. {.Math.}_{h, i .Math. .Math. 2} \\ .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. & .Math. \\ 0 & 0 & \frac{T_{{il}_{i}}}{\begin{matrix} \frac{1}{2} \\ .Math. .Math. p_{0} .Math. .Math. .Math. {Sb}_{w} \end{matrix}} & - \frac{\begin{matrix} m_{MTOW} \\ b_{w}^{2} .Math. N_{s, {il}_{i}}^{2} .Math. .Math. \end{matrix}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{w}} & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 & M_{{il}_{i}}^{2} & M_{{il}_{i}}^{2} .Math._{{il}_{i}} & M_{{il}_{i}}^{3} & M_{{il}_{i}}^{2} .Math. {.Math.}_{h, {il}_{i}} \end{matrix} .Math.] .Math. .Math. z = [\begin{matrix} E_{1} \\ E_{2} \\ E_{3} \\ C_{M_{T}} \\ r_{0} \\ r_{1} \\ r_{M .Math. .Math. 1} \\ r_{{.Math.}_{h} .Math. 1} \\ p_{0} \\ p_{1} \\ p_{M .Math. .Math. 1} \\ p_{{.Math.}_{h} .Math. 1} \\ y_{0} \\ y_{1} \\ y_{M .Math. .Math. 1} \\ y_{{.Math.}_{h} .Math. 1} \end{matrix}] .Math. .Math. \dim (z) = 16 1 .Math. .Math. O_{i} = [\begin{matrix} 0 \\ 0 \\ .Math. \\ 0 \end{matrix}] .Math. .Math. \dim (O_{i}) = 3 .Math. .Math. l_{i} 1 .Math. .Math. H = [\begin{matrix} H_{1} \\ H_{2} \\ .Math. \\ H_{k} \end{matrix}] .Math. .Math. \dim (H) = 3 .Math. .Math. {kl}_{i} 16 .Math. .Math. O = [\begin{matrix} O_{1} \\ O_{2} \\ .Math. \\ O_{k} \end{matrix}] .Math. .Math. \dim (O) = 3 .Math. .Math. {kl}_{i} 1 .Math. .Math. z = {(H^{T} .Math. H)}^{- 1} .Math. H^{T} .Math. O .Math. .Math. E_{1} = z [1] .Math. .Math. E_{2} = z [2] .Math. .Math. E_{3} = z [3] .Math. .Math. C_{M_{T}} = z [4]$

Step 250 for Determination of the Aerodynamic Moments Model

[0119] Once the parameters of the propulsive models have been identified, expression [53] can be used again, this time as a direct observer for the aerodynamic moment coefficients as follows:

[00042] $\begin{matrix} [\begin{matrix} \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. {Sb}_{w} .Math. C_{R} \\ \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. {Sc}_{w} .Math. C_{P} \\ \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. {Sb}_{w} .Math. C_{Y} \end{matrix}] = [\begin{matrix} I_{xx} & 0 & - I_{xz} \\ 0 & I_{yy} & 0 \\ - I_{xz} & 0 & I_{zz} \end{matrix}] .Math. [\begin{matrix} \dot{p} - D_{9} .Math. pq + D_{10} .Math. qr \\ \dot{q} + D_{11} (p^{2} - r^{2}) + D_{12} .Math. pr \\ \dot{r} - D_{13} .Math. pq + D_{14} .Math. qr \end{matrix}] - [\begin{matrix} R_{P} \\ P_{P} \\ Y_{P} \end{matrix}] & [201] \\ .Math. [\begin{matrix} R_{P} \\ P_{P} \\ Y_{P} \end{matrix}] = T [\begin{matrix} E_{1} \\ E_{2} \\ E_{3} \end{matrix}] + M_{T} [\begin{matrix} 1 \\ - \\ - .Math. \end{matrix}] & [202] \end{matrix}$

[0120] Expression [201] allows the estimation of the aerodynamic moment coefficients C.sub.R, C.sub.P and C.sub.Y in terms of their dependency variables through flight testing, for which the respective domains have to be swept, which can be done manually and or with the help of control loops such as AOS-on-rudder, bank-on-ailerons and q-on-elevator.

[0121] In the general case:

C.sub.R=C.sub.R(,,M,{dot over ({circumflex over ()})},{dot over ({circumflex over ()})},{circumflex over (p)},{circumflex over (q)},{circumflex over (r)},.sub.h,.sub.a,.sub.r,.sub.e)[203]

C.sub.P=C.sub.P(,,M,{dot over ({circumflex over ()})},{dot over ({circumflex over ()})},{circumflex over (p)},{circumflex over (q)},{circumflex over (r)},.sub.h,.sub.a,.sub.r,.sub.e)[204]

C.sub.Y=C.sub.Y(,,M,{dot over ({circumflex over ()})},{dot over ({circumflex over ()})},{circumflex over (p)},{circumflex over (q)},{circumflex over (r)},.sub.h,.sub.a,.sub.r,.sub.e)[205]

[0122] Typical symmetry assumptions render:

C.sub.R=C.sub.R(,,M,{dot over ({circumflex over ()})},{dot over ({circumflex over ()})},{circumflex over (p)},{circumflex over (r)},.sub.a,.sub.r)[206]

C.sub.P=C.sub.P(,,M,{dot over ({circumflex over ()})},{dot over ({circumflex over ()})},{circumflex over (q)},.sub.h,.sub.e,.sub.r)[207]

C.sub.Y=C.sub.Y(,,M,{dot over ({circumflex over ()})},{dot over ({circumflex over ()})},{circumflex over (p)},{circumflex over (r)},.sub.a,.sub.r)[208]

[0123] Additionally, for uncompressible aerodynamics and small AOS the main dependencies of the aerodynamics coefficients under the quasi-steady state and typical symmetry assumptions remain:

C.sub.R=C.sub.R(,{dot over ({circumflex over ()})},{circumflex over (p)},{circumflex over (r)},.sub.a,.sub.r)[209]

C.sub.P=C.sub.P(,{dot over ({circumflex over ()})},{circumflex over (q)},.sub.h,.sub.e,.sub.r)[210]

C.sub.Y=C.sub.Y(,{dot over ({circumflex over ()})},{circumflex over (p)},{circumflex over (r)},.sub.a,.sub.r)[211]

[0124] Finally, for quasi-steady flight:

C.sub.R=C.sub.R(,{circumflex over (p)},{circumflex over (r)},.sub.a,.sub.r)[212]

C.sub.P=C.sub.P(,{circumflex over (q)},.sub.h,.sub.e,.sub.r)[213]

C.sub.Y=C.sub.Y(,{circumflex over (p)},{circumflex over (r)},.sub.a,.sub.r)[214]

[0125] For fixed-wing UAVs flying at moderate speeds, in quasi-steady coordinated flight, the aerodynamic model of interest is:

C.sub.L=C.sub.L(,{circumflex over (q)},.sub.h,.sub.e)[215]

C.sub.Q=C.sub.Q(,{circumflex over (p)},{circumflex over (r)},.sub.a,.sub.r)[216]

C.sub.D=C.sub.D(,{circumflex over (q)},.sub.h,.sub.e)[217]

C.sub.R=C.sub.R(,{circumflex over (p)},{circumflex over (r)},.sub.a,.sub.r)[218]

C.sub.P=C.sub.P(,{circumflex over (q)},.sub.h,.sub.e,.sub.r)[219]

C.sub.Y=C.sub.Y(,{circumflex over (p)},{circumflex over (r)},.sub.a,.sub.r)[220]

[0126] The propulsive model is characterized by:

C.sub.T(,,M,.sub.T)[221]

{,}[222]

{E.sub.1,E.sub.2,E.sub.3}[223]

C.sub.M.sub.T[224]

[0127] The mass and inertia properties are characterized by:

C.sub.F=C.sub.F(,,.sub.T)[225]

I.sub.xx=I.sub.xx(m)[226]

I.sub.yy=I.sub.yy(m)[227]

I.sub.zz=I.sub.zz(m)[228]

I.sub.xz=I.sub.xz(m)[229]

[0128] Based on the afore described approach, an automated flight test procedure may be implemented that obtains an APM with the help of several basic control loops: [0129] AOS-on-rudder [0130] Altitude-on-bank [0131] Speed-on-elevator [0132] AOA-on-elevator [0133] Pitch rate-on-elevator

[0134] Once the propulsive forces and moments model are well identified, the aerodynamic model may be continuously identified while-on-the-flight, populating the respective n-dimensional hypercubes allocated to each aerodynamic coefficient (cell-mapping) continuously averaging the new value estimated for each cell with the previously existing one. To that end, an acceptable tradeoff between the range and discretization step of each dependency variable and available runtime memory must be achieved.

[0135] A measure of the goodness of the LS fit in every linear observation problem considered above can be obtained through the RMS (Root Mean Square) error:

[00043] $\begin{matrix} e = Hz - O & [230] \\ SSE = e^{T} .Math. e .Math. .Math. Sum .Math. .Math. of .Math. .Math. Square .Math. .Math. Errors & [231] \\ RMS = \sqrt{\frac{SSE}{\dim (e)}} .Math. .Math. Root .Math. .Math. Mean .Math. .Math. Square & [232] \end{matrix}$

Linearization of the Equations of Motion (Linear Motion of the CoG)

Linear Motion Equations in BFS:

[0136] [00044] $\begin{matrix} [\begin{matrix} L .Math. .Math. \sin .Math. .Math. + Q .Math. .Math. \sin .Math. .Math. .Math. .Math. \cos .Math. .Math. - D .Math. .Math. \cos .Math. .Math. .Math. .Math. \cos .Math. .Math. \\ - D .Math. .Math. \sin .Math. .Math. - Q .Math. .Math. \cos .Math. .Math. \\ - L .Math. .Math. \cos .Math. .Math. + Q .Math. .Math. \sin .Math. .Math. .Math. .Math. \sin .Math. .Math. - D .Math. .Math. \cos .Math. .Math. .Math. .Math. \sin .Math. .Math. \end{matrix}] + T [\begin{matrix} 1 \\ - \\ - .Math. \end{matrix}] = m [\begin{matrix} a_{1}^{BFS} \\ a_{2}^{BFS} \\ a_{3}^{BFS} \end{matrix}] & [233] \end{matrix}$

Longitudinal Motion with AOS-On-Rudder to Ensure Coordinated Flight:

[00045] $\dot{} 0$ $p r 0$ $L .Math. .Math. \sin .Math. .Math. - D .Math. .Math. \cos .Math. .Math. + T = {ma}_{1}^{BFS} .Math. - L .Math. .Math. \cos .Math. .Math. - D .Math. .Math. \sin .Math. .Math. - T .Math. .Math. .Math. = {ma}_{3}^{BFS}$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S [C_{L} .Math. \sin .Math. .Math. - C_{D} .Math. \cos .Math. .Math.] + W_{MTOW} .Math. .Math. .Math. C_{T} = {ma}_{1}^{BFS} .Math. - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S [C_{L} .Math. \cos .Math. .Math. + C_{D} .Math. \sin .Math. .Math.] - W_{MTOW} .Math. .Math. .Math. C_{T} .Math. .Math. = {ma}_{3}^{BFS}$

[0137] Close enough to the balanced flight condition and, as long as E, is held null, .sub.a and .sub.r are close to null and .sub.h,0 is held constant, the aerodynamic lift and drag coefficients can be approximated by the respective Taylor expansions in the form:

[00046] $\begin{matrix} C_{L} (, 0, M, \hat{\dot{}}, 0, 0, \hat{q}, 0, {.Math.}_{h, 0}, 0, 0, 0) = C_{L} (, M, \hat{\dot{}}, \hat{q}) = C_{L, 0} + C_{L,} .Math. .Math. .Math. + C_{L, M} .Math. .Math. .Math. M + C_{L, \hat{\dot{}}} .Math. \hat{\dot{}} + C_{L, \hat{q}} .Math. \hat{q} = C_{L, 0} + .Math. .Math. C_{L} & [238] \\ C_{D} (, 0, M, \hat{\dot{}}, 0, 0, \hat{q}, 0, {.Math.}_{h, 0}, 0, 0, 0) = C_{D} (, M, \hat{\dot{}}, \hat{q}) = C_{D, 0} + C_{D,} .Math. .Math. .Math. + C_{D, M} .Math. .Math. .Math. M + C_{D, \hat{\dot{}}} .Math. \hat{\dot{}} + C_{D, \hat{q}} .Math. \hat{q} = C_{D, 0} + .Math. .Math. C_{D} & [239] \\ .Math. C_{L, \hat{\dot{}}} .Math. \hat{\dot{}} = C_{L, \hat{\dot{}}} .Math. \frac{\dot{} .Math. .Math. c_{w}}{2 .Math. .Math. v_{TAS}} = C_{L, \hat{\dot{}}} .Math. \frac{c_{w}}{2 .Math. a_{0}} .Math. \frac{\dot{a}}{M} = C_{L, \dot{a}} .Math. \frac{\dot{a}}{M} .Math. .Math. .Math. C_{L, \dot{}} = C_{L, \hat{\dot{}}} .Math. \frac{c_{w}}{2 .Math. a_{0}} & [240] \\ .Math. C_{L, \hat{q}} .Math. \hat{q} = C_{L, \hat{q}} .Math. \frac{{qc}_{w}}{2 .Math. v_{TAS}} = C_{L, \hat{q}} .Math. \frac{c_{w}}{2 .Math. a_{0}} .Math. \frac{q}{M} = C_{L, q} .Math. \frac{q}{M} .Math. .Math. .Math. C_{L, q} = C_{L, \hat{q}} .Math. \frac{c_{w}}{2 .Math. a_{0}} & [241] \\ .Math. C_{D, \hat{\dot{}}} .Math. \hat{\dot{}} = C_{D, \hat{\dot{}}} .Math. \frac{\dot{} .Math. .Math. c_{w}}{2 .Math. v_{TAS}} = C_{D, \hat{\dot{}}} .Math. \frac{c_{w}}{2 .Math. a_{0}} .Math. \frac{\dot{}}{M} = C_{D, \dot{}} .Math. \frac{\dot{}}{M} .Math. .Math. .Math. C_{D, \dot{}} = C_{L, \hat{\dot{}}} .Math. \frac{c_{w}}{2 .Math. a_{0}} & [242] \\ .Math. C_{D, \hat{q}} .Math. \hat{q} = C_{D, \hat{q}} .Math. \frac{{qc}_{w}}{2 .Math. v_{TAS}} = C_{D, \hat{q}} .Math. \frac{c_{w}}{2 .Math. a_{0}} .Math. \frac{q}{M} = C_{D, q} .Math. \frac{q}{M} .Math. .Math. .Math. C_{D, q} = C_{D, \hat{q}} .Math. \frac{c_{w}}{2 .Math. a_{0}} & [243] \\ C_{L} (, M, \hat{\dot{}}, \hat{q}) = C_{L, 0} + C_{L,} .Math. .Math. .Math. + C_{L, M} .Math. .Math. .Math. M + C_{L, \dot{}} .Math. \frac{\dot{}}{M} + C_{L, q} .Math. \frac{q}{M} & [244] \\ C_{D} (, M, \hat{\dot{}}, \hat{q}) = C_{D, 0} + C_{D,} .Math. .Math. .Math. + C_{D, M} .Math. .Math. .Math. M + C_{D, \dot{}} .Math. \frac{\dot{}}{M} + C_{D, q} .Math. \frac{q}{M} & [245] \end{matrix}$

[0138] As for the thrust coefficient, assuming level flight at the trim condition and the fact that the throttle level does not change:

C.sub.T(,,M,.sub.T,0)=C.sub.T,0(,,M.sub.0,.sub.T,0)+C.sub.T,M(,,M.sub.0,.sub.T,0)M[246]

[0139] At the balanced flight condition:

[00047] $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M_{0}^{2} .Math. S [C_{L, 0} .Math. \sin .Math. .Math._{0} - C_{D, 0} .Math. \cos .Math. .Math._{0}] + W_{MTOW} .Math. .Math. .Math. C_{T, 0} = m_{0} .Math. a_{1, 0}^{BFS} .Math. - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M_{0}^{2} .Math. S [C_{L, 0} .Math. \cos .Math. .Math._{0} + C_{D, 0} .Math. \sin .Math. .Math._{0}] - W_{MTOW} .Math. .Math. .Math. C_{T, 0} .Math. .Math. = m_{0} .Math. a_{3, 0}^{BFS}$

[0140] And close enough to the balanced flight condition:

[00048] $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S [C_{L} .Math. \sin .Math. .Math. - C_{D} .Math. \cos .Math. .Math.] + W_{MTOW} .Math. .Math. .Math. C_{T} = m_{0} (a_{1, 0}^{BFS} + .Math. .Math. a_{1}^{BFS}) - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S [C_{L} .Math. \cos .Math. .Math. + C_{D} .Math. \sin .Math. .Math.] - W_{MTOW} .Math. .Math. .Math. C_{T} .Math. .Math. = m_{0} (a_{3, 0}^{BFS} + .Math. .Math. a_{3}^{BFS})$ $.Math. .Math. .Math. a_{1}^{BFS} = a_{1}^{BFS} - a_{1, 0}^{BFS}$ $.Math. .Math. .Math. a_{3}^{BFS} = a_{3}^{BFS} - a_{3, 0}^{BFS}$ $C_{L} .Math. \sin (_{0} +) - C_{D} .Math. \cos (_{0} +) = C_{L} (\sin .Math. .Math._{0} .Math. \cos .Math. .Math. + \cos .Math. .Math._{0} .Math. \sin .Math. .Math.) - C_{D} (\cos .Math. .Math._{0} .Math. \cos .Math. .Math. - \sin .Math. .Math._{0} .Math. \sin .Math. .Math.)$ $C_{L} .Math. \cos (_{0} +) + C_{D} .Math. \sin (_{0} +) = C_{L} (\cos .Math. .Math._{0} .Math. \cos .Math. .Math. - \sin .Math. .Math._{0} .Math. \sin .Math. .Math.) + C_{D} (\sin .Math. .Math._{0} .Math. \cos .Math. .Math. + \cos .Math. .Math._{0} .Math. \sin .Math. .Math.)$ $.Math. 1$ $C_{L} .Math. \sin (_{0} +) - C_{D} .Math. \cos (_{0} +) C_{L} .Math. (\sin .Math. .Math._{0} + \cos .Math. .Math._{0} .Math.) - C_{D} (\cos .Math. .Math._{0} - \sin .Math. .Math._{0} .Math.) = C_{L} .Math. \sin .Math. .Math._{0} - C_{D} .Math. \cos .Math. .Math._{0} + (C_{L} .Math. \cos .Math. .Math._{0} + C_{D} .Math. \sin .Math. .Math._{0}) .Math.$ $C_{L} .Math. \cos (_{0} +) + C_{D} .Math. \sin (_{0} +) C_{L} (\cos .Math. .Math._{0} - \sin .Math. .Math._{0} .Math.) + C_{D} (\sin .Math. .Math._{0} + \cos .Math. .Math._{0} .Math.) == C_{L} .Math. \cos .Math. .Math._{0} + C_{D} .Math. \sin .Math. .Math._{0} - (C_{L} .Math. \sin .Math. .Math._{0} - C_{D} .Math. \cos .Math. .Math._{0}) .Math.$ $C_{L} .Math. \sin (_{0} +) - C_{D} .Math. \cos (_{0} +) + C_{L} .Math. \cos (_{0} +) + C_{D} .Math. \sin (_{0} +) == C_{L} .Math. \sin .Math. .Math._{0} - C_{D} .Math. \cos .Math. .Math._{0} + (C_{L} .Math. \cos .Math. .Math._{0} + C_{D} .Math. \sin .Math. .Math._{0}) .Math. + C_{L} .Math. \cos .Math. .Math._{0} + C_{D} .Math. \sin .Math. .Math._{0} - (C_{L} .Math. \sin .Math. .Math._{0} - C_{D} .Math. \cos .Math. .Math._{0}) .Math.$ $C_{L} .Math. \sin (_{0} +) - C_{D} .Math. \cos (_{0} +) - C_{L} .Math. \cos (_{0} +) - C_{D} .Math. \sin (_{0} +) == C_{L} .Math. \sin .Math. .Math._{0} - C_{D} .Math. \cos .Math. .Math._{0} + (C_{L} .Math. \cos .Math. .Math._{0} + C_{D} .Math. \sin .Math. .Math._{0}) .Math. - C_{L} .Math. \cos .Math. .Math._{0} - C_{D} .Math. \sin .Math. .Math._{0} + (C_{L} .Math. \sin .Math. .Math._{0} - C_{D} .Math. \cos .Math. .Math._{0}) .Math.$ $C_{L} .Math. \sin (_{0} +) - C_{D} .Math. \cos (_{0} +) + C_{L} .Math. \cos (_{0} +) + C_{D} .Math. \sin (_{0} +) == \sin .Math. .Math._{0} .Math. C_{L} (1 -) - \cos .Math. .Math._{0} .Math. C_{D} (1 -) + \cos .Math. .Math._{0} .Math. C_{L} (1 +) + \sin .Math. .Math._{0} .Math. C_{D} (1 +)$ $C_{L} .Math. \sin (_{0} +) - C_{D} .Math. \cos (_{0} +) - C_{L} .Math. \cos (_{0} +) - C_{D} .Math. \sin (_{0} +) == \sin .Math. .Math._{0} .Math. C_{L} (1 +) - \cos .Math. .Math._{0} .Math. C_{D} (1 +) - \cos .Math. .Math._{0} .Math. C_{L} (1 -) - \sin .Math. .Math._{0} .Math. C_{D} (1 -)$ $.Math. 1$ $C_{L} .Math. \sin (_{0} +) - C_{D} .Math. \cos (_{0} +) + C_{L} .Math. \cos (_{0} +) + C_{D} .Math. \sin (_{0} +) \sin .Math. .Math._{0} .Math. C_{L} - \cos .Math. .Math._{0} .Math. C_{D} + \cos .Math. .Math._{0} .Math. C_{L} + \sin .Math. .Math._{0} .Math. C_{D} .Math.$ $C_{L} .Math. \sin (_{0} +) - C_{D} .Math. \cos (_{0} +) - C_{L} .Math. \cos (_{0} +) - C_{D} .Math. \sin (_{0} +) \sin .Math. .Math._{0} .Math. C_{l} - \cos .Math. .Math._{0} .Math. C_{D} - \cos .Math. .Math._{0} .Math. C_{L} - \sin .Math. .Math._{0} .Math. C_{D}$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S (\sin .Math. .Math._{0} .Math. C_{L} - \cos .Math. .Math._{0} .Math. C_{D} + \cos .Math. .Math._{0} .Math. C_{L} + \sin .Math. .Math._{0} .Math. C_{D}) + W_{MTOW} .Math. .Math. .Math. C_{T} (1 + .Math.) = m_{0} (a_{1, 0}^{BFS} - a_{3, 0}^{BFS} + .Math. .Math. a_{1}^{BFS} - .Math. .Math. a_{3}^{BFS})$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S (\sin .Math. .Math._{0} .Math. C_{L} - \cos .Math. .Math._{0} .Math. C_{D} - \cos .Math. .Math._{0} .Math. C_{L} - \sin .Math. .Math._{0} .Math. C_{D}) + W_{MTOW} .Math. .Math. .Math. C_{T} (1 - .Math.) = m_{0} (a_{1, 0}^{BFS} + a_{3, 0}^{BFS} + .Math. .Math. a_{1}^{BFS} + .Math. .Math. a_{3}^{BFS})$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S (\sin .Math. .Math._{0} .Math. C_{L} - \cos .Math. .Math._{0} .Math. C_{D}) + W_{MTOW} .Math. .Math. .Math. C_{T} = m_{0} (a_{1, 0}^{BFS} + .Math. .Math. a_{1}^{BFS}) - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S (\cos .Math. .Math._{0} .Math. C_{L} + \sin .Math. .Math._{0} .Math. C_{D}) - W_{MTOW} .Math. .Math. .Math. C_{T} .Math. .Math. = m_{0} (a_{3, 0}^{BFS} + .Math. .Math. a_{3}^{BFS})$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S [\sin .Math. .Math._{0} (C_{L, 0} + .Math. .Math. C_{L}) - \cos .Math. .Math._{0} (C_{D, 0} + .Math. .Math. C_{D})] + W_{MTOW} .Math. .Math. .Math. (C_{T, 0} + C_{T, M} .Math. .Math. .Math. M) = m_{0} (a_{1, 0}^{BFS} + .Math. .Math. a_{1}^{BFS}) - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S [\cos .Math. .Math._{0} (C_{L, 0} + .Math. .Math. C_{L}) + \sin .Math. .Math._{0} (C_{D, 0} + .Math. .Math. C_{D})] - W_{MTOW} .Math. .Math. .Math. (C_{T, 0} + C_{T, M} .Math. .Math. .Math. M) .Math. .Math. = m_{0} (a_{3, 0}^{BFS} + .Math. .Math. a_{3}^{BFS})$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S (\sin .Math. .Math._{0} .Math. C_{L, 0} - \cos .Math. .Math._{0} .Math. C_{D, 0}) + \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S (\sin .Math. .Math._{0} .Math. .Math. .Math. C_{L} - \cos .Math. .Math._{0} .Math. .Math. .Math. C_{D}) + W_{MTOW} .Math. (C_{T, 0} + C_{T, M} .Math. .Math. .Math. M) = m_{0} (a_{1, 0}^{BFS} + .Math. .Math. a_{1}^{BFS}) - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S (\cos .Math. .Math._{0} .Math. C_{L, 0} + \sin .Math. .Math._{0} .Math. C_{D, 0}) - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S (\cos .Math. .Math._{0} .Math. .Math. .Math. C_{L} + \sin .Math. .Math._{0} .Math. .Math. .Math. C_{D}) - W_{MTOW} .Math. (C_{T, 0} + C_{T, M} .Math. .Math. .Math. M) .Math. .Math. = m_{0} (a_{3, 0}^{BFS} + .Math. .Math. a_{3}^{BFS})$ $M^{2} = {(M_{0} + .Math. .Math. M)}^{2} = M_{0}^{2} + 2 .Math. .Math. M_{0} .Math. .Math. .Math. M + .Math. .Math. M^{2} = M_{0}^{2} [1 + 2 .Math. \frac{.Math. .Math. M}{M_{0}} + {(\frac{.Math. .Math. M}{M_{0}})}^{2}]$ $.Math. \frac{.Math. .Math. M}{M_{0}} 1$ $.Math. M^{2} M_{0}^{2} (1 + 2 .Math. \frac{.Math. .Math. M}{M_{0}})$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. S (M_{0}^{2} + 2 .Math. .Math. M_{0} .Math. .Math. .Math. M) .Math. (\sin .Math. .Math._{0} .Math. C_{L, 0} - \cos .Math. .Math._{0} .Math. C_{D, 0}) + \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {SM}^{2} (\sin .Math. .Math._{0} .Math. .Math. .Math. C_{L} - \cos .Math. .Math._{0} .Math. .Math. .Math. C_{D}) + W_{MTOW} .Math. (C_{T, 0} + C_{T, M} .Math. .Math. .Math. M) == m_{0} (a_{1, 0}^{BFS} + .Math. .Math. a_{1}^{BFS}) - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. S (M_{0}^{2} + 2 .Math. .Math. M_{0} .Math. .Math. .Math. M) .Math. (\cos .Math. .Math._{0} .Math. C_{L, 0} + \sin .Math. .Math._{0} .Math. C_{D, 0}) - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {SM}^{2} (\cos .Math. .Math._{0} .Math. .Math. .Math. C_{L} + \sin .Math. .Math._{0} .Math. .Math. .Math. C_{D}) - W_{MTOW} .Math. (C_{T, 0} + C_{T, M} .Math. .Math. .Math. M) .Math. .Math. == m_{0} (a_{3, 0}^{BFS} + .Math. .Math. a_{3}^{BFS})$

[0141] Bearing in mind [247]:

[00049] $.Math. .Math. p_{0} .Math. .Math. .Math. {SM}_{0} .Math. .Math. .Math. M (\sin .Math. .Math._{0} .Math. C_{L, 0} - \cos .Math. .Math._{0} .Math. C_{D, 0}) + \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {SM}^{2} (\sin .Math. .Math._{0} .Math. .Math. .Math. C_{L} - \cos .Math. .Math._{0} .Math. .Math. .Math. C_{D}) + W_{MTOW} .Math. .Math. .Math. C_{T, m} .Math. .Math. .Math. M = m_{0} .Math. .Math. .Math. a_{1}^{BFS} - .Math. .Math. p_{0} .Math. .Math. .Math. {SM}_{0} .Math. .Math. .Math. M (\cos .Math. .Math._{0} .Math. C_{L, 0} + \sin .Math. .Math._{0} .Math. C_{D, 0}) - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {SM}^{2} (\cos .Math. .Math._{0} .Math. .Math. .Math. C_{L} + \sin .Math. .Math._{0} .Math. .Math. .Math. C_{D}) - W_{MTOW} .Math. .Math. .Math. C_{T, m} .Math. .Math. .Math. M .Math. .Math. .Math. = m_{0} .Math. .Math. .Math. a_{3}^{BFS}$ $.Math. .Math. p_{0} .Math. .Math. .Math. {SM}_{0} (\sin .Math. .Math._{0} .Math. C_{L, 0} - \cos .Math. .Math._{0} .Math. C_{D, 0}) .Math. .Math. .Math. M ++ .Math. \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S [\sin .Math. .Math._{0} .Math. (C_{L,} .Math. + C_{L, M} .Math. .Math. .Math. M + C_{L, \dot{}} .Math. \frac{\dot{}}{M} + C_{L, q} .Math. \frac{q}{M}) - \cos .Math. .Math._{0} (C_{D,} .Math. + C_{D, M} .Math. .Math. .Math. M + C_{D, \dot{}} .Math. \frac{\dot{}}{M} + C_{D, q} .Math. \frac{q}{M})] ++ .Math. W_{MTOW} .Math. .Math. .Math. C_{T, M} .Math. .Math. .Math. M = m_{0} .Math. .Math. .Math. a_{1}^{BFS} - .Math. .Math. p_{0} .Math. .Math. .Math. {SM}_{0} (\cos .Math. .Math._{0} .Math. C_{L, 0} + \sin .Math. .Math._{0} .Math. C_{D, 0}) .Math. .Math. .Math. M -- .Math. \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S [\cos .Math. .Math._{0} (C_{L,} .Math. + C_{L, M} .Math. .Math. .Math. M + C_{L, \dot{}} .Math. \frac{\dot{}}{M} + C_{L, q} .Math. \frac{q}{M}) + \sin .Math. .Math._{0} (C_{D,} .Math. + C_{D, M} .Math. .Math. .Math. M + C_{D, \dot{}} .Math. \frac{\dot{}}{M} + C_{D, q} .Math. \frac{q}{M})] - W_{MTOW} .Math. .Math. .Math. C_{T, M} .Math. .Math. .Math. M .Math. .Math. .Math. = m_{0} .Math. .Math. .Math. a_{3}^{BFS}$ $M^{2} .Math. .Math. .Math. M = {(M_{0} + .Math. .Math. M)}^{2} .Math. .Math. .Math. M = (M_{0}^{2} + 2 .Math. .Math. M_{0} .Math. .Math. .Math. M + .Math. .Math. M^{2}) .Math. .Math. .Math. M = M_{0}^{2} [1 + 2 .Math. \frac{.Math. .Math. M}{M_{0}} + {(\frac{.Math. .Math. M}{M_{0}})}^{2}] .Math. .Math. .Math. M M_{0}^{2} (1 + 2 .Math. \frac{.Math. .Math. M}{M_{0}}) .Math. .Math. .Math. M$ $C_{X, M} = \frac{2}{M_{0}} .Math. (\sin .Math. .Math._{0} .Math. C_{L, 0} - \cos .Math. .Math._{0} .Math. C_{D, 0}) + (\sin .Math. .Math._{0} .Math. C_{L, M} - \cos .Math. .Math._{0} .Math. C_{D, M}) + \frac{W_{MTOW} .Math. C_{T, M}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. {SM}_{0}^{2}}$ $C_{Z, M} = \frac{2}{M_{0}} .Math. (\cos .Math. .Math._{0} .Math. C_{L, 0} + \sin .Math. .Math._{0} .Math. C_{D, 0}) + (\cos .Math. .Math._{0} .Math. C_{L, M} + \sin .Math. .Math._{0} .Math. C_{D, M}) + \frac{W_{MTOW} .Math. C_{T, M}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. {SM}_{0}^{2}} .Math. .Math.$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {SM}_{0}^{2} .Math. C_{X, M} .Math. .Math. .Math. M ++ .Math. \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. S [M^{2} (\sin .Math. .Math._{0} .Math. C_{L,} - \cos .Math. .Math._{0} .Math. C_{D,}) .Math. + 2 .Math. .Math. M_{0} (\sin .Math. .Math._{0} .Math. C_{L, M} - \cos .Math. .Math._{0} .Math. C_{D, M}) .Math. .Math. .Math. M^{2} + M (\sin .Math. .Math._{0} .Math. C_{L, \dot{}} - \cos .Math. .Math._{0} .Math. C_{D, \dot{}}) .Math. \dot{} + M (\sin .Math. .Math._{0} .Math. C_{L, q} - \cos .Math. .Math._{0} .Math. C_{D, q}) .Math. q] == m_{0} .Math. .Math. .Math. a_{1}^{BFS} - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {SM}_{0}^{2} .Math. C_{Z, M} .Math. .Math. .Math. M -- .Math. \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. S [M^{2} (\cos .Math. .Math._{0} .Math. C_{L,} + \sin .Math. .Math._{0} .Math. C_{D,}) .Math. + 2 .Math. .Math. M_{0} (\cos .Math. .Math._{0} .Math. C_{L, M} + \sin .Math. .Math._{0} .Math. C_{D, M}) .Math. .Math. .Math. M^{2} + M (\cos .Math. .Math._{0} .Math. C_{L, \dot{}} + \sin .Math. .Math._{0} .Math. C_{D, \dot{}}) .Math. \dot{} + M (\cos .Math. .Math._{0} .Math. C_{L, q} + \sin .Math. .Math._{0} .Math. C_{D, q}) .Math. q] == m_{0} .Math. .Math. .Math. a_{3}^{BFS}$ $.Math. C_{X,} = \sin .Math. .Math._{0} .Math. C_{L,} - \cos .Math. .Math._{0} .Math. C_{D,}$ $.Math. C_{Z,} = \cos .Math. .Math._{0} .Math. C_{L,} + \sin .Math. .Math._{0} .Math. C_{D,}$ $.Math. C_{X, .Math. .Math. M} = \sin .Math. .Math._{0} .Math. C_{L, M} - \cos .Math. .Math._{0} .Math. C_{D, M}$ $.Math. C_{Z, .Math. .Math. M} = \cos .Math. .Math._{0} .Math. C_{L, M} + \sin .Math. .Math._{0} .Math. C_{D, M}$ $.Math. C_{X, \dot{}} = \sin .Math. .Math._{0} .Math. C_{L, \dot{}} - \cos .Math. .Math._{0} .Math. C_{D, \dot{}}$ $.Math. C_{Z, \dot{}} = \cos .Math. .Math._{0} .Math. C_{L,} + \sin .Math. .Math._{0} .Math. C_{D, \dot{}}$ $.Math. C_{X, q} = \sin .Math. .Math._{0} .Math. C_{L, q} - \cos .Math. .Math._{0} .Math. C_{D, q}$ $.Math. C_{Z, q} = \cos .Math. .Math._{0} .Math. C_{L, q} + \sin .Math. .Math._{0} .Math. C_{D, q}$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. S (M_{0}^{2} .Math. C_{X, M} .Math. .Math. .Math. M + M^{2} .Math. C_{X,} .Math. + 2 .Math. .Math. M_{0} .Math. C_{X, .Math. .Math. M} .Math. .Math. .Math. M^{2} + {MC}_{X, \dot{}} .Math. \dot{} + {MC}_{X, q} .Math. q) = m_{0} .Math. .Math. .Math. a_{1}^{BFS} - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. S (M_{0}^{2} .Math. C_{Z, M} .Math. .Math. .Math. M + M^{2} .Math. C_{Z,} .Math. + 2 .Math. .Math. M_{0} .Math. C_{Z, .Math. .Math. M} .Math. .Math. .Math. M^{2} + {MC}_{Z, \dot{}} .Math. \dot{} + {MC}_{Z, q} .Math. q) = m_{0} .Math. .Math. .Math. a_{3}^{BFS}$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. S [M_{0}^{2} .Math. C_{X, M} .Math. .Math. .Math. M + M_{0}^{2} (1 + 2 .Math. \frac{.Math. .Math. M}{M_{0}}) .Math. C_{X,} .Math. + 2 .Math. .Math. M_{0} .Math. C_{X, .Math. .Math. M} .Math. .Math. .Math. M^{2} + (M_{0} + .Math. .Math. M) .Math. C_{X, \dot{}} .Math. \dot{} + (M_{0} + .Math. .Math. M) .Math. C_{X, q} .Math. q] = m_{0} .Math. .Math. .Math. a_{1}^{BFS} - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. S [M_{0}^{2} .Math. C_{Z, M} .Math. .Math. .Math. M + M_{0}^{2} (1 + 2 .Math. \frac{.Math. .Math. M}{M_{0}}) .Math. C_{Z,} .Math. + 2 .Math. .Math. M_{0} .Math. C_{Z, .Math. .Math. M} .Math. .Math. .Math. M^{2} + (M_{0} + .Math. .Math. M) .Math. C_{Z, \dot{}} .Math. \dot{} + (M_{0} + .Math. .Math. M) .Math. C_{Z, q} .Math. q] = m_{0} .Math. .Math. .Math. a_{3}^{BFS}$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {SM}_{0}^{2} [C_{X,} .Math. + C_{X, M} .Math. .Math. .Math. M + \frac{C_{X, \dot{}}}{M_{0}} .Math. \dot{} + \frac{C_{X, q}}{M_{0}} .Math. q + (2 .Math. \frac{C_{X,}}{M_{0}} .Math. + 2 .Math. \frac{C_{X, .Math. .Math. M}}{M_{0}} .Math. .Math. .Math. M + \frac{C_{X, \dot{}}}{M_{0}^{2}} .Math. \dot{} + \frac{C_{X, q}}{M_{0}^{2}} .Math. q) .Math. .Math. .Math. M] = m_{0} .Math. .Math. .Math. a_{1}^{BFS} - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {SM}_{0}^{2} [C_{Z,} .Math. + C_{Z, M} .Math. .Math. .Math. M + \frac{C_{Z, \dot{}}}{M_{0}} .Math. \dot{} + \frac{C_{Z, q}}{M_{0}} .Math. q + (2 .Math. \frac{C_{Z,}}{M_{0}} .Math. + 2 .Math. \frac{C_{Z, .Math. .Math. M}}{M_{0}} .Math. .Math. .Math. M + \frac{C_{Z, \dot{}}}{M_{0}^{2}} .Math. \dot{} + \frac{C_{Z, q}}{M_{0}^{2}} .Math. q) .Math. .Math. .Math. M] = m_{0} .Math. .Math. .Math. a_{3}^{BFS}$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {Sc}_{Z, \dot{}} (M_{0}^{2} .Math. C_{X, M} .Math. .Math. .Math. M + M^{2} .Math. C_{X,} .Math. + 2 .Math. .Math. M_{0} .Math. C_{X, .Math. .Math. M} .Math. .Math. .Math. M^{2} + {MC}_{X, \dot{}} .Math. \dot{} + {MC}_{X, q} .Math. q) = m_{0} .Math. C_{Z, \dot{}} .Math. .Math. .Math. a_{1}^{BFS} - \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {SC}_{X, \dot{}} (M_{0}^{2} .Math. C_{Z, M} .Math. .Math. .Math. M + M^{2} .Math. C_{Z,} .Math. + 2 .Math. .Math. M_{0} .Math. C_{Z, .Math. .Math. M} .Math. .Math. .Math. M^{2} + {MC}_{Z, \dot{}} .Math. \dot{} + {MC}_{Z, q} .Math. q) = m_{0} .Math. C_{X, \dot{}} .Math. .Math. .Math. a_{3}^{BFS}$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. S [M_{0}^{2} (C_{Z, \dot{}} .Math. C_{X, M} - C_{X, \dot{}} .Math. C_{Z, M}) .Math. .Math. .Math. M + M^{2} (C_{Z, \dot{}} .Math. C_{X,} - C_{X, \dot{}} .Math. C_{Z,}) .Math. + 2 .Math. .Math. M_{0} (C_{Z, \dot{}} .Math. C_{X, .Math. .Math. M} - C_{X, \dot{}} .Math. C_{Z, .Math. .Math. M}) .Math. .Math. .Math. M^{2} + M (C_{Z, \dot{}} .Math. C_{X, q} - C_{X, \dot{}} .Math. C_{Z, q}) .Math. q] == m_{0} (C_{Z, \dot{}} .Math. .Math. .Math. a_{1}^{BFS} + C_{X, \dot{}} .Math. .Math. .Math. a_{3}^{BFS})$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. S [M_{0}^{2} (C_{Z, \dot{}} .Math. C_{X, M} + C_{X, \dot{}} .Math. C_{Z, M}) .Math. .Math. .Math. M + M^{2} (C_{Z, \dot{}} .Math. C_{X,} + C_{X, \dot{}} .Math. C_{Z,}) .Math. + 2 .Math. .Math. M_{0} (C_{Z, \dot{}} .Math. C_{X, .Math. .Math. M} + C_{X, \dot{}} .Math. C_{Z, .Math. .Math. M}) .Math. .Math. .Math. M^{2} + 2 .Math. .Math. {MC}_{X, \dot{}} .Math. C_{Z, \dot{}} .Math. \dot{} + M (C_{Z, \dot{}} .Math. C_{X, q} + C_{X, \dot{}} .Math. C_{Z, q}) .Math. q] == m_{0} (C_{Z, \dot{}} .Math. .Math. .Math. a_{1}^{BFS} - C_{X, \dot{}} .Math. .Math. .Math. a_{3}^{BFS})$ $.Math. C_{1} = C_{Z, \dot{}} .Math. C_{X, M} - C_{X, \dot{}} .Math. C_{Z, M}$ $.Math. C_{2} = C_{Z, \dot{}} .Math. C_{X,} - C_{X, \dot{}} .Math. C_{Z,}$ $.Math. C_{3} = C_{Z, \dot{}} .Math. C_{X, .Math. .Math. M} - C_{X, \dot{}} .Math. C_{Z, .Math. .Math. M}$ $.Math. C_{4} = C_{Z, \dot{}} .Math. C_{X, q} - C_{X, \dot{}} .Math. C_{Z, q}$ $.Math. D_{1} = C_{Z, \dot{}} .Math. C_{X, M} + C_{X, \dot{}} .Math. C_{Z, M}$ $.Math. D_{2} = C_{Z, \dot{}} .Math. C_{X,} + C_{X, \dot{}} .Math. C_{Z,}$ $.Math. D_{3} = C_{Z, \dot{}} .Math. C_{X, .Math. .Math. M} + C_{X, \dot{}} .Math. C_{Z, .Math. .Math. M}$ $.Math. D_{4} = C_{Z, \dot{}} .Math. C_{X, q} + C_{X, \dot{}} .Math. C_{Z, q}$ $.Math. D_{5} = 2 .Math. .Math. C_{X, \dot{}} .Math. C_{Z, \dot{}}$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. S (M_{0}^{2} .Math. C_{1} .Math. .Math. .Math. M + M^{2} .Math. C_{2} .Math. + 2 .Math. .Math. M_{0} .Math. C_{3} .Math. .Math. .Math. M^{2} + {MC}_{4} .Math. q) = m_{0} (C_{Z, \dot{}} .Math. .Math. .Math. a_{1}^{BFS} + C_{X, \dot{}} .Math. .Math. .Math. a_{3}^{BFS})$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. S (M_{0}^{2} .Math. D_{1} .Math. .Math. .Math. M + M^{2} .Math. D_{2} .Math. + 2 .Math. .Math. M_{0} .Math. D_{3} .Math. .Math. .Math. M^{2} + {MD}_{4} .Math. \dot{} + {MD}_{5} .Math. q) = m_{0} (C_{Z, \dot{}} .Math. .Math. .Math. a_{1}^{BFS} - C_{X, \dot{}} .Math. .Math. .Math. a_{3}^{BFS})$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {SD}_{3} (M_{0}^{2} .Math. C_{1} .Math. .Math. .Math. M + M^{2} .Math. C_{2} .Math. + 2 .Math. .Math. M_{0} .Math. C_{3} .Math. .Math. .Math. M^{2} + {MC}_{4} .Math. q) = m_{0} .Math. D_{3} (C_{Z, \dot{}} .Math. .Math. .Math. a_{1}^{BFS} + C_{X, \dot{}} .Math. .Math. .Math. a_{3}^{BFS})$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {SC}_{3} (M_{0}^{2} .Math. D_{1} .Math. .Math. .Math. M + M^{2} .Math. D_{2} .Math. + 2 .Math. .Math. M_{0} .Math. D_{3} .Math. .Math. .Math. M^{2} + {MD}_{4} .Math. q + {MD}_{5} .Math. \dot{}) = m_{0} .Math. C_{3} (C_{Z, \dot{}} .Math. .Math. .Math. a_{1}^{BFS} - C_{X, \dot{}} .Math. .Math. .Math. a_{3}^{BFS})$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. S [M_{0}^{2} (C_{1} .Math. D_{3} - C_{3} .Math. D_{1}) .Math. .Math. .Math. M + M^{2} (C_{2} .Math. D_{3} - C_{3} .Math. D_{2}) .Math. + M (C_{4} .Math. D_{3} - C_{3} .Math. D_{4}) .Math. q - {MC}_{3} .Math. D_{5} .Math. \dot{}] = m_{0} [.Math. .Math. a_{1}^{BFS} .Math. C_{Z, \dot{}} (D_{3} - C_{3}) + .Math. .Math. a_{3}^{BFS} .Math. C_{X, \dot{}} (D_{3} + C_{3})]$ $.Math. F_{1} = C_{1} .Math. D_{3} - C_{3} .Math. D_{1}$ $.Math. F_{2} = C_{2} .Math. D_{3} - C_{3} .Math. D_{2}$ $.Math. F_{3} = C_{4} .Math. D_{3} - C_{3} .Math. D_{4}$ $.Math. F_{4} = C_{3} .Math. D_{5}$ $.Math. F_{5} = C_{Z, \dot{}} (D_{3} - C_{3}) = 2 .Math. .Math. C_{Z, \dot{}} .Math. C_{X, \dot{}} .Math. C_{Z, .Math. .Math. M}$ $.Math. F_{6} = C_{X, \dot{}} (D_{3} + C_{3}) = 2 .Math. .Math. C_{X, \dot{}} .Math. C_{Z, \dot{}} .Math. C_{X, .Math. .Math. M}$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. S (M_{0}^{2} .Math. F_{1} .Math. .Math. .Math. M + M^{2} .Math. F_{2} .Math. + {MF}_{3} .Math. q - {MF}_{4} .Math. \dot{}) = m_{0} (F_{5} .Math. .Math. .Math. a_{1}^{BFS} + F_{6} .Math. .Math. .Math. a_{3}^{BFS})$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. {SMF}_{4} .Math. \dot{} .Math. \frac{C_{P, \dot{}}}{C_{P, \dot{}}} = \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. S (M_{0}^{2} .Math. F_{1} .Math. .Math. .Math. M + M^{2} .Math. F_{2} .Math. + {MF}_{3} .Math. q) = m_{0} (F_{5} .Math. .Math. .Math. a_{1}^{BFS} + F_{6} .Math. .Math. .Math. a_{3}^{BFS}) = \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. S [M_{0}^{2} .Math. F_{1} .Math. .Math. .Math. M + M^{2} .Math. F_{2} .Math. + {MF}_{3} .Math. q - \frac{m_{0}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. S} .Math. (F_{5} .Math. .Math. .Math. a_{1}^{BFS} + F_{6} .Math. .Math. .Math. a_{3}^{BFS})]$ $M .Math. .Math. \dot{} = \frac{1}{F_{4}} [M_{0}^{2} .Math. F_{1} .Math. .Math. .Math. M + M^{2} .Math. F_{2} .Math. + {MF}_{3} .Math. q - \frac{m_{0}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. S} .Math. (F_{5} .Math. .Math. .Math. a_{1}^{BFS} + F_{6} .Math. .Math. .Math. a_{3}^{BFS})]$

[0142] Lateral-directional motion with pitch rate-on-elevator to ensure q{dot over ()}0:

D sin Q cos T=ma.sub.2.sup.BFS[287]

Close enough to the balanced flight condition and, as E, is held close to null, .sub.a and .sub.r are held null and .sub.h,0 is held constant, the aerodynamic drag and side-force coefficients can be approximated by the respective Taylor expansions in the form:

[00050] $\begin{matrix} C_{D} (,, M, 0, \hat{\dot{}}, \hat{p}, 0, \hat{r}, {.Math.}_{h, 0}, 0, 0, 0) = C_{D} (,, M, \hat{\dot{}}, \hat{p}, \hat{r}) = C_{D, 0} + C_{D,} .Math. .Math. .Math. + C_{D,} .Math. + C_{D, M} .Math. .Math. .Math. M + C_{D, \hat{\dot{}}} .Math. \hat{\dot{}} + C_{D, \hat{p}} .Math. \hat{p} + C_{D, \hat{r}} .Math. \hat{r} = C_{D, 0} + .Math. .Math. C_{D} & [288] \\ C_{D} (,, M, 0, \hat{\dot{}}, \hat{p}, 0, \hat{r}, {.Math.}_{h, 0}, 0, 0, 0) = C_{Q} (,, M, \hat{\dot{}}, \hat{p}, \hat{r}) = C_{Q, 0} + C_{Q,} .Math. .Math. .Math. + C_{Q,} .Math. .Math. .Math. + C_{Q, M} .Math. .Math. .Math. M + C_{Q, \hat{\dot{}}} .Math. \hat{\dot{}} + C_{Q, \hat{p}} .Math. \hat{p} + C_{Q, \hat{r}} .Math. \hat{r} = C_{Q, 0} + .Math. .Math. C_{Q} & [289] \\ .Math. C_{D, \hat{\dot{}}} .Math. \hat{\dot{}} = C_{D, \hat{\dot{}}} .Math. \frac{\dot{} .Math. .Math. b_{w}}{2 .Math. .Math. v_{TAS}} = C_{D, \hat{\dot{}}} .Math. \frac{b_{w}}{2 .Math. a_{0}} .Math. \frac{\dot{}}{M} = C_{D, \dot{}} .Math. \frac{\dot{}}{M} .Math. .Math. .Math. C_{D, \dot{}} = C_{D, \hat{\dot{}}} .Math. \frac{b_{w}}{2 .Math. a_{0}} & [290] \\ .Math. C_{Q, \hat{\dot{}}} .Math. \hat{\dot{}} = C_{Q, \hat{\dot{}}} .Math. \frac{\dot{} .Math. .Math. b_{w}}{2 .Math. v_{TAS}} = C_{Q, \hat{\dot{}}} .Math. \frac{b_{w}}{2 .Math. a_{0}} .Math. \frac{\dot{}}{M} = C_{Q, \dot{}} .Math. \frac{\dot{}}{M} .Math. .Math. .Math. C_{Q, \dot{}} = C_{Q, \hat{\dot{}}} .Math. \frac{b_{w}}{2 .Math. a_{0}} & [291] \\ .Math. C_{D, \hat{p}} .Math. \hat{p} = C_{D, \hat{p}} .Math. \frac{q .Math. .Math. b_{w}}{2 .Math. v_{TAS}} = C_{D, \hat{p}} .Math. \frac{b_{w}}{2 .Math. a_{0}} .Math. \frac{p}{M} = C_{D, p} .Math. \frac{p}{M} .Math. .Math. .Math. C_{D, p} = C_{D, \hat{p}} .Math. \frac{b_{w}}{2 .Math. a_{0}} & [292] \\ .Math. C_{D, \hat{r}} .Math. \hat{r} = C_{D, \hat{r}} .Math. \frac{{qb}_{w}}{2 .Math. v_{TAS}} = C_{D, \hat{r}} .Math. \frac{b_{w}}{2 .Math. a_{0}} .Math. \frac{r}{M} = C_{D, r} .Math. \frac{r}{M} .Math. .Math. .Math. C_{D, r} = C_{D, \hat{r}} .Math. \frac{b_{w}}{2 .Math. a_{0}} & [293] \\ .Math. C_{Q, \hat{p}} .Math. \hat{p} = C_{Q, \hat{p}} .Math. \frac{q .Math. .Math. b_{w}}{2 .Math. v_{TAS}} = C_{Q, \hat{p}} .Math. \frac{b_{w}}{2 .Math. a_{0}} .Math. \frac{p}{M} = C_{Q, p} .Math. \frac{p}{M} .Math. .Math. .Math. C_{Q, p} = C_{Q, \hat{p}} .Math. \frac{b_{w}}{2 .Math. a_{0}} & [294] \\ .Math. C_{Q, \hat{r}} .Math. \hat{r} = C_{Q, \hat{r}} .Math. \frac{q .Math. .Math. b_{w}}{2 .Math. v_{TAS}} = C_{Q, \hat{r}} .Math. \frac{b_{w}}{2 .Math. a_{0}} .Math. \frac{r}{M} = C_{Q, r} .Math. \frac{r}{M} .Math. .Math. .Math. C_{Q, r} = C_{Q, \hat{r}} .Math. \frac{b_{w}}{2 .Math. a_{0}} & [295] \\ C_{D} (,, M, \hat{\dot{}}, \hat{p}, \hat{r}) = C_{D, 0} + C_{D,} .Math. .Math. .Math. + C_{D,} .Math. + C_{D, M} .Math. .Math. .Math. M + C_{D, \dot{}} .Math. \frac{\dot{}}{M} + C_{D, p} .Math. \frac{p}{M} + C_{D, r} .Math. \frac{r}{m} & [296] \\ C_{Q} (,, M, \hat{\dot{}}, \hat{p}, \hat{r}) = C_{Q, 0} + C_{Q,} .Math. .Math. .Math. + C_{Q,} .Math. + C_{Q, M} .Math. .Math. .Math. M + C_{Q, \dot{}} .Math. \frac{\dot{}}{M} + C_{Q, p} .Math. \frac{p}{M} + C_{Q, r} .Math. \frac{r}{m} & [297] \end{matrix}$

[0143] As for the thrust coefficient, expression [246] still applies.

[0144] At the balanced flight condition:

p.sub.0M.sub.0.sup.2SC.sub.Q,0W.sub.MTOWC.sub.T,0=m.sub.0a.sub.2,0.sup.BFS[298]

[0145] And close enough to the balanced flight condition:

[00051] $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S (- C_{D} .Math. \sin .Math. .Math. - C_{Q} .Math. \cos .Math. .Math.) - W_{MTOW} .Math. .Math. .Math. C_{T} .Math. = m_{0} (a_{2, 0}^{BFS} + .Math. .Math. a_{2}^{BFS})$ $.Math. .Math. .Math. a_{2}^{BFS} = a_{2}^{BFS} - a_{2, 0}^{BFS}$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S [- (C_{D, 0} + .Math. .Math. C_{D}) .Math. \sin .Math. .Math. - (C_{Q, 0} + .Math. .Math. C_{Q}) .Math. \cos .Math. .Math.] - W_{MTOW} .Math. (C_{T, 0} + C_{T, M} .Math. .Math. .Math. M) .Math. = m_{0} (a_{2, 0}^{BFS} + .Math. .Math. a_{2}^{BFS})$

[0146] Bearing in mind [298]:

[00052] $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S (- .Math. .Math. C_{D} .Math. \sin .Math. .Math. - .Math. .Math. C_{Q} .Math. \cos .Math. .Math.) - W_{MTOW} .Math. .Math. .Math. C_{T, M} .Math. .Math. .Math. M = m_{0} .Math. .Math. .Math. a_{2}^{BFS}$ $.Math. 1$ $.Math. \frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S (- .Math. .Math. C_{D} .Math. - .Math. .Math. C_{Q}) - W_{MTOW} .Math. .Math. .Math. C_{T, M} .Math. .Math. .Math. M = m_{0} .Math. .Math. .Math. a_{2}^{BFS}$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S [- (C_{D,} .Math. + C_{D,} .Math. + C_{D, M} .Math. .Math. .Math. M + C_{D, \dot{}} .Math. \frac{\dot{}}{M} + C_{D, p} .Math. \frac{p}{M} + C_{D, r} .Math. \frac{r}{M}) .Math. - (C_{Q,} .Math. + C_{Q,} .Math. + C_{Q, M} .Math. .Math. .Math. M + C_{Q, \dot{}} .Math. \frac{\dot{}}{M} + C_{Q, p} .Math. \frac{p}{M} + C_{Q, r} .Math. \frac{r}{M})] - W_{MTOW} .Math. .Math. .Math. C_{T, M} .Math. .Math. .Math. M .Math. .Math. = m_{0} .Math. .Math. .Math. a_{2}^{BFS}$ $\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. .Math. .Math. M^{2} .Math. S [- \frac{\dot{}}{M} .Math. (C_{D, \dot{}} .Math. + C_{Q, \dot{}}) - (C_{D,} .Math. + C_{D,} .Math. + C_{D, M} .Math. .Math. .Math. M + C_{D, p} .Math. \frac{p}{M} + C_{D, r} .Math. \frac{r}{M}) .Math. - (C_{Q,} .Math. + C_{Q,} .Math. + C_{Q, M} .Math. .Math. .Math. M + C_{Q, p} .Math. \frac{p}{M} + C_{Q, r} .Math. \frac{r}{M})] - W_{MTOW} .Math. .Math. .Math. C_{T, M} .Math. .Math. .Math. M .Math. .Math. = m_{0} .Math. .Math. .Math. a_{2}^{BFS}$ $M .Math. .Math. \dot{} = - \frac{1}{C_{D, \dot{}} .Math. + C_{Q, \dot{}}} .Math. [M^{2} [(C_{D,} .Math. + C_{D,} .Math. + C_{D, M} .Math. .Math. .Math. M + C_{D, p} .Math. \frac{p}{M} + C_{D, r} .Math. \frac{r}{M}) .Math. + (C_{Q,} .Math. + C_{Q,} .Math. + C_{Q, M} .Math. .Math. .Math. M + C_{Q, p} .Math. \frac{p}{M} + C_{Q, r} .Math. \frac{r}{M})] + \frac{W_{MTOW} .Math. C_{T, M}}{\frac{1}{2} .Math. .Math. .Math. p_{o} .Math. S} .Math. .Math. .Math. M + \frac{m_{0}}{\frac{1}{2} .Math. .Math. .Math. p_{o} .Math.} .Math. .Math. .Math. a_{2}^{BFS}] .Math. .Math. .Math. \frac{1}{C_{D, \dot{}} .Math. + C_{Q, \dot{}}} = \frac{1}{C_{Q, \dot{}} (1 + \frac{C_{D, \dot{}}}{C_{Q, \dot{}}} .Math.)} .Math. .Math. .Math. .Math. k .Math. = .Math. \frac{C_{D, \dot{}}}{C_{Q, \dot{}}} .Math. < 1 .Math. .Math. \begin{matrix} Side .Math. .Math. force .Math. .Math. changes .Math. .Math. faster .Math. .Math. with .Math. .Math. a \\ variation .Math. .Math. of .Math. .Math. AOS .Math. .Math. than .Math. .Math. drag .Math. .Math. does \end{matrix} .Math. .Math. \frac{1}{C_{D, \dot{}} .Math. + C_{Q, \dot{}}} = \frac{1}{C_{Q, \dot{}} (1 + k .Math. .Math.)} \frac{1}{C_{Q, \dot{}}} .Math. (1 - k .Math. .Math.) \frac{1}{C_{Q, \dot{}}} .Math. .Math. Taylor .Math. .Math. expansion .Math. .Math. M .Math. .Math. \dot{} = - \frac{1}{C_{Q, \dot{}}} [M^{2} [(C_{D,} .Math. + C_{D,} .Math. + C_{D, M} .Math. .Math. .Math. M + C_{D, p} .Math. \frac{p}{M} + C_{D, r} .Math. \frac{r}{M}) .Math. + (C_{Q,} .Math. + C_{Q,} .Math. + C_{Q, M} .Math. .Math. .Math. M + C_{Q, p} .Math. \frac{p}{M} + C_{Q, r} .Math. \frac{r}{M})] + \frac{W_{MTOW} .Math. C_{T, M}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. S} .Math. .Math. .Math. M + \frac{m_{0}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math.} .Math. .Math. .Math. a_{2}^{BFS}] .Math. .Math. M .Math. .Math. \dot{} = - \frac{1}{C_{Q, \dot{}}} [(M^{2} .Math. C_{D,} .Math. + M^{2} .Math. C_{D,} .Math. + M^{2} .Math. C_{D, M} .Math. .Math. .Math. M + {MC}_{D, p} .Math. p + {MC}_{D, r} .Math. r) .Math. + (M^{2} .Math. C_{Q,} .Math. + M^{2} .Math. C_{Q,} .Math. + M^{2} .Math. C_{Q, M} .Math. .Math. .Math. M + {MC}_{Q, p} .Math. p + {MC}_{Q, r} .Math. r) + \frac{W_{MTOW} .Math. C_{T, M}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. S} .Math. .Math. .Math. M + \frac{m_{0}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math.} .Math. .Math. .Math. a_{2}^{BFS}]$

[0147] Taking into account [263]:

[00053] $M .Math. .Math. \dot{} = - \frac{1}{C_{Q, \dot{}}} [M^{2} .Math. C_{D,} .Math. + M^{2} .Math. C_{D,} .Math.^{2} + M^{2} .Math. C_{D, M} .Math. .Math. .Math. M .Math. .Math. + {MC}_{D, p} .Math. p .Math. .Math. + {MC}_{D, r} .Math. r .Math. .Math. + M^{2} .Math. C_{Q,} .Math. + M^{2} .Math. C_{Q,} .Math. + M_{0}^{2} (1 + 2 .Math. \frac{.Math. .Math. M}{M_{0}}) .Math. C_{Q, M} .Math. .Math. .Math. M + {MC}_{Q, p} .Math. p + {MC}_{Q, r} .Math. r + \frac{W_{MTOW} .Math. C_{T, M}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math.} .Math. .Math. .Math. a_{2}^{BFS}]$ $.Math. C_{y, M} = C_{Q, M} + \frac{W_{MTOW} .Math. C_{T, M}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math. {SM}_{0}^{2}}$ $M .Math. .Math. \dot{} = - \frac{1}{C_{Q, \dot{}}} .Math. (C_{D,} .Math. M^{2} .Math. + C_{D,} .Math. M^{2} .Math.^{2} + C_{D, M} .Math. M^{2} .Math. .Math. .Math. .Math. M + C_{D, p} .Math. Mp .Math. .Math. + C_{D, r} .Math. Mr .Math. .Math. + C_{Q,} .Math. M^{2} .Math. + C_{Q,} .Math. M^{2} .Math. + M_{0}^{2} .Math. C_{y, M} .Math. .Math. .Math. M + 2 .Math. .Math. M_{0} .Math. C_{Q, M} .Math. .Math. .Math. M^{2} + C_{Q, p} .Math. Mp + C_{Q, r} .Math. Mr + \frac{m_{0}}{\frac{1}{2} .Math. .Math. .Math. p_{0} .Math.} .Math. .Math. .Math. a_{2}^{BFS})$

[0148] As these teachings show, a manufacturer of a fixed-wings AV can inexpensively generate an accurate APM of its aircraft. This APM may be employed for flight planning, simulation or design. Advantageously, the APM may continually be fine-tune and refine, should the payload or any other aspect of the AV geometry or mass configuration change.

[0149] These and other features, functions, and advantages that have been discussed can be achieved independently in various embodiments or may be combined in yet other embodiments.

Computer-Implemented Method and System for Modelling Performance of a Fixed-Wing Aerial Vehicle with Six Degrees of Freedom

Inventors

Cpc classification

Classification Explorer

G06F30/15

PHYSICS

Classification Explorer

G06F30/20

PHYSICS

Classification Explorer

B64D43/00

PERFORMING OPERATIONS; TRANSPORTING

Classification Explorer

G01M9/08

PHYSICS

Classification Explorer

B64C17/00

PERFORMING OPERATIONS; TRANSPORTING

Classification Explorer

G05D1/0825

PHYSICS

Classification Explorer

G07C5/085

PHYSICS

International classification

Classification Explorer

G06F17/50

PHYSICS

Classification Explorer

G07C5/08

PHYSICS

Abstract

Claims

Description