About this Journal Submit a Manuscript Table of Contents
Journal of Applied Mathematics

Volume 2014 (2014), Article ID 701545, 11 pages

http://dx.doi.org/10.1155/2014/701545
Research Article

Analysis of a Single Species Model with Dissymmetric Bidirectional Impulsive Diffusion and Dispersal Delay

College of Mathematics and System Sciences, Xinjiang University, Urumqi 830046, China

Received 19 November 2013; Accepted 25 February 2014; Published 9 April 2014

Academic Editor: Yongkun Li

Copyright © 2014 Haiyun Wan 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

In most models of population dynamics, diffusion between two patches is assumed to be either continuous or discrete, but in the real natural ecosystem, impulsive diffusion provides a more suitable manner to model the actual dispersal (or migration) behavior for many ecological species. In addition, the species not only requires some time to disperse or migrate among the patches but also has some possibility of loss during dispersal. In view of these facts, a single species model with dissymmetric bidirectional impulsive diffusion and dispersal delay is formulated. Criteria on the permanence and extinction of species are established. Furthermore, the realistic conditions for the existence, uniqueness, and the global stability of the positive periodic solution are obtained. Finally, numerical simulations and discussion are presented to illustrate our theoretical results.

1. Introduction

In the last few years, mathematicians and ecologists have been actively investigating the dispersal of populations, a ubiquitous phenomenon in population dynamics. Levin [1] showed that both spatial dispersal of populations and population dynamics are much affected by spatial heterogeneity. In real life, dispersal often occurs among patches in ecological environments; because of the ecological effects of human activities and industries, such as the location of manufacturing industries and the pollution of the atmosphere, soil, and rivers, reproduction- and population-based territories and other habitats have been broken into patches. Thus, realistic models should include dispersal processes that take into consideration the effects of spatial heterogeneity.

In recent years, increasing attention has been paid to the dynamics of a large number of mathematical models with diffusion, and many nice results have been obtained. The persistence and extinction for ordinary differential equation and delayed differential equation models were investigated in [26]. Global stability of equilibrium and periodic solution for diffusing model were studied in [712]. However, in all of above population dispersing systems, it is always assumed that the dispersal occurs at every time. For example, in [7], Beretta and Takeuchi proposed the following single-species diffusion Volterra models with continuous time delays: where is the number of patches, and is the population density in the th patch. The form of the dispersal established in this model is continuous; that is, the dispersal is always happening at any time.

Actually, real dispersal behavior is very complicated and is always influenced by environmental change and human activities. In many practical situations, it is often the case that maybe one of the species suffers a significant loss or increase in density for some reason at some transitory time slots. These short-term perturbations are often assumed to be in the form of impulses in the modeling process. For example, when winter comes, birds will migrate between patches in search for a better environment, whereas they do not diffuse in other seasons, and the excursion of foliage seeds occurs during a fixed period of time every year. Therefore, impulsive differential equations [13] provide a natural description of such system. With the developments and applications of impulsive differential equations, theories of impulsive differential equations have been introduced into population dynamics, and many important studies have been performed [1420].

In [14], the authors studied the following autonomous single-species model with impulsively bidirectional diffusion: where , are the intrinsic growth rate and density-dependent parameters of the population and is the dispersal rate in the th patch. Consider , where represents the density of population in the th patch immediately after the th diffusion pulse at time , while represents the density of population in the th patch before the th diffusion pulse at time ( the period of dispersal between any two pulse events is a positive constant, ). It is assumed here that the net exchange from the th patch to th patch is proportional to the difference of population densities. The dispersal behavior of populations between two patches occurs only at the impulsive instants . Obviously, in this model, species inhabits, respectively, two patches before the pulse appears, and when the time at the pulse comes, species in two patches disperses from one patch to the other. The boundedness and global stability of positive periodic solution were obtained.

Time delay often appears in many control systems (such as aircraft, chemical, or process control systems) either in the state, the control input, or the measurements. In order to reflect the dynamical behaviors of models that depend on the past history of system, it is often necessary to incorporate time delays into systems [21]. There have been extensive theoretical works on delay differential equations in the past three decades. The research topics include global asymptotic stability of equilibria, existence of periodic solutions, complicated behavior, and chaos (e.g., [8, 2224]).

Takeuchi et al. in [25] studied the following population model with time delays that introduced the dispersal time for individuals to move from one patch to other patches: where , , are positive constants, , , , are nonnegative constants, satisfy , , , are constants, and some of them may be negative. In this paper, the authors took account of dispersal delay; however, they assumed that the dispersal is continuous.

It is well known that the application of impulsive delay differential equations to population dynamics has been an interesting topic since it is reasonable and correct in modelling the evolution of population, such as pest management [26].

However, in all of the impulsive dispersal models studied up till now, there are few papers considering the dispersal delay, which is really a pity. Actually, in the real world, the migration between patches is usually not immediate; that is, dispersal processes often involve time delay. For example, elks move from higher to lower elevations to escape cold in winter, and ungulates migrate annually among grazing areas to follow spatiotemporal changes in rainfall. Obviously, this kind of dispersal delay between patches extensively exists in the real world. Therefore, it is a very basilic problem to research this kind of population dynamic systems.

Moreover, in the above impulsive dispersal models, it is assumed that the dispersal occurs between homogeneous habitat patches; that is, the dispersal rate between any two patches is equal or symmetrical [1, 11], which is really too idealized for a real ecosystem. Actually, in the real world, due to the heterogeneity of the spatiotemporal distributions in nature, movement between fragments of patches is usually not the same rate in both directions. In addition, once the individuals leave their present habitat, they may not successfully reach a new one, due to predation, harvesting, or other reasons, so that there are traveling losses. Thus, the dispersal rates among these patches are not always the same. Rather, in real ecological situations they are different (or dissymmetrical) [27, 28].

Therefore, it is our basilic goal to investigate a single species model with dissymmetric impulse dispersal and dispersal delay. Motivated by the calculation hereinbefore, in this paper, we extend system (2) with dispersal delay and dispersal loss and consider it where is the rate of population emigrating from the th patch and is the rate of population immigrating from the th patch. Here we assume , which means that there possibly exists mortality during migration between two patches. stands for the time delay; that is, a period of time of species dispersing between patches.

The organization of this paper is as follows. In Section 2, as preliminaries, the definition of permanence and some useful lemmas are introduced. From discrete dynamic system theory, we establish the stroboscopic map of system (4), by which we can obtain the dynamical behaviors of it. In Section 3, the results of permanence and extinction for the system are presented. The existence and the uniqueness of the positive periodic solution for system (4) are established in Section 4. In Section 5, using the discrete dynamic system theory in [29], we can get the global stability of the positive periodic solution for the system. Finally, we give a brief discussion and our theoretical results are conformed by numerical simulations.

2. Preliminaries

In this section we introduce a definition and some notations and state some results which will be useful in subsequent sections.

For any fixed , let be a piecewise continuous function such that is continuous in , and exist. For any fixed , let denote the Banach space of all such piecewise continuous functions with the norm . Further, let = for all and for .

Motivated by the biological background of system (4), in this paper we always assume that all solutions of system (4) satisfy the following initial conditions: where . For any , by the fundamental theory of impulsive functional differential equations [30, 31], system (4) has a unique solution satisfying the initial conditions (5).

Definition 1. System (4) is said to be permanent, if there are positive constants and such that for any positive solution of system (4).

We consider the following scalar impulsive differential equation: where , , is an impulsive time sequence, , is continuous, and is a nondecreasing function. We have the following comparison theorem [30, 31] for (7).

Lemma 2. Let be a solution of system (7) defined on and satisfy

If , then for all .

Lemma 3 (see [32]). Consider the following nonlinear impulsive system: where and are bounded and continuous -periodic functions defined on , for all and impulsive coefficients , is a fixed positive integer, and . If , then there exists a unique positive periodic solution of system (9), which is globally asymptotically stable.

Remark 4. If system (9) degenerates into the following autonomous impulsive differential equation: where , , are positive constants, . As a consequence of Lemma 3, we have the following result: if , then system (10) has a unique positive periodic solution , which is globally asymptotically stable. In fact, here .

Lemma 5 (see [29]). Let be continuous, in int , and suppose exists with . In addition, assume(a) , if ;(b) , if .

If , let . If , then, for every as ; if , then either as for every or there exists a unique nonzero fixed point of . In the latter case, and, for every , as .

If , then either as for every or there exists a unique fixed point of . In the latter case, and, for every , as .

Next, we analyze system (4). Integrating and solving the first two equations of system (4) between pulses, we have

Similarly, considering the last two equations of system (4), we obtain the following stroboscopic map: here , , , .

Remark 6. System (12) is a difference system, which means that densities of population in two patches have values at the previous pulse. We are, in other words, stroboscopically sampling at its pulsing period. The dynamical behavior of system (12), coupled with (11), determines the dynamical behavior of system (4). In the following sections, we will focus our attention on system (12) and investigate various aspects of its dynamical behavior.

To write system (12) as a map, we define the map :

The set of all iterations of the map is equivalent to the set of all density sequences generated by system (12); is the map evaluated at the point . Consequently, in system (12), describes the population densities in the time .

On the positivity of solutions of system (4) we have the following result.

Lemma 7. The solution of system (4) with initial condition (5) is positive, that is, on the interval of the existence.

The proof of Lemma 7 is simple; we hence omit it here.

3. Permanence and Extinction

In this section, we present conditions to ensure that system (12) is permanent and extinct which will imply the permanence and extinction of system (4). The permanence plays an important role in mathematical ecology since the criterion of permanence for ecological systems is a condition ensuring the long-term survival of all species. So, we firstly prove system (12) is permanent.

Theorem 8. Suppose hold; then system (12) is permanent.

Proof. Let be the solution of system (4) satisfying the initial conditions (5). From the first equation of system (12), we have

Similarly, we have

Hence, by (14) and (15) we know that system (12) has an ultimately upper bound.

Next, we prove that all the solutions of system (4) are ultimately below bounded. Since , from the third equation of system (4), we have

Similarly,

Thus, system (4) becomes

From (18), we find that there is no relation between and . Therefore, we will discuss them, respectively;

If holds, from Remark 4, we can obtain that the auxiliary system has a unique positive periodic solution which is globally asymptotically stable.

Let be the solution of system (20) with initial value . By Lemma 2, we have

Hence, for any sufficiently small, there exists a such that

Similarly, if holds, for above , there exists a such that

Denote and ; then we have and , . Finally, we can determine that there exist constants , , such that , and . The proof of Theorem 8 is completed.

Next, we present condition to ensure that system (12) is extinct.

Theorem 9. System (12) is extinct if

Proof. Let us consider the system (13). Obviously, is continuous in , and . We obtain

Obviously, ; if ; if , . We have the characteristic equation of :

Let ; then we have

Assume ; then by (27) we can obtain that is, which contradicts with (24). Therefore we have . By Lemma 5, we can get as , which means that system (12) is extinct. This completes the proof.

4. Existence and Uniqueness of Positive Periodic Solution

In this part, we will prove the existence and uniqueness of the fixed points of system (12), which means that system (4) has a uniquely positive periodic solution.

Theorem 10. If holds, then there exists a unique positive fixed point of system (12).

Proof. Corresponding to (12), let us consider the following system:

From (30), we have hence

From (30), we also obtain

Thus

Let ; then . And is an increasing function on the interval . Since , so there exists such that . We can easily find that . By the zero theory of continuous function, there exists such that

Next, we will prove the uniqueness of the fixed point.

It follows from (33) that we obtain

Let then

By (31), we have ; since , so . Similarly, we have . Therefore, we obtain , which implies that is a decreasing function on the interval .

Since using the zero theory of continuous function, there exists a unique point such that . Besides, thus which, together with , leads to . By , we have that there exists a unique point such that . The proof is completed.

5. Global Stability

Now, we prove that the positive fixed points of system (30) are globally stable by using Lemma 5, which means that the positive periodic solution of system (4) is globally stable.

Theorem 11. If holds, then there exists a unique positive fixed point of the map , and, for every as .

Proof. For any small , , we make the change of variable

By (13), we get the map ; that is,

Now, we show that satisfies the hypotheses of Lemma 5. It is easy to see that is continuous, in , and .

Since so, . Obviously, if , then ; if , then . It satisfies all the conditions of Lemma 5; then, for every , we have as . Corresponding to coordinate, this means, for , , the system (30) tends to the unique fixed point.

It follows from the permanence of system (12) that we have , for any initial value .

From the above analysis, we can know that, for every , , the trajectory of system (12) will tend to . The proof is completed.

6. Numerical Simulation and Discussion

In order to test the validity of our results, first, for (4) we use the parameters values (Val. 1) in Table 1. We can easily test that the assumptions in Theorems 8 and 10 hold, which means the populations in the two patches are permanent and have a unique periodic solution which is globally stable (see Figure 1(a)). Moreover, if, in Table 1, we consider the influence of time delay, then we can see that the permanence and stability for species unchanged. The details are given in Table 2. However, the longer the duration of the time delay, the lower the limit inferior and the limit superior of (see Figures 1(b), 1(c), and 1(d)). This implies that the case with dispersal delay is harmful to live for species .

tab1
Table 1: Parameter values used in the simulations of Model (4).
tab2
Table 2: Simulations of Model (4).
fig1
Figure 1: (a), (b), (c), and (d) Dynamical behavior of system (4). Here, we take three sets of initial values (1, 1.2), (0.6, 0.8), and (0.2, 1).

Next, we take the parameters values (Val. 2) in Table 1. We can easily test that assumption in Theorem 9 holds, which implies the populations in the two patches are extinct (see Figure 2). Comparing Figure 2(a) and Figure 2(b), we realize that system (4) with time delay accelerates the extinction comparing with no delay (see Figure 2). This is reasonable from a biological point of view. Without delay means less loss during dispersion, which implies that more members can arrive to other patches. Otherwise, populations go extinct due to much loss. The details are given in Table 3.

tab3
Table 3: Simulations of Model (4).
fig2
Figure 2: (a) and (b) Dynamical behavior of system (4). Here, we take three sets of initial values (1.7, 2.1), (1.9, 1.5), and (1.4, 2.2).

Furthermore, we take , , , , and keep other parameters unchanged from in Val. 1. Here we can see that the period of individual intrinsic growth rate and density dependence rate is equal to the period of migration . In this case, we have let the period of the environment match the period of migration. If we take , from simulation (see Figure 3(a)), we can see the populations are permanent and have a uniquely periodic solution which is globally stable. However, if we take the migration period , respectively, with other parameters unchanged, from numerical simulations (see Figures 3(b) and 3(c)), we can see that population dynamics change from almost periodic to chaotic. For Case 4 of Table 4, we consider the influence of time delay (see Figure 3(d)). Comparing Figures 3(c) and 3(d), we realize that system (4) with time delay is more complicated than without. The details are given in Table 4.

tab4
Table 4: Simulations of Model (4).
fig3
Figure 3: (a), (b), (c), and (d) Dynamical behavior of system (4). Here, we take the initial value = (1.735, 1.873).

Lastly, if we take and keep other parameters unchanged with Figure 3, by numerical simulations (see Figure 4), we find that all of the solutions of system (4) which through the initial points will converge to the positive periodic solution . Therefore, we can guess that under the assumptions of Theorem 10 system (4) has a unique positive periodic solution which is globally asymptotically stable. In addition, the periodic solution with time delay is larger than without delay which indicates that the duration of the time delay is beneficial to live for species compared with species .

fig4
Figure 4: (a) and (b) Dynamical behavior of system (4). We take a series of initial points, such as (1.935,2.3), (1.939,2.34), and (1.943,2.38).

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (nos. 10901130, 11361059, 10961022, and 11271312), the China Scholarship Council, the Natural Science Foundation of Xinjiang Province of China (2012211B07), and the Scientific Research Programmes of Colleges in Xinjiang (XJEDU2013I03).

References

  1. S. A. Levin, “Dispersion and population interactions,” The American Naturalist, vol. 108, pp. 207–228, 1974.
  2. Z. Teng and Z. Lu, “The effect of dispersal on single-species nonautonomous dispersal models with delays,” Journal of Mathematical Biology, vol. 42, no. 5, pp. 439–454, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  3. J. Cui, Y. Takeuchi, and Z. Lin, “Permanence and extinction for dispersal population systems,” Journal of Mathematical Analysis and Applications, vol. 298, no. 1, pp. 73–93, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  4. J. Cui and L. Chen, “Permanence and extinction in logistic and Lotka-Volterra systems with diffusion,” Journal of Mathematical Analysis and Applications, vol. 258, no. 2, pp. 512–535, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  5. R. Xu and Z. Ma, “The effect of dispersal on the permanence of a predator-prey system with time delay,” Nonlinear Analysis. Real World Applications, vol. 9, no. 2, pp. 354–369, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  6. Y. Takeuchi, J. Cui, R. Miyazaki, and Y. Saito, “Permanence of delayed population model with dispersal loss,” Mathematical Biosciences, vol. 201, no. 1-2, pp. 143–156, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  7. E. Beretta and Y. Takeuchi, “Global stability of single-species diffusion Volterra models with continuous time delays,” Bulletin of Mathematical Biology, vol. 49, no. 4, pp. 431–448, 1987. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  8. E. Beretta and Y. Takeuchi, “Global asymptotic stability of Lotka-Volterra diffusion models with continuous time delay,” SIAM Journal on Applied Mathematics, vol. 48, no. 3, pp. 627–651, 1988. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  9. E. Beretta, P. Fergola, and C. Tenneriello, “Ultimate boundedness for nonautonomous diffusive Lotka-Volterra patches,” Mathematical Biosciences, vol. 92, no. 1, pp. 29–53, 1988. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  10. H. I. Freedman, J. B. Shukla, and Y. Takeuchi, “Population diffusion in a two-patch environment,” Mathematical Biosciences, vol. 95, no. 1, pp. 111–123, 1989. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  11. A. Hastings, “Dynamics of a single species in a spatially varying environment: the stabilizing role of high dispersal rates,” Journal of Mathematical Biology, vol. 16, no. 1, pp. 49–55, 1982/83. View at Publisher · View at Google Scholar · View at MathSciNet
  12. W. Wang, L. Chen, and Z. Lu, “Global stability of a population dispersal in a two-patch environment,” Dynamic Systems and Applications, vol. 6, no. 2, pp. 207–215, 1997. View at Zentralblatt MATH · View at MathSciNet
  13. L. Dong, L. Chen, and P. Shi, “Periodic solutions for a two-species nonautonomous competition system with diffusion and impulses,” Chaos, Solitons & Fractals, vol. 32, no. 5, pp. 1916–1926, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  14. J. Hui and L.-S. Chen, “A single species model with impulsive diffusion,” Acta Mathematicae Applicatae Sinica, vol. 21, no. 1, pp. 43–48, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  15. L. Wang, Z. Liu, Jinghui, and L. Chen, “Impulsive diffusion in single species model,” Chaos, Solitons & Fractals, vol. 33, no. 4, pp. 1213–1219, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  16. G. Ballinger and X. Liu, “Permanence of population growth models with impulsive effects,” Mathematical and Computer Modelling, vol. 26, no. 12, pp. 59–72, 1997. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  17. L. Zhang and Z. Teng, “N-species non-autonomous Lotka-Volterra competitive systems with delays and impulsive perturbations,” Nonlinear Analysis: Real World Applications, vol. 12, no. 6, pp. 3152–3169, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  18. J. Vandermeer, L. Stone, and B. Blasius, “Categories of chaos and fractal basin boundaries in forced predator-prey models,” Chaos, Solitons & Fractals, vol. 12, no. 2, pp. 265–276, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  19. L. Zhang, Z. Teng, D. L. DeAngelis, and S. Ruan, “Single species models with logistic growth and dissymmetric impulse dispersal,” Mathematical Biosciences, vol. 241, no. 2, pp. 188–197, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  20. Z. Zhao, X. Zhang, and L. Chen, “The effect of pulsed harvesting policy on the inshore-offshore fishery model with the impulsive diffusion,” Nonlinear Dynamics, vol. 63, no. 4, pp. 537–545, 2011. View at Publisher · View at Google Scholar · View at MathSciNet
  21. X.-P. Yan and W.-T. Li, “Hopf bifurcation and global periodic solutions in a delayed predator-prey system,” Applied Mathematics and Computation, vol. 177, no. 1, pp. 427–445, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  22. L. Zhang and Z. Teng, “Permanence for a delayed periodic predator-prey model with prey dispersal in multi-patches and predator density-independent,” Journal of Mathematical Analysis and Applications, vol. 338, no. 1, pp. 175–193, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  23. H. L. Smith, Monotone Dynamical Systems: An Introduction to the Theory of Competitive and Cooperative Systems, vol. 41 of Mathematical Surveys and Monographs, American Mathematical Society, Providence, RI, USA, 1995. View at MathSciNet
  24. X.-Q. Zhao, “Global attractivity in a class of nonmonotone reaction-diffusion equations with time delay,” Canadian Applied Mathematics Quarterly, vol. 17, no. 1, pp. 271–281, 2009. View at Zentralblatt MATH · View at MathSciNet
  25. Y. Takeuchi, W. Wang, and Y. Saito, “Global stability of population models with patch structure,” Nonlinear Analysis: Real World Applications, vol. 7, no. 2, pp. 235–247, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  26. X. Meng and L. Chen, “Permanence and global stability in an impulsive Lotka-Volterra n-species competitive system with both discrete delays and continuous delays,” International Journal of Biomathematics, vol. 1, no. 2, pp. 179–196, 2008. View at Publisher · View at Google Scholar · View at MathSciNet
  27. H. I. Freedman, B. Rai, and P. Waltman, “Mathematical models of population interactions with dispersal. II. Differential survival in a change of habitat,” Journal of Mathematical Analysis and Applications, vol. 115, no. 1, pp. 140–154, 1986. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  28. R. H. MacArthur and E. O. Wilson, The Theory of Island Biogeography, Princeton University Press, 1967.
  29. H. L. Smith, “Cooperative systems of differential equations with concave nonlinearities,” Nonlinear Analysis: Theory, Methods & Applications, vol. 10, no. 10, pp. 1037–1052, 1986. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  30. V. Lakshmikantham, D. D. Baĭnov, and P. S. Simeonov, Theory of Impulsive Differential Equations, vol. 6 of Series in Modern Applied Mathematics, World Scientific, Teaneck, NJ, USA, 1989. View at MathSciNet
  31. D. Bainov and P. Simeonov, Impulsive Differential Equations: Periodic Solutions and Applications, Longman, London, UK, 2003.
  32. F. Zhang, S. Gao, and Y. Zhang, “Effects of pulse culling on population growth of migratory birds and economical birds,” Nonlinear Dynamics, vol. 67, no. 1, pp. 767–779, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet