# Damping Coefficient Induces Stochastic Multiresonance in Bistable System with Asymmetric Dichotomous Noise

Authors: Juan Wu; Yong Xu
Journal: Discrete Dynamics in Nature and Society
 (2014)

---

## Abstract

Stochastic resonance (SR) and stochastic multiresonance (SMR) phenomena as a function of the underdamping and overdamping coefficients in bistable system with asymmetric dichotomous noise are investigated numerically. By the efficient numerical simulation of the asymmetric dichotomous noise and the fourth-order Runge-Kutta algorithm, we calculate the system responses, the averaged power spectrum, and the signal-noise-ratio (SNR) that can be a measure of the existence of SR and SMR phenomenon. And the effects of damping coefficients on the three characteristics are analyzed. Firstly, it is found that the periodic asymmetric distribution of the particle’s hopping between two potential wells in the system response is gradually weakened as underdamping coefficient is increased to overdamping coefficient. And it also displays the periodic asymmetric distribution under the circumstance of overdamping coefficient. Then the averaged power spectrum exhibits multiple sharp peaks, and the highest peak increases and decreases for underdamping coefficient which is added to overdamping coefficient. Finally, SNR versus the damping coefficient for the system parameters and the noise parameters are acquired and they show multiple peaks and valleys, which illustrates the obvious SMR phenomena in bistable system with asymmetric dichotomous noise.


---

## Body

## 1. Introduction

Previously, noise is regarded as an ingredient that has a negative effect on system. Nevertheless the appearance of SR phenomenon alters people’s attitude, which assists the sequence of the disordered system to be well organized. Given this, SR has been investigated profoundly by a number of scientists. SR phenomenon was firstly proposed to explain the periodicity in the earth’s ice ages [1]. Then SR phenomenon was observed in various fields in nature, for instance, physics, chemistry, biology, ecology, graphics, and so on [2]. Up to now, there are abundant publications on SR based on the original discovery of SR, and the conception of SR was developed in a broad sense, for example, SMR, quantum SR, autonomous SR, aperiodic SR, coherence resonance, and logical SR [3–16].SR and SMR are thought to be a kind of the coherent manner among nonlinear system, random noise, and periodic force. And various characteristics can be used to measure the emerging of SR and SMR phenomena, such as the amplitude of system response, the output amplitude gain, the averaged power spectrum, the spectral power amplification, the signal-noise-ratio, the residence time distributions, and the information and probability of detection, by which SR and SMR phenomena can be easily explored [2]. For example, Ray and Sengupta measured SR in underdamped bistable system by the power spectrum [17], Xu et al. found SR phenomena in a bistable system with Lévy noise by the SNR [18], and Zhang et al. demonstrated SMR phenomena in a linear system driven by multiplicative polynomial dichotomous noise by the SNR [19]. In this paper, the amplitude of system response, the averaged power spectrum, and the SNR are the main characteristic indicators to demonstrate the existence of SR and SMR phenomena.As we know, SR and SMR phenomena in overdamped bistable system and underdamped bistable system have been researched widely both theoretically and experimentally for their applications [17, 20–23]. However, fewer investigators pay attention to SR and SMR phenomena in the overdamped and underdamped bistable system simultaneously. Thus, we will focus on SR and SMR phenomena versus the damping coefficients. Besides, a mass of the researches focus on Gaussian noise or white noise just for its simplicity. But they are just a sort of ideal model of the actual noises, and they cannot represent noise with the exponential relevance. Nevertheless asymmetric dichotomous noise is a non-Gaussian color noise and is widespread for its manageability in some fields, which leaps between two fixed points (a,-b,(a,b>0)); its waiting time is submitted to the exponential distribution, and this kind of the leaps inspires the nonequilibrium more easily than Gaussian noise [24]. In addition, up to now SR and SMR phenomena about the dichotomous noise have been studied largely in theory. Jin, Li, and some other experts have done some great work theoretically on SR in typical systems with the dichotomous noise [25–27]. And the research of the numerical simulation of SR and SMR phenomena is less relative. Consequently in this paper we will explore SR and SMR phenomena versus the damping coefficients in bistable system with asymmetric dichotomous noise numerically.This paper is organized as the following. In Section2, the bistable system and asymmetric dichotomous noise are introduced. In Section 3, SR and SMR phenomena are researched in three aspects. In Section 3.1, SR phenomenon is investigated by the transition between the two potential wells of the response of the system driven by the asymmetric dichotomous noise. In Section 3.2, the averaged power spectrum is computed numerically, and the signature of SR phenomenon is obviously displayed. In Section 3.3, we obtain the SNR versus the damping coefficient. And there are multiple peaks in the SNR which demonstrates the existence of SMR phenomena in the bistable system with asymmetric dichotomous noise.


## 2. Bistable System with Asymmetric Dichotomous Noise

We consider a bistable system with a periodic signal driven by asymmetric dichotomous noise, which can be described by the following Langevin equation:(1)d2xdt2=-ηdxdt-dV(x)dx+F(t)+φ(t),
where η is the damping coefficient. It is an underdamped bistable system when η<1; it is an overdamped bistable system when η>1. V(x) is the double potential well function defined as -mx2/2+nx4/4 (m>0,n>0), which has two stable points at V±=±m/n, and the height of potential barrier is ΔV(x)=m2/4n. When we choose the parameters as m=n=1, it is the standard double potential well function. And the two stable fixed points are V±=±1, ΔV(x)=1/4. F(t) is the periodic signal, which can be described as Acos⁡(ωt+θ), and A and ω, respectively, are the amplitude and forcing frequency of the periodic signal.φ
(
t
) is asymmetric dichotomous noise which jumps between two values a and b with mean waiting times ta and tb, and the rates of the switching can be obtained from the mean waiting; that is μa=1/ta, μb=1/tb.The master equation of this noise can be described as(2)∂tP(a,t∣x,t0)=-μaP(a,t∣x,t0)+μbP(b,t∣x,t0),∂tP(b,t∣x,t0)=μaP(a,t∣x,t0)-μbP(b,t∣x,t0),
with the initial condition and the total probability condition:
(3)P(x′,t∣x,t0)=δx′x,P(a,t∣x,t0)+P(b,t∣x,t0)=1.
Then the solution of the master equation is
(4)P(a,t∣x,t0)=μbμb+μa+(μaμb+μaδax-μbμb+μaδbx)P(a,t∣x,t0)=×exp⁡(-(μa+μb)(t-t0)),P(b,t∣x,t0)=μaμb+μa-(μaμb+μaδax-μbμb+μaδbx)P(b,t∣x,t0)=×exp⁡(-(μa+μb)(t-t0)).
And the stationary solution of (2) can be easily obtained as
(5)Ps(a)=μbμa+μb,Ps(b)=μaμa+μb.
Then the stationary mean which can be obtained by using (4) is
(6)〈ξ(t)〉s=aμb+bμaμa+μb,
and the stationary correlation function is
(7)〈ξ(t)ξ(t′)〉=(bμa+aμbμb+μa)2+μaμb(a-b)2(μb+μa)2×exp⁡(-(μa+μb)(t-t′)).
Here the mean function and the correlation function satisfy the following conditions:
(8)〈ξ(t)〉=0,〈ξ(t)ξ(t′)〉=λDexp⁡(-λ|t-t′|),
where D is the noise intensity and λ is the noise correlation time. Thus, the noise intensity of the asymmetric dichotomous noise can be computed as
(9)D=12∫-∞∞(〈ξ(τ)ξ(0)〉-〈ξ(τ)〉2)dτ=(a-b)2μaμb(μa+μb)3.
Finally, the rates of the switching μa and μb of the asymmetric dichotomous noise are computed in terms of the state values and the noise intensity. Also the conditional probabilities Paa and Pba can also be computed as follows:
(10)Paa=P(a,tn+1∣a,tn)=μbμb+μa+μaμb+μa×exp⁡⁡(-(μa+μb)dt),Pba=P(a,tn+1∣b,tn)=μbμb+μa-μbμb+μa×exp⁡(-(μa+μb)dt).And the numerical series of the asymmetric dichotomous noise can be obtained by the above formulas and the following relevant algorithm. The algorithm procedure of the asymmetric dichotomous noise can be described as follows. Firstly, the initial state of the asymmetric dichotomous noise can be supposed to bex0=a, (x0=b); a series of random numbers Rn (n=0,1,2,…) in the interval [0,1] are generated in computer, which are compared with the conditional probability Paa or Pba. Then we consider that if R0<Paa (R0<Pba), we will ascertain x1=a, or else x1=b. Next, we also should decide that if R1<Paa (R1<Pba), we will ascertain x2=a, or else x2=b. And keep doing this. Finally, a series of random numbers Rn are got in computer. Moreover, the timestep dt should be much lesser [28].


## 3. Stochastic Multiresonance

We are devoted to researching SR and SMR phenomena of this bistable system induced by the change of the damping coefficients in this section. According to this complex two-dimensional system, numerical simulation method is a good research approach, so the specific numerical simulation program is shown in regard to SR and SMR in bistable system with asymmetric dichotomous noise. And we study the system responses, the averaged power spectrum, and the SNR which can be used to reflect SR phenomenon.For the system responses, the particle oscillates at the bottom of the one potential well at a level of the certain damping coefficients, and then it oscillates between the two potential wells with the increasing and decreasing of the damping coefficients. That is to say, SR phenomenon can be discovered with the damping coefficients of the system. And the averaged power spectrum displays the relationship between the system, the cosine signal, and the asymmetric dichotomous noise for different underdamping and overdamping coefficients, from which we can find SR phenomenon. Moreover, we also devote ourselves to computing the SNR as one important symbol of SR phenomenon. If the SNR gives rise to one or more extreme values as we modulate the damping coefficients from the underdamping to overdamping at a certain range of parameters, it shows that the bistable system with asymmetric dichotomous noise has presented SR phenomenon or SMR phenomena.


### 3.1. System Responses

In order to obtain the responses of this system, the bistable system is transformed into two one-dimensional systems as(11)dxdt=y,dydt=-ηy-dV(x)dx+F(t)+φ(t),
and then the system responses are calculated by the discretization of the above form and the fourth-order Runge-Kutta algorithm.We fix the parameters of the system, the signal, and the asymmetric dichotomous noise asA=0.2, θ=0.03, ω=0.04, D=0.01, a=3.0, b=-1.0. And the effects of the different damping coefficients including the underdamping coefficients and overdamping coefficients on the SNR are investigated. In Figure 1(a), the underdamping coefficients of this system is 0.2, and we can find that the particle oscillates in the V+=1 around, surmounts the potential barrier, jumps into and oscillates around the V-=-1. Then the particle oscillates in the V-=-1 around, surmounts the potential barrier, and oppositely jumps into and oscillates around the V+=1. For some time, the periodic transition of the particle between the two potential wells and the uniform asymmetry are obviously shown. With the increasing of the underdamping coefficients to 0.6, Figure 1(b) also clearly reflects the periodic transition of the particle between the two potential wells and the uniform asymmetry. When the underdamping coefficients of this system is increased to 1.0, the underdamped bistable system transforms into the overdamped one, the transition of the particle between the two potential wells and the asymmetry of the distribution can also be found but the clear periodicity is not occurred, and it is important that the particle remains mostly in the V-=-1. Finally, when the overdamping coefficients are continually increased to 1.4, compared with the case of Figure 1(c), the transition between the two potential wells still remains, and the most important thing is that the asymmetry of the distribution is more severe and the vast majority of the states are located in the V-=-1, yet the initial state is located in the V+=1 in Figure 1(d).The system responses of the underdamped and overdamped bistable system (a)η=0.2, (b) η=0.6, (c) η=1.0, and (d) η=1.4.
(a)
(b)
(c)
(d)According to three subfigures of Figure2, we pay attention to the overdamped bistable system and the asymmetry of the system response. In Figure 2(a), the overdamping coefficient is adjusted to 1.3, and the asymmetric dichotomous noise parameter is increased to a=4.0, b=-2.0. It is illustrated that the particle oscillates between the two potential wells, the intense asymmetry of the distribution exists, and most states are located in the V-=-1; meanwhile the random vibration is intense compared with Figure 1(d).   In Figure 2(b), we interchange the state value of the asymmetric dichotomous noise a=2.0,b=-4.0; it is clear that the asymmetry of the distribution and the action of the random vibration still exist; nevertheless most states are located in the V+=1 conversely, compared with Figure 2(a). Moreover, with the asymmetric dichotomous noise parameter decreased to a=3.0,b=-1.0, Figure 2(c) shows us the more regular asymmetry of the distribution as the action of the random vibration is weak.The system responses of the overdamped bistable system for different states values of the asymmetric dichotomous noise (a)a=4.0,b=-2.0, (b) a=2.0,b=-4.0, and (c) a=3.0,b=-1.0.
(a)
(b)
(c)


### 3.2. Averaged Power Spectrum

Power spectrum can reflect a kind of the coordination between the signal and the noise. But it includes a large number of random factors. So we make use of the method of average to eliminate the random factors and then obtain the averaged power spectrum which reflects the coefficient properties of the cosine signal and the asymmetric dichotomous noise. Accordingly, whether the damping coefficient including the underdamping and the overdamping can induce SR phenomenon in terms of the averaged power spectrum is the cure of our research in this section.The power spectrum density can be obtained by the following formula of the Fourier transform of the autocorrelation function:(12)P(ω)=∫-∞∞〈x(t+τ)x(t)〉exp⁡(-iωτ)dτ.
Next, it is found that ensemble averaging on 600 power spectrum trajectories is more sufficient to achieve the averaged power spectrum. Thus several averaged power spectrum figures are obtained as the following for the different damping coefficient.In Figure3(a), we choose the damping coefficient η=0.05, which means that the system becomes the underdamped bistable system with the fixed parameters A=1.0, a=3.0, b=-1.0, D=0.001, ω=1.0, θ=0.03. And it is easily observed that three distinct peaks appear on the averaged power spectrum due to the effect of the periodicity of the periodic signal and the asymmetry of the asymmetric dichotomous noise on the system, although the underdamping coefficient is smaller. They are marked P1, P2, and P3 from left to right so as to reveal the meaning that the figure contains conveniently and particularly. Between the three peaks, the value of the middle peak P2 is the highest; it is almost double P3’s, and the value of the left peak P1 is the lowest. Then the underdamping coefficient is added to 0.1 as other parameters are invariant, and three distinct peaks are still observed in Figure 3(b). Here the value of the highest peak P2 is rapidly increased, but the value of P3 is decreased and the value of P1 is decreased slightly. And as it is increased to 0.3, 0.6, 0.9, 1.0, and 1.3, the similar increase and decrease are clearly present in Figures 3(c), 3(d), 3(e), 3(f), and 3(g).The averaged power spectrum of the system for different damping coefficient (a)η=0.05, (b) η=0.1, (c) η=0.3, (d) η=0.6, (e) η=0.9, (f) η=1.0, and (g) η=1.3.
(a)
(b)
(c)
(d)
(e)
(f)
(g)In general, when the damping coefficient is increased gradually from the underdamping to overdamping, the value ofP2 that is the highest peak increases sharply and then decreases slowly, the value of P3 that is the second highest peak decreases all the time, and the value of P1 that is the lowest peak decreases firstly and then increases slowly and disappears last. And all those variations of the three peaks with the change of damping coefficient reveal SR phenomena.


### 3.3. Signal-Noise-Ratio

SNR is a typical method to measure SR and SMR phenomena. And in this section as a function of damping coefficient it displays the conspicuous SR and SMR phenomena. At present, there are many several numerical simulation methods about SNR. We employ the following formula [29]:
(13)SNR=10log⁡PS(ω)PN(ω),
where PS(ω) and PN(ω) are the output power spectrum of the periodic signal and the asymmetric dichotomous noise, respectively. And in the following figures, the peaks of the SNR phenomenon are marked as S1, S2, S3, and S4 similarly in order to represent those figures conveniently and particularly.In Figure4, SNR as the functions of the damping coefficient for the different amplitudes of signal and the fixed system A=1.0,0.7,0.4, a=3.0, b=-1.0, D=0.001, ω=1.0, θ=0.03, displays the conspicuous SR and SMR phenomena. It is clearly showed that, as the damping coefficient increases from 0.04 to 1.4, the SNR firstly increases sharply, next decreases slowly, and then increases more slowly. The nonmonotonic behaviors of the SNR obviously reveal the occurrence of SMR phenomenon. At the same time, the value of the peak descends slowly and then rapidly; also the value of the valley descends always and it moves towards the left, when the amplitude of the signal decreases from 1.0 to 0.7 and then to 0.4.Figure 4
SNR versus the damping coefficient for the different amplitudes of the signal.In Figure5, we observe the SNR versus the damping coefficient for the different forcing frequency of the periodic signal, and four subfigures represent the changes of SNR with the fixed system parameter A=1.0, a=3.0, b=-1.0, D=0.001, θ=0.03. In Figure 5(a), the obvious SR phenomenon is easily observed, but SR phenomenon demonstrates to us that there are two peaks; in other words, it is SMR phenomenon. And as the forcing frequency of the periodic signal ω is reduced from 1.5 to 1.4 and to 1.3, the peaks S1 and S2 heighten and move right gradually and they become more and more evident, while it barely heighten when ω is reduced from 1.4 to 1.3 according to S2. And when it is decreased from 1.3 to 0.7, the similar move, increase, and decrease of the peaks are present, respectively, in Figures 5(b), 5(c), and 5(d).SNR versus the damping coefficient for the different forcing frequencies of the periodic signal.
(a)
(b)
(c)
(d)In short, as the forcing frequency of the periodic signal constantly lessens from 1.5 to 0.7, the SNR demonstrate the variation between two peaks and one peak, which can declare that the noteworthy SMR phenomena exist in the bistable system with asymmetric dichotomous noise.In Figure6, similarly we focus on the SNR versus the damping coefficient for the different state a of the asymmetric dichotomous noise with the fixed system parameters A=1.0, θ=0.03, ω=1.3, b=-1.0, D=0.001. According to the states of the asymmetric dichotomous noise, the state b is fastened and the state a is changed continually. When the state a is selected as a=4.0 and the other state is b=-1.0, the noise is still the asymmetric dichotomous noise. It can be easily found that there are two peaks in the SNR, which go up increasingly but do not move toward left and right, when the state a is reduced to 3.0 and then to 2.0. However, three distinct peaks turn up as the state a is changed to 0.5. During the process of the decrease of the state a, the SNR show us two peaks and three peaks, which demonstrates SMR phenomena.Figure 6
SNR versus the damping coefficient for the different state values of the asymmetric dichotomous noise.The effect of the different asymmetric states of dichotomous noise on the SNR versus the damping coefficient is just investigated for SMR phenomena. Now we research how the symmetric states of dichotomous noise influence SMR phenomena versus the damping coefficient. And the parameter of the periodic signal and the noise is fastenedA=1.0, θ=0.03, ω=1.3, D=0.001. By the clear manifestation of Figure 7, it is found that the SNR versus the damping coefficient present multiple peaks with the different symmetric states of dichotomous noise. When it is a=3.0b=-3.0, there are a peak and a valley. As the states are lessened to a=2.0b=-2.0, the peak and the valley rise up but their positions do not shift to the left or the right; meanwhile another peak appears in the SNR. When the states are continuously reduced to a=1.0b=-1.0, the peak S1 and the valley go up and similarly their positions do not shift left and right; the peak S2 yet develops into two peaks. And then the three peaks S1, S2, and S3 get more remarkable, when the states are lessened lower a=0.5b=-0.5. So SMR phenomena are always present, when the state of the symmetric dichotomous noise is increased constantly.Figure 7
SNR versus the damping coefficient for the different state values of the symmetric dichotomous noise.Finally, whether the noise intensity of the asymmetric dichotomous noise can cause the prominent influence on the peak of the SNR is researched. Figure8 shows us how the SNR changes based on the fixed parameters A=1.0, θ=0.03, ω=1.3, a=0.5, b=-1.0 with the increase of the noise intensity. It is demonstrated that three noteworthy peaks S1, S2, and S3 appear in the SNR, when the noise intensity is 0.001. And as it is increased to 0.01, the peak S2 disappears, and the peaks S1, S3 descend. There are two peaks in the SNR, which are marked again as S1, S2. Then we continue to improve the noise intensity to 0.2 and find that the two peaks S1 and S2 both decline gradually. Until the noise intensity is increased to 0.3, the peaks S1 and S2 decline gradually once again, but it is particular that the peak S2 disappears. There are a peak and a valley in the SNR with D=0.03. And all give evidence of the existence of SMR phenomena.Figure 8
SNR versus the damping coefficient with the different noise intensity of the asymmetric dichotomous noise.Furthermore, some situations in Figures6, 7, and 8 are different from Figures 4 and 5’s, which demonstrate to us that the amplitude and the forcing frequency of the periodic signal can change the position toward left and right of the peaks of SNR versus the damping coefficient, while the states and the intensity of the asymmetric dichotomous noise do not change the position of the peaks toward left or right but the high of the peaks.


## 3.1. System Responses

In order to obtain the responses of this system, the bistable system is transformed into two one-dimensional systems as(11)dxdt=y,dydt=-ηy-dV(x)dx+F(t)+φ(t),
and then the system responses are calculated by the discretization of the above form and the fourth-order Runge-Kutta algorithm.We fix the parameters of the system, the signal, and the asymmetric dichotomous noise asA=0.2, θ=0.03, ω=0.04, D=0.01, a=3.0, b=-1.0. And the effects of the different damping coefficients including the underdamping coefficients and overdamping coefficients on the SNR are investigated. In Figure 1(a), the underdamping coefficients of this system is 0.2, and we can find that the particle oscillates in the V+=1 around, surmounts the potential barrier, jumps into and oscillates around the V-=-1. Then the particle oscillates in the V-=-1 around, surmounts the potential barrier, and oppositely jumps into and oscillates around the V+=1. For some time, the periodic transition of the particle between the two potential wells and the uniform asymmetry are obviously shown. With the increasing of the underdamping coefficients to 0.6, Figure 1(b) also clearly reflects the periodic transition of the particle between the two potential wells and the uniform asymmetry. When the underdamping coefficients of this system is increased to 1.0, the underdamped bistable system transforms into the overdamped one, the transition of the particle between the two potential wells and the asymmetry of the distribution can also be found but the clear periodicity is not occurred, and it is important that the particle remains mostly in the V-=-1. Finally, when the overdamping coefficients are continually increased to 1.4, compared with the case of Figure 1(c), the transition between the two potential wells still remains, and the most important thing is that the asymmetry of the distribution is more severe and the vast majority of the states are located in the V-=-1, yet the initial state is located in the V+=1 in Figure 1(d).The system responses of the underdamped and overdamped bistable system (a)η=0.2, (b) η=0.6, (c) η=1.0, and (d) η=1.4.
(a)
(b)
(c)
(d)According to three subfigures of Figure2, we pay attention to the overdamped bistable system and the asymmetry of the system response. In Figure 2(a), the overdamping coefficient is adjusted to 1.3, and the asymmetric dichotomous noise parameter is increased to a=4.0, b=-2.0. It is illustrated that the particle oscillates between the two potential wells, the intense asymmetry of the distribution exists, and most states are located in the V-=-1; meanwhile the random vibration is intense compared with Figure 1(d).   In Figure 2(b), we interchange the state value of the asymmetric dichotomous noise a=2.0,b=-4.0; it is clear that the asymmetry of the distribution and the action of the random vibration still exist; nevertheless most states are located in the V+=1 conversely, compared with Figure 2(a). Moreover, with the asymmetric dichotomous noise parameter decreased to a=3.0,b=-1.0, Figure 2(c) shows us the more regular asymmetry of the distribution as the action of the random vibration is weak.The system responses of the overdamped bistable system for different states values of the asymmetric dichotomous noise (a)a=4.0,b=-2.0, (b) a=2.0,b=-4.0, and (c) a=3.0,b=-1.0.
(a)
(b)
(c)


## 3.2. Averaged Power Spectrum

Power spectrum can reflect a kind of the coordination between the signal and the noise. But it includes a large number of random factors. So we make use of the method of average to eliminate the random factors and then obtain the averaged power spectrum which reflects the coefficient properties of the cosine signal and the asymmetric dichotomous noise. Accordingly, whether the damping coefficient including the underdamping and the overdamping can induce SR phenomenon in terms of the averaged power spectrum is the cure of our research in this section.The power spectrum density can be obtained by the following formula of the Fourier transform of the autocorrelation function:(12)P(ω)=∫-∞∞〈x(t+τ)x(t)〉exp⁡(-iωτ)dτ.
Next, it is found that ensemble averaging on 600 power spectrum trajectories is more sufficient to achieve the averaged power spectrum. Thus several averaged power spectrum figures are obtained as the following for the different damping coefficient.In Figure3(a), we choose the damping coefficient η=0.05, which means that the system becomes the underdamped bistable system with the fixed parameters A=1.0, a=3.0, b=-1.0, D=0.001, ω=1.0, θ=0.03. And it is easily observed that three distinct peaks appear on the averaged power spectrum due to the effect of the periodicity of the periodic signal and the asymmetry of the asymmetric dichotomous noise on the system, although the underdamping coefficient is smaller. They are marked P1, P2, and P3 from left to right so as to reveal the meaning that the figure contains conveniently and particularly. Between the three peaks, the value of the middle peak P2 is the highest; it is almost double P3’s, and the value of the left peak P1 is the lowest. Then the underdamping coefficient is added to 0.1 as other parameters are invariant, and three distinct peaks are still observed in Figure 3(b). Here the value of the highest peak P2 is rapidly increased, but the value of P3 is decreased and the value of P1 is decreased slightly. And as it is increased to 0.3, 0.6, 0.9, 1.0, and 1.3, the similar increase and decrease are clearly present in Figures 3(c), 3(d), 3(e), 3(f), and 3(g).The averaged power spectrum of the system for different damping coefficient (a)η=0.05, (b) η=0.1, (c) η=0.3, (d) η=0.6, (e) η=0.9, (f) η=1.0, and (g) η=1.3.
(a)
(b)
(c)
(d)
(e)
(f)
(g)In general, when the damping coefficient is increased gradually from the underdamping to overdamping, the value ofP2 that is the highest peak increases sharply and then decreases slowly, the value of P3 that is the second highest peak decreases all the time, and the value of P1 that is the lowest peak decreases firstly and then increases slowly and disappears last. And all those variations of the three peaks with the change of damping coefficient reveal SR phenomena.


## 3.3. Signal-Noise-Ratio

SNR is a typical method to measure SR and SMR phenomena. And in this section as a function of damping coefficient it displays the conspicuous SR and SMR phenomena. At present, there are many several numerical simulation methods about SNR. We employ the following formula [29]:
(13)SNR=10log⁡PS(ω)PN(ω),
where PS(ω) and PN(ω) are the output power spectrum of the periodic signal and the asymmetric dichotomous noise, respectively. And in the following figures, the peaks of the SNR phenomenon are marked as S1, S2, S3, and S4 similarly in order to represent those figures conveniently and particularly.In Figure4, SNR as the functions of the damping coefficient for the different amplitudes of signal and the fixed system A=1.0,0.7,0.4, a=3.0, b=-1.0, D=0.001, ω=1.0, θ=0.03, displays the conspicuous SR and SMR phenomena. It is clearly showed that, as the damping coefficient increases from 0.04 to 1.4, the SNR firstly increases sharply, next decreases slowly, and then increases more slowly. The nonmonotonic behaviors of the SNR obviously reveal the occurrence of SMR phenomenon. At the same time, the value of the peak descends slowly and then rapidly; also the value of the valley descends always and it moves towards the left, when the amplitude of the signal decreases from 1.0 to 0.7 and then to 0.4.Figure 4
SNR versus the damping coefficient for the different amplitudes of the signal.In Figure5, we observe the SNR versus the damping coefficient for the different forcing frequency of the periodic signal, and four subfigures represent the changes of SNR with the fixed system parameter A=1.0, a=3.0, b=-1.0, D=0.001, θ=0.03. In Figure 5(a), the obvious SR phenomenon is easily observed, but SR phenomenon demonstrates to us that there are two peaks; in other words, it is SMR phenomenon. And as the forcing frequency of the periodic signal ω is reduced from 1.5 to 1.4 and to 1.3, the peaks S1 and S2 heighten and move right gradually and they become more and more evident, while it barely heighten when ω is reduced from 1.4 to 1.3 according to S2. And when it is decreased from 1.3 to 0.7, the similar move, increase, and decrease of the peaks are present, respectively, in Figures 5(b), 5(c), and 5(d).SNR versus the damping coefficient for the different forcing frequencies of the periodic signal.
(a)
(b)
(c)
(d)In short, as the forcing frequency of the periodic signal constantly lessens from 1.5 to 0.7, the SNR demonstrate the variation between two peaks and one peak, which can declare that the noteworthy SMR phenomena exist in the bistable system with asymmetric dichotomous noise.In Figure6, similarly we focus on the SNR versus the damping coefficient for the different state a of the asymmetric dichotomous noise with the fixed system parameters A=1.0, θ=0.03, ω=1.3, b=-1.0, D=0.001. According to the states of the asymmetric dichotomous noise, the state b is fastened and the state a is changed continually. When the state a is selected as a=4.0 and the other state is b=-1.0, the noise is still the asymmetric dichotomous noise. It can be easily found that there are two peaks in the SNR, which go up increasingly but do not move toward left and right, when the state a is reduced to 3.0 and then to 2.0. However, three distinct peaks turn up as the state a is changed to 0.5. During the process of the decrease of the state a, the SNR show us two peaks and three peaks, which demonstrates SMR phenomena.Figure 6
SNR versus the damping coefficient for the different state values of the asymmetric dichotomous noise.The effect of the different asymmetric states of dichotomous noise on the SNR versus the damping coefficient is just investigated for SMR phenomena. Now we research how the symmetric states of dichotomous noise influence SMR phenomena versus the damping coefficient. And the parameter of the periodic signal and the noise is fastenedA=1.0, θ=0.03, ω=1.3, D=0.001. By the clear manifestation of Figure 7, it is found that the SNR versus the damping coefficient present multiple peaks with the different symmetric states of dichotomous noise. When it is a=3.0b=-3.0, there are a peak and a valley. As the states are lessened to a=2.0b=-2.0, the peak and the valley rise up but their positions do not shift to the left or the right; meanwhile another peak appears in the SNR. When the states are continuously reduced to a=1.0b=-1.0, the peak S1 and the valley go up and similarly their positions do not shift left and right; the peak S2 yet develops into two peaks. And then the three peaks S1, S2, and S3 get more remarkable, when the states are lessened lower a=0.5b=-0.5. So SMR phenomena are always present, when the state of the symmetric dichotomous noise is increased constantly.Figure 7
SNR versus the damping coefficient for the different state values of the symmetric dichotomous noise.Finally, whether the noise intensity of the asymmetric dichotomous noise can cause the prominent influence on the peak of the SNR is researched. Figure8 shows us how the SNR changes based on the fixed parameters A=1.0, θ=0.03, ω=1.3, a=0.5, b=-1.0 with the increase of the noise intensity. It is demonstrated that three noteworthy peaks S1, S2, and S3 appear in the SNR, when the noise intensity is 0.001. And as it is increased to 0.01, the peak S2 disappears, and the peaks S1, S3 descend. There are two peaks in the SNR, which are marked again as S1, S2. Then we continue to improve the noise intensity to 0.2 and find that the two peaks S1 and S2 both decline gradually. Until the noise intensity is increased to 0.3, the peaks S1 and S2 decline gradually once again, but it is particular that the peak S2 disappears. There are a peak and a valley in the SNR with D=0.03. And all give evidence of the existence of SMR phenomena.Figure 8
SNR versus the damping coefficient with the different noise intensity of the asymmetric dichotomous noise.Furthermore, some situations in Figures6, 7, and 8 are different from Figures 4 and 5’s, which demonstrate to us that the amplitude and the forcing frequency of the periodic signal can change the position toward left and right of the peaks of SNR versus the damping coefficient, while the states and the intensity of the asymmetric dichotomous noise do not change the position of the peaks toward left or right but the high of the peaks.


## 4. Discussion and Conclusion

Stochastic resonance (SR) and stochastic multiresonance (SMR) phenomena versus the damping coefficient in bistable system with asymmetric dichotomous noise have been researched numerically in this paper. The system response, the averaged power spectrum, and the signal-noise-ratio (SNR) have been applied to investigate and demonstrate SR and SMR phenomena.Firstly, by the fourth-order Runge-Kutta numerical algorithm, it is found that the asymmetric dichotomous noise can induce the uniform asymmetry and the irregular asymmetry of the system response in the bistable system with the appropriately fixed parameters, as the damping coefficient is increased gradually from the underdamping 0.2 to the overdamping 1.4. Also the uniform asymmetry can be shown in the system response of the overdamped system, after the states values parameters of the asymmetric dichotomous noise are adjusted properly. Then we obtain the averaged power spectrum by the Fourier transform of the autocorrelation function. It is observed that there are three obvious peaks in the averaged powering spectrum. And when the damping coefficient is increased gradually from the underdamping to overdamping, the three peaks generate various transformations. The two parts above reveal SR phenomena in bistable system with asymmetric dichotomous noise.Finally, the SNR versus the damping coefficient is researched. It is found that several peaks appear in the SNR under some circumstances. The gradual augment of the amplitude of the periodical signal can make those peaks increase and the position toward the right. The decrease of the forcing frequency of the periodical signal can induce that the two peaks of the SNR rise and then decline and one peak disappears in the end; meanwhile the position of the peaks move toward the right. And one peak of the SNR develops to two and three peaks with the decrease of the states of the asymmetric and symmetric dichotomous noise. Also the reduction of the noise intensity of the asymmetric dichotomous noise can arouse two or three peaks of the SNR. Furthermore, the states and the intensity of the asymmetric dichotomous noise do not make the position of the peaks of the SNR move toward left and right. And this various situations of the peaks of the SNR versus the damping coefficient demonstrate the existence of SMR phenomena in the bistable system with asymmetric dichotomous noise.



---

*Source: 850361-2014-10-19.xml*