This paper deals with the problem of limit cycles for the whirling pendulum equation x=y,y=sin x(cosx-r)under piecewise smooth perturbations of polynomials of cos x,sin x and y of degree n with the switching line x=0....This paper deals with the problem of limit cycles for the whirling pendulum equation x=y,y=sin x(cosx-r)under piecewise smooth perturbations of polynomials of cos x,sin x and y of degree n with the switching line x=0.The upper bounds of the number of limit cycles in both the oscillatory and the rotary regions are obtained using the Picard-Fuchs equations,which the generating functions of the associated first order Melnikov functions satisfy.Furthermore,the exact bound of a special case is given using the Chebyshev system.At the end,some numerical simulations are given to illustrate the existence of limit cycles.展开更多
The maximal number of limit cycles for a particular type Ⅲ system x = -y + lx2 + mxy, y =x(1 + ax + by) is studied and some errors that appeared in the paper by Suo Mingxia and Yue Xiting (Annals of Differential Equa...The maximal number of limit cycles for a particular type Ⅲ system x = -y + lx2 + mxy, y =x(1 + ax + by) is studied and some errors that appeared in the paper by Suo Mingxia and Yue Xiting (Annals of Differential Equations, 2003,19(3):397-401) are corrected. By translating the system to be considered into the Lienard type and by using some related properties, we obtain several theorems with suitable conditions coefficients of the system, under which we prove that the system has at most two limit cycles. The conclusions improve the results given in Suo and Yue's paper mentioned above.展开更多
Aris and Amundson studied a chemical reactor and obtained the two equationsDaoud showed that at most one limit cycle may exist in the region of interest. Itis showed in this paper that other singular points exist and ...Aris and Amundson studied a chemical reactor and obtained the two equationsDaoud showed that at most one limit cycle may exist in the region of interest. Itis showed in this paper that other singular points exist and that a stable limitt cycle existsaround the singularity (1/2, 2) when K∈(9-δ, 9).展开更多
This paper presents a novel mechanical attachment, i.e., nonlinear energy sink (NES), for suppressing the limit cycle oscillation (LCO) of an airfoil. The dynamic responses of a two-degree-of-freedom (2-DOF) air...This paper presents a novel mechanical attachment, i.e., nonlinear energy sink (NES), for suppressing the limit cycle oscillation (LCO) of an airfoil. The dynamic responses of a two-degree-of-freedom (2-DOF) airfoil coupled with an NES are studied with the harmonic balance method. Different structure parameters of the NES, i.e., mass ratio between the NES and airfoil, NES offset, NES damping, and nonlinear stiffness in the NES, are chosen for studying the effect of the LCO suppression on an aeroelastic system with a supercritical Hopf bifurcation or subcritical Hopf bifurcation, respectively. The results show that the structural parameters of the NES have different influence on the supercritical Hopf bifurcation system and the subcritical Hopf bifurcation system.展开更多
Based on the piston theory of supersonic flow and the energy method, the flutter motion equations of a two-dimensional wing with cubic stiffness in the pitching direction are established. The aeroelastic system contai...Based on the piston theory of supersonic flow and the energy method, the flutter motion equations of a two-dimensional wing with cubic stiffness in the pitching direction are established. The aeroelastic system contains both structural and aerodynamic nonlinearities. Hopf bifurcation theory is used to analyze the flutter speed of the system. The effects of system parameters on the flutter speed are studied. The 4th order Runge-Kutta method is used to calculate the stable limit cycle responses and chaotic motions of the aeroelastic system. Results show that the number and the stability of equilibrium points of the system vary with the increase of flow speed. Besides the simple limit cycle response of period 1, there are also period-doubling responses and chaotic motions in the flutter system. The route leading to chaos in the aeroelastic model used here is the period-doubling bifurcation. The chaotic motions in the system occur only when the flow speed is higher than the linear divergent speed and the initial condition is very small. Moreover, the flow speed regions in which the system behaves chaos axe very narrow.展开更多
This paper applies washout filter technology to amplitude control of limit cycles emerging from Hopf bifurcation of the van der Pol-Duffing system. The controlling parameters for the appearance of Hopf bifurcation are...This paper applies washout filter technology to amplitude control of limit cycles emerging from Hopf bifurcation of the van der Pol-Duffing system. The controlling parameters for the appearance of Hopf bifurcation are given by the Routh-Hurwitz criteria. Noticeably, numerical simulation indicates that the controllers control the amplitude of limit cycles not only of the weakly nonlinear van der Pol-Duffing system but also of the strongly nonlinear van der Pol-Duffing system. In particular, the emergence of Hopf bifurcation can be controlled by a suitable choice of controlling parameters. Gain-amplitude curves of controlled systems are also drawn.展开更多
In this article,the authors consider a class of Kukles planar polynomial differential system of degree three having an invariant parabola.For this class of second-order differential systems,it is shown that for certai...In this article,the authors consider a class of Kukles planar polynomial differential system of degree three having an invariant parabola.For this class of second-order differential systems,it is shown that for certain values of the parameters the invariant parabola coexists with a center.For other values it can coexist with one,two or three small amplitude limit cycles which are constructed by Hopf bifurcation.This result gives an answer for the question given in[4],about the existence of limit cycles for such class of system.展开更多
To continue the discussion in (Ⅰ ) and ( Ⅱ ),and finish the study of the limit cycle problem for quadratic system ( Ⅲ )m=0 in this paper. Since there is at most one limit cycle that may be created from critic...To continue the discussion in (Ⅰ ) and ( Ⅱ ),and finish the study of the limit cycle problem for quadratic system ( Ⅲ )m=0 in this paper. Since there is at most one limit cycle that may be created from critical point O by Hopf bifurcation,the number of limit cycles depends on the different situations of separatrix cycle to be formed around O. If it is a homoclinic cycle passing through saddle S1 on 1 +ax-y = 0,which has the same stability with the limit cycle created by Hopf bifurcation,then the uniqueness of limit cycles in such cases can be proved. If it is a homoclinic cycle passing through saddle N on x= 0,which has the different stability from the limit cycle created by Hopf bifurcation,then it will be a case of two limit cycles. For the case when the separatrix cycle is a heteroclinic cycle passing through two saddles at infinity,the discussion of the paper shows that the number of limit cycles will change from one to two depending on the different values of parameters of system.展开更多
This paper studies a class of quartic system which is more general and realistic than the quartic accompanying system.x'=-y+ex+lx^2+mxy+ny^2,y'=x(1-Ay)(1+Cy^2),(*)where C 〉 0. Sufficient conditions are ...This paper studies a class of quartic system which is more general and realistic than the quartic accompanying system.x'=-y+ex+lx^2+mxy+ny^2,y'=x(1-Ay)(1+Cy^2),(*)where C 〉 0. Sufficient conditions are obtained for the uniqueness of limit cycle of system (*) and some more in-depth conclusion such as Hopf bifurcation.展开更多
This paper is concerned with a cubic Kolmogorov system with a solution of central quadratic curve which neither contacts with the coordinate axes, nor passes through the origin. The conclusion is that such a system ma...This paper is concerned with a cubic Kolmogorov system with a solution of central quadratic curve which neither contacts with the coordinate axes, nor passes through the origin. The conclusion is that such a system may possess limit cycles.展开更多
It is proved that the quadratic system with a weak focus and a strong focus has at most one limit cycle around the strong focus, and as the weak focus is a 2nd order(or 3rd order) weak focus the quadratic system ha...It is proved that the quadratic system with a weak focus and a strong focus has at most one limit cycle around the strong focus, and as the weak focus is a 2nd order(or 3rd order) weak focus the quadratic system has at most two(one) limit cycles which have (1,1) distribution ((0,1) distribution).展开更多
In this paper, we discuss the limit cycles of the systemdx/dt=y·[1+(A(x)]oy/dt=(-x+δy+α_1x^2+α_2xy+α_5x^2y)[1+B(x)] (1)where A(x)=sum form i=1 to n(a_ix~), B(x)=sum form j=1 to m(β_jx^j) and 1+B(x)>0. We ...In this paper, we discuss the limit cycles of the systemdx/dt=y·[1+(A(x)]oy/dt=(-x+δy+α_1x^2+α_2xy+α_5x^2y)[1+B(x)] (1)where A(x)=sum form i=1 to n(a_ix~), B(x)=sum form j=1 to m(β_jx^j) and 1+B(x)>0. We prove that (1) possesses at most one limit cycle and give out the necessary and sufficient conditions of existence and uniqueness of limit cycles.展开更多
A quantitative analysis of limit cycles and homoclinic orbits, and the bifurcation curve for the Bogdanov-Takens system are discussed. The parameter incremental method for approximate analytical-expressions of these p...A quantitative analysis of limit cycles and homoclinic orbits, and the bifurcation curve for the Bogdanov-Takens system are discussed. The parameter incremental method for approximate analytical-expressions of these problems is given. These analytical-expressions of the limit cycle and homoclinic orbit are shown as the generalized harmonic functions by employing a time transformation. Curves of the parameters and the stability characteristic exponent of the limit cycle versus amplitude are drawn. Some of the limit cycles and homoclinic orbits phase portraits are plotted. The relationship curves of parameters μ and A with amplitude a and the bifurcation diagrams about the parameter are also given. The numerical accuracy of the calculation results is good.展开更多
The dynamics of the confinement transition from L mode to H mode(LH) is investigated in detail theoretically via the extended three-wave coupling model describing the interaction of turbulence and zonal flow(ZF) f...The dynamics of the confinement transition from L mode to H mode(LH) is investigated in detail theoretically via the extended three-wave coupling model describing the interaction of turbulence and zonal flow(ZF) for the first time.Thereinto, turbulence is divided into a positive-frequency(PF) wave and a negative-frequency(NF) one, and the gradient of pressure is added as the auxiliary energy for the system. The LH confinement transition is observed for a sufficiently high input energy. Moreover, it is found that the rotation direction of the limit cycle oscillation(LCO) of PF wave and pressure gradient is reversed during the transition. The mechanism is illustrated by exploring the wave phases. The results presented here provide a new insight into the analysis of the LH transition, which is helpful for the experiments on the fusion devices.展开更多
In this paper we consider a class of polynomial planar system with two small parameters,ε and λ,satisfying 0<ε《λ《1.The corresponding first order Melnikov function M_(1) with respect to ε depends on λ so tha...In this paper we consider a class of polynomial planar system with two small parameters,ε and λ,satisfying 0<ε《λ《1.The corresponding first order Melnikov function M_(1) with respect to ε depends on λ so that it has an expansion of the form M_(1)(h,λ)=∑k=0∞M_(1k)(h)λ^(k).Assume that M_(1k')(h) is the first non-zero coefficient in the expansion.Then by estimating the number of zeros of M_(1k')(h),we give a lower bound of the maximal number of limit cycles emerging from the period annulus of the unperturbed system for 0<ε《λ《1,when k'=0 or 1.In addition,for each k∈N,an upper bound of the maximal number of zeros of M_(1k)(h),taking into account their multiplicities,is presented.展开更多
Bifurcations of subharmonic solutions of order m of a planar periodic perturbed system near a hyperbolic limit cycle are discussed. By using a Poincare map and the method of rescaling a discriminating condition for th...Bifurcations of subharmonic solutions of order m of a planar periodic perturbed system near a hyperbolic limit cycle are discussed. By using a Poincare map and the method of rescaling a discriminating condition for the existence of subharmonic solutions of order m is obtained. An example is given in the end of the paper.展开更多
This paper presents the research on the laws of systematic-parameter dependent variation in the vibration amplitude of drum-brake limit cycle oscillations (LCO). We established a two-degree non-linear dynamic model to...This paper presents the research on the laws of systematic-parameter dependent variation in the vibration amplitude of drum-brake limit cycle oscillations (LCO). We established a two-degree non-linear dynamic model to describe the low-frequency vibration of the drum brake, applied the centre manifold theory to simplify the system, and obtained the LCO amplitude by calculating the normal form of the simplified system at the Hopf bifurcation point. It is indicated that when the friction coefficient is smaller than the friction coefficient at the bifurcation point, the amplitude decreases; whereas with a friction coefficient larger than the friction coefficient of bifurcation point, LCO occurs. The results suggest that it is applicable to suppress the LCO amplitude by changing systematic parameters, and thus improve the safety and ride comfort when applying brake. These findings can be applied to guiding the design of drum brakes.展开更多
This paper is concerned with the number and distributions of limit cycles of a cubic Z2-symmetry Hamiltonian system under quintic perturbation.By using qualitative analysis of differential equation,bifurcation theory ...This paper is concerned with the number and distributions of limit cycles of a cubic Z2-symmetry Hamiltonian system under quintic perturbation.By using qualitative analysis of differential equation,bifurcation theory of dynamical systems and the method of detection function,we obtain that this system exists at least 14 limit cycles with the distribution C91 [C11 + 2(C32 2C12)].展开更多
In this paper, we have proved several theorems which guarantee that the Lienard equation has at least one or n limit cycles without using the traditional assumption G Thus some results in [3-5] are extended. The limit...In this paper, we have proved several theorems which guarantee that the Lienard equation has at least one or n limit cycles without using the traditional assumption G Thus some results in [3-5] are extended. The limit cycles can be located by our theorems. Theorems 3 and 4 give sufficient conditions for the existence of n limit cycles having no need of the conditions that the function F(x) is odd or 'nth order compatible with each other' or 'nth order contained in each other'.展开更多
Influences of prey refuge on the dynamics of a predator-prey model with ratio-dependent functional response are investigated. The local and global stability of positive equilibrium of the system are considered. Theore...Influences of prey refuge on the dynamics of a predator-prey model with ratio-dependent functional response are investigated. The local and global stability of positive equilibrium of the system are considered. Theoretical analysis indicates that constant refuge leads to the system undergo supercritical Hopf bifurcation twice with the birth rate of prey species changing continuously.展开更多
基金supported by the Natural Science Foundation of Ningxia(2022AAC05044)the National Natural Science Foundation of China(12161069)。
文摘This paper deals with the problem of limit cycles for the whirling pendulum equation x=y,y=sin x(cosx-r)under piecewise smooth perturbations of polynomials of cos x,sin x and y of degree n with the switching line x=0.The upper bounds of the number of limit cycles in both the oscillatory and the rotary regions are obtained using the Picard-Fuchs equations,which the generating functions of the associated first order Melnikov functions satisfy.Furthermore,the exact bound of a special case is given using the Chebyshev system.At the end,some numerical simulations are given to illustrate the existence of limit cycles.
文摘The maximal number of limit cycles for a particular type Ⅲ system x = -y + lx2 + mxy, y =x(1 + ax + by) is studied and some errors that appeared in the paper by Suo Mingxia and Yue Xiting (Annals of Differential Equations, 2003,19(3):397-401) are corrected. By translating the system to be considered into the Lienard type and by using some related properties, we obtain several theorems with suitable conditions coefficients of the system, under which we prove that the system has at most two limit cycles. The conclusions improve the results given in Suo and Yue's paper mentioned above.
文摘Aris and Amundson studied a chemical reactor and obtained the two equationsDaoud showed that at most one limit cycle may exist in the region of interest. Itis showed in this paper that other singular points exist and that a stable limitt cycle existsaround the singularity (1/2, 2) when K∈(9-δ, 9).
基金Project supported by the National Natural Science Foundation of China(No.11172199)the KeyProgram of Tianjin Natural Science Foundation of China(No.11JCZDJC25400)
文摘This paper presents a novel mechanical attachment, i.e., nonlinear energy sink (NES), for suppressing the limit cycle oscillation (LCO) of an airfoil. The dynamic responses of a two-degree-of-freedom (2-DOF) airfoil coupled with an NES are studied with the harmonic balance method. Different structure parameters of the NES, i.e., mass ratio between the NES and airfoil, NES offset, NES damping, and nonlinear stiffness in the NES, are chosen for studying the effect of the LCO suppression on an aeroelastic system with a supercritical Hopf bifurcation or subcritical Hopf bifurcation, respectively. The results show that the structural parameters of the NES have different influence on the supercritical Hopf bifurcation system and the subcritical Hopf bifurcation system.
基金supported by the National Natural Science Foundation of China and China Academy of Engineering Physics(No. 10576024).
文摘Based on the piston theory of supersonic flow and the energy method, the flutter motion equations of a two-dimensional wing with cubic stiffness in the pitching direction are established. The aeroelastic system contains both structural and aerodynamic nonlinearities. Hopf bifurcation theory is used to analyze the flutter speed of the system. The effects of system parameters on the flutter speed are studied. The 4th order Runge-Kutta method is used to calculate the stable limit cycle responses and chaotic motions of the aeroelastic system. Results show that the number and the stability of equilibrium points of the system vary with the increase of flow speed. Besides the simple limit cycle response of period 1, there are also period-doubling responses and chaotic motions in the flutter system. The route leading to chaos in the aeroelastic model used here is the period-doubling bifurcation. The chaotic motions in the system occur only when the flow speed is higher than the linear divergent speed and the initial condition is very small. Moreover, the flow speed regions in which the system behaves chaos axe very narrow.
基金Project supported by the National Natural Science Foundation of China (Grant No 10672053)
文摘This paper applies washout filter technology to amplitude control of limit cycles emerging from Hopf bifurcation of the van der Pol-Duffing system. The controlling parameters for the appearance of Hopf bifurcation are given by the Routh-Hurwitz criteria. Noticeably, numerical simulation indicates that the controllers control the amplitude of limit cycles not only of the weakly nonlinear van der Pol-Duffing system but also of the strongly nonlinear van der Pol-Duffing system. In particular, the emergence of Hopf bifurcation can be controlled by a suitable choice of controlling parameters. Gain-amplitude curves of controlled systems are also drawn.
基金NNSF of China(10671211)NSF of Hunan Province(07JJ3005)USM(120628 and 120627)
文摘In this article,the authors consider a class of Kukles planar polynomial differential system of degree three having an invariant parabola.For this class of second-order differential systems,it is shown that for certain values of the parameters the invariant parabola coexists with a center.For other values it can coexist with one,two or three small amplitude limit cycles which are constructed by Hopf bifurcation.This result gives an answer for the question given in[4],about the existence of limit cycles for such class of system.
基金Project supported by the National Natural Science Foundation of China (10471066).
文摘To continue the discussion in (Ⅰ ) and ( Ⅱ ),and finish the study of the limit cycle problem for quadratic system ( Ⅲ )m=0 in this paper. Since there is at most one limit cycle that may be created from critical point O by Hopf bifurcation,the number of limit cycles depends on the different situations of separatrix cycle to be formed around O. If it is a homoclinic cycle passing through saddle S1 on 1 +ax-y = 0,which has the same stability with the limit cycle created by Hopf bifurcation,then the uniqueness of limit cycles in such cases can be proved. If it is a homoclinic cycle passing through saddle N on x= 0,which has the different stability from the limit cycle created by Hopf bifurcation,then it will be a case of two limit cycles. For the case when the separatrix cycle is a heteroclinic cycle passing through two saddles at infinity,the discussion of the paper shows that the number of limit cycles will change from one to two depending on the different values of parameters of system.
基金Supported by the Natural Science Foundation of Fujian Province(Z0511052,2006J0209)the Foundation of Fujian Education Department(JA04158,JA04274)and the Foundation of Developing ScienceTechnology of Fuzhou University(2005-QX-20)
文摘This paper studies a class of quartic system which is more general and realistic than the quartic accompanying system.x'=-y+ex+lx^2+mxy+ny^2,y'=x(1-Ay)(1+Cy^2),(*)where C 〉 0. Sufficient conditions are obtained for the uniqueness of limit cycle of system (*) and some more in-depth conclusion such as Hopf bifurcation.
基金The NSF of Liaoning provinceFoundation of returned doctors and Foundation of LiaoningEducation Committee.
文摘This paper is concerned with a cubic Kolmogorov system with a solution of central quadratic curve which neither contacts with the coordinate axes, nor passes through the origin. The conclusion is that such a system may possess limit cycles.
文摘It is proved that the quadratic system with a weak focus and a strong focus has at most one limit cycle around the strong focus, and as the weak focus is a 2nd order(or 3rd order) weak focus the quadratic system has at most two(one) limit cycles which have (1,1) distribution ((0,1) distribution).
文摘In this paper, we discuss the limit cycles of the systemdx/dt=y·[1+(A(x)]oy/dt=(-x+δy+α_1x^2+α_2xy+α_5x^2y)[1+B(x)] (1)where A(x)=sum form i=1 to n(a_ix~), B(x)=sum form j=1 to m(β_jx^j) and 1+B(x)>0. We prove that (1) possesses at most one limit cycle and give out the necessary and sufficient conditions of existence and uniqueness of limit cycles.
基金the National Natural Science Foundation of China (No.10672193)
文摘A quantitative analysis of limit cycles and homoclinic orbits, and the bifurcation curve for the Bogdanov-Takens system are discussed. The parameter incremental method for approximate analytical-expressions of these problems is given. These analytical-expressions of the limit cycle and homoclinic orbit are shown as the generalized harmonic functions by employing a time transformation. Curves of the parameters and the stability characteristic exponent of the limit cycle versus amplitude are drawn. Some of the limit cycles and homoclinic orbits phase portraits are plotted. The relationship curves of parameters μ and A with amplitude a and the bifurcation diagrams about the parameter are also given. The numerical accuracy of the calculation results is good.
基金supported by the National Natural Science Foundation of China(Grant Nos.11305010 and 11475026)the Joint Foundation of the National Natural Science FoundationChina Academy of Engineering Physics(Grant No.U1530153)
文摘The dynamics of the confinement transition from L mode to H mode(LH) is investigated in detail theoretically via the extended three-wave coupling model describing the interaction of turbulence and zonal flow(ZF) for the first time.Thereinto, turbulence is divided into a positive-frequency(PF) wave and a negative-frequency(NF) one, and the gradient of pressure is added as the auxiliary energy for the system. The LH confinement transition is observed for a sufficiently high input energy. Moreover, it is found that the rotation direction of the limit cycle oscillation(LCO) of PF wave and pressure gradient is reversed during the transition. The mechanism is illustrated by exploring the wave phases. The results presented here provide a new insight into the analysis of the LH transition, which is helpful for the experiments on the fusion devices.
基金The first author is supported by the National Natural Science Foundation of China(11671013)the second author is supported by the National Natural Science Foundation of China(11771296).
文摘In this paper we consider a class of polynomial planar system with two small parameters,ε and λ,satisfying 0<ε《λ《1.The corresponding first order Melnikov function M_(1) with respect to ε depends on λ so that it has an expansion of the form M_(1)(h,λ)=∑k=0∞M_(1k)(h)λ^(k).Assume that M_(1k')(h) is the first non-zero coefficient in the expansion.Then by estimating the number of zeros of M_(1k')(h),we give a lower bound of the maximal number of limit cycles emerging from the period annulus of the unperturbed system for 0<ε《λ《1,when k'=0 or 1.In addition,for each k∈N,an upper bound of the maximal number of zeros of M_(1k)(h),taking into account their multiplicities,is presented.
文摘Bifurcations of subharmonic solutions of order m of a planar periodic perturbed system near a hyperbolic limit cycle are discussed. By using a Poincare map and the method of rescaling a discriminating condition for the existence of subharmonic solutions of order m is obtained. An example is given in the end of the paper.
基金the Natural Science Foundation of China (No. 50075029)
文摘This paper presents the research on the laws of systematic-parameter dependent variation in the vibration amplitude of drum-brake limit cycle oscillations (LCO). We established a two-degree non-linear dynamic model to describe the low-frequency vibration of the drum brake, applied the centre manifold theory to simplify the system, and obtained the LCO amplitude by calculating the normal form of the simplified system at the Hopf bifurcation point. It is indicated that when the friction coefficient is smaller than the friction coefficient at the bifurcation point, the amplitude decreases; whereas with a friction coefficient larger than the friction coefficient of bifurcation point, LCO occurs. The results suggest that it is applicable to suppress the LCO amplitude by changing systematic parameters, and thus improve the safety and ride comfort when applying brake. These findings can be applied to guiding the design of drum brakes.
基金Supported by the Natural Science Foundation of China(10802043 10826092) Acknowledgements We are grateful to Prof Li Ji-bin for his kind help and the referees' valuable suggestions.
文摘This paper is concerned with the number and distributions of limit cycles of a cubic Z2-symmetry Hamiltonian system under quintic perturbation.By using qualitative analysis of differential equation,bifurcation theory of dynamical systems and the method of detection function,we obtain that this system exists at least 14 limit cycles with the distribution C91 [C11 + 2(C32 2C12)].
文摘In this paper, we have proved several theorems which guarantee that the Lienard equation has at least one or n limit cycles without using the traditional assumption G Thus some results in [3-5] are extended. The limit cycles can be located by our theorems. Theorems 3 and 4 give sufficient conditions for the existence of n limit cycles having no need of the conditions that the function F(x) is odd or 'nth order compatible with each other' or 'nth order contained in each other'.
基金Supported by the NNSF of China(11126284)Supported by the NSF of Department of Education of Henan Province(12A110012)Supported by the Young Scientific Research Foundation of Henan Normal University(1001)
文摘Influences of prey refuge on the dynamics of a predator-prey model with ratio-dependent functional response are investigated. The local and global stability of positive equilibrium of the system are considered. Theoretical analysis indicates that constant refuge leads to the system undergo supercritical Hopf bifurcation twice with the birth rate of prey species changing continuously.