Next Article in Journal
Signal Source Localization with Long-Term Observations in Distributed Angle-Only Sensors
Previous Article in Journal
Intra Prediction Method for Depth Video Coding by Block Clustering through Deep Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Single-Channel Blind Signal Separation of the MHD Linear Vibration Sensor Based on Singular Spectrum Analysis and Fast Independent Component Analysis

1
School of Instrument Science and Opto-Electronics Engineering, Hefei University of Technology, Hefei 230009, China
2
State Key Laboratory of Precision Measuring Technology and Instruments, Tianjin University, Tianjin 300072, China
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(24), 9657; https://doi.org/10.3390/s22249657
Submission received: 29 October 2022 / Revised: 4 December 2022 / Accepted: 8 December 2022 / Published: 9 December 2022
(This article belongs to the Section Physical Sensors)

Abstract

:
An MHD vibration sensor, as a new type of sensor used for vibration measurements, meets the technical requirements for the low-noisy measurement of acceleration, velocity, and micro-vibration in spacecraft during their development, launch, and orbit operations. A linear vibration sensor with a runway type based on MHD was independently developed by a laboratory. In a practical test, its output signal was mixed with a large amount of noise, in which the continuous narrowband interference was particularly prominent, resulting in the inability to efficiently carry out the real-time detection of micro-vibration. Considering the high interference of narrowband noise in linear vibration signals, a single-channel blind signal separation method based on SSA and FastICA is proposed in this study, which provides a new strategy for linear vibration signals. Firstly, the singular spectrum of the linear vibration signal with noise was analyzed to suppress the narrowband interference in the collected signal. Then, a FastICA algorithm was used to separate the independent signal source. The experimental results show that the proposed method can effectively separate the useful linear vibration signals from the collected signals with low SNR, which is suitable for the separation of the MHD linear vibration sensor and other vibration measurement sensors. Compared with EEMD, VMD, and wavelet threshold denoising, the SNR of the separated signal is increased by 10 times on average. Through the verification of the actual acquisition of the linear vibration signal, this method has a good denoising effect.

1. Introduction

With the development of space technology, the field of space engineering demands increasingly high standards for the structural design of spacecraft, among which, the analysis of the mechanical environmental conditions of spacecraft is key to its overall design [1]. The mechanical environment of spacecraft refers to the environment of vibration, acceleration, noise, and microgravity that the spacecraft is subjected to during operation. Wide-band and low-noise vibration sensors can meet the needs of the wide-band, high-precision mechanical environment testing of spacecraft [2].
In this paper, a gallium–indium–tin alloy with good conductivity and no bias is used as the conductive fluid material for the MHD linear vibration sensors. A linear vibration sensor based on the MHD effect has the advantages of a high yield strength, good dynamic performance, and no mechanical wear phenomenon between internal firmware. It can be easily used in the scenarios of wide-band frequency and high-precision measurements, and can be used for effective in-orbit space vibration measurements. In the actual vibration measurement, the measurement system contains multiple error sources inside, and the weak broadband signals detected by the MHD linear vibration sensors are often mixed with multiple noise sources. At the same time, the lack of accuracy in the dynamic measurement process remains a key concern for engineering measurement personnel and instrument designers. We need to take effective and efficient measures to eliminate or reduce interference information, and thus improve the accuracy of the dynamic measurement system.
The requirement for the MHD linear vibration sensors in aerospace technology to have amplitudes within millivolts and vibration frequencies within kilohertz, coupled with the complex external environment [1], makes signal acquisition and extraction separation difficult. Therefore, the extraction and separation of weak signals output from the MHD linear vibration sensors is one of the key techniques used to perform micro-vibration measurements [3].
Blind source separation (BSS) is a blind signal extraction technology, which is often used to separate and extract independent sources in mixed signals [4,5,6]. Independent component analysis (ICA), as a blind source separation technology, is widely used to separate covered source signals from mixed signals [7]. An independent component analysis cannot directly solve the underdetermined blind source separation problem of the MHD line vibration sensor output signal presented in this paper. Single-channel signals can be converted into multi-channel signals by some methods, and then the converted multi-channel signals can be used for ICA [8]. For example, in 2016, Lu et al. proposed a new wavelet threshold function and its denoising application [9]. From the SNR comparison of the processed signals, using this method it can be seen that the SNR of the output signal is the largest; 1–2.7 times higher than that of the traditional threshold denoising method. In 2018, Ajay Kumar Maddirala and Rafi Ahamed Shaik proposed a single-channel EEG signal source separation method combining SSA and ICA [10]. Through the simulation of analog signals, it was verified that this method can decompose the hidden source in a single-channel EEG signal and has a high source separation efficiency. In 2019, Meng Zhou et al. proposed a photoacoustic sensing imaging denoising algorithm based on empirical mode decomposition (EMD) and ICA [11], which can be applied to most pulse signal denoising problems. In the same year, Tao Jin et al. proposed a novel adaptive integrated empirical mode decomposition (EEMD) method based on partial discharge [12] and applied it to the simulated pulse signal. This method has a good denoising effect, and the SNR improves by about 1.2 times. In 2020, Dariusz Mika et al. proposed a single-channel source separation algorithm based on ICA time frequency decomposition [13], which used an ICA method to separate single-channel source signals from single mixed signals. The main limitation of this method is the lack of generality. The choice of distance measure and clustering algorithm depends on the time frequency structure of the mixed signal. As the number of mixed signals increases, the quality of the signal of separated components decreases. In 2020, Liu et al. proposed the vibration signal analysis of water seal blasting based on wavelet threshold denoising and HHT transform [14]. In 2022, Yufei Zhou et al. proposed an improved single-channel BSS algorithm based on variational mode decomposition (VMD) [15]. This algorithm is applied to the simulated signal, which effectively avoids the endpoint effect and mode mixing, and makes decomposition more stable. The decomposition layers of different signals are determined by setting different spectral correlation coefficient thresholds, and the sparse frequency division of the received mixed signals is realized. In summary, modal aliasing occurs when the signals are decomposed by EMD [16,17,18]. EEMD can solve the pattern aliasing phenomenon in the EMD method [19,20,21,22], but it needs to increase the ensemble average to compensate for the added white noise. White noise and the reference signal are not easy to apply when a real-time vibration measurement is carried out by an MHD linear vibration sensor. Moreover, the EMD-like method is used to process the output signal of the MHD linear vibration sensor, with a long calculation time and large reconstruction error, so the processed signal is likely to have a low SNR. Therefore, the EMD-like method is not suitable for processing the output signal of an MHD linear vibration sensor.
In order to separate independent source signals from mixed signals with a low SNR, this paper proposes a method of MHD linear vibration signal separation based on SSA and FastICA. In this paper, SSA was first used to suppress the periodic narrowband interference and white noise in an acquired signal, and a noise-reduced mixed acceleration signal was obtained. Then, the noise-reduced signal was separated by a modified, fast, independent component analysis to obtain an unmixed signal. The unmixed signal was analyzed and identified to determine the useful and noisy components of the acquired magnetohydrodynamic sensor signal. Experimental analysis shows that this method can achieve the effective separation of independent line vibration signals with reasonable suppression of linear vibration signal admixture noise.
In the second section of this paper, the working principle of the MHD linear vibration sensor is described. In the third section, the blind source separation algorithm based on SSA and FastICA is proposed. In the fourth section, the experimental results are compared and verified.

2. The Working Principle of the MHD Linear Vibration Sensor

When the MHD linear vibration sensor is functioning, the conductive fluid moves in a circle in the magnetic field inside the channel. Figure 1(a) shows the structure diagram of the MHD linear sensor. The principle of the MHD linear vibration sensor is shown in Figure 1(b), where h is the height of the channel, d is the thickness of the channel, the entire working area is in a constant magnetic field B, F0 is the external inertia force on the measured carrier, and F1 is the Lorentz force on the conducting fluid in the magnetic field. Based on the principle of electromagnetic induction, the conductive fluid moves along the cutting magnetic inductive v under the superposition of inertia force and Lorentz force. According to the laws of electromagnetic induction, i.e., mechanical energy into electrical energy, due to the upper and lower wall conduction, an electric potential difference between the two walls is produced on the MHD channel between the internal and external wall, and the relative velocity and magnetic induction intensity forms a certain relationship between motional and electromotive force Uout, Uout = v·B (v is the relative velocity between the conductive fluid and the wall of the fluid channel).
According to the working principle of the MHD linear vibration sensor, the input–output model of the sensor was theoretically derived. The MHD explores the interaction between the motion field of conductive fluid and the magnetic field of the medium. It not only observes the influence of the electromagnetic field on the motion of conductive fluid, but also explores the reaction of fluid motion on the electromagnetic field. Without considering the influence of temperature change, according to the theory of MHD, the motion of Newtonian fluid and the basic size diagram of the sensor fluid channel are combined, as shown in Figure 2.
Without considering the influence of temperature change, according to the theory of MHD [23], by combining continuity equation and Navier–Stokes equation, we can obtain:
v t + ( v ) v = 1 ρ p f + μ 2 v 1 ρ J × B + a
where ρ is the fluid density, p f is the pressure, μ is the dynamic viscosity, F1 = −J × B (F1 is the Lorentz force [24]), a is acceleration in the inertial force F0 (F0 = m × a).
In the coordinate system, expanding along the X-axis, the main direction of movement is as follows:
v x t + ( v x v x x + v y v x y + v z v x z ) = 1 ρ p f x + μ ( 2 v x x 2 + 2 v x y 2 + 2 v x z 2 ) 1 ρ J × B + a
The following assumptions are made about the sensor during the model construction:
(1) Ignoring Secondary Flow Assumption:
In order to minimize the impact of secondary flows, the following should be carried out: v y v x , v z v x .
( v x v x x + v y v x y + v z v x z ) ~ ( v x v x x )
(2) Poiseuille Flow Hypothesis:
Assuming that the fluid flow in the working area along the X direction is a time-invariant one-dimensional linear flow [25], and conforms to the Hartmann flow in the two-dimensional straight channel, the Navier–Stokes equation and the magnetic diffusion equation need to be used. Since the coupling between the velocity vector and the magnetic field in the MHD can be expressed through Maxwell’s equations [26,27], the magnetic diffusion equation is composed of Maxwell’s equations, which are transformed by Faraday’s law, Ohm’s law, and Ampere’s law through differentiation, and its expression form is as follows:
B t = × ( v × B ) + η 2 B
In the working area, the conductive fluid moving in the X-axis direction will have a flow velocity distribution in the Z-axis direction, which follows the law of the Poisson lobe flow between the two plates [28], and its steady-state solution flow velocity is trapezoidal.
When z = ± k ( k = 1 / 2 h , h = h 2 h 1 ) , v x = 0 , and since V = v x v , the maximum velocity in X direction can be deduced as follows:
V x f = 1 2 μ f p f x ( z 2 k 2 ) = ( v x v ) ( 1 z 2 k 2 )
(3) Assumption from Ignoring the Induced Magnetic Field:
When the magnetic Reynolds number of the model is small, the magnetic field generated by the internal induced current (induced magnetic field) is much smaller than that of the external applied magnetic field. Therefore, the generation of the coupled magnetic field is not considered here, and the induced magnetic field is ignored.
(4) Ideal Fluid Hypothesis:
The conductive fluid is assumed to be an ideal fluid with a small viscosity and incompressibility. According to Bernoulli’s equation, we can conclude that the sum of the pressure of potential energy, upper potential energy, and kinetic energy is equal to the constant.
By partial differential treatment, the following results can be obtained:
p f x = 1 2 ρ 2 v x v x x = ρ v x v x x
Combined with the above assumptions, the Navier–Stokes equation can be simplified as follows:
v x t + v x l = 1 ρ ( ρ d h v x l ) μ v x v h 2 1 ρ σ B 2 v + a
The simplified Navier–Stokes equations are obtained by the Laplace transformation as:
a ( s ) = s 2 + ( μ h 2 + 1 d h l ) s s + μ h 2 + σ B 2 ρ v x ( s )
From the law of electromagnetic induction, the formula for the induced electric potential can be derived as:
U = 0 k B v x d h
Therefore, the transfer function of the runway-type fluid ring linear vibration sensor based on the MHD is:
| U ( s ) V o u t ( s ) | = B h ( s + μ h 2 + σ B 2 ρ ) s + ( μ h 2 + 1 d h ρ l )

3. Single-Channel Blind Signal Separation Based on SSA and FastICA

3.1. Principle of SSA

Based on subspace technology [29], SSA can identify and enhance the characteristics of periodic signals in a time series without prior information, and is not subject to the advantages of sine wave constraints. It combines the empirical orthogonal function to reconstruct the time series, and can extract the signal from the time series data with limited length, which is especially suitable for studying the periodic oscillation behavior of a time series. Currently, SSA is also widely used to decompose single-channel signals into trend, oscillation, and noise components.
A singular analysis is based on the practice of constructing the observed sequence path matrix, and the trajectory of matrix decomposition and reconstruction. The aim is to extract signals on behalf of the original time series signal of different components, such as long-term trend, and periodic and noise signals, and the structure of time series analysis can be used for further predictions. In this paper, SSA denoising actually disregards the main components of the original data and considers the unimportant components as noise, so as to remove this noise. The basic steps are to first arrange the original sequence into a matrix form with a delay, construct the trajectory matrix, and then use singular value decomposition. Then, a new matrix formed by the obtained principal components is grouped and diagonally averaged [19]. The main calculation steps of SSA technology are given as follows:
(1) Decomposition
The decomposition stage consists of two sub-steps, namely embedding and SVD. In the embedding step, the N sampled signal vectors y = [y(1), y(2), …, y(N)] are mapped into the M × L trajectory matrix Y:
Y = [ y ( 1 ) y ( 2 ) y ( L ) y ( 2 ) y ( 3 ) y ( L + 1 ) y ( M ) y ( M + 1 ) y ( L + M ) ]
The trajectory matrix Y is called the Hankel matrix, and the value of each element is equal to the value of the adjacent elements in the diagonal direction of the matrix. Where M is the decomposition length, its value can be flexibly selected according to the actual processing of the signal, usually not more than half of the signal length.
In the SVD step, a singular value decomposition (SVD) is performed on the trajectory matrix Y, that is, Y = UΣVT, where U is called the left matrix, only has values on the main diagonal, which is the singular value, and all other elements are zero. V represents the right matrix. In addition, both U and V are unit orthogonal matrices, which satisfy UUT = E, VVT = E. As it is difficult to directly decompose the trajectory matrix, the covariance matrix of the trajectory matrix, S = YYT, was first calculated and then the eigenvalue decomposition of S was carried out to obtain eigenvalues λ1 > λ2 > … > λL ⩾ 0 and the corresponding eigenvectors U1, U2, …, Un; then, U = [U1, U2, …, UL], λ11/2 > λ21/2 > … > λL1/2 = 0 represents the singular spectrum of the original sequence.
(2) Reconstruction
Based on the singular value decomposition (SVD), the original time series are divided into several disjoint groups, and the original data are divided into several series according to their degree of importance, so as to prepare for the next reconstruction. Finally, the first M (decomposition length) main components are selected to form a new time series to achieve reconstruction.

3.2. FastICA Algorithm

FastICA is a fast optimization iterative algorithm derived from the principle of non-Gaussianity maximization [30]. Considering that an algorithm to maximize the objective function is needed in practical applications, the FastICA algorithm is adopted in this paper, which adopts a fixed point and Newton iteration; therefore, it has a fast convergence speed.
FastICA assumes that the acquired signal X is preprocessed to obtain a new signal z, where w represents the row vectors in the separation matrix W. The FastICA algorithm maximizes the non-Gaussianity of the vector WTX (where T denotes the transpose) by a Newton iteration. From the theory of information theory, we know that negative entropy remains unchanged during invertible linear transformation. In order to obtain a measure for the Gaussian distribution that is zero and always non-negative, negative entropy is the best measure.
In general, the observation data are first centralized and whitened to obtain the observation signal with zero mean and unit variance, and then FastICA analysis is performed. The steps of the FastICA algorithm are as follows: (1) Firstly remove the mean of the observed signal X and then whiten it to obtain the variable z. The maximum iteration number N and tolerance error ε are set. (2) Set the initial weight vector w and make the iteration number i = 1. (3) Update W through a Newton iteration and normalize it after each iteration. (4) Suppose i = i + 1, if it does not converge, i.e., i ≤ M and w(i + 1) k − w(i)k > ε, then go back to step 3. The flow chart of the FastICA algorithm is shown in Figure 3.

3.3. The Single-Channel MHD Linear Vibration Signal Blind Source Separation Method Based on SSA and FastICA

In real experiments, the output signal of the shaking table is disturbed by a variety of noises when it is received and collected by the MHD linear vibration sensor. Since the noise sources are unknown, in this study, we assume that there are N noise sources, namely X(t) = S(t) + Ni(t), i = 1, 2, 3, …, n, where, X(t) is the acquisition signal, S(t) is the source signal, and N(t) is the noise signal.
The collected linear vibration signal is a mixed signal containing multiple noise signals. The signal-to-noise ratio of such mixed signals is low, and the effect of using ICA directly is not good; thus, the mixed signals need to be denoised first. In this paper, a single-channel line vibration signal separation method based on SSA-FastICA is proposed. Firstly, the output signal of the MHD linear vibration sensor is collected, and the collected signal is a noise-laden mixed signal. Then, the SSA method is preliminarily used to denoise the acquired signal with noise. Finally, FastICA is used to separate the collected signal after denoising, and the unmixed signal was obtained. The separation process of mixed signals based on SSA and FastICA is shown in Figure 4.
The specific implementation steps of single-channel MHD linear vibration signal blind source separation based on SSA-FastICA are as follows:
(1) The MHD linear vibration sensor prototype is placed on the vibration table for vibration detection and its output line vibration signal is collected.
(2) The collected MHD linear vibration signal is processed by SSA, the contribution rate of the eigenvalues of each order component is calculated, and according to the contribution rate of the eigenvalues, the useful signals and noise signals are determined from the components.
(3) The correlation coefficient between each component and the original collected signal is determined, and the threshold of the correlation coefficient is calculated. The relationship between the correlation coefficient and its threshold is used as the criterion to reconstruct the signal components. The reconstructed signal is combined with the original ultrasonic echo signal to form a new multidimensional observation signal, and its dimension is consistent with the estimated number of source signals.
(4) The FastICA algorithm is used to process the newly constructed observation signal to obtain the separated noise signal and useful linear vibration signal.

4. Experimental Analysis and Verification

4.1. Signal Collected of the MHD Linear Vibration Sensor

The MHD linear vibration sensor is a vibration sensor independently developed by the laboratory. In order to measure the accuracy of its actual output linear vibration signal, we chose a quartz flexible accelerometer as the comparison sensor, and its physical diagram is shown in Figure 5. The quartz flexible accelerometer selected in this paper has a high reliability and can be used for dynamic testing. It is a standard vibration sensor with an analog output similar to an MHD line vibration sensor. The parameters of the quartz flexible accelerometer are shown in Table 1.
The experimental wiring diagram is shown in Figure 6. An MHD linear vibration sensor and quartz flexible accelerometer were placed on the vibration table, and the acquisition card was connected. Finally, the acquired line vibration signal is displayed through the upper computer. Before the experiment, the sensor was fixed with a fixture to prevent position deviation caused by the vibration of the shaking table and affecting the accuracy of the sensor output signal.
Four groups of experimental data were collected, including 0.6 V, 10 Hz; 0.6 V, 60 Hz; 0.8 V, 5 Hz; and 0.8 V, 60 Hz, denoted as Collected Signal 1, Collected Signal 2, Collected Signal 3, and Collected Signal 4, respectively. The corresponding spectrum analysis of the four groups of collected signals is shown in Figure 7.

4.2. The Collected MHD Linear Vibration Signal Is Denoised

For the MHD linear vibration collected signals doped with a variety of noises, their SNR is low and not suitable for direct signal separation. Therefore, the collected linear vibration signals (Collection Signal 1–4) should be preliminarily denoised to obtain a higher signal SNR.
The Collected Signals 1–4 were processed by EEMD, VMD, wavelet threshold denoising, and SSA, respectively, and the useful linear vibration signal components and noise components under the four denoising methods were obtained, respectively.
In order to evaluate the denoising performance of EEMD, VMD, wavelet threshold denoising, and SSA, SNR was selected as the evaluation index of signal denoising. The formula of SNR after the preliminary noise reduction in the collected signal is:
S N R = 10 lg i = 1 N s i 2 ( t ) i = 1 N [ s i ( t ) X i ( t ) ] 2
where si(t) represents the source signal, where denoised linear vibration signal is substituted, and Xi(t) indicates the collected signals. When the SNR after noise reduction is higher, the separation effect is greater. The SNR analysis of these four groups of useful linear vibration signals was carried out, as shown in Table 2.
As seen from the above table, among the collected signals of group 1, 2, 3, and 4, the signal after SSA processing has a higher SNR than the signal after EEMD, VMD, and wavelet threshold denoising processing, and the SNR of the signal after SSA processing is 1.6 to 32 times that of the other three kinds. Therefore, it can be seen that SSA technology is a better choice for pre-denoising the collected signals when comparing the SNR improvement of the four groups of collected signals after processing.

4.3. Collection Signal Noise Reduction and Re-Separation

For mixed signals with a low SNR, after SSA is first used to reduce noise and improve the SNR, FastICA can be directly used to separate the collected signals with high SNR, and independent source signals, i.e., useful linear vibration signals, can be obtained. Blind signal separation was performed in four groups of collected signals to obtain six separated signals: separated signal 1, separated signal 2, separated signal 3, separated signal 4, separated signal 5, and separated signal 6. The similarity coefficient between each group of separated signals and the corresponding collected signals is shown in Table 3.
The four groups of signals with the highest similarity coefficients among the decomposed signals are found, as shown in Table 4. These four groups of signals are the line vibration signals obtained by processing the four groups of collected signals by the SSA-FastICA algorithm. It is observed that the output signals of the quartz flexible accelerometer are essentially the same as the decomposed linear vibration signals waveform, as shown in Figure 8. In contrast, in the quartz flexible acceleration sensor and our experimental test sensor, the output voltage signals amplitudes are different. Therefore, it is more appropriate to measure the output signal of the MHD linear vibration sensor and quartz flexible accelerometer by signal similarity. The similarity of the output signals is more than 98%. It can be seen that the single-channel linear vibration signal blind source separation method based on SSA-FastICA can effectively extract standard linear vibration signals. In order to test the repeatability of the sensor, multiple sets of data of the MHD linear vibration sensors at different times on different dates are collected. The formula for the standard error (SE) [31] is shown below:
S E = 1 n 1 i = 1 n ( X i x ¯ ) 2
x ¯ = 1 n i = 1 n X i
where Xi represents the collected signal segment and x ¯ represents the average of all signal segments.
The SE of the data segments of 10 Hz and 60 Hz is calculated, respectively. Eight data segments of 10 Hz and 60 Hz are randomly selected, as shown in Figure 9. The signal average SE of the quartz flexible accelerometer is 0.02348 mV at 10 Hz and 0.56456 mV at 60 Hz. The signal average SE of the MHD linear vibration sensor is 0.00620 mV at 10 Hz and 0.03684 mV at 60 Hz. After SSA-FastICA algorithm processing, the average SE of signals of the MHD linear vibration sensor is 0.00420 mV at 10 Hz and 0.00392 mV at 60 Hz. Therefore, it can be considered that the MHD linear vibration sensor has good repeatability.

5. Conclusions

Based on the theoretical basis of MHD, this paper designs and prepares a prototype of a runway-type fluid ring type linear vibration sensor, and samples and extracts the output linear vibration velocity signal. Since the acquired signals belong to mixed signals containing multiple complex noises, a single-channel linear vibration signal separation method combining SSA and FastICA is proposed through the study of blind signal extraction techniques, and the effectiveness of the method is verified by experimental results. The following conclusions are obtained:
  • The self-made MHD linear vibration sensor possesses the advantages of no mechanical wear between internal firmware and no additional power supply. However, its limitation is that it is not suitable for constant speed measurement.
  • For mixed signals containing multiple narrowband noise signals, there is a large difference between the signals directly separated by the ICA algorithm and the source signal. For the collection signal with more noise, noise reduction is needed to improve the SNR.
  • Compared with the processing methods of EEMD, VMD, and wavelet threshold noise reduction, the SSA algorithm can effectively suppress the narrowband periodic interference and white noise in the mixed signal, improve the SNR of the acquired signal, and effectively retain the useful information of the signal while reducing noise.
  • The proposed SSA-FastICA algorithm can effectively solve the problem of blind source separation and the extraction of the single-channel MHD linear vibration signals to minimize the dynamic measurement error. Thus, it shows that the SSA-FastICA method can be effective in single-channel MHD linear vibration signal separation and processing.

Author Contributions

M.X. and X.L. conceived and designed the study. M.X. and F.J. developed the sensor in the experiment. J.W. and J.M. carried out the sensor experiment, J.W. and L.Y. completed the data processing. M.X. and J.W. wrote the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research and the APC was funded by the National Natural Science Foundation of China (Grant No. 61903116) and the Natural Science Foundation of Anhui Province (Grant No. 1908085QE232).

Acknowledgments

The study was supported by the National Natural Science Foundation of China (No. 61903116) (applied January 2020), and the Natural Science Foundation of Anhui Province (No. 1908085QE232) (applied May 2019).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Addari, D.; Aglietti, G.S.; Remedia, M. Experimental and numerical investigation of coupled microvibration dynamics for satellite reaction wheels. J. Sound. Vib. 2017, 386, 225–241. [Google Scholar] [CrossRef]
  2. Xu, M.; Li, X.; Wu, T.; Chen, C.; Ji, Y. Error analysis of theoretical model of angular velocity sensor based on magnetohydrodynamics at low frequency. Sens. Actuators A Phys. 2015, 226, 116–125. [Google Scholar] [CrossRef]
  3. Ji, Y.; Li, X.; Wu, T.; Chen, C. Modified analytical model of magnetohydrodynamics angular rate sensor for low-frequency compensation. Sens. Rev. 2016, 36, 193–199. [Google Scholar] [CrossRef]
  4. Ji, Y.; Li, X.; Wu, T.; Chen, C. Theoretical and Experimental Study of Radial Velocity Generation for Extending Bandwidth of Magnetohydrodynamic Angular Rate Sensor at Low Frequency. Sensors 2015, 15, 31606–31619. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Rejer, I.; Gorski, P. MAICA: An ICA-based method for source separation in a low-channel EEG recording. J. Neural Eng. 2019, 16, 056025. [Google Scholar] [CrossRef] [PubMed]
  6. Albera, L.; Kachenoura, A.; Comon, P.; Karfoul, A.; Wendling, F.; Senhadji, L.; Merlet, I. ICA-based EEG denoising: A comparative analysis of fifteen methods. B. Pol. Acad. Sci-Tech. 2012, 60, 407–418. [Google Scholar] [CrossRef] [Green Version]
  7. Debals, O.; De Lathauwer, L. Stochastic and Deterministic Tensorization for Blind Signal Separation. In International Conference on Latent Variable Analysis and Signal Separation; Springer: Cham, Switzerland, 2015; pp. 3–13. [Google Scholar]
  8. Kopsinis, Y.; McLaughlin, S. Development of EMD-Based Denoising Methods Inspired by Wavelet Thresholding. IEEE Trans. Signal Process. 2009, 57, 1351–1362. [Google Scholar] [CrossRef]
  9. Jing-yi, L.; Hong, L.; Dong, Y.; Yan-sheng, Z. A New Wavelet Threshold Function and Denoising Application. Math. Probl. Eng. 2016, 2016, 3195492. [Google Scholar] [CrossRef] [Green Version]
  10. Maddirala, A.K.; Shaik, R.A. Separation of Sources from Single-Channel EEG Signals Using Independent Component Analysis. IEEE Trans. Instrum. Meas. 2018, 67, 382–393. [Google Scholar] [CrossRef]
  11. Zhou, M.; Zhao, H.; Xia, H.; Zhang, J.; Liu, Z.; Liu, C.; Gao, F. De-noising of photoacoustic sensing and imaging based on combined empirical mode decomposition and independent component analysis. J. Biophotonics 2019, 12, e201900042. [Google Scholar] [CrossRef]
  12. Jin, T.; Li, Q.; Mohamed, M.A. A Novel Adaptive EEMD Method for Switchgear Partial Discharge Signal Denoising. IEEE Access 2019, 7, 58139–58147. [Google Scholar] [CrossRef]
  13. Mika, D.; Budzik, G.; Jozwik, J. Single Channel Source Separation with ICA-Based Time-Frequency Decomposition. Sensors 2020, 20, 2019. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Liu, J.; Gao, W. Vibration Signal Analysis of Water Seal Blasting Based on Wavelet Threshold Denoising and HHT Transformation. Adv. Civ. Eng. 2020, 2020, 4381480. [Google Scholar] [CrossRef] [Green Version]
  15. Zhou, Y.; Fu, W.; Shi, F.; Zhang, X.; Zheng, X. Improved Single-Channel Blind Source Separation Algorithm based on Variational Mode Decomposition. Res. Sq. 2021. [Google Scholar] [CrossRef]
  16. Zou, L.; Chen, X.; Dang, G.; Guo, Y.; Wang, Z.J. Removing Muscle Artifacts from EEG Data via Underdetermined Joint Blind Source Separation: A Simulation Study. IEEE Trans. Circuits Syst. II Express Briefs 2020, 67, 187–191. [Google Scholar] [CrossRef]
  17. Guo, Y.; Huang, S.; Li, Y. Single-Mixture Source Separation Using Dimensionality Reduction of Ensemble Empirical Mode Decomposition and Independent Component Analysis. Circ. Syst. Signal. Pr. 2012, 31, 2047–2060. [Google Scholar] [CrossRef]
  18. He, P.; She, T.; Li, W.; Yuan, W. Single channel blind source separation on the instantaneous mixed signal of multiple dynamic sources. Mech. Syst. Signal. Pr. 2018, 113, 22–35. [Google Scholar] [CrossRef]
  19. Fu, J.; Cai, F.; Guo, Y.; Liu, H.; Niu, W. An Improved VMD-Based Denoising Method for Time Domain Load Signal Combining Wavelet with Singular Spectrum Analysis. Math. Probl. Eng. 2020, 2020, 1485937. [Google Scholar] [CrossRef]
  20. Li, B.; Zhang, L.; Zhang, Q.; Yang, S. An EEMD-Based Denoising Method for Seismic Signal of High Arch Dam Combining Wavelet with Singular Spectrum Analysis. Shock. Vib. 2019, 2019, 4937595. [Google Scholar] [CrossRef]
  21. Yang, Q.; Wang, K.; Ma, Y.; Wu, J.; Zhu, Q.; Gu, Y.; Qi, H. Design parameters analysis and verification of angular vibration sensor based on magnetohydrodynamics. Int. J. Metrol. Qual. Eng. 2021, 12, 2. [Google Scholar]
  22. Cheng, J.; Sacchi, M.D. Separation and reconstruction of simultaneous source data via iterative rank reduction. Geophysics 2015, 80, V57–V66. [Google Scholar] [CrossRef]
  23. Gómez, D.O.; Mininni, P.D.; Dmitruk, P. MHD simulations and astrophysical applications. Adv. Space Res. 2005, 35, 899–907. [Google Scholar] [CrossRef]
  24. Koochi, A.; Abadyan, M.; Gholami, S. Electromagnetic instability analysis of nano-sensor. Eur. Phys. J. Plus 2021, 136, 44. [Google Scholar] [CrossRef]
  25. Smolentsev, S. MHD duct flows under hydrodynamic “slip” condition. Theor. Comput. Fluid Dyn. 2009, 23, 557–570. [Google Scholar] [CrossRef] [Green Version]
  26. Koochi, A.; Abadian, F.; Rezaei, M.; Abadyan, M. Electromagnetic instability of electromechanical nano-bridge incorporating surface energy and size dependency. Phys. E Low-Dimens. Syst. Nanostructures 2021, 129, 114643. [Google Scholar] [CrossRef]
  27. Koochi, A.; Goharimanesh, M.; Gharib, M.R. Nonlocal electromagnetic instability of carbon nanotube-based nano-sensor. Math. Methods Appl. Sci. 2021. [Google Scholar] [CrossRef]
  28. Daprà, I.; Scarpi, G. Couette–Poiseuille flow of the Giesekus model between parallel plates. Rheol. Acta 2008, 48, 117–120. [Google Scholar] [CrossRef]
  29. Cheng, J.; Li, L.; Li, C.; Liu, Y.; Liu, A.; Qian, R.; Chen, X. Remove Diverse Artifacts Simultaneously from a Single-Channel EEG Based on SSA and ICA: A Semi-Simulated Study. IEEE Access 2019, 7, 60276–60289. [Google Scholar] [CrossRef]
  30. He, X.; He, F.; He, A. Super-Gaussian BSS Using Fast-ICA with Chebyshev–Pade Approximant. Circ. Syst. Signal. Pr. 2017, 37, 305–341. [Google Scholar] [CrossRef]
  31. Chai, T.; Draxler, R.R. Root mean square error (RMSE) or mean absolute error (MAE)?—Arguments against avoiding RMSE in the literature. Geosci. Model Dev. 2014, 7, 1247–1250. [Google Scholar] [CrossRef]
Figure 1. (a) The physical map of the quartz flexible accelerometer; (b) schematic diagram of the MHD linear vibration sensor principle.
Figure 1. (a) The physical map of the quartz flexible accelerometer; (b) schematic diagram of the MHD linear vibration sensor principle.
Sensors 22 09657 g001
Figure 2. Schematic diagram of basic dimensions of the MHD linear vibration sensor fluid channel.
Figure 2. Schematic diagram of basic dimensions of the MHD linear vibration sensor fluid channel.
Sensors 22 09657 g002
Figure 3. Flowchart of FastICA algorithm.
Figure 3. Flowchart of FastICA algorithm.
Sensors 22 09657 g003
Figure 4. Flow chart of mixed signal separation process based on SSA and FastICA.
Figure 4. Flow chart of mixed signal separation process based on SSA and FastICA.
Sensors 22 09657 g004
Figure 5. Physical graph of quartz flexible accelerometer.
Figure 5. Physical graph of quartz flexible accelerometer.
Sensors 22 09657 g005
Figure 6. Experimental connection diagram of signal collection of the MHD linear vibration sensor.
Figure 6. Experimental connection diagram of signal collection of the MHD linear vibration sensor.
Sensors 22 09657 g006
Figure 7. Frequency spectrum of Collected Signal 1–4.
Figure 7. Frequency spectrum of Collected Signal 1–4.
Sensors 22 09657 g007aSensors 22 09657 g007b
Figure 8. Comparison of output signals of quartz flexible accelerometer and linear vibration signals of MHD after SSA-FastICA processing.
Figure 8. Comparison of output signals of quartz flexible accelerometer and linear vibration signals of MHD after SSA-FastICA processing.
Sensors 22 09657 g008
Figure 9. (a) Random four signal segments at 10 Hz; (b) Random four signal segments at 60 Hz.
Figure 9. (a) Random four signal segments at 10 Hz; (b) Random four signal segments at 60 Hz.
Sensors 22 09657 g009
Table 1. Parameter list of quartz flexible accelerometer.
Table 1. Parameter list of quartz flexible accelerometer.
ParametersParameters Setting
0 g 4-h short-term stability≤20 μg
1 g 4-h short-term stability≤20 ppg
Partial value synthetic repeatability≤40 μg
Scale factor synthesis repeatability≤50 ppg
Comprehensive repeatability of nonlinear coefficient≤±20 μg/g2
Partial temperature coefficient≤±30 μg/°C
Scale factor temperature coefficient≤50 ppg/°C
Noise≤4 mV
Table 2. Table of signals SNR comparison after EEMD, VMD, wavelet threshold denoising, and SSA processing.
Table 2. Table of signals SNR comparison after EEMD, VMD, wavelet threshold denoising, and SSA processing.
Denoising AlgorithmCollected Signal 1Collected Signal 2Collected Signal 3Collected Signal 4
EEMD10.71210.71964.02920.8049
VMD34.990710.176917.560711.3256
Wavelet Threshold Denoising33.164511.318132.244611.1439
SSA58.642823.113752.441221.5621
Table 3. Similarity factor between the separated signals and the collection signals 1–4.
Table 3. Similarity factor between the separated signals and the collection signals 1–4.
Serial NumberCollected Signal 1Collected Signal 2Collected Signal 3Collected Signal 4
Separated Signal 12.61%13.67%1.68%9.02%
Separated Signal 298.45%98.82%2.59%99.26%
Separated Signal 33.20%4.60%16.31%4.53%
Separated Signal 411.98%1.54%−6.35%0.88%
Separated Signal 511.86%−3.50%6.97%3.89%
Separated Signal 61.54%2.20%98.15%−4.93%
Table 4. Similarity factor (max) between the separated signals and the Collected Signals 1–4.
Table 4. Similarity factor (max) between the separated signals and the Collected Signals 1–4.
Collected Signal 1Collected Signal 2Collected Signal 3Collected Signal 4
Similarity Factor (max)98.45%98.82%98.15%99.26%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xu, M.; Wang, J.; Mo, J.; Li, X.; Yang, L.; Ji, F. Single-Channel Blind Signal Separation of the MHD Linear Vibration Sensor Based on Singular Spectrum Analysis and Fast Independent Component Analysis. Sensors 2022, 22, 9657. https://doi.org/10.3390/s22249657

AMA Style

Xu M, Wang J, Mo J, Li X, Yang L, Ji F. Single-Channel Blind Signal Separation of the MHD Linear Vibration Sensor Based on Singular Spectrum Analysis and Fast Independent Component Analysis. Sensors. 2022; 22(24):9657. https://doi.org/10.3390/s22249657

Chicago/Turabian Style

Xu, Mengjie, Jianhan Wang, Jiahui Mo, Xingfei Li, Lei Yang, and Feng Ji. 2022. "Single-Channel Blind Signal Separation of the MHD Linear Vibration Sensor Based on Singular Spectrum Analysis and Fast Independent Component Analysis" Sensors 22, no. 24: 9657. https://doi.org/10.3390/s22249657

APA Style

Xu, M., Wang, J., Mo, J., Li, X., Yang, L., & Ji, F. (2022). Single-Channel Blind Signal Separation of the MHD Linear Vibration Sensor Based on Singular Spectrum Analysis and Fast Independent Component Analysis. Sensors, 22(24), 9657. https://doi.org/10.3390/s22249657

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop