Research on Stochastic Resonance Detection Method for Periodic Signals under Low SNR and α-Stable Noise
In Dual Sequence Frequency Hopping (DSFH) communication mode, aiming at improving the detection performance to weak signal under low signal-to-noise ratio (SNR) conditions, stochastic resonance (SR) detection method is proposed. First, the α-stable distribution is used as the impulsive noise model and the influence of α value on the properties of α-stable noise is analyzed. Second, the transmitting and receiving signal model of DSFH communication system is introduced. The SR method is used to detect DSFH signal. In order to analyze the output signal, the fractional Fokker–Planck equation (FFPE) is established, and a new simplified solution method based on sampling decision time is proposed to solve the time-varying fractional differential equation. Base on the theoretical solution of FFPE, a binary hypothesis test statistic is constructed to quantify the signal detection probability and false alarm probability, and the detection performance is analyzed. Finally, simulation experiments verify the theoretical conclusions. The minimum effective SNR for SR detection is obtained, and it is about −20 dB, which provides a theoretical basis for the application of SR in the DSFH communication system.
High frequency (HF)/very high frequency (VHF) hopping communication is the primary means in the campaign and tactical communication. In the HF and VHF band, there are many kinds of radiation sources, including communication, radar, navigation, identification of friend or foe (IFF), and other equipments. The electromagnetic environment is very bad. Under the condition of communication countermeasure, the communication signal will be seriously suppressed by the enemy’s malicious interference. People have devoted themselves to research an emergency communication means suitable for low signal-to-noise ratio (SNR). Dual Sequence Frequency Hopping (DSFH) is a similar kind of communication method. The DSFH communication signal receiving can be attributed to the signal detection problem under very low SNR. Stochastic resonance (SR) method have been used to detect weak signals, but the noise background is often Gaussian white noise, which does not match the battlefield impulsive noise environment, resulting in large deviation of the research results. For DSFH communication mode, there is no research on the application of SR method for communication signal detection under impulse noise. The motivation of the research is to study an effective communication method for low SNR and impulsive noise.
SR detection of DSFH communication signal under impulse noise is the main content of this paper. DSFH is a kind of frequency hopping communication mode, based on “the medium is the message” . DSFH was first proposed in [2, 3], the basic transmitting and receiving models were constructed, and its antijamming performance was analyzed. In the DSFH model, both the transmitter and the receiver have two channels which represent two states of 0 and 1. Different from the traditional communication mode, the existence state of radio frequency (RF) sequences represents the symbols 0 and 1. DSFH signal does not need to be demodulated, which reduces the quality requirements of the received signal. However, the antijamming ability is naturally stronger. Aiming at the poor antijamming ability of dual channels, linear frequency modulated (LFM) wave and pseudorandom code were modulated in DSFH RF signal to enhance the system's antijamming ability in [4–6].
Stochastic resonance (SR) theory is a powerful tool for weak signal detection under low SNR conditions. SR is applied to detect a weak harmonic signal or weak fault feature in [7–9]. The optimal local detector for the SR output based on the offset coefficient is studied in . The detection of BPAM, BFSK, and BPSK signal using SR is realized in [11, 12]. In the DSFH communication mode, the feasibility of SR is demonstrated and the quantitative results to a suitable range of SNR and BER are given in [13, 14].
The SR principle is that the weak signal uses noise energy drive particles to resonance between two potential wells. Its essence is a kind of energy transfer from noise to weak signal, which increases the SNR. Therefore, noise is an essential condition for SR system, and selecting a suitable noise model is the key to study DSFH communication performance. Gaussian white noise is used in classical SR, and many researchers study SR under Gaussian white noise [15–18]. Although the Gaussian distribution has a solid theoretical foundation and mature application, many investigations have shown that it is not suitable to describe impulsive noise. This mismatch will lead to performance degradation or failure of the signal detection algorithm.
Battlefield wireless communication environment contains many random intense pulses, which show typical non-Gaussian property. The common non-Gaussian mathematical models mainly include Gaussian mixture distribution, t-distribution, generalized Gaussian distribution, and α-stable distribution. Alpha-stable distribution is a typical thick-tailed distribution, which is successfully used in modeling atmospheric and marine noise, and it is also widely used in other impulsive noise situations. It has a complete mathematical definition and characteristic function. It is convenient for theoretical analysis and closely related to the statistical physical model [19–23]. Therefore, in this paper, α-stable distribution is used as the impulsive noise model, so SR detection is studied under α-stable noise.
There are two main ways to study SR. One is the signal simulation method based on the Langevin equation (LE), the other is the theoretical method by solving probability density function (PDF) based on Fokker–Planck equations (FPE), and the SR output signal is analyzed based on the solution.
In [24, 25], the solution process of fractional Fokker–Planck Equations (FFPE) is studied and the stationary solution of the SR system is obtained by the approximate method. The path integral solution (PIS) is adopted for determining the PDF and mean first-passage time (MFPT) in [26, 27]. The bistable system SR phenomenon with the periodic input signal and α-stable noise is demonstrated by the numerical simulation method in . The SR phenomenon under α-stable noise is verified, and the bistable system’s response excited by constant signal and periodic square wave signal are studied in [29–31]. The parameters induced SR are studied and the relationship between system parameters , , and SNR with α-stable noise is analyzed by simulation experiments in . In recent years, many domestic and foreign scholars have utilized different potential functions to realize the SR detection under α-stable noise [32–37]. The problem of weak signal detection using a cascade SR system is studied in . The study of SR in α-stable noise has attracted more and more attention [39–44].
In this paper, the SR method is applied to detect DSFH signal under impulsive noise, and a new simplified solution method based on sampling decision time is proposed to solve the time-varying fractional differential equation. The accuracy of theoretical solution is demonstrated by simulation experiments, and some quantitative performances are obtained for the first time. By analyzing parameter correlation between bistable system, periodic signal and noise, the influence of parameter changes on SNR, suitable SNR range, and achievable BER can be obtained quantitatively. The detection probability and false alarm probability can be described qualitatively.
2. Alpha-Stable Noise
2.1. Definition and Mathematical Expression
α-stable distribution is a mathematical empirical model suitable for impulsive noise modeling. It satisfies the generalized central limit theorem. The sum of infinite independent and identically distributed random variables is α-stable distribution. Unlike the central limit theorem, the random variables may not be finite variance, which makes it more widely used. Alpha-stable distribution does not have a unified PDF, but it has a characteristic function as follows:where and Different types of impulsive noise can be described by adjusting its parameters.
is the characteristic parameter, which determines the impulsivity of α-stable distribution. The smaller the value of is, the stronger the impulsivity is, the higher the probability of the abnormal value is, and the thicker the tail is in PDF. Since α-stable distribution turns into Gaussian distribution when , α-stable distribution is a generalized form of Gaussian distribution.
is the skew parameter which determines the skewness of α-stable distribution. When , it is called symmetric α stable distribution.
is the scale parameter, which is used to describe the random variables dispersion from the mean or median value.
is the shift parameter, which is used to describe the position of the PDF of α-stable distribution. When and , it is called the standard α stable distribution. When and , the α-stable distribution becomes Gaussian distribution. When and , the α-stable distribution becomes Cauchy distribution. When and , the α-stable distribution becomes Levy distribution.
2.2. Generation of α-Stable Distribution Random Variables
Simulation is a main mean to study α-stable noise and signal processing in the noise. The generation of α-stable distribution random variables is an important condition. If is α-stable distribution random numbers, which can be generated by the Janicki Weron (JW) algorithm :where , , is a random variable uniform distributed over , and is a random variable exponentially distributed with a unit mean. The PDF curve of α-stable distribution with a different α is shown in Figure 1. It can be seen that the smaller α is, the stronger the impulsivity is.
3. Theoretical Analysis of SR Signal Detection in DSFH Communication System
3.1. Analysis of DSFH Receiving Signal Property
3.1.1. The Transmitting Signal Model of DSFH
In DSFH model, when sending data “0,” channel 0 is used as the data channel to transmit frequency signal, and channel 1 as the dual channel does not transmit frequency signal, as shown in Figure 2.
The carrier frequency of channel 0 is which is controlled by pseudorandom (PN) sequence , and the carrier frequency of channel 1 is which is controlled by the PN sequence . When the transmitting symbol is 0, the carrier frequency controlled by the PN sequence is transmitted. When the transmitting symbol is 1, the carrier frequency controlled by the PN sequence is transmitted. At the time of , the sinusoidal signal with the frequency of is transmitted when the transmitted symbol is 0. Otherwise, the sinusoidal signal with the frequency of is transmitted when the transmitted symbol is 1. Moreover, the signals and are both sinusoidal. After the channel switch, the final transmitted signal of DSFH is the combination of and .
Assume that the transmitted data is , for the main analysis, and the baseband filter is neglected. Then, the transmitted signal is the sine carrier with a frequency of , and the RF signal of DSFH can be represented aswhere is the hop duration, is the step signal, is the RF of the nth hop, and and .
3.1.2. The Receiving Signal Model of DSFH
The super-heterodyne receiver is adopted in DSFH mode, which is depicted in Figure 3. The receiver receives the signal of two channels according to the same frequency hopping pattern of the transmitter. After receiving the signal , the RF front-end carries out mixing and band-pass filtering with the carrier signal controlled by the corresponding PN sequence. The waveform and frequency of the intermediate frequency (IF) signal obtained by the two branches are the same, which can be shown aswhere is the preset frequency of the IF signal.
By analyzing the DSFH model, we can see that the received signal after frequency mixing is a single frequency sinusoidal signal. So, the problem of DSFH signal receiving becomes the problem of judging whether the sinusoidal signal exists.
3.2. Theoretical Analysis of SR Signal Detection under α-Stable Noise
Driven by the α-stable noise and periodic signal, the output of the bistable SR system can be described by LE, which is shown aswhere is a weak external sinusoidal signal and is asymmetric unbiased α stable distribution white noise with characteristic parameter and noise intensity . In this paper, zero-skewness and zero-shift will be considered, so , , and according to formula (6), . A typical nonlinear potential system is utilized in the bistable system as
LE describes amplitude property of the output signal. The FFPE corresponding to the statistical property of the output signal is shown aswhere is the PDF of the output signal, which represents the amplitude statistical property at time . is assumed to be expressed as , and formula (7) contains a time-varying term without stationary solution. Since the signal transformation is instantaneous, it can be assumed that the signal reaches a stable state instantaneously and no longer changes with time as . The sampling value on signal decision time in (7) can be expressed aswhere is regarded as a constant , and (8) can be rewritten aswhere , , , and are constants.
According to the definition of Gyunwald–Letnikov (GL) fractional calculus, after the differential processing, the equation iswhere , , , and and are left and right bounds of the signal amplitude boundary, respectively. The and space domain is divided into divisions in the differential processing.
Let and , with
Let , and is a matrix with the size of , wherewithwhere , , and .
The Gauss elimination method is used to obtain the solution of equation (13), which is also the probability density of SR output signal. The PDF without weak sinusoidal signal input can be also obtained by the same method. The difference between the two states can be analyzed theoretically.
Due to the bistable system’s effect, the PDF has changed from a single peak structure to an asymmetric bimodal structure, indicating that particles gather in two potential wells. The SNR was calculated and SNR peak appears obviously. It proves the occurrence of SR phenomenon. With sinusoidal signal driving, the peak value is higher than that without driving, which indicates that the weak periodic force affects particles’ aggregation. The difference is amplified by the SR system, which proves the effectiveness of SR to detect weak sinusoidal signal.
3.3. Performance Analysis of Detection and Reception
In the DSFH system, two channels receive signals separately. Since the two channels cannot receive two signals at the same time, signal reception can be considered as binary hypothesis test problem.where is SR output driven by only impulsive noise, is the SR output driven by weak sinusoidal signal and impulsive noise of DSFH system received, and is the sampling decision time. Based on Bayesian criterion, maximum a posterior probability (MAP) can be used for binary communication problems, and the test statistic can be obtained aswhere is the decision threshold and is expressed as
The detection probability is the probability that the real state is and it is judged as . The false alarm probability is the probability that the state is , but it is judged as . The PDF with or without periodic signal, which has been solved by theoretical method, is shown in Figure 4. The decision area and decision probability of SR is marked, where is the decision threshold.
Since the motion orbit of the SR particles is symmetrical, the PDF of the SR system output signal and the decision domain are also symmetrical. The decision domain of is . The decision domain of is .
can be obtained as
can be obtained as
Based on the PDF of the SR system output signal, we can obtain and .
4. Simulation Results and Discussion
To the typical bistable system represented by equation (6), according to the adiabatic approximation theory, there are two identical potential wells and one potential barrier. The potential well position is located at , the barrier center is at 0, and the barrier height is . When there is only noise and no sinusoidal signal excitation, the particles will vibrate in the potential wells with small noise, while the particles will transition between the two potential wells according to Kramers escape rate when the noise increases, and the transition is random with no periodicity. However, SR will occur in the system if the sinusoidal signal is added, and the transition frequency of the particles is the same as the sinusoidal signal.
4.1. Verification of SR Phenomenon
For the bistable system in equation (6), based on the requirements of a low-frequency small signal for the adiabatic approximation theory, sine wave is selected as the periodic input signal. The frequency of the sinusoidal signal is 10 Hz, and the amplitude is 0.1. The signal is in the α-stable noise with , , and SNR = −15 dB. The sampling frequency is 200 KHz. To increase the adaptability and dynamic of the bistable system, let , and b ranges from to . The simulation results is shown in Figure 5(a), where the vertical axis is the output SNR, and the horizontal axis is the parameter , and the actual value is . It can be seen that when is higher than , the SNR reaches the maximum region.
Let and , and the simulation results are shown in Figure 5(b), where the vertical axis is SNR, and the horizontal axis is the noise intensity . The amplitude of the input period weak signal is 0.1. The output SNR increases with the increase of noise intensity obviously, and the peak value appears when is about 0.3. It indicates that the energy of the noise transfers to the signal during the process, and the SR effect is better near . The adiabatic approximation theory gives an analytical expression for the relationship between noise intensity and SNR. There is no analytical expression for impulsive noise, but the phenomenon is the same as SR under Gaussian noise. The results are the same as the previous theoretical calculation results. Figure 5(c) shows the probability density curve of the output signal. Red indicates the output without sinusoidal signal, and black represents the output after adding sinusoidal signal. It can be seen that the signal probability density is more concentrated after adding the signal, and the difference is more obvious.
The SR phenomenon can also be proved by the time-domain output signal of the bistable system. As shown in Figure 6, the first line shows the input signal, the second line shows the output signal of the SR system, and the third line shows the periodic input signal. As shown in Figure 6(a), in the case of no sinusoidal signal input, although the output signal of the bistable system is also two-state switching, the periodicity cannot be seen. As shown in Figure 6(b), after adding a sin signal with amplitude of 0.1, the periodic property of the input signal cannot be seen because the signal is covered by noise. However, the periodicity of the output is more obvious after the SR system, which also proves the SR phenomenon.
Figure 7 shows the PDF of the bistable system output signal under different noise intensities. Red is the SR output without a sinusoidal signal, and black is the output after adding a sinusoidal signal. Figure 7(a) shows the PDF of the bistable system output signal when the noise intensity is 0.1. Due to the low noise intensity, the system does not have a potential well transition, and the output signal is mainly in a single well. When the sinusoidal signal is added, the SR phenomenon occurs. Figure 7(b) shows the PDF of the bistable system output signal when the noise intensity is 0.15. There is no sinusoidal signal input in the system, but the distribution of the two potential wells is not uniform. Increasing the noise intensity further, as shown in Figure 7(c), the transition between two potential wells can occur with or without sinusoidal signal, but SR only occurs when sinusoidal signal input.
4.2. Influence of Noise Parameter on SR
Alpha is an important parameter of α-stable distribution, which determines the impulsivity of noise. The smaller is, the stronger the impulse and the thicker the tail of PDF will be, which indicates that there are more abnormal values. Different impulse degrees will have a significant impact on the signal reception, especially the SR. The following experiment mainly focuses on the influence of on SR.
First, the stronger the impulse of noise, the more likely it is to cause SR. As shown in Figures 8(a)–8(d), there are the input and output SNR curves of the SR system when the noise parameter is 2.0, 1.5, 0.85, and 0.5, respectively. The red line is the input SNR curve calculated according to the 10 Hz sinusoidal signal and noise amplitude in time domain. The blue line is the mixed signal and noise SNR curve calculated by FFT in frequency domain. It can be seen that the SNR curves of the two methods are the same. The black line is the SNR curve of SR output signal. Compared with the FFT method, SR can generate SNR gain, which is about 22 dB, 19 dB, 17 dB, and 14 dB. It indicates that SR occurs, and the stronger the impulsivity is, the smaller the gain is. However, in the SR stage, the stronger the noise impulsivity is, the more likely it is to induce SR. When the noise intensity is 0.5, 0.25, 0.02, and near-zero, the SNR output peak is reached. It shows that the larger the abnormal value of the noise, the easier the energy transfer and the more likely to generate SR.
Under α-stable noise, the theoretical and simulation values of PDF after SR processing to the same SNR and a different are shown in Figure 9. With the same SNR, it can be seen that the smaller the is, the more different the PDF of SR output signal is. The reason is that can be used to characterize the impulsivity of α-stable noise. The smaller the is, the stronger the color property of noise is. The large the difference between α-stable distribution and Gaussian white noise is, the more favorable it is for SR to start resonance and achieve the best resonance state and the better the resonance effect is. When , α-stable noise becomes Gaussian white noise, and the PDF difference is the smallest. It is proved that the driving impact of Gaussian white noise on SR is worse than that of α-stable noise, that is, color noise is more conducive to SR.
4.3. Influence of Input SNR to SR
Under α-stable noise, the theoretical and simulation value of PDF after SR processing to the same and different SNR are shown in Figure 10. It can be seen that when is the same, the larger the SNR is, the greater the sinusoidal signal amplitude and the driving force are, which means the better the starting resonance and SR effect. Under α-stable noise, the difference of output PDF with or without sinusoidal signal makes it possible for signal detection and processing.
4.4. Analysis of System Receiving Performance
and can be obtained based on PDF. The receiver operation characteristic (ROC) curve of sinusoidal signal after SR processing in α-stable noise is shown in Figure 11, reflecting the relationship between and . A different α corresponds to different ROC curves when SNR = −8 dB. Regardless of the value of 0.5, 1.5, or 2.0, the ROC curve is far above the line, which represents good detection performance.
The BER can be obtained based on and . The BER curve of DSFH sinusoidal signal with α-stable noise after SR processing is shown in Figure 12. It can be seen that the theoretical value and simulation value of decrease with the increase of SNR. It is because the improvement of SNR can effectively improve the detection probability. Meanwhile, the smaller the is, the smaller the is. The reason is that the noise power is constant when the SNR is fixed, and the smaller represents the more substantial noise impact, which is more conducive to electromagnetic particles in SR potential wells. Therefore, the detection performance is better. It can be proved that the DSFH signal can be detected and processed by SR under α-stable noise, and the application range of SNR is further extended to about −20 dB.
According to the characteristics of the DSFH signal under α-stable noise, the SR method is introduced to detect the sinusoidal signal. Based on the instantaneity of electrical signal transformation, the sampling decision time is used to solve FFPE and the theoretical solution of SR output signal PDF under α stable noise is obtained and verified by simulation. The stages of starting resonance and stable resonance are proposed, and the characteristic is analyzed according to the stage in different noises. Through SNR, noise intensity, and PDF difference analysis, the influence of noise on staring resonance and SNR gain is quantitatively obtained. Noise impulsivity is beneficial to start SR, but it will reduce SNR gain. Based on PDF and MAP criterion, the detection probability and false alarm probability are calculated, and the ROC and Pe-SNR curves are obtained. The availability of the detection method is verified, and the lowest applicable SNR of the system is about −20 dB. The quantitative conclusion provides the basis for SR to detect and receive DSFH signal under α-stable noise and provides a reference for SR application in electromagnetic signal processing. In the future, we will continue to study how to further reduce the BER by designing different receiving structures.
The data used to support the findings of this study are available from the corresponding or submitting author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
The research was funded by the project “Research on anti jamming technology of dual sequence frequency hopping communication”, Natural Science Foundation of Hebei Province of China (Grant nos. F2019506031, F2019506037, and F2020506036), Frontier Innovation Program of Army Engineering University under Grant no. KYSZJQZL2005, Basic Frontier Science and Technology Innovation Program of Army Engineering University under Grant no. KYSZJQZL2020, and Research Innovation Development Funding of Army Engineering University Shijiazhuang Campus under Grant no. KYSZJQZL1902.
F. H. P. Fitzek, “The medium is the message,” in Proceedings of the IEEE International Conference on Communications, pp. 5016–5021, Istanbul, Turkey, August 2006.View at: Google Scholar
H. Zhao, H. D. Quan, and P. Z. Cui, “Follower-jamming resistible multi-sequence frequency hopping wireless communication system,” Systems Engineering and Electronics, vol. 37, no. 3, pp. 671–678, 2015.View at: Google Scholar
Z. Q. Tang, H. D. Quan, and H. X. Sun, “Anti-jamming method of communication based on binary- sequence frequency hopping,” Journal of Detection & Control, vol. 42, no. 1, pp. 56–61, 2020.View at: Google Scholar
R. Gao, “Research on signal detection technology in complex electromagnetic environment,” Xidian University, Xi’an, China, 2015, Doctoral Thesis.View at: Google Scholar
L. Liang, “Research on digital signal detection technology based on nonlinear stochastic resonance,” Xidian University, Xi’an, China, 2018, Doctoral Thesis.View at: Google Scholar
G. K. Liu, H. D. Quan, H. X. Sun, P. Z. Cui, C. Chi, and S. L. Yao, “Stochastic resonance detection method for the dual-sequence frequency hopping signal under extremely low signal-to-noise radio,” Journal of Electronics and Information Technology, vol. 41, no. 10, pp. 2342–2349, 2019.View at: Google Scholar
G. Hu, Random Forces and Nonlinear Systems, Shanghai Education Press, Shanghai, China, 1994.
X. Q. Hu, Stochastic Resonance Weak Characteristic Signal Detection Theory and Method, National Defense Industry Press, Beijing, China, 2012.
F. B. Duan, “Application of parameter-adjusted stochastic resonance in digital signal transmission,” Zhejiang University, Zhejiang, China, 2002, Doctoral Thesis.View at: Google Scholar
N. H. Nedev, “Analysis of the impact of impulse noise in digital subscriber line systems univers,” University of Edinburgh, Edinburgh, Scotland, 2003, Ph.D. Thesis.View at: Google Scholar
B. Dybiec and E. Gudowska-Nowak, “Stochastic resonance: the role of alpha-stable noise,” Acta Physica Polonica Seriries B, vol. 37, no. 5, pp. 1479–1490, 2006.View at: Google Scholar
W. Y. Zhang, “Research on detection of weak periodic signal in levy noise based on stochastic resonance,” Shanghai Jiaotong University, Shanghai, China, 2009, Master Thesis.View at: Google Scholar
J. M. Huang, “Research on stochastic resonance phenomenon of nonlinear system under the background of stable distributed noise,” Zhejiang University, Zhejiang, China, 2012, Doctoral Thesis.View at: Google Scholar
G. Zhang, Y. Song, and T. Q. Zhang, “Characteristics analysis of exponential type monostable stochastic resonance under Levy noise,” Journal of Electronics & Information Technology, vol. 39, no. 4, pp. 893–900, 2017.View at: Google Scholar
S. B Jiao, J. Li, Q. Zhang, and G. Xie, “Stochastic resonance in time-delayed asymmetric monostable system with stable noise,” Journal of System Simulation, vol. 28, no. 1, pp. 139–146, 2016.View at: Google Scholar
S. B Jiao, D. Sun, J. Li, G. Xie, Y. L. Wu, and Q. Zhang, “Vibrational resonance in a periodic potential system with stable noise,” Acta Physica Sinica, vol. 66, no. 10, p. 100501, 2017.View at: Google Scholar
L. F. He, X. C. Zhou, G. Zhang, and T. Q. Zhang, “Stochastic resonance characteristics analysis of the new potential function under Levy noise and bearing fault detection,” Journal of Vibration and Shock, vol. 38, no. 12, pp. 53–62, 2019.View at: Google Scholar
Y. J. ZhangDuan, F. Z. Wang, and L. Liu, “Extraction of weak signal in cascaded system under Levy noise,” Computer Measurement & Control, vol. 27, no. 1, pp. 190–194, 2019.View at: Google Scholar
G. L. Zhang, X. L. Lv, and Y. M. Kang, “Parameter-induced stochastic resonance in overdamped system with stable noise,” Acta Physica Sinica, vol. 61, no. 4, pp. 1–8, 2012.View at: Google Scholar