Next Article in Journal
Corrosion Monitoring of Reinforced Steel Embedded in Cement Mortar under Wet-And-Dry Cycles by Electrochemical Impedance Spectroscopy
Previous Article in Journal
A Hydrodynamic Model for Measuring Fluid Density and Viscosity by Using Quartz Tuning Forks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Function Extension Based Real-Time Wavelet De-Noising Method for Projectile Attitude Measurement

1
School of Automation, Beijing Institute of Technology, Beijing 100081, China
2
School of Automation, Beijing Information Science and Technology University, Beijing 100192, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(1), 200; https://doi.org/10.3390/s20010200
Submission received: 26 November 2019 / Revised: 24 December 2019 / Accepted: 26 December 2019 / Published: 30 December 2019
(This article belongs to the Section Physical Sensors)

Abstract

:
The real-time measurement of the projectile attitude is the key to realize the whole process guidance of the projectile. Due to the high dynamic characteristics of the projectile motion, the attitude measurement is affected by the real-time and accuracy of the gyro signal de-noising. For the nonlinear discontinuity of the conventional extension method in real-time wavelet de-noising, a function extension real-time wavelet de-noising method is proposed in this paper. In this method, a prediction model of gyro signal is established based on the Roggla formula. According to the model, the signal is fitted in the sliding window, and the signal of the same length is predicted to realize the right boundary extension. The simulation and experiment results show that compared with the traditional extension method, the proposed method can in-crease the signal-to-noise ratio (SNR) and the smoothness, and can decrease the attitude mean absolute error (AMAE) and the attitude root mean square error (ARMSE). Moreover, the time delay caused by signal de-noising can be effectively solved. The real-time performance of the attitude measurement can be ensured.

1. Introduction

Gyro is the core component of the projectile attitude measurement. Its accuracy and real-time performance affect the flight and control safety of the projectile. Due to the high dynamic characteristics of projectile in flight, such as high overload and high speed [1], there are many uncertain factors, such as the complex working environment, which cause interference to the measurement of gyro and affect its measurement accuracy [2,3]. Figure 1 shows a group of actual measurement data of gyro. It can be seen that the dynamic characteristics of the projectile flight environment are strong, the data fluctuate greatly, and the introduced nonlinear noise has a great influence on the attitude measurement. Therefore, in the high dynamic environment, it is the key to remove the gyro noise effectively for accurate measurement of the projectile attitude.
For the gyro random noise processing, the main methods at present include: low-pass, high-pass, band-pass or band-stop filter algorithm [4,5], Kalman filter and its derivative algorithm [6]. The wavelet de-noising algorithm overcomes the defects of traditional filter methods, and is with good local characteristics, adaptive characteristics and multi-resolution characteristics [7]. Therefore, wavelet algorithm is widely used in signal de-noising processing.
At present, the wavelet de-noising method mainly focuses on offline batch processing of measurement data [8]. In practical engineering, there is boundary distortion phenomenon in real-time wavelet de-noising of gyro signal [9,10,11], and boundary extension processing is needed for data window. Jianing S proposed the periodic extension technique to eliminate the boundary issue inherent in kernel regression means [12]. Yang A proposed a construction method of wavelet transform matrices with arbitrary decomposition depth while the signal was symmetrically extended [13]. Wei J proposed a novel method for reducing these tile-boundary artifacts. The method employed odd tile sizes (2N + 1 samples) rather than the conventional even tile sizes (2N samples) [14]. Qian et al. proposed neural network prediction model and wavelet packet transform technology into the empirical mode decomposition method to improve the border effect and to enhance the ability of signal de-noising [15]. Common extension methods mainly include: zero extension, periodic extension, symmetric extension and linear extension [8,9,16]. These methods can achieve better results in extension processing of signals with low dynamic performance, but in the high dynamic environment, the output signal of gyro has strong nonlinear characteristics, so the results of these methods are not ideal. Therefore, it is urgent to find a wavelet de-noising boundary extension method for high dynamic application environment.
In this paper, the Roggla formula is introduced to establish the gyro signal prediction model of the projectile attitude measurement. According to this model, a boundary extension method based on function fitting is proposed. The least square method is used to fit the data of the current sliding window and the data of the next sliding window is predicted. A new sliding window with the fitting and prediction results is formed, and the wavelet threshold de-noising is realized. The signal-to-noise ratio (SNR), root mean square error (RMSE), smoothness, attitude mean absolute error (AMAR) and attitude root mean square error (ARMSE) are used as the evaluation indexes of the algorithm. Several common extension methods and the function extension method proposed in this paper are compared and analyzed through simulation and experiment. The results show that the effects of the five evaluation indexes are improved by the function extension method, and the processing results are stable to follow the real signal, and the problem of time delay caused by signal de-noising is solved.
The rest of the paper is organized as follows. In Section 2, we introduce the principle of real-time wavelet threshold de-noising. In Section 3, we establish the prediction model of gyro signal, and propose the function extension method. In Section 4, we analyze the simulation and experiment results. Section 5 draws the conclusions.

2. Principle of Real-Time Wavelet Threshold De-Noising

The real-time wavelet threshold de-noising process is shown in Figure 2 [17].
In the process of real-time wavelet threshold de-noising, the current signal is located at the right end of the sliding window, and the optimal position of de-noising is located in the middle of the sliding window. Therefore, in this paper only the right boundary extension is performed to shorten the sliding window length and improve the real-time performance of the algorithm. Common boundary extension methods include:
(1)
Zero extension. As shown in Figure 3a, this method is easy to produce large step change at the boundary, and introduce large errors.
(2)
Periodic extension. As shown in Figure 3b, this method is not applicable when the data changes sharply, which is easy to cause boundary discontinuity.
(3)
Symmetric extension. As shown in Figure 3c, this method is not applicable when the data changes violently.
(4)
Linear extension. As shown in Figure 3d, this method is not applicable for nonlinear systems.
In the high dynamic environment, the gyro output data changes violently. In order to avoid the problems introduced by the above method, in this paper a boundary extension method is proposed based on function fitting.

3. Real-Time Wavelet De-Noising Method Based on Function Extension

3.1. Prediction Model of Gyro Signal for Projectile Attitude Measurement

The yaw angle ψ , pitch angle θ and roll angle γ of the projectile attitude are a set of Euler angles. The angular position relationship of the body coordinate system ( O x b y b z b ) relative to the navigation coordinate system ( O x n y n z n ) are described by the Euler angles. The definition of coordinate system is shown in Figure 4.
Since the latitude and longitude of the projectile from launching point to the landing point do not change much, it can be considered that the gyro output angular rate ω b = [ ω b x ω b y ω b z ] T is equal to the attitude rate ω n b , and the relationship between the Euler angles and ω b can be obtained as follows:
[ ω b x ω b y ω b z ] = [ sin ( θ ) 0 1 cos ( θ ) sin ( γ ) cos ( γ ) 0 cos ( θ ) cos ( γ ) sin ( γ ) 0 ] [ ψ ˙ θ ˙ γ ˙ ]
Thus:
{ ω b x = sin ( θ ) ψ ˙ + γ ˙ ω b y = cos ( θ ) sin ( γ ) ψ ˙ + cos ( γ ) θ ˙ ω b z = cos ( θ ) cos ( γ ) ψ ˙ + sin ( γ ) θ ˙
Due to the high dynamic characteristics of the projectile, the distance spread is much larger than the direction spread, and the change of yaw angle is very little, the change is generally only about 2° [18]. Therefore, it can be assumed that ψ ˙ = 0 . Thus the prediction model of gyro signal for projectile attitude measurement can be obtained by the simplified formula (2):
{ ω b x = γ ˙ ω b y = cos ( γ ) θ ˙ ω b z = sin ( γ ) θ ˙
After the projectile is ejected, the angular rate of the projectile self-rotation is exponentially attenuated, and the process can be described by the Roggla formula [19]:
γ ˙ ( t ) = ω g e ( 0.075 k L d D 3 A t )
where, γ ˙ ( t ) is the angular rate of the projectile self-rotation, ω g is the angular rate of initial self-rotation, L d is the length of the projectile, D is the diameter of the projectile, A is the polar moment of inertia of the projectile, and k is the correction coefficient.

3.2. Function Extension Method

It can be known from formulas (3) and (4) that the angular motion of the gyro-sensitive projectile has a certain regularity. Therefore, the boundary extension method can be designed according to this regularity.
Set: a = ω g , b = 0.075 k L d D 3 A , then formula (4) can be abbreviated as:
γ ˙ ( t ) = a e b t
By integral of formula (5), the following formula can be obtained:
γ ( t ) = γ ˙ ( t ) d t = a e b t d t = a b e b t + c
Combining formulas (3), (5) and (6), a function extension method can be designed and its basic idea is shown in Figure 5.
It can be seen from Figure 5 that the data in the sliding window is fitted by the method. The prediction data of the right boundary extension can be obtained.
In the process of function extension, the sliding window length is L. In the sliding window, the change of pitch angle is small, and it can be considered that d = θ ˙ is constant. Therefore, simplified formula (3) can obtain as follows:
{ ω b x = a e b t ω b y = d cos ( a b e b t + c ) = d cos ( ω b x b + c ) ω b z = d sin ( a b e b t + c ) = d sin ( ω b x b + c )
According to formula (7), the function is fitted and the right boundary data is predicted. The specific operations are as follows:
For the x-axis, the logarithm on both sides as follows:
ln ( ω b x ) = ln ( a ) b t
For the sliding window with length L, the least square method is used for curve fitting. W x = [ ln ( ω b x , 1 ) ln ( ω b x , L ) ] T , T x = [ 1 t 1 1 t L ] , x = [ ln ( a ) b ] T . According to the least square method:
x = ( T x T T x ) 1 T x T W x
For the original data sliding window with length L, the extension length is also L, then a new sliding window data are formed by x-axis prediction and fitting results as follows:
[ ω ^ b x , 1 ω ^ b x , L ω ^ b x , L + 1 ω ^ b x , 2 L ] T = [ a e b t 1 a e b t L a e b t L + 1 a e b t 2 L ] T
For the y-axis, the parameter b can be obtained by fitting the x-axis data. Set ω x = ω b x b , then ω b y = d cos ( ω x + c ) . According to the formula of trigonometric function, the formula (7) can be expanded as follows:
ω b y = d cos ( ω x + c ) = d cos ( c ) cos ( ω x ) d sin ( c ) sin ( ω x )
The least square method is also used for curve fitting, W y = [ ω b y , 1 ω b y , L ] T , T y = [ cos ( ω x , 1 ) sin ( ω x , 1 ) cos ( ω x , L ) sin ( ω x , L ) ] , y = [ d cos ( c ) d sin ( c ) ] T . According to the least square method:
y = ( T y T T y ) 1 T y T W y
For the original data sliding window with length L, the extension length is also L, then a new sliding window data are formed by y-axis prediction and fitting results as follows:
[ ω ^ b y , 1 ω ^ b y , L ω ^ b y , L + 1 ω ^ b y , 2 L ] T = [ d cos ( a e b t 1 b + c ) d cos ( a e b t L b + c ) d cos ( a e b t L + 1 b + c ) d cos ( a e b t 2 L b + c ) ] T
For the z-axis, the y-axis curve can be used directly to fit the parameters. For z-axis, y-axis curve fitting parameters can be used directly.
The flowchart of real-time wavelet threshold de-noising based on function extension method is shown in Figure 6.
In Figure 6, ca represents low frequency wavelet coefficient, cd represents high frequency wavelet coefficient, and threshold processing mainly deals with high frequency coefficients, and j is the number of decomposition layers. The z-axis de-noising process is the same as the y-axis.
After the boundary extension, the wavelet decomposition, threshold processing, wavelet reconstruction and de-noising results are extracted by combining Figure 2 and Figure 6.

4. Algorithm Verification and Result Analysis

4.1. Evaluation Indexes of Algorithm

For the simulation, the SNR and root mean square error (RMSE) are used as evaluation indexes because of the true value [20].
(1) SNR:
S N R = 10 l o g 10 ( N ω 2 ( t ) N ( ω ( t ) ω ˜ ( t ) ) 2 )
where, ω ( t ) represents the true value of the signal, ω ˜ ( t ) represents the filtering result, N represents the signal length, and the SNR unit is decibel (db). The larger the SNR is, the more useful signals are, and the better the de-noising effect is.
(2) RMSE:
R M S E = 1 N N ( ω ( t ) ω ˜ ( t ) ) 2
The smaller the RMSE is, the closer the filtering result is to the real signal, and the better the de-noising effect is.
For the experiment, it is difficult to obtain the true value. The smoothness is used as the evaluation index of the de-noising effect, and the smoothness can reflect the smoothness of the de-noising signal. The smaller the smoothness, the smoother the curve and the better the de-noising effect. Its expression is as follows:
r = N 1 ( ω ˜ ( t + 1 ) ω ˜ ( t ) ) 2 N 1 ( ω ¯ ( t + 1 ) ω ¯ ( t ) ) 2
where, ω ¯ ( t ) represents the original signal.
When de-noising based on the function extension wavelet is completed, the single-sample algorithm is used to the attitude calculation, and the AMAE and the ARMSE are defined as the evaluation indexes.
(1) AMAE:
A M A E = 1 3 N ( N | θ ( t ) θ L ( t ) | + N | γ ( t ) γ L ( t ) | + N | ψ ( t ) ψ L ( t ) | )
where, [ θ ( t ) γ ( t ) ψ ( t ) ] is the result of attitude calculation using the true value of the signal; [ θ L ( t ) γ L ( t ) ψ L ( t ) ] is the result of attitude calculation after de-noising by window length L function extension.
(2) ARMSE:
A R M S E = 1 3 N ( N ( θ ( t ) θ L ( t ) ) 2 + N ( γ ( t ) γ L ( t ) ) 2 + N ( ψ ( t ) ψ L ( t ) ) 2 )

4.2. Simulation Comparison of Different Extension Methods

Simulation conditions: under standard meteorological conditions, the full ballistic attitude data is simulated. The initial rotation rate of the projectile is 10 r/s, the sampling frequency is 1000 Hz, and the gyro bias stability is 10°/h. According to multiple simulations, the relevant parameters are selected as shown in Table 1.
According to the characteristics of the wavelet transform, the length of sliding window L = 4, 8, 16, 32, 64, 128 are selected. The real-time wavelet threshold de-noising is performed on the three axial gyro data by using the methods of no extension, zero extension, periodic extension, symmetric extension, linear extension and function extension respectively. The simulation results of the six extension methods and the original signal SNR and RMSE are shown in Figure 7 and Figure 8.
It can be seen from Figure 7 and Figure 8 that the function extension method is superior to other methods, and with the increase of window length, the SNR increases, and RMSE decreases. For the three axial gyros, the results of zero extension and linear extension are not as good as the original signal, which indicates that the two methods are not applicable in high dynamic environment. For the y-axis and the z-axis, the effect of no extension and periodic extension is poor, which indicates that the real-time wavelet threshold de-noising for these two axes needs to be extended, and the extension method has a great influence on the de-noising effect. The performance improvement (PI) is defined as the increase percentage of SNR for different extension methods relative to the original signal. The window length L = 128 is taken as an ex-ample. The SNR and PI results of different methods are shown in Table 2.
It can be seen from Table 2 that compared with the original signal, for the x-axis, the zero extension and the linear extension do not improve the SNR, and the performance of the function extension method increases the most. For the y-axis, only symmetric extension and function extension have PI. Due to the high dynamic environment, the y-axis data changes drastically in a short time, other extension methods are not suitable, and the PI of the function extension method is 127.69%, which proves the practicability of this method. For the z-axis, no extension, symmetric extension and function extension have PI.
In order to further illustrate the applicability of the pro-posed method, taking the z-axis signal as an example, the de-noising results using different extension methods are shown in Figure 9.
It can be seen from Figure 9 that for the z-axis data change drastically, there will be signal delay phenomenon by using other extension methods, and the de-noising result cannot follow the original signal change normally. Because the high-frequency wavelet coefficients are set to zero by the wavelet threshold. The high-frequency component of the signal has a large attenuation amplitude, the time delay is gradually increased, and the signal amplitude is gradually decreased, and the dynamic response performance is obviously reduced. The function extension method can follow the original signal normally. Because the function fitting result is used to the wavelet threshold de-noising. In the process of extension processing, some high-frequency signals have been removed, which ensures the accuracy of wavelet threshold de-noising.

4.3. Experiment Comparison of Different Extension Methods

According to the actual experiment data of projectile, the performance of different methods are verified, and the parameter setting is the same as the simulation part. Smoothness is used as an evaluation index for de-noising effects. The real-time wavelet threshold de-noising is performed on the three axial gyro data by using the six extension methods respectively. The simulation results of the six extension methods and the original signal smoothness are shown in Figure 10.
It can be seen from Figure 10 that the function extension method is superior to other methods, and with the window length increases, the smoothness decreases within 64. Because the actual projectile motion is more complicated than the simulation motion, it may contain regular precession, general angular velocity motion and random angular vibration. Therefore, the simulation results are different from the experiment results, but in terms of the overall trend, the simulation results and experiment results have certain commonality. For the three axial gyro de-noising results, the linear extension results are not as good as the original signal, and the periodic extension and zero extension are not good, which is similar to the simulation results.
In order to further illustrate the applicability of the proposed method, taking the z-axis signal as an example, the de-noising results using different extension methods are shown in Figure 11.
It can be seen from Figure 11 that although there are some differences between the experiment and the simulation, the effect of wavelet threshold de-noising using different extension methods is similar.

4.4. Real-Time Verification Comparison of Different Extension Methods

The AMD A8-4599M processor with the main frequency of 1.90 GHz is used for simulation in MATLAB R2014b. The actual experiment data is used as the simulation object, and the single calculation period is used as the evaluation index. The specific results are shown in Table 3.
It can be seen from Table 3 that with the window length increases, the single calculation period of the different extension methods does not increase significantly.
The results of simulation and experiment show that the function extension method is superior to other extension methods.

4.5. Accuracy of Attitude Calculation Comparison of Different Window Length

In order to facilitate the description, the function extension of six kinds of length windows L = 4, 8, 16, 32, 64, 128 are defined as “function-4”, “function-8”, “function-16”, “function-32”, “function-64”, “function-128”; the original noise signal is de-fined as “noise-data”, and the simulation results are shown in Figure 12.
It can be seen from Figure 12 that the noise has a great influence on the attitude calculation result. De-noising by function extension method can weaken the influence of noise, and the attitude calculation results are different with different window length selection. According to Figure 12, the AMAE and ARMSE are shown in Table 4.
It can be seen from Table 4 that AMAE and ARMSE of the original signal attitude calculation are relatively large, indicating that the noise has a great influence on the attitude calculation, and it is necessary to de-noise the gyro signal. Compared with the original signal, AMAE and ARMSE are significantly reduced by wavelet de-noising based on function extension, which indicates that the function extension method has practicality. The window length affects the function extension processing effect, and when L = 16, AMAE and ARMSE are the smallest and the de-noising effect is the best.

5. Conclusions

In this paper, the problem of poor wavelet threshold de-noising effect of gyro signal is studied by the conventional extension method in high dynamic environment. According to the flight characteristics of the projectile, the gyro signal prediction model of the projectile attitude measurement is established by introducing the Roggla formula. A boundary extension method of function fitting is proposed. The least square method is used to fit the data of the sliding window, and the prediction of the next sliding window data is realized. The fitting and prediction results form a new sliding window, which is de-noised by wavelet threshold. Simulation and experiment results show that:
(1)
Compared with other extension methods, the gyro signal output characteristics of the projectile attitude measurement is combined, the SNR is improved greatly, the RMSE of the signal is reduced, the smoothness of the signal increased, and the AMAE and ARMSE are reduced, and the real-time performance of the algorithm can be guaranteed by the function extension method, which has great research significance in practical applications.
(2)
For the gyro y-axis and z-axis data, other extension methods will produce time delay, while the function extension method can stably follow the signal, effectively removing random noise and ensuring the accuracy of gyro data, which have a strong practical application value.

Author Contributions

J.W. and Z.D. originated the research and designed the method. X.L. and N.L. helped in results analysis and paper revision. All authors contributed to the written manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China under Grant, grant number 61801032.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shang, J.; Deng, Z.; Fu, M.; Wang, S. Advance and Perspective on Spin Rate Measurement Technology for Guided Projectile. Acta Autom. Sin. 2016, 42, 1620–1629. [Google Scholar]
  2. Ahmed, H.; Tahir, M. Accurate Attitude Estimation of a Moving Land Vehicle Using Low-Cost MEMS IMU Sensors. IEEE Trans. Intell. Transp. Syst. 2017, 18, 1723–1739. [Google Scholar] [CrossRef]
  3. Gong, X.; Zhang, J.; Fang, J. A Modified Nonlinear Two-Filter Smoothing for High-Precision Airborne Integrated GPS and Inertial Navigation. IEEE Trans. Instrum. Meas. 2015, 64, 3315–3322. [Google Scholar] [CrossRef]
  4. Agrawal, N.; Kumar, A.; Bajaj, V.; Singh, G.K. Design of Bandpass and Bandstop Infinite Impulse Response Filters Using Fractional Derivative. IEEE Trans. Ind. Electron. 2019, 66, 1285–1295. [Google Scholar] [CrossRef]
  5. Mallahzadeh, A.R.; Bod, M. Method for designing low-pass filters with a sharp cut-off. IET Microw. Antennas Propag. 2014, 8, 10–15. [Google Scholar] [CrossRef]
  6. Shmaliy, Y.S.; Lehmann, F.; Zhao, S.; Ahn, C.K. Comparing Robustness of the Kalman, H_∞, and UFIR Filters. IEEE Trans. Signal Process. 2018, 66, 3447–3458. [Google Scholar] [CrossRef]
  7. Mallat, S.G.; Heil, C.; Walnut, D.F. A Theory for Multiresolution Signal Decomposition: The Wavelet Representation. Fundam. Pap. Wavelet Theory 2009, 11, 494–513. [Google Scholar]
  8. Chen, X.; Huang, L.; Fu, P.; Gao, G.; He, S.; Shen, J.; Zhu, L.; Dong, L. The Design on the Real-Time Wavelet Filter for ITER PF AC/DC Converter Control System. IEEE Trans. Plasma Sci. 2016, 44, 1178–1186. [Google Scholar] [CrossRef]
  9. Kharitonenko, I.; Zhang, X.; Twelves, S. A wavelet transform with point-symmetric extension at tile boundaries. IEEE Trans. Image Process. 2002, 11, 1357–1364. [Google Scholar] [CrossRef] [PubMed]
  10. Mallat, S.; Zhong, S. Characterization of signals from multiscale edges. IEEE Trans. Pattern Anal. Mach. Intell. 1992, 14, 710–732. [Google Scholar] [CrossRef] [Green Version]
  11. Lei, S.; Jian-Shu, L. Research of point-symmetric boundary extension in wavelet transform. J. Comput. Appl. 2008, 28. [Google Scholar] [CrossRef]
  12. Jianing, S.; Jiawei, Q.; Guodong, X. Modified kernel regression method for the de-noising of X-ray pulsar profiles. Adv. Space Res. 2018, 62, 683–691. [Google Scholar]
  13. Yang, A.; Hou, Z.; Wang, C.; Ding, X.; Gao, Z. Construction of Wavelet Transform Matrices with Symmetric Boundary-Extension. In Proceedings of the 8th international Conference on Signal Processing, Beijing, China, 16–20 November 2006; Volume 2, pp. 1–4. [Google Scholar]
  14. Wei, J.; Pickering, M.R.; Frater, M.R.; Arnold, J.F.; Boman, J.A.; Zeng, W. Tile-boundary artifact reduction using odd tile size and the low-pass first convention. IEEE Trans. Image Process. 2005, 14, 1033–1042. [Google Scholar] [PubMed]
  15. Qian, L.; Xu, G.; Tian, W.; Wang, J. A novel hybrid EMD-based drift de-noising method for a dynamically tuned gyroscope (DTG). Measurement 2009, 42, 927–932. [Google Scholar] [CrossRef]
  16. Ishibashi, R.; Suzuki, T.; Kyochi, S.; Kudo, H. Image Boundary Extension with Mean Value for Cosine–Sine Modulated Lapped/Block Transforms. IEEE Trans. Circuits Syst. Video Technol. 2019, 29, 1–11. [Google Scholar] [CrossRef]
  17. Gaci, S. The Use of Wavelet-Based De-noising Techniques to Enhance the First-Arrival Picking on Seismic Traces. IEEE Trans. Geosci. Remote Sens. 2014, 52, 4558–4563. [Google Scholar] [CrossRef]
  18. Shang, J.; Deng, Z.; Fu, M.; Wang, S. A High-Spin Rate Measurement Method for Projectiles Using a Magnetoresistive Sensor Based on Time-Frequency Domain Analysis. Sensors 2016, 16, 894. [Google Scholar] [CrossRef] [PubMed]
  19. Wang, Y. Half-Experiential Formulas for Calculating Decreasing Angular Velocity of Projectile in Trajectory. J. Detect. Ctrl. 2003, 25, 1–6. [Google Scholar]
  20. Yi, T.H.; Li, H.-N.; Zhao, X.-Y. Noise Smoothing for Structural Vibration Test Signals Using an Improved Wavelet Thresholding Technique. Sensors 2012, 12, 11205–11220. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Actual measurement data of gyro.
Figure 1. Actual measurement data of gyro.
Sensors 20 00200 g001
Figure 2. Flowchart of real-time wavelet threshold de-noising.
Figure 2. Flowchart of real-time wavelet threshold de-noising.
Sensors 20 00200 g002
Figure 3. Common boundary extension methods.
Figure 3. Common boundary extension methods.
Sensors 20 00200 g003
Figure 4. Definition of projectile attitude measurement coordinate system.
Figure 4. Definition of projectile attitude measurement coordinate system.
Sensors 20 00200 g004
Figure 5. Function extension method.
Figure 5. Function extension method.
Sensors 20 00200 g005
Figure 6. Flowchart of real-time wavelet threshold de-noising based on function extension.
Figure 6. Flowchart of real-time wavelet threshold de-noising based on function extension.
Sensors 20 00200 g006
Figure 7. Simulation results of SNR for different extension methods.
Figure 7. Simulation results of SNR for different extension methods.
Sensors 20 00200 g007
Figure 8. Simulation results of RMSE for different extension methods.
Figure 8. Simulation results of RMSE for different extension methods.
Sensors 20 00200 g008
Figure 9. De-noising results of different extension methods.
Figure 9. De-noising results of different extension methods.
Sensors 20 00200 g009
Figure 10. Experiment results of smoothness for different extension methods.
Figure 10. Experiment results of smoothness for different extension methods.
Sensors 20 00200 g010
Figure 11. De-noising results of different extension methods for the actual data.
Figure 11. De-noising results of different extension methods for the actual data.
Sensors 20 00200 g011
Figure 12. Attitude calculation error of different window length function extension.
Figure 12. Attitude calculation error of different window length function extension.
Sensors 20 00200 g012
Table 1. Simulation parameters.
Table 1. Simulation parameters.
ParametersAngle Random WalkWavelet BaseDecomposition LayersThreshold PrincipleThreshold Function
Axial
x 1.5 ° / h Bior 2.65Fixed thresholdSoft threshold
y 1 ° / h Coif 54Unbiased Risk Estimate
z 2 ° / h
Table 2. SNR and PI results of different extension methods.
Table 2. SNR and PI results of different extension methods.
MethodsAxial
xyz
SNR/dbPISNR/dbPISNR/dbPI
Original signal172.81160%31.35780%19.08730%
No extension195.084612.89%26.69530%26.266937.61%
Zero extension133.32050%11.83730%11.90500%
Periodic extension192.080811.15%5.57270%5.22390%
Symmetric extension200.498316.02%32.87804.85%30.932662.06%
Linear extension151.97840%16.20210%6.95400%
Function extension203.843517.96%71.3984127.69%63.2918231.59%
Table 3. The single calculation period results of different extension methods (ms).
Table 3. The single calculation period results of different extension methods (ms).
MethodsWindow Length
48163264128
No extension4.7834.5014.2354.0594.1024.414
Zero extension4.3354.2664.2594.2624.4874.876
Periodic extension4.9794.1894.3284.6824.5984.341
Symmetric extension4.3584.1334.1624.6084.4524.629
Linear extension4.8605.4795.3775.4145.7255.767
Function extension4.0984.9914.6375.9235.6746.108
Table 4. AMAE and ARMSE results of different window length function extension (°).
Table 4. AMAE and ARMSE results of different window length function extension (°).
LengthNoise-DataFunction-4Function-8Function-16Function-32Function-64Function-128
Indexes
AMAE6.70385.56815.19574.64104.68114.68484.8271
ARMSE8.07666.92996.52715.90795.93795.91636.0618

Share and Cite

MDPI and ACS Style

Deng, Z.; Wang, J.; Liang, X.; Liu, N. Function Extension Based Real-Time Wavelet De-Noising Method for Projectile Attitude Measurement. Sensors 2020, 20, 200. https://doi.org/10.3390/s20010200

AMA Style

Deng Z, Wang J, Liang X, Liu N. Function Extension Based Real-Time Wavelet De-Noising Method for Projectile Attitude Measurement. Sensors. 2020; 20(1):200. https://doi.org/10.3390/s20010200

Chicago/Turabian Style

Deng, Zhihong, Jinwen Wang, Xinyu Liang, and Ning Liu. 2020. "Function Extension Based Real-Time Wavelet De-Noising Method for Projectile Attitude Measurement" Sensors 20, no. 1: 200. https://doi.org/10.3390/s20010200

APA Style

Deng, Z., Wang, J., Liang, X., & Liu, N. (2020). Function Extension Based Real-Time Wavelet De-Noising Method for Projectile Attitude Measurement. Sensors, 20(1), 200. https://doi.org/10.3390/s20010200

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