Next Article in Journal
Potentiometric Aptasensing of Vibrio alginolyticus Based on DNA Nanostructure-Modified Magnetic Beads
Next Article in Special Issue
Cramer-Rao Lower Bound Evaluation for Linear Frequency Modulation Based Active Radar Networks Operating in a Rice Fading Environment
Previous Article in Journal
A Touch Sensing Technique Using the Effects of Extremely Low Frequency Fields on the Human Body
Previous Article in Special Issue
Development of Cloud-Based UAV Monitoring and Management System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Method for Vertical Acceleration Noise Suppression of a Thrust-Vectored VTOL UAV

1
Department of Automation, Tsinghua University, Beijing 100084, China
2
Department of Energy Resources Engineering, Stanford University, Stanford, CA 94305, USA
*
Author to whom correspondence should be addressed.
Sensors 2016, 16(12), 2054; https://doi.org/10.3390/s16122054
Submission received: 7 September 2016 / Revised: 24 November 2016 / Accepted: 29 November 2016 / Published: 2 December 2016
(This article belongs to the Special Issue UAV-Based Remote Sensing)

Abstract

:
Acceleration is of great importance in motion control for unmanned aerial vehicles (UAVs), especially during the takeoff and landing stages. However, the measured acceleration is inevitably polluted by severe noise. Therefore, a proper noise suppression procedure is required. This paper presents a novel method to reduce the noise in the measured vertical acceleration for a thrust-vectored tail-sitter vertical takeoff and landing (VTOL) UAV. In the new procedure, a Kalman filter is first applied to estimate the UAV mass by using the information in the vertical thrust and measured acceleration. The UAV mass is then used to compute an estimate of UAV vertical acceleration. The estimated acceleration is finally fused with the measured acceleration to obtain the minimum variance estimate of vertical acceleration. By doing this, the new approach incorporates the thrust information into the acceleration estimate. The method is applied to the data measured in a VTOL UAV takeoff experiment. Two other denoising approaches developed by former researchers are also tested for comparison. The results demonstrate that the new method is able to suppress the acceleration noise substantially. It also maintains the real-time performance in the final estimated acceleration, which is not seen in the former denoising approaches. The acceleration treated with the new method can be readily used in the motion control applications for UAVs to achieve improved accuracy.

1. Introduction

Vertical takeoff and landing (VTOL) unmanned aerial vehicles (UAV) have attracted the attention of many researchers [1,2]. For VTOL UAVs, the vertical motion control design is of great importance, especially during the takeoff and landing stages [3,4]. This is because a motion control failure can result in severe damage to the UAV. To achieve well-performed vertical motion control, high-quality vertical acceleration feedback should be used because it enables the controllers to take action before the unexpected dynamic manifests itself in velocity and height [5,6,7].
However, the acceleration measurement can be severely polluted by noise in real applications, which is a major challenge of acceleration feedback control [8]. A lot of research has been done to reduce the acceleration noise.
Sun et al. developed an infinite impulse response (IIR) digital lowpass filter to remove the high frequency noise in acceleration measurement [9]. The parameters of the lowpass filter were adopted according to the characteristics of the acceleration noise. An experiment was conducted to verify the performance of the lowpass filter. Because of its simplicity, a lowpass filter has been widely used in denoising applications [10]. Lu et al. developed a procedure which consists of an IIR lowpass filter and a Kalman filter to attenuate the influence of sensor noise and external disturbance [11]. The lowpass filter is applied at first to eliminate the high frequency noise of the measured acceleration. The Kalman filter is then used to further reduce the noise. The noise reduction performance was verified by experiment. The standard deviation of acceleration noise was reduced from 0.33 m/s 2 to 0.02 m/s 2 . In addition, El-Sheimy et al. introduced the wavelet filter to reduce the acceleration noise [12]. The six-level wavelet decomposition is used to eliminate the high frequency noise from the low frequency signal. A significant improvement in the quality of the acceleration signal was achieved. The standard deviation of the acceleration noise was reduced by about 99%.
Though the above methods are able to reduce acceleration noise effectively, these studies were all taken in a constrained environment in which the acceleration is very stable with only small variations. For acceleration, which varies significantly with time, the noise reduction becomes much more challenging. Kownacki investigated such problems using the Kalman filter, though the model was built with the assumption of constant acceleration [13]. A conflict between the output signal noise level and filter response rate was discussed. This discussion indicates that the filter response rate has to be sacrificed in order to achieve effective noise suppression performance. In fact, this dilemma also exists in all of the above noise reduction methods. Hebbale and Ghoneim exploited constant jerk (the derivative of acceleration) model in the design of a Kalman filter to reduce the acceleration noise [14]. Nevertheless, this method also suffers from the time delay just like the other approaches.
The above-mentioned methods have been used for UAV acceleration noise reduction. For example, Abellanosa et al. used lowpass filter to suppress the acceleration noise for position estimation of a Quadcopter UAV [15]. Quadri et al. applied a Kalman filter to reduce acceleration noise for attitude estimation of a multi-rotor UAV [16]. Related work can also be found in [17,18]. In those studies, the accelerations were mainly used for position and attitude estimation for UAVs. The presence of time delay in the existing denoising methods prevents the usage of acceleration feedback control for UAVs.
To overcome the problems in the previous denoising methods, we develop a new frame for acceleration noise suppression of a thrust-vectored UAV. Instead of the direct filter used in previous work, the new procedure combines the thrust information with the accelerometer measurement to obtain a more accurate acceleration estimate. The new method is able to suppress the acceleration noise without introducing time delay. A thrust-vectored tail-sitter UAV takeoff experiment is conducted to verify the effectiveness of our method. To the best of our knowledge, this is the first attempt to estimate UAV acceleration in this way.
This paper is organized as follows. In Section 2, the configuration of the thrust-vectored VTOL UAV and the experiment platform is introduced. Section 3 presents the detailed analysis of the statistical characteristics of the acceleration noise, which is the prerequisite for the noise suppression approach. In Section 4, we build an engine thrust model and a thrust vector model based on ground experiments which are used to calculate the vertical thrust when UAV is flying. The noise reduction procedure is then described in Section 5, which includes a Kalman filter to estimate UAV mass and a data fusion step to obtain the minimum variance estimate of acceleration. In Section 6, we present the results for the UAV takeoff experiment. Two existing methods are tested as well for comparison purpose. A section to summarize the work is finally presented.

2. UAV Configuration and Experiment Platform

Tsinghua University has developed a new prototype tail-sitter VTOL UAV, which is equipped with a thrust-vectoring engine system as shown in Figure 1. The UAV has a length of 2 m, a wingspan of 1.6 m, and the gross weight of the UAV (including the thrust vectoring system) is about 24 kg. The thrust vectoring system includes two micro turbine engines with independent vectoring nozzles. The turbine engine used is Jetcat P200. P200 is a single-spool axial flow turbine engine with a net weight of 2.4 kg and a maximum thrust of 230 N which corresponds to a maximum rotor speed of 112,000 revolutions per minute (RPM). With two turbine engines, this UAV has a maximum thrust of 460 N which enables an approximately 46 kg maximum takeoff mass. The detailed technical parameters of the engine are available in [19].
The electronic system of the UAV includes a microcomputer, an inertial measurement unit (IMU), a Hall sensor, a laser range finder and other auxiliary communication and measurement instruments. The microcomputer board with STM32F407 microprocessor is used to perform data processing and give commands to the thrust vectoring system in order to keep the UAV stable. The IMU (ADIS16488A from Analog Devices) provides the attitude and acceleration information to the microcomputer with a sample frequency of 25 Hz. It is installed in the center of gravity of the UAV to prevent the rotation of the UAV body from disturbing acceleration measurement. The Hall sensor is used to measure the rotor speed with very high accuracy (±10 RPM for rotor speed at over 85,000 RPM). The laser range finder (DLSB-15 from Dimetix) with a sample frequency of 5 HZ and a measurement error within ±1.5 mm is used to measure the height from the ground for the UAV.
To ensure the safety of the experiment (for both the staff and UAV body) and also for the convenience of the experiment, the UAV does not directly take off and land on the ground. Similarly to the UAV VTOL experiment study in [20], it is hung with a cable under a gantry crane which is shown in Figure 2a. The height of the gantry crane is 7.2 m, and the width is 6.8 m. To avoid the ground effect influencing the thrust performance, the tail of the UAV is hung at about 1.5 m high. As mentioned previously the UAV itself has a length of 2 m, therefore the operating window for the UAV is about 3.5 m. During the experiment, the UAV height is monitored closely to ensure it does not hit the crane. The UAV hanging under the crane is shown in Figure 2b. The experiment is conducted in calm days with wind speed less than 1 km/h to minimize the external disturbance.

3. Analysis of Acceleration Noise

The acceleration measured with IMU is inevitably polluted by noise. The noise can be caused by many reasons, for example, the external disturbances (such as wind), and the vibration of the turbine engines. To develop a noise reduction procedure, the characteristics of the noise need firstly to be investigated.
In this work, only the vertical acceleration is analyzed because it is of greater importance for UAV motion control. The vertical acceleration is not directly measured in the strap-down inertial system. Instead, the acceleration is measured in the UAV body coordinate system. This measured acceleration needs to be transformed from the UAV body coordinate system to the inertial coordinate system. Previous work on coordinate transformation for tail-sitter UAVs has been done by other researchers [21]. It is directly used in this paper without repetitive description. In this work, all the acceleration data presented is in the inertial coordinate system (after transformation).
The measurements of UAV vertical acceleration at three different engine states are shown in Figure 3. Figure 3a shows the acceleration before the engines start, Figure 3b shows the acceleration when the engines work at an idle speed (33,000 RPM), while Figure 3c shows the acceleration when the engines work close to the takeoff speed (85,000 RPM). It should be noted that the engine system does not reach the takeoff thrust at the rotor speed of 85,000 RPM, therefore the UAV remains stationary for all three states.
Because the IMU measures the proper acceleration (physical acceleration) instead of coordinate acceleration (rate of change of velocity), the measured acceleration should be the gravitational acceleration when the UAV stays stationary. As shown in Figure 3a the acceleration measurement is very close to 9.8 m/s 2 with small-scale noise. In contrast, the acceleration noise becomes more significant after the engines start and it grows as the rotor speed increases from 33,000 to 85,000 RPM. This is because of the vibration in the engines due to the rotation of the rotor and the disturbance of the gas flow. Note that the rotor speed of 85,000 RPM is close to the critical rotor speed for UAV takeoff.
The mean and variance for each of the three acceleration measurements are shown in Table 1. It can be seen that the mean values are very close to each other. The slight difference between them is because of the finite amount of data measured during the experiments. The largest difference is only 0.009 m/s 2 , which is negligible compared with the gravitational acceleration. In this work, gravitational acceleration is adopted as 9.797 m/s 2 which is the average of the three mean values in Table 1. In contrast, the computed variance values differ from each other significantly. When engines are shut down (before start), the variance is very small and it increases dramatically when engine rotor speed is of 33,000 RPM or 85,000 RPM. This is consistent with the observation in Figure 3, since the variance in acceleration represents the intensity of the noise.
Many filtration and data fusion algorithms have intrinsic assumptions on the statistical characteristics of the noise. More detailed investigation of the noise is therefore needed. In this work, the acceleration noise in measured data is defined as:
w = a m a ¯ m
where w represents the noise, a m is the measured acceleration and a ¯ m is the mean of the measured acceleration.
We first compute the normalized autocorrelation, designated R ( k ) , of the acceleration noise using Equation (2):
R ( k ) = 1 N σ 2 i = 1 N w ( i ) w ( i + k )
where N is the number of sampled data points, σ 2 is the variance of the noise, w ( i ) and w ( i + k ) are the noise values at sampled data points i and i + k respectively, and k is the lag of the autocorrelation. The acceleration noise for engine rotor speed of 85,000 RPM (as shown in Figure 3c) is used for analysis because it corresponds to a state of greater importance (near the takeoff speed).
The result of Equation (2) is shown in Figure 4 with the x-axis indicating the lag and the y-axis indicating the normalized autocorrelation of the sampled data. According to Figure 4, the autocorrelation of the acceleration noise can be regarded as an impulse function with high fidelity, which indicates that the acceleration noise can be treated as white noise.
In addition, Figure 5 shows the histogram of the acceleration noise with a normal distribution curve fit (purple line). It can be seen that the normal distribution fits the experimental results very well, which suggests that the noise is of the Gaussian type.
In addition to the above analysis, Woodman [22] and Titterton [23] also stated that it is standard practice to assume that random errors follow a Gaussian distribution when modelling errors in acceleration measurement. Therefore, it is reasonable to assume the accelerometer measurement noise is a white Gaussian noise. This assumption is widely used in the literature, for example, in [11,13,16] mentioned before. The same assumption is used in this work as well.

4. Modeling of Engine Thrust System

During the takeoff and landing stages, the aerodynamic forces can be ignored because of the low air speed relative to the UAV. The engine thrust is the only non-gravitational force exerts on the UAV. The UAV motion dynamics in the vertical direction can therefore be written as:
F v m g = m a
where F v denotes the vertical thrust of the engine system, m is the mass of the UAV, g is the gravitational acceleration, and a is the vertical coordinate acceleration with the positive direction upward. By introducing the proper acceleration a p , Equation (3) can be rewritten as:
a p = g + a = F v / m
This proportional relationship (Newton’s Second Law) between the proper acceleration and the vertical thrust is demonstrated in Equation (4). It inspires us to incorporate the thrust information into the acceleration estimation. However, because the engine thrust cannot be directly measured during flight, mathematical models are needed to calculate the thrust. In this section, we develop models to compute the thrust and thrust deflection angle based on measurements that can be obtained when the UAV is flying.
Note that the noise reduction method developed in this work (presented in Section 5) is independent of the engine thrust system model built here. The method can be readily used to reduce acceleration noise for other types of UAV platforms with different engines for which different thrust system models are needed. The thrust system model described in this section is for Jetcat P200 turbine engine used in our UAV platform.

4.1. Engine Thrust Model

Previous work has demonstrated that a single input single output model can be used to characterize the relationship between the thrust and rotor speed for micro turbine engines [24]. In this work, a similar relationship is built to model the P200 turbine engine. Since the engine thrust cannot be measured during flight, the ground experiments are conducted.
Figure 6 shows the settings for engine thrust ground experiment. A turbine engine and a vectoring nozzle are fixed on a precise balance which can measure the triaxial forces. The detailed descriptions for the turbine engine ground experiment can be found at [24]. For the thrust experiment, the deflection angle of the vectoring nozzle is set to be zero, meaning the nozzle axis is aligned with the engine axis. During the experiment, the rotor speed is fixed at a constant value and the corresponding thrust is measured at the steady state. The rotor speed is then adjusted to different values and the thrust measurements are taken to build the model.
The thrust experiment results for a single P200 engine at steady state are shown in Figure 7 as the blue points for rotor speed ranges from 60,000 RPM to 94,000 RPM. A cubic polynomial function is used to fit the experimental data, and it gives us the red curve as shown in Figure 7. The mathematical expression for the cubic polynomial relationship is written as:
F t = 439 . 61 Ω 3 499 . 46 Ω 2 + 293 . 67 Ω 55 . 26
where F t represents the engine thrust with the unit of Newton (N) and Ω represents the rotor speed with the unit of 10 5 RPM. It can be seen that the polynomial function matches the experiment data closely. The mean of the regression residuals is 1 . 2 × 10 4 N and the standard deviation of the residuals is 0.67 N.
The above polynomial relationship is validated by comparing the experimental thrust with the computed thrust (using Equation (5)). The comparison is shown in Table 2 for the rotor speeds of 96,000, 98,000 and 100,000 RPM which are not used in obtaining Equation (5). The relative error, designated E F , is computed as:
E F = F t m F t e F t e
where F t m is the engine thrust computed using the model, and F t e is the engine thrust measured in the experiment. The close agreement between the measured and computed thrusts and the very small (0.3%) relative errors shown in Table 2 demonstrate high accuracy of the thrust model. Different experimental data (rather than the blue points shown in Figure 7) can be used to fit the cubic equation as long as the relative errors shown in Table 2 are within an acceptable range (0.5% used in this work). Detailed discussion about the cubic thrust model can be found in an earlier work by Tsinghua University [24].
For real problems, the engines typically operate with time-varying rotor speed, especially during the takeoff and landing stages. Therefore, the polynomial thrust model needs to be validated by a dynamic thrust experiment as well. The time-varying rotor speed as shown in Figure 8a is used in the dynamic experiment. The range for this time-varying rotor speed is from 45,000 to 95,000 RPM which covers the rotor speed for UAV takeoff and landing. The corresponding thrust computed using the polynomial model is then compared with the thrust measured during the experiment in Figure 8b. A close match is observed which again confirms the accuracy and effectiveness of the thrust model.
As mentioned before, the engine thrust cannot be measured during flight, but the rotor speed can be measured by Hall sensor with very high accuracy. The measured rotor speed will serve as the input in the engine thrust model (Equation (5)) to compute the thrust.

4.2. Thrust Vector Model

During flight, the attitude of the UAV is controlled by the engine thrust vectoring system to ensure the stability of the UAV body (with the UAV body nose pointing up vertically). The thrust vectoring system generates triaxial moments by deflecting the vectoring nozzles. Consequently, the engine thrust is not exactly aligned with the engine axis due to the nozzle deflection. The engine thrust model (Equation (5)) developed in the previous section is used to compute the overall thrust. To obtain the vertical component of the overall thrust, the deflection angle of the engine thrust is needed.
We are not able to measure the deflection angle of the thrust directly when the UAV is flying. However, the deflection angle is controlled by the microcomputer and the command of the deflection angle can be recorded. The thrust vectoring system is well calibrated to ensure the thrust deflection angle is the same as the deflection command. It is validated by performing a thrust vectoring experiment on the ground.
The same experimental settings shown in Figure 6 are used. The engine rotor speed for this experiment is fixed at 84,140 RPM which corresponds to a total thrust of 100 N. Since the nozzle is axisymmetric, only the experimental results for yaw movement are presented. During the experiment, the sinusoidal deflection angle command with an amplitude of 10 degree (shown in Figure 9 as the purple curve) is used and the corresponding thrust vector is measured. The deflection angle of the thrust, designated α F , is calculated using Equation (7) as:
α F = t a n 1 F t y b F t x b
where F t x b and F t y b are the measured thrusts in the x and y directions in the UAV body coordinate system. The UAV coordinate is defined as: the x-axis points out of the tailsitter nose, the y-axis points out the right wing and the z-axis points out the bottom of the fuselage.
The actual thrust deflection angle ( α F ) derived from the thrust vector measurement using Equation (7) is compared with the deflection command in Figure 9 using the blue curve. A close match between the two curves is observed. This suggests that the deflection command can be used as the actual thrust deflection angle with sufficient accuracy. Therefore, the thrust vector model can be written as:
F t x b F t y b F t z b = F t 1 + t a n 2 α + t a n 2 β 1 t a n β t a n α
where α and β are the pitch and yaw deflection commands, F t is the total engine thrust which is estimated using engine thrust model described in the previous section, F t x b , F t y b and F t z b are the engine thrusts in the x, y and z directions respectively. Note that the x, y and z directions are in the UAV coordinate system, not in the inertial coordinate system.
By combining the engine thrust model (Equation (5)) and thrust vector model (Equation (8)), we are able to obtain the thrust vector using the measured rotor speed and the deflection commands. The coordinate transformation is again used to obtain the vertical thrust, which is indicated by F v c .
Note that there are two engines in the propulsion system. During the flight, the two engines receive the same rotor speed command, which means they have the same overall thrust. But the deflection angles are controlled independently. Therefore, the vertical thrusts generated by the two engines, designated F v c 1 and F v c 2 , need to be computed separately. The total vertical thrust of the engine system is the sum of F v c 1 and F v c 2 :
F v c = F v c 1 + F v c 2
Though the modeled results are of good quality, small modeling errors are inevitable. In addition, the measurement noise is also present. To incorporate these effects into the vertical thrust computation, we assume the error of the computed vertical thrust is white Gaussian noise. The genuine vertical thrust is thus expressed as:
F v = F v c + w F
where F v is the genuine vertical thrust, F v c is the computed vertical thrust, and w F is the white Gaussian noise for vertical thrust with the variance of σ F 2 .

5. Methodology for Acceleration Noise Suppression

By defining λ = 1 m in which m denotes the UAV mass, we can rewrite Equation (4) as:
a m = λ F v + w m
where F v is the vertical thrust force and a m = a + g + w m represents the measurement of the IMU with w m is the vertical acceleration noise which is white Gaussian noise with a variance of σ a m 2 . By combining with Equation (10), the above equation can be rewritten as:
a m = λ ( F v c + w F ) + w m = λ F v c + ( λ w F + w m )
Since λ is a constant (the inverse of the UAV mass), λ w F is white Gaussian noise with the variance of λ 2 σ F 2 .
We then define a new term, w, as:
w = λ w F + w m
where w represents the overall noise which combines the measured acceleration noise and the thrust modeling error. Because both λ w F and w m are white Gaussian noise, w is white Gaussian noise as well. In addition, λ w F is independent of w m due to the independence between w F and w m . According to the above analysis, w can be treated as white Gaussian noise with the variance, designated σ 2 , written as:
σ 2 = λ 2 σ F 2 + σ a m 2
We take λ as the state variable, a m as the observation variable and F v c as the input parameter, the discrete-time state space model can be written as:
λ ( k + 1 ) = λ ( k ) a m ( k ) = λ ( k ) F v c ( k ) + w ( k )
The standard Kalman filter can be used to estimate the λ by combining the computed thrust and measured acceleration information. The updating process is as follows:
K λ ( k ) = P λ ( k 1 ) F v c ( k ) [ F v c ( k ) P λ ( k 1 ) F v c ( K ) + σ 2 ( k ) ] 1 P λ ( k ) = [ 1 K λ ( k ) F v c ( k ) ] P λ ( k 1 ) λ ^ ( k ) = λ ^ ( k 1 ) + K λ ( k ) [ a m ( k ) F v c ( k ) λ ^ ( k 1 ) ]
where λ ^ ( k 1 ) and λ ^ ( k ) are the estimates of λ at sampled data points k 1 and k respectively, K λ ( k ) is the Kalman filter gain, and P λ ( k ) is the variance of λ ^ ( k ) . Note that the variance of w ( k ) , designated σ 2 ( k ) , is estimated as follows in the calculation process:
σ 2 ( k ) = λ ^ 2 ( k 1 ) σ F 2 + σ a m 2
Here we assume the variance of the measured acceleration noise during UAV takeoff and landing is the same as the acceleration variance when the engine rotor speed is at 85,000 RPM. This is because that the rotor speed of 85,000 RPM is very close to speed of the engines during takeoff and landing. Therefore, σ a m = 0 . 351 (from Table 1) is used in this work. The variance of vertical thrust error is assumed to be σ F 2 = 10 N 2 .
With the estimated λ and calculated vertical thrust F v c , vertical acceleration can be obtained as:
a c ( k ) = λ ^ ( k ) F v c ( k )
where a c ( k ) is the calculated acceleration at sampled data point k. By combing with Equation (10), the above equation can be rewritten as:
a c ( k ) = λ ^ ( k ) [ F v ( k ) w F ( k ) ] = λ ^ ( k ) F v ( k ) λ ^ ( k ) w F ( k )
We then investigate the expectation of a c ( k ) as:
E { a c ( k ) } = E { λ ^ ( k ) F v ( k ) } E { λ ^ ( k ) w F ( k ) } = λ F v ( k ) ρ λ ^ ( k ) , w F ( k ) σ λ ^ ( k ) σ F
where ρ λ ^ ( k ) , w F ( k ) represents the correlation coefficient between λ ^ ( k ) and w F ( k ) , and the derivation of ρ λ ^ ( k ) , w F ( k ) is shown in Equation (A1). Note that the absolute value of the correlation coefficient does not exceed 1. Because σ λ ^ ( k ) approaches zero as k tends to infinity, the expectation of a c ( k ) approaches the genuine vertical acceleration, λ F v ( k ) , as k tends to infinity according to Equation (20). Therefore a c ( k ) is the asymptotically unbiased estimate of the vertical acceleration.
According to Equation (19), the variance of a c ( k ) is then given by:
σ a c 2 ( k ) = F v 2 ( k ) σ λ ^ ( k ) 2 + λ ^ ( k ) 2 σ F 2 2 ρ λ ^ ( k ) , w F ( k ) λ ^ ( k ) F v ( k ) σ λ ^ ( k ) σ F
To obtain the minimum variance estimate of the vertical acceleration, the calculated acceleration a c ( k ) and the measured acceleration a m ( k ) are fused as :
a f ( k ) = η ( k ) a c ( k ) + [ 1 η ( k ) ] a m ( k )
where a f ( k ) is the acceleration fusion result, η ( k ) is the weighting coefficient with 0 < η ( k ) < 1 . It is obvious that a f ( k ) is asymptotically unbiased due to the asymptotically unbiasedness of a c ( k ) and the unbiasedness of a m ( k ) . The variance of a f ( k ) is given by:
σ a f ( k ) 2 = η ( k ) 2 σ a c ( k ) 2 + [ 1 η ( k ) ] 2 σ a m 2 + 2 η ( k ) [ 1 η ( k ) ] ρ a c ( k ) , a m ( k ) σ a c ( k ) σ a m
where ρ a c ( k ) , a m ( k ) is the correlation coefficient between a c ( k ) and a m ( k ) which is derived in Equation (A2). The variable η ( k ) is adopted as follows to achieve the minimum variance estimate:
η ( k ) = σ a m 2 ρ a c ( k ) , a m ( k ) σ a c ( k ) σ a m σ a c ( k ) 2 + σ a m 2 2 ρ a c ( k ) , a m ( k ) σ a c ( k ) σ a m
The flowchart of the noise reduction procedure developed in this work is shown in Figure 10. At first, the engine system model (including engine thrust model and thrust vector model in Section 4) needs to be built. It is used to obtain the UAV vertical thrust when it is in flight. Then the Kalman filter (Equation (16)) is applied to estimate UAV mass by combining the information of engine vertical thrust and the IMU measured acceleration. The estimated UAV mass can then be used together with the vertical thrust to compute acceleration using Newton’s Second Law. A data fusion step (Equation (22)) is finally applied to obtain the minimum variance estimate of the UAV acceleration. This overall procedure can be used for any type of UAVs for noise reduction purpose. The only change will be in the engine system model if the method is applied to different UAV platforms.

6. Experiment Result

To assess the performance of the new method for UAV vertical acceleration noise suppression, a vertical takeoff experiment is conducted. The data record began at the time when the UAV has positive (upward) elevation change which is indicated by the height measurement from the laser range finder as shown in Figure 11a. The UAV ascended from the initial position of about 1.5 m to around 3.2 m at t = 2 . 8 s. This corresponds to a maximum height gain of about 1.7 m. After reaching the highest point, the UAV started to descend, though it had not returned to its original position for the period that we measured. The rotor speed measured during the experiment is shown in Figure 11b. As can be seen, the rotor speed increased for t < 1 s, and dropped quickly during 1 < t < 2 s. The rotor speed reached a stable value around 86,000 RPM eventually. The corresponding vertical acceleration measurement from IMU is shown in Figure 11c. The same as in Figure 3, the acceleration measurement is heavily polluted by the noise, though the similar trend with the rotor speed is observed.
We then subtract the gravitational acceleration from the IMU measured acceleration shown in Figure 11c. This gives us the coordinate acceleration (rate of change of velocity) which is displayed in Figure 11d. It can be clearly seen that this UAV experiment includes both the upward acceleration period (positive coordinate acceleration) and downward acceleration period (negative coordinate acceleration). Therefore, it is sufficient to validate the noise reduction method developed in this work.
We apply the engine thrust model and thrust vector model developed in Section 3 to the experimental data. The calculated vertical thrust is shown in Figure 12. It is then used to estimate λ using the Kalman filter (Equation (16)). As shown in Figure 13, the value of λ quickly converges to 0.0404 kg 1 , which means that the mass of the UAV (1/λ) is about 24.75 kg. Note that, the variable λ is unknown during flight, though it can be considered as a constant within a short period of time (for example, during the takeoff and landing stages). This is because that the UAV mass is impacted by the fuel consumption, which is unknown when the UAV is flying. Therefore, the Kalman filter (Equation (16)) is required to estimate λ.
The final acceleration fusion result is shown in Figure 14 as the red line. For comparison purpose, the measured vertical acceleration is shown again in Figure 14. It can be seen clearly that the acceleration fusion result suppresses the noise in the measured acceleration and meanwhile preserves the trend in measurements. To quantify that, we calculate the variance of the acceleration computed with data fusion method using Equation (23). The variance is shown in Figure 15. Initially the variance is of similar magnitude as the variance of the measured acceleration (as shown in Table 1). It drops quickly as more data is used in the fusion calculation and becomes stable at a value close to 0.02 m 2 /s 4 after about 2 s.
For the purpose of comparison, an infinite impulse response (IIR) lowpass filter developed in [9] and the Kalman filter developed in [13] are also used to denoise the measured acceleration. The parameters of the IIR filter and the Kalman filter are adopted to achieve similar denoising performance to the data fusion method developed in this work.
The IIR lowpass filter is shown as follows:
a l ( k ) = b 0 a m ( k ) + b 1 a m ( k 1 ) + b 2 a m ( k 2 ) + b 3 a m ( k 3 ) + b 4 a m ( k 4 ) + b 5 a m ( k 5 ) c 1 a l ( k 1 ) c 2 a l ( k 2 ) c 3 a l ( k 3 ) c 4 a l ( k 4 ) c 5 a l ( k 5 )
where a l ( k ) and a m ( k ) are the filtered and measured accelerations for sampled data point k, respectively. The coefficients used here are:
b 0 = 0 . 0013 , b 1 = 0 . 0064 , b 2 = 0 . 0128 , b 3 = 0 . 0128 , b 4 = 0 . 0064 , b 5 = 0 . 0013 c 1 = 2 . 9754 , c 2 = 3 . 8060 , c 3 = 2 . 5453 , c 4 = 0 . 8811 , c 5 = 0 . 1254
The state space model used in the Kalman filter developed in [13] is shown as:
a s ( k + 1 ) = a s ( k ) + w s ( k ) a m ( k ) = a s ( k ) + w m ( k )
where a k ( k ) and a k ( k + 1 ) are the acceleration states at sampled point k and k + 1 , w s ( k ) is the process noise, a m ( k ) is the observed acceleration with observation noise indicated by w m ( k ) . Note that w s ( k ) and w m ( k ) are mutually independent white Gaussian noise. Their variances used here are adopted as σ w s 2 = 0 . 01 and σ w m 2 = 0 . 351 , respectively.
The comparison is shown in Figure 16. At a very early time value ( t < 0 . 2 s), the results using the Kalman filter and the new data fusion method are very close. As more data is included into the calculation, they start to deviate and an obvious time delay is seen in the Kalman filter result for t < 2 . 5 s. In contrast, the lowpass result differs from the fusion result from the very beginning, and the time delay in the lowpass result is very obvious as well. At the late stage ( t > 2 . 5 s), when the acceleration becomes relatively stable (due to stable engine rotor speed shown in Figure 11b), the results using the three approaches become close.
The difference between each of the three computed accelerations and the measured acceleration, designated a d , is computed as:
a d = a c a m
where a c is the computed acceleration, and a m is the measured acceleration. The acceleration differences are shown in Figure 17. Dashed lines are used for Kalman filter and data fusion results for better visualization. We can see that the acceleration differences are of high amplitude due to the significant noise in the measured acceleration.
A nine-order mean filter is used to extract the low frequency component of the acceleration differences shown in Figure 17. The nine-order mean filter is expressed as:
a ¯ d ( k ) = 1 9 i = k 4 k + 4 a d ( i ) , for 4 < k < N 4
where a ¯ d ( k ) is the low frequency component result at sampled data point k. The results are shown in Figure 18, from which we can clearly see the difference between the new method and the two existing methods. The lowpass filter and Kalman filter provide similar results and both of them have apparent non-zero values in their low frequency components (greater than zero for t < 0 . 9 s and less than zero during 1 < t < 2 . 2 s). This is because of the time delay introduced by the filtering process. The time delay is due to the weighting of the past time acceleration measurements in the filtered results in both the lowpass and Kalman filters. In contrast, the results using the new method closely distribute around zero without any obvious non-zero trends. This is because that the new method uses the past time information only for the estimation of the variable λ. As shown in Figure 13, the λ converges very quickly to a constant value during the takeoff stage. Therefore, the new method does not introduce any time delay in acceleration computation.
Note that the time delay can be reduced for lowpass filter and Kalman filter by adjusting their parameters, but the noise suppression performance of the two methods will be severely degraded. Compared with the lowpass filter and Kalman filter, the newly developed method aided with the thrust information have much better real time performance with strong ability of noise suppression.

7. Conclusions

In this work, we developed a novel method to reduce the noise in acceleration measurement for a thrust-vectored VTOL UAV. The method combines the information in engine thrust with the information in the measured acceleration. To do that, a Kalman filter is first applied to estimate the UAV mass, which is used to compute acceleration together with the vertical thrust. This estimated acceleration is further fused with the measured acceleration to obtain the minimum variance acceleration estimate.
To test the effectiveness of the newly developed approach, a vertical takeoff experiment for the UAV was performed. The new method was used to compute the vertical acceleration using the data collected from the experiment. The results demonstrated that the new method has very good performance in terms of noise reduction. The variance of the vertical acceleration is reduced by about 95%. In addition, the acceleration calculated using the new approach shows no time delay, which enables its usage in acceleration feedback control for UAV takeoff and landing (or for any other situations with fast-changing accelerations). In contrast, two existing approaches based on IIR and Kalman filter suffer from significant time delay which limits their usability for real problems, although they can achieve similar denoising performance.
In addition, we developed numerical models for engine thrust and thrust vector which are important components in the entire noise suppression procedure. A polynomial function is used to characterize the relationship between the rotor speed and the engine thrust. Though the model was built based on steady-state ground experiment, it is able to predict accurately the thrust in dynamic condition (with time-varying rotor speed). The thrust deflection angle was found to be consistent with the deflection angle command, which is then used to compute the engine thrust in the vertical direction.
We analyzed the statistical characteristics of the acceleration noise. It was found that the amplitude of the noise grows with the increasing rotor speed. We also demonstrated that the acceleration noise can be treated as white Gaussian noise based on the autocorrelation and histogram of the sampled noise data.
The noise suppression method developed here can be readily applied into the acceleration feedback motion control for VTOL UAVs during a critical period, for example, takeoff and landing stages. Improved precision in motion control is expected, which will be tested in future work.

Acknowledgments

The authors thank Chao Zhang for his experiment advices.

Author Contributions

The method was proposed by Huanyu Li. Huanyu Li and Linfeng Wu designed the experiment. The experiments were performed by Huanyu Li, Linfeng Wu, Yingjie Li. Data analysis was done by Huanyu Li and Hangyu Li, and the paper was mainly written by them as well. Chunwen Li supervised the work. All authors contributed to writing the article.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
UAVUnmanned aerial vehicle
VTOLVertical takeoff and landing
IMUInertial measurement unit
RPMRevolutions per minute
IIRInfinite impulse response

Appendix A

ρ λ ^ ( k ) , w F ( k ) = E { w F ( k ) ( λ ^ ( k ) λ ) } / ( σ λ ^ ( k ) σ F ) = E { w F ( k ) [ λ ^ ( k 1 ) + K λ ( k ) ( a m ( k ) F v c ( k ) λ ^ ( k 1 ) ) λ ] } / ( σ λ ^ ( k ) σ F ) = E { w F ( k ) [ λ ^ ( k 1 ) + K λ ( k ) ( λ F v c ( k ) + λ w F ( k ) + w m ( k ) F v c ( k ) λ ^ ( k 1 ) ) λ ] } / ( σ λ ^ ( k ) σ F ) = E { w F ( k ) [ ( λ λ ^ ( k 1 ) ) ( K λ ( k ) F v c ( k ) 1 ) + w F ( k ) K λ ( k ) ( λ w F ( k ) + w m ( k ) ) ] } / ( σ λ ^ ( k ) σ F ) = E { K λ ( k ) w F ( k ) [ λ w F ( k ) + w m ( k ) ] } / ( σ λ ^ ( k ) σ F ) = K λ ( k ) λ σ F / σ λ ^ ( k )
ρ a c ( k ) , a m ( k ) = E { w m ( k ) [ λ ^ ( k ) w F ( k ) + λ ^ ( k ) F v ( k ) λ F v ( k ) ] } / ( σ a c ( k ) σ a m ) = E { λ ^ ( k ) w m ( k ) w F ( k ) + w m ( k ) ( λ ^ ( k ) λ ) F v ( k ) } / ( σ a c ( k ) σ a m ) = F v ( k ) E { w m ( k ) [ λ ^ ( k ) λ ] } / ( σ a c ( k ) σ a m ) = F v ( k ) E { w m ( k ) [ λ ^ ( k 1 ) + K λ ( k ) ( a m ( k ) F v c ( k ) λ ^ ( k 1 ) ) λ ] } / ( σ a c ( k ) σ a m ) = F v ( k ) E { w m ( k ) [ ( λ λ ^ ( k 1 ) ) ( K λ ( k ) F v c ( k ) 1 ) + K λ ( k ) w m ( k ) + K λ λ w F ( k ) ] } / ( σ a c ( k ) σ a m ) = F v ( k ) K λ σ a m / σ a c ( k )

References

  1. Johnson, E.N.; Turbe, M.A. Modeling, control, and flight testing of a small-ducted fan aircraft. J. Guid. Control Dyn. 2006, 29, 769–779. [Google Scholar] [CrossRef]
  2. Kubo, D.; Suzuki, S. Tail-sitter vertical takeoff and landing unmanned aerial vehicle: Transitional flight analysis. J. Aircr. 2008, 45, 292–297. [Google Scholar] [CrossRef]
  3. Tomic, T.; Schmid, K.; Lutz, P.; Domel, A.; Kassecker, M.; Mair, E.; Grixa, I.L.; Ruess, E.; Suppa, M.; Burschka, D. Toward a fully autonomous UAV: Research platform for indoor and outdoor urban search and rescue. IEEE Robot. Autom. Mag. 2012, 19, 46–56. [Google Scholar] [CrossRef]
  4. Lei, X.; Li, J. An adaptive altitude information fusion method for autonomous landing processes of small unmanned aerial rotorcraft. Sensors 2012, 12, 13212–13224. [Google Scholar] [CrossRef] [PubMed]
  5. Peddle, I.K. Acceleration Based Manoeuvre Flight Control System for Unmanned Aerial Vehicles; Stellenbosch University: Stellenbosch, South Africa, 2008. [Google Scholar]
  6. Boyle, D.P.; Chamitoff, G.E. Autonomous maneuver tracking for self-piloted vehicles. J. Guid. Control Dyn. 1999, 22, 58–67. [Google Scholar] [CrossRef]
  7. Wang, X. Takeoff/landing control based on acceleration measurements for VTOL aircraft. J. Frankl. Inst. 2013, 350, 3045–3063. [Google Scholar] [CrossRef]
  8. Gan, Y.; Sui, L.; Wu, J.; Wang, B.; Zhang, Q.; Xiao, G. An EMD threshold de-noising method for inertial sensors. Measurement 2014, 49, 34–41. [Google Scholar] [CrossRef]
  9. Sun, F.; Sun, W. Mooring alignment for marine SINS using the digital filter. Measurement 2010, 43, 1489–1494. [Google Scholar] [CrossRef]
  10. Parks, T.W.; Burrus, C.S. Digital Filter Design; Wiley-Interscience: Hoboken, NJ, USA, 1987. [Google Scholar]
  11. Lu, S.; Xie, L.; Chen, J. New techniques for initial alignment of strapdown inertial navigation system. J. Frankl. Inst. 2009, 346, 1021–1037. [Google Scholar] [CrossRef]
  12. El-Sheimy, N.; Nassar, S.; Noureldin, A. Wavelet de-noising for IMU alignment. IEEE Aerosp. Electron. Syst. Mag. 2004, 19, 32–39. [Google Scholar] [CrossRef]
  13. Kownacki, C. Optimization approach to adapt Kalman filters for the real-time application of accelerometer and gyroscope signals’ filtering. Dig. Signal Process. 2011, 21, 131–140. [Google Scholar] [CrossRef]
  14. Hebbale, K.V.; Ghoneim, Y.A. A speed and acceleration estimation algorithm for powertrain control. In Proceedings of the 1991 American Control Conference, Boston, MA, USA, 26–28 June 1991; pp. 415–420.
  15. Abellanosa, C.B.; Lugpatan, R.P.J.; Pascua, D.A.D. Position Estimation using Inertial Measurement Unit (IMU) on a Quadcopter in an Enclosed Environment. Int. J. Comput. Commun. Instrum. Eng. 2016, 3, 332–336. [Google Scholar]
  16. Quadri, S.A.; Sidek, O.; Bin Abdullah, A. A Study of State Estimation Algorithms in an OktoKopter. Int. J. u- e-Serv. Sci. Technol. 2014, 7, 247–266. [Google Scholar] [CrossRef]
  17. Edu, I.R.; Adochiei, F.C.; Grigorie, T.L.; Botez, R.M. Tuning of a Wavelet Filter for Miniature Accelerometers Denoising based Joint Symbolic Dynamics (JSD) Method. INCAS Bull. 2015, 7, 71–81. [Google Scholar]
  18. Nebylov, A.; Sukrit, S.; Arifuddin, F. Perspectives for Development of an Autonomous and Intelligent WIG-Craft and Its Peculiar Control Problems. In Proceedings of the 2009 IFAC “AGNFC” Workshop Proceedings, Samara, Russia, 30 June–2 July 2011; pp. 1–6.
  19. Li, H.; Wu, L.; Li, Y.; Li, C. Identification of turbine engine dynamics with the governor in the loop. In Proceedings of the 2016 IEEE International Conference on Systems, Man, and Cybernetics (SMC), Budapest, Hungary, 9–12 October 2016.
  20. Stone, R.H.; Anderson, P.; Hutchison, C.; Tsai, A.; Gibbens, P.; Wong, K.C. Flight testing of the t-wing tail-sitter unmanned air vehicle. J. Aircr. 2008, 45, 673–685. [Google Scholar] [CrossRef]
  21. Beach, J.M.; Argyle, M.E.; McLain, T.W.; Beard, R.W.; Morris, S. Tailsitter heading estimation using a magnetometer. In Proceedings of the 2014 American Control Conference, Portland, OR, USA, 4–6 June 2014; pp. 91–96.
  22. Woodman, O.J. An Introduction to Inertial Navigation; Technical Reports; UCAMCL-TR-696; University of Cambridge, Computer Laboratory: Cambridge, UK, 2007. [Google Scholar]
  23. Titterton, D.; Weston, J. Strapdown Inertial Navigation Technology, 2nd ed.; The Institution of Engineering and Technology: London, UK, 2004. [Google Scholar]
  24. Yang, J.; Zhu, J. Dynamic modelling of a small scale turbojet engine. In Proceedings of the 2015 European Control Conference (ECC), Linz, Austria, 15–17 July 2015; pp. 2750–2755.
Figure 1. The thrust-vectored tail-sitter unmanned aerial vehicle (UAV).
Figure 1. The thrust-vectored tail-sitter unmanned aerial vehicle (UAV).
Sensors 16 02054 g001
Figure 2. The experiment platform. (a) Gantry crane; (b) Suspended UAV.
Figure 2. The experiment platform. (a) Gantry crane; (b) Suspended UAV.
Sensors 16 02054 g002
Figure 3. Vertical acceleration measurements at different engine states. (a) Vertical acceleration measurement before the engines start; (b) Vertical acceleration measurement when the engines work at idle (33,000 RPM); (c) Vertical acceleration measurement when the engines work close to takeoff (85,000 RPM).
Figure 3. Vertical acceleration measurements at different engine states. (a) Vertical acceleration measurement before the engines start; (b) Vertical acceleration measurement when the engines work at idle (33,000 RPM); (c) Vertical acceleration measurement when the engines work close to takeoff (85,000 RPM).
Sensors 16 02054 g003
Figure 4. Normalized autocorrelation of the noise in measured acceleration.
Figure 4. Normalized autocorrelation of the noise in measured acceleration.
Sensors 16 02054 g004
Figure 5. Histogram of the acceleration noise and normal distribution curve fit.
Figure 5. Histogram of the acceleration noise and normal distribution curve fit.
Sensors 16 02054 g005
Figure 6. Engine thrust ground experiment platform.
Figure 6. Engine thrust ground experiment platform.
Sensors 16 02054 g006
Figure 7. Relationship between the rotor speed and engine thrust in steady state conditions.
Figure 7. Relationship between the rotor speed and engine thrust in steady state conditions.
Sensors 16 02054 g007
Figure 8. Thrust model validation using time-varying rotor speed. (a) Time-varying rotor speed; (b) Measured thrust and modeled thrust for time-varying rotor speed.
Figure 8. Thrust model validation using time-varying rotor speed. (a) Time-varying rotor speed; (b) Measured thrust and modeled thrust for time-varying rotor speed.
Sensors 16 02054 g008
Figure 9. Comparison of the deflection angle command and the actual deflection angle of thrust.
Figure 9. Comparison of the deflection angle command and the actual deflection angle of thrust.
Sensors 16 02054 g009
Figure 10. Flowchart of the noise reduction procedure developed in this work.
Figure 10. Flowchart of the noise reduction procedure developed in this work.
Sensors 16 02054 g010
Figure 11. Experiment results. (a) Experiment height measurement; (b) Experiment rotor speed; (c) Experiment acceleration measurement; (d) Coordinate acceleration (computed from Figure 11c).
Figure 11. Experiment results. (a) Experiment height measurement; (b) Experiment rotor speed; (c) Experiment acceleration measurement; (d) Coordinate acceleration (computed from Figure 11c).
Sensors 16 02054 g011
Figure 12. Calculated vertical thrust.
Figure 12. Calculated vertical thrust.
Sensors 16 02054 g012
Figure 13. Estimate of λ.
Figure 13. Estimate of λ.
Sensors 16 02054 g013
Figure 14. Acceleration fusion result.
Figure 14. Acceleration fusion result.
Sensors 16 02054 g014
Figure 15. Variance of computed acceleration.
Figure 15. Variance of computed acceleration.
Sensors 16 02054 g015
Figure 16. Acceleration fusion result compared with the IIR lowpass filter and Kalman filter results.
Figure 16. Acceleration fusion result compared with the IIR lowpass filter and Kalman filter results.
Sensors 16 02054 g016
Figure 17. Acceleration differences.
Figure 17. Acceleration differences.
Sensors 16 02054 g017
Figure 18. Low frequency component of acceleration differences.
Figure 18. Low frequency component of acceleration differences.
Sensors 16 02054 g018
Table 1. Mean and variance of acceleration measurement.
Table 1. Mean and variance of acceleration measurement.
Engine StatesMean (m/s 2 )Variance (m 2 /s 4 )
Before start9.8010.00022
Idle speed (33,000 RPM)9.7920.193
Rotor speed (85,000 RPM)9.7970.351
Table 2. Comparison of the experimental thrust and the modeled thrust.
Table 2. Comparison of the experimental thrust and the modeled thrust.
Rotor Speed (RPM)96,00098,000100,000
Experimental thrust (N)154.8167.3177.1
Computed thrust (N)154.6166.8177.4
Relative error (%) 0 . 13 0 . 30 0.17

Share and Cite

MDPI and ACS Style

Li, H.; Wu, L.; Li, Y.; Li, C.; Li, H. A Novel Method for Vertical Acceleration Noise Suppression of a Thrust-Vectored VTOL UAV. Sensors 2016, 16, 2054. https://doi.org/10.3390/s16122054

AMA Style

Li H, Wu L, Li Y, Li C, Li H. A Novel Method for Vertical Acceleration Noise Suppression of a Thrust-Vectored VTOL UAV. Sensors. 2016; 16(12):2054. https://doi.org/10.3390/s16122054

Chicago/Turabian Style

Li, Huanyu, Linfeng Wu, Yingjie Li, Chunwen Li, and Hangyu Li. 2016. "A Novel Method for Vertical Acceleration Noise Suppression of a Thrust-Vectored VTOL UAV" Sensors 16, no. 12: 2054. https://doi.org/10.3390/s16122054

APA Style

Li, H., Wu, L., Li, Y., Li, C., & Li, H. (2016). A Novel Method for Vertical Acceleration Noise Suppression of a Thrust-Vectored VTOL UAV. Sensors, 16(12), 2054. https://doi.org/10.3390/s16122054

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