Research article

Optimal strategies for a fishery model applied to utility functions

  • Received: 05 September 2020 Accepted: 22 November 2020 Published: 11 December 2020
  • This work examines aquaculture-related activities in the commercial exploitation of fish reproduction. Fisheries' problem of maximizing utility is modeled for the state of Puebla, Mexico, to determine optimal fish production. The problem of maximizing utility subject to the fish production function is solved using an approach based on Euler's equation. The theoretical results are then applied, using data on aquaculture production and tilapia sales prices in the state of Puebla, Mexico. A logarithmic regression is used to approximate the utility function. The optimal fishing production and utility functions are thus explicitly obtained. Furthermore, this work shows how to obtain greater profits from the amount of fish that can be extracted without reducing the fish population.

    Citation: Carlos Camilo-Garay, R. Israel Ortega-Gutiérrez, Hugo Cruz-Suárez. Optimal strategies for a fishery model applied to utility functions[J]. Mathematical Biosciences and Engineering, 2021, 18(1): 518-529. doi: 10.3934/mbe.2021028

    Related Papers:

    [1] Hua Guo, Yuan Tian, Kaibiao Sun, Xinyu Song . Dynamic analysis of two fishery capture models with a variable search rate and fuzzy biological parameters. Mathematical Biosciences and Engineering, 2023, 20(12): 21049-21074. doi: 10.3934/mbe.2023931
    [2] Yuan Tian, Yan Gao, Kaibiao Sun . A fishery predator-prey model with anti-predator behavior and complex dynamics induced by weighted fishing strategies. Mathematical Biosciences and Engineering, 2023, 20(2): 1558-1579. doi: 10.3934/mbe.2023071
    [3] Ali Moussaoui, Pierre Auger, Christophe Lett . Optimal number of sites in multi-site fisheries with fish stock dependent migrations. Mathematical Biosciences and Engineering, 2011, 8(3): 769-783. doi: 10.3934/mbe.2011.8.769
    [4] Xin Gao, Yue Zhang . Bifurcation analysis and optimal control of a delayed single-species fishery economic model. Mathematical Biosciences and Engineering, 2022, 19(8): 8081-8106. doi: 10.3934/mbe.2022378
    [5] Andras Balogh, Tamer Oraby . Stochastic games of parental vaccination decision making and bounded rationality. Mathematical Biosciences and Engineering, 2025, 22(2): 355-388. doi: 10.3934/mbe.2025014
    [6] Summer Atkins, William W. Hager, Maia Martcheva . The switch point algorithm applied to a harvesting problem. Mathematical Biosciences and Engineering, 2024, 21(5): 6123-6149. doi: 10.3934/mbe.2024269
    [7] Shixuan Yao, Xiaochen Liu, Yinghui Zhang, Ze Cui . An approach to solving optimal control problems of nonlinear systems by introducing detail-reward mechanism in deep reinforcement learning. Mathematical Biosciences and Engineering, 2022, 19(9): 9258-9290. doi: 10.3934/mbe.2022430
    [8] Wei Wu, Zhiyi Lin, Ming Wei . Research on a collaborative evolution model of multi-airport route network considering subsidy strategies. Mathematical Biosciences and Engineering, 2023, 20(11): 19808-19838. doi: 10.3934/mbe.2023877
    [9] Pierre Auger, Tri Nguyen-Huu, Doanh Nguyen-Ngoc . On the impossibility of increasing the MSY in a multisite Schaefer fishing model. Mathematical Biosciences and Engineering, 2025, 22(2): 415-430. doi: 10.3934/mbe.2025016
    [10] Yuqing Lu . A knowledge-driven decision support system for remote medical management. Mathematical Biosciences and Engineering, 2023, 20(2): 2732-2749. doi: 10.3934/mbe.2023128
  • This work examines aquaculture-related activities in the commercial exploitation of fish reproduction. Fisheries' problem of maximizing utility is modeled for the state of Puebla, Mexico, to determine optimal fish production. The problem of maximizing utility subject to the fish production function is solved using an approach based on Euler's equation. The theoretical results are then applied, using data on aquaculture production and tilapia sales prices in the state of Puebla, Mexico. A logarithmic regression is used to approximate the utility function. The optimal fishing production and utility functions are thus explicitly obtained. Furthermore, this work shows how to obtain greater profits from the amount of fish that can be extracted without reducing the fish population.


    Aquaculture has had great success in the fishing subsector in recent years, using controlled reproduction through breeding and cultivation techniques that are amenable to commercial exploitation. Aquaculture has great potential in Mexico, since it has the necessary features for generating jobs and commercially exploiting different species of fish, given rising demand that is expected to increase even further over the following decades [1]. In order of importance, the main aquaculture species harvested in Mexico are mojarra, shrimp, oyster, carp, and trout.

    In particular, Puebla has been one of the most successful states in aquaculture development. This progress has led to the consolidation of aquaculture and new opportunities for companies and individuals dedicated to tilapia production. The state of Puebla has no coastline, another factor that increases the importance of tilapia production, since tilapia easily adapt to captivity. In addition, aquaculture activity is also being used in to reduce poverty and develop new sources of capital income.

    Given the importance of aquaculture in Mexico, this study models the fishery problem as a Markov decision process (MDP) [2,3], [4] with a state space related to the fish production, a set of actions involving the quantity of fish removed from the pond, a set of admissible actions that are a function of fish production, the transition law, and the utility function. In the previously approach, the main objective of a country, region, or company was to maximize the profits from its fish production, that is, profit maximization subject to a transition law. Furthermore, to obtain greater benefits from the amount of fish that can be extracted without reducing the population of fish. In this document, a discrete-type dynamic is proposed for the study of the fish population. On the other hand, in the literature on the subject, there exist references about continuous time dynamic population, for example, see [5] and [6]. The fisheries problem has been previously addressed in the context of dynamic programming [7,8]; however, those works do not use the Euler's Equation (EE) to find the optimal solution [9]. In the present work, the optimal solution is characterized by and derived explicitly from the EE.

    After modeling the fisheries problem, this study demonstrates that the dynamic programming conditions are satisfied, guaranteeing the existence of an optimal solution. In this case, each of the conditions necessary for the existence of the optimal solution is demonstrated in detail. However, this only proves the existence of an optimal policy; next, the EE is used to determine this optimal policy and the optimal value function.

    Finally, the parameters are estimated using fish production statistics for the state of Puebla, Mexico. To this end, a logarithmic regression is used to approximate the utility function, following the methodology exposed in [10]. Thus, this model is used to determine the optimal policy and optimal value solution numerically and thus make comparisons with previous work.

    This paper is organized as follows. First, in section materials and method is presented the statistical analysis and later, the topics of Markov decision processes and EE are commented in a general way. In section 3, the main results of the manuscript are presented, in this section the existence of optimal extraction strategies is guaranteed and they are characterized using the Euler equation. Finally, the work presents some numerical results and a discussion section of the results obtained.

    This work uses the tilapia catch data from CONAPESCA [11] between 2006 and 2014 for the state of Puebla, Mexico. Figure 1 plots the weight (in kilograms) of fish caught versus total sales (in Mexican pesos) per month, over 102 months. Following CONAPESCA, the total sales data coincide with the demand for tilapia in the region. Therefore, fish demand is related to sales data.

    Figure 1.  Amount of fish caught versus total sales per month.

    In this case, the plot in Figure 1 has been adjusted to fit a logarithmic utility. Table 1 presents the results.

    Table 1.  Descriptive statistics of the estimation.
    Parameter Estimator Standard error Confidence Interval of 95%
    b1 –4,697,590 1,149,049 (–6,977,269, –2,417,909)
    b2 567,437 109,491 (350,210,784,665)

     | Show Table
    DownLoad: CSV

    Figure 2 shows that the data fit a logarithmic function, which is a utility function that expresses the total amount sold in Mexican pesos, U(a)=b1ln(a)b2 (with conditions satisfied by a utility function [12,13], where b1:=567,437, b2:=4,697,590, and a is the amount (in kilograms) of fish harvested. Note that, in this case, a is always a non-negative number.

    Figure 2.  Logarithmic regression: kilograms of fish caught versus sales (in pesos).

    In the next section, the problem of fisheries is presented in the context of an MDP, where the state space is related to fish production. The set of actions is the amount of fish extracted, and the set of admissible actions is a function of fish production. Then, the main objective of a country, region, or company is to maximize its profits, according to its fish production, that is, profit maximization subject to a transition law. To this end, consider h as the fish production function, which denotes the quantity of fish produced in a certain period (for the conditions that must satisfy a production function see [12]). In particular, let h(x):=xγ, where 0<γ<1 is the production factor. This class of function is widely used in fish production, see e.g., [7,8]. Then, to solve this optimal control problem, a discrete-time version of the Euler-Lagrange equation (see chapter 6 of [14]) is deduced. In the context of MDP, this equation is simply known as EE, a deterministic version of EE can be consulted in [15] and a stochastic version for linear quadratic models is proposed in [16].

    The fishery model can be interpreted as follows. Let xt be the amount of fish at time t, where xt denotes the state of the system at time t. Thus, X=[0,) is the state space and at denotes the amount of fish caught. In this case, the action and the state spaces are the same, that is, A=[0,). The production function, h:X[0,), is given by h(x)=xγ. This implies that the set of admissible actions, when the system is in state x, is A(x)=[0,xγ] for each state xX. Furthermore, K:={(x,a)xX,aA(x)}. The evolution equation of the fishery process thus has the following form:

    xt+1=(xtat)γξt (3.1)

    where, for t=0,1,, ξ0,ξ1,, are nonnegative independent and identically distributed random variables, independent of the initial stock x0, that represent the fish mortality rate, and their probability distribution function is denoted by F, that is, F(s)=P(ξs) for all sR and density ΔC2([0,)). The expected value E[ξ]:=μ<+, and ξ represents a generic element of the sequence {ξt}. Finally, the utilities are represented by the following function U:AR:

    U(a)={0ifexp(b2/b1)>a,b1ln(a)b2ifexp(b2/b1)a.

    Remark 3.1. In this case, the transition law Q is determined for xX, aA(x), and BB(X) (where B(X) is the Borel σ-algebra of X) in the following form:

    Q(Bxt=x,at=a)=Bω(x,y,a)dy,

    where the function ω:[0,)3[0,) is defined by

    ω(x,y,a):=Δ(y(xa)γ)1(xa)γ

    for x,yX and a[0,xγ).

    The MDP M:=(X,A,{A(x)xX},Q,U) is thus defined [3], and M is referred to as the fishery model.

    Definition 3.2. For each t=0,1,, define the space Ht of admissible histories up to time t as H0=X, and Ht=K×Ht1 for t=1,2,. A policy is defined as a sequence π={πt,t=0,1,} of stochastic kernels, defined on A given Ht [3]. In this paper, the set of policies is denoted by Π. In particular, a stationary policy is of the form π={f,f,}, where f is defined as a measurable function f:XA such that f(x)[0,xγ] for all xX. The set of all stationary policies is denoted by F.

    The sequence of consecutive states and corresponding actions are denoted by {xt} and {at}, respectively, and Eπx denotes the expected value with respect to the probability measure Pπx, which is defined on a canonical measurable space (Ω,F) according to the Ionescu–Tulcea theorem [3].

    For each πΠ and initial state xX, let

    V(π,x)=Eπx[t=0αtU(at)]

    be the total expected discounted utility when policy π is applied, given initial state x. The constant α(0,1) is a fixed discount factor.

    Definition 3.3. A policy π is optimal if V(π,x)=υ(x) for all xX, where

    υ(x):=supπΠV(π,x),    xX

    and υ is called the optimal value function.

    Let M be the fishery model described in the previous section. In this section, the necessary conditions to guarantee the existence of an optimal policy in the fishery model M will be proven. It is worth mentioning that the proof is carried out theoretically, so that, in the next sections, this policy can be calculated explicitly and, consequently, the optimal value function can be determined.

    Lemma 3.4. There exists a sequence {Xj}jN0 (where N0 is the set of natural numbers union zero) of non-empty Borel subsets X such that XjXj+1 for all jN0 and

    X=j0int(Xj),

    where int(Xj) denotes the interior of set Xj for each jN0. In addition, for any jN0, let mj be

    mj=supxXjsupaA(x)U+(a), (3.2)

    where U+(a)=max{U(a),0}. Assume that m0>0 and mj< for every jN0 and

    αlim supjmj+1mj<1. (3.3)

    Finally, for each jN0 and xXj, aA(x),

    Q(Xj+1x,a)=1. (3.4)

    Proof. Let ε>0 be a fixed number and define Xj:=[0,j+ε], for each jN0. Note that XjXj+1 and

    X=j0int(Xj).

    Furthermore,

    mj=supxXjsupaA(x)U+(a)=supx[0,j+ε]supa[0,xγ]b1ln(a)b2=supx[0,j+ε]γb1ln(x)b2=γb1ln(j+ε)b2

    for each jN0 and aexp(b2/b1). On the other hand, substituting mj and mj+1 in Eq. (3.3) and applying L'Hˆopital's rule yields

    αlim supjmj+1mj=α1+ε2+ε<α<1.

    Finally, for t=0,1,2,, if xt=xX and atA(x), there exists jN0 such that xXj, and, since ξ(0,1), thus xt+1Xj+1 almost surely with respect to the probability measure Pπx, that is, Q(Xj+1x,a)=1.

    Definition 3.5. The value iteration functions are defined as

    Vn(x)=maxaA(x){b1ln(a)b2+αE[Vn1((xa)γξ)]} (3.5)

    for all xX and n=1,2,, with V0()=0.

    Using Lemma 3.4, the proof of the following lemma can be adapted from the proof of Proposition 1 and Theorem 1 by Jaskiewicz and Nowak [17].

    Lemma 3.6. The fishery problem satisfies the following conditions.

    a) The optimal value function υ:XR is a measurable function and a solution of

    υ(x)=maxaA(x){b1ln(a)b2+αE[υ((xa)γξ)]} (3.6)

    for each xX.

    b) There exists a selector fF such that, in Eq. (3.6), the maximum is attained, that is,

    υ(x)=b1ln(f(x))b2+αE[υ((xf(x))γξ)] (3.7)

    for all xX and f is optimal.

    c) Vn:XR is a measurable function, n=1,2,, and {Vn} converges pointwise to υ. Moreover, for each n=1,2,, there exists fnF such that

    Vn(x)=b1ln(fn(x))b2+αE[Vn1((xfn(x))γξ)]

    for all xX, where fn is the maximizer of the value iteration functions.

    We define the functions

    G(x,a):=b1ln(a)b2+αE[υ(((xa)γ)ξ)]

    and

    Gn(x,a):=b1ln(a)b2+αE[Vn1(((xa)γ)ξ)],

    n=1,2,, for each (x,a)K.

    Theorem 3.7. The iterative value functions satisfy the functional equation

    Vn(x)=αγE[Vn1((xVn(x)b1Vn(x))γξ)ξ](xVn(x)b1Vn(x))γ1 (3.8)

    for each n2 and xint(X), where ξ is a generic element of the sequence {ξt} (see Eq.3.1). This equation is known in the literature of MDPs as Euler equation.

    Proof. Deriving the function Gn with respect to a yields

    Gna(x,a)=b1aαγE[Vn1((xa)γξ)ξ](xa)γ1,

    (x,a)K. Evaluating the preceding equation in the maximizer fnF results in

    b1fn(x)αγE[Vn1((xfn(x))γξ)ξ](xfn(x))γ1. (3.9)

    Furthermore,

    Vn(x)=αγE[Vn1((xfn(x))γξ)ξ](xfn(x))γ1, (3.10)

    xint(X), and substituting Eq. (3.9) in Eq. (3.10) yields

    b1fn(x)=Vn(x),

    xint(X). Then,

    fn(x)=b1Vn(x). (3.11)

    Finally, substituting Eq. (3.11) in Eq. (3.10) yields Eq. (3.8).

    Theorem 3.8. The optimal value function and the optimal policy are given by

    υ(x)=b1Klnx+C

    and

    f(x)=(1αγ)x,

    respectively, for each xX, where

    K=11αγ,
    C=11α[b1αγ1αγln(αγ)+b1ln(1αγ)+b1αμ1αγb2].

    Proof. First, note that

    V1(x)=b1γlnxb2,

    xX, it follows that

    V1(x)=b1γx,

    x(0,+). For n=2, using Eq. (3.8) and the preceding equation, yields

    V2(x)=αE[V1((xV2(x)b1V2(x))γξ)ξ]γ(xV2(x)b1V2(x))γ1=b1αγ2V2(x)xV2(x)b1.

    Therefore,

    V2(x)=b1(αγ2+1)x,

    x(0,+). It will be proven by induction that, for all n2, Vn(x)=b1Knx. To this end, suppose that

    Vn1(x)=b1Kn1x, (3.12)

    x(0,+), where Kn1=n2i=0(αγ)i+αn2γn1. Thus, substituting Eq. (3.12) into Eq. (3.8) yields

    Vn(x)=b1αγKn1Vn(x)xVn(x)1,

    from which it follows that

    Vn(x)=b1αγKn1+1x=b1Knx,

    x(0,+), where Kn=n1i=0(αγ)i+αn1γn. Then

    Vn(x)=b1Knlnx+Cn

    for all n1, xX. Taking the limit when n,

    Vn(x)=b1Knlnx+Cnυ(x)=b1Klnx+C,

    xX, where

    K=limnKn=11αγ. (3.13)

    Therefore, for all xX,

    b1Klnx+C=maxa[0,xγ]{b1ln(a)b2+αE[υ((xa)γξ)]}=maxa[0,xγ]{b1[ln(a)+αγKln(xa)+αμK]+αCb2}. (3.14)

    The first-order condition of Eq. (3.14) is given by

    b1ab1αγKxa=0,

    and, using Eq. (3.13), this yields

    f(x)=(1αγ)x,

    xX. Thus, for xX,

    b1Klnx+C=b1(1+αγK)ln(x)+b1αγKln(αγ)+b1ln(1αγ)+b1αμK+αCb2.

    Then, K=1+αγK=11αγ and

    C=11α[b1αγ1αγln(αγ)+b1ln(1αγ)+b1αμ1αγb2].

    In the example presented above, with catch data from CONAPESCA [11] between 2006 and 2014 for the state of Puebla, Mexico, the production of mojarra is obtained for different values of the production factor γ. For this numerical case, the discount factor α=0.5 is used and π={a0,a1,a2,,a102}, where 102 is the number of months for which information is available between 2006 and 2014. Note that each ai, i=1,2,,102, corresponds to the weight of fish caught monthly, starting with a0=0 (for more details, see CONAPESCA [11]). In the first case, the production factor is γ=0.85, with an initial state x0=436,250, which is the initial amount of fish; the production is thus calculated for the first and subsequent months. Using the logarithmic utility function U(a)=b1ln(a)b2, aA(x) from Section 1, one can conclude that V102(x,π)=1,251,934.12. Analogously, when the production factor is greater than γ=0.90 but the initial stock is smaller, x0=220,000, the value function increases to V102(x,π)=1,252,077.98 with respect to the previous case, or γ=0.95 and x0=115,000, such that the value function is V102(x,π)=1,235,629.78. Table 2 summarizes the results.

    Table 2.  Utility estimations with discount factor α=0.5.
    γ x0 (kg of fish) V102(x,π) π
    0.85 436,250 1,251,934.12 {a0,a1,a2,,a102}
    0.90 220,000 1,252,077.98 {a0,a1,a2,,a102}
    0.95 115,000 1,235,629.78 {a0,a1,a2,,a102}

     | Show Table
    DownLoad: CSV

    The following tables are similar to Table 2, but with a different discount factor.

    According to Tables 2, 3, and 4, note that if α increases, then V102(x,π) increases. This means that as the discount factor increases, profits from fish sales will increase.

    Table 3.  Utility estimations with discount factor α=0.2.
    γ x0 (kg of fish) V102(x,π) π
    0.85 436,250 365,030.39 {a0,a1,a2,,a102}
    0.90 220,000 368,347.63 {a0,a1,a2,,a102}
    0.95 115,000 366,801.64 {a0,a1,a2,,a102}

     | Show Table
    DownLoad: CSV
    Table 4.  Utility estimations with discount factor α=0.9.
    γ x0 (kg of fish) V102(x,π) π
    0.85 436,250 7,616,999.13 {a0,a1,a2,,a102}
    0.90 220,000 8,857,481.60 {a0,a1,a2,,a102}
    0.95 115,000 12,433,010.56 {a0,a1,a2,,a102}

     | Show Table
    DownLoad: CSV
    Table 5.  Utility estimations with discount factor α=0.5.
    γ x0 (kg of fish) V102(x,π) π
    0.85 436,250 1,251,934.12 {a0,a1,a2,,a102}
    0.90 436,250 1,663,042.04 {a0,a1,a2,,a102}
    0.95 436,250 2,119,445.87 {a0,a1,a2,,a102}

     | Show Table
    DownLoad: CSV

    The following table shows a comparison between the value functions when the initial state is the same and the production factor is different.

    In this case, when the initial state is the same for each production factor, observe that V increases as the production factor increases. That is, if production increases, then profits will also increase.

    Table 6 presents the results for the optimal value function and optimal policy given in Theorem 3.8. In this case, the same initial values and production factors are used.

    Table 6.  Utility estimations for the optimal policy when α=0.5.
    γ x0 (kg of fish) υ(x) f(x)
    0.85 436,250 2,074,213.42 250,843.75
    0.90 220,000 1,876,290.76 121,000.00
    0.95 115,000 1,703,773.25 60,375.00

     | Show Table
    DownLoad: CSV

    Note that, in the three cases presented, υ evaluated in the initial state is always greater than the respective cases find in the last section. Furthermore, with this optimal policy, not higher profits are obtained, but also the fish population is in equilibrium. Additionally, a production factor equal to 0.85 increases profits by 66%, whereas production factors of 0.90 and 0.95 increase profits by 50% and 38%, respectively.

    Unexploited fish stocks are in equilibrium (zero population growth), the causes of increment only have to balance the losses produced by natural mortality. However, once over-fishing begins, more fish is then caught than the population can replace through natural reproduction, consequently, the population level is reduced. If this reduction is not drastic, less than the initial level, and if the exploitation remains stable and below f levels given in Tables 6 and 7, the most probable is that the fish stock increases or a decrease in natural mortality will allow the number of fish remain stable.

    Table 7.  Utility estimations for the optimal policy when α=0.2.
    γ x0 (kg of fish) υ(x) f(x)
    0.85 436,250 2,616,400.34 362,087.50
    0.90 220,000 2,232,771.20 180,400.00
    0.95 115,000 1,865,407.39 93,150.00

     | Show Table
    DownLoad: CSV

    Following the above, the conclusion is that the fish stock will remain in equilibrium as long as there is a reasonable fish exploitation, in this sense, it is proposed to follow the strategy f as shown in Tables 6 and 7. Thereby, it is guaranteed that both the population size and the catch amounts are kept at optimal levels, guaranteeing a sustainable performance. In other words, when f is applied, the fish stock and captures remain at a level that allows the population to fully develop its natural growth capacity. Additionally, from an economic point of view, the optimal profit is obtained. Conversely, when the fishing is greater than f, the stock can be reduced to lower levels than the initial state x0, which would reduce the natural capacity for population growth. In this case, there is a period of over-fishing that excessively reduces the stock and can lead to the collapse of the fishery.

    This work models the fisheries problem using CONAPESCA's mojarra extraction data from 2006 to 2014 for the state of Puebla, Mexico. Furthermore, concludes that the utility function is logarithmic. This fisheries model is studied as a MDP, allowing for the demonstration of the existence of an optimal policy. Using a dynamic programming equation and a version of the EE for the iterative value functions, the optimal value function and optimal policy are calculated for the case of an infinite horizon. This study shows how to obtain greater profits from the quantity of fish that can be extracted without reducing its population.

    Subsequently, particular values are assigned to the production factor to obtain numerical results for which the value function is calculated for the state of Puebla. The utility function is calculated with data available from CONAPESCA. An analysis is then presented that compares the optimal value function to increase utility, applying the optimal policy found in this work. The study concludes both theoretically and numerically that the optimal policy is indeed the best strategy to follow to maximize profits.

    The authors declare that they have no potential conflicts of interest.



    [1] D. Ngom, A. Iggidir, A. Guiro, A. Ouahbi, An observer for a nonlinear age-structured model of a harvested fish population, Math. Biosci. Eng., 5 (2008), 337–354. doi: 10.3934/mbe.2008.5.337
    [2] D. P. Bertsekas, Dynamic Programming and Stochastic Control, 1st edition, Academic Press, New York, 1976.
    [3] O. Hernández-Lerma, J. B. Lasserre, Discrete-Time Markov Control Processes: Basic Optimality Criteria, 1st edition, Springer-Verlag, New York, 1996.
    [4] M. L. Puterman, Markov Decision Processes: Discrete Stochastic Dynamic, 1st edition, Wiley-Interscience, California, 2005.
    [5] B. Liao, Q. Liu, X. Wang, K. Zhang, J. Zhang, K. H. Memon, M. A. Kalhoro, Asymptotic behaviour of the n-species stochastic Gilpin-Ayala Cooperative model, Stoch. Environ. Res. Risk Assess, 30 (2015), 39-45.
    [6] B. Liao, X. Shan, C. Zhou, Y. Han, Y. Chen, Q. Liu, A dynamic energy budget–integral projection model (DEB-IPM) to predict population-level dynamics based on individual data: A case study using the small and rapidly reproducing species engraulis japonicus, Mar. Freshwater Res., 71 (2019), 461–468.
    [7] D. Levhari, L. J. Mirman, The great fish war: An example using a dynamic Cournot-Nash solution, Bell J. Econ., 11 (1980), 322–334. doi: 10.2307/3003416
    [8] L. J. Mirman, Dynamic models of fishing: A heuristic approach, Control Theory Math. Econ., (1979), 39–73.
    [9] H. Cruz-Suárez, R. Montes-de-Oca, An envelope theorem and some applications to discounted Markov decision processes, Math. Methods Oper. Res., 67 (2008), 299–321. doi: 10.1007/s00186-007-0155-z
    [10] Y. K. Kwan, G. C. Chow, Estimating economic effects of political povements in China, J. Comp. Econ., 23 (1996), 192–208. doi: 10.1006/jcec.1996.0054
    [11] Información Estadística por Especie y Entidad. March 10th, 2020, de Gobierno de México Sitio web, CONAPESCA. (2006-2014). https://www.conapesca.gob.mx/wb/cona/informacion_estadistica_por_especie_y_entidad
    [12] H. Cruz-Suárez, R. Montes-de-Oca, G. Zacarías, A consumption-investment problem modelled as a discounted Markov decision process, Kybernetika (Prague), 67 (2011), 909–929.
    [13] A. de la Fuente, Mathematical Methods and Models for Economists, 1st edition, Cambridge University Press, 2000.
    [14] J. L. Troutman, Variational calculus and optimal control: optimization with elementary convexity, 2nd edition, Springer Science and Business Media, 2012.
    [15] H. Cruz-Suárez, R. Montes-de-Oca, Discounted Markov control processes induced by deterministic systems, Kybernetika, 42(6) (2006), 647–664.
    [16] H. Cruz-Suárez, G. Zacarías-Espinoza, V. Vázquez-Guevara, A version of the Euler equation in discounted Markov decision processes, J. Appl. Math., (2012), 1–16.
    [17] A. Jaskiewicz, A. S. Nowak, Discounted dynamic programming with unbounded returns: Application to economic models, J. Math. Anal. Appl., 378 (2011), 450–462. doi: 10.1016/j.jmaa.2010.08.073
  • Reader Comments
  • © 2021 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(2712) PDF downloads(151) Cited by(0)

Figures and Tables

Figures(2)  /  Tables(7)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog