Research article Special Issues

Global analysis and optimal harvesting for a hybrid stochastic phytoplankton-zooplankton-fish model with distributed delays

  • Received: 15 June 2020 Accepted: 24 August 2020 Published: 14 September 2020
  • In this paper, we formulate a phytoplankton-zooplankton-fish model with distributed delays and hybrid stochastic noises involving Brownian motion and Markov chain, and propose an optimal harvesting problem pursuing the maximum of total economic income. By global analysis in terms of some system parameters, we investigate the dynamical behaviors on the well-posedness, bounded- ness, persistence, extinction, stability and attractiveness of the solutions for the stochastic delayed system. Moreover, we provide sufficient and necessary condition ensuring the existence of the optimization solution for the optimization problem and obtain the optimal harvesting effect and the maximum of sustainable yield. Lastly, two numerical examples and their simulations are given to illustrate the effectiveness of our results.

    Citation: Yuanpei Xia, Weisong Zhou, Zhichun Yang. Global analysis and optimal harvesting for a hybrid stochastic phytoplankton-zooplankton-fish model with distributed delays[J]. Mathematical Biosciences and Engineering, 2020, 17(5): 6149-6180. doi: 10.3934/mbe.2020326

    Related Papers:

    [1] Ruhma Mahmood, Mahmood Shaukat, Mahmood S Choudhery . Biological properties of mesenchymal stem cells derived from adipose tissue, umbilical cord tissue and bone marrow. AIMS Cell and Tissue Engineering, 2018, 2(2): 78-90. doi: 10.3934/celltissue.2018.2.78
    [2] Ryota Tamura, Masahiro Toda . Stem cell-based therapies for neurological disorders. AIMS Cell and Tissue Engineering, 2018, 2(1): 24-46. doi: 10.3934/celltissue.2018.1.24
    [3] Robert J. Henning, Qing Zhu, Xiao Wang . Therapeutic effects of paracrine factors secreted by human umbilical cord blood mononuclear cells in myocardial infarctions
    Paracrine effects of cord blood cells. AIMS Cell and Tissue Engineering, 2018, 2(4): 220-237. doi: 10.3934/celltissue.2018.4.220
    [4] Edgar J. Cubillo, Sang M. Ngo, Alejandra Juarez, Joshuah Gagan, Gisel D. Lopez, David A. Stout . Embryonic stem cell therapy applications for autoimmune, cardiovascular, and neurological diseases: A review. AIMS Cell and Tissue Engineering, 2017, 1(3): 191-223. doi: 10.3934/celltissue.2017.3.191
    [5] Lauren N. Randolph, Yuqian Jiang, Xiaojun Lian . Stem Cell Engineering and Differentiation for Disease Modeling and Cell-based Therapies. AIMS Cell and Tissue Engineering, 2017, 1(2): 140-157. doi: 10.3934/celltissue.2017.2.140
    [6] Shihori Tanabe, Ryuichi Ono . The gene and microRNA networks of stem cells and reprogramming. AIMS Cell and Tissue Engineering, 2018, 2(4): 238-245. doi: 10.3934/celltissue.2018.4.238
    [7] Shihori Tanabe, Kazuhiko Aoyagi, Hiroshi Yokozaki, Hiroki Sasaki . Molecular pathway network of EFNA1 in cancer and mesenchymal stem cells. AIMS Cell and Tissue Engineering, 2018, 2(2): 58-77. doi: 10.3934/celltissue.2018.2.58
    [8] Kiran Dhaliwal, Muhammad Junejo . Penile Tissue Engineering: a Review of the Current Progress of Penile Reconstruction. AIMS Cell and Tissue Engineering, 2017, 1(3): 180-190. doi: 10.3934/celltissue.2017.3.180
    [9] Vijay Kumar Kuna, Sanchari Paul, Bo Xu, Robert Sjöback, Suchitra Sumitran-Holgersson . Human fetal kidney cells regenerate acellular porcine kidneys via upregulation of key transcription factors involved in kidney development
    . AIMS Cell and Tissue Engineering, 2019, 3(1): 26-46. doi: 10.3934/celltissue.2019.1.26
    [10] Cindy E. McKinney, Katherine M. Baumgarner . CRISPR engineering cardiometabolic disease models using human iPSC. AIMS Cell and Tissue Engineering, 2018, 2(3): 185-202. doi: 10.3934/celltissue.2018.3.185
  • In this paper, we formulate a phytoplankton-zooplankton-fish model with distributed delays and hybrid stochastic noises involving Brownian motion and Markov chain, and propose an optimal harvesting problem pursuing the maximum of total economic income. By global analysis in terms of some system parameters, we investigate the dynamical behaviors on the well-posedness, bounded- ness, persistence, extinction, stability and attractiveness of the solutions for the stochastic delayed system. Moreover, we provide sufficient and necessary condition ensuring the existence of the optimization solution for the optimization problem and obtain the optimal harvesting effect and the maximum of sustainable yield. Lastly, two numerical examples and their simulations are given to illustrate the effectiveness of our results.


    Can cats or dogs barter? Trade is one of the essential features of human intelligence. Brosnan et al., (2008) reported that chimpanzees, who are an intelligent species like humans, can trade but are reluctant to trade. This subsequently led to major divergence in the fates of the two species. The market is nothing but an expression of human intelligence. As Malone and Bernstein (2015) mentioned, intelligence does not arise only in individual brains; it also arises in groups of individuals. The securities market, including crypto-assets, is the ultimate expression of human intelligence. Fama (1965) stated that the efficiency of the securities market in the strong or semi-strong sense causes instantaneous changes in traders' subjective equilibrium with the incorporation of securities information, and the speed at which they respond to this information determines the winners and losers. In the weak sense, market equilibrium follows changes in traders' subjective equilibria due to noise. The efficiency of the securities market can be summarized as the correct and rapid incorporation of information into prices, the presence of arbitrage opportunities, countless traders always searching for such opportunities, the gradual loss of arbitrage opportunities, and the market becoming more efficient. However, it should be noted that without noise, transactions would not occur, and the price distortion would persist in the presence of strong bias.

    In an efficient market where noise and bias have no effect, and information is perfectly symmetrical, security prices should accurately reflect only information. However, if traders had perfectly and simultaneously symmetrical knowledge and information, including the asset valuation model, the transaction would not function because the traders' valuation of the asset would be the same. They must be doing some rational calculations via intelligence in natural and digital computing. In this sense, a rational representative agent in macroeconomics corresponds to the perfect symmetry of information (Lucas (1976), Kirman (1992), Hartley (1996)).

    While Black (1986) treated this "symmetry breaking" as noise, the effect of noise on a security's price is expected to be symmetrical based on its nature. However if there exist so many irrational noise traders synchronizing erroneous stochastic beliefs that both affect prices and earn higher expected returns, the unpredictability of noise traders' beliefs creates a risk in the price of the asset that deters rational arbitrageurs from aggressively betting against them. As a result, prices can diverge significantly from fundamental values even in the absence of fundamental risk. Moreover, bearing a disproportionate amount of risk that they themselves create enables noise traders to earn a higher expected return than rational investors (De Long et al., (1990)).

    On the other hand, Tversky and Kahneman (1974) found that the effect of bias is asymmetric. They described three heuristics that are employed in making judgments under uncertainty: 1) representativeness, 2) availability of instances or scenarios, and 3) adjustment from an anchor. These heuristics are highly economical and usually effective, but they lead to systematic and predictable errors. However, the effect of these heuristics has not been detected as a global bias in the securities market. Consequently, identifying the specific effects of noise and bias on security prices is challenging.

    If attention is paid to any statistical property in any complex system, the log-normal distribution is the most natural and appropriate among the standard or "normal" statistics to overview the whole system (Kobayashi et al., (2011)). Log-normality emerges as a familiar and typical example of statistical aspects in various complex systems. Since every member of any complex system has its own history, each member is in the process of growth (or retrogression). The log-normal distribution is realized because of Gibrat's law, or the Matthew effect. It is applied to cities' size and growth rate, where the proportionate growth process may cause a distribution of city sizes that is log-normal. When considering the entire size distribution, not just the largest cities, the city size distribution is log-normal (Samuels (1965)). However, it has been argued that it is problematic to define cities through their fairly arbitrary legal boundaries. According to Gabaix (1999), Zipf's law is a very tight constraint on the class of admissible models of local growth. It says that for most countries, the size distribution of cities strikingly fits a power law: the number of cities with populations greater than S is proportional to 1/S. Suppose that, at least in the upper tail, all cities follow some proportional growth process (this appears to be verified empirically). This automatically leads their distribution to converge to Zipf's law.

    According to Gibrat's law of proportionate effect, the relative change in a firm's size is not affected by its absolute size, which implies that both small and large firms experience similar average rates of growth. However, Samuels (1965) provides evidence that contradicts this law by showing that larger firms grow significantly faster. Aoki and Nirei (2017) developed a neoclassical growth model that produces Pareto's law of income distribution and Zipf's law of firm size distribution based on firm-level productivity shocks. Executives and entrepreneurs invest in both risk-free assets and their own firms' risky stocks, which affect their income and wealth. They used the model to investigate how changes in tax rates can explain the evolution of top incomes in the US. The model accounts for the recent decline in the Pareto exponent of the income distribution and the trend of the top 1 percent income share. In the same research direction, Nirei and Aoki (2016) developed a neoclassical growth model that incorporates heterogeneous households and explains the Pareto distributions of income and wealth in the upper tail. They introduced households' business productivity risks and borrowing constraints in a standard Bewley (1977) model to generate Pareto distributions. Low-productivity households rely on wages and safe asset returns, whereas high-productivity households do not diversify their business risks. The model can accurately explain the observed income distribution in the United States with reasonable parameter values. The authors conducted comparative statics to investigate how changes in parameters affect Pareto distributions and found that the increase in top income dispersion in recent decades may be attributed to changes in the top tax rates in the 1980s. Their analytical finding offers a consistent explanation for numerical comparative statics.

    Few studies are available on global bias in the stock market. The fundamental indicators of stocks include information about the effects of noise and bias on stock prices; however, distinguishing between them is generally hard. In this article, I present the fundamentals hypothesis based on rational expectations (Muth (1961)) and, using a log-normal distribution model, detect global bias components from the price-earnings (P/E), price-to-book (P/B) and price-to-cash flow (P/CF) ratios. The traditional theory of the firm assumes that the firm acts in the stockholders' interests and that stockholders are interested in profits such that the object of the firm is to maximize profit. However, there is a range in the profit concept (Bodenhorn (1964)). The analysis results support our fundamentals hypothesis as the detected biases show similar characteristics. Additionally, the results show that the cash flow indicators contain relatively few bias components and are closer to the fundamentals. I further demonstrate and examine why the positive P/IC ratio among the indicators analyzed is a proxy for the fundamentals that do not include bias components. Globally, this implies that the efficient market hypothesis holds. The shape of the curve, and in particular the strength of the bias, is stable throughout time and independent of whether the economy is good or bad. The answer is simple: Cash is a fact, and profit is an opinion. Namely, opinions of management and accountants are added as noise to fundamentals.

    When Xt denotes the true fundamentals of listed companies at time t, and Rt denotes their growth rates, the following Gibrat process represents the growth of those companies

    Xt=RtXt1. (1)

    It is important to note that I assume that the growth rates Rt are mutually independent random variables that follow the same distribution with finite variance. The initial value of the fundamentals being set as X0 yields

    XT=X0Tt=1Rt (2)

    at time T. Taking the log of both sides of the equation results in

    logXT=logX0+logR1++logRT. (3)

    Therefore,

    logXTLN(μ,σ2) (4)

    would hold for a sufficiently large T based on the central limit theorem. Essentially, the true fundamental XT of listed companies follows the log-normal distribution.

    Furthermore, by assuming rational expectations through the future point in time T as of the present point in time 0 on the premise of a going concern, the following equation becomes true:

    X0=E[XT]Tt=1E[R1t]. (5)

    Therefore, the rational expectations X0 for fundamentals follow the log-normal distribution. In other words, the fundamentals Xt at time t follow the log-normal distribution.

    Among the actual fundamental indicators, we analyze the price-earnings (P/E), forward price-earnings (P/FE), price-to-book (P/B), price-to-cash flows from operating activities (P/OC), price-to-cash flows from investing activities (P/IC) price-to-cash flows from financing activities (P/FC), and the price-to-cash equivalents at the year-end (P/CE) ratios. We divide the cash flow ratios into positive and negative data. We use daily data with absolute values less than 1,000 for all companies listed in Japan for the period spanning 1,817 business days from January 2007 to May 2014 (data source: Nikkei NEEDS https://nkbb.nikkei.co.jp/en/service/nikkei-needs/).

    If our fundamentals hypothesis is true, the actual fundamental indicators should follow the log-normal distribution. To test our hypothesis, we calculate the p-values using the Kolmogorov-Smirnov test, Pearson's χ2 test and Anderson-Darling test, to test the goodness of fit to the log-normal distribution.

    Our tests can be strongly affected by the sample size that differs by indicator, so we extract samples at the 300 quartiles from each data set to calculate the p-value for testing the goodness of fit of the overall average of each indicator to the log-normal distribution with variance. For all cases, the null hypothesis is that "the indicator follows the log-normal distribution, " or the otherwise worded "the indicator is a proxy of true fundamentals." The indicator reflects the fundamentals more when the p-value is closer to 1, and when p-value is closer to 0, the bias is stronger. We show the lowest p-value for each test in Table 1 for the data spanning 1,817 business days.

    Table 1.  The lowest p-values for the data spanning 1,817 business days. Samples at the 300 quartiles. Significant (bias) levels: 1%(***), 3%(**), 5%(*).
    Kolmogorov-Smirnov Pearson's χ2 Anderson-Darling
    P/OC+ 0.0986 0.2843 0.0398
    P/OC 0.0634 0.0045 0.0850
    P/IC+ 0.6205 0.4441 0.7201
    P/IC 0.2787 0.7012 0.1849
    P/FC+ 0.2042 0.2287 0.1093
    P/FC 0.2522 0.5763 0.1269
    P/CE 0.4611 0.8740 0.2376
    P/E 0.0170 0.0000 0.0018
    P/FE 0.0107 0.0040 0.0019
    P/B 0.0032 0.0023 0.0005

     | Show Table
    DownLoad: CSV

    The test results indicate that bias is strong on the P/E, P/FE, and P/B ratios. The forward price-earnings ratio (P/FE) is the most biased. Among all price-to-cash flow ratios, the bias is significant only for the P/OC ratio. It can be easier to manipulate than others. I illustrate the time series of p-values since the p-value is the indicator of the strength of bias. From Table 1, Figure 1, and Figure 2, the positive P/IC ratio is the best proxy for the fundamentals that do not include bias components. Probably, it should be hard to manipulate according to the opinions of management and accounting.

    Figure 1.  Time series of the goodness of fit test to the log-normal distribution. Blue: Kolmogorov-Smirnov, Yellow: Pearson's χ2, Blue: Anderson-Darling. The higher the value, the less bias there is, and vice versa.
    Figure 2.  The more the data deviates from a straight line showing a log-normal distribution, the stronger the bias. These graphs are drawn with data from the first day of the time series. The shape of the bias varies proportionally to Figure 1, but is stable throughout the time. https://figshare.com/projects/Intelligence_and_Global_Bias_in_the_Stock_Market/131711.
    Figure 3.  Comparison of three distributions: Log-Normal (blue), GDP-1 (yellow), GDP-2 (green).
    Figure 4.  Left: GDP-1 to fit P/FE (κ=18.82,α=1,γ=0.385,μ=0.993), Right: GDP-2 to fit P/FE (κ=13.70,α=0.515,γ=0.238,μ=0.993). Note that P/FE has the strongest bias.
    Figure 5.  Number of all listed companies in the Japanese stock market (data source: Nikkei NEEDS https://nkbb.nikkei.co.jp/en/service/nikkei-needs/).

    The detected biases exhibit similar characteristics. Each indicator is shown in Figure 2 and compared to the log-normal distribution. While the figures only show the data for the first day of all periods, January 4, 2007, the biases remain unchanged. It is worth noting that the shape and strength of the curve are independent of whether the economy is good or bad. The artificial stock market appears to have natural intelligence.

    The results indicate that the fundamental indicators, including the P/E, P/FE, and P/B ratios, are strongly affected by bias. Bias also has a significant effect on the positive and negative P/OC ratios. Additionally, there is a weak bias on the negative P/IC, positive and negative P/FC, and P/CE ratios.

    When we compare the test results, the positive P/IC ratio is the stable proxy of fundamentals among all these indicators. Positive cash flows from investing activities represent the realized gain or loss from past investments such as marketable securities, tangible fixed assets, the sales of investment securities, and income from the collection on loans declared at the end of the accounting period. In other words, positive cash flow is the indicator that most directly reflects past business decisions. We can interpret the year-end cash equivalent ratio and the other cash flow ratios reflecting the fundamentals because they have less bias of the components.

    On the other hand, although the P/E ratio and P/B ratio are definite values, investors might not view them as indicators that reflect the fundamentals because there is a high degree of freedom in accounting. The P/OC ratio might also have a lower credibility than other cash flow indicators.

    The answer is simple: "Cash is a fact, and profit is an opinion." Namely, opinions of management and accountants are added as noise to true fundamentals. As a result, the Kesten (1973) process

    Xt=RtXt1+ϵt,E[ϵt]>0, (6)

    is realized so that the Pareto distribution is obtained. This result means that a positive bias accompanied their opinions.

    In fact, these biases fit the Pareto distribution quite well. The following functions represent generalized Pareto distribution (GPD).

    F(x)=1(1+(xμκ)1/γ)α (7)
    f(x)=αγκ1/γ(xμ)1+1γ(1+(xμκ)1/γ)α1 (8)

    What are the implications of the results? Does the existence of bias negate the efficient market hypothesis? From the viewpoint of intelligence in natural and digital computing, how should we deal with this problem? The artificial stock market, where countless traders use digital computers to trade, appears to have natural intelligence. Intelligence does not arise only in individual brains but also arises in groups of individuals. The interest in whether markets are efficient is based on the idea that inefficiencies are the reason active returns are possible. If our fundamentals hypothesis is true, the fundamental indicators, including the P/E, P/FE, and P/B ratios, are strongly affected by bias. This result means that investors might not view them as indicators that reflect the fundamentals since they have a high degree of freedom in accounting. Can we conclude that the market is inefficient because of this? This issue has important economic implications. Globally, our result implies that the efficient market hypothesis holds. The shape of the curve, and in particular the strength of the bias, is stable throughout time and independent of whether the economy is good or bad. It is the path of the ants (Kirman (1993); Sano (2022)) to avoid obstacles, and traders move back and forth along it in response to noise (Black (1986); De Long et al., (1990)). This noise, among other things, plays a role in activating the market system, which knows the fundamentals of firms and represents a stable global bias in the stock market. Nevertheless, when the biases revealed by this study disappear, the stock market will become more efficient. It is the responsibility of firms and financial authorities to disclose fair information.

    The author declares he has not used Artificial Intelligence (AI) tools in the creation of this article.

    The author is grateful to the Department of Economics of Fukui Prefectural University. This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

    The author declares no conflicts of interest in this paper.



    [1] J. Steele, Stability of Plankton Ecosystem, Chapman and Hall, London, 1974.
    [2] T. Saha, M. Bandyopadhyay, Dynamical analysis of toxin producing phytoplankton-zooplankton interactions, Nonlinear Anal. Real World Appl., 10 (2009), 314-332.
    [3] Y. Lv, Y. Pei, S. Gao, C. Li, Harvesting of a phytoplankton-zooplankton model, Nonlinear Anal. Real World Appl., 11 (2010), 3608-3619.
    [4] C. Liu, L. Wang, Q. Zhang, Y. Yun, Dynamical analysis in a bioeconomic phytoplankton zooplankton system with double time delays and environmental stochasticity, Phys. A, 482 (2017), 682-698.
    [5] X. Meng, J. Li, Stability and Hopf bifurcation analysis of a delayed phytoplankton-zooplankton model with Allee effect and linear harvesting, Math. Biosci. Eng., 17 (2019), 1973-2002.
    [6] G. Denaro, D. Valenti, B. Spagnolo, G. Basilone, S. Mazzola, S. W. Zgoz, et al., Dynamics of two picophytoplankton groups in mediterranean sea: Analysis of the deep chlorophyll maximum by a stochastic advection-reaction-diffusion model, PLoS One, 8 (2013), e66765.
    [7] G. Denaro, D. Valenti, A. La Cognata, B. Spagnolo, A. Bonanno, G. Basilone, et al., Spatio-temporal behaviour of the deep chlorophyll maximum in Mediterranean Sea: Development of a stochastic model for picophytoplankton dynamics, Ecol. Complex, 13 (2013), 21-34.
    [8] D. Valenti, G. Denaro, A. L. Cognata, B. Spagnolo, Picophytoplankton dynamics in noisy marine environment, Acta. Phys. Pol B, 43 (2012), 1227-1240.
    [9] A. S. Heiskanen, T. Tamminen, K. Gundersen, Impact of planktonic food web structure on nutrient retention and loss from a late summer pelagic system in the coastal northern Baltic Sea, Mar. Ecol. Prog. Ser., 145 (1996), 195-208.
    [10] M. Scheffe, Fish and nutrients interplay determines algal biomass: a minimal model, Oikos, 62 (1991), 271-282.
    [11] A. B. Medvinsky, S. V. Petrovskii, I. A. Tikhonova, H. Malchow, B. L. Li, Spatio-temporal complexity of plankton and fish dynamics, SIAM Rev., 44 (2002), 311-370.
    [12] P. Panja, S. K. Mondal, Stability analysis of coexistence of three species prey-predator model, Nonlinear Anal., 81 (2015), 373-382.
    [13] A. Sharma, A. K. Sharma, K. Agnihotri, Complex dynamic of plankton-fish interaction with quadratic harvesting and time delay, Model Earth Syst. Environ., 2 (2016), 1-17.
    [14] X. Y. Meng, Y.Q. Wu, Bifurcation and control in a singular phytoplankton-zooplankton-fish model with nonlinear fish harvesting and taxation, Int. J. Bifurcat. Chaos, 28 (2018), 1850042.
    [15] Z. Wei, J. Sugie, Global asymptotic stability and equiasymptotic stability for a time-varying phytoplankton-zooplankton-fish system, Nonlinear Anal. Real World Appl., 46 (2019), 116-136.
    [16] B. Buonomo, M. Cerasuolo, The effect of time delay in plant-pathogen interactions with host demography, Math. Biosci. Eng., 12 (2015), 473-490.
    [17] B. Tian, Y. Qiu, N. Chen, Periodic and almost periodic solution for a non-autonomous epidemic predator-prey system with time-delay, Appl. Math. Copmut., 215 (2009), 779-790.
    [18] O. A. Chichigina, A. A. Dubkov, D, Valenti, B. Spagnolo, Stability in a system subject to noise with regulated periodicity, Phys. Rev. E, 84 (2011), 021134.
    [19] D. Valenti, L. Tranchina, M. Brai, A. Caruso, C. Cosentino, B. Spagnolo, Environmental metal pollution considered as noise: Effects on the spatial distribution of benthic foraminifera in two coastal marine areas of Sicily (Southern Italy), Ecol. Model, 213 (2008), 449-462.
    [20] A. A. Dubkov, B. Spagnolo, Verhulst model with Lévy white noise excitation, Eur. Phys. J. B, 65 (2008), 361-367.
    [21] Q. Luo, X. Mao, Stochastic population dynamics under regime switching, J. Math. Anal. Appl., 355 (2009), 577-593.
    [22] M. Slatkin, The dynamics of a Population in a Markovian environment, Ecology, 59 (1978), 249-256.
    [23] H. Qiu, W. Deng, Stationary distribution and global asymptotic stability of a three-species stochastic food-chain system, Turk. J. Math., 41 (2017), 1292-1307.
    [24] Y. Ma, Q. Zhang, L. Wang, T. Kang, Dissipative control of a three-species food chain stochastic system with a hidden Markovchain, Adv. Differ. Equ-Ny., 2017 (2017), 1-22.
    [25] Y. Lin, D. Jiang, Long-time behavior of a stochastic predator-prey model with modified Leslie-Gower and Holling-type II schemes, Int. J. Biomath, 9 (2016), 1650039.
    [26] J. Lv, K. Wang, Asymptotic properties of a stochastic predator-prey system with Holling II functional response, Commun. Nonlinear. Sci. Numer Simulat., 16 (2011), 4037-4048.
    [27] Z. Liu, N. Shi, D. Jiang, C. Ji, The Asymptotic behavior of a stochastic predator-prey system with Holling II functional response, Abstr. Appl. Anal., 2012 (2012), 1-14.
    [28] G. Gilioli, S. Pasquali, F. Ruggeri, Nonlinear functional response parameter estimation in a stochastic predator-prey model, Math. Biosci. Eng., 9 (2012), 75-96.
    [29] C. S. Reynolds, The Ecology of Freshwater Phytoplankton, Cambridge University Press, Cambridge, 1984.
    [30] C. W. Clark, Mathematical Bio-Economics: The Optimal Management of Renewable Resources, Wiley, New York, 1976.
    [31] C. W. Clark, Bioeconomic Modeling and Resource Management, in Applied Mathematical Ecology (eds. S. A. Levin), Springer, 1989, 11-57.
    [32] M. Mesterton-Gibbons, On the optimal policy for combining harvesting of predator and prey, Nat. Resour. Model, 3 (1988), 63-90.
    [33] M. Mesterton-Gibbons, A technique for finding optimal two-species harvesting policies, Ecol. Model, 92 (1996), 235-244
    [34] S. Wang, L. Wang, T. Wei, Optimal harvesting for a stochastic predator-prey model with S-type distributed time delays, Methodol Comput. Appl., 20 (2016), 37-68.
    [35] M. Liu, C. Bai, Analysis of a stochastic tri-trophic food-chain model with harvesting, J. Math. Biol., 73 (2016), 597-625.
    [36] M. Liu, X. He, J. Yu, Dynamics of a stochastic regime-switching predator-prey model with harvesting and distributed delay, Nonlinear Anal. Hybrid Syst., 28 (2018), 87-104.
    [37] W. J. Anderson, Continuous-Time Markov Chains, Springer, New York, 1991.
    [38] R. M. May, Stability and Complexity in Model Ecosystems, Princeton University Press, Princeton, 1975.
    [39] X. Zhang, W. Li, M. Liu, K. Wang, Dynamics of a stochastic Holling II one-predator two-prey system with jumps, Phys. A, 421 (2015), 571-582.
    [40] X. Mao, C. Yuan, Stochastic Differential Equations with Markovian Switching, Imperial College Press, London, 2006.
    [41] M. Liu, J. Yu, P. S. Mandal, Dynamics of a stochastic delay competitive model with harvesting and Markovian switching, Appl. Math. Comput., 337 (2018), 335-349.
    [42] V. M. Popov, Hyperstability of Control Systems, Springer-Verlag, New York, 1973.
    [43] D. Prato, J. Zabczyk, Ergodicity for Infinite Dimensional Systems, Cambridge University Press, Cambridge, 1996.
    [44] L. Thomas, Estimating Phytoplankton Growth Rates from Compositional Data, in Oceanography/Biological Oceanography Massachusetts Institute of Technology and Woods Hole Oceanographic Institution, Massachusetts Institute Of Technology, 2008.
    [45] T. Nanazato, M. Yasuno, Population dynamics and production of cladoceran zooplankton in the highly eutrophic Lake Kasumigaura, Hydrobiologia, 124 (1981), 13-22.
    [46] Y. Wang, Q. Liu, Estimating natural mortality from stock size and catch data (in Chinese), Period. Ocean Univ. China, 35 (2005), 020-024.
  • This article has been cited by:

    1. Muhammad Waseem Sajjad, Fatima Muzamil, Maida Sabir, Usman Ali Ashfaq, Regenerative Medicine and Nanotechnology Approaches against Cardiovascular Diseases: Recent Advances and Future Prospective, 2025, 20, 1574888X, 50, 10.2174/011574888X263530230921074827
  • Reader Comments
  • © 2020 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(4654) PDF downloads(113) Cited by(8)

Figures and Tables

Figures(4)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog