Next Article in Journal
Developing a Decision Tree Algorithm for Wind Power Plants Siting and Sizing in Distribution Networks
Next Article in Special Issue
A Multilayer Brushless DC Motor for Heavy Lift Drones
Previous Article in Journal
Performance of the Condensation Process for Water Vapour in the Presence of a Non-Condensable Gas on Vertical Plates and Horizontal Tubes
Previous Article in Special Issue
Assessment of All-Electric General Aviation Aircraft
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Control and Estimation of the Condition of a Small Unmanned Aircraft Using a Kalman Filter

Department of Avionics, Faculty of Aeronautics, Technical University of Košice, 040 01 Košice, Slovakia
*
Author to whom correspondence should be addressed.
Energies 2021, 14(8), 2292; https://doi.org/10.3390/en14082292
Submission received: 16 March 2021 / Revised: 16 April 2021 / Accepted: 17 April 2021 / Published: 19 April 2021
(This article belongs to the Special Issue Algorithms and Aircraft Electric Power Systems)

Abstract

:
The article was motivated by the design of adaptive control algorithms for the control of a a fixed wing unmanned aerial vehicle (UAV). An adaptive system is a system that, with its structure or parameters, adapts to changes in the behavior of the object and based on the knowledge of variable properties, maintains the quality of its regulatory transition. The knowledge gained on this small UAV can be applied to larger aircraft. The creation of the proposed adaptive control into the UAV consisted of the creation of a simulation model of the aircraft based on known physical laws, the properties of the aircraft and a mathematical description. An adaptive PID controller for stabilization with changing coefficients based on the airspeed of the aircraft was designed and simulated. A validated control of the mathematical model of an unmanned aircraft was designed and simulated using the methods of estimation and identification of the UAV model parameter based on measured data from flight tests. Identifying dynamic parameters is a challenging task due to several factors, such as random vibration noise, interference, and sensor measurement uncertainty. The designed adaptive UAV control provides very promising results in improving the controllability of the aircraft while reducing the effect of speed changes on the stability and controllability of the system compared to the conventional PID controller. The comparison was performed on three selected types of PID controllers. The first type had fixed coefficients for the entire range of speeds calculated using the Control toolbox in MatLab. The second type also had constant coefficients over the entire range of speeds calculated using the Naslin method. The third adaptive type of PID controller had variable coefficients based on approximate polynomials dependent on the change in flight speed. The reason for the comparison was to show an increase in margin of stability using the method of variable coefficients of the PID controller based on the change of flight speeds. The obtained results show that the proposed adaptive control algorithm is robust enough to control the movement of the aircraft in the longitudinal plane and due to the introduced process and measurement errors, while the used Kalman filter effectively eliminates these errors.

Graphical Abstract

1. Introduction

At present, aircraft are equipped with various systems for flight automation, their main tasks being to alleviate the load on the aircraft crew, increase safety, comfort, and flight economy. Aircraft flight control is a closed-loop human-aircraft process or an automatic-aircraft control system. The task of suppressing unwanted oscillating movement of the aircraft is transferred to assistance systems that can work simultaneously with the pilot. These systems affect the rudder deflection, which is set by the pilot based on the measured angular velocities of the aircraft around the center of gravity. PID control methods, robust control and progressive fuzzy logic methods are currently used for aircraft control. When applying robust controllers for aircraft control systems, relatively widespread methods of H and μ —synthesis are used. A useful publication in the field of robust management is [1]. The feedback system is robust if its basic qualitative properties such as stability, quality of regulation under the influence of various faults are maintained. When there are changes in the parameters of the object, we need to use adaptive control. The theory of robust control of dynamic systems develops methods for the analysis of robust properties of real objects and proposes controllers that ensure the robustness of the feedback system. The peculiarity of this control theory of dynamical systems is that modelling, analysis of object properties, as well as synthesis of the controller is performed using an incomplete and inaccurate mathematical description of the controlled process [2]. In the past, PID controllers have been used in autopilots from simple control laws to advanced control systems, which are currently being developed into structures with variable coefficients according to flight conditions. In practice, this problem is solved by designing several controllers for linearized models at different operating points, the controllers are subsequently switched according to the current state variables, method (Gain scheduling) [3]. This process of changing the properties of the controller is called adaptation, by which we mean the process of automatic tuning of PID controllers depending on a certain parameter. This automatic tuning can be performed at the request of the operator or automatically when the parameter is changed. Although Gain scheduling is often used in aviation, it also has its disadvantages, which are manifested when the change of the adaptation parameter is significantly fast. The calculation of the PID controller coefficients itself can be a relatively lengthy process in terms of computational time. If the value of the adaptation parameter changes more than once during this process, the quality of the control may decrease [4].
During the design phase of a suitable control system for a fixed wing unmanned aerial vehicle (UAV), it was necessary to create a mathematical-simulation model that served as a test platform for the proposed control algorithms to assess the impact of these algorithms on the properties of the controlled object. The development of a mathematical-simulation model was based on physical quantities obtained during experimental flights and by measuring its aerodynamic and mechanical properties. The created mathematical-simulation model had to meet the specified conditions for accuracy and reliability of simulated quantities with a comparison of quantities obtained by measurements during experimental flights. An approach known as the “six degrees of freedom” (6DOF) was used to create a mathematical-simulation model, which is now widely used in practice. This model serves as a simulation and test platform for developed flight control and stabilization systems. Simulation testing has several advantages, which include mainly low experimental costs, repeatability, and testing critical conditions in which there is no danger. The disadvantage of simulations is their limited the number of processed quantities and inputs, which results in differences between simulated ones and real measured waveforms. Specifically in-flight simulations, the simulation algorithm does not count with environmental changes. Also entered parameters obtained by measurement, calculation or other simulations may not be highly accurate and thus, deviations from the actual state occur in flight simulations. To reduce these and more errors of the simulation model, it is necessary to create the most accurate model of the investigated object (aircraft). It is ensured by entering the most accurate coefficients into the model airplane [5]. The implementation of the mathematical-simulation model was performed in the MatLab/Simulink environment. In the next stage, a control system was created, the task of which was to stabilize the angular velocities of the aircraft. The adaptive control algorithm was designed and tested on the created simulation model longitudinal motion of airplane, verification of the operation to control and stabilize the angular velocity of angle of longitudinal inclination was realized using the simulation model of longitudinal motion [6]. After the design, testing and tuning stage, a comparison of the dynamic properties of the proposed adaptive control algorithm with a simple type of PID controller was performed. The use of a steering assistance system, where a controller is used, whose task is to control angular velocities based on the requirements of the pilot, eliminates the effort and necessary experience of the pilot. Such a control system eliminates the influence of flight parameters on the control because the rudder deflection and angular velocity are no longer directly controlled by the parameter and measured by gyroscopic sensors. Similar systems exist in practice and are described in [7].
In 2008, Leong [5] implemented a nonlinear 6-DOF model for the Yak-54 RC aircraft using the Aerospace Blockset. Advanced Airplane Analysis (AAA) software was used to obtain aerodynamic derivatives. The dynamics of the actuator was modelled by a first-order system. Engine dynamics have been simplified by using a first order time delay. The contact force was calculated based on the dynamic characteristics of the propeller, which were estimated using a computer program called JavaProp.
In 2015, Kamal [8] presented the process of compiling a detailed nonlinear 6-DOF model DOF for a small Tiger-Trainer RC aircraft. The propulsion model with propeller characteristics was investigated based on experimental results of wind tunnel tests, the dynamics of the engine and its servomotor were identified as a first-order model. The complete aircraft simulation model was compiled in the MatLab/Simulink environment, the procedure of tuning the model parameters was carried out using aircraft system identification techniques.
In his work in 2014, Hoffer [9] used system identification using recursive least squares methods with online filtering error as an alternative to analytical methods and wind tunnel testing of a small low-cost fixed-wing UAV. The author analyzed his identification method on UAV systems using MEMS sensors and GPS and their error signals.
In his 2009 work [6] Bazezit determined the stability coefficients of the UAV Trainer60 aircraft model using Tornado software, identified the moments of inertial forces using a special design of pendulum systems. Designed PID controllers for the internal and external loops of the Trainer60 model were simulated for modified modes in the Hardware-in-the-Loop simulator. To achieve accurate low altitude evaluation, an ultrasonic and laser UAV landing meter was used in an urban environment.
In his contribution [10] Al-Mashhadani presented in 2019 an optimal estimation and state estimation approach for controlling UAV parameters under the conditions of random vibration noise and uncertainty. The results showed that to minimize the failure rate of UAVs, the optimal estimation approach for nonlinear parameters should be used rather than a linear deterministic model.
In a 2015 paper Ravi Kumar [11] proved a model for controlling the inclination, tilt and lateral slip of an aircraft using an LQR controller and an LQG controller. The simulation and analysis results showed that the LQR controller gave relatively better performance. The LQG controller achieved very good results with the output of a device with a limited constant shift error, and the Kalman filter was the optimal estimator when working with Gaussian white noise.
Based on [12,13,14] it can be stated that a PID, LFC or LRQ controller can be used to increase the stabilization of the longitudinal inclination angle. The use of one of the controller types will increase the overall stability of the UAV and reduce its susceptibility to noise that may occur in the control process. The system’s response to the required change is shortened and unwanted oscillations are also eliminated. According to [8], in the process of estimating the optimal state of the system, it is also necessary to consider random vibrations caused either by the UAV engine or by external influences, among which the crosswind can be included.
For the needs of simulation of the selected aircraft, we must know its physical properties and be able to describe them mathematically. These are converted into simulation schemes in the form of equations and mathematical operations, the outputs of which correspond to the behavior of the real object. We will use the required rudder deflection as the input variable. The output variable will be the angular velocity of the longitudinal inclination of the aircraft. The creation of an aircraft simulation model is described in the literature [5,8,10]. During the flight, several different types of forces act on the aircraft, the interaction of which affects the movement and position of the aircraft in space: the lift, thrust, aerodynamic drag and weight. These forces are calculated based on the various characteristics of the aircraft. In the case of a small unmanned aircraft, the mass is constant during the flight and the change in gravitational acceleration due to the change in altitude is negligible. Equation (1) describes the dynamic system described by the UAV mathematical model without the influence of interfering signals with the moment of inertia included for the longitudinal movement:
W ϑ ˙ δ V s = ϑ ˙ s δ V s = n V s + n 22 s 2 + 2 d 0 ω 0 s + ω 0 2
where:
2 d 0 ω 0 = n 0 + n 22 + n 33
ω 0 2 = n 32 + n 22 n 33
The introduced constants of the transfer function have the form:
n 33 = μ m z ω z τ a
n 22 = Y + P c o s α m z V g s i n θ V
n 32 = μ m z α
n 0 = P P O M c o s α 0 S V 2 ρ
n V = M Z δ V J Z
M Z δ V = m Z b A S ρ V 2 2
where δ v is the rudder deflection angle, ϑ is the angle of longitudinal inclination, m z is the moment coefficients of all forces acting per plane equal to f V , α , ϑ , δ v , μ is the auxiliary coefficient, τ a is the delay time constant, Y is the lift, P is the engine thrust, α is the angle of attack, g is the gravitational acceleration, θ is the aerodynamic angle of longitudinal inclination, V is the speed of flight, P P O M is the engine control lever deflection, α 0 is the zero angle of attack, S is the wing area, ρ is the air density, J Z is the moment of inertia in the Z axis, M Z δ V is the inclination moment of aerodynamic forces, b A is the mean length of the aerodynamic chord of the wing. The kinematic relationship between the angles is:
ϑ = θ + α
The characteristic equation of the longitudinal inclination of an aircraft usually has two pairs of complexly associated roots, their real parts are usually negative. It is clear from the transfer function that in relation to the longitudinal inclination, the behavior of the aircraft is expressed by the combination of the frequency, proportional derivative, and integration term. Due to the angle of attack, when the rudder is deflected, the aircraft behaves as an oscillating member with its own frequency ω 0 and damping coefficient d 0 . Figure 1 shows angles and basic coordinate systems used in longitudinal motion of UAV.
The basic coordinate systems describing the movement of the aircraft include the body frame (index b), aerodynamic (index a) and earth (index e) coordinate system. The aerodynamic coordinate system represented by the X axis with the designation Xa is chosen as the basic coordinate system. The angle between the axis of the aerodynamic and body frame coordinate systems is defined as the angle of attack α . The angle between the axis of the aerodynamic and earth coordinate systems is defined as the aerodynamic angle of longitudinal inclination of the airplane θ . Between the body frame and earth coordinate system, we can detect the angular velocity of the longitudinal inclination of the aircraft.
The UAV management system is based on micro-electro-mechanical systems (MEMS) used for UAV guidance without pilot assistance, consisting of hardware supported by software. The UAV autopilot can navigate the UAV along a reference route or through waypoints. More powerful autopilot systems can control UAVs at every stage of the flight, such as take-off, climb, flight, descent, and landing. Figure 2 shows basic UAV flight control scheme.
The UAV autopilot is a closed loop control system consisting of two parts. The first is the observation part, which is represented by an inertial measurement unit (IMU) containing the sensors needed to determine the position of the UAV. The second part of closed loop control system is state estimate part, where the sensor data along with the GPS position information is passed to a filter to make an estimate of the current state, which is needed for later use in driving. This closed-loop system can be presented using the block diagram in Figure 3 [15,16,17].
Primary stabilization performed by an inner loop. The outer loop performs the activity associated with maneuvering and navigating the UAV. The outer loop is also responsible for capturing waypoints. To perform a fully automated UAV flight, it is necessary to use an autopilot that will control this flight. For its proper functionality, the autopilot needs to know the current position of the UAV in space. To determine this position, it is appropriate to use an inertial navigation system (INS), which also includes calculations and corrections. The basic element of every INS is an inertial measurement unit, called IMU. The IMUs have a compact size that allows mounting in a small UAV, but the MEMS-based sensors installed in the IMU do not have such high accuracy. Therefore, further processing of the measured IMU data and correction of the noise entered in the measurement is necessary. The accuracy of the inertial sensors plays a key role in autonomous navigation. The errors of current inertial sensors are approximately 0.01°/h for angular velocity sensors and 100 μg for acceleration sensors. These errors integrate in time and cause a positioning error, which is expressed by the measurement uncertainty per hour, for modern IMU this error is approximately 2 km/h. The measured data from angular velocity sensors represent the vector of angular velocity of navigated object to the inertial coordinate system indicated by the index i and measured in the individual axes of the Cartesian coordinate system of object navigation indicated by the index b:
ω i b b = ω i b b x ,   ω i b b y , ω i b b z T
Angular velocities are transformed between the inertial and body frame coordinate systems (subscript ib) with a final transformation into the body frame coordinate system (superscript b).
By mathematically adjusting and integrating angular velocities, we get information about the angular position of the object’s navigation relative to the reference system, such as, tilt, and slope. The data measured by a three-way accelerometer represent a vector of a specific force, resp. vector of acceleration measured in the axes of navigation of the object:
f b = f b x , f b y , f b z T
The ability of the algorithm to determine changes in navigation parameters with the required accuracy in a vibrating environment is a critical factor in determining the navigation parameters (angular position) while maintaining the required high accuracy. The source of errors in each IMU sensor can be modelled to create a sensor model. To be able to estimate any significant sources of error whose value cannot be measured, calculated or otherwise determined, parameters describing the behavior of the errors must be included in the sensor model. In practice, simplified error models of inertial sensors are most often used. The simplified error model of the uniaxial angular velocity sensor assumes that the measurable sensor output ω ˜ ib b (in the x, y or z axis) consists of the actual angular velocity ω ˜ ib b , deterministic bias d , bias-drift (randomly fluctuating bias is often referred to only as drift) d ω and broadband noise sensor w ω [18].
The paper deals with the creation of a model with an adaptive PID controller, which adjusts its coefficients based on a change in the speed of the UAV and thus increases the stock of stability. A Kalman filter is also created, which by its estimation eliminates interfering influences from the output signals and thus increases the overall stability of the UAV.

2. Materials and Methods

The Department of Avionics, Faculty of Aviation, Technical University in Košice has a “SkyDog” fixed wing unmanned aerial vehicle, which is mostly used for recreational and sport model flying. Figure 4 shows UAV “SkyDog”. The aircraft is remotely controlled by a pilot from the ground using a walkie-talkie. Its control and flight mechanics correspond to large crew aircraft and is therefore suitable for conducting experiments in the field of autonomous aircraft control. It is assumed that the measured results and the proposed systems may be applicable in practice on large commercial aircraft. The “SkyDog” UAV has a wingspan of 2 m and its take-off weight is approximately 8 kg. These parameters allow him to carry a wide range of measuring equipment or small experimental equipment.
Based on the knowledge of the equations of motion, it is necessary to determine the aerodynamic coefficients for the given type of UAV. The XFLR5 simulation program was used to determine the aerodynamic coefficients of the “SkyDog” UAV. To design systems for this aircraft, it is necessary to mathematically identify the system. In our case, we will consider the movement of the “SkyDog” UAV in the longitudinal plane under the term system. A 6DoF “SkyDog” UAV model will be used for identification. By system identification we mean the search for the relevant transfer function, which most accurately describes the movement of the aircraft in the longitudinal plane. For the possibility of further simulation, it is also necessary to consider action interventions on the control surfaces of the aircraft. In our case, the action intervention is converted by actuators, which convert the electrical signal into an action variable. The action variable is the angle of deflection of the control surface and the actuator is the servo that performs this movement. By obtaining the transfer functions of the actuator and the aircraft itself, it is possible to calculate the PID controller coefficients for the system. Figure 5 shows a control block diagram with a PID controller.
The coefficients of the regulators were calculated by using the Control toolbox in the MatLabprogram. It is characterized by the short rise time and relatively low overshoots. The mathematical form of the PID controller has the form:
F R s = K P + K I s + K D N 1 + N 1 s
The values of the PID controller coefficients depend on the coefficients of the system transfer function. The transfer function of the actuator is unchanged unless we consider a system failure or change in its properties. The transfer function of the aircraft changes with respect to the current flight speed and the current altitude. In our case, the effect of height is negligible, as we consider a small UAV, whose access is in the tens of meters. However, the effect of speed is considerable because the range of speeds at which the aircraft can perform its flight is in the range of approximately 50 to 150 km/h. If we want to maintain the quality of control in the entire operating range of aircraft speeds and thus increase the overall stability of the system during the entire flight, it is necessary to use one of the methods of adapting the controller to a variable parameter, which in this case it is airspeed [19,20].
To design an optimal “SkyDog” UAV motion control process, model blocks such as the PID controller, actuator transfer function and airplane transfer function need to be integrated into a single transfer function block. Model feedback must also be included in this block. The feedback transfer function has following form:
F F = H s 1 + H s
The Kalman filter is used to estimate the state variables of the system for a given measurement of the output signals of the system containing the measured noise. Figure 6 shows a block diagram of the Kalman filter connection.
The filter estimation is optimal in the sense that the sum of the errors of the estimated state variables at any time has the smallest possible value. We consider the difference between the estimate processed by the filter and the actual value of the system state variable to be an estimation error. Thus, the Kalman filter is used for the best determination of variables, in other words the optimal suppression of the measured noise. When estimating the state vector, we start from the equations characterizing the discrete state space. The considered n-dimensional Gaussian white noise w k has a mean value equal to zero. It is necessary to perform an IMU analysis installed in the UAV, as its sensors introduce errors into the measurement process. These errors can cause a decrease in the stability and manageability of the UAV. In terms of sensor error analysis, Allan variance method is one of the basic tools for analysis. Its use is determined by the IEEE standard [21], which defines not only the application of this tool, but also the procedure of analysis of results. The result is the representation of noise in the output of the IMU inertial sensors in terms of long-term development introduced into the UAV control.
Figure 7 shows the types of noise that can be identified from the Allan variation graph.
Process noise is created by disturbances that act on the UAV during flight and landing, these are moments caused by horizontal and vertical air flow. The action of the wind causes a change in the modulus of the airspeed vector, the angle of attack, the tilt angle, and the tilt speed. Another component of the process failure is the pilot himself in the control of the UAV as a member of the control complex. The pilot’s sensory organs perform the function of sensors, the pilot uses the central nervous system to process information. It uses executive bodies for the implementation of operations, the interconnection of the above elements acts in unity as a regulatory member. The simplified theory of the pilot-operator model, designed to control the longitudinal motion of the aircraft, is a rough empirical value of the real values needed to mathematically express the task of controlling the angle of inclination θ and angle of attack α at different speeds of flight UAV. At present, experimental methods for determining the pilot’s ability to control a given type of UAV are well developed, using special simulators, together with mathematical models of the properties of a “piloted UAV”. For the spectral density of random perturbations of air flow functions, for the longitudinal axis we start from the relation:
S x ω = 2 σ x 2 L x π 1 1 + L x V ω 2
where σ x 2 is turbulence intensity value, ω is angular velocity, V . is air speed, Lx is a measure of the lengths of turbulence acting on the individual components of the velocity vector which has the form:
L x = H 0.177 + 0.000823 H 1 2
where H is the altitude. The intensity of turbulence acting on the UAV can be calculated according to the relation:
σ x = 0.1 * W 20
where W 20 is wind speed at the altitude of 6 m. For slight turbulence, this value is given by wind speed 7.72   ms 1 , for medium turbulence 15.43   ms 1 and for strong turbulence 23.15   ms 1 .
Based on these process noise components, the process error vector w k was determined w k is the vector of process errors (Gaussian white noise with zero mean value):
w k = 0.07 0 0 0 0 0.07 0 0 0 0 0.07 0 0 0 0 0.07
Measurement noise is generated by the UAV sensors themselves and their error models. For the analysis of the sensors, a static measurement of the IMU UAV was performed, focused on the measurement errors of the IMU PixHawk PX-4 gyroscope. Static measurements were performed in a time window of 3 h with a sampling frequency of 100 Hz. The IMU was stored in a horizontal position isolated from the surroundings, thus ensuring the stability of the measurement and the elimination of interfering vibrations. As little as the angular velocity of the Earth’s rotation is negligible, we can assume that the data measured by the gyroscope in the Y-axis show the magnitude of the individual errors affecting the IMU measurements. For recognizing the measuring noise, the measurement results were analyzed by Allan variation. The coefficients for the noise mechanisms angle random walk and bias drift were determined from the Allan variation graph. The bias drift value for the Y-axis gyroscope is B Y = 8.0712 * 10 5 ° / s . The standard deviation of the wideband noise representing the angle random walk is σ A R W 2 = 4.8147 * 10 9 ° / s . It should be noted that although the static measurement was performed for 3 h, the “SkyDog” UAV can fly approximately 40 min with a maximum flight altitude of approximately 150 m. Based on these components of the measuring noise, the measurement error vector v k was determined. The measurement error vector is (Gaussian white noise with zero mean value):
v k = 0.0054
The connected Kalman filter estimates the non-measurable output of the longitudinal slope of the UAV and filters the noisy signals from the other outputs of airspeed, angle of attack and roll rate. In the design of Kalman filter, the covariance of the Q and R matrix proved to be optimal. In determining the covariance matrices of Q and R, it was necessary to find a compromise between the accuracy of the estimate and the smoothing of noisy signals. The diagonal elements of the matrix represent the scatter of the individual sensors and in the simulations, they had Q = d i a g 0.001   0.001   0 . 001   0 . 001 . Optimization of Q and R matrices is essential for the best possible KF response to noisy signals, wind gusts, pilot characteristics or an inaccurate aircraft model.
To use the UAV model in the Kalman filter, it is necessary to convert the created state space of the UAV model into its discrete form. Discrete state space is characterized by the following relationships:
x k + 1 = A D x k + B D u k + G w k
y k + 1 = C D x k + 1 + D D v k + 1
where x k + 1 is an n-dimensional state vector, y k + 1 is an m-dimensional measurement vector, v k + 1 is an m-dimensional measurement error vector (Gaussian white noise), u k is the input signal, w k is an n-dimensional Gaussian white noise with zero mean, A D is a system matrix of n × n , B D is an excitation matrix of n × m , C D is an output matrix of m × n , D D is a direct feedback matrix bonds if k = 1 , 2 , 3 , is the sample number. To calculate the matrices of the discrete state space, it is possible to fully use the matrices of the continuous state space A , B , C and D . The value of the matrix A D is determined from the matrix A for continuous state space model using the following relation:
A D = A D T = e A T
We determine the matrix B D from the matrix B for continuous state space model using the relation:
B D = A 1 e A T I B
where T is the discretization time and I is the unit matrix. The matrices C D and D D remain unchanged as in the continuous state space, as the output equation does not change with the initial conditions of the state vector when converting from the continuous shape to the discrete shape.

3. Results

The values of forces and moments acting on the aircraft also affect the geometric coefficients. These coefficients are affected by the dimensions of the aircraft and were obtained from the technical documentation and the measurement of the aircraft itself. Table 1 shows the identified geometric coefficients that were used in the next step of creating the “SkyDog” UAV simulation model.
Using the XFLR5 program, aerodynamic coefficients for UAV “SkyDog” were obtained. Their values are given in Table 2.
Based on the obtained aerodynamic coefficients shown in Table 2, a mathematical simulation model of the dynamic movement of the “SkyDog” UAV was created as shown in the Figure 8. Simulations were performed at a sampling time of T = 0.01   s .
To verify it, a comparison was made in the longitudinal plane. During the comparison, the angle of longitudinal inclination of the aircraft θ was observed in a few seconds because the deviation of the simulated data from the measured ones increases with time. This phenomenon is since the model does not include forces acting on the aircraft during real flight. This force group can include wind forces and the like. The simulation model is not expected to deviate zero from a real aircraft. The essence of creating the model is to interpret as accurately as possible realistically the outputs measured during the flight of the aircraft in a simulation environment under the same initial conditions. Figure 9 shows a comparison of the angle of longitudinal inclination during real and simulation flight.
As can be seen, the simulated output course of the pitch angle is almost identical to the course during the experimental flight. In the monitored section, the absolute difference between these two values is 0.046 series. The accuracy of position angle simulation is important for the further design of flight stabilization and control systems, as position angle output information enters control systems. Therefore, great emphasis is placed on the accuracy of the simulation in position angles of the aircraft. Using the 6-DOF model UAV, the waveforms of the input signal in the form of rudder deflection and the output signals in the form of the angular velocity of the longitudinal inclination ϑ ˙ were obtained. Using the Gauss-Neton iteration method, a transfer function describing the motion of the UAV in the longitudinal plane was obtained. Its mathematical form is as follows:
F ϑ ˙ δ v = b 1 s + b 0   a 2 s 2 + a 1 s + a 0  
Based on the knowledge of the transfer function, mathematical forms of approximate polynomials for individual members of the transfer function were created, the shape of which is:
b 1 = 0.0954 v 2 13.8807 v + 538.22
b 0 = 0.0221 v 2 + 3.8803 v 141.61
a 2 = 1
a 1 = 0.0034 v 2 0.5181 v + 23.85
a 0 = 0.0068 v 2 + 0.1018 v 2.2698
The transfer function was calculated for a constant flight speed of 90 km/h. This flight speed was chosen because it is the “SkyDog” UAV cruising speed and so a substantial part of the experimental flights take place at approximately this speed. The same method was used to determine the transfer functions at the other two speeds, namely 60 km/h and 120 km/h. The transfer functions for the given flight speeds have the following form:
F ϑ ˙ 60 δ v = 48.82 s + 11.57 s 2 + 4.881 s + 1.401
F ϑ ˙ 90 δ v = 61.7 s + 28.43 s 2 + 4.482 s + 1.41
F ϑ ˙ 120 δ v = 246.3 s + 5.47 s 2 + 10.14 s + 0.2014
As mentioned in the previous chapter, it is also necessary to know the transfer function of the actuator that controls the control surfaces of the aircraft. The measurement therefore obtained the transfer function of the actuator for its mathematical model, which has the form:
F δ v δ a c t = 1 0.125 s + 1
The input to the actuator block is the required rudder deflection and the output is the real rudder deflection. The transfer function of the actuator is already included in the transfer functions of the airplane, insofar as these transfer functions were calculated from the data measured during the flight. After obtaining mathematical models of the aircraft and actuator, it was possible to design a PID controller. Since the transfer function of the aircraft also changes with the change of air speed, it is necessary to change the PID controller coefficients so that its properties are preserved together with the stability and controllability of the object. For this reason, the approximation polynomials of the PID controller coefficients were mathematically expressed by tuning the parameters for the variable flight speed in the MatLab program. Equations (34)–(37) describe the mathematical form of the approximate polynomial for the calculation of the individual PID controller coefficients:
K P = 9.756 * 10 5 v 3 0.0294 v 2 + 2.6253 v 57.5736
K I = 0.00059 v 3 0.1763 v 2 + 15.4961 v 326.5208
K D = 4.0255 * 10 6 v 3 0.0012 v 2 + 0.1106 v 0.317593
N = 0.002666
Based on the knowledge of the approximate polynomials, it is possible to calculate the PID controller coefficients for all speeds in the “SkyDog” UAV operating speed range. Table 3 shows the coefficients for the individual speeds.
After the identification of mathematical models for individual blocks simulating the longitudinal movement of the “SkyDog” UAV, it is possible to combine them into a complex transfer function, which is needed for further processing. We consider the transfer function to be a complex transfer function, which includes the transfer functions of the designed PID controller, the servo mechanism, and the “SkyDog” UAV transfer function itself. Complex transfer functions for simulations converted at flight speeds of 60, 90 and 120 km/h have the following form:
F c ϑ ˙ 60 δ v = 9640 s 3 + 2.493 * 10 5 s 2 + 1.642 * 10 6 s + 3.752 * 10 5 s 4 + 1.002 * 10 4 s 3 + 2.512 * 10 5 s 2 + 1.642 * 10 6 s + 3.752 * 10 5
F c ϑ ˙ 90 δ v = 9654 s 3 + 2.457 * 10 5 s 2 + 1.618 * 10 6 s + 6.944 * 10 5 s 4 + 1.003 * 10 4 s 3 + 2.474 * 10 5 s 2 + 1.619 * 10 6 s + 6.944 * 10 5
F c ϑ ˙ 120 δ v = 9548 s 3 + 2.994 * 10 5 s 2 + 2.261 * 10 6 s + 5.006 * 10 4 s 4 + 9933 s 3 + 3.032 * 10 5 s 2 + 2.261 * 10 6 s + 5.006 * 10 4
Bode diagrams can be used to evaluate the adaptive controller. Logarithmic amplitude-frequency characteristics are an expression of the properties of transfer functions and angular velocity. Analysis of the tilting properties of the aircraft shows that the parameters of its longitudinal movement, i.e., its natural frequency and the coefficient of relative damping, do not differ from the required values. The optimal values of the natural frequency and the coefficients of ratio damping d 0 and ω 0 from the pilot’s point of view were determined correctly.
The Bode diagrams presented below were created based on the general “SkyDog” UAV transfer function given in Equation (24), which is bound to the approximation polynomials of the individual members of the transfer function given in Equations (25)–(29) for variable flight speed. The transfer function thus created was multiplied by the transfer function of the created PID controller and analyzed by Bode diagrams. Two types of robust controllers were used for comparison. The first robust controller used coefficients calculated using the Control Toolbox in MatLab for a speed of 90 km/h, while the second type of robust controller uses coefficients calculated by Naslin’s method. These coefficients were subsequently applied to the entire range of “SkyDog” UAV flight speeds. The third type of controller comparison is an adaptive controller, which calculates the values of the coefficients based on the change in UAV flight speed, which are shown in Table 3. In each experiment, the margin of stability of the overall system in the transmission and phase characteristics was monitored. The first simulation experiment followed the properties of a system that used a classical PID controller with constant coefficients for a range of speeds. The coefficients have the following values: K P = 10.25 ; K I = 65.12 ; K D = 0.3898 and N = 0.002666 . The graphical result of the experiment is shown in Figure 10.
The first condition for maintaining the stability of the system is that at the moment when the transmission characteristic reaches the value 0 dB, the phase shift does not exceed the value −180°. The second condition says that, when the phase characteristic exceeds the value −180°, the transmission system must be less than 0 dB. The difference between these values is called the margin of stability. The larger the margin of stability, the more stable the system. From the simulation graph the first conditions are met to what extent the phase shift of the system is less than −180°. Since the system never reaches a phase shift of −180°, the system is always stable from the point of view of the second condition. Failure to exceed the phase shift value of −180° means that the supply of stability in the transmission characteristic is infinite. In the second experiment, a system using a PID controller with coefficients determined by Naslin’s method was monitored.
The coefficients were constant for the whole range of velocities and had the following values: K P = 3.4741 ; K I = 13.6215 ; K D = 0.3364 and N = 0.002666 . The graphical result of the experiment is shown in Figure 11.
This type of PID controller meets the stability conditions as well as the previous type of controller. The third type of PID controller investigated is an adaptive controller using the Gain scheduling method. The controller coefficients are calculated based on the approximate polynomials dependent on the UAV airspeed. The value of the coefficients for each speed is shown in Table 3. The graphical result is shown in Figure 12.
As can be seen from the graphs, the use of such a controller ensures the same result at almost the entire range of speeds. The system again meets both conditions and therefore the supply of stability in the transmission part is infinite. At 150 km/h the system has different transmission and frequency characteristics from other simulated speeds. This phenomenon may be since the UAV can fly at a speed of 150 km/h, but its stability is beginning to decline. From the created complex transfer function for the longitudinal movement of the “SkyDog” UAV it is possible to create a model in the state space. By performing the conversion to the state space, we get three models in the state space, which the individual matrices have the following form:
Model at 60 km/h:
A 60 = 1.002 * 10 4 1 0 0 2.512 * 10 5 0 1 0 1.642 * 10 6 0 0 1 3.752 * 10 5 0 0 0
B 60 = 9640 2.493 * 10 5 1.642 * 10 6 3.752 * 10 5
Model at 90 km/h:
A 90 = 1.003 * 10 4 1 0 0 2.474 * 10 5 0 1 0 1.619 * 10 6 0 0 1 6.944 * 10 5 0 0 0
B 90 = 9654 2.457 * 10 5 1.618 * 10 6 6.944 * 10 5
Model at 120 km/h:
A 120 = 9933 1 0 0 3.032 * 10 5 0 1 0 2.261 * 10 6 0 0 1 5.006 * 10 4 0 0 0
B 120 = 9548 2.994 * 10 5 2.261 * 10 6 5.006 * 10 4
The matrices C and D are the same for all speeds and have the form:
C = 1 0 0 0
D = 0
Simulations were performed with the UAV model for three selected speeds. Process noise was estimated based on the assumed pilot control errors, the effect of atmospheric errors up to a height of 150 m. This process noise is then transformed into velocity, angle of attack, angle of longitudinal inclination and angular velocity of UAV in longitudinal inclination.
To the extent that the Kalman filter works in discrete time, it is necessary to convert the created models in state space into a discrete form. Based on Equations (22) and (23), we create models in a discrete state space with a discretization time T = 0.01 s. Matrices C and D remain the same even in the discrete form of state space:
C 60 = C 90 = C 120 = C D 60 = C D 90 = C D 120
D 60 = D 90 = D 120 = D D 60 = D D 90 = D D 120
Matrices A and B change their shape in a discrete form and have the following shape. Discrete time model at 60 km/h:
A D 60 = 0.002086 7.734 * 10 5 8.746 * 10 7 4.515 * 10 9 20.86 0.7729 0.008841 4.612 * 10 5 127.4 1.438 0.9926 0.009975 29.02 0.3282 0.001694 1
B D 60 = 0.9711 90.92 682.1 156.8
Discrete time model at 90 km/h:
A D 90 = 0.002058 7.757 * 10 5 8.752 * 10 7 4.515 * 10 9 20.61 0.7762 0.008858 4.618 * 10 5 126.2 1.42 0.9927 0.009975 53.87 0.6078 0.003135 1
B D 90 = 0.9711 90.63 672.2 289.8
Discrete time model at 120 km/h:
A D 120 = 0.00245 7.368 * 10 5 8.592 * 10 7 4.476 * 10 9 24.28 0.7294 0.008609 4.532 * 10 5 166.6 1.943 0.9899 0.009966 3.688 0.04302 0.0002241 1
B D 120 = 0.9708 95.78 952.8 21.11
Since the proposed model of the longitudinal inclination introduces process noise and measurement noise into the process, it is necessary to define these two parameters, for which a static measurement of the IMH PixHawk PX-4 was performed. The IMU is installed on board the “SkyDog” UAV and has the task of measuring and storing flight data. During the static measurement, the IMU was statically stored on a flat surface, which means that the output data of the accelerometers and gyroscopes correspond to their noise, which they introduce into the measurement. In the case of the longitudinal movement of the UAV, the simulations require data from a gyroscope, which has its axis of sensitivity identical to the y-axis of the aircraft coordinate system. The Allan variation method was used for sensor analysis. Graphs of the Allan variation for the acceleration and angular velocity sensors of all three aircraft axes are shown in Figure 13. The initial analysis was performed for all three axes of the IMU PixHawk PX-4 accelerometers and gyroscopes. To what extent we simulate the longitudinal movement of the aircraft, the Y-axis gyroscope data is authoritative for us. The basic parameter by which it is possible to determine the accuracy of the sensor is the standard deviation of broadband noise, known as angular random walk (ARW). In the Allan variation graph, this noise mechanism is defined by the standard by a curve with a slope of 1 2 . For our chosen gyroscope, the standard deviation value of the wideband noise is 7.3548 × 10−4. The second parameter determining the accuracy of the sensor is the bias instability, which is defined by a curve with a slope of 0. The size of the bias instability in the gyroscope’s Y axis is 8.0878 × 10−5.
We defined the standard deviation of the measured noise signal of the IMU gyroscope as process noise. The value of measurement uncertainty type A was chosen for the measurement noise. The resulting noise values are:
σ G y 2 = 0.00070910
u A G y = 0.0054
The UAV model has thus been created, which is shown in Figure 14.
The output of the models created in this way is very noisy compared to the ideal model. The output of the model with the noises is shown in Figure 15. The figure shows the change in the angular velocity of the longitudinal inclination over time.
Two types of Kalman filters have been developed to effectively filter such a signal as well as to estimate the state. For its correct functionality, it is necessary to determine the covariance of the matrix Q and R , which significantly affect the calculation of the gain of the Kalman filter. The values of covariance matrices were determined by estimation and their size is:
Q = 0.001 0 0 0 0 0.001 0 0 0 0 0.001 0 0 0 0 0.001
R = 0.5
For the Q and R was determined, a gain was calculated for each of the three simulated velocities, the value of which is:
K 60 = 4.1427 * 10 6 0.0414 0.2533 0.0577
K 90 = 4.0862 * 10 6 0.0409 0.2509 0.1071
K 120 = 4.8670 * 10 6 0.0482 0.3315 0.0073
Created models of Kalman filters are shown in Figure 16.
Figure 17 shows the real state of the system and an estimate of the state of the system for the calculated gain and their statistical data at a UAV speed of 90 km/h.
Figure 18 shows that the proposed adaptive controller loaded with process noise and measurement errors using a Kalman filter optimally monitors the angular velocity UAV with minimal errors. The conclusions resulting from the plots shown in Figure 17 confirm the values of covariance errors. The figure shows the change in the angular velocity of the longitudinal inclination over time.
The values of the covariance errors of the unfiltered signal and the signal after filtering at selected values of the Q and R matrices are given in Table 4.

4. Discussion

The stability and controllability of UAVs in critical flight situations can be reduced by induced illicit, but also accidental stimuli or conditions. The mentioned properties of the control complex can be the cause of dangerous conditions that require risky decisions of the operator—pilot or the intervention of assistance systems. This paper describes the physical properties of the UAV control object on a mathematical model in the form of nonlinear equations, these equations are linearized and then implemented in the software environment. Applied methods of aircraft control in longitudinal motion were analyzed on classical approaches and on new progressive methods. The PID parameters of a robust and adaptive controller for longitudinal control of a fixed-wing unmanned aircraft were determined by the method of analysis and synthesis. The proposed optimal discrete Kalman filter allows to determine the optimal discrete output signal (angular velocity information), while very effectively eliminating modulated process and measurement noise for different UAV speeds (Figure 15 and Figure 16). The function of the longitudinal velocity angular velocity estimation algorithm has been successfully applied not only to the model UAV but also its involvement with adaptive flight control systems. The requirements for the created control system were determined by the criteria contained in the identification of the control object so that the control system was robust enough to ensure the stability of the system at different speeds of the UAV. At the same time, the proposed control system is characterized by a certain degree of adaptation, the advantages of adaptation were visible compared to the use of the classic robust controller in Figure 11. To meet the main goal, it was necessary to sunk one of the sub-goals. The created simulation model works based on known mathematical relations describing the motion of an aircraft in three-dimensional space with six degrees of freedom. The construction of the model was based on the technical knowledge of the selected UAV and simulation experiments. This model is characterized by relative simplicity with relatively high accuracy of simulation. An adaptive controller was designed together with an aircraft model for a selected speed range which was burdened by process and measurement errors. The designed controller using Kalman filter optimally monitored the angular velocity of the UAV in the longitudinal plane with minimal errors, which is documented in Figure 16. The conclusions from Table 4 confirm the allowable values of covariation errors of our proposed adaptive UAV control system.

Author Contributions

Conceptualization, D.M., R.B. and M.S.; methodology, R.B.; validation, D.M., R.B. and M.S.; formal analysis, D.M., R.B. and M.S.; investigation, D.M. and R.B.; data curation, D.M. and M.S.; writing—original draft preparation, D.M. and M.S.; writing—review and editing, D.M. and R.B.; visualization, D.M.; supervision, R.B.; funding acquisition, R.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This work was supported by project no. KEGA 044TUKE-4/2019—a small, unmanned airplane—the platform for education in the area of intelligent avionics.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Szászi, I.; Kulcsár, B. Robust Control and Fault Detection Filter Design for Aircraft Pitch Axis. Period. Polytech. Transp. Eng. 2001, 29, 83–100. [Google Scholar]
  2. Liu, M.; Egan, G.K.; Santoso, F. Modeling, Autopilot Design, and Field Tuningof a UAV With Minimum Control Surfaces. IEEE Trans. Control Syst. Technol. 2015, 23, 2353–2360. [Google Scholar] [CrossRef]
  3. Raju, S.S.; Srividya, A.; Sudarsan, S.; Kiran, R. Design and Implementation of Model Predictive Control for Non-Linear System and Performance Comparison with Conventional PID. In Proceedings of the IX Control Instrumentation System Conference, Manipal, Karnataka, India, 16–17 November 2012. [Google Scholar]
  4. Aishah, J.; Fitri, Y.; Hatta, A.; Zainudin, A.R.; Shamsul, S.; Rudzidatulm, A.D.; Ab Rashid, M.Z.; Azizan, A.; Mori, Y. Improvement of Pitch Motion Control of an Aircraft Systems. TELKOMNIKA 2018, 16, 2263–2274. [Google Scholar]
  5. Leong, H.I. Development of a 6DoF Nonlinear Simulation Model Enhanced with Fine Tuning Procedures. Ph.D. Thesis, University of Kansas, Lawrence, KS, USA, 2008. [Google Scholar]
  6. Bayezit, I. The Autonomy of Fixed-Wing Aerial Vehicles and Experimental Design Steps on Implementing Autonomous Navigation, Landing And Take Off Operations For The Trainer60 Model Aircraft. In Proceedings of the 5th Ankara International Aerospace Conference, METU, Ankara, Turkey, 17–19 August 2009. [Google Scholar]
  7. Vural, M.; Nicolai, L.M. Estimating R/C Model Aerodynamics and Performance; Illinois Institute of Technology: Chicago, IL, USA, 2009. [Google Scholar]
  8. Kamal, A.M.; Bayoumy, A.A.; Elshabka, A. Modeling, Analysis and Validation of a Small Airplane Flight Dynamics. In Proceedings of the AIAA Modeling and Simulation Technologies Conference, Kissimmee, FL, USA, 5–9 January 2015. [Google Scholar]
  9. Hoffer, N.V. System Identification of a Small Low-Cost Unmanned Aerial Vehicle Using Flight Data from Low-Cost Sensors. Ph.D. Thesis, Utah State University, Logan, UT, USA, 2014. [Google Scholar]
  10. Al-Mashhadani, M.A. Optimal control and state estimation for unmanned aerial vehicle under random vibration and uncertainty. Meas. Control. 2019, 52, 1264–1271. [Google Scholar] [CrossRef] [Green Version]
  11. Ravi Kumar, M.G.; Ramamurthy, T.V. Fault Detection of Aircraft Plant Using KALMAN FILTER. IJCTA 2015, 8, 923–932. [Google Scholar]
  12. Ingabire, A.; Sklyarov, A.A. Control of longitudinal flight dynamics of a fixedwing UAV using LQR, LQG and nonlinear control. E3S Web Conf. 2019, 104, 02001. [Google Scholar] [CrossRef]
  13. Hajiyev, C.; Vural, S.Y. LQR Controller with Kalman Estimator Applied to UAV Longitudinal Dynamics. Positioning 2013, 4, 36–41. [Google Scholar] [CrossRef] [Green Version]
  14. Chrif, L.; Kadda, Z.M. Aircraft Control System Using LQG and LQR Controller with Optimal Estimation-Kalman Filter Design. Procedia Eng. 2014, 80, 245–257. [Google Scholar] [CrossRef] [Green Version]
  15. Ahmed, W.; Li, Z.; Maqsood, H.; Anwar, B. System Modeling and Controller Design for Lateral and Longitudinal Motion of F-16. Autom. Control Intell. Syst. 2019, 4, 39–45. [Google Scholar] [CrossRef] [Green Version]
  16. Chaoqiang, L. An Estimation Method for Aircraft Longitudinal Static Stability Prior to Flight Test. In Proceedings of the 2nd International Conference on Computer Application and System Modeling, Taiyuan, Shanxi, China, 27–29 July 2012. [Google Scholar]
  17. Chao, H.Y.; Cao, Y.C.; Chen, Y.Q. Autopilots for Small Unmanned Aerial Vehicles: A Survey. Int. J. Control Autom. Syst. 2010, 8, 36–44. [Google Scholar] [CrossRef]
  18. Titterton, D.; Weston, J. Strapdown Inertial Technology; Institution of Engineering and Technology: Stevenage, UK, 2004. [Google Scholar]
  19. Noshahri, H.; Kharrati, H. PID controller design for unmanned aerial vehicle using genetic algorithm. In Proceedings of the 23rd International Symposium on Industrial Electronics, Istanbul, Turkey, 1–4 June 2014. [Google Scholar] [CrossRef]
  20. Mehdi, H. Unmanned Aerial Systems, Chapter 22—UAV geofencing: Navigation of UVAs in constrained environments. In Theoretical Foundation and ApplicationsAdvances in Nonlinear Dynamics and Chaos (ANDC); Academic Press: Cambridge, MA, USA, 2021; pp. 567–594. [Google Scholar]
  21. Freescale Semiconductor, Inc. Allan Variance: Noise Analysis for Gyroscopes 2015. Available online: https://telesens.co/wp-content/uploads/2017/05/AllanVariance5087-1.pdf (accessed on 10 February 2015).
Figure 1. Longitudinal movement of the aircraft.
Figure 1. Longitudinal movement of the aircraft.
Energies 14 02292 g001
Figure 2. UAV flight control system.
Figure 2. UAV flight control system.
Energies 14 02292 g002
Figure 3. Block diagram of the UAV control loop.
Figure 3. Block diagram of the UAV control loop.
Energies 14 02292 g003
Figure 4. “SkyDog” fixed wing UAV.
Figure 4. “SkyDog” fixed wing UAV.
Energies 14 02292 g004
Figure 5. Block diagram of PID controller.
Figure 5. Block diagram of PID controller.
Energies 14 02292 g005
Figure 6. Block diagram of Kalman filter.
Figure 6. Block diagram of Kalman filter.
Energies 14 02292 g006
Figure 7. Identification of noises from the Allan variation graph.
Figure 7. Identification of noises from the Allan variation graph.
Energies 14 02292 g007
Figure 8. Mathematical simulation 6DoF model of the dynamic movement of the “SkyDog” UAV.
Figure 8. Mathematical simulation 6DoF model of the dynamic movement of the “SkyDog” UAV.
Energies 14 02292 g008
Figure 9. Comparison of real and simulation flight in the longitudinal plane.
Figure 9. Comparison of real and simulation flight in the longitudinal plane.
Energies 14 02292 g009
Figure 10. Bode diagrams for a classic PID controller.
Figure 10. Bode diagrams for a classic PID controller.
Energies 14 02292 g010
Figure 11. Bode diagrams for PID controller using Naslin’s method.
Figure 11. Bode diagrams for PID controller using Naslin’s method.
Energies 14 02292 g011
Figure 12. Bode diagrams for adaptive PID controller.
Figure 12. Bode diagrams for adaptive PID controller.
Energies 14 02292 g012
Figure 13. Graphs of Allan variation of IMH PixHawk PX-4 sensors.
Figure 13. Graphs of Allan variation of IMH PixHawk PX-4 sensors.
Energies 14 02292 g013
Figure 14. Created model of longitudinal movement of UAV including noise.
Figure 14. Created model of longitudinal movement of UAV including noise.
Energies 14 02292 g014
Figure 15. Graphical outputs of created UAV models including noises for a speed of 90 km/h.
Figure 15. Graphical outputs of created UAV models including noises for a speed of 90 km/h.
Energies 14 02292 g015
Figure 16. Created model of Kalman filters for UAV models.
Figure 16. Created model of Kalman filters for UAV models.
Energies 14 02292 g016
Figure 17. The course of the system state and its estimation for the calculated gain.
Figure 17. The course of the system state and its estimation for the calculated gain.
Energies 14 02292 g017
Figure 18. Graphical course of angular velocity for UAV speed 90 km/h.
Figure 18. Graphical course of angular velocity for UAV speed 90 km/h.
Energies 14 02292 g018
Table 1. “SkyDog” geometric coefficients.
Table 1. “SkyDog” geometric coefficients.
AbbreviationName of VariableValue
W s u r f Wing surface0.68 m2
W s p a n Wing span2.04 m
W c Middle wing chord0.32 m
C G X Centre of gravity position on axis X0.08 m
C G Y Centre of gravity position on axis Y0 m
C G Z Centre of gravity position on axis Z0.02 m
C P X Pressure centre position on axis X0.12 m
C P Y Pressure centre position on axis Y0 m
C P Z Pressure centre position on axis Z0 m
Table 2. “SkyDog” aerodynamic coefficients.
Table 2. “SkyDog” aerodynamic coefficients.
EffectsShort CutDependency on VariableValue
Aerodynamic’s Coefficients of ForcesAxis Z C L _ α Dependency of lift on α 4.56
C L _ α 0 Dependency of lift on α = 0 0.51
C L _ Δ EL Dependency of lift on Δ E L −0.85
CL_qDependency of lift on q7.2
Axis X C D _ α Dependency of drag on α 0.39
C D _ α 0 Dependency of drag on α = 0 0.03
Axis Y C Y _ β Dependency of side force on β −0.18
C Y _ r Dependency of side force on r 0.314
C Y _ Δ R U D Dependency of side force on Δ R U D 0.2
Aerodynamic’s Coefficients of MomentsAxis X C l _ Δ R U D Dependency of roll moment on Δ R U D 0.004
C l _ Δ A L Dependency of roll moment on Δ A L   −0.1
C l _ β Dependency of roll moment on β 0.012
C l _ p Dependency of roll moment on p −0.32
C l _ r Dependency of roll moment on r 0.028
Axis Y C m _ Δ E L Dependency of pitch moment on Δ E L −2.442
C m _ q Dependency of pitch moment on q −48
C m _ α Dependency of pitch moment on α −1.13
C m _ α 0 Dependency of pitch moment on α = 0 −0.04
Axis Z C n _ r Dependency of yaw moment on r −0.09
C n _ p Dependency of yaw moment on p −0.07
C n _ Δ A L Dependency of yaw moment on Δ A L −0.12
C n _ Δ R U D Dependency of yaw moment on Δ R U D −0.04
C n _ β Dependency of yaw moment on β 0.089
Table 3. PID controller coefficients for individual airspeeds.
Table 3. PID controller coefficients for individual airspeeds.
Air Speed [km/h] K P K I K D N
508.51959.530.2920.002666
6013.2686.460.49110.002666
7018.17112.50.70850.002666
801698.10.62980.002666
9010.2565.120.38980.002666
1006.40843.070.22940.002666
1104.32731.060.14280.002666
1203.17324.40.094890.002666
1302.498220.067020.002666
1402.05220.310.049110.002666
1502.77631.730.053320.01917
Table 4. Error covariance of Kalman filter at selected values of Q and R matrices.
Table 4. Error covariance of Kalman filter at selected values of Q and R matrices.
Matrix ValuesBefore FiltrationAfter Filtration
Q = 0.001
R = 0.5
0.0860350.030103
Q = 0.001
R = 0.05
0.0860350.030125
Q = 0.01
R = 0.05
0.0860350.032823
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Megyesi, D.; Bréda, R.; Schrötter, M. Adaptive Control and Estimation of the Condition of a Small Unmanned Aircraft Using a Kalman Filter. Energies 2021, 14, 2292. https://doi.org/10.3390/en14082292

AMA Style

Megyesi D, Bréda R, Schrötter M. Adaptive Control and Estimation of the Condition of a Small Unmanned Aircraft Using a Kalman Filter. Energies. 2021; 14(8):2292. https://doi.org/10.3390/en14082292

Chicago/Turabian Style

Megyesi, Dávid, Róbert Bréda, and Martin Schrötter. 2021. "Adaptive Control and Estimation of the Condition of a Small Unmanned Aircraft Using a Kalman Filter" Energies 14, no. 8: 2292. https://doi.org/10.3390/en14082292

APA Style

Megyesi, D., Bréda, R., & Schrötter, M. (2021). Adaptive Control and Estimation of the Condition of a Small Unmanned Aircraft Using a Kalman Filter. Energies, 14(8), 2292. https://doi.org/10.3390/en14082292

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop