The Adaptive Analysis of Shock Signals on the Basis of Improved Morlet Wavelet Clusters

Morlet wavelets do not satisfy the permissibility condition of wavelet analysis, and there are therefore no inverse transformations for Morlet wavelet transforms. In this paper, we put forward the Yang and Pan transform (YPT), which is an adaptive discrete analysis method for shock signals. First, we improved the Morlet wavelet so that the centre and radius of the frequency window can be easily adjusted in the frequency domain. Second, we proposed the extremum frequency concept and analysed the extremum situation of the improved Morlet wavelet. /ird, combining the improved Morlet wavelet and extremum frequency, we advanced the theory of the YPT, which does not need to satisfy the permissibility condition. We then continued by using a smoothing operator that can smooth the potentially distorted signal reconstructed after being analysed by the YPTand filtered by using the threshold filtering theory. /is operator proved to be simple and efficient. Finally, a noisy signal was reconstructed after being analysed and filtered using the YPTand threshold filtering, respectively, to verify the validity of the theory, and the YPT was compared with the discrete wavelet transform (DWT). As a supplement to the theory in engineering, the shock signals about a gun automatic mechanism were also analysed using the theory in this paper. Good results were obtained, thereby demonstrating that the YPT can be helpful to further extract the features of shock signals in pattern recognition and fault diagnosis.


Introduction
Wavelet transforms are a significant method for analysing signals, and Morlet wavelets are often used as the kernel functions of wavelet transforms.In 1982, Morlet first used them to analyse seismic signals, employing a complex wavelet using a Gauss envelope and a special case of Gabor wavelets as well [1].e continuous wavelet transform of a signal by Morlet wavelets can achieve arbitrary high resolution in the time or frequency domains [2,3].ere is a 90-degree phase shift between its real and imaginary parts, which makes it easy to obtain the instantaneous frequency and phase of the signal, and it therefore has a wide range of applications.For example, Morlet wavelets have been applied to signal filtering and denoising [4], mechanical fault diagnosis [5][6][7][8], the analysis of medical signals [9], research on river runoff in natural environments [10], research on rainwater evaporation [11], problems of polymer pollution [12], the atmospheric system [13], the influence of cosmic rays on organisms [14], and the motion of celestial bodies [15].
However, the Morlet wavelet does not satisfy the permissibility condition of wavelets, and wavelet reconstruction therefore cannot be realised.To obtain the inverse transformation of the Morlet wavelet, many scholars have improved the Morlet wavelet and obtained corresponding reconstruction transforms.Grossmann et al. [16] added proper correction terms to the Morlet wavelet to satisfy the permissibility condition and properly set the parameters for the tuning signal so that the rounding error of the computer had the same order of magnitude as the correction terms so that the correction terms could be omitted [17].Ji and Yan [17], undertaking research at Northwestern Polytechnical University in China, improved the Morlet wavelet.ey fixed some parameters in the Morlet wavelet and put forward the Morlet and Ji transform (MJT) as well.e MJTcan reconstruct a signal without needing to meet the admissible condition.It is a continuous transformation of a signal and has a good time-frequency localisation property.Partially referencing the MJT and on the basis of the Morlet wavelet, we improved the Morlet wavelet in another way, further advancing a method for adaptively and discretely analysing signals.
In this paper, we first summarise the MJT and then propose another method to improve the Morlet wavelet cluster.Subsequently, the time and frequency windows of the Morlet wavelet, which were improved by Ji and Yan [17], are analysed, as is the Morlet wavelet improved by us.ird, the concept of extremum frequency is proposed in this paper, and the extremum properties of the Morlet wavelet cluster as improved by us are then analysed.Based on the improved Morlet wavelet cluster and extremum frequency, this paper then presents a new method that adaptively and discretely analyses a vibration signal and simultaneously yields the YPT, which can completely reconstruct the original signal.Additionally, we continue by putting forward a smooth operator that can smooth the potentially distorted signal reconstructed after analysis using the YPT and filtered employing threshold filtering theory.Finally, through the inverse transformation of the YPT, a filter is constructed to verify the correctness and practicability of the YPT by combining a signal with high noise.At the same time, the YPT is compared with the DWT.As a supplement to the theory in engineering, the shock signals about a gun automatic mechanism are also analysed by using the theory in this paper, which can provide help for further extracting the features of shock signals in pattern recognition and fault diagnosis.

MJT and Improved Morlet Wavelet Cluster
is section first introduces the content of the MJT, and we then improve the Morlet wavelet cluster in the MJT to obtain our improved Morlet wavelet cluster and analyse the properties of the time and frequency windows for the latter improved wavelet cluster.

MJT.
e analytical expression of the Morlet wavelet is e Fourier transform of (1) is ( Given a wavelet ψ(t) ∈ L 1 ∩ L 2 , the permissibility condition of the wavelet transform is We know ψ(0) � 0 by (3), whereas m(0) ≠ 0; thus, the Morlet wavelet does not satisfy the permissibility condition of the wavelet transform.Following Ji and Yan [17], let )e −(t 2 /2)+j2πt , and the modified Morlet wavelet cluster is then obtained by stretching and translating ψ(t).
eir improved Morlet wavelet cluster is Assuming signal h(t) ∈ L 2 , the expression for the MJT is [17] ( According to (5), the Fourier transformation of J(f, τ) with respect to τ can be obtained: erefore, according to (6), the inverse transformation of h(t) can be obtained: dτ e j2πft df.(7) Formulas ( 5) and ( 7) are, respectively, the positive transformation and inverse transformation formulas of the MJT.According to (5), the signal h(t) can be transformed into the f -τ domain for analysis.e MJT is a continuous transformation that has good time-frequency localisation properties [17], and it is a redundant transformation that requires large amounts of calculations and storage space.

Improved Morlet Wavelet
Cluster.In this section, we improve the Morlet wavelet cluster in (4) and analyse the properties of the time and frequency windows for the improved Morlet wavelet cluster.e improved Morlet wavelet cluster, denoted as φ f,δ (t), is expressed as (8).Formula ( 9) is the Fourier transform of (8):

2
Shock and Vibration Cauchy's theorem on the two-dimensional complex plane can be used to compute the definite integral in (9), and it can also be obtained directly from the Poisson integral formula.φ f,0 (t) � ψ f,0 (t) can be obtained by using ( 4) and (8).We will now analyse the properties of the time and frequency windows for φ f,δ (t), the centre of which in the time domain is denoted as t * , which can be obtained by the following equation: In the process of calculating the definite integral in (10), we obtain )sgn(|f| + δ) by utilising the properties of the Gamma function.Of course, as far as (10) is concerned, we do not have to calculate the latter in detail if we obtain the former.We can conclude that t * is independent of f and τ via (10); therefore, the centre of the time domain for ψ f,0 (t) is also 0. Furthermore, the radius of the time window, which is denoted as Δt, can be obtained by the following equation: Obviously, the radius of the time window for ψ f,0 |f|, and f together with δ jointly affects the radius of the time window for φ f,δ (t).When the absolute value of f is large, the radii of the time windows for ψ f,0 (t) and φ f,δ (t) decrease, but the latter can be adjusted by using the parameter δ.Similarly, the centre c * and radius Δc of the frequency window for φ f,δ (t) in the frequency domain are obtained as We can conclude Δt × Δc � 1/4π via ( 11) and ( 13), which as a constant cannot be affected by δ, but it can be separately adjusted by δ, and the width of the frequency window can be linearly adjusted by f and δ; what is more, f determines the centre of the frequency window.

Extremum Frequency
In this section, the concept of extremum frequency, which reflects the average rate of change for a sequence and provides a method for adaptively analysing signals, is first proposed.Based on this concept, the extremum and extremum frequency of the improved Morlet wavelet cluster are analysed, which paves the way for the adaptive analysis of signals in the next section.

3.1.
e Extremum Frequency of a Sequence.
To go a step further, assuming that the sequence x(n)|n � 1, 2, . . ., M { } contains P h amounts of maximum aggregates and P l amounts of minimum aggregates as well, the extremum frequency of x(n)|n � 1, 2, . . ., M { }, denoted as f e , can be defined as

e Extremum of the Improved Morlet Wavelet Cluster.
In this section, we let |f| + δ ≠ 0 and analyse the real part of φ f,δ (t), which is equal to Without loss of generality, let |f| + δ > 0. Some properties regarding the maximum of φ f,δ (t) are now described as follows: (1) ere is only one point , and there is therefore only one point, When f < 0, similarly, there exists only one point (ii) Let n be an odd number; similarly, ξ n is the only minimum point in the interval )) ≠ 0, and with the information given above, the property of (2) has been proved.

e Discrete Sampling of the Improved Morlet Wavelet
Cluster.We can obtain a sequence, which is denoted as from the start time t s to the end time t e and by sampling interval t 0 .Figure 1 shows a schematic diagram for sampling φ R (t).
According to the parameters for sampling φ R (t), there exists the following formula: According to the parameters for sampling φ R (t) and the analysis of the extremum for the improved Morlet wavelet cluster in Section 3.2, we can obtain the numbers of the extreme aggregates for u(n)|n � 1, 2, . . ., N 0  .Let the numbers be P, which can be expressed as Combining ( 14) and ( 20) with (21), the extremum frequency of u(n)|n � 1, 2, . . ., N 0  , denoted as f em , can be obtained: e effect of analysing a sequence whose extremum frequency is f e may be better when we utilise the improved Morlet wavelet, if we make f em � f e .Because the Morlet wavelet, on the condition of f em � f e , will have better adaptability to the sequence, the following formula is therefore established: For a continuous signal h(t) defined in the interval [t s , t e ] and containing P c amounts of extreme points, we can determine that the extremum frequency of the sequence obtained by sampling h(t) is f ec � (P c t 0 )/((t e − t s ) + t 0 ) if setting the sampling interval equal to t 0 .When we utilise the improved Morlet wavelet to analyse h(t), we set f em f ec .Formula (25) is therefore established: It is equal to On the basis of (26), the improved Morlet wavelet, whose parameter f is equal to (P c )/(2(t e − t s )), will have better adaptability when analysing the signal h(t).According to (23)-(25), we obtain f f e /2t 0 P c /2(t e − t s ) on the conditions (t e − t s ) → +∞ and t e − t s ≫ t 0 .  2 and 3.

Analysing a Vibration Signal
Let |f| + δ ≥ 0. According to Figures 2 and 3, the following formula must be true if we make the radii of the two adjacent Morlet wavelets just connect with each other in the frequency domain: It is equal to According to Section 3.3, we know we should make f f e /2t 0 P c /2(t e − t s ), when analysing the signal h 0 (t) by utilising the improved Morlet wavelet.We now transform the signal h 0 (t) by using the following equations, which are collectively named the "YPT": -1

Shock and Vibration
where k ∈ N + in the formulas in (29).On their basis, to go a step further, we take the Fourier transforms of the formulas in (29): e equations in (29) dispersedly and adaptively transform the signal h 0 (t) into the time domain, and on the contrary, the formulas in (30) dispersedly and adaptively transform the signal h 0 (t) into the frequency domain.Taken together, both constitute the positive transformation of the YPT.f k and δ k are adaptive to h k−1 (t).We can also obtain the inverse transformation of the YPT according to the formulas in (29) and (30), which can be described by the following equations: For the convenience of our expression in this paper, Y pk (t; f k , δ pk ) and Y nk (t; −f k , δ nk ) are uniformly named Y(t; f, δ); similarly, h k (t) is designated as h(t), Ŷpk (c; f k , δ pk ) and Ŷnk (c; −f k , δ nk ) are designated as Ŷ(c; f, δ), and ĥk (c) is designated as ĥ(c).According to (31), we can directly reconstruct the original signal h 0 (t) by adding signals Y(t; f, δ) and h(t), having been analysed and processed, in the time domain.In addition, we can continue to further transform signals Y(t; f, δ) and h(t) into signals Ŷ(c; f, δ) and ĥ(c) in the frequency domain and then continue to process signals Ŷ(c; f, δ) and ĥ(c) and reconstruct the original signal h 0 (t) through (32).On a similar principle, according to (32), we can reconstruct the original signal h 0 (t) through signals Ŷ(c; f, δ) and ĥ(c), having been analysed and processed, in the frequency domain.In addition, we can continue to transform signals Ŷ(c; f, δ) and ĥ(c) into signals Y(t; f, δ) and h(t) in the time domain and then continue to process signals Y(t; f, δ) and h(t) and reconstruct the original signal h 0 (t) through (31).erefore, (31) and (32) provide us 4 types of YPT inverse transforms, namely, 2 types of YPT inverse transforms from the frequency domain to the time domain and 2 types of YPT inverse transforms from the time domain to the time domain.

Smoothing Operator.
A reconstructed vibration signal can be seriously distorted in some local areas after being analysed and processed.It is usually necessary to properly smooth the signal for the speci c transformation, which is helpful to re ect the change law of the original signal.In this section, we propose a smooth operator S that can smooth the sequence corresponding to the signal.Before introducing the operator S, three operators, A, B, and L, are introduced rst.Assuming that there is a nite sequence X x 1 , x 2 , . . ., x n , the operator A can perform the following operation on the sequence X: In particular, e operator B performs the following operation: 6 Shock and Vibration

B(X) ≜ A(A(A(. . . A(X, 1) . . .), n − 1), n). (35)
e operator L performs the following operation: (36) Obviously, L(X, n) � L(X, −n) � X.On the basis of operators A, B, and L, we can define the operator S: In particular, S(X, 0) ≜ X. rough a large number of experiments on sequences obtained by the YPT and the filtering theory that we utilised in this paper, we found that the operator S can effectively smooth sequences.However, we cannot excessively smooth a signal.

e Process of Analysing a Vibration Signal.
From the analyses in Sections 4.1 and 4.2, if we analyse the signal h 0 (t) through the YPT, the kernel is determining the extremum frequency of the signal and controlling the radius of the frequency window.
which has been discussed through (28) in Section 4.1, and we can let δ 1 be 0. e radius of the frequency window will be zero if |f k | + δ k � 0 in the process of analysing the signal.However, this matter will make the process of the YPT stagnate.We therefore let δ k � 0 again if |f k | + δ k � 0 to continue transforming the signal until the transformation reaches the required number of the transformation.We can denote this required number as K 0 .e flow chart of the YPT can be divided into two parts, the first of which concerns analysing the signal, and the second concerns reconstructing the signal; both parts are shown in Figures 4 and 5, respectively.e parameter model determines how to reconstruct a signal.For this paper, we chose model � 4 to reconstruct a signal.

Analysis and Reconstruction of a Signal
In this section, the correctness and practicability of the YPT are examined from two aspects: signal analysis and reconstruction.Based on the above research, we analysed a signal according to the flow chart in Figure 4 and reconstructed it according to the flow chart in Figure 5, using model � 4. We first analysed a signal containing high amounts of white noise, and we then reconstructed the signal after the process of filtering it.reshold filtering was adopted as filtering theory to remove the noise in the signal, while the threshold filtering established by us is neither hard threshold filtering nor conventional soft threshold filtering but is rather a theory established on the probability that an interval signal belongs to noise.In summary, first, the probability function that judges the probability that an interval signal belongs to noise is established on the basis of the characteristics that the expectation of the noise is zero and that the noise has a uniform distribution; the threshold is then determined by the probability function, and at last, the interval signals belonging to noise are filtered according to the threshold function.In this paper, this filtering theory is named "probability filtering" and is not discussed in detail.

Analysis of a Signal. Noise was added to the ideal signal
), by utilising the random function rand(t) in the Matlab numerical calculation software environment.
From Figure 6, it is difficult to observe the change regulation of the noisy signal when we set the maximum noise amplitude to 3.After calculation, the signal-to-noise ratio (SNR) of the noisy signal was 3.13 dB. is noise can therefore be regarded as high noise.
e signal h(t) was transformed 6 times using the YPT, and K 0 was therefore equal to 6.For the convenience of calculation, we equally divided the transformed signal h(K 0 ) into two parts when we stored the transformed signals in memory.ere were a total of 14 transformed signals resulting from the YPT.To conveniently observe the transformed signals in the time domain, the real parts of the signals in three-dimensional space are plotted in Figure 7.
Because the noisy signal is real, the two transformed signals whose serial numbers are opposite are conjugates and have the same real part.According to Figure 7, the signals e real parts of the transformed signals in threedimensional space are plotted in Figure 9. From Figure 9, we can clearly that the signals with serial numbers 1 and −1 contain large amounts of noise, which accords with the characteristics that the expectation of the noise is zero and is uniformly distributed in the frequency domain.For the signals with serial numbers 2 and −2, each of them has a positive impact amplitude and a lower negative impact amplitude in the frequency domain.
e real parts of the signals with serial numbers 3 and −3 are separately plotted on two-dimensional planes in Figures 10 and 11    Shock and Vibration and negative signals."eir real parts have symmetry if the signal h(t) is real.

Reconstruction of a Signal.
For this section, we adopted model 4 to reconstruct the signal h(t) from Section 5.1.First, the signals in the frequency domain were ltered employing probability ltering to remove a part of the noise, and we then transformed the ltered signals in the frequency domain into the time domain.After that, we continued to again lter the signals in the time domain using probability ltering to remove a part of the noise in the di erent space.We nally reconstructed the original signal h(t) using the ltered signals in the time domain.e process of reconstructing the original signal employed the YPT inverse transform, the ltering theory, and the smoothing operator.In Figure 12, the real parts of the new ltered signals in the frequency domain are shown in threedimensional space.
Compared with those shown in Figure 9, the signals with serial numbers 1 and −1 Figure 12 signi cantly reduced the noise and retained a part of the useful signals, as was the case for the signals with serial numbers 2 and −2.However, for the serial number 2 signal, the positive impact amplitude was reduced slightly because the positive impact amplitude was a ected by the noise and the di erent ltering order when the threshold function readjusted the values of the interval signals.To compare with the signals in Figures 10  and 11 and to also observe the e ects of the YPT and ltering theory, the real parts of the new ltered signals with serial numbers 3 and −3 are separately plotted on two-dimensional planes in Figures 13 and 14.From Figures 13 and 14 15 generally has less noise.On the contrary, the new ltered signal with serial number 2 lost a large amount of energy in the ltering process, which exposes the lack of this ltering theory but does not in uence the validity and practicability of the YPT.
Finally, according to (31), we reconstructed the original signal h(t) through the new ltered signals in Figure 15.

Shock and Vibration
We then properly smoothed the signal using the smoothing operator S(X, w), and the parameter w was 5. e new signal processed by the YPT and ideal signal without noise correspond to the B curve and the C curve in Figure 16, respectively.From Figure 16, we can determine that the amplitude variation of the new signal was not as high as the ideal signal.Combining Figures 6 and 16, we can conclude that the e ect of reconstructing a signal is ideal when utilising the YPT, probability ltering, and the smoothing operator.Compared with the noisy signal in Figure 6, the noise in the new signal was very low, and it clearly re ected its own law of change.
To further verify the validity of the YPT, it is compared with the DWT.We also decomposed the noisy signal into 4 layers by the DWT, and the wavelets for the DWT are 3 order Daubechies wavelets.For the noisy signal decomposed by the DWT, after removing the rst three layers of detail signals, we can then obtain the best reconstructed signal by the inverse discrete wavelet transform (IDWT), which corresponds to the A curve in Figure 16.After calculation, the SNR of the noisy signal, A curve, and B curve is 3.13 dB, 10.61 dB, and 16.39 dB, respectively, and as far as denoising and the noisy signal in this paper are concerned, the YPT is therefore better than the DWT.

The Application of the YPT in Engineering
In this section, two types of serious fault cracks on the latch sheets of a gun automatic mechanism were set by us, and the shock signals that are obtained by sampling the data, respectively, when the gun automatic mechanism works in normal and two fault patterns were analysed by utilising the method in this paper.e ring frequency of the gun automatic mechanism is 10Hz, and the main parameters of the sensors that we used in this paper are shown in Table 1.To obtain more details in the time domain, we set the sampling

Shock and Vibration
frequency to the maximum sampling frequency (204.8kHz) of the signal collector, and the sampling time is 5s.On the contrary, for a shock signal, the sampling frequency should be as high as possible to utilise the signal near the resonant frequency of the sensor, although the output of the sensor may be distorted at high frequency.e environment of the experiment is shown in Figure 17, and the main measuring point is shown in Figure 17(d).
Due to the bad working environment or high overload of the equipment, a lot of noise is usually contained in the shock signals.e noise reduction for the signals therefore needs to be done to eliminate the e ect of the noise if we do pattern recognition or fault diagnosis for the equipment.For the gun automatic mechanism, there are two main sources of noise, one of which is at ammunition feeding mechanism that is numbered (1) in Figure 17(a) and another is at the base of the equipment, numbered (2) in Figure 17(a).
e former is produced by the irregular friction and collision between the bullet chain and the body of the equipment, and the latter is produced by self-excited vibration on the base that is with the complex structure or the random vibration of the base such as a gun automatic mechanism that works on a tank or ghter plane.ese two kinds of noise have the property that the frequency of them is lower when compared with that of the e cient signals  18(a)-18(c) and are ltered by the method, respectively.To further observe the e ects of the method, we transformed them into the frequency domain, which is shown in Figure 19.From Figure 19, we can conclude that contrast to original signals, the amplitude of the part in low frequency is signi cantly reduced, relatively to the part in high frequency, and the shape of the part in high frequency almost does not change.
e results of the analysis show that the noise in shock signals is ltered e ectively.e YPT therefore shows good performance in the analysis of the shock signal in engineering, which can provide help for further extracting features for shock signals in pattern recognition and fault diagnosis.

Summary and Outlook
Because the Morlet wavelet does not satisfy the permissibility condition of the wavelet transform, it has no inverse 12 Shock and Vibration wavelet transform.On the contrary, the Morlet wavelet has good time-frequency characteristics.In this paper, it was easy to adaptively adjust the centre and radius of the frequency window for the Morlet wavelet after its improvement.e essence of the wavelet transform is that it is a type of filtering.For this paper, the essence of transforming a signal into the time and frequency domains is also a filter.e YPT in this paper can separate noise from the signal, the essential aspect of which is that noise and useful signals have different extreme frequencies.On the principle of extreme frequency, the original signal is decomposed into different components using the YPT.
In this paper, the noisy signal and the shock signals in engineering could also be complex or real.e extremum frequency is determined by the real part of the signal when analysing a complex signal.ere will be two ways for us to conduct research if further analysing a multidimensional signal.For the first, we can analyse each dimension of the multidimensional signal using the YPT alone.For the other, we can extend the dimension of the improved Morlet wavelet so that it has the ability to analyse multidimensional signals.
For the latter method, how to determine the centre and radius of the Morlet wavelet in the frequency domain needs to be researched more deeply.Matching a multidimensional signal to the Morlet wavelet in phase is worthy of study.

4. 1 .
Two Adjacent Morlet Wavelets.In this section, we analyse two adjacent Morlet wavelets improved by us.e centres of the frequency windows separately are f k−1 and f k , and the radii are Δc k−1 and Δc k , respectively, where k ≥ 2 and k ∈ N. We made their radii just connect with each other in the frequency domain, the situations of which are shown in the following two cases and seen in Figures

Figure 2 :
Figure 2: e rst situation of two adjacent Morlet wavelets.

Figure 3 :
Figure 3: e second situation of two adjacent Morlet wavelets.

Figure 4 :
Figure 4: e chart for analysing a signal.
. e signals in Figures 10 and 11 and those in Figure 9 have similar properties.For the convenience of expression in this paper, the signals whose serial numbers are opposite to each other are, respectively, named "positive

Figure 5 :
Figure 5: e ow chart for reconstructing a signal.

8
, they have properties similar to those of the new ltered signals with serial numbers 2 and −2.What follows is to transform the new ltered signals in the frequency domain into the time domain, and we can then continue to lter them in the time domain by probability ltering.e new ltered signals in three-dimensional space are plotted in the time domain in Figure 15.Compared with the signals in Figure 7, each of the signals in Figure

Figure 9 :Figure 10 :Figure 11 :Figure 12 :
Figure 9: e real part of the signals in the frequency domain.

Figure 13 :Figure 14 :
Figure 13: e real part of the new positive ltered signal in the frequency domain.

Figure 15 :Figure 16 :
Figure 15: e real part of the new ltered signals in the time domain.

Figure 17 :Figure 18
Figure 17: e environment of the experiment.(a) e positions prone to noise; (b) the place for the experiment; (c) signal collector; (d) the main measuring point; (e) the rst kind of cracks; (f ) the second kind of cracks.
that are produced by the free vibrations of each part of the gun automatic mechanism.e actual e ective time period of original signals was intercepted to be analysed.e original signals and the signals ltered by the YPTare plotted in Figure 18.Figures 18(a)-18(c) are the original signals, and Figures 18(d)-18(f) are the signals that correspond to Figures

Figure 19
Figure 19: e original signals and the signals ltered by the YPT in the frequency domain.

Table 1 :
e parameters of the sensor.