Stochastic Hopf–Hopf bifurcation of two-species discrete coupling logistic system with symbiotic interaction

In this paper, stochastic Hopf–Hopf bifurcation of the discrete coupling logistic system with symbiotic interaction is investigated. Firstly, orthogonal polynomial approximation of discrete random function in the Hilbert spaces is applied to reduce the discrete coupling logistic system with random parameter to the deterministic equivalent system. Then, it is concluded that Hopf–Hopf bifurcation exists in the equivalent deterministic system according to the principle of algebraic criteria. Numerical simulations show that the bifurcation critical value varies with the intensity of random parameter, and Hopf–Hopf bifurcation and period-doubling bifurcation behavior exist. In particular, Hopf–Hopf bifurcation can be drift with the change of random intensity, and frequency locking phenomenon occurs in the stochastic system.


Introduction
In 1798, Malthus proposed the population growth model which stated that the world population grew geometrically as the food and space resources became more abundant. The mathematical ecologist May suggested in one of his influential articles published in 1976 by Nature [1] that the non-heterogamous insect population model in ecology can be explained in terms of the nonlinear difference equation We call Eq. (1) logistic map or logistic system. As we all know, Eq. (1) and the form of generalization were studied in detail in reference [2]. After that many researchers have studied the logistic system in mathematical ecology, social science, economic net, and so on [3][4][5][6][7]. At the same time people studied the logistic system further in more aspects. So, some references [8][9][10] have paid more attention to the dynamical behavior in the logistic system, such as bifurcation, chaos, and stability. Recently, the dynamical behavior in many complicated logistic systems, such as coupling logistic system, has been explored.
Masoaller, Cavalcante, and Leite studied the delayed coupling of logistic maps [11]. The reference [12] analyzed synchronous chaos in the coupled system of two logistic maps. The symmetry breaking bifurcations in a circular chain of N coupled logistic maps have been presented in reference [13]. Zhang and Wang studied the spatiotemporal chaos in mixed linear-nonlinear coupled logistic map [14]. Leonel Rocha and Abdel-Kaddous Taha studied the Allee effect on the dynamical behavior of a new class of generalized logistic maps [15]. Duygu Arugaslan and Leyla Guzel considered the logistic equation with piecewise constant argument of generalized type and analyzed the stability of the trivial fixed point and the positive fixed point after reducing the equation into a nonautonomous difference equation [16]. The dynamical properties of a stochastic susceptible-infected epidemic model with logistic growth were investigated in reference [17]. Qifa Lin performed stability analysis of a single species logistic model with Allee effect and feedback control [18]. On the other hand, these results are limited to a deterministic logistic system. As is known to all, in the real world, uncertainty of the system is inevitable, which can affect the dynamical behavior of the nonlinear system. Therefore the influence of random factors on the logistic system are worthy of exploring. The effects of randomness on chaos and order of coupled logistic maps have been stated by Marcelo [19]. The stochastic resonance problem for a class of time-delay logistics systems has been discussed [20]. Xing and Peng studied the boundedness, persistence, and extinction of a stochastic nonautonomous logistic system with time delays [21]. The reference [22] studied the dynamical behaviors of stochastic delay logistic system with impulsive toxicant in polluted environment. The stability of a stochastic logistic model with distributed delay has been explored in reference [23]. Xu Yong and Ma Shaojuan studied Hopf bifurcation control in a nonlinear stochastic dynamical system with nonlinear random feedback method [24]. Xu Yong et al. performed stability analysis in a stochastic Brusselator system with random parameter [25]. The reference [26] discussed a sliding mode control (SMC) for a wing model excited by harmonic forces and colored Gaussian noise together. JinZhong MA et al. introduced a periodically driven bistable eutrophication model with Gaussian white noise as a prototype class of real systems; they found that stochastic perturbations and periodic excitations were also able to slow down an imminent critical transition [27]. Qi Liu et al. analyzed the influences of random fluctuations on a two-degree-of-freedom (TDOF) airfoil model with viscoelastic terms [28]. Zhang Xiaoyu et al. considered the tipping behaviors of the thermoacoustic system with the time-varying parameters and the combined excitations of additive and multiplicative colored noises [29]. Although the research results on stochastic logistics system have been obtained, there are few studies on the influence of uncertainty on bifurcation. In addition, as far as we know, the influence of internal random parameters on the dynamic behavior of logistics system has not been found in the existing literature. Motivated by the above discussions, this paper takes the coupled logistic system, which has a role of cohesion between one-dimensional and high-dimensional logistic model, as an example. By using the statistical characteristic of random variable, we build a stochastic coupled logistic system with random parameter. The influence of a random parameter in the two-species coupled logistic system with symbiotic interaction on the Hopf-Hopf bifurcation is studied by orthogonal polynomial approximation [30][31][32][33][34][35][36].
This paper is organized as follows. In Sect. 2, we transform the stochastic two-species coupling logistic system with symbiotic interaction with random parameter into its equivalent deterministic one by orthogonal polynomial approximation. We analyze the Hopf-Hopf bifurcation of the stochastic coupling logistic system in Sect. 3. The numerical simulations of Hopf-Hopf bifurcation about the two-species coupling logistic system with symbiotic interaction are shown in Sect. 4. Finally, conclusions are drawn in Sect. 5.

Stochastic coupling logistic system with random parameter and its orthogonal polynomial approximation
Consider a discrete coupling logistic model for the symbiotic interaction of the dynamic system [37] ⎧ ⎨ ⎩ x n+1 = μ(y n )x n (1x n ), where x n , y n are symbiotic species, the symbiotic interaction leads to both species that the growth rate μ(z) is varying with time. As it is seen in the equations, we call Eq. (2) symmetrical interaction. And it is shown that the particular dynamics of each species is a logistic map whose parameter μ n is not fixed, x n+1 = μx n (1x n ), but it is forced to remain in the interval [38][39][40][41]. The existence of a nontrivial fixed point at each step n ensures the nontrivial evolution of the system [42] where u is a positive constant, and we call that the mutual benefit. At the same time,ū is a random parameter which can be described as where u is the deterministic system parameter ofū, δ is regarded as the intensity of random disturbance, k is a random variable which obeys density function p k with standard deviation λ. By Eq. (2) and Eq. (3), we obtain Obviously, there are many fixed points in this discrete coupling logistic system as we choose reasonable system variables. Origin P 0 (0, 0) is a fixed point, other fixed points are: . In order to compute simply, we select origin P 0 (0, 0) to analyze the Hopf-Hopf bifurcation. So it follows from the orthogonal polynomial approximation that the response of Eq. (5) can be expressed by the following Fourier series under the condition of convergence in mean square: where is the ith standard orthogonal polynomial, M represents the largest order of the polynomial we have taken.
Substituting Eq. (6) and Eq. (4) into Eq. (5), we obtain With the aid of a cycle recurrence formula of orthogonal polynomial [43] kP λ the nonlinearity term in the right equation of Eq. (7) can be respectively written as follows: and and The stochastic term in the right equation of Eq. (7) can be respectively written as follows: and and and and and and where S i (n), W i (n), X i (n), Z i (n), and Y i (n) stand for the linear combination of nonlinearity terms calculated by computer algebraic, such as Maple. According to Eqs. (9)- (20), we obtain that Eq. (7) can be further reduced to where after standardization of polynomial. According to the statistical characteristics and orthogonal polynomial approximation of discrete random function in the Hilbert spaces and the orthogonality of Charlier orthogonal polynomials, we can finally get the equivalent deterministic logistic equation. As M → ∞, the discrete coupling stochastic logistic system is strictly equivalent to system (21) in the sense of mean square convergence. In order to facilitate the numerical analysis of this paper, we select M = 1, λ = 0.5 and approximately obtain the equivalent deterministic system , y 0 (n + 1) = 3uS 0 (n) + uy 0 (n) + 3 2 δ[S 1 (n) + S 0 (n)] + 1 2 δ[y 1 (n) + y 0 (n)] -3uZ 0 (n) -uY 0 (n) - 3 2 Then the approximate random response of the original stochastic logistic system can be expressed as and as k = 0, the sample response of mean parameter system (SMR) and the ensemble mean response of it (EMR) are calculated as follows: In this paper, we take the initial conditions of deterministic equivalent system (22) and the same initial conditions of the deterministic system So we take

Hopf bifurcation analysis
We firstly introduce the Hopf-Hopf bifurcation conditions [44] about the deterministic discrete system as follows. F(x, μ), x ∈ R n , μ ∈ R m , if the eigenvalues of a Jacobian matrix at the bifurcation parameter point μ = μ c , we have the following properties:
Let us assume that the solutions to the deterministic equivalent discrete coupling logistic system undergo a Hopf-Hopf bifurcation on some submanifold in parameter space corresponding to a critical valueμ =μ c .

The numerical simulation and numerical analysis
In this section, when the strength of random intensity δ = 0.000, the discrete coupling stochastic logistic system with symbiotic interaction is investigated. System (5) can be turned into a deterministic original coupling logistic system. We can know this bifurcation parameter u c = 1, system (5) undergoes a Hopf-Hopf bifurcation at the fixed point (0, 0). Figure 1 shows each of the bifurcation phenomena about the coupling stochastic logistic system with symbiotic interaction, when the random intensity δ = 0.000, 0.005, 0.020, 0.240. By the numerical analysis, the critical value for Hopf-Hopf bifurcation of the coupling stochastic logistic system with symbiotic interaction is u c = -δ + 1 -1 4 δ 2 . At the same time, the Hopf-Hopf bifurcation under the influence of different random strength is in advance.
When u = 0.965 and random intensity δ is equal to 0.000, 0.002, 0.006, 0.010, 0.050, 0.090, 0.150, 0.200, 0.240, respectively, the phase trajectories of DR and EMR all converge to zero. It is shown in Fig. 2(a) and Fig. 2(b) that fixed points are stable with discrete two-species coupled logistic system with symbiotic interaction, mean as δ = 0.000 with two-species discrete coupling stochastic logistic system with symbiotic interaction, and Fig. 2(b) is a partial enlarged effect of Fig. 2(a). We know that the two-species discrete coupling stochastic logistic system with symbiotic interaction does not undergo Hopf-Hopf bifurcation, which is illustrated by Fig. 2(a). Increasing the random intensity to δ = 0.002,

Figure 1
The bifurcation diagrams of deterministic and coupling stochastic logistic system the phase trajectories of system (22) that converge to limit circle are shown in Fig. 2(c), (d), and Fig. 2(d) is a partial enlarged effect of Fig. 2(c). From Fig. 2(c), (e), (g), (i), (j), we find that with the change of random intensity, the amplitude of limit circle is increased. From Fig. 2(k) we can see that the limit circle happened to intersect when random intensity increased to δ = 0.150. At the same time, from Fig. 2(l), (m), (n), Fig. 2(f ) is a partial enlarged effect of Fig. 2(e), and Fig. 2(h) is a partial enlarged effect of Fig. 2(g), Fig. 2(o) is a partial enlarged effect of Fig. 2(n). We can know that with the increase of random intensity the phase trajectories show random motion, and frequency lock phenomena appear at last. We find that with the change of bifurcation parameter the phase trajectories of the deterministic system accord with the phase trajectories of two-species discrete coupling stochastic logistic system with symbiotic interaction, the bifurcation happens in both systems. According to the above numerical analysis, we discover that the critical value for Hopf-Hopf bifurcation in the two-species discrete coupling stochastic logistic system with symbiotic interaction is varying from random intensity. Based on the theoretical analysis and numerical simulations, we find that the Hopf-Hopf bifurcation occurs in two-species discrete coupling stochastic logistic system with symbiotic interaction with the variation of bifurcation parameter. Compared with the deterministic system, the random intensity has an obvious influence on the bifurcation critical value of its coupled random logic system, and the bifurcation critical value decreases with the increase of random intensity, indicating that the Hopf-Hopf bifurcation point appears in advance and frequency locking phenomenon occurs. In other words, within a certain range of random intensity, two-species coupled logistic system with symbiotic interaction will never die out.

Conclusions
The orthogonal polynomial approximation theory of discrete random function is applied to research the Hopf-Hopf bifurcation of the two-species discrete coupling stochastic logistic system with symbiotic interaction. Analysis shows that the orthogonal polynomial the discussion of numerical analysis shows that Hopf-Hopf bifurcation point appears in advance with the change of random intensity.