# An Extraction Method of Weak Low-Frequency Magnetic Communication Signals Based on Multisensor

Authors: Chao Huang; Xin Xu; Dunge Liu; Wanhua Zhu; Xiaojuan Zhang; Guangyou Fang
Journal: Journal of Electrical and Computer Engineering
 (2015)

---

## Abstract

It is a technical challenge to effectively remove the influence of magnetic noise from the vicinity of the receiving sensors on low-frequency magnetic communication. The traditional denoising methods are difficult to extract high-quality original signals under the condition of low SNR (the signal-to-noise ratio). In this paper, we analyze the numerical characteristics of the low-frequency magnetic field and propose the algorithms of the fast optimization of blind source separation (FOBSS) and the frequency-domain 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.


---

## Body

## 1. Introduction

Because the electromagnetic wave of RF (radio frequency) attenuates in the water or on land, the application of low-frequency (LF, 30 Hz~300 kHz), ultralow-frequency (ULF, 300~3000 Hz), and extremely low-frequency (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 low-frequency 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 signal-to-noise ratio or existing complex noise components, these methods are difficult to work effectively.In the application of low-frequency 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 low-frequency magnetic field is analyzed, and the algorithms of the fast optimization of blind source separation and the frequency-domain 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. Section1 introduces the background of the low-frequency field application and the present situation of low-frequency magnetic signal processing methods; Section 2 introduces the basic theory of blind source separation; Section 3 analyzes the numerical characteristics of low-frequency 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 thatx
=
x
1
,
x
2
,
…
,
x
n
T is the observation signal vector of n sensors, and each observation element of x is the linear mixture of m signal sources s
=
s
1
,
s
2
,
…
,
s
m
T. The linear model of x
i is written as follows:(1)
x
i
=
∑
j
=
1
m
a
i
j
s
j
,
i
=
1,2
,
…
,
n
,where x
i is the measured signal of the ith sensor. a
i
j is the unknown mixing coefficient between the jth source and the ith sensor. s
j is the jth signal source. Usually, m and n are equal. Equation (1) is expressed as matrix form:(2)
x
=
A
·
s
,
A
∈
a
i
j
n
×
m
.Generally, before signal separation, x is preprocessed by spheroidizing decomposition. Blind source separation is to find the solution of inverse mixed matrix W and estimate the optimal source signals. The process can be shown as (3)
y
t
=
W
·
x
t
=
W
·
A
·
s
t
=
s
^
t
.The basic requirements should be satisfied: (
1
) signal sources are statistically independent of each other; (
2
) the rank of the covariance matrix x is full. The results of separation are evaluated by the coherence. The time-domain coherence function between y
i and s
j is given by(4)
ρ
y
i
s
j
=
cov
⁡
y
i
t
,
s
j
t
cov
⁡
y
i
t
cov
⁡
s
j
t
,where cov
(
·
,
·
) 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 Low-Frequency Magnetic Field Signal

### 3.1. The Correlation Characteristics of Low-Frequency Magnetic Field

The propagations of low-frequency and high-frequency electromagnetic field are in different ways (e.g., the conduction current mainly contributes to the low-frequency field, and the displacement current plays a main role in high-frequency field). The propagation mode of low-frequency field is extremely complex. However, the near field of low-frequency 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 J
(
r
′
,
t
) are, respectively, written as(5)
A
r
,
t
=
μ
0
4
π
∫
V
′
J
r
′
,
t
-
r
/
c
r
d
V
′
=
A
r
e
-
j
ω
t
,
(6)
B
r
=
∇
×
A
r
,
t
=
∇
×
A
r
e
-
j
ω
t
,where μ
0 is the permeability of air, c is the speed of light, and A
(
r
) is expressed as A
(
r
)
=
(
μ
0
/
4
π
)
∫
V
′
(
J
r
′
e
i
k
r
/
r
)
d
V
′. ω is the frequency. k is the wave number. r is the observation point. r
′ is the source point. r
=
|
r
-
r
′
|. The frequency band of magnetic source is finite (ω
l
≤
ω
≤
ω
h). Frequency is described as discrete form through the Fourier analysis. The field of observation point r
i can be given by(7)
B
r
i
=
∇
×
A
′
r
i
∑
ω
=
ω
l
ω
h
e
-
j
ω
t
-
k
r
i
,where A
′
(
r
i
) is a factor related to the distance. Considering the wavelength of low-frequency 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 as(8)
ω
t
-
k
r
i
=
ω
t
-
2
π
r
i
λ
ω
⟶
ω
t
,where λ
ω is the wavelength of ω. So the ratio of magnetic field at different observation points is (9)
B
r
i
B
r
j
=
∇
×
A
′
r
i
∇
×
A
′
r
j
.And the coherence between the magnetic field of B
(
r
i
) and B
(
r
j
) is (10)
ρ
B
r
i
B
r
j
=
1
.In theory, the magnetic fields of different observation points are highly relevant. The covariance matrix of B
(
r
i
)  (i
=
1,2
,
…
,
m) is E
(
B
B
T
) (B
=
B
r
1
,
B
r
2
,
…
,
B
r
m
T). The singular value decomposition of E
(
B
B
T
) is represented as E
(
B
B
T
)
=
Q
Σ
Q
T. 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 matrixW (see Section 2) contains m
2 unknown elements, BSS is the process of W rotating x (see also Section 2). W is a unit orthogonal matrix. Now there is a method that reduces the numbers of the unknown elements of W, so the iteration time of finding W can be shortened in calculation. According to the law of Givens rotation transformation [14], the mathematical process is given by (11)
T
ξ
k
=
T
1
m
T
1
m
-
1
⋯
T
12
ξ
k
=
ξ
k
e
j
,where T
1
i is the Givens matrix, ξ
k is an arbitrary nonzero vector, and e
j is an arbitrary standard unit vector. If ξ
k is an arbitrary unit vector, there is a result: ξ
k
=
T
-
1
e
j. Assume that T
-
1
=
T
′ (T
-
1 denotes reverse rotation). The rotation process of e (e is a standard orthogonal basis) can be expressed as T
′
e
=
T
12
′
T
13
′
⋯
T
1
m
′
=
ξ (ξ is an arbitrary orthogonal basis). So the arbitrary orthogonal matrix of m-dimensional space can be acquired through the rotation of a standard unit matrix, and the number of rotations is at most m
-
1 times. Therefore W can be expressed as W
(
θ
1
,
θ
2
,
…
,
θ
m
-
1
) by m
-
1 unknown elements. The iteration time of W
(
θ
1
,
θ
2
,
…
,
θ
m
-
1
) can be defined as t
(
θ
1
,
θ
2
,
…
,
θ
m
-
1
) (θ
i is the ith rotation angle). The iteration time of basic BSS is defined as t
(
w
11
,
w
12
,
…
,
w
i
j
…
,
w
m
m
) (w
i
j is the element of W). Under the condition of the same computational resource, the iteration time of m
-
1 unknown elements is statistically shorter than that of m
2 unknown elements. Therefore the ratio of mean iteration time between t
(
θ
1
,
θ
2
,
…
,
θ
m
-
1
) and t
(
w
11
,
w
12
,
…
,
w
i
j
…
,
w
m
m
) is (12)
E
t
θ
1
,
θ
2
,
…
,
θ
m
-
1
E
t
w
11
,
w
12
,
…
,
w
m
m
≈
m
-
1
m
2
.


### 3.3. The Frequency-Domain 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 frequency-domain extraction contains the process of power spectrum estimation. The correlation reflects in each frequency point. u
(
t
) and v
(
t
) are, respectively, defined as the time series of signals. u
(
t
) consists of the convolutions of various signals. v
(
t
) is the known waveform similar to a component of u
(
t
). The basic model is simply described as (13)
u
t
=
∑
i
=
1
n
h
i
t
∗
u
i
t
,
v
t
=
u
k
t
,
k
=
1,2
,
…
,
n
,where different components of u
(
t
) are statistically independent and h
i
(
t
) (i
=
1,2
,
…
,
n) is an unknown kernel function of linear convolution. Now in order to extract h
k
(
t
)
∗
u
k
(
t
) of u
(
t
), we get the cross-correlation function R
u
v
(
τ
) between u
(
t
) and v
(
t
). R
u
v
(
τ
) can be transformed to the corresponding frequency domain as S
u
v
(
ω
) by the Fourier transform. The mathematical process is described as (14)
R
u
v
τ
=
h
k
τ
∗
R
u
k
u
k
τ
,
S
u
v
ω
=
h
k
ω
S
u
k
u
k
ω
,
h
k
ω
=
S
u
v
ω
S
u
k
u
k
ω
=
S
u
v
ω
S
v
v
ω
,where S
u
v
(
ω
) is the cross power spectrum density function between u
(
t
) and v
(
t
), h
k
(
ω
) is the Fourier transform of h
k
(
t
), and S
u
k
u
k
(
ω
) is the power spectrum density function of u
k
(
t
). Therefore the result of extraction is (15)
h
k
t
∗
u
k
t
=
F
-
1
h
k
ω
V
ω
,where V
(
ω
) is the Fourier transform of v
(
t
). We give a deviation ratio λ
e of the extracted signal and the original signal such that(16)
λ
e
=
E
z
^
t
-
z
t
E
z
t
×
100
%
,where z
(
t
) is the original signal, z
^
(
t
) is the extracted signal, and E(·) is the mean value. If the mean values of z
(
t
) and z
^
(
t
) are zero, z
(
t
) and z
^
(
t
) are, respectively, replaced as the peak-to-peak values z
p
-
p
(
t
) and z
^
p
-
p
(
t
).


## 3.1. The Correlation Characteristics of Low-Frequency Magnetic Field

The propagations of low-frequency and high-frequency electromagnetic field are in different ways (e.g., the conduction current mainly contributes to the low-frequency field, and the displacement current plays a main role in high-frequency field). The propagation mode of low-frequency field is extremely complex. However, the near field of low-frequency 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 J
(
r
′
,
t
) are, respectively, written as(5)
A
r
,
t
=
μ
0
4
π
∫
V
′
J
r
′
,
t
-
r
/
c
r
d
V
′
=
A
r
e
-
j
ω
t
,
(6)
B
r
=
∇
×
A
r
,
t
=
∇
×
A
r
e
-
j
ω
t
,where μ
0 is the permeability of air, c is the speed of light, and A
(
r
) is expressed as A
(
r
)
=
(
μ
0
/
4
π
)
∫
V
′
(
J
r
′
e
i
k
r
/
r
)
d
V
′. ω is the frequency. k is the wave number. r is the observation point. r
′ is the source point. r
=
|
r
-
r
′
|. The frequency band of magnetic source is finite (ω
l
≤
ω
≤
ω
h). Frequency is described as discrete form through the Fourier analysis. The field of observation point r
i can be given by(7)
B
r
i
=
∇
×
A
′
r
i
∑
ω
=
ω
l
ω
h
e
-
j
ω
t
-
k
r
i
,where A
′
(
r
i
) is a factor related to the distance. Considering the wavelength of low-frequency 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 as(8)
ω
t
-
k
r
i
=
ω
t
-
2
π
r
i
λ
ω
⟶
ω
t
,where λ
ω is the wavelength of ω. So the ratio of magnetic field at different observation points is (9)
B
r
i
B
r
j
=
∇
×
A
′
r
i
∇
×
A
′
r
j
.And the coherence between the magnetic field of B
(
r
i
) and B
(
r
j
) is (10)
ρ
B
r
i
B
r
j
=
1
.In theory, the magnetic fields of different observation points are highly relevant. The covariance matrix of B
(
r
i
)  (i
=
1,2
,
…
,
m) is E
(
B
B
T
) (B
=
B
r
1
,
B
r
2
,
…
,
B
r
m
T). The singular value decomposition of E
(
B
B
T
) is represented as E
(
B
B
T
)
=
Q
Σ
Q
T. 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 matrixW (see Section 2) contains m
2 unknown elements, BSS is the process of W rotating x (see also Section 2). W is a unit orthogonal matrix. Now there is a method that reduces the numbers of the unknown elements of W, so the iteration time of finding W can be shortened in calculation. According to the law of Givens rotation transformation [14], the mathematical process is given by (11)
T
ξ
k
=
T
1
m
T
1
m
-
1
⋯
T
12
ξ
k
=
ξ
k
e
j
,where T
1
i is the Givens matrix, ξ
k is an arbitrary nonzero vector, and e
j is an arbitrary standard unit vector. If ξ
k is an arbitrary unit vector, there is a result: ξ
k
=
T
-
1
e
j. Assume that T
-
1
=
T
′ (T
-
1 denotes reverse rotation). The rotation process of e (e is a standard orthogonal basis) can be expressed as T
′
e
=
T
12
′
T
13
′
⋯
T
1
m
′
=
ξ (ξ is an arbitrary orthogonal basis). So the arbitrary orthogonal matrix of m-dimensional space can be acquired through the rotation of a standard unit matrix, and the number of rotations is at most m
-
1 times. Therefore W can be expressed as W
(
θ
1
,
θ
2
,
…
,
θ
m
-
1
) by m
-
1 unknown elements. The iteration time of W
(
θ
1
,
θ
2
,
…
,
θ
m
-
1
) can be defined as t
(
θ
1
,
θ
2
,
…
,
θ
m
-
1
) (θ
i is the ith rotation angle). The iteration time of basic BSS is defined as t
(
w
11
,
w
12
,
…
,
w
i
j
…
,
w
m
m
) (w
i
j is the element of W). Under the condition of the same computational resource, the iteration time of m
-
1 unknown elements is statistically shorter than that of m
2 unknown elements. Therefore the ratio of mean iteration time between t
(
θ
1
,
θ
2
,
…
,
θ
m
-
1
) and t
(
w
11
,
w
12
,
…
,
w
i
j
…
,
w
m
m
) is (12)
E
t
θ
1
,
θ
2
,
…
,
θ
m
-
1
E
t
w
11
,
w
12
,
…
,
w
m
m
≈
m
-
1
m
2
.


## 3.3. The Frequency-Domain 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 frequency-domain extraction contains the process of power spectrum estimation. The correlation reflects in each frequency point. u
(
t
) and v
(
t
) are, respectively, defined as the time series of signals. u
(
t
) consists of the convolutions of various signals. v
(
t
) is the known waveform similar to a component of u
(
t
). The basic model is simply described as (13)
u
t
=
∑
i
=
1
n
h
i
t
∗
u
i
t
,
v
t
=
u
k
t
,
k
=
1,2
,
…
,
n
,where different components of u
(
t
) are statistically independent and h
i
(
t
) (i
=
1,2
,
…
,
n) is an unknown kernel function of linear convolution. Now in order to extract h
k
(
t
)
∗
u
k
(
t
) of u
(
t
), we get the cross-correlation function R
u
v
(
τ
) between u
(
t
) and v
(
t
). R
u
v
(
τ
) can be transformed to the corresponding frequency domain as S
u
v
(
ω
) by the Fourier transform. The mathematical process is described as (14)
R
u
v
τ
=
h
k
τ
∗
R
u
k
u
k
τ
,
S
u
v
ω
=
h
k
ω
S
u
k
u
k
ω
,
h
k
ω
=
S
u
v
ω
S
u
k
u
k
ω
=
S
u
v
ω
S
v
v
ω
,where S
u
v
(
ω
) is the cross power spectrum density function between u
(
t
) and v
(
t
), h
k
(
ω
) is the Fourier transform of h
k
(
t
), and S
u
k
u
k
(
ω
) is the power spectrum density function of u
k
(
t
). Therefore the result of extraction is (15)
h
k
t
∗
u
k
t
=
F
-
1
h
k
ω
V
ω
,where V
(
ω
) is the Fourier transform of v
(
t
). We give a deviation ratio λ
e of the extracted signal and the original signal such that(16)
λ
e
=
E
z
^
t
-
z
t
E
z
t
×
100
%
,where z
(
t
) is the original signal, z
^
(
t
) is the extracted signal, and E(·) is the mean value. If the mean values of z
(
t
) and z
^
(
t
) are zero, z
(
t
) and z
^
(
t
) are, respectively, replaced as the peak-to-peak values z
p
-
p
(
t
) and z
^
p
-
p
(
t
).


## 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 than1
×
1
0
6
Ω. 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 24-bit 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 f
s 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 in-band 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.Figure 1
Field experiment.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 Figure1, 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 Section2, 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:(17)
x
i
=
∑
j
=
1
n
a
i
j
s
j
+
N
i
,
i
=
1,2
,
…
,
n
,where a
i
j is the mixing coefficient, s
j is the jth signal source, and N
i is the micro noise of the ith sensor (N
i 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 band-pass 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.Figure 2
The procedure of signal extraction.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 frequency-domain coherence function): (18)
ρ
y
i
s
j
ω
=
S
y
i
s
j
ω
S
y
i
y
i
ω
S
s
i
s
i
ω
,where S
y
i
s
j
(
ω
) is the cross power spectrum density function between y
i and s
j, S
y
i
y
i
(
ω
) is the power spectrum density function of y
i, and S
s
i
s
i
(
ω
) is the power spectrum density function of s
j.As shown in Figure1, 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 frequency-domain coherence, this feature meets the correlation characteristics of low-frequency 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.Figure 3
The frequency-domain coherence between magnetic noise signals of different receivers.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. Figure4 shows the time-domain 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 low-frequency 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.Figure 4
Time-domain waveforms of the received signals (three noise sources).Figure 5
Frequency-domain waveforms of the received signals (three noise sources).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. Figure6 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.Figure 6
The separated results.
(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 time-domain coherence is given by (4), and the deviation ratio is given by (16).Table 1
The results of signal processing.
Experiment site
Number of noise sources
Time-domain coherenceρ
Deviation ratioλ
e
An Guli Prairie
1
0.9924
3.42%
An Guli Prairie
1
0.9533
4.18%
An Guli Prairie
1
0.9501
4.22%
Kangxi Prairie
3
0.8945
7.51%
Kangxi Prairie
3
0.8184
9.03%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 time-frequency domain is difficult to filter out in-band noise; the ability of empirical mode decomposition to overcome in-band 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 frequency-domain 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 Figure7. 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.Figure 7
The comparison of frequency-domain coherence.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.


## 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 than1
×
1
0
6
Ω. 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 24-bit 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 f
s 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 in-band 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.Figure 1
Field experiment.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 Figure1, 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 Section2, 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:(17)
x
i
=
∑
j
=
1
n
a
i
j
s
j
+
N
i
,
i
=
1,2
,
…
,
n
,where a
i
j is the mixing coefficient, s
j is the jth signal source, and N
i is the micro noise of the ith sensor (N
i 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 band-pass 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.Figure 2
The procedure of signal extraction.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 frequency-domain coherence function): (18)
ρ
y
i
s
j
ω
=
S
y
i
s
j
ω
S
y
i
y
i
ω
S
s
i
s
i
ω
,where S
y
i
s
j
(
ω
) is the cross power spectrum density function between y
i and s
j, S
y
i
y
i
(
ω
) is the power spectrum density function of y
i, and S
s
i
s
i
(
ω
) is the power spectrum density function of s
j.As shown in Figure1, 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 frequency-domain coherence, this feature meets the correlation characteristics of low-frequency 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.Figure 3
The frequency-domain coherence between magnetic noise signals of different receivers.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. Figure4 shows the time-domain 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 low-frequency 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.Figure 4
Time-domain waveforms of the received signals (three noise sources).Figure 5
Frequency-domain waveforms of the received signals (three noise sources).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. Figure6 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.Figure 6
The separated results.
(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 time-domain coherence is given by (4), and the deviation ratio is given by (16).Table 1
The results of signal processing.
Experiment site
Number of noise sources
Time-domain coherenceρ
Deviation ratioλ
e
An Guli Prairie
1
0.9924
3.42%
An Guli Prairie
1
0.9533
4.18%
An Guli Prairie
1
0.9501
4.22%
Kangxi Prairie
3
0.8945
7.51%
Kangxi Prairie
3
0.8184
9.03%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 time-frequency domain is difficult to filter out in-band noise; the ability of empirical mode decomposition to overcome in-band 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 frequency-domain 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 Figure7. 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.Figure 7
The comparison of frequency-domain coherence.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 low-frequency 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 low-frequency magnetic field, this paper has made preliminary progress in narrowband signal processing. In order to expand the application range of low-frequency 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.



---

*Source: 312347-2015-05-10.xml*