Next Article in Journal
Short Review of EMB Systems Related to Safety Concepts
Next Article in Special Issue
A Novel Radially Closable Tubular Origami Structure (RC-ori) for Valves
Previous Article in Journal
Research on the Influence of Friction Pairs on the Output Characteristics of the Piezoelectric Ultrasonic Actuator
Previous Article in Special Issue
Design of an Intuitive Master for Improving Teleoperation Task Performance Using the Functional Separation of Actuators: Movement and Gravity Compensation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Flying State Sensing and Estimation Method of Large-Scale Bionic Flapping Wing Flying Robot

School of Mechanical Engineering and Automation, Harbin Institute of Technology, Shenzhen 518055, China
*
Author to whom correspondence should be addressed.
Actuators 2022, 11(8), 213; https://doi.org/10.3390/act11080213
Submission received: 22 June 2022 / Revised: 24 July 2022 / Accepted: 25 July 2022 / Published: 31 July 2022
(This article belongs to the Special Issue Mechanism Design and Control for Robotics)

Abstract

:
A large bionic flapping wing robot has unique advantages in flight efficiency. However, the fluctuation of fuselage centroid during flight makes it difficult for traditional state sensing and estimation methods to provide stable and accurate data. In order to provide stable and accurate positioning and attitude information for a flapping wing robot, this paper proposes a flight state sensing and estimation method integrating multiple sensors. Combined with the motion characteristics of a large flapping wing robot, the autonomous flight, including the whole process of takeoff, cruise and landing, is realized. An explicit complementary filtering algorithm is designed to fuse the data of inertial sensor and magnetometer, which solves the problem of attitude divergence. The Kalman filter algorithm is designed to estimate the spatial position and speed of a flapping wing robot by integrating inertial navigation with GPS (global positioning system) and barometer measurement data. The state sensing and estimation accuracy of the flapping wing robot are improved. Finally, the flying state sensing and estimation method is integrated with the flapping wing robot, and the flight experiments are carried out. The results verify the effectiveness of the proposed method, which can provide a guarantee for the flapping wing robot to achieve autonomous flight beyond the visual range.

1. Introduction

Flying in the blue sky like a bird has always been a dream of people. Learning about and utilizing birds’ body structures and flight mechanisms so as to develop bionic flapping aircraft with high maneuverability and low energy consumption has broad application prospects [1,2]. Large-scale bionic flapping wing robots developed by imitating the flight patterns of birds in nature have unique advantages in flight efficiency, wind resistance, and bionic concealment [3,4]. In order to improve the survivability of the flapping wing robot in complex environments, it is extremely important for the position and attitude calculation and state prediction of the flapping wing robot to perform beyond visual range and long-distance flight missions. Since a single sensor is affected by environmental disturbances such as gusts and magnetic field pulses during the flight of the flapping wing robot, the position and attitude estimation information of the flapping wing robot is seriously inconsistent with the actual situation [5,6]. The use of the flying state sensing and estimation method can reasonably use the measurement information of each sensor to complement the advantages and disadvantages, thereby improving the reliability of the flapping wing robot in performing tasks.
At present, there are two main methods for solving the attitude of flapping wing robots. One is to use MEMS (Micro-Electro-Mechanical System, MEMS) (Norwood, MD, USA) sensors (including accelerometers, gyroscopes and magnetometers) to estimate the attitude. This method is suitable for large- and medium-sized flapping wing robots with strong load capacity, such as Bat Bot flapper from University of Illinois at Urbana-Champaign, and pigeon flapper from Northwestern Polytechnic University [7,8]. At the same time, Kamankesh et al. and Yu et al. [9,10] studied the stability characteristics of periodic flapping wing structures and attitude stabilization of the aircraft. Verboom et al. and Hao et al. proposed a method based on periodic filtering to estimate the attitude of the flapping wing robot [11,12]. Yang et al. and Tu et al. [13,14] used Kalman filtering and adaptive attitude update methods to study the attitude information of UAVs (unmanned aerial vehicles). He et al. [15,16] proposed several adaptive control schemes and position estimation methods. In addition, Taha, Bialy, Dou, Zhang and other researchers [17,18,19,20] have also conducted in-depth research on the attitude and control of flapping wing robots. The other is the visual image pose solution method using the ground station for image transmission. This method is relatively common in small flapping wing robots, mainly including RoboBee developed by Harvard University and Delfly micro flapping wing robot developed by Karasek [21,22].
Solving the position and attitude based on flying state sensing and estimation is an important prerequisite for the flapping wing robot to achieve attitude control. Only by accurately obtaining the state parameters of the prototype can the position and attitude be controlled, and then the prototype can be guided to fly according to the predetermined trajectory by adjusting the attitude [23,24]. In a word, solving the pose of the flapping wing robot based on the flying state sensing and estimation method can improve the adaptability of the flapping wing robot to complex environments and improve the reliability of task execution.
However, different from small- and micro-sized flapping wing flying robots, the two wings of the large flapping wing robots have a large proportion of mass and inertia, and the fluctuation of the fuselage center of mass during flight cannot be ignored. The asymmetric flapping or unbalanced force leads to serious attitude deflection, and the existence of delay aggravates the inaccuracy and time inconsistency of perceptual information. Therefore, the traditional flying state sensing and estimation methods for fixed wing and rotor are difficult to provide stable and accurate data, which seriously limits the autonomous flight ability of large bionic flapping wing flying robots. Combined with the long-term accumulation of research on large-scale bionic flapping wing flying robot in the laboratory, and based on its flight characteristics in different stages of takeoff, cruise and landing, this paper proposes a method of flying state sensing and estimation, which realizes its autonomous navigation in the whole process of flight. It was tested in the actual flight experiment to verify the effectiveness of the proposed method.

2. Flying State Sensing and Estimation Scheme Design and Sensor Calibration

2.1. Mechanical Structure Analysis and Hardware Connection Design

In this paper, the bionic flapping wing robot “Phoenix” independently developed by the laboratory is taken as the research object. The characteristic parameters and flight control hardware connection of the prototype are shown in Figure 1. The design of the transmission mechanism and steering mechanism of the flapping wing robot is shown in Figure 2. The wingspan of the prototype is 2.2 m, and the average chord length is 40 cm. The driving element is a brushless DC (direct current) motor (specification: T-motor AT2308, kv2600), which drives the crank to create circular motion after deceleration by a two-stage gear set (the two-stage gear set is composed of gear 1, gear 2 and gear 3, and gear 1 is fixed on the motor shaft). Then the linkage mechanism is used to control the flutter of the wings up and down. The transmission ratio of the first stage gear (composed of gear 1 and gear 2) is 7, the transmission ratio of the second stage gear (composed of gear 2 and gear 3) is 11. The addendum circle diameter of gear 1 is 7.5 mm, the addendum circle diameter of gear 2 is 48 mm, and the addendum circle diameter of gear 3 is 73 mm. The lift and thrust generated by the periodic flapping of wings are the main source of power for the “Phoenix” flapping robot.
The steering mechanism is composed of two parallelogram mechanisms, which are driven by two steering engines, respectively, and are responsible for steering adjustment of the flapping wing robot. The specification of the pitching steering engine is GDW DS298MG, the torque is 6.5 kg.cm, and the speed is 0.12 s/60°. The specification of the yaw steering engine is Emax ES08MD, the torque is 2.0 kg.cm, and the speed is 0.1 s/60°. The flight control module is powered by the signal line of the electronic speed controller. After completing the tasks of sensor data acquisition and processing, position and attitude solution, and position and attitude control, it outputs three control signals to control the adjustment of motor speed, pitch and roll attitude adjustment, respectively.

2.2. Design of Flying State Sensing and Estimation Scheme

The design of the flying state sensing and estimation function is shown in Figure 3, and the black solid arrows indicate the flow direction of data. First, the RC (remote control) sends control commands, including the selection of flight mode, sensor calibration, etc. After receiving the sensor calibration and filtering command, the flight control module will run the calibration task to calibrate the original data of IMU (inertial measurement unit), barometer and magnetometer. The calibrated data will enter the flying state sensing and estimation task module after filtering.
In the flying state sensing and estimation task module, the attitude of the “Phoenix” flapping wing robot is solved by fusing IMU and magnetometer data, and the position of the “Phoenix” flapping wing robot is solved by fusing GPS, a barometer and inertial navigation. When the position and attitude information of the “Phoenix” flapping wing robot is obtained, the flapping wing robot will run the corresponding guidance algorithm according to its own position information and enter the position and attitude control link. The position control outputs the desired attitude, and then combines the estimated attitude information of the flapping wing robot into the position and attitude control links of autonomous takeoff, cruise and landing flight to output three-way attitude control quantities. It can be seen from Figure 2 that the roll and yaw attitudes of the “Phoenix” flapping wing robot are coupled with each other. The yaw control quantities and the roll control quantities need to be converted into a control quantities output through the preset logic of the mixer, and then used for the control of the roll servo. The pitch control quantities is used to control the pitch servo.

2.3. Sensor Calibration and Filtering

Due to an assembly error and the influence of the surrounding environment, the original data output by the sensor is not accurate enough. Therefore, it is necessary to calibrate and filter the sensor data to obtain relatively accurate sensor data.
(1)
Calibration of the sensor. The calibration of the sensor mainly includes the calibration of the gyroscope, magnetometer and accelerometer. The gyroscope is a device that detects the angular motion of an object, and its error can be calibrated by calculating the average value and Kalman filter by collecting the angular velocity data of three-axis output of gyroscope in static state. The magnetometer is an electronic device that obtains the heading angle by measuring the strength of the geomagnetic field. Here, the ellipsoid fitting based on the least squares method is used to calibrate the magnetometer. The calibration principle of the magnetometer can be found in the literature [12,25]. The effect comparison after calibration is shown in Figure 4. It can be seen that the center of the sphere formed by the raw magnetometer data is not at the origin, which indicates that there is an obvious offset error in the magnetometer. The fitted data are more evenly distributed on the spherical surface and the center of the sphere is located at the origin, which indicates that the ellipsoid fitting method has achieved better results.
For the accelerometer, the ellipsoid fitting method based on nine parameters is used to calibrate it, and the comparison of calibration results is shown in Figure 5. The calibration principle of accelerometer can be found in the literature of Verboom and Hao [11,26]. It can be seen that after the accelerometer is calibrated, the data envelope forms a sphere, and the center of the sphere is closer to the origin, which shows that the ellipsoid fitting method is effective in correcting the systematic error of the accelerometer. Then, keeping the prototype in a horizontal state, the raw data is collected and the calibrated data output by the IMU and the magnetometer when stationary are shown in Table 1. It can be seen that the influence of the sensor system error is basically eliminated by the calibration.
(2)
Sensor filtering. During the flapping process of the flapping wing robot, a lot of vibration noise will be mixed with the flapping of the wings, which will seriously affect the attitude and position estimation, so filtering processing is required.
The comparisons of the accelerometer and gyroscope signals before and after calibration filtering are shown in Figure 6 and Figure 7, respectively. It can be seen that the fluctuation of the horizontal acceleration and the three-axis angular velocity data curve after calibration filtering is significantly reduced, and the signal-to-noise ratio and anti-interference ability are enhanced. At the same time, because the magnetometer signal is not easily affected by mechanical vibration, the sliding mean filter is used to suppress the magnetic field pulse noise.

3. Solution Method of Position and Attitude for Flying State Sensing and Estimation

3.1. Coordinate System Definition and Attitude Kinematics

In order to describe the state of the flapping wing robot in three-dimensional space, including the position of its own center of mass and its attitude in space, the coordinate system fixed with flapping wing robot is defined as the body system {OBXBYBZB}. The origin of the body system is selected at the position of the centroid of the flapping wing robot. The reference coordinate system that characterizes the attitude change of the flapping wing robot is called the inertial frame {OIXIYIZI}, as shown in Figure 8.
The rotation matrix R B I is used to describe the attitude change of the body system {B} relative to the inertial frame {I}, which is called the rotation transformation matrix of the body system {B} relative to the inertial frame {I}. Similarly, the rotation matrix R I B is used to describe the attitude of the inertial frame relative to the body system. Due to the large range of attitude changes of the flapping wing robot, in order to avoid the appearance of attitude jumping, the unit quaternion q = [q0 q1 q2 q3] is used to represent the attitude of the flapping wing robot. The quaternion represents the rotation matrix R B I of the body system relative to the inertial frame as shown in Equation (1)
R B I = [ 2 ( q 0 2 + q 1 2 ) 1 2 ( q 1 q 2 q 0 q 3 ) 2 ( q 1 q 3 + q 0 q 2 ) 2 ( q 1 q 2 + q 0 q 3 ) 2 ( q 0 2 + q 2 2 ) 1 2 ( q 2 q 3 q 0 q 1 ) 2 ( q 1 q 3 q 0 q 2 ) 2 ( q 2 q 3 + q 0 q 1 ) 2 ( q 0 2 + q 3 2 ) 1 ]
Quaternion was used to solve the differential equation of the attitude of flapping wing robot, as shown in Equation (2), and the first-order Runge–Kutta method was also used to update the quaternion differential equation. Then the quaternion expression of the current moment was obtained as shown in Equation (3). Due to the instability of the gyroscope’s bias, the bias error changes gradually over time. Therefore, the key to attitude calculation is to compensate the change of gyroscope bias in real time. In this paper, an explicit complementary filtering algorithm is used to filter the gyroscope data to eliminate the interference of the deviation instability on the attitude solution.
q ^ ˙ = 1 2 q ^ ( ω ^ b + δ )
where q ^ ˙ is the derivative of quaternion, q ^ is the estimated value of the quaternion at any time; w ^ b is the angular velocity vector under the body system; δ is the error compensation term; is the Quaternion multiplication.
q ^ = q ^ init + q ^ ˙ Δ t
where q ^ init is the quaternion at the initial moment; Δ t is the integration time.

3.2. Periodic Equivalent Strategy

On the basis of the coordinate system in Figure 8, the whole-body coordinate system {OsysXsysYsysZsys} is added, as shown in Figure 9. rsys is the position coordinate vector of the center of mass of the whole-body system relative to the inertial reference frame. r0 is the position vector of the center of mass of the fuselage relative to the center of mass of the whole-body system. r1 and r2 are the position vectors of the left and right wings relative to the center of mass of the whole system, respectively. The mass of the fuselage is m0, and the masses of the left and right wings are m1 and m2, respectively.
Because the center of mass of the fuselage fluctuates relative to the center of mass of the whole system during the flapping process of the aircraft, and the sensor is fixedly connected to the fuselage, the height measured by the sensor is not the actual flying height of the aircraft. The relationship between the center of mass of the fuselage and the center of mass of the whole system is:
m 0 r 0 + m 1 r 1 + m 2 r 2 = ( m 0 + m 1 + m 2 ) r sys
The position vector of the body’s center of mass can be simplified as:
r 0 = ( m 0 + m 1 + m 2 ) r sys m 1 r 1 m 2 r 2 m 0 = m 0 + m 1 + m 2 m 0 r s y s m 1 m 0 r 1 m 2 m 0 r 2
Equation (5) can be rewritten as
r 0 = r sys r 1 r 2
where
{ r sys = m 0 + m 1 + m 2 m 0 r sys r 1 = m 1 m 0 r 1 r 2 = m 2 m 0 r 2
Then the fluctuation equation of the center of mass of the fuselage relative to the center of mass of the whole system is:
r 0 r sys = r 1 r 2 = Δ hsin ( 2 π f t )
The fluctuation curve of the center of mass of the fuselage relative to the center of mass of the whole system is obtained, as is shown in Figure 10.
Then the change equation of the center of mass of the fuselage is:
h ( t ) = H ( t ) ± Δ hsin ( 2 π f t )
where H(t) is the height change of the center of mass of the whole system. Δ h is the fluctuation range of the center of mass of the fuselage relative to the center of mass of the whole system. f is the flapping frequency of the wings.
(1)
When the aircraft is in cruise flight
| r sys | = H(t) is a constant value, then the variation curve of the center of mass of the fuselage relative to the center of mass of the whole system is shown in Figure 11.
(2)
When the aircraft is in maneuvering flight
| r sys | = H(t) a is a variable value, then the variation curve of the center of mass of the fuselage relative to the center of mass of the whole system is shown in Figure 12.

3.3. Attitude Solution Based on Explicit Complementary Filtering

Because the traditional complementary filtering algorithm is not suitable for the flapping wing system with severe vibration [9,27], a more widely practical display complementary filtering algorithm was used to solve the attitude of the “Phoenix” flapping wing robot. It combines the data of the accelerometer, magnetometer and gyroscope, and uses the estimation error of the accelerometer and the magnetometer to a certain constant vector to correct the zero-drift change of the three axes of the gyroscope so that the angular velocity data is closer to the real value, thereby suppressing the attitude divergent.

3.3.1. Gravitational Acceleration Compensation

The unit vector of the estimated gravitational acceleration in the inertial frame is shown in Equation (10).
g ¯ I = g I g I = [ 0 0 1 ] T
where g I and g I are the gravitational acceleration vector and modulo length (m/s2) in the inertial frame, respectively. With the help of the rotation matrix (1), the representation of the estimated unit vector of the gravitational acceleration in the body system can be obtained, which is given by the Equation (11).
v ¯ = R I B g ¯ = R B I T g ¯ = [ 2 ( q 1 q 3 q 0 q 2 ) 2 ( q 2 q 3 + q 0 q 1 ) 2 ( q 0 2 + q 3 2 ) 1 ] T
The vector of the specific acceleration measured by the accelerometer under the body system is shown in Equation (12),
f B = a B g B
where Ba is the inertial acceleration vector under the body system; Bg is the gravitational acceleration vector under the body system.
When the flapping robot is in equilibrium, a B 0 ,then the output signal of the accelerometer is the estimation of the gravitational acceleration, that is, g ^ B = f B . After normalization, the Equation (13) will be obtained,
v ^ = g ^ B g ^ B = f B f B
Then the unit vector v ¯ of the estimated gravitational acceleration under the body system is cross-multiplied with the unit vector v ^ of the measured gravitational acceleration to obtain the gravity error correction vector, as shown in Equation (14), which is used to compensate the zero-drift change of the gyroscope.
e 1 = v ^ × v ¯

3.3.2. Magnetic Deflection Compensation

In this paper, the magnetometer is used to correct the yaw angle, and the divergence of the yaw angle is suppressed by the invariance of the magnetic field resultant vector in the horizontal plane. Firstly, the geomagnetic field signal measured under the body system is converted to the inertial frame through the rotation matrix, as shown in Equation (15).
m I = R B I × m B
where Bm is the geomagnetic vector under the body system; Im is the geomagnetic vector in the inertial frame.
The equation for solving the geomagnetic declination in the inertial frame is shown in Equation (16).
ψ ^ 0 = arctan ( m y I / m x I )
The geomagnetic declination of Shenzhen is ψ ^ 0 = 0.052 . Thus, it can be obtained that the difference between the measured geomagnetic declination and the predicted geomagnetic declination is Δ ψ 0 = ψ ^ 0 ψ ¯ 0 , which is rewritten into a vector form, as shown in Equation (17).
e = [ 0 0 Δ ψ 0 ] T
Convert the geomagnetic error vector into the body system, as shown in Equation (18).
e 2 = R I B e
Adding the gravitational acceleration error vector and the geomagnetic error vector, and then the error compensation vector was obtained through the proportional–integral link, as shown in Equation (19). The error compensation vector and the angular velocity vector output by the gyroscope are added to obtain the corrected angular velocity vector, which is then substituted into the quaternion differential equation for integration. The quaternion is then iteratively updated to finally obtain the precise attitude information of the “Phoenix” flapping wing robot.
δ = K p ( e 1 + e 2 ) + K I ( e 1 + e 2 )
where Kp and KI are proportional coefficient and integral coefficient, respectively.
When the “Phoenix” flapping wing robot performs the attitude calculation task, it needs to obtain the accelerometer and magnetometer data to estimate the initial attitude of the “Phoenix” flapping wing robot. In the equilibrium state of the “Phoenix” flapping wing robot, the three-axis accelerometer data is collected to calculate the initial roll angle and pitch angle of the flapping wing robot, as shown in Equations (20) and (21).
φ = arctan ( a y b / a x b )
θ = arcsin ( a x b )
where φ is the roll angle of flapping wing robot, and θ is the pitch angle.
The geomagnetic vector obtained by the magnetometer under the body system is projected into the horizontal plane of the inertial frame, as shown in Equations (22) and (23).
m x = m x b cos θ + m y b sin θ sin φ m z b sin θ cos φ
m y = m y b cos φ + m z b sin φ
Initialize the yaw angle (looking from the tail to the nose, the right is positive) using the magnetometer data in the north (X-axis) and east (Y-axis) directions, as shown in Equation (24).
ψ = atan 2 ( - m y , m x ) mag _ decline
where mag_decline is the magnetic declination angle, which is about −0.052 rad in Shenzhen.
Then, the initial Euler angle is converted into a quaternion by Equation (25), and then the quaternion vector at the initial moment can be obtained by normalizing.
{ q 0 = cos ( ψ / 2 ) cos ( θ / 2 ) cos ( φ / 2 ) + sin ( ψ / 2 ) sin ( θ / 2 ) sin ( φ / 2 ) q 1 = cos ( ψ / 2 ) cos ( θ / 2 ) sin ( φ / 2 ) sin ( ψ / 2 ) sin ( θ / 2 ) cos ( φ / 2 ) q 2 = cos ( ψ / 2 ) sin ( θ / 2 ) cos ( φ / 2 ) + sin ( ψ / 2 ) cos ( θ / 2 ) sin ( φ / 2 ) q 3 = cos ( ψ / 2 ) sin ( θ / 2 ) sin ( φ / 2 ) + sin ( ψ / 2 ) cos ( θ / 2 ) cos ( φ / 2 )
After the quaternion initialization is completed, the offset of the three-axis gyroscope is compensated by displaying the complementary filtering algorithm to obtain the error compensation vector δ , and then the quaternion is updated by using Equations (2) and (3). If the solved quaternion is a finite real number, the solution is considered successful, otherwise the solution result is discarded. After the quaternion solution is completed, the normalization processing is carried out, and the corresponding rotation matrix R B I is updated. In this way, the attitude solution task of one cycle (0.01 s) is completed.

3.4. Position Solution Based on Kalman Filter

The position solution of the “Phoenix” flapping wing robot in the horizontal direction requires the use of sensor data. The GPS sensor output data has long-term stability and does not diverge over time, but the data refresh rate is low and is easily affected by surrounding obstructions. The INS (inertial navigation system) can be used in any complex environment, with high accuracy in a short time, but the position measured by INS will produce integral drift for a long time. A barometer and GPS can be used to estimate altitude information. In a windless and stable flight environment, the barometer is used to perceive the altitude information of the “Phoenix” flapping wing robot. When the flapping robot is disturbed by strong winds or its own flapping frequency is too fast, it will cause large fluctuations in the surrounding air pressure, resulting in inaccurate data output of the barometer [20,23]. When the output data measured by the barometer fluctuates unreasonably, the altitude data output by GPS is selected as the height measurement information of the flapping robot.
In summary, it is difficult for each navigation system to be used independently to meet the reliability requirements. A feasible way to improve the accuracy of the navigation system is the integrated navigation technology, that is, two or more navigation systems are used to measure and solve the same navigation information to form a more accurate measurement value.

3.4.1. Design of Kalman Filter

In this study, the Kalman filter is used to fuse multi-sensor data to solve the position and velocity information of the “Phoenix” flapping wing robot [13,28]. The current state can be calculated only by inputting the measured value of the state vector in the current solution cycle (which may contain multiple sensor data) and the estimated value of the previous solution cycle.
Setting x = [p v b]T as the northward state vector of the flapping wing robot, where p represents the position; v represents the velocity; b represents the inertial acceleration bias. The solution using the Kalman filter includes prediction process and update process. The prediction process is shown in Equations (26) and (27).
x ^ k = F x ^ k 1 + B a k + Γ ζ k 1
where F is the state transition matrix; B is the control matrix; Γ is the process noise driving matrix. x ^ k and x ^ k 1 are the prior state vector at time k and the posterior state vector at time (k–1), respectively; ak is the inertial acceleration in the inertial frame at time k; ζ k 1 is the zero mean white noise.
P ^ k = F P ^ k 1 F T + Q
where P ^ k is the prior estimated state error covariance matrix at time k; P ^ k 1 is the posterior estimated state error covariance matrix at time k–1; Q is the process excitation noise covariance matrix.
The update process of the Kalman filter is shown in Equations (28)–(31). Equations (26) and (28) are collectively referred to as the Kalman filter system equation of the “Phoenix” flapping wing robot.
z k = H x ^ k + v k
x ^ k = x ^ k + K ( z k H x ^ k )
K = P ^ k H T H P k H T + R
P ^ k = ( I K H ) P ^ k
where z k is measurement vector, z k = [ p k m   v k m ] T ; p k m and v k m are the position and velocity measured by the sensor, respectively; vk is the measurement error vector, v k = [ p e v e ] T ; p e and v e represent the measurement error of position and velocity, respectively. I is the three order unit matrix; K is the Kalman gain matrix; H is the measurement matrix; R is the measurement noise covariance matrix.
The design of the Kalman filter focuses on three covariance matrices, namely the initial state error covariance matrix P ^ 0 , the process excitation noise covariance matrix Q and the measurement noise covariance matrix R. Here, the initial state error covariance matrix P ^ 0 is a third-order square matrix, the elements on the diagonal are the variance of the initial estimation error of the state vector, and the elements on the non-diagonal line are the covariance of the initial estimation error of the state vector.
When setting the process excitation noise covariance matrix, it can be seen from the system model that the errors of position and speed are related to the change of acceleration, and the relationship is shown in Equations (32) and (33). Assuming that the variance of the acceleration is D(a), the expressions of the position and velocity variances are shown in (34) and (35) according to the properties of the variance, thus determining the elements on the diagonal of the Q matrix. The calculation of covariance is given by Equation (36), and the correlation coefficient is taken as ρ = ± 1 . When the direction of position and velocity is the same as the direction of acceleration, the correlation coefficient ρ = 1 , otherwise ρ = 1 . In summary, under the condition that the state quantity is positively correlated with the driving quantity, the final expression of the process excitation noise Q matrix is shown in Equation (37).
Δ v = a Δ t
Δ p = v 0 Δ t + 1 2 a Δ t 2
D ( v ) = D ( a t ) = Δ t 2 D ( a )
D ( p ) = D ( 1 2 a Δ t 2 ) = 1 4 Δ t 4 D ( a )
cov ( X , Y ) = ρ D ( X ) D ( Y )
where D(X) and D(Y) represent the variance of X variable and Y variable, respectively; ρ represents the correlation coefficient, 1 ρ 1 ;
Q = [ 1 4 Δ t 4 1 2 Δ t 3 1 2 Δ t 2 1 2 Δ t 3 Δ t 2 Δ t 1 2 Δ t 2 Δ t 1 ] D ( a )
The measurement noise covariance is mainly related to the measurement accuracy of the sensor, and the measurement error can be determined according to the horizontal (vertical) positioning accuracy factor output by the GPS sensor. It is also assumed that the GPS measurements of position and velocity are independent of each other, then the non-diagonal elements of the measurement noise covariance matrix R are set to 0, and the diagonal elements are related to the positioning accuracy factor. The final expression of the measurement noise covariance moment RH in the horizontal direction is shown in Equation (38).
R H = [ D ( h d o p ) 0 0 D ( v d o p ) ]
where D(hdop) and D(vdop) represent the variance of the horizontal position and horizontal velocity positioning accuracy factor, respectively. The measurement noise covariance Rv in the vertical direction is related to the vertical positioning factor of the barometer or GPS. When the barometer is used to measure the height information of the “Phoenix” flapping wing robot, the expression of Rv is shown in Equation (39).
R V = D ( b a r o )
where D(baro) is the variance of the height measured by the barometer.
When the error of the barometer output data is large, GPS is used to sense the height change of the flapping wing robot, and the expression of Rv is shown in Equation (40).
R V = D ( v d o p )
where D(vdop) is the variance of the vertical positioning precision factor of GPS.

3.4.2. Improvement of the Kalman Filter Algorithm

(1)
H∞ filter
If the prior information of the covariance matrices Q, R, and P0 in the Kalman filter equation is uncertain, when the flapping robot is disturbed by gusts, it will easily lead to inaccuracy of the system model and unclear statistical characteristics of noise. At this time, it is necessary to appropriately increase the value on the diagonal of the process noise covariance matrix Q to increase the utilization weight of the real-time measured value and reduce the utilization weight of the predicted value. This process has great randomness, and it is impossible to determine the appropriate value of the process noise covariance matrix Q to optimize the performance of the filter. H∞ filtering is an effective method to solve the above problems [19,22], which is of great help to improve the robustness of the system. The discrete H∞ filter system equations and algorithms are shown in Equations (41) and (42), respectively.
x ^ k + 1 = F x ^ k + B a k + w k z k = H x ^ k + v k y k = L k x ^ k
where yk is the estimated vector, which can be a linear combination of state variables, and it is necessary to ensure that Lk is full rank; wk, vk are system noise and measurement noise, respectively.
x ^ k + 1 = F x ^ k + F K ( z k H x ^ k ) K = P k ( I + H T R 1 H P k ) 1 H R 1 P k + 1 = F P k ( I + H T R 1 H P k ) 1 F T + Q
After many experiments, the values of each covariance matrix are as follows. Among them, initial state error covariance matrix and process noise matrix in the X direction:
P ^ 0 x = [ 8 0 0 0 0.18 0 0 0 6.4 × 10 6 ] , Q x = [ 1.15 0 0 0 0.8 0 0 0 1 × 10 6 ]
Initial state error covariance matrix and process noise matrix in Y direction:
P ^ 0 y = [ 32 0 0 0 0.5 0 0 0 6.4 × 10 6 ] , Q y = [ 5.8 0 0 0 5 0 0 0 1 × 10 6 ]
Horizontal measurement noise matrix:
R H = [ 1 0 0 0.04 ]
Initial state error covariance matrix and process noise matrix in Z direction:
P ^ 0 z = [ 12.5 0 0 0 3.38 0 0 0 1 × 10 4 ] , Q z = [ 3.5 0 0 0 2.0 0 0 0 1 × 10 6 ]
When using barometer data to represent height measurements, the measurement noise in the Z direction is:
R V = 1.524
When using the GPS output data to represent the height measurement value, the measurement noise in the Z direction is:
R V = 9
(2)
Suppression of filter divergence
In the actual flight experiment, it was found that with the prolongation of flight time, the navigation parameters output by the Kalman filter tended to be stable, but the deviation from the real value became larger and larger. This is because the torsion, bending and other deformations will occur during the flapping process of the wings, which makes it difficult to establish an accurate mechanism model to analyze the force of the flapping robot during flight. Therefore, the physical model error is the main factor leading to the divergence of the filter. The basic idea of model error correction is to change the value of the elements in the covariance matrix to affect the size of the elements of the Kalman gain matrix, and then reasonably assign the weights of the predicted value and measured value of the state quantity. Here, the attenuation memory method is used to suppress the filter divergence.
The attenuation memory method filtering equations are shown in Equations (43)–(45). Compared with Equation (27), the filtering equation only has one more scalar factor λ in Equation (43). Since λ > 1 , the prior state error covariance matrix P ˜ k is always larger than P k , so there is always a Kalman gain matrix K > K . This shows that the use of the attenuation filter algorithm for the new measurement value is more weighted than the Kalman filter, which will cause the filter gain matrix to converge to a larger value, thus reducing the impact of model error. The larger the value of the scalar factor λ is, the higher the trust in the latest measurement value and the lower the trust in the historical measurement value. Through actual debugging, it is determined that λ = 1.02 .
P ˜ k = λ 2 F P ˜ k 1 F T + Q
where λ is the scalar factor, λ > 1 .
K = P ˜ k H T ( H P ˜ k H T + R ) 1
P ˜ k = P ˜ k K H P ˜ k

4. Flight Test Verification

4.1. Verification Experiment of Attitude Solution

In order to verify the actual effect of the complementary filtering algorithm, the flight control module is installed at the center of mass of the “Phoenix” flapping wing robot, as shown in Figure 13. The flight control module also integrates the industrial-grade attitude calculation module Xsens-MTi3 and the MPU6050 chip at the same time. Among them, the MPU6050 is a traditional motion tracking chip, which is packaged with a three-axis gyroscope, a three-axis accelerometer and a DMP (Motion Digital Processor, DMP). The performance comparison between the chip and the Xsens-MTi3 module is shown in Table 2. It can be seen that the Xsens-MTi3 module has higher attitude estimation capability and attitude calculation accuracy, but it is expensive. Therefore, from the perspective of cost, MPU6050 is more suitable for attitude estimation of the “Phoenix” flapping wing robot.
After the flight control module is installed, the flapping wing robot is controlled by the remote controller to fly for a period of time. The attitude information solved by the display complementary filtering algorithm and the attitude data output by Xsens-MTi3 are recorded by the SD (Secure Digital, SD) card. The comparison of attitude changes over time is shown in Figure 14.
It can be seen that the initial attitude solution by complementary filtering gradually converges with time. During the flight process of the “Phoenix” flapping wing robot, the results of the display complementary filtering algorithm and the output of Xsens-MTi3 are basically consistent, and the attitude estimation is always convergent during the whole flight process. The results show that the attitude calculation accuracy is less than 5 degrees, which basically meets the stability and rapidity indicators of the system design.

4.2. Verification Experiment of Position Solution

The flight test was conducted in an open sports field with no obstructions and strong magnetic field interference. The SD card was used to record the flight trajectory of the “Phoenix” flapping wing robot, and the accuracy and stability of the position calculation algorithm were evaluated. First of all, in order to verify the accuracy of the position calculation, let the “Phoenix” flapping robot flap at a certain frequency, making a circle around the middle of the sports field (the length of the middle ring of the sports field is about 158 m and the width is about 78 m). Then, walk around the outside of the sports field for half a circle, and the obtained plane trajectory is shown in Figure 15 (left). From the figure, the northward distance between points 1 and 2 is 158.13 m, and the eastward distance between points 3 and 4 is 78.98 m. It can be seen that the estimation accuracy of the horizontal position is less than 1.5 m, which meets the relevant requirements. The variation of height with time is shown in Figure 15 (right). Taking the height at the beginning of the detour as the benchmark, it can be seen that the fluctuation of the vertical position is within 1 m, which also meets the requirements.
Then, the “Phoenix” flapping wing robot is controlled by the remote control to fly over the playground for a period of time, and the position and speed information during the flight are recorded. The variation of speed with time during flight is shown in Figure 16. It can be seen that the maximum horizontal resultant velocity speed of the “Phoenix” flapping wing robot is about 10 m/s, and the average horizontal resultant velocity speed is 7 m/s, which has strong maneuverability. In the vertical direction, due to the reciprocating vibration of the fuselage, the speed changes very drastically. The maximum climbing speed is about 2.5 m/s, and the maximum descending speed is about 5 m/s.
Based on the above two experiments, it can be seen that in the whole flight process, the Kalman filter algorithm does not have any divergence in the estimation of the horizontal position, and the stability is strong. The estimation of the flying height of the flapping wing robot tends to be stable on the whole. the Kalman filter algorithm is used to solve the position and speed of the “Phoenix” flapping wing robot, which meets the requirements of accuracy and stability, and has received better results.

5. Conclusions

In view of the large mass ratio of the wings of the large bionic flapping wing flying robot and the fluctuation of the fuselage centroid during the flight, the traditional state perception and estimation methods struggle to provide stable and accurate data. This paper takes the “Phoenix” large-scale bionic flapping wing flying robot with a wingspan of 2.2 m as the research object, studies its integrated navigation method of attitude and position, and develops a corresponding flight controller to carry out flight experiment verification. A three-axis attitude estimation method for flapping wing flying robot based on the fusion of inertial sensor and magnetometer data is proposed, which solves the problem of attitude divergence caused by periodic heaving and pitching motions. The developed explicit complementary filtering algorithm fuses the IMU and magnetometer data to solve the attitude information of the flapping wing robot. The results obtained by the explicit complementary filtering algorithm are compared with the high-precision attitude module Xsens-MTi3. The results show that the attitude calculation accuracy is less than 5 degrees, which indicates that the use of low-cost sensors has achieved high-precision position and attitude estimation.
The position and velocity of Phoenix flapping wing robot in three-dimensional space are solved by using integrated navigation technology and the Kalman filter to fuse inertial navigation, GPS and barometer data. The Kalman filter is optimized by using the filtering algorithm, which improves the robustness of the Kalman filter. Aiming at the problem of filter divergence, the attenuation memory method is used to correct the model error, which improves the accuracy of position estimation. Finally, a complete position calculation process is planned and flight experiments are carried out. The results show that the accuracy and stability of the position solution meet the requirements when the sensor signal is good. At present, this method is mainly studied from the perspective of multi rigid bodies, and satisfactory results are obtained. In the future, the flexible vibration characteristics and fluid structure coupling characteristics of the wing will be considered, and on this basis, the method will be further optimized to improve the accuracy of state perception and estimation, ensuring real-time performance.

Author Contributions

G.L. contributed the central idea and wrote the initial draft of the paper. S.W., and W.X. contributed to refining the ideas, and revised this paper. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the Shenzhen excellent scientific and technological innovation talents training project (basic research on outstanding young people). Research on flight dynamics and cooperative control of large Bionic Flapping Wing robot formation, No. RCJC2020071414436040.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ajanic, E.; Feroskhan, M.; Mintchev, S.; Noca, F.; Floreano, D. Bioinspired wing and tail morphing extends drone flight capabilities. Sci. Robot. 2020, 5, eabc2897. [Google Scholar] [CrossRef] [PubMed]
  2. Lau, G.A. Stunt flying hawk-inspired drone. Sci. Robot. 2020, 5, eabe8379. [Google Scholar] [CrossRef] [PubMed]
  3. Wang, S.P.; Du, C.P.; Zheng, Y. Path planning algorithm of flapping wing aircraft based on reinforcement learning. J. Control Decis. 2022, 37, 851–860. [Google Scholar]
  4. Wang, Y.P. Research on Autonomous Formation Flight Control Method of Large Bionic Flapping Wing Flying Robot; Harbin University of Technology: Harbin, China, 2021; pp. 12–23. [Google Scholar]
  5. Li, G.; Chen, X.; Zhou, F.; Liang, Y.; Xiao, Y.; Cao, X.; Zhang, Z.; Zhang, M.; Wu, B.; Yin, S.; et al. Self-powered soft robot in the mariana trench. Nature 2021, 591, 66–71. [Google Scholar] [CrossRef] [PubMed]
  6. Di Luca, M.; Mintchev, S.; Su, Y.; Shaw, E.; Breuer, K. A bioinspired separated flow wing provides turbulence resilience and aerodynamic efficiency for miniature drones. Sci. Robot. 2020, 5, eaay8533. [Google Scholar] [CrossRef]
  7. Xuan, J.L.; Song, B.F.; Song, W.P.; Yang, W.Q.; Xue, D.; Liang, S.R. Progress of chinese dove and future studies on flight mechanism of birds and application system. Trans. Nanjing Univ. Aeronaut. Astronaut. 2020, 37, 663–675. [Google Scholar]
  8. Ramezani, A.; Chung, S.J.; Hutchinson, S. A biomimetic robotic platform to study flight specializations of bats. Sci. Robot. 2017, 2, eaal2505. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Kamankesh, Z.; Banazadeh, A. Stability analysis for design improvement of bio-inspired flapping wings by energy method. Aerosp. Sci. Technol. 2021, 111, 106558. [Google Scholar] [CrossRef]
  10. Yu, B.; Cheng, Y.Y.; Jin, X.Z.; Du, H.B. Robust finite time attitude stabilization of rigid spacecraft with angular velocity constraints. Control. Decis. 2021; in press. [Google Scholar]
  11. Verboom, J.L.; Tijmons, S.; Wagter, C.D.; Remes, B.; Babuska, R.; de Croon, G.C.H.E. Attitude and altitude estimation and control on board a flapping wing micro air vehicle. In Proceedings of the 2015 IEEE International Conference on Robotics and Automation (ICRA), Seattle, WA, USA, 26–30 May 2015; pp. 5846–5851. [Google Scholar]
  12. Zhu, X.; Zhao, T.; Cheng, D.; Zhou, Z. A three-step calibration method for tri-axial field sensors in a 3D magnetic digital compass. Meas. Sci. Technol. 2017, 28, 055106. [Google Scholar] [CrossRef]
  13. Yang, W.; Wang, L.; Song, B. Dove: A biomimetic flapping-wing micro air vehicle. Int. J. Micro Air Veh. 2018, 10, 70–84. [Google Scholar] [CrossRef] [Green Version]
  14. Tu, Z.; Yang, Y.; Fei, F.; Zhang, J.; Deng, X. Realtime on-board attitude estimation of high-frequency flapping wing mavs under large instantaneous oscillation. In Proceedings of the 2018 IEEE International Conference on Robotics and Automation (ICRA), Brisbane, QLD, Australia, 21–25 May 2018; pp. 6806–6811. [Google Scholar]
  15. He, W.; Mu, X.; Zhang, L.; Zou, Y. Modeling and trajectory tracking control for flapping-wing micro aerial vehicles. IEEE/CAA J. Autom. Sin. 2021, 8, 148–156. [Google Scholar] [CrossRef]
  16. He, W.; Ding, S.Q.; Sun, C.Y. Research progress on modeling and control of flapping wing aircraft. Acta Autom. Sin. 2017, 43, 685–696. [Google Scholar]
  17. Taha, H.E.; Kiani, M.; Hedrick, T.L.; Greeter, J.S. Vibrational control: A hidden stabilization mechanism in insect flight. Sci. Robot. 2020, 5, eabb1502. [Google Scholar] [CrossRef] [PubMed]
  18. Bialy, B.J.; Chakraborty, I.; Cekic, S.C.; Dixon, W.E. Adaptive boundary control of store induced oscillations in a flexible aircraft wing. Automatica 2016, 70, 230–238. [Google Scholar] [CrossRef]
  19. Dou, L.Q.; Mao, Q.; Su, P.H. Reentry attitude control of reusable launch vehicle based on adaptive fuzzy h∞ control. J. Control Decis. 2018, 33, 1181–1189. [Google Scholar]
  20. Zhang, J.; Cheng, B.; Yao, B.; Deng, X. Adaptive robust wing trajectory control and force generation of flapping wing MAV. In Proceedings of the 2015 IEEE International Conference on Robotics and Automation (ICRA), Seattle, WA, USA, 26–30 May 2015; pp. 5825–5857. [Google Scholar]
  21. Jafferis, N.T.; Helbling, E.F.; Karpelson, M.; Wood, R.J. Untethered flight of an Insect-sized flapping-wing microscale aerial vehicle. Nature 2019, 570, 491–495. [Google Scholar] [CrossRef] [PubMed]
  22. Karasek, M.; Percin, M.; Cunis, T.; van Oudheusden, B.W.; De Wagter, C.; Remes, B.D.; de Croon, G.C. Accurate position control of a flapping-wing robot enabling free-flight flow visualisation in a wind tunnel. Int. J. Micro Air Veh. 2019, 11, 1756829319833683. [Google Scholar] [CrossRef] [Green Version]
  23. Bailing, T.; Pinpin, L.; Hanchen, L.; Qun, Z. Trajectory and attitude cooperative control of multiple uavs in complex environment. Acta Aeronaut. Et Astronaut. Sin. 2020, 41, 36–43. [Google Scholar]
  24. Wang, Z.W.; Wang, F.J.; Di, C.C.; Shi, Z.; Yang, G. Nonlinear alignment method of land inertial navigation system. Acta Aeronaut. Et Astronaut. Sin. 2018, 39, 283–293. [Google Scholar]
  25. Wu, H.; Pei, X.; Li, J.; Gao, H.; Bai, Y. An improved magnetometer calibration and compensation method based on Levenberg–Marquardt algorithm for multi-rotor unmanned aerial vehicle. Meas. Control. 2020, 53, 276–286. [Google Scholar] [CrossRef]
  26. Zhao, G.L.; Xu, D.F.; Ge, L.Z.; Guo, Q.R. Research on calibration method of MEMS accelerometer based on kalman filter. Transducer Microsyst. Technol. 2021, 40, 25–27. [Google Scholar]
  27. Zhang, Z.; Wang, J.B.; Song, B.; Tong, G.-f. Adaptive complementary filtering algorithm for imu based on mems. In Proceedings of the 2020 Chinese Control And Decision Conference (CCDC), Hefei, China, 22–24 August 2020; pp. 5409–5416. [Google Scholar]
  28. Khamseh, H.B.; Ghorbani, S.; Janabi, F. Unscented kalman filter state estimation for manipulating unmanned aerial vehicles. Aerosp. Sci. Technol. 2019, 92, 446–463. [Google Scholar] [CrossRef]
Figure 1. The characteristic parameters and flight control hardware connection of the ”Phoenix” flapping wing robot. (a) Flapping wing robot and its characteristic parameters. (b) flight control hardware connection.
Figure 1. The characteristic parameters and flight control hardware connection of the ”Phoenix” flapping wing robot. (a) Flapping wing robot and its characteristic parameters. (b) flight control hardware connection.
Actuators 11 00213 g001
Figure 2. Transmission mechanism and steering mechanism of flapping wing robot. (a) Transmission mechanism. (b) steering mechanism.
Figure 2. Transmission mechanism and steering mechanism of flapping wing robot. (a) Transmission mechanism. (b) steering mechanism.
Actuators 11 00213 g002
Figure 3. Functional block diagram of flying state sensing and estimation.
Figure 3. Functional block diagram of flying state sensing and estimation.
Actuators 11 00213 g003
Figure 4. Ellipsoid fitting effect of magnetometer.
Figure 4. Ellipsoid fitting effect of magnetometer.
Actuators 11 00213 g004
Figure 5. Ellipsoid fitting effect of accelerometer.
Figure 5. Ellipsoid fitting effect of accelerometer.
Actuators 11 00213 g005
Figure 6. Comparison diagram of accelerometer signal before and after calibration and filtering. (a) Specific acceleration in X-axis. (b) Specific acceleration in Y-axis. (c) Specific acceleration in Z-axis.
Figure 6. Comparison diagram of accelerometer signal before and after calibration and filtering. (a) Specific acceleration in X-axis. (b) Specific acceleration in Y-axis. (c) Specific acceleration in Z-axis.
Actuators 11 00213 g006
Figure 7. Comparison diagram of gyro signal before and after calibration and filtering. (a) Angular velocity of X-axis. (b) Angular velocity of Y-axis. (c) Angular velocity of Z-axis.
Figure 7. Comparison diagram of gyro signal before and after calibration and filtering. (a) Angular velocity of X-axis. (b) Angular velocity of Y-axis. (c) Angular velocity of Z-axis.
Actuators 11 00213 g007
Figure 8. Schematic diagram of inertial frame and body system orientation. (a) Schematic diagram of inertial system. (b) Schematic diagram of body system.
Figure 8. Schematic diagram of inertial frame and body system orientation. (a) Schematic diagram of inertial system. (b) Schematic diagram of body system.
Actuators 11 00213 g008
Figure 9. Vector position diagram of aircraft centroid.
Figure 9. Vector position diagram of aircraft centroid.
Actuators 11 00213 g009
Figure 10. Height fluctuation curve of fuselage centroid relative to the centroid of the whole machine system.
Figure 10. Height fluctuation curve of fuselage centroid relative to the centroid of the whole machine system.
Actuators 11 00213 g010
Figure 11. Variation curve of centroid height and periodic equivalent height during cruise flight.
Figure 11. Variation curve of centroid height and periodic equivalent height during cruise flight.
Actuators 11 00213 g011
Figure 12. Centroid height variation curve and periodic equivalent height during maneuvering flight.
Figure 12. Centroid height variation curve and periodic equivalent height during maneuvering flight.
Actuators 11 00213 g012
Figure 13. Installation diagram of flight control module.
Figure 13. Installation diagram of flight control module.
Actuators 11 00213 g013
Figure 14. Comparison diagram of attitude solution. (a) Comparison of pitch angle changes. (b) Comparison of rolling angle changes.
Figure 14. Comparison diagram of attitude solution. (a) Comparison of pitch angle changes. (b) Comparison of rolling angle changes.
Actuators 11 00213 g014
Figure 15. Plane track of detour playground and variation of height. (a) Plane track of detour playground. (b) Variation of height.
Figure 15. Plane track of detour playground and variation of height. (a) Plane track of detour playground. (b) Variation of height.
Actuators 11 00213 g015
Figure 16. Flight speed variation diagram. (a) Speed variation in X direction. (b) Speed variation in Y direction. (c) Speed variation of horizontal resultant speed. (d) Speed variation in Z direction.
Figure 16. Flight speed variation diagram. (a) Speed variation in X direction. (b) Speed variation in Y direction. (c) Speed variation of horizontal resultant speed. (d) Speed variation in Z direction.
Actuators 11 00213 g016
Table 1. Data comparison table before and after sensor calibration.
Table 1. Data comparison table before and after sensor calibration.
Static MeasurementRaw DataCalibrated Data
Accelerometer
(m/s2)
X0.420.15
Y0.04−0.20
Z10.19−9.69
Module10.219.70
Gyroscope
(°/s)
X2.38−0.04
Y1.700.005
Z1.340.05
Module3.170.06
MagnetometerX0.150.16
Y0.490.47
Z0.140.14
Module0.530.51
Table 2. Comparison between MPU6050 and Xsens-MTi3.
Table 2. Comparison between MPU6050 and Xsens-MTi3.
MPU6050Xsens-MTi3
size4 × 4 × 0.9 mm12 × 12 × 2 mm
weight<1 g8 g
Calculating power200 Hz1 kHz
accuracyØ 10°<1°
cost8~12¥1000~9000¥
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, G.; Wang, S.; Xu, W. Flying State Sensing and Estimation Method of Large-Scale Bionic Flapping Wing Flying Robot. Actuators 2022, 11, 213. https://doi.org/10.3390/act11080213

AMA Style

Liu G, Wang S, Xu W. Flying State Sensing and Estimation Method of Large-Scale Bionic Flapping Wing Flying Robot. Actuators. 2022; 11(8):213. https://doi.org/10.3390/act11080213

Chicago/Turabian Style

Liu, Guangze, Song Wang, and Wenfu Xu. 2022. "Flying State Sensing and Estimation Method of Large-Scale Bionic Flapping Wing Flying Robot" Actuators 11, no. 8: 213. https://doi.org/10.3390/act11080213

APA Style

Liu, G., Wang, S., & Xu, W. (2022). Flying State Sensing and Estimation Method of Large-Scale Bionic Flapping Wing Flying Robot. Actuators, 11(8), 213. https://doi.org/10.3390/act11080213

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