Mobile Information Systems
Volume 2020 (2020), Article ID 8845160, 8 pages
https://doi.org/10.1155/2020/8845160

Research Article

The Mean First Passage Time of the Stochastic Resonance Driven by Dual-Sequence-Frequency-Hopping Signal and Noise

Zhaorui Li,1,2 Xiaobei Wu,¹ Peizhang Cui,² and Guangkai Liu³

¹School of Automation, Nanjing University of Science and Technology, Nanjing 210094, China
²Department of Electronics and Optical Engineering, Army Engineering University, Shijiazhuang 050003, China
³Beijing Institute of Tracking and Telecommunications Technology, Beijing 100094, China

Correspondence should be addressed to Peizhang Cui; dreamer_gkc@163.com

Received 21 March 2020; Revised 17 August 2020; Accepted 29 August 2020; Published 15 September 2020

Academic Editor: Adrian Kliks

Copyright © 2020 Zhaorui Li 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.

Abstract.

The mean first passage time (MFPT) represents the dynamic characteristic of stochastic resonance (SR). The study focuses on how can the Dual-Sequence-Frequency-Hopping (DSFH) signal influence the MFPT and any difficulty in solving the MFPT problem considering the DSFH signal. In this current study, the SR system driven by DSFH signal and Gaussian white noise is described with the parameters of the signal amplitude, the frequency of Intermediate Frequency (IF) of the receptive DSFH signal, the SR system parameter, scale transformation coefficient, the noise intensity, and the sampling multiple, firstly. Secondly, under the assumption that MFPT is small aqueous about the domain of 0, the nonautonomous differential equation with MFPT is transformed to a nonhomogeneous differential equation with one unknown variable coefficient of second order. Finally, the numerical solution of MFPT can be obtained by the method of Runge–Kutta. Theoretical and simulation results are shown as below: (1) the effect of the signal amplitude, the IF frequency, the noise intensity, the SR system parameter, and the scale transformation coefficient, for decrease the MFPT, are positive; however, the effect of the sampling multiple is negative; (2) the MFPT cannot follow the dynamic period of the SR controlled by the IF frequency, when SNR is low; (3) when SNR = −12 dB, the sampling multiple is 200, the IF frequency is 2100, and the duty cycle reaches 25% (available for DSFH signal detection with peek or valley decision Liu et al. (2019)), so we need to decrease the IF frequency or increase the SNR for further availability.

1. Introduction

The Dual-Sequence-Frequency-Hopping (DSFH) communication mode draws lessons from “the medium is the message” [1, 2]. Its typical character is that the message is modulated in the radio frequency, which is the hopping frequency controlled by the PN sequence. When the symbol 1 is transmitted, the hopping frequency controlled by PN 1 is transmitted. And when 0 is transmitted, the hopping frequency controlled by PN 0 is transmitted. The radio frequency of DSFH is combined by hopping frequency controlled by PN 1 and PN 0. The chosen medium is the communication channel; meanwhile, the unchosen one is the dual channel controlled by two PN sequences representing the symbol 0 or 1. Then, the receptive symbol is decided by whether the channel is occupied or not [3–5]. Due to the character of DSFH, the receptive signal of DSFH is a sinusoidal signal that is the typical input signal of stochastic resonance (SR). And the DSFH processed by stochastic resonance (SR) can work under strong color noise as a military emergency communication mode [6, 7]. SR is a nonlinear physical phenomenon that goes against the conventional concept that noise is always harmful for signal detection. Noise can enhance the detection performance of nonlinear SR systems when the signal, noise, and the SR parameters are matched. The SR theory is first proposed in the research of the earth glacier by Benzi et al. [8] and then demonstrated in physics, biology, and electronics [9–12]. And the application of the SR to signal procession is reviewed and discussed [13]. However, SR has a dynamic nature as a nonlinear system. The signal processing needs a steady state of the SR, not a dynamic one. So, when it is applied to signal detection, the researchers must ensure that the output of SR is already steady at the signal decision time during the bit period of DSFH. In other words, the SR system particles are jumping in some well instead of between wells at the decision time.

The mean first passage time (MFPT) is used to describe the dynamic characteristic of SR. The MFPT is defined as the time that a walker or a particle spends before crossing a given position for the first time [14]. As for the MFPT, most researchers mainly focus on the relationship between noise and the SR system. Guardia and San Miguel firstly study the escape time of bistable SR caused by multiple and additive Gaussian noise and obtain the analytic expression of MFPT related with potential function and noise intensity [15], which initiates the research of MFPT. Kang researches the MFPT in underdamped bistable Duffing oscillator [16] and the phenomenological gene transcriptional regulatory model [17] by the method of moments and the period steady solution of Fokker–Planck equation (FPE). Wang studies the MFPT in the Levine model, focusing on the white and color noise in the biological population system [18, 19]. Xu and Jin studies the MFPT in the tristable SR driven by correlated multiplicative and additive white noises [20, 21]. However, most research studies focus on the effect of MFPT caused by the SR system and noise. The adiabatic elimination theory is mostly used to obtain the MFPT with the assumption that the signal amplitude and frequency are small. And there has been no result which can be used to analyze the influence caused by the input signal frequency so far.

When the DSFH signal is considered in the MFPT problem, the nonautonomous item caused by the time-vary signal is introduced in the differential equation related to MFPT, which causes many difficulties in solving the MFPT. Aimed at solving this problem, we assume that MFPT is small aqueous about the domain of 0 based on the experience of communication signals detection and simulation results. And we try to research the relationship between MFPT and the DSFH signal, noise, and SR system. The paper is structured as follows. In Section 2, the SR system driven by DSFH signal and Gaussian white noise is described. In Section 3, the MFPT affected by DSFH signal and Gaussian white noise are analyzed, and the nonhomogeneous differential equation with MFPT of second order is solved with the method of Runge–Kutta. In Section 4, the relationship between MFPT and the DSFH signal, noise, and SR system is discussed in the theoretical and simulation view. Conclusions are given in Section 5.

2. The System Model of the SR Driven by DSFH Signal and Gaussian White Noise

2.1. The Transmitted Signal of the DSFH

The communication and dual carrier controlled by PN sequences are chosen by the transmitted symbol in the DSFH mode and described as Figure 1.

Figure 1: The transmitted structure of the DSFH.

Channel 0 and 1 are, respectively, the carrier and controlled by PN sequences and . When the send symbol is 0, the carrier controlled by PN sequences is transmitted. When the send symbol is 1, the carrier controlled by PN sequences is transmitted. At the time of t, the sine carrier with the frequency of is transmitted if the transmitted symbol is 0. Otherwise, the sine carrier represented of symbol 1 is transmitted. Furthermore, the signals and are sinusoidal. The final transmitted signal of DSFH is the combination of the and the , after the channel switch.

Assume that the transmitted data is , for the main analysis, the baseband filter is neglect. Then, the transmitted signal is the sine carrier with a frequency of . So, the radio frequency (RF) signal of DSFH can be represented as

where is the hop duration, is the step signal, is the radio frequency of the nth hop, and .

2.2. The Receptive Signal of the DSFH

The superheterodyne receiver is adopted in the DSFH mode, as depicted in Figure 2.

Figure 2: The received structure of the DSFH.

The signal is receipted under noise at the RF front end. The noise mostly is white Gaussian noise, the channel fading is neglect, and the receptive signal is ideal sinusoidal signal. Then, it will be mixed with the superheterodyne carriers controlled by two PN sequences. Subsequently, the analogy beat signal as the intermediate frequency (IF) signal be converted to digital one by the A/D convertor. The frequency and waveform of the IF signal obtained by the two receptive branches are the same but not at the same time, which can be descripted as

where is the preset frequency of the IF signal.

Noting that, in the DSFH system, the signal is just exit in one branch and the other branch is zero at the same time.

Due to the superheterodyne reception and the special modulation of DSFH, the received signal of the DSFH is the simple sinuous signal, which is the typical input signal of SR. With the reason that the and the sampling rate are set as 1 kHz and 200 kHz, the IF signal can be viewed as low-frequency signal, which can be processed by low-pass filter (LF). In the DSFH communication system, part information of the received signal such as the frequency is known in the reception end. So, we can design the LF based on this. Then, the scale transact unit (ST) completes the parameters of IF signal to fit the demand of SR.

When the receptive signal is deciding, the output of SR must be already steady. So, we need to analyze the dynamic characteristic of SR driven by the receptive signal of DSFH and Gaussian white noise.

2.3. The SR Processing of the Receptive Signal

The IF signal of the two receptive branches is the sine wave as depicted in equation (2). The overdamped bistable SR enforced by sine wave and noise can be depicted as follows:

where a and b are the SR parameters, is the periodic signal, is the white Gaussian noise satisfied with , , and D is the intensity of the noise.

With the reason that SR unit can only settle the small signal with small frequency and amplitude, so the IF signal with significant frequency and amplitude must be transacted to a small one, which is carried out by the ST unite. Therefore, we introduce the variable substitution [9] ; equation (3) can be transacted to

So, the frequency scale transacted equation is , where the frequency is transacted to the times of the original. The amplitude scale transacted equation is . When the parameter of a is large enough and b is small enough, the large IF signal can be transacted to small signal after ST. At the same time, the noise intensity becomes .

Therefore, the corresponding FPE of equation (4) can be described as

where is the PDF of the particle locating at time t.

As for the formal unification, we rewrite equation (5) as

where and .

Next, we will analyze the MFPT with time-vary item of equation (6).

3. The MFPT Influenced by DSFH Signal and Gaussian White Noise

3.1. The Nonautonomous Effect of MFPT Caused by a Time-Vary Input Signal

Let represent the probability of the first passage time from the well to the well is T and the initial position is x. So, MFPT is the average of T as

For simplicity, we use to represent the . So, the probability of the first passage time among can be described as

where is the probability of the particles positing at and at time t.

And equation (8) can be simplified as

Substituting equation (9) to equation (7), we obtain

So, using the FPE, we can obtain

Integrating t on both sides of equation (11), it shows

When the input signal is a sinusoidal wave, the time-vary item can be decomposed into the nontime-vary item representing the SR system and time-vary item representing the input signal as

Then, equation (12) can be expressed as

Substituting equation (7) to equation (14), we can obtain that

Next, we emphasize on the nonautonomous item and obtain the solution of .

3.2. The Solution of MFPT with the Transaction from Autonomous Item to Nonautonomous One

The output of SR and the sketch of MFPT are described in Figure 3. We can see that the SR system reaches a steady state quickly in one period. The transaction between wells is so quick, that is, the MFPT is much smaller than the period.

Figure 3: The output of SR and the sketch of MFPT (SNR = −18 dB, a = 1e4, b = 2.59e12, f₀ = 1e3, and R = 200).

Based on this phenomenon that particles of the SR system jump in a potential well after completing the barrier transition and the electromagnetic particles jump so quickly, a given time of relaxation can represent the drift effect of . And as for the cosine form, the given time is . So, the nonautonomous item can be described as

Due to ( is the period of the output signal of the SR system), so . And in consideration of the experience of communication signals detection and simulation results, which can be viewed as small aqueous about the domain of 0. So, we apply the Taylor expansion to the with the whole of as (neglecting the phase effect)

So, equation (15) can be described as

Equation (18) is the nonhomogeneous differential equation with one unknown variable coefficient of second order, and there is no analytical solution. So, let , and equation (18) can be described as

So, we can obtain the numerical solution of with the method of Runge–Kutta.

As a reminder that when the input signal is a sine wave, the nonautonomous item is . And the steady well of SR lies in the peak or valley of the sine wave. And the given time of relaxation which can represent the drift effect of is the . So, the drift sine wave can be represented by . And we apply the Taylor expansion to the with the whole of as

The form of the differential equation drifted by sine wave is the same with the one drifted by cosine wave. The drift effect caused by the cosine wave is the same as the sine wave demonstrated in the mathematical analytic form.

4. Simulated Analysis

In this section, we set the optimal parameter of SR drawing lessons from the optimal matching SR and statistics for the first passage time of the particles jumping from one well to the other as the simulation results in the Simulink model. And we will analyze the theory and simulation in three aspects. (1) Parameters of a and b of the SR system and the scale transformation coefficient , which represent the intrinsic system nature together. (2) The signal amplitude A and the signal frequency of the receptive signal of DSFH. (3) The noise intensity D and the multiple sampling R.

4.1. The Relationship between MFPT and Signal Amplitude A

The theoretical and simulation results of the MFPT varying with A are described in Figure 4. We can see that the MFPT decreases as the SNR increases no matter in theoretical or simulation curve. This is because that when the signal frequency , sampling multiple R, and noise power (corresponding to noise variance) are constant, the frequency scale transformation a (one of the parameters of the SR system) and noise intensity D is constant. So, the three factors listed above which can affect the particles to jump vary as follows: the scale transformation coefficient and signal amplitude A increase, while the noise intensity D is constant as the SNR increases. And the scale transformation coefficient representing the system characteristic can promote the particles to jump between wells in the intrinsic nature of the SR system. Furthermore, the signal can pull the particles to jump between wells in external factors. So, it is more conducive to particle transitions that and A are increase, and that is why MFPT decreases as the SNR increases.

Figure 4: The MFPT vary with A (f₀ = 1e3, R = 200, a = 1e4, and  = 4).

4.2. The Relationship between MFPT and Noise Intensity D

The theoretical and simulation results of the MFPT varying with D are described in Figure 5. We can see that the MFPT decreases as the SNR increases no matter in theoretical or simulation curve. Because when the signal frequency , sampling multiple R, and signal amplitude A are constant, the frequency scale transformation a is constant. So, the three factors listed above which can affect the particles to jump vary as follows: the scale transformation coefficient increases and the noise intensity D decreases, while signal amplitude A is constant as the SNR increases. And representing the system characteristic can promote the particles to jump between wells in intrinsic nature of the SR system. Furthermore, the noise can pull the particles to jump between wells in external factors. So, increase has a positive effect on the decrease of MFPT, while D decrease has a negative effect on the dynamic characteristic of SR leading decrease of MFPT. And the overall effect depends on the relative result caused by system and noise. Furthermore, the MFPT decreases when the system parameters and noise intensity are listed in our experiment.

Figure 5: The MFPT vary with D (f₀ = 1e3, R = 200, a = 1e4, and A = 1).

4.3. The Relationship between MFPT and Sampling Multiple R

The theoretical and simulation results of the MFPT varying with R are described in Figure 6. We can see that there are two results. One is that the MFPT increases as the R increases when SNR is constant no matter in the theoretical or simulation curve. Because when SNR, , A and are constant, the parameters a, are constant, and D decreases with the increase of R. So, the positive effect for transaction bring by noise decreases becomes weak as the increase of R, leading the increase of MFPT. The other result is that the MFPT decreases as the SNR increases when R is constant no matter in the theoretical or simulation curve. Because when , R, and are constant, the parameters a and D are constant, while and A increase with increase of SNR. The effect caused by and A are all positive. So, the MFPT decreases. Furthermore, the bigger the R which leads higher quality of A/D converter, the worse the MFPT. However, the SR system vibrate needs a certain R, so for the dynamic system performance, we should decrease R meeting the demand of SR vibration as possible.

Figure 6: The MFPT vary with R (f₀ = 1e3, a = 1e4, b = (2.59e12 (SNR = −18 dB), 7.42e12 (SNR = −12 dB), and 18.47e12 (SNR = −6 dB)), and  = 4).

As a reminder, the bigger the , the bigger the error between theoretical and simulation for a constant SNR, which also verify our assumption that the smaller the , the more accuracy of the Taylor expansion to the with the whole of . And different SNRs need different system parameters a and b, so the error between theoretical and simulation at different SNRs is mainly due to the system characteristic.

4.4. The Relationship between MFPT or Duty Cycle and the IF Frequency

The theoretical and simulation results of the MFPT varying with are described in Figure 7(a) and the duty cycle varying with are described in Figure 7(b).

(a)

(b)

Figure 7: The MFPT vary with f₀ (R = 200 and  = 4).

Two results are shown in Figure 7(a). One is that the MFPT decreases as the increases when SNR is constant, no matter in the theoretical or simulation curve. This is because that when SNR, R, A, and are constant, the parameters a and D decrease with the increase of . The overall effect of the positive one for transaction caused by system parameters and a negative one for transaction caused by noise is positive, leading to a decrease of MFPT. The other one is that the MFPT decreases as the SNR increases when is constant no matter in the theoretical or simulation curve. This is because that when , R, and are constant, the parameters a and D are constant, while and A increase with increase of SNR. The effect of and A are all positive. So, the MFPT decreases.

In Figure 7(b), we give the normalization of MFPT as varies. Because the output of SR has a distinct periodic characteristic that reflects the frequency of the input signal, the normalization is the duty cycle which is . We can also see two results. One is that the duty cycle increases as the increases when SNR is constant no matter in the theoretical or simulation curve because the period of the SR system at constant SNR mainly depends on the input signal period. And its periodicity can react quickly following with the input signal. However, the MFPT reacts slower than periodicity, which leading the result that the duty cycle increases as the increases when SNR is constant. And the parameters of SR varies more, following the .

Furthermore, the tracking performance of SR following the , which explains why the SR system is limited to be applied to communication systems with varying. That is to say when SNR is low, the MFPT cannot follow the dynamic period of SR controlled by the . The other result is that the duty cycle decreases as the SNR increases when is constant in the theoretical or simulation curve. The reason is the same to the one in Figure 7(a), so there is no more detailed description. And related to the sample frequency, the frequency of the received signal locates at the low region of the frequency band, which can be processed by SR.

5. Conclusion

MFPT can represent the dynamic characteristic of SR, and the effect on MFPT caused by the DSFH signal is studied in this paper. Based on the assumption that MFPT is small aqueous about the domain of 0, the nonautonomous differential equation with MFPT is transformed to a nonhomogeneous differential equation with one unknown variable coefficient of second order. The numerical solution of MFPT is obtained by Runge–Kutta. The quantitative relationship between MFPT and DSFH signal, noise, and parameters of SR is obtained. Furthermore, the conclusion that “when SNR = −12 dB, R = 200, and  = 2100, the duty cycle reaches 25% (available for DSFH signal detection), so it is suggested to decrease the IF frequency of the receptive DSFH or increase the SNR for further availability.” The research method illuminates the way of effect on MFPT caused by input signal. Furthermore, the conclusion provides reference for the dynamic characteristic and signal detection of the SR system driven by the DSFH signal and Gaussian white noise.

Data Availability

The data used to support the findings of the study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

References

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.

X. Zhou, P. Kyritsi, P. C. F. Eggers, and F. H. P. Fitzik, ““The medium is the message”: secure communication via waveform coding in MIMO systems,” in Proceedings of the 2007 IEEE 65th Vehicular Technology Conference—VTC2007-Spring, pp. 491–495, Dublin, Ireland, April 2007.

H. Quan, H. Zhao, and P. Cui, “Anti-jamming frequency hopping system using multiple hopping patterns,” Wireless Personal Communications, vol. 81, no. 3, pp. 1159–1176, 2015.

H. Zhao, H. D. Quan, and C. Pei-Zhang, “Follower-jamming resistible multi-sequence frequency hopping wireless communication,” Systems Engineering and Electronics, vol. 3, pp. 671–678, 2015.

C. Du, H. Quan, P. Cui, W. Liang, P. Zhou, and J. Dou, “Carrier sense random packet CDMA protocol in dual-channel networks,” Radioengineering, vol. 24, no. 2, pp. 507–517, 2015.

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.

G. Liu, Y. Kang, H. Quan, H. Sun, P. Cui, and C. Guo, “The detection performance of the dual-sequence-frequency-hopping signal via stochastic resonance processing under color noise,” Radioengineering, vol. 27, no. 3, pp. 618–626, 2019.

R. Benzi, A. Sutera, and A. Vulpiani, “The mechanism of stochastic resonance,” Journal of Physics A: Mathematical and General, vol. 14, no. 11, pp. 453–457, 1981.

G. K. Liu, H. D. Quan, Y. M. Kang, H. X. Sun, P. Z. Cui, and Y. M. Han, “A quadratic polynomial receiving scheme for sine signals enhanced by stochastic resonance,” Acta Physica Sinica, vol. 68, no. 21, Article ID 210501, 2019.

L. Gammaitoni, P. Hänggi, P. Jung, and F. Marchesoni, “Stochastic resonance,” Reviews of Modern Physics, vol. 70, no. 1, pp. 223–287, 1998.

F. Chapeau-Blondeau and X. Godivier, “Theory of stochastic resonance in signal transmission by static nonlinear systems,” Physical Review E, vol. 55, no. 2, pp. 1478–1495, 1997.

J. Tougaard, “Signal detection theory, detectability and stochastic resonance effects,” Biological Cybernetics, vol. 87, no. 2, pp. 79–90, 2002.

Z. Qiao, Y. Lei, and N. Li, “Applications of stochastic resonance to machinery fault detection: a review and tutorial,” Mechanical Systems and Signal Processing, vol. 122, pp. 502–536, 2019.

C. W. Gardiner, Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences, Springer-Verlag, Berlin, Germany, 1985.

E. Guardia and M. San Miguel, “Escape time and state dependent fluctuations,” Physics Letters A, vol. 109, no. 1-2, pp. 9–12, 1985.

Y. M. Kang, J. X. Xu, and Y. Xie, “Observing stochastic resonance in an underdamped bistable Duffing oscillator by the method of moments,” Physical Review E, vol. 68, no. 3, Article ID 036123, 2003.

Y. M. Kang, X. Chen, X. D. Lin et al., “Mean first passage time and stochastic resonance in a transcriptional regulatory system with non-Gaussian noise,” Fluctuation and Noise Letters, vol. 16, no. 1, pp. 201–207, 2017.

K. K. Wang, Stochastic Dynamical Characteristics for a Few Class of Biological Population Models, Nanjing University of Aeronautics and Astronautics College of Aerospace Engineering, Nanjing, China, 2014.

K. K. Wang, X.-B. Liu, and J.-H. Yang, “The mean extinction time and stability for a metapopulation system driven by colored cross-correlated noises,” Acta Physica Sinica, vol. 63, no. 10, Article ID 100502, 2013.

P. Xu and Y. Jin, “Mean first-passage time in a delayed tristable system driven by correlated multiplicative and additive white noises,” Chaos, Solitons & Fractals, vol. 112, pp. 75–82, 2018.

P. Xu and Y. Jin, “Stochastic resonance in an asymmetric tristable system driven by correlated noises,” Applied Mathematical Modelling, vol. 77, pp. 408–425, 2020.