Research Article  Open Access
An Extraction Method of Weak LowFrequency Magnetic Communication Signals Based on Multisensor
Abstract
It is a technical challenge to effectively remove the influence of magnetic noise from the vicinity of the receiving sensors on lowfrequency magnetic communication. The traditional denoising methods are difficult to extract highquality original signals under the condition of low SNR (the signaltonoise ratio). In this paper, we analyze the numerical characteristics of the lowfrequency magnetic field and propose the algorithms of the fast optimization of blind source separation (FOBSS) and the frequencydomain correlation extraction (FDCE). FOBSS is based on blind source separation (BSS). Signal extraction of low SNR can be implemented through FOBSS and FDCE. This signal extraction method is verified in multiple field experiments which can remove the magnetic noise by about 25 dB or more.
1. Introduction
Because the electromagnetic wave of RF (radio frequency) attenuates in the water or on land, the application of lowfrequency (LF, 30 Hz~300 kHz), ultralowfrequency (ULF, 300~3000 Hz), and extremely lowfrequency (ELF, 3~30 Hz) electromagnetic spectrum receives widespread attention [1, 2]. The electromagnetic wave of low frequency is mainly used in submarine communication [2] and resource exploration [3].
Many scholars have made great efforts toward the theory of propagation of lowfrequency electromagnetic field [4, 5]. They have created magnetic communication channel models [6] and analyzed the physical and statistical characteristics of channel noise [7]. Actually, many excellent signal processing methods of RF can be extended to that of low frequency. For the engineering application, numbers of signal denoising methods are proposed, such as the matched filtering [8], the optimization of receiving antenna [9], and denoising based on the noise distribution and empirical mode decomposition [10]. However, in the case of low signaltonoise ratio or existing complex noise components, these methods are difficult to work effectively.
In the application of lowfrequency magnetic field, not only the signal waveform should be obtained, but also the original signal should be acquired for locating magnetic source [11]. Magnetic communication receivers are often affected by the surrounding noise and signal distorts seriously. In this paper, the correlation of lowfrequency magnetic field is analyzed, and the algorithms of the fast optimization of blind source separation and the frequencydomain correlation extraction are proposed. The extraction of signal can be implemented through these algorithms. FOBSS can obtain the signal similar to the original signal waveform. FDCE recovers the original signal combining with FOBSS. In essence, we acquire more information (or more relationships) through multisensor comparing with single sensor. We associate the advantages of the algorithms and the inner relationships to achieve the signal extraction.
In this paper, contents are arranged as follows. Section 1 introduces the background of the lowfrequency field application and the present situation of lowfrequency magnetic signal processing methods; Section 2 introduces the basic theory of blind source separation; Section 3 analyzes the numerical characteristics of lowfrequency magnetic field sources and proposes FOBSS and FDCE; Section 4 presents the experimental design process and verifies algorithms through multiple experiments; finally, Section 5 summarizes the full text.
2. BSS
When signal is interfered by multiple noise sources under the condition of low SNR, the methods of signal decomposition (Fourier analysis, wavelet analysis, empirical mode decomposition, etc.) are not effective. BSS [12] can work well in some cases, but the requirements of BSS must be met.
In the signal processing of multisensor, assume that is the observation signal vector of sensors, and each observation element of is the linear mixture of signal sources . The linear model of is written as follows:where is the measured signal of the th sensor. is the unknown mixing coefficient between the th source and the th sensor. is the th signal source. Usually, and are equal. Equation (1) is expressed as matrix form:Generally, before signal separation, is preprocessed by spheroidizing decomposition. Blind source separation is to find the solution of inverse mixed matrix and estimate the optimal source signals. The process can be shown as The basic requirements should be satisfied: signal sources are statistically independent of each other; the rank of the covariance matrix is full. The results of separation are evaluated by the coherence. The timedomain coherence function between and is given bywhere is the covariance. The sequence and amplitude of the separated signal through BSS are uncertain. When the coherence is larger than 0.8, signal separation is successful. High coherence denotes that the waveform similarity of the separated result and the original signal is high.
3. Characteristics and Extraction Method of LowFrequency Magnetic Field Signal
3.1. The Correlation Characteristics of LowFrequency Magnetic Field
The propagations of lowfrequency and highfrequency electromagnetic field are in different ways (e.g., the conduction current mainly contributes to the lowfrequency field, and the displacement current plays a main role in highfrequency field). The propagation mode of lowfrequency field is extremely complex. However, the near field of lowfrequency magnetic field is approximately equivalent to the radiation field of magnetic source in electromagnetic field theory [13]. The vector potential and magnetic field excited by the changing current density source are, respectively, written aswhere is the permeability of air, is the speed of light, and is expressed as . is the frequency. is the wave number. is the observation point. is the source point. . The frequency band of magnetic source is finite (). Frequency is described as discrete form through the Fourier analysis. The field of observation point can be given bywhere is a factor related to the distance. Considering the wavelength of lowfrequency field that is much longer than the distance between the observation point and the source point (the ratio of wavelength to distance tends to be infinite), the phase of (7) can be expressed aswhere is the wavelength of . So the ratio of magnetic field at different observation points is And the coherence between the magnetic field of and is In theory, the magnetic fields of different observation points are highly relevant. The covariance matrix of () is (). The singular value decomposition of is represented as . If the rank of matrix is full, the locations of the observation points meet the requirements of signal separation. If the rank of matrix is not full, the locations of the observation points need to be adjusted (usually, placing the sensors randomly is appropriate; they need not to be relocated).
3.2. The Fast Optimization of Blind Source Separation
Because matrix (see Section 2) contains unknown elements, BSS is the process of rotating (see also Section 2). is a unit orthogonal matrix. Now there is a method that reduces the numbers of the unknown elements of , so the iteration time of finding can be shortened in calculation. According to the law of Givens rotation transformation [14], the mathematical process is given by where is the Givens matrix, is an arbitrary nonzero vector, and is an arbitrary standard unit vector. If is an arbitrary unit vector, there is a result: . Assume that ( denotes reverse rotation). The rotation process of ( is a standard orthogonal basis) can be expressed as ( is an arbitrary orthogonal basis). So the arbitrary orthogonal matrix of dimensional space can be acquired through the rotation of a standard unit matrix, and the number of rotations is at most times. Therefore can be expressed as by unknown elements. The iteration time of can be defined as ( is the th rotation angle). The iteration time of basic BSS is defined as ( is the element of ). Under the condition of the same computational resource, the iteration time of unknown elements is statistically shorter than that of unknown elements. Therefore the ratio of mean iteration time between and is
3.3. The FrequencyDomain Correlation Extraction
FOBSS can only acquire the waveforms similar to the original signals. The recovery of original signals requires excellent algorithms of signal amplitude extraction. In signal extraction, the parameters of the adaptive filtering [15] are difficult to be adjusted and the adaptive filter cannot work effectively for nonperiodic signal. The signal extraction can be implemented in the frequency domain based on the correlation function. The advantage of frequencydomain extraction contains the process of power spectrum estimation. The correlation reflects in each frequency point. and are, respectively, defined as the time series of signals. consists of the convolutions of various signals. is the known waveform similar to a component of . The basic model is simply described as where different components of are statistically independent and () is an unknown kernel function of linear convolution. Now in order to extract of , we get the crosscorrelation function between and . can be transformed to the corresponding frequency domain as by the Fourier transform. The mathematical process is described as where is the cross power spectrum density function between and , is the Fourier transform of , and is the power spectrum density function of . Therefore the result of extraction is where is the Fourier transform of . We give a deviation ratio of the extracted signal and the original signal such thatwhere is the original signal, is the extracted signal, and () is the mean value. If the mean values of and are zero, and are, respectively, replaced as the peaktopeak values and .
4. Processing of the Experimental Data
4.1. Experiment Design
We have carried out multiple sets of magnetic communication experiment. Experimental equipment consists of the transmitter, the receivers, and the magnetic noise sources. The transmitter is an audio power amplifier AE7224 and the output impedance of the amplifier is less than . The transmitter antenna is a coil loaded by a magnetic core that it looks like a cylinder which is about 1 meter long. Besides, the diameter of antenna is 0.07 meter and the coil has 300 windings. The transmitting current is 1.5 A. The receiving antenna is also a magnetic induced coil and its conversion factor is 100 mV/nT. The AD of the receiver is a 24bit converter and its LSB (Least Significant Bit) is 20 bits. The quantization noise is about 0.048 pT. The experiment sites are Kangxi Prairie in Beijing and An Guli Prairie in Hebei province, China. The scene is shown in Figure 1: the magnetic signal (the signal bandwidth is about 15 Hz) is transmitted with the carrier frequency of 480 Hz, and the transmitter is located at 100 meters in the Cartesian coordinate system; the receivers (multisensor) distribute in the range of 2 meters~4 meters away from the original point, and the sampling rate is 2400 Hz. The actual sampling rate should be higher than Nyquist sampling rate at a suitable value; the magnetic noise is from the car engines, and the engines are in the range of 2 meters~4 meters away from the receivers. Magnetic noise is stochastic, and its frequency bandwidth covers the signal’s bandwidth. So the signal contains a certain amount of inband noise. The received signals are intensely interfered by the magnetic noise in time domain and frequency domain. Only the signals of the multisensor close to the original point are applied to signal separation and signal extraction.
We have verified that the signal of the sensor at 50 meters is not affected by magnetic noise. When the noise sources and the transmitter are open in Figure 1, the coherence between the signal of sensor at 75 meters and that of sensor at 50 meters is 0.98. The coherence between the signal of sensor at 75 meters and that of sensor at original point nearby is 0. The coherence of the signal of sensor at 50 meters and that of sensor at original point nearby is 0. But when only the transmitter is open, the coherence of different sensors is 0.98. This phenomenon verifies that the reference sensor is not interfered by magnetic noise sources. The signal of the sensor (the reference sensor) at 50 meters is to prove the effectiveness of algorithms.
4.2. Data Processing
In Section 2, the signal model is theoretical. Under the condition of reality, the sensor itself contains micro noise, and the signal model is modified as another form:where is the mixing coefficient, is the th signal source, and is the micro noise of the th sensor ( is statistically independent). The procedure of signal extraction is shown in Figure 2, mainly including the preprocessing, the spectral subtraction [16], FOBSS, and FDCE. Firstly, whether the rank of the covariance matrix of the received signals is full or not should be confirmed; secondly, remove the micro noise of sensors through spectral subtraction; thirdly, signal separation is executed by FOBSS; finally, the corresponding signal waveforms are obtained by component matching, and the original signal is obtained by FDCE from the row data (the signals after bandpass filtering). The main parts are FOBSS and FDCE. Component matching is to select the required component after FOBSS (not matched filter). Besides, computing correlation function through the waveform of code is to decide which component is the required component. The required component has a similar waveform to the original signal.
Experiments are carried out by the combinations which consist in whether transmitter and magnetic noise sources are open or not, and the reference sensor at 50 meters is used to validate the correctness of algorithms. In order to illustrate data, we transform (4) to frequency domain (the frequencydomain coherence function): where is the cross power spectrum density function between and , is the power spectrum density function of , and is the power spectrum density function of .
As shown in Figure 1, magnetic noise sources are within the range of 2 meters~5 meters away from receivers. When only one magnetic noise is open, the coherence between the signals of different receivers is above 0.96 and there is no time delay between the signals of different receivers (these characteristics are consistent with the analysis of Section 3.1). Figure 3 shows the frequencydomain coherence, this feature meets the correlation characteristics of lowfrequency magnetic field. The reference sensor at 50 meters cannot be affected by magnetic noise, and we only use it to verify the correctness of algorithms (the reference sensor is not used for signal separation and signal extraction). The coherence between different noise sources is zero. There is no correlation between different physical sources, so components separation is possible.
When one noise source and the transmitter are open, the SNR is about −9 dB. When three noise sources and the transmitter are open, the SNR is about −20 dB. Figure 4 shows the timedomain waveforms of different receivers in the case that three magnetic noise sources and the transmitter are open, and the receiver at 50 meters is not influenced by noise sources. The sequence of sensors is arranged in Figure 1. The reference sensor is at 50 meters. Figure 5 shows the frequency spectra of Figure 4. The actual locations of sensors 1, 2, 3, and 4 have a little deviation compared with the locations shown in Figure 4. According to basic electromagnetic theory, the lowfrequency magnetic field value decreases inversely proportional to distance when the source can be deemed as one point source. Because three noise sources are distributed near the sensors, three sources are not equivalent to one point source when the distance and the range of source are approximately at the same level. So field values in Figure 4 are not strictly consistent with that of one point source. There maybe exists offset when peaks and troughs appear at the same time. Since the signals are intensely interfered in time domain and frequency domain, the signal processing methods should be based on components separation but not on the synthetic filtering (it maybe contains the matched filtering, the optimal notch filtering, etc.) based on signal decomposition. The signal separation performance is generally influenced by the noise of sensor itself, industrial interference, and the experiment errors.
After preprocessing and independent noise removing, the received signals will be processed by FOBSS. We make a comparison between the separated results and the reference signal. Figure 6 shows a comparison between the separated results and reference signal (the separation results and reference signal are adjusted to the same scale), and the coherence calculated by (4) is 0.89 (in the case that three magnetic noise sources and the transmitter are open). The coherence of the separated signal under multiple sets of experiments is higher than 0.8. The communication channel may be influenced slightly by industrial frequency interference. Since the data of Figures 6(a) and 6(b) are obtained at different time and places, it is reasonable that their envelope causes slight difference in a larger image display. Both the separated signal and the signal of the reference sensor can be decoded correctly in BPSK modulation.
(a) The comparison of separated results in time domain (one noise source)
(b) The comparison of separated results in time domain (three noise sources)
(c) The comparison of separated results in frequency domain (one noise source)
(d) The comparison of separated results in frequency domain (three noise sources)
The signal separation algorithm (FOBSS) can only acquire the waveforms similar to the original signals. The signal extraction is realized by (15) (FDCE). The processed results of multiple experiments are shown in Table 1; the timedomain coherence is given by (4), and the deviation ratio is given by (16).

The key purpose of this paper is to extract weak true signal masked by strong noise. The “weak” (SNR is too low) is the main problem. The extracted signal through FOBSS and FDCE can be decoded errorless and it has the same waveform as the reference signal. They also obtained the true amplitude. In the actual environment, the signal in the transmitter has minute difference from the signal in the air. In order to verify algorithms, the reference signal only can be the signal in the air but not that in the transmitter. Even though there is a slight influence factor, both the receiving end and the reference at 50 meters are affected (other algorithms can easily solve usually). The effectiveness cannot be influenced.
4.3. Data Analysis
The coherence larger than 0.8 indicates high correlation in statistical theory. The synthetic filtering based on the timefrequency domain is difficult to filter out inband noise; the ability of empirical mode decomposition to overcome inband noise is limited and it is difficult to identify the mode; when the situation is more complex, it needs signal separation through building conditions. We regard the similarity (the frequencydomain coherence) to the reference signal as the standard to evaluate the ability of different signal processing methods. There is a comparison between FOBSS and the synthetic filter in Figure 7. The coherence of FOBSS tends to 1, and the coherence of the synthetic filter is low. Simultaneously, because design of standard narrowband filter is not easy, the synthetic filter is difficult to solve the problem.
The signal preprocessing and independent noise removal are important work. Compared with basic BSS, FOBSS can effectively improve the computation time in multiple experiments. Even though (12) is a statistical result, FOBSS is more suitable for engineering application. In the original signal extraction, the deviation ratio of FDCE is about 6 percent. In short, according to the performance of FOBSS and FDCE, the proposed signal extraction method is effective.
5. Conclusion
In this paper, we have analyzed the characteristics of lowfrequency magnetic field and have proposed FOBSS and FDCE. FOBSS improves processing efficiency. Even though FDCE is essentially from Fourier transform, it can be simply implemented and works well. A systemic method of extraction has been formed through the combination of FOBSS and FDCE.
This method has been effectively verified in multiple field experiments. Its advantage is that it can extract the waveform and the true amplitude under the condition of low SNR. It can remove the magnetic noise by about 25 dB or more. In future work, we will study the statistic characteristics of magnetic noise and give the mathematical relationship between the performance of the algorithms and independent noise.
Weak signal extraction is a difficult problem in signal processing. For the engineering application of lowfrequency magnetic field, this paper has made preliminary progress in narrowband signal processing. In order to expand the application range of lowfrequency signal, such as oil exploration and mining exploration, we need to increase the signal bandwidth. Because the processing of wideband signal has more difficulties, we will make further study.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This research is sponsored by the National High Technology Research and Development Program of China (Grant no. 2014AA093407) and the National Natural Science Foundation of China (Grant no. 41374186).
References
 A. Abubakar, T. M. Habashy, V. L. Druskin, L. Knizhnerman, and D. Alumbaugh, “2.5D forward and inverse modeling for interpreting lowfrequency electromagnetic measurements,” Geophysics, vol. 73, no. 4, pp. F165–F177, 2008. View at: Publisher Site  Google Scholar
 S. A. Wolf, J. R. Davis, and M. Nisenoff, “Superconducting extremely low frequency (ELF) magnetic field sensors for submarine communications,” IEEE Transactions on Communications, vol. 22, no. 4, pp. 549–554, 1974. View at: Publisher Site  Google Scholar
 R. E. Grimm, “Lowfrequency electromagnetic exploration for groundwater on Mars,” Journal of Geophysical Research: Planets, vol. 107, no. E2, p. 1, 2002. View at: Google Scholar
 W. Pan and K. Li, Propagation of SLF/ELF Electromagnetic Waves, Springer, 2014.
 L. O. Løseth, H. M. Pedersen, B. Ursin, L. Amundsen, and S. Ellingsrud, “Lowfrequency electromagnetic fields in applied geophysics: waves or diffusion?” Geophysics, vol. 71, no. 4, pp. W29–W40, 2006. View at: Publisher Site  Google Scholar
 H. E. Rowe, “Extremely low frequency (ELF) communication to submarines,” IEEE Transactions on Communications, vol. 22, no. 4, pp. 371–385, 1974. View at: Publisher Site  Google Scholar
 E. C. Field Jr. and M. Lewinstein, “Amplitudeprobability distribution model for VLF/ELF atmospheric noise,” IEEE Transactions on Communications, vol. 26, no. 1, pp. 83–87, 1978. View at: Publisher Site  Google Scholar
 J. E. Evans and A. S. Griffiths, “Design of a sanguine noise processor based upon worldwide extremely low frequency (ELF) recordings,” IEEE Transactions on Communications, vol. 22, no. 4, pp. 528–539, 1974. View at: Publisher Site  Google Scholar
 P. R. Bannister, “Orbiting transmitter and antenna for spaceborne communications at ELF/VLF to submerged submarines,” in AGARD, ELF/VLF/LF Radio Propagation and Systems Aspects 14 p (SEE N9330727 1132), vol. 1, 1993. View at: Google Scholar
 Y. Jeng, M.J. Lin, C.S. Chen, and Y.H. Wang, “Noise reduction and data recovery for a VLFEM survey using a nonlinear decomposition method,” Geophysics, vol. 72, no. 5, pp. F223–F235, 2007. View at: Publisher Site  Google Scholar
 A. P. Nickolaenko and I. G. Kudintseva, “A modified technique to locate the sources of ELF transient events,” Journal of Atmospheric and Terrestrial Physics, vol. 56, no. 11, pp. 1493–1498, 1994. View at: Publisher Site  Google Scholar
 M. T. Akhtar, T.P. Jung, S. Makeig, and G. Cauwenberghs, “Recursive independent component analysis for online blind source separation,” in Proceedings of the IEEE International Symposium on Circuits and Systems (ISCAS '12), pp. 2813–2816, IEEE, May 2012. View at: Publisher Site  Google Scholar
 C. H. Papas, Theory of Electromagnetic Wave Propagation, Courier Dover Publications, 2013.
 F. Ling, “Givens rotation based least squares lattice and related algorithms,” IEEE Transactions on Signal Processing, vol. 39, no. 7, pp. 1541–1551, 1991. View at: Publisher Site  Google Scholar
 S. T. Smith, Geometric optimization methods for adaptive filtering [PhD thesis], Harvard University, 1993.
 T. Inoue, H. Saruwatari, Y. Takahashi et al., “Theoretical analysis of musical noise in generalized spectral subtraction based on higher order statistics,” IEEE Transactions on Audio, Speech, and Language Processing, vol. 19, no. 6, pp. 1770–1779, 2011. View at: Google Scholar
Copyright
Copyright © 2015 Chao Huang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.