Application of Wavelet Analysis in The Prediction of Telemetry Data

## Publications

/ Export Citation / / / Text size:

#### International Journal of Advanced Network, Monitoring and Controls

Xi'an Technological University

Subject: Computer Science , Software Engineering

eISSN: 2470-8038

8
29
Visit(s)
0
Comment(s)
0
Share(s)

SEARCH WITHIN CONTENT

FIND ARTICLE

Volume / Issue / page

Archive
Volume 5 (2020)
Volume 4 (2019)
Volume 3 (2018)
Volume 2 (2017)
Volume 1 (2016)
Related articles

VOLUME 4 , ISSUE 2 (Oct 2019) > List of articles

### Application of Wavelet Analysis in The Prediction of Telemetry Data

Citation Information : International Journal of Advanced Network, Monitoring and Controls. Volume 4, Issue 2, Pages 28-34, DOI: https://doi.org/10.21307/ijanmc-2019-044

Published Online: 08-October-2019

### ARTICLE

#### ABSTRACT

With the rapid development of space technology, the increasing number of spacecraft, in-orbit risk also increases, how to ensure that the spacecraft safety and reliability is particularly important. Prediction technology can predict the failure of the spacecraft in advance, and it has won valuable time for the fault of the spacecraft troubleshooting, thereby increasing the safety and reliability of spacecraft operation. In this paper, based on the non-stationary and periodicity of telemetry data. Based on the wavelet analysis, the prediction of the data is introduced, the establishment of a short-term forecasting model based on Mallat algorithm. The experimental results show that the prediction curve is basically consistent with the actual curve.

## I. INTRODUCTION

The prediction technology of spacecraft fault has been a hot research field. After 20 years development of prediction theory, until the discrete parameters of the linear model of a finite parameter linear model is proposed, and it is possible to combine the prediction theory with the computer. According to the different properties of the forecast, the forecasting methods are generally divided into two categories: time series forecasting and causal prediction. Time series prediction is made by the past predict the future value of the prediction, and causal forecasting is through the known variables to predict the values of other variables. In this paper, the time series forecasting method is used to forecast the future development trend of the telemetry data.

## II. WAVELET ANALYSIS THEORY

### A. Wavelet analysiss

The wavelet analysis method has the characteristics of low frequency and high frequency of the non-stationary signal that change with the low-frequency information signals using a wide time window, high frequency information using a narrow time window. Wavelet is a small area of the wave, waveform with special length, average of 0. Wavelet are defined as follows[1].

Set ψ(t) to one square integrable function, namely ψ(t) ∈ L2(R), if the Fourier transform to meet the conditions:

##### (1)
$Cψ=∫R|ψ^(ω)||ω|dω<∞$

(1) formula called ψ(t) is a basic wavelet or wavelet generating function. When the generating function ψ(t) is expanding and translating, it can get function ψa,τ(t):

##### (2)
$ψa,τ(t)=1|a|ψ(t−τa),amongthema,τ∈R;a≠0$

In (2) formula, a is the scaling factor, t is the translation factor, Because the value of scale factor and translation factor is continuously changing, and depends on the parameters, it is a set of sequence of functions which are obtained by the expansion and translation of the generating function, also called sub-wavelet.

### B. Mallat algorithm

The basic idea of the Mallat algorithm is as follows: Let Hjf as the approximation of the energy limited signal fL2(R) in the resolution 2j, Then the Hjf is further decomposed into the approximation of Hj−1f under the f resolution 2j−1, and the details of Dj−1f between 2j−1 and 2j.

#### 1) Mallat algorithm based on wavelet decomposition

From Multi-resolution analysis: L2(R) = ⊕jZWj, To arbitrary function f(t) ∈ L2(R), get

##### (3)
$f(t)=Σj,k∈Zckjψj,k(t)$

Take the inner product in the side of the equation with ψj,k, because {ψj,k(t)}j,kZ is the orthonormal basis of L2(R), get $ckj=〈f,ψj,k〉$, thus to be

##### (4)
$f(t)=Σj,k∈Z〈f,ψj,k〉ψj,k(t)$

From multi-resolution analysis, we can know that any function fj of Vj, can be expressed as the following form L2(R),

$fj=fj−1+dj−1=fj−2+dj−2+dj−1=…=fM+dM+dM+1+…+dj−1$

Among

##### (5)
$fj(t)=Σkckjφj−1,k(t)+Σkdkj−1ψj−1,k(t)$

fj represents the low frequency components of fM(t), while dl(t), l = M,…,j−1 indicates the high frequency components of fj at different resolutions. Because of $fj(t)=Σkckjφj−1,k(t)=Σkckj−1φj−1,k(t)+Σkdkj−1ψj−1,k(t)$ and ϕ, ψ binary translation and scalability of orthogonality, Can be obtained

##### (6)
$ckj−1=Σncnj<φj,n,φj−1,k>=Σncnjhn−2k*$
##### (7)
$dkj−1=Σncnj<φj,n,ψj−1,k>=Σncnjgn−2k*$

The formula (6) and (7) called Wavelet decomposition algorithm of Mallat algorithm, among wherein {hk}kZ is a filter coefficient sequence by a two-scale equation corresponding orthogonal scaling functions.

#### 2) Reconstruction algorithm of mallat algorithm

The reconstruction algorithm of mallat algorithm is the inverse process of its decomposition algorithm. the convolution of mallat algorithm is represented:

##### (8)
${cj−1=D(cj*h¯*)dj−1=D(cj*g¯*)cj=(Ucj−1)*h+(Udj−1)*g$

Among $h¯*$, is represented conjugate inversion of filter h; $cj*h¯*$ represent conjugate of cj and $h¯*$; $D(cj*h¯*)$ represent Under the dual sampling of conjugate $cj*h¯*$.

## III. THE RESEARCH OF TELEMETRY DATA TIME SERIES PREDICTION BASED ON MALLAT ALGORITHM

### A. The characteristics of telemetry data

Telemetry data has the characteristics of non-stationary variation, commonly used statistics of the telemetry data (such as the mean and autocorrelation function, etc.) often varies with time changing, it bring very great difficulty to the telemetry data forecast. Through the telemetry data 1 and 2 (table 1, 2) statistics, difference is very big, every stage of the statistical parameters show that the sequence of non-stationary time series. Wavelet analysis to deal with this kind of data has a great advantage.

##### TABLE I.

THE TEST RESULTS OF A REMOTE SENSING DATA 1 STATIONARITY

##### TABLE II.

THE TEST RESULTS OF A REMOTE SENSING DATA 2 STATIONARITY

Figures 1, 2 is telemetry data 1 and 2 for four hours of change curve, it can be seen that the output power sequence is periodicity, as well as randomness. The coexistence of periodicity and randomness, the result can be seen as the superposition of different frequency components, these frequency components superimposed on each other in the interior have similar frequency characteristics and the same variation. If the subsequence to establish a prediction model for single change, due to the change of the data characteristics of a single, reduce the difficulty of forecasting model selection.

##### Figure 1.

Data of 1 consecutive 4 hours curve graph

##### Figure 2.

Data of 2 consecutive 4 hours curve graph

### B. The choice of wavelet function

There are some mutations in the trend of spacecraft telemetry data, and these mutations reflect the actual state of the satellite. In order to accurately capture the point of the mutation, the wavelet function is usually required to have a fast convergence, which can quickly attenuate to zero [5]. In this paper, db3 wavelet is chosen as the wavelet base for the different scale of a certain remote data sequence. DbN wavelet in N=1, 3, 4, 1, the results of 2 scale decomposition of a telemetry data (Figure 3).

##### Figure 3.

Comparison results of dbN wavelet based N=1,2,3,4

The choice of wavelet function can satisfy the following 3 conditions except that the condition and the regularity condition[2].

• 1) Good compact support;

• 2) ψ(t) has vanishing moments; vanishing moments can make the wavelet function has a good locality in the frequency domain;

• 3) Satisfy orthogonality.

### C. Wavelet decomposition scale study

Due to the telemetry data is not stable, its change cycle is difficult to see, and its change is slow and fast change together, that is, the change of telemetry data cycle is the size of the cycle of nested together. Therefore, Separating different frequency component can make its change rule is more intuitive, and can also improve data non-stationary[6]. The figure 4 Show the results that comparison db3 wavelet approximation part aN decomposed at different scales and a0 of the original sequence.

##### Figure 4.

Comparison of different decomposition scale approximation section

It can be seen that the decomposition scale is 2, the curve of the approximate part a2 has smooth enough, and basically maintained the shape of the original curve, while A3, A4, because of the increase in the number of points, the sampling point is reduced, the approximate partial curve is too smooth, the sequence of changes in the trend has been distorted, so this paper chose the decomposition of 2.

### D. Forecasting model of time series

Time series forecasting[3][7] is one of the methods of statistical analysis. Its modeling idea is the basic assumption that the change of the past of the telemetry data will continue into the future, that is, the future is a continuation of the past. In this paper, the use of the periodic autoregressive model (PAR model)[8][9] is as follows:

If there is a time series X, the expression is

##### (9)
$Xt=a0t+a1,tXt−1+a2,tXt−2+…+ap,txt−p+ϵt$

Meet the following conditions:

• 1) εt is independent sequence, Expected Value is Eεt = 0, variance is $Eϵt2=σt2$;

• 2) For any i = 0,1,…,p, ait = ait+T, $σt2=σit+T2$, t = 0,±1,…, among T is a positive number, and the model is model of PAR, T is the Cycle length of PAR model, t is phase of PAR model.

The forecast model of the telemetry data is:

$XKT+1=a01+a11XKT+…+aγ1XKT+1−γ+ϵKT+1XKT+2=a2+a2XKT+1+…+aγ2XKT+2−γ+ϵKT+2……XKT+T=a0T+a1TXKT+T−1+…+aγTXKT+T−γ+ϵKT+T$

Set X1, X2,…, Xn is telemetry data samples of per minute, the value of Xt(h) in the future h is the Xt+h in the condition of t, and thus as its predictive value, denoted as $X^n(t)$, according to the definition of:

Decomposition for telemetry data sequence prediction model is established for an hour, the selection of the length of the cycle is 60, namely, p=T=60. This load sequence PAR prediction model can be represented by the expression of the following:

$X^n(k)=akT,k+akT+1,kX^n(k−1)+…+akT+60,kX^n(k−60)akT+i,k(i,k=1,2,…,60)$

### E. Predicted results analysis

The predicted values of the reconstructed sub sequences of the scales are compared with the original output power trends (Figure 5).

##### Figure 5.

Comparison results between predicted values and actual results

By analyzing the comparison between an actual telemetry data and the predicted value, it can be seen that the boundary of the predicted results and the trend of mutation is not very ideal. In this paper, the extension of the periodic continuation to the boundary of the sequence is extended. The idea of periodic continuation is that the signal is considered as a periodic signal, and the extension process is as follows[10][11]:

${xn=xn+M,n<0xn=xn−M,n

Among, M is the length of sequence. The comparison result of the prediction curve and the actual curve after eliminating the boundary error is shown in Figure 6.

##### Figure 6.

The comparison results between the predicted values and the actual values of the modified boundary

From Figure 6, we can see that the results of the prediction of the sequence are better than the results obtained by the wavelet transform.

## IV. THE PERFORMANCE EVALUATION OF THE TELEMETRY DATA FORECASTING MODEL

Thought of optimal decision method is using linear transform to normalize the attribute value, and at the same time using the ideal point and negative ideal point, compared with the traditional method has more rationality and reliability. The ideal point is the best solution, and its target value is the best, the worst is the worst. The worst is the worst. The optimal solution algorithm is as follows:

### A. Set decision matrix of A is

##### (10)
$A=f1(x)f2(x)⋮fp(x)[f1(x1)f1(x2)…f1(xn)f2(x1)f2(x2)…f2(xn)⋮⋮⋮⋮fP(x1)fP(x2)…fP(xn)]$

to

##### (11)
$U−i=min⁡(fi(x1),fi(x2),…,fi(xn)},i=1,2,…,k$
$U¯i=max⁡(fi(x1),fi(x2),…,fi(xn)},i=k+1,k+2,…,PV¯i=max⁡(fi(x1),fi(x2),…,fi(xn)},i=1,2,…,kV−i=min⁡(fi(x1),fi(x2),…,fi(xn)},i=k+1,k+2,…,P$

### B. To determine the ideal point and negative ideal point

The ideal point:

##### (12)
$x*=(U_1,U_2,…,U_k,U¯k+1,U¯k+2,…,U¯P)T$

Negative ideal point:

$x¯=(V¯1,V¯2,…,V¯k,V_k+1,…,V_P)T$

### C. Calculate the close degree

The proximity of the ideal point

##### (13)
$Ri=1P[Σj=1kU_jfj(xi)+Σj=k+1pfj(xi)U¯j],i=1,2,…,n$

The proximity of the negative ideal point

##### (14)
$ri=1P[Σj=1kfj(xi)V¯j+Σj=k+1pV_jfj(xi)],i=1,2,…,n$

### D. Calculate the relative close degree

##### (15)
$ϵi=Riri+Ri,0≤ϵi≤1,i=1,2,…,n$

By calculating the formula (1) (11) (12) is as follows:

The ideal point:

$x*=(U_1,U_2,…,U_k,U¯k+1,U¯k+2,…,U¯P)T=(1.23,0.978,4.219,1.491)$

Negative ideal point:

$x¯=(V¯1,V¯2,…,V¯k,V_k+1,…,V_P)T=(3.76,3.679,7.896,4.118)$

Calculate the relative close degree by the formula (13) (14) (15): ε =0.34

From the relative closeness of view, the time series forecasting model program has a higher rationality and reliability.

## V. CONCLUSION

This paper studies the telemetry data forecasting method based on wavelet analysis. Through the analysis of the characteristics of the telemetry data, the characteristics of the non - stationary and certain periodicity of the telemetry data are established, and the prediction algorithm based on wavelet analysis is established. By choosing different wavelet bases and the decomposition scale, the decomposition results show that 2, db3 wavelet decomposition scale is the best. Based on the Mallat algorithm, the time series forecasting model is established according to the characteristics of detail data and approximate data. The experimental results show that the predicted values are in good agreement with the actual values. Finally, through the analysis of the performance of the forecasting model, the forecasting model is reasonable and reliable.

## References

1. Li-Zhi Cheng, Hong-Xia Wang, Yong Luo. The theory and application of the wavelet[M]. beijing: Science press, 2004:75-116
2. Xiang-Bing Meng. Research and implementation of short term load forecasting model based on wavelet analysis[D]. dalian: Computer application technology of Dalian University of Technology, 2008.
3. Jian Tang, Jia-Hui Luan, Chen Lv. Interval forecasting technique of remote sensing data for small satellite power supply system[J]. Journal of Huazhong University of Science and Technology (NATURAL SCIENCE EDITION), 2009(1):209-212.
4. Jia-Lin Zhang, Xiao-Jun Wei. Based on the optimization model of target decision method and its application research[J]. Financial theory and practice. 2007(145):114-119.
5. Zhen-Ming Sun. Forecasting theory and technology in the application of the spacecraft[D]. haerbing: Harbin institute of technology, 2004.
6. Zhen-Ming Sun, Wei-Guang An, Hui Zhang. Spacecraft data to predict the causal relation adjustment technology application research[J]. Journal of aerospace, 2006(6):1323-1327.
7. Keon-Tae SOHN, Deuk-Kyun RHA. The 3-Hour-Interval Prediction of Ground-Level Temperature[J]. Advanced in Atmosphric Sicences. 2003(4):576-583.
8. Chang-il K, In-keun Y, Y.H. Song. Kohonen neural network and wavelet transform based approach to short-trem load forecasting [J]. Electric Power Systems Research, 2002(3):167-176.
9. Soltani Skander. On the use of the wavelet decmpositoin for time series prediction[J]. Neurocomotuing, 2002(9):267-277.
[CROSSREF]
10. Zheng Hua, Zhang Lizi. The factor Analysis of short-trem laod forecast based on wavelet transform[J]. IEEE, 2002:165-169.
11. Zhang B L, Coggins R. Multiresolution forecasting for futures trading using wavelet decmpositions[J]. IEEE Transactions on Neural Netoworks, 2001, 12(4):765~774.
[CROSSREF]

### FIGURES & TABLES

Figure 1.

Data of 1 consecutive 4 hours curve graph

Figure 2.

Data of 2 consecutive 4 hours curve graph

Figure 3.

Comparison results of dbN wavelet based N=1,2,3,4

Figure 4.

Comparison of different decomposition scale approximation section

Figure 5.

Comparison results between predicted values and actual results

Figure 6.

The comparison results between the predicted values and the actual values of the modified boundary

### REFERENCES

1. Li-Zhi Cheng, Hong-Xia Wang, Yong Luo. The theory and application of the wavelet[M]. beijing: Science press, 2004:75-116
2. Xiang-Bing Meng. Research and implementation of short term load forecasting model based on wavelet analysis[D]. dalian: Computer application technology of Dalian University of Technology, 2008.
3. Jian Tang, Jia-Hui Luan, Chen Lv. Interval forecasting technique of remote sensing data for small satellite power supply system[J]. Journal of Huazhong University of Science and Technology (NATURAL SCIENCE EDITION), 2009(1):209-212.
4. Jia-Lin Zhang, Xiao-Jun Wei. Based on the optimization model of target decision method and its application research[J]. Financial theory and practice. 2007(145):114-119.
5. Zhen-Ming Sun. Forecasting theory and technology in the application of the spacecraft[D]. haerbing: Harbin institute of technology, 2004.
6. Zhen-Ming Sun, Wei-Guang An, Hui Zhang. Spacecraft data to predict the causal relation adjustment technology application research[J]. Journal of aerospace, 2006(6):1323-1327.
7. Keon-Tae SOHN, Deuk-Kyun RHA. The 3-Hour-Interval Prediction of Ground-Level Temperature[J]. Advanced in Atmosphric Sicences. 2003(4):576-583.
8. Chang-il K, In-keun Y, Y.H. Song. Kohonen neural network and wavelet transform based approach to short-trem load forecasting [J]. Electric Power Systems Research, 2002(3):167-176.
9. Soltani Skander. On the use of the wavelet decmpositoin for time series prediction[J]. Neurocomotuing, 2002(9):267-277.
[CROSSREF]
10. Zheng Hua, Zhang Lizi. The factor Analysis of short-trem laod forecast based on wavelet transform[J]. IEEE, 2002:165-169.
11. Zhang B L, Coggins R. Multiresolution forecasting for futures trading using wavelet decmpositions[J]. IEEE Transactions on Neural Netoworks, 2001, 12(4):765~774.
[CROSSREF]