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.
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
(1) formula called
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.
The basic idea of the Mallat algorithm is as follows: Let
From Multi-resolution analysis:
Take the inner product in the side of the equation with
From multi-resolution analysis, we can know that any function
Among
The formula (6) and (7) called Wavelet decomposition algorithm of Mallat algorithm, among wherein {
The reconstruction algorithm of mallat algorithm is the inverse process of its decomposition algorithm. the convolution of mallat algorithm is represented:
Among
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.
THE TEST RESULTS OF A REMOTE SENSING DATA 1 STATIONARITY
Time | 24 | 60 | 120 | 240 |
---|---|---|---|---|
Mean Value | 608.618 | 620.9572 | 622.6287 | 622.6287 |
Variance | 198.2080 | 604.1601 | 611.2902 | 674.6010 |
THE TEST RESULTS OF A REMOTE SENSING DATA 2 STATIONARITY
Time | 24 | 60 | 120 | 240 |
---|---|---|---|---|
Mean Value | 29.2479 | 35.413 | 37.1234 | 39.2378 |
Variance | 10.3366 | 30.9019 | 37.4902 | 42.5010 |
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.
Data of 1 consecutive 4 hours curve graph
Data of 2 consecutive 4 hours curve graph
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).
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].
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.
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.
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
Meet the following conditions:
The forecast model of the telemetry data is:
Set
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:
The predicted values of the reconstructed sub sequences of the scales are compared with the original output power trends (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]:
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.
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.
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:
to
The ideal point:
Negative ideal point:
The proximity of the ideal point
The proximity of the negative ideal point
By calculating the formula (1) (11) (12) is as follows:
The ideal point:
Negative ideal point:
Calculate the relative close degree by the formula (13) (14) (15):
From the relative closeness of view, the time series forecasting model program has a higher rationality and reliability.
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.