Noise Reduction for Modal Parameter Identification of the Measured FRFsUsing theModal Peak-BasedHankel-SVDMethod

(e measured frequency response functions (FRFs) in the modal test are usually contaminated with noise that significantly affects the modal parameter identification. In this paper, a modal peak-based Hankel-SVD (MPHSVD) method is proposed to eliminate the noise contaminated in the measured FRFs in order to improve the accuracy of the identification of modal parameters.(is method is divided into four steps. Firstly, themeasured FRF signal is transferred to the impulse response function (IRF), and the Hankel-SVDmethod that works better in the time domain rather than in the frequency domain is further applied for the decomposition of component signals. Secondly, the iteration of the component signal accumulation is conducted to select the component signals that cover the concerned modal features, but some component signals of the residue noise may also be selected.(irdly, another iteration considering the narrow frequency bands near the modal peak frequencies is conducted to further eliminate the residue noise and get the noise-reduced FRF signal. Finally, themodal identificationmethod is conducted on the noise-reduced FRF to extract themodal parameters. A simulation of the FRF of a flat plate artificially contaminatedwith the randomGaussian noise and the randomharmonic noise is implemented to verify the proposedmethod. Afterwards, amodal test of a flat plate under the high-temperature conditionwas undertaken using scanning laser Doppler vibrometry (SLDV). (e noise reduction and modal parameter identification were exploited to the measured FRFs. Results show that the reconstructed FRFs retained all of the modal features we concerned about after the noise elimination, and the modal parameters are precisely identified. It demonstrates the superiority and effectiveness of the approach.


Introduction
In modal testing, the measured FRF signals are usually contaminated with noise [1,2], and the noise will interfere with the accurate extraction of modal features. To eliminate the noise and purify the measured FRF signals, the noise reduction method should be exploited. In the past decades, many noise reduction methods such as the traditional low/ high-pass filter [3], singular value decomposition (SVD) [4,5], minimum mean square error (MMSE) [6,7], Wiener filter [8,9], wavelet transform (WT) [10,11], empirical mode decomposition (EMD) [12,13], independent component analysis (ICA) [14,15], and deep recurrent neural network (DRNN) [16,17] have been developed. Among these methods, the SVD-related filters draw the most noticeable attention for their convenience, simplicity, and nonparametric properties and have been widely used in the noise reduction of various signals in engineering.
Generally speaking, the noised high-dimensional data, like the image signal [18] and multisensor signals [19], can be directly handled by the SVD filters without any preprocessing. But, difficulty exists for the low-dimensional data such as in the case of modal testing, in which only a few sensors are used to measure the frequency response functions (FRFs). In order to purify the measured FRF signals, the preprocessing techniques are conducted to improve their dimension before the SVD filter is subsequently carried out. e commonly used technique is the Hankel-SVD filter [20] that constructs a Hankel matrix from the separated single FRF signal to improve its dimension and then uses the SVD to purify the constructed Hankel matrix to reduce the noise of the correlated single FRF signal. e traditional Hankel-SVD filter mainly contains two steps [21]. e first step is the decomposition of the noise signal into the component signals and their correlated singular values using the Hankel-SVD method. e next is to select the useful component signals and the corresponding singular values to form the noise-reduced signal. In both steps, the decisive accuracy for the noise removal in the single FRF signal is how to select the useful singular values or the component signals.
In recent years, Yang and Tse [22] proposed a singular entropy selection method for singular value selection, which is to select the saturation point of the singular entropy calculated from the singular value sequences. Fan et al. [23] proposed a cluster method to select proper singular values, which is to select the big difference point of the standard deviation sequence calculated from the singular value clusters. Zhao and Ye [24] argued that the singular entropy is just like an inverse transform of the singular value sequence and proposed a difference spectrum of the singular value (DSSV) selection method to select the biggest point of the difference spectrum of the singular values. Bao et al. [25] introduced the model order indicator (MOI) to the singular value selection of the IRF signal for its noise reduction, which is to select the biggest point of the MOI spectrum calculated from singular values, and the formula of MOI is like an extension of DSSV. Although the above methods achieved good results in many cases, some drawbacks still exist. e above methods are based on the assumption that there will be a huge difference of the singular values produced by the noise and by the real signal. e big difference point between them can be used to select the number of singular values for denoising the signal. However, this assumption does not make any sense when dealing with the weak features in the signal because the big singular values do not represent the weak features sometimes and they do not contain a clear physical meaning.
Instead of using the biggest different point of the singular values to select the number of useful components, some new feature selection methods are based on the decomposed component signals from Hankel-SVD. Qiao and Pan [26] proposed a correlation coefficient selection method, which is to select the biggest difference point of the correlation coefficient spectrum calculated from the component signals and the original signal. e method shows a better result than the DSSV selection method. Zhao and Jia [27] proposed a periodic modulation intensity (PMI) method to select and enhance the fault feature of the rolling machine signal, which is to select the component signals whose impulse energy is bigger than noise energy for noise reduction. e result shows the method well enhances the impulse feature. Recently, Zhao and Ye [28] proved that the SVD has a feature to recover certain frequency components, and Cheng et al. [29] proposed a resonance enhancement SVD for feature extraction by adding the excitation sinusoidal signals with concerned frequencies to the original signal to form the enhanced signal. en, a rectangle Hankel matrix with a certain dimension is constructed from the enhanced signal.
e Hankel matrix is further decomposed into component signals using the SVD technique, and the first component signal, together with the excitation sinusoidal signals subtracted, is extracted to realize the feature extraction.
In this paper, a novel modal peak-based Hankel-SVD (MPHSVD) method is proposed to eliminate the noise and retain the whole concerned modal features of the FRF signals whose modal peak frequencies and modal peak amplitudes can be selected by observation. e MPHSVD is an iteration Hankel-SVD filter, which contains the accumulation iteration and reselection iteration. e accumulation is to get the component signals of the transformed IRF signal that covers the whole concerned modal features, but these component signals will also cover some noise. e reselection is based on the accumulation result, which aims to separate the component signals which belong to the modal modes from the component signals of the noise and add them up to complete noise reduction thoroughly. Afterwards, the modal identification based on the rational fraction polynomial (RFP) method is conducted to the noise-reduced FRF signal to extract the modal parameters. e simulated and two experimental FRF signals of a flat plate with heavy noise are utilized to verify the approach. e results show the validity and superiority of the proposed method. e rest of the paper is divided into four sections. In Section 2, the modal identification method based on the MPHSVD technique for noise reduction is represented. In Section 3, a simulated case is studied. In Section 4, the experimental cases are further used to verify the proposed method. e conclusion is drawn in the last section.

Modal Parameter Identification of the Measured FRFs with Noise Reduction by MPHSVD
In this section, the modal identification based on MPHSVD is described in detail, and the MPHSVD method is the Hankel-SVD filter with the modal peak-based component signal selection method, which is used for eliminating the noise of the single FRF signal.

2.1.
Hankel-SVD of the Transformed IRF Signal. As the transformed impulse response function (IRF) signal is better in the noise reduction than the FRF signal, the FRF signal H(ω) is therefore transformed to the IRF signal by inverse Fourier transform: where h(t) represents the transformed IRF signal and iFT represents the inverse Fourier transform. Besides, the inverse Fourier transform is conducted in the fast Fourier transform framework, so the transformed IRF signal correlated with the N/2 length long FRF signal H(ω) contains N points.
After the transformed IRF signal is obtained, the Hankel-SVD is conducted to decompose the IRF signal into component signals. For the length N transformed IRF signal h(t), its Hankel matrix A is where m is the row number of the Hankel matrix A, n is the column number, and m + n − 1 � N. Here, m is set as floor(N/2) to make the Hankel matrix nearly square, where the symbol "floor" means rounding down to the nearest integer. e reason of setting the Hankel matrix nearly square is that a nearly square Hankel matrix can achieve the maximum noise reduction quality of the Hankel-SVD filter.
en, matrix A can be decomposed by the SVD technique: where U is the left orthogonal matrix, u i is the column vector of U, V is the right orthogonal matrix, v i is the column vector of V, S is the singular value matrix, and σ i is the i th singular value. e dimension of S is q, and q is equal to min(m, n). en, by expanding equation (4), we get the following equation: where A is decomposed to the summation of q component matrices and A i is the i th component matrix of the component matrix set. e dimension of A i is equal to the dimension of A, which is m × n. Besides, we assume the component matrix still contain the Hankel structure, and matrix A i is shown in the following equation: where s i a,b represents the element of the a th row and the b th column in matrix A i and A i contains mn elements. Due to the Hankel structure of A i , the values of the elements in the antidiagonal line are actually the same, so the elements in each antidiagonal line can be represented by one element. us, the whole mn elements in A i can be represented by N elements because there are totally N antidiagonal lines. Besides, N is equal to m + n − 1, and h i j represents the element of the j th antidiagonal line of matrix A i . en, we define the first row of A i as L i and define the last column without the first element of A i as R i . e elements of L i and R i can form a component signal h i (t), which is correlated with its singular value σ i : en, because Hankel matrix A can be divided into the summation of the component matrices, A � A 1 + A 2 + · · · + A q , the first row L and the final column R of matrix A can also be divided to the sum of the first rows and the last columns of the component matrices according to the properties of matrix addition: Finally, the IRF signal h is decomposed to the component signals formed by L i and R i because the signal h can be represented by the combination of L and R: e above derivation is based on the assumption that component matrices still contain the Hankel structure after SVD, but the decomposition result of SVD does not actually follow the assumption. us, we usually average the antidiagonal lines of the component matrices such that they contain the Hankel structure and get the component signals. e antidiagonal averaging of the component matrix A i is shown in the following equation: where the averaging of the l th antidiagonal line means adding up the elements whose coordinate summation in A i is equal to l + 1 and dividing the total element number d l to get the element h i l . e red lines in the equation mean averaging, Shock and Vibration 3 and the equation means averaging each antidiagonal line of A i to get the elements of h i .

Selection of the Number of the Accumulated Component
Signals. Based on the Hankel-SVD technique, the frequencies and the amplitudes at modal peaks of the FRF signal can be decomposed to the component signals [28]. e component signals are arranged in orders corresponding to the singular values sorted from the largest to the smallest. Generally, one modal frequency mode can be decomposed into one or several pairs of adjacent component signals, and one pair of adjacent component signals correlated with the largest pair of singular values of this mode is likely dominated. For the experimental FRF signal H(ω), the concerned number of modal peaks is set as k.
en, the concerned modal peak frequencies, f � f 1 f 2 · · · f k , and their amplitudes, R � R 1 R 2 · · · R k , can be selected. Considering all the decomposed component signals of the IRF signal, the iteration of the component signal accumulation is conducted. For the i th iteration, the i th component signal h i is added to the (i − 1) th accumulated signal h i−1 acu to form the i th accumulated signal h i acu : Afterwards, the energy spectrum of the accumulated signal h i acu is calculated by Fourier transform: where E i represents the energy spectrum of the accumulated signal h i acu . en, the differences between modal peak amplitudes R and the amplitudes of energy spectrum E i at modal peak frequencies f are calculated: where the symbol "abs" represents the absolute value, the symbol "./" means the vector point division, D i is the i th difference vector, and D i l is the difference between the l th modal peak amplitude R l and the energy spectrum E i at the modal peak frequency f l .
Finally, D i is judged by the threshold T. If the arbitrary component of vector D of the i th iteration is below or equal to the threshold T, the iteration stops, and the rest of the component signals are treated as noise. If not, the next iteration will be conducted. After the accumulation ends, the whole modal modes will recover, and the number of the selected component signals is recorded as m. However, some residue noise may also be recovered in the component signal accumulation. To eliminate the residue noise, and accomplish the noise elimination thoroughly, a further reselection method will be conducted.

Component Signal Reselection for Eliminating the Residue
Noise.
e peak frequencies of the component signals correlated with the modes are near the modal peak frequencies, while the peak frequencies of the component signals correlated with the residue noise will be away from the modal peak frequencies in some degrees. us, if some narrow frequency bands are set near the peak frequencies, the residue noises can be eliminated by reselecting the component signals whose peak frequencies are in the bands from the accumulation result for noise reduction.
For there are k cared modal frequencies [f 1 , f 2 , . . . , f k ] in H, and k frequency bands with the bandwidth d are set as the passband P: Since there are m component signals extracted from the component signal accumulation, another iteration with the max iteration number m is conducted to reduce the residual noise, and the j th iteration is shown below.
Firstly, for the j th component signal h j extracted from the accumulation, its correlated energy spectrum is calculated by the Fourier transform, which is shown in the following equation: where E j c is the energy spectrum of the j th component signal h j and FT represents the Fourier transform. en, the frequency of the peak amplitude in the energy spectrum E j c is extracted: where f j max is the frequency of the max amplitude of E j c . If the frequency f j max is in the passband P, the component signal is added to the second accumulated signal: If f j max is not in the band, the component signal h j is abandoned, and the next iteration is conducted. When the iteration turn is m, the iteration ends, and the second accumulated signal h acu2 is signed as the noise-reduced IRF signal h. e noise-reduced FRF H is transformed from it:

Modal Parameter Identification with the Noise-Reduced FRF.
After the noise reduction, the noise-reduced FRF signal H is acquired, and the modal parameters are subsequently extracted by the modal identification. e rational fraction polynomial (RFP) method that is widely used in the modal test package is selected here. One benefit of using the RFP method is that it is suitable for a single FRF to extract modal parameters. Since the method is well developed, and the main formulas of the method are described thoroughly in reference [30], more detailed description of the method is not given here.
In summary, the framework of the modal identification by MPHSVD is given in Figure 1. It is divided into four steps as follows: Step 1: transform the noised FRF signal to noised IRF signal by IFFT and decompose it to n component signals and the corresponding n singular values by Hankel-SVD Step 2: accumulate the component signals in the order of the corresponding singular values from the largest to the smallest to get the first m component signals that cover the modal modes and some residue noise Step 3: reselect q component signals of the modal modes from the first m component signals and add them up to get the noise-reduced IRF signal and the corresponding FRF signal Step 4: identify the modal parameters from the noisereduced FRF signal by modal identification

e Simulated FRF Signal with the Added Noise.
To verify the proposed method, a simulation of a steel flat plate is carried out using the finite element (FE) analysis. An FRF signal of one point at the middle left of the plate is simulated. e calculated frequency band is [0∼3200] Hz with an interval of 0.5 Hz. erefore, the FRF signal contains 6400 points with a maximum frequency of 3200 Hz. e corresponding impulse response function (IRF) is attained by the inverse Fourier transform from the FRF signal and contains 12800 points. e simulated FRF signal and its corresponding IRF signal are shown in Figure 2.
In Figure 2(a), we can clearly see seven frequencies at the modal peaks which are 64 Hz, 399 Hz, 428 Hz, 1119 Hz, 1334 Hz, 2195 Hz, and 2381 Hz, respectively. Figure 2(b) is the transformed IRF signal. It can be seen that the amplitude of the transformed IRF signal does not ideally decrease to zero since the IRF signal is taken from the real part of the discrete inverse Fourier transform of the simulated FRF signal.
is will cause some digital errors and make the energy spectrum of the transformed IRF signal slightly different from the energy spectrum of the original FRF signal. Besides, the proposed method is mainly conducted on the transformed IRF signal, so the frequencies and the amplitudes at the modal peaks mentioned in this paper are selected from the energy spectrum of the transformed IRF signal instead of the energy spectrum of the original FRF signal for accuracy. en, to simulate the IRF signal contaminated with noise, the additive noise is added to the IRF signal by the following equation: (20) where h noise (t) is the noised IRF signal, h(t) is the noise-free IRF signal, n(t) is the random Gaussian noise, A n is the amplitude of the random noise, which is set as two point five percent of the max amplitude of the noise IRF signal h(t), and A ni cos(2πf ni t) is the random harmonic noise; its frequency is obtained from the frequency band [100,400], [500,1000], [1500,2000], and [2500,3000] randomly. Its amplitude is randomly obtained from the amplitude band [0,0.05A max ], where A max is the max amplitude of the noisefree IRF signal h(t). After the noise signals are added to the noise-free IRF signal, the SNR between the original IRF signal and the added noise is 0.4691, and the noise IRF signal with its correlated noised FRF signal is shown in Figure 3.
Obviously, the FRF signal in Figure 3(a) is polluted heavily by the noise. e energy spectrum of the noised FRF signal severely distorts from the noise-free one, and the phase spectrum is totally drown into the noise, where many of the phase changes are hard to be seen directly. From Figure 3(b), it can be seen that the noise pollutes most parts of the IRF signal, which makes the noised IRF signal deeply distort from the simulated IRF signal. us, the qualities of the FRF and IRF signals are seriously damaged by the noise. en, the real modal parameters are extracted from the original noise-free FRF signal using the RFP method, and the noised modal parameters are extracted from the noised FRF signal using the same curve fitting frequency intervals of the noise-free FRF signal. e curve fittings are plotted in Figure 4, and the identified modal parameters are listed in Table 1. Figures 4(a) and 4(b) show the curve fitting of the original FRF signal, and the fitted curve matches the original FRF signal well. e curve fitting of the noised FRF signal is shown in Figures 4(c) and 4(d). It can be found that there are some differences in the modal peak frequencies 1 and 5 between the fitted curve and the noised FRF in the energy spectrum. Besides, by comparing the fitted curve and the noise-free FRF signal, there is some difference between the fitted curve and the noise-free FRF near modes 2 and 3.
us, the fitted curve is inaccurate. From Table 1, the modal parameters of the noised FRF and real modal parameters extracted from the original noise-free FRF signal are given. Besides, the differences between them are also calculated. e frequency differences are the absolute differences between the modal frequencies of the two FRF signals. e damping factor differences are calculated by the following equation: where D damp i is the damping factor difference of the i th mode, "abs" means the absolute value, ξ noisefree i is the damping factor from the noise-free FRF signal, and ξ noised i is the damping factor of the noised FRF signal. Similarly, there is Step 1 Step 2 Step 4 Step 3 Hankel-SVD IFFT Noised IRF signal is decomposed to n component signals  nearly a 8 Hz difference in mode 1 between the extracted modal frequency and the real modal frequency. Besides, there is a 6 Hz difference in mode 5 between the extracted and the real modal frequency. e extracted damping factor of mode 1 is negative, and the difference between damping factors of mode 1 extracted from the noise-free and noised FRFs is 11 times bigger than the real damping factor of mode 1. It is noticed that the damping factor difference of mode 1 is far beyond the damping factor differences of the other modes. It is because that the mode 1 has the lowest energy of all concerned modes, and it is heavily affected by the noise. Besides, due to the inaccurate curve fitting of mode 2, the damping factor extracted from the noised FRF signal is nearly 2.6 times bigger than the factor of the noise-free FRF signal. en, the analysis of the noised FRF signal shows that the quality of the FRF signal is  Shock and Vibration 7 damaged by the noise seriously, and the extracted modal parameters are inaccurate.

e Noise Reduction of the Simulated FRF Signal Using the MPHSVD Method.
To eliminate the noise contaminated in the FRF signal, the corresponding noised IRF signal in Figure 3(b) is taken into consideration. Firstly, a nearly square Hankel matrix with 6400 rows and 6401 columns is constructed from the noised IRF signal. e Hankel matrix is subsequently decomposed into 6400 component signals using the SVD technique. e seven modal frequencies and their amplitudes are selected for noise elimination. e amplitude difference threshold T is set as 0.15, and the iteration of the component signal accumulation is conducted to recover the selected modal peak amplitudes and modal peak frequencies.
e iteration stops at the 80 th iteration turn, and the first 80 component signals are accumulated to form the accumulated signal. en, the accumulated signal is transformed to the frequency domain to acquire the correlated FRF signal, and the overlay between the FRF signal from the accumulated signal and the FRF signal from the transformed noise-free IRF signal is given in Figure 5. Obviously, seven modal frequencies and their amplitudes are recovered from Figure 5(a), which is marked in red circles in the figure. Clearly, the random Gaussian noise is mostly eliminated, but many random harmonic noises are simultaneously recovered after the accumulation, which makes it hard to distinguish which peaks are correlated with the real modal modes. e phase spectrum after accumulation is shown in Figure 5(b). Similarly, the phase changes correlated with the modal modes are recovered, and the changes of the harmonic noises are also recovered, which makes the phase spectrum still distort from the noise-free one a lot. en, from Figure 5, the accumulation removes most parts of the random Gaussian noise, but it recovers most parts of the random harmonic noises with the recovery of the modal modes. To eliminate the residue harmonic noises and complete the noise reduction thoroughly, the further component signal reselection is conducted.
Because seven selected modal peaks are supposed to be kept after noise recovery, seven frequency bands with the bandwidth 35 Hz near the seven modal peak frequencies are set as the threshold for the further noise removal. e component signals, that the frequency of the biggest peak amplitude is in the set frequency bands, are reselected for the further noise reduction. en, for the 80 component signals selected from accumulation, the component signal reselection is conducted for them to separate the component signals correlated with the modal modes for recovery and abandon the component signals of the random harmonic noises. Finally, 20 component signals are reselected, and they are added up to form the noise-reduced IRF signal. e overlays of the noise-reduced IRF signal against the noised IRF signal and the noise-reduced IRF signal against the noise-free IRF signal are given in Figures 6(a)  and 6(b).
It can be found that most parts of the noise are removed, and there is mostly no difference between the noise-reduced and the noise-free IRF signals, except the part of the red circle in the figure. en, a noise removal rate is calculated by the following equation where h denoised is the noise-reduced IRF signal and NR is the noise removal rate. e calculation result from Figure 6(b) is 0.9365, which indicates that 93.65% of the noise is removed. en, the noise-reduced IRF signal is transformed to the noise-reduced FRF signal, and the comparison between the noise-reduced FRF signal and the FRF signal from the noisefree IRF signal is shown in Figure 7.
e noise-reduced energy spectrum comparison is shown in Figure 7(a). It is obvious that the random harmonic noises are eliminated, where the seven modal peaks are the only peaks that can be seen directly. Besides, the noise-reduced spectrum fits the noise-free one well, which shows the validity of the noise elimination. Figure 7(b) gives the noise-reduced phase spectrum of the reselected 22 component signals, the phase changes correlated with the modal features are recovered, and they fit the noise-free one well. From the analysis of the noise reduction result, the MPHSVD filter successively eliminates the random Gaussian noise and the random harmonic noises thoroughly. Besides, the final noise reduction result fits the noise-free one well, which shows the validity of the proposed method.      Figure 9, with the overlay of the noise reduction result of the pass bandwidth 35 Hz. It can be found that the results of the reselection of bandwidths 5 Hz, 10 Hz, 20 Hz, 35 Hz, 45 Hz, and 60 Hz are the same, which means the reselection iteration with the arbitrary bandwidth between 60 Hz and 5 Hz for the noised IRF here will achieve the same noise reduction quality. en, when the bandwidth improves to 80 Hz, one harmonic noise at 104.5 Hz is recovered. When the bandwidth decreases to 2 Hz, the recovery of mode 1 at 65 Hz and mode 7 at 2381 Hz is incomplete. us, the reselection with too big bandpass width will recover some additional noise, and the reselection with too small bandpass width will lose some modal features. en, since the least frequency interval between the modal peaks and the harmonic noise in the real modal test is unknown, it is difficult to decide the biggest frequency bandwidth that the residue noise will not be recovered. So, the frequency pass bandwidth of the reselection should be small to avoid the recovery of the noise, but not too small that the modal features are abandoned. One frequency bandwidth like 5 Hz or 10 Hz is enough for the recovery of the whole modal modes and the removal of the noise.

Modal Parameter Identification with the Noise-Reduced FRF.
After the MPHSVD is implemented, the curve fitting with the rational fraction polynomial (RFP) method is conducted to extract the modal parameters. Modal modes are identified piecewise, and the modal frequencies with the correlated damping factors are extracted for each mode. e curving fitting of the noise-reduced FRF signal is plotted in Figure 10, and the identified modal parameters are listed in Table 2. From Figure 10, the fitted curve fits the noise-reduced FRF signal well, which shows the validity of curve  fitting. From Table 2, it can be seen that there is little difference between identified modal frequencies and the original modal frequencies, which is more precise than the modal frequencies extracted from the noised FRF signal.
en, compared with the modal damping factors from the noised FRF signal, the damping factor of mode 1 is positive, and the difference between the damping factors of mode 1 significantly decreases from 1100% to 60%. In practical applications, the first mode is usually considered as the fundamental frequency, and the 60% difference of the damping factor still looks high. e better way to deal with it maybe is to pick up another measurement location in order to increase the energy of mode 1. e damping factor of mode 2 from the noise-reduced FRF signal is much closer to the correlated damping factor of the noise-free FRF signal, where the difference of damping factors of mode 2 decreases from 160% to 32%. e other damping factors extracted from the noise-reduced FRF signal are also closer to the original damping factors extracted from the noisefree FRF signal compared with those extracted from the noised FRF signal. us, the analysis shows that the curve fitting of the noise-reduced FRF signal is precise, and the modal parameters extracted from the noise-reduced FRF signal are more precise than those from the noised FRF signal, which further proves that the noise reduction is successful.
In order to further investigate the effectiveness of the MPHSVD to deal with various degrees of the noise, four degrees of the random Gaussian noise, listed as 3%, 5%, 7%, and 9% ratios between the amplitudes of the added Gaussian noise and the max amplitude of the noise-free IRF signal, respectively, are added to the simulated noise-free IRF signal. Five simulations of each ratio are conducted, and their corresponding mean SNRs are 7.68 dB, 3.25 dB, 0.31 dB, and −1.89 dB, respectively. en, the noise removal rates are calculated by equation (22), and the results with the corresponding noise reduction successful rate are shown in Table 3. It shows that as the noise amplitude increases, the noise removal rate decreases and the noise reduction successful rate also decreases. e overlays of the noise-reduced FRF signals against the noised FRF signals of the simulations given in italics in Table 3 are conducted, and the result is shown in Figure 11. Figures 11(a) and 11(b) show that the noise removal rates can reach over 90% when SNRs of the noising FRFs are 7.68 dB and 3.25 dB, respectively. e noise removal rate decreases when the SNR becomes worse. In Figure 11(d), the FRF in green color contaminated with the heaviest noise with an SNR of −1.89 dB can hardly distinguish the vibration modes. After reducing the noise, most of concerned modes in black color except for mode 5 can be clearly seen although the noise-reduced FRF looks a bit rougher than the noise-free FRF.

Experiment Configuration.
To further verify the proposed method, an experimental test is undertaken. e force  hammer typed as PCB086C03, manufactured by PCB, is used for excitation. e UM-40-AB-HF electromagnetic heater and the heating coil, produced by UIHM, are used for heating the plate. e 3D scanning laser Doppler vibrometers (SLDVs) typed with PSV-400-3D for response measurement are manufactured by Polytec. e experimental setup is shown in Figure 12(a), and the measuring points are shown in Figure 12(b).
In the modal test, the flat plate is firstly grounded at the one end to a metal baffle tightened with four bolts, and the baffle is connected to the test bench to monitor the rigid support. en, the heating coil is connected to the bench behind the plate and connected to the electromagnetic heater for heating the plate. e force hammer is then connected to the computer for excitation and measuring the input force. Subsequently, the SLDVs are placed in front of the test plate for the noncontact measurement. e measuring points are evenly distributed on the plate, and totally 45 points are selected. e test starts when the temperature of the plate is heated and stabled at the 200 centigrade degrees. e analysis frequency range is 6400 Hz, and it takes 2 seconds once for each measuring point. After testing, the FRF signals are obtained for these 45 measuring points.

e Analysis of the Noised FRF Signals.
e measured FRF signals of Point 40 and Point 1, marked as the red circle in Figure 12(b), at 200 centigrade degrees are taken as the example. e overlay of the energy and the phase spectra of the FRF signals of Points 1 and 40 in the heated temperature is plotted in Figure 13.
It can be found from Figure 13(a) that the energy spectrum of Point 1 is clearer, but the spectrum of Point 40 contains heavy noise and cannot be directly used for modal parameter identification. One contributing factor is thought that the speckle noise of the FRF signal is enhanced by the high temperatures, since the speckle noise is usually dominant in the measurement of the SLDV [31]. e ten peak amplitudes and frequencies in the red and blue circles are the ten modal modes we concerned about. ese peak frequencies are selected from the clear energy spectrum of point 1 by observation. e modal peaks in the blue circles can only be observed in the energy spectrum of Point 1, and the correlated modal peaks of Point 40 are totally drown into the noise, which means they cannot be recovered by the MPHSVD filter. us, the modal peak frequencies and amplitudes of the red circles in the light blue energy spectrum are the concerned modal features of Point 40 and those of the whole red and blue circles in the black curve are the concerned modal features of Point 1. e other peaks in the signal are thought to be induced by the added mass of the vibrator, and they are not selected for recovery Figure 13(b) gives the phase spectra of the two FRF signals. Similarly, the phase spectrum of Point 40 contains heavy noise, and some of the modal phase changes cannot be seen clearly. Finally, the overlay between the denoised FRF signals and the FRF signals from the transformed IRF signals is given in Figure 14.

e Noise Reduction of the Measured FRF
From Figures 14(a) and 14(b), the noise reduction result of Point 1 can be seen. It can be seen that the concerned modes are retained, and the spectra of the energy and phase are much smoother and clearer after noise reduction, which shows the validity of the noise reduction method. Figures 14(c) and 14(d) give the noise reduction result of Point 40. It shows that the spectra are much smoother than the original spectra of the transformed IRF signal, and most concerned modal modes are recovered completely. However, the amplitudes of the modal modes at 1127 and 1149 Hz are not recovered completely, and the modal mode at 1014 Hz is nearly eliminated after the noise reduction, which is marked in the red circles in Figure 14(a). e main reason is that these modes are contaminated with the noise heavily, and their correlated component signals are aliased with the noise. e modes at 1014, 1127 Hz, and 1149 Hz with their correlated decomposed component signals are shown in Figures 15 and 16. Figure 15(a) gives the energy spectra of modal mode at 1127 and 1149 Hz, and it is obvious that these modal modes are contaminated with the noise heavily, since their amplitudes fluctuate in the whole observed frequency range. Figure 15(

Modal Identification of the Noise-Reduced Measured FRF
Signal. After the MPHSVD is done, the RFP method is conducted to extract the modal parameters of the two denoised FRF signals. e curve fitting result and the extracted modal parameters of the concerned modes are shown in Figure 17 and Table 4. Figure 17 gives the curve fitting result of the two noise-reduced FRF signals. It can be found that the modal modes are fitted well, and the fitted curves have a good coincidence with the denoised FRF signals. Besides, the mode at 1014 Hz of Point 40 is nearly eliminated since it is contaminated with the noise heavily, and it is not fitted to extract the modal parameters. Table 4 gives the extracted modal parameters of the two denoised FRF signals and the differences between the modal parameters. e modal frequency difference is the absolute difference between the modal frequencies of the two denoised FRF signals. e modal damping factor difference is calculated by equation (21). e extracted modal frequencies of the two signals are nearly the same, and the max difference between the extracted modal frequencies of the two denoised FRF signals is 1.21 Hz, which shows the validity of extracted modal frequencies. For the damping factors, the max difference between them is the difference at mode 5, which is 68%. e main reason is thought that the mode 5 of Point 40 loses a large part of its amplitude in the noise elimination, and the correlated damping factor is estimated to be higher than the real one. en, for the other damping factors, the differences between them are not big, and the max difference is 32%, which shows that the extracted damping factors are also validated. en, the analysis shows that the curve fitting of the two denoised FRF signals is precise, and the extracted modal parameters match each other, which further proves the validity of the noise reduction.

Conclusion
In this paper, a novel noise reduction method based on the modal peak-based Hankel-SVD technique (MPHSVD) is proposed, which aims to acquire the precise FRF signal from the noised FRF signal for the modal identification. e proposed method contains one accumulation iteration to select the component signals which cover the whole modal modes and some residue noise and one reselection iteration to separate the component signals of the modal modes from the residue noise for completing the noise reduction. e advantage of the method is that the component signals correlated with the modal modes can be determined by the two iterations. e simulated FRF of a plate with the added heavy noise is taken as an example to show how the denoising method works. e noise removal rate reaches above 93%, and the further modal analysis shows that the modal frequencies and modal damping factors are extracted precisely. Besides, the detailed analysis of the correlation between the component signals and the modal modes and effect of different bandwidths for the noise reduction are also investigated. e simulation analysis of FRFs contaminated with various degrees of Gaussian noise shows that the noise reductions are fully successful for FRFs with SNRs of 7.68 dB and 3.25 dB, respectively. e application of two experimental FRF signals of a flat plate in the high-temperature environment also shows the superiority and effectiveness of the proposed method. However, some amplitudes of the modal modes contaminated with too much noise decrease after the noise elimination. Nevertheless, the method shows the practical potential to engineering application in the modal test.

Data Availability
Data are available on request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.