About this Journal Submit a Manuscript Table of Contents
Discrete Dynamics in Nature and Society
Volume 2008 (2008), Article ID 793512, 21 pages
Research Article

Model-Matching-Based Control of the Beverton-Holt Equation in Ecology

Department of Electricity and Electronics, Faculty of Science and Technology, Institute of Research and Development of Processes, Campus of Leioa, Leioa Bilbao 48940, Spain

Received 24 July 2007; Revised 24 September 2007; Accepted 16 November 2007

Copyright © 2008 M. De La Sen and S. Alonso-Quesada. 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.


This paper discusses the generation of a carrying capacity of the environment so that the famous Beverton-Holt equation of Ecology has a prescribed solution. The way used to achieve the tracking objective is the design of a carrying capacity through a feedback law so that the prescribed reference sequence, which defines the suitable behavior, is achieved. The advantage that the inverse of the Beverton-Holt equation is a linear time-varying discrete dynamic system whose external input is the inverse of the environment carrying capacity is taken in mind. In the case when the intrinsic growth rate is not perfectly known, an adaptive law implying parametrical estimation is incorporated to the scheme so that the tracking property of the reference sequence becomes an asymptotic objective in the absence of additive disturbances. The main advantage of the proposal is that the population evolution might behave as a prescribed one either for all time or asymptotically, which defines the desired population evolution. The technique might be of interest in some industrial exploitation problems like, for instance, in aquaculture management.

1. Introduction

The nonautonomous discrete so-called Beverton-Holt equation (BHE) is very common in Ecology and, in particular, in studying the growth population dynamics (see, e.g., [19]). The equation is of great importance in the fishery industry concerning the growth and exploitation of species like, for instance, plaice, haddock and coho salmon, as well as other bottom feeding fish populations [1]. It is proven that, for subcritical biomass, according to a low value of the intrinsic growth rate of the studied population, zero is a stable attractor if the carrying capacity is not independent of the intrinsic growth rate in order that the model is well posed. In particular, the carrying capacity should decrease as the spawning stock decreases converging to zero. The BHE is a two-parameter widely used nonlinear equation of the form [19] ????+1=????????????????+????????-1??,?????0?0?:=???,(1.1)with ?? being the set of natural numbers, where ??0>0, ???????+0:=??+?{0} (??+ being the set of positive real numbers) is the intrinsic growth rate of the population, determined by life cycle and demographic properties (specie growth rate, survivorship rate, etc.), and ???????+, so-called the environment carrying capacity, is a characteristic of the habitat being dependent on resources availability, temperature, humidity, and so on. Typically, ????>1 and ????=????+?????+0, for some ??(=2)???, due to periodic fluctuations. It has been reported that a negative carrying capacity is nonsense since a small carrying capacity is sufficient to interpret a very unfavorable situation for the population [6], even to model the so-called Allee effect, that is, the difficulty of finding mates due to the very scarce number of population [7]. Note that ????=0 is sufficiently unfavorable to extinguish the population since, from (1.1), a recruitment ????+1=0 would occur for any spawning stock. The population extinction has been dealt with in [10] for intrinsic growth rates less than unity by discussing the spawning stock biomass conditions for the zero equilibrium to be a stable attractor. It has been concluded that the environment carrying capacity is not independent of the intrinsic growth rate for subcritical spawning stock biomass. Generalizations of the BHE including bounded discontinuities at sampling points, which can be interpreted as the consequence of impulsive effects [11], in the generating continuous logistic equation, have been also studied.

This paper is devoted to the model matching of a prescribed reference model which is also defined in practice by a BHE. The standard case of intrinsic growth-rate sequences being greater than unity is considered. The environment carrying capacity is locally modified around its reference values to achieve the prescribed behavior. Its inverse plays the role of the control of the IBHE, namely, the inverse of the BHE which is a linear dynamic system [9], and then much easier to deal with than the BHE. The control law consists of calculating the appropriate carrying capacity inverse so that the solution of the IBHE (namely, the inverse of the population levels) coincides with the solution of the reference IBHE. In that way, the population levels also coincide with the reference-suited ones. In the case when the intrinsic growth rate is not exactly known, the scheme is extended by incorporating parametrical estimation and a related adaptive control law. The tracking objective becomes asymptotic as a result of the incorporation of the estimation process. Since, in many practical situations, the reference model and the current one are locally deviated from each other, the solution matches very closely the reference solution for all time. The proposed technique is very feasible in close or semi-open environments like, for instance, aquaculture industry, for instance, by local modification of the temperature. In the case when the BHE parameters are not exactly known, an adaptive extension of the method incorporating parameter estimation is used instead. The paper is organized as follows. Section 2 is devoted to the statement of a control law for the case of known parameters. The BHE is generalized to include, if suited, additive disturbances which may be useful to describe uncertainties in the parameters of the BHE or phenomena which are not included in the nominal standard equation like, for instance, local or global migrations of the population. The reference model to be tracked (also often referred to as “matched” in the literature) is also a BHE which potentially includes additive disturbances corresponding to those present in the given BHE being defined as that possessing as its solution, the IBHE of the given BHE. Such an inverse has a linear structure by nature so that it is easy to deal with, and it describes the appropriate system behavior. The control law is of a feedback type and consists of designing appropriately the environment carrying capacity, within some prescribed margins, so that the reference model is perfectly tracked (i.e., the solution of the controlled equation and that of the reference model coincide for all samples). This is very feasible in closed environments like, for instance, fisheries. It turns out that the achievement of a perfect tracking of the IBHE corresponds to achieving the objective of a perfect tracking of the corresponding reference BHE. Section 3 extends the method to the adaptive case by incorporating an adaptive scheme. It is assumed that certain a priori knowledge on the additive disturbances is available so that a relative dead zone is used in the adaptation algorithm to freeze the adaptation under small-adaptation errors to prevent against stability caused by the disturbances when the parameters are not fully known [12]. The adaptive control law is implemented by replacing the controller parameters in the case of perfect knowledge by those computed based on the estimates of the IBHE. The closed-loop system is proven to be stable if the additive disturbances grow, at most linearly, with the maximum of the solution of the IBHE. In the adaptive case, the tracking objective is only achievable asymptotically provided that the contribution of the additive disturbances is zero or becomes extinguished asymptotically. The method to discuss the stability and convergence of estimates relies on the use of a Lyapunov sequence, [8, 13, 14]. Section 4 presents two numerical examples and, finally, conclusions end the paper.

2. Control Law for the Case of Known Parameters

Define an inverse system of (1.1) through the change of variable ????=????-1 resulting directly from (1.1) in the time-varying discrete linear system ????+1=????????+????????,??0=??0-1>0,(2.1)where ????=????-1, ????=1-????-1=1-????, and ????=????-1, that is, the inverse carrying capacity is the control action. If there exists an additive disturbance sequence ???? in (2.1), one gets, instead of (2.1) [9], ????+1=????????+????????+????.(2.2)The IBHE (2.2) corresponds, through the change of variable ????=????-1, to an extended BHE which includes a sequence of additive disturbances: ????+1=????????????????+?????+?????????????=????????????????+?????-1+?????????????????=????????????????+????????-1??+??0??=????????????????+????????-1??+????,(2.3)where ??0??=????????????????,??????=-????????????0???????+????????-1????????+?????-1+????????????????????=-2????2????2???????????+????????-1????????+?????-1+??????????????????.(2.4)Note that the sequence {????}80 is an additive disturbance of the BHE which is identically zero if the additive disturbance sequence {????}80 of its inverse is identically zero. The disturbance may include the effects in the solution of parametrical uncertainties, for instance, in the intrinsic growth rate, or effects like, for instance, migrations or local migrations which are not taken into account in the standard BHE. The following assumptions are concerned with the stability and controllability of the IBHE and with the knowledge of an upper bound of the disturbance term if such a disturbance is unknown.

Assumption 2.1. 8>????=1+???? and 8>????=????, for all ?????0, and some ???????+,???????+.

Assumption 2.2. ????(???+0)=????, for all ?????0, and {????}80 is known.

Note that Assumption 2.1 implies that the BHE is stable since ????-1<1, and also controllable since ????=?????:=????/(1+????)>0 so that (1-????-1), for all ?????0, so that ???? in (2.2), and then the BHE solution ????=????-1 may be driven to sequences of prescribed values. These two properties are crucial for the solution to track a prescribed reference sequence. Note also that ???? may be chosen in the real interval (0,1) without any loss in generality since, if Assumption 2.1 holds for any given ????>0, then it always holds for some ?????(0,1). This fact is used in the proof of the subsequent result, which is direct and related to the positivity and boundedness of the solution.

Assertion 2.3. If Assumptions 2.1 and 2.2 hold and {????}80 is uniformly bounded, then all solutions of the BHE and the IBHE are uniformly bounded and positive, provided that 8>??0>0.

Proof. ????=1+??????????[?????,1-?????]. Define ??:=Max0=??<8(????-1)?(0,(1/(1+?????))]?(0,1). Since ??0=??0-1>0, it follows from (2.2) that 0<????-1=????=??????0+??-1???=0????-??-1?1-???????????-1+?????=??????0+1-?????1-??1-??????Max0=??<8?????-1+?????=??0+?1-?????2??????????Max0=??<8??????<8????????0,8,??????01,since=1-??1-??????????,Max0=??<8?????-1?=1????.(2.5)

The subsequent set of structures is considered in order to be then able to formulate the control law.

(1) A Reference
BHE ??*??+1=??*????*????*??/(??*??+(??*??-1)??*??), for all ?????0, which defines the suitable solution through the appropriate reference values of the intrinsic growth rate and the environment carrying capacity of its corresponding reference IBHE ??*??+1=??*????*??+??*??????, with reference input ????=(??*??)-1>0, and parameter sequences ??*??=(??*??)-1<1 and ??*??=(??*??)-1>0.

(2) The Current
BHE ????+1=????????????/(????+(????-1)????)+???? and its corresponding IBHE ????+1=????????+????????+????,(2.6) whose parameters are ????=????-1, ????=1-????-1 are deviated from the suited ones ??*??=(??*??)-1 and ??*??=1-(??*??)-1, which may include disturbance terms grouped in a sequence {????}80 with control input ????=????-1 obtained according to a feedback law.

(3) The Control Law:
????=????????+????????-????,(2.7) provided that ???? is known, which is a relaxation of Assumption 2.2 which is parameterized by ????=??*??-????????,????=??*??????,????=????????,(2.8)with ????=1-????=1-????-1. Note that the control law generates the inverse carrying capacity from the reference inverse carrying capacity and the solution of the IBHE with two parametrical sequences {????}80 and {????}80 and a correcting sequence {????}80 calculated from (2.8). The subsequent constraint ????=???*???-1>????????-????????????=1?????????1-????+????-??*???????(2.9)must hold for all samples to guarantee that the environment carrying capacity is positive. The following result establishes that the tracking error in between the reference sequence and the current solution for the given control law is zero. In that way, the perfect tracking objective is achieved if the disturbance sequence is known.

Proposition 2.4. If ??0=??*0, then the control law (2.7)-(2.8) achieves an identically zero tracking error {????:=????-??*??}80.

Proof. It follows, by complete induction, by assuming ????=??*???????=0, for all ?????0, that ????+1=?????+????????-??*?????*??+?????????-??*???????=0,??????0.(2.10)

If the disturbance is unknown, but Assumption 2.2 still holds, then the last identity in the parameterization of the control law in (2.8) is replaced, with similar expressions, with the disturbances absolute upper bounds, namely, ????=????/????, where |????|=????. Thus the tracking error is not identically zero but it is uniformly bounded. It is also very close to zero, provided that the deviation of the upper bound of the disturbance with respect to such a disturbance is very close to zero. Then one has the subsequent result.

Proposition 2.5. If Assumption 2.1 holds and |????|=???? with ???? being unknown, but ???? being known, then the control law (2.7)-(2.8) of Proposition 2.4, with the modification ????=????/????, achieves a uniformly bounded tracking error {????:=????-??*??}80, provided that the disturbance sequence is upper bounded by a bounded sequence.

Proof. Direct calculations yield ????+1=????????+????????+????,??*??+1=??*????*??+??*??????(2.11)so that, from (2.7)-(2.8), ????+1=??*??????+?????-?????=???*???+1??0+????1-??*<8,??????0,(2.12)where 0<??*??=??*<1, for all ?????0, and ????:=Max0=??<8?|||????-????|||?=2Max0=??<8??????<8.(2.13)

Since the inverse of the environment carrying capacity is the control action, a large deviation from its nominal values for tracking purposes may be not admissible. Note that Assumption 2.1 also applies to the reference to be tracked by the control action consisting of appropriate achievement of an environment carrying capacity, that is, 8>??*??=1+????>1, and 8>??*??=????>0, for all ?????0. That means that the maximum deviation of the suited reference sequence, related to the current solution, depends on the maximum allowed variations of the environment carrying capacity with respect to its nominal value. More precisely, the subsequent result follows.

Proposition 2.6. Assume that the maximum allowed absolute variation of the carrying capacity with respect to its nominal value is ???>0. Then the following properties hold.
(i) The control parameter sequence {????}80 is subject to the constraint ?????[1-???/??*??,1+???/??*??], for all ?????0, which induces the parametrical constraints ???????[*????*??-?????*??-???,??*????*??+?????*??+???],???????[*????*????*??+???,??*????*????*??-???].(2.14)(ii) The tracking-error sequence {????:=????-??*??}80 is subject to the constraint (????+1-????????+???*??-???????*??)?[????-??*???????*??2,????+??*???????*??2],(2.15) for all ?????0, which is trivially bounded, provided that the disturbance sequence is bounded.(iii) A positive carrying capacity is obtained from the controller synthesis if (??*??)-1>(????????-????)/????????=(1/????)(????/(1-????)+(????-??*??)/????), when ???? is unknown but ????=|????| is known. The control parameter ?????(0,1] and the calculated ???? from the control law are positive if ???? is calculated satisfying ????[0,??*??-????/(1-????)+(??*??-????)/????) with ?????[0,(1-????)(??*??+(????-??*??)/????)).

Proof. (i) Note that ??????[1-????*???,1+????*???]?1-????????[*???1-??*?????*??+???,??*???1-??*?????*??-???]???1-*???1-??*?????*??+???=??????=1-*???1-??*?????*??-???,(2.16)which yields directly Property (i) since ????=1-????.
(ii)One gets from (2.6)–(2.8) and Property (i) that ????+1=????+1-??*??+1=????????+?????-??*?????*??+????????-1*??????+?????[????????+?????-??*?????*??-?????*????*??????+????,????????+?????-??*?????*??+?????*????*??????+????],(2.17)and the result follows.(iii) The first part follows with the replacement of?????????=|????| in (2.9). The second part follows directly by replacing the admissibility domain of the controller parameter ???? in the above constraint.

3. Adaptive Control

For the case when the parameters of the BHE are unknown, an estimation scheme with adaptation dead zone for robust closed-loop stabilization is incorporated (see, e.g., [12, 13, 15, 16]). Such a mechanism governs the evolution of the IBHE and it is now discussed from a theoretical point of view. It will then corroborate, through numerical simulation, the potential usefulness of numerical and adaptive techniques as previously tested in some ecological models (see, e.g., [1719]). One takes the advantage that, since the parameters fulfill the constraint ????+????=1, for all ?????0, the same constraint might be introduced in the estimation scheme, what reduces it to the estimation of one parameter only. A modification would consist of estimating both sequences of parameters separately. In that case, the estimation is performed without taking into account the parametrical constraint at the expense of an increase in computation time and memory storage. However, it may be proven that all the properties of boundedness and convergence of the estimates still hold if existing parametrical constraints are not taken into account, provided that each parameter is updated with its corresponding component in the set of measured data (see, e.g., [12, 13, 15]). In a first step, a simple procedure is developed for the nominal BHE through controlling its inverse, and then an extension is made to prevent a good operation in the presence of additive disturbances. A relative dead zone to the estimation scheme is added to prevent potential instability caused by those disturbances.

3.1. Adaptive Control approach Based on the Intrinsic Growth Rate Estimate

Additive Disturbances
The deviations of the intrinsic growth rate with respect to a certain unknown constant value ??=??-1, implying ??=1-??, are incorporated within a disturbance sequence {????}80 which can also include other unstructured disturbance contributions in a sequence {??0??}80: ???????:=?????-????+????????-????+??0??=???????-??????-?????+??0??,(3.1)and ????=????-1. In this way, the nominal parameter vector of the IBHE, namely, ??=(??,??)??, is constant, which then facilitates an updating algorithm that has an identical structure for the parameter estimates as for the parametrical errors. The IBHE, its estimate, the identification error, the reference IBHE, the tracking error, and the parameter-adaptive law are, respectively, as follows.

????+1=??????+??????+???????=????-?????+????+????.(3.2)If ????=0, the resulting particular case of (3.2) is called the nominal IBHE, and its corresponding inverse equation, whose solution is the sequence {????=????-1}80, is called the nominal BHE.

Estimate of the IBHE
????+1=????????+????????=?????????-?????+????,(3.3)where ???? and ????=1-???? are the estimates of ?? and ?? at the ??th sample.

Identification Error
????+1:=????+1-????+1=-?????????-?????????+????=-??????????-?????+????,(3.4) where ?????:=????-??, ??????:=??-??=??-????=-????? are the parametrical errors. Note that although ???? is unknown, the identification error is available for measurement through the first identity in (3.4).

Reference IBHE
It is assumed to be time invariant as the nominal IBHE, namely, ??*??+1=??*??*??+??*????=??*???*??-?????+????(3.5)with ??*=1-??* and ????=(??*??)-1, is the reference environment carrying capacity inverse.

Tracking Error

Adaptation Algorithm 1:
????+1=????+?????????-?????????+11+?????????-?????2,??0???0,1,??????0,(3.7)is a recursive estimation of the intrinsic growth rate, where {????}80 and {????}80 are real sequences satisfying ????=0, ????>0 to be specified later for purposes of convergence of the algorithm. The control law has the same structure as (2.7)–(2.9) by replacing the true parameters by its estimates, and by deleting the correcting sequence ????. Instability drawbacks due to disturbances will be dealt with by using a specific dead-zone updating rule for the sequence ????=0 in the parameter-adaptive algorithm.

Adaptive Control Law:
????=????????+????????=11-??????1-??*???????+???*??-??????????(3.8)with ????=(??*??)-1 with the sample-dependent controller parameters ????=??*??-????1-????,????=1-??*??1-????.(3.9)The following result concerning useful technical relationships in between the IBHE solution, its estimated IBHE, and the reference IBHE independent of the parameter adaptive algorithm is immediate.

Proposition 3.1. The adaptive control law (3.8)-(3.9) leads to ????=??*?? (i.e., the estimated solution equalizes the reference solution), ????=-???? (i.e., the tracking error equalizes the minus-identification error), for all ?????, for any finite ??0, ??0=??*0. The above results are independent of the parameter-adaptive algorithm.

Proof. Proceed, by induction, by assuming ????=??*?? for any given ?????0. Then one gets, from (3.8)-(3.9) into (3.3) by using ????=1-????, ????+1=????????+?1-??????????????+?????????=??????1-?????+?????????+?1-?????????????=??*??????+??*??????=??*??+1(3.10)if ????=??*??, and then using (3.6), ????+1=???*??+1-????+1?+?????+1-????+1?=?????+1-????+1?=-????+1,??????0,(3.11)what completes the proof.

3.2. Boundedness and Convergence Results of the Adaptation Algorithm

In the following, the boundedness and convergence properties of the adaptation algorithm are investigated in the subsequent result.

Theorem 3.2. Assume that a sequence {????}80 satisfying ????=|????|, for all ?????0, is known that the adaptation algorithm (3.7) is updated by choosing ????=??0, and ????=???????|||??0??????+1|||=??????,2???????-1-??1???|||????????+1|||>??????,(3.12)for all ?????0, for some prefixed real constants ??>1, ??-1>??0>0, and ??1?[??0,??-1). Then the following properties hold irrespective of the control law.
(i) The sequences {????}80, {(????/(1+????(????-????)2))1/2|????+1|}80, and {(????/(1+????(????-????)2))??2??+1}80 are bounded and have limits ??8<8, and zero, respectively.(ii)There exist the subsequent limits lim???8(?????=0??????2??1+?????????-?????2)=lim???8(?????=0??????2??1+?????????-?????2)=lim???8(?????=0??????2??+11+?????????-?????2)<8,lim???8(?????=0(??????2??1+?????????-?????2)1/2)=lim???8(?????=0(??????2??1+?????????-?????2)1/2)=lim???8??((????2??+11+?????????-?????2)1/2)<8.(3.13)

Proof. Let {????}80 be defined by ????:=???2??. Note, by subtracting ?? in both sides of (3.7), that ?????+1=?????+????????????+1/(1+??????2??), where ????:=????-????. Since ?????????=????-????+1=??????-????+1,(3.14)from (3.4), for any absolute upper-bound ????=|????|, one gets ?????:=???2??+1-???2??=??2????2????2??+1?1+??????2???2+2????|||????+1|||?????-|||????+1|||?1+??????2??=????|||????+1|||1+??????2??(??????2??|||????+1|||1+??????2???|||??-2??+1|||-?????).(3.15)From (3.15), ????=0??????=0. Assume |????+1|>?????? so that for any real ??>1, one gets, from (3.15), ?????=??????2??+11+??????2??(??????2??1+??????2???-21-??-1?)=-2??????1??2??+1???1+??????2???=0(3.16)for all ?????0, provided that ????>0 and ???? satisfies 0=??????2??/(1+??????2??)=????/????=2(??-1-??1)/??0 for any real ??1?[??0,??-1). Then ????=??0<8, what implies that the sequences {|?????|}80, {?????}80 and {????}80, and ??{(??1+?????????-?????2)1/2|||????+1|||}80??,{(??1+?????????-?????2)1/2|||????|||}80(3.17)are uniformly bounded and there exist limits lim???8(|?????|)=lim???8(??????2??+1/(1+????(????-????)2))=0, lim???8(????)=??8=??0<8, lim???8(????)=??8<8, and lim???8(?????)=???8<8. lim???8(??????2??+1/(1+????(????-????)2))=0?lim???8((????/(1+????(????-????)2))1/2|????+1|)=0 so that lim???8((????/(1+????(????-????)2))1/2|????|)=0 since either ????=2??0(??-1-??1)/?? or ????=0, for all ?????0, and Property (i) is proven. Since ???? converges to a finite limit, Property (ii) follows from ?????=0??????2??1+?????????-?????2=?????=0??????2??+11+?????????-?????2=?????=0||?????||=??0-????+1=??0-??8<8(3.18)as well as the boundedness of the square of each element of the sequence associate with the above series.

The boundedness of the estimates and estimation and tracking errors as well as the convergence of the estimates to finite limits are crucial issues to formulate a well-posed problem. The relative dead zone of the algorithm (3.12) is relevant for adaptive stabilization under external disturbances. The interpretation is direct; if the identification error is small related to the known upper bound of the disturbances, the estimation is frozen (i.e., stopped). The reason is that the contribution to the identification error of the parameterized part, which is being updated, may be smaller than that of the disturbance and then to maintain the estimation, updating may lead to instability. The incorporation of relative adaptation dead zones into the estimation scheme is a standard mechanism to cope against adaptive instability since the sequence {????}80 is bounded, see, for instance, [12, 13, 15]. Otherwise, the estimates computations and then the control law are unfeasible. If Assumption 2.1 is satisfied, then a projection method may be used for the estimates to guarantee that they lie inside its definition domain. Theorem 3.2 is then extended as follows.

Theorem 3.3. If both the nominal and current intrinsic growth rates satisfy Assumption 2.1, then there exist real constants ???????+0, ??1>1, such that ??,??????11+??1????,11+?????,??????0,(3.19)and the above adaptation algorithm is replaced with the set of modified recursive equations (3.20)–(3.22) below as to include projection of the estimates on an admissibility domain.
Adaptation Algorithm 2
It consists of two steps, namely, a priori and a posteriori estimations as follows.

A Priori Estimation:
??0??+1=????+1+?????????-?????????+11+?????????-?????2,??0???0,1,??????0,(3.20)being updated by choosing ????=??0, and ????=???????|||??0if??+1|||=????0??,2???????-1-??1???|||??if??+1|||>????0???:??????+?1-??????0?????0???????-????,??(3.21)for all ?????0, for some prefixed real constants ??>1, ??0>0, and ??1?[??0,??-1).

A Posteriori Estimation Via Projection on the Interval (1/(1+??1????),1/(1+????)):
????+1=?????????11+??1????if??0??+1=11+??1????,??0??+1if??0??+1??11+??1????,11+?????,11+????if??0??+1=11+????.(3.22)Then Theorem 3.2(i)-(ii) holds. In addition, similar properties hold by replacing ???????0??=????+(1-??????0??)??0???????-?????,???????0??=????+(1-??????0??)??0???????-?????,????+1??Aposterioriidenti?cationerror:=????+1-????+1=????+1-?????????-????????0??+1??Aprioriidenti?cationerror:=????+1-??0??+1=????+1-??0???????-?????.(3.23)

Proof. Define the sequence 0=???0??:=(???0??+1)2-???2??=????|||????+1|||1+??????2??(??????2??|||????+1|||1+??????2??|||??-2(??+1|||-??0??))=0(3.24)from (3.20)-(3.21) using a close development to that in Theorem 3.2 so that it is monotonically decreasing ?????:=???2??+1-???2??=???2??+1-????0??+1?2-|||???0??|||=-(1-(????+1??0??+1)2)????0??+1?2-|||???0??|||=0(3.25)for all ?????0 since 0<????+1/??0??+1=1, for all ?????0, from the projection part (3.22) of the estimation algorithm. Then {????}80 is bounded and converges to a finite limit, and {?????}80 is also bounded and converges to zero. Also |???0??|=(1-????+1/??0??+1)<1 for all ?????0, and |???0??|=(1-????+1/??0??+1)<1. Since ???0??:=(???0??+1)2-???2??=0, for all ?????0, then there exists a finite limit lim???8(?????=0|???0??|)<8, then ((???0??+1)2-???2??) converges to a finite limit and |???0??| converges since |?????| converges (since ???? converges) to finite limits as ???8. Thus Theorem 3.2(i)-(ii) follows. On the other hand, note that ????+1??Aposterioriidenti?cationerror:=????+1-????+1=????-??????????-?????=??0??+1+???0??-??????????-?????=??0??-???0???????-?????,??0??+1??Aprioriidenti?cationerror:=????+1-??0??+1=????-???0???????-?????=????+1+?????-??0????????-?????=??0??+?1-??????0????2??0????-??????-?????.(3.26)Since the normalized sequences {(????/(1+??????2??))1/2|????+1|}80 and {(????/(1+??????2??))1/2??0??}80 are bounded, it follows that the reported replacements may be performed leading to Theorem 3.3(iii).

3.3. Some Practical Design Constraints

The tracking error of the solutions of the IBHE and the reference IBHE are given by the set of equivalent expressions: ????+1=??*??+1-????+1=??*??+1-????+1-????+1=??*??+1-????+1+??????????-?????-????.(3.27)The solution of the IBHE is equivalently calculated either as a function of the nominal parameters ?? and ??=1-?? or as a function of their estimates as follows, by considering the adaptive control law (3.8): ????+1???=????-?????+????+????=????+1+????+1=?????????-?????+????+????-??????????-?????(3.28a)=????+1+????+1=?????????-?????+????+????-??????????-?????(3.28b)=????????+??1-??*???????+???*??-??????????+????+1(3.28c)=??*??????+?1-??*???????+????+1.(3.28d)It is important to discuss when the adaptive control law has no division by zero with the property being numerically robust. This translates, in practice, to ?????1 being sufficiently deviated from unity and {????}80 with all its elements positive and sufficiently deviated from zero. This guarantees that the BHE does not diverge when driven by an appropriate nonzero finitely upper-bounded environment carrying capacity. More formally, a reasonable assumption for practical design purposes is to extend the assumption of Theorem 3.3 to also consider the same admissibility range referred to for the reference intrinsic growth-rate inverse and the estimated intrinsic growth rate.

Assumption 3.4. ??,????,????,??*???[1/(1+??1????),1/(1+????)], for all ?????0, and some prefixed real constants ????>0 and ??1>1. A direct calculation with Assumption 3.4 yields ?1-????????[??1+????,??1????1+??1??????1]?0,1,(1-????1)?[1+??1????1,1+????|||??],*??-????|||=???????1?-1?1+??1??????1+?????,???*??-????????=[-1???-1???1+??1??????1+?????,???1???-1???1+??1??????1+?????],(3.29)while similar expressions follow for replacements involving any two elements of the quadruplet (??,????,????,??*??) as follows: ?1+????1+??1?????[??1????+??1-11+????????]=????=11-????[?1-??*???????+?????-??*???????]=?1+????1+??1?????[??1?????-11+?????????+11+?????????]=1+????1+??1????????1?????-?????+11+?????????.(3.30)The above expressions (3.29)-(3.30) lead to the following result referring to explicit positive lower and upper bounds of ????=????-1=(??*??+?????)-1 so that ????-1=??*??+?????=??*??+????-1 since ????=(??*??)-1.

Proposition 3.5. Assume ????<??1(1+????)/(??1-1)??*??, for all ?????0. Then the following constraints hold: 1+??????1?1+????????*???-1+???1???-1??=1????=1+??????1?1+?????????*???-1-?????+????,(3.31)1+??????1?1+????????*???-1+???1???-1??-??*??=?????=1+??????1?1+????????*???-1-???1???-1??-??*??.(3.32)

Proof. It follows directly from (3.29)-(3.30). The necessary condition ????<??1(1+????)/(??1-1)??*?? is needed to keep (??1(1+????)(??*??)-1-(??1-1)????)>0 in the right-hand-side upper bound of ?????, for all ?????0 in (3.32).

3.4. Closed-Loop Stability

The boundedness and convergence properties of the parameter estimates have been proven. In the following, the closed-loop stability is proven under a condition of slow growing of the disturbances with respect to the solution of the IBHE if such a disturbance is unknown. In particular, it is assumed that ????=??[Max0=??=??(????)] with known constants of the related linear upper bound and with sufficiently slow growth of ???? with respect to Max0=??=??(????). This is a reasonable assumption used in the background literature since a complete lack of knowledge of disturbances makes impossible the stabilization in the general case (see, e.g., [12, 13, 15] and references therein).

Assumption 3.6. There exist known finite nonnegative real constants ??1,2 such that ????=??1+??2Max0=??=????????.(3.33)

Then, from (3.28b), (3.31)-(3.32), and Assumptions 3.63.4, one gets ????+1=?????????-?????+????+????+1.(3.34)The following results are preparatory for the stability theorem to be stated later on. They basically establish that the absolute value of the identification error grows not faster than linearly with the upper-bounding sequence {????}80 of the additive disturbances.

Assertion 3.7. If the Adaptation Algorithm 1 is used, then |????+1|=??1????+??2, for all ?????0, for some real constants ??1(=??)???+ and ??2???+0.

Proof. If |????+1|=??????, the assertion is true with ??1=?????+ and ??2=0. If |????+1|>??????, then, from Theorem 3.2, (i) 0<(|????+1|-??????)/(1+????(????-????)2)1/2=??<8, for all ?????0. Thus ||????+1||=???????v+??1+????||????-????||?=?v??+?????????v+????||????||+??(3.35)and the assertion follows with ??1v=??+????,??2?v=??1+???Max0=??<8?||????||?<8(3.36)since Min0=??<8(????)>0.

Assertion 3.8. If the Adaptation Algorithm 2 is used and Assumption 3.4 holds, then |????+1|=??1????+??2, for all ?????0, with ??12???=1+1???-1???1+??1??????1+?????,??2=Max0=??<8?|||?????|||?Max0=??<8?|||????|||?.(3.37)

Proof. It follows from Assumption 3.4, (3.4) and (3.29), and the fact that {????}80 is positively lower bounded, which yield directly ||????+1||=||??????????-?????||+????=(1+Max0=??<8?||?????||?)????+Max0=??<8?||?????||?Max0=??<8?||????||?=??1????+??2,??????0.(3.38)

From Assertions 2.3 and 3.7, for the Adaptation Algorithms 1 and 2, respectively, Assumption 3.6 and (3.28c), ????+1=|||??*??????|||+|||1-??*??|||????+??1????+??2=??????0+1-???????1-??1??1Max0=??=??+1??????+??1??2+??1??2+??2?=??1??11-??Max0=??=??+1??????+1???1-??1??2+??1??2+??2?+??0,??????0,(3.39)for any ??0>0 and some positive finite constants ??,??1, and ??2 satisfying 1>??=Max0=??<8(??*??), ??1=1-Min0=??<8(??*??)<1, and ??2=Max0=??<8(????). The last two inequalities of (3.39) follow directly by proceeding recursively with the preceding one by using ???-1??=0????=(1-????)/(1-??)=1/(1-??) since ???(0,1). It is obvious that, for each ?????0, it exists a unique ???(??)???0, dependent in general on ??, such that ???(??)=?? such that ????'=Max0=??=??'(????)=Max0=??=??(????) and there is no other ?????0 such that ??=??>???(??), where the same maximum is reached. A precise definition is ???????:=Arg(?????0???:0=??=????'=Max0=??=??'??????????=Max0=??=????????,???¬?????????<??=?????0:????=Max0=??=????????)(3.40)with “?” and “¬” being the conjunction and negation logic symbols. The substitution of (3.40) into (3.39) yields Max0=??=????????=??1-??-1??<8,??????0,(3.41)provided that ??1??1/(1-??)<1, which holds if and only if 0=??1<(1-??)/??1 (i.e., if Assumption 3.6 holds with sufficiently small ??1), where ????:=1??111-????:=???1-??1??2+??1??2+??2?+??0(3.42)with ?? being a finite positive real, which depends on the initial conditions of the IBHE (and thus on the initial conditions of the BHE) on the available upper bound of the additive disturbances, on the reference, and on the free-design parameters of the adaptation algorithm, ????, ??>1, ??0, ??1, and so on. The above result also holds for the Adaptation Algorithm 2 with specific constants redefined by 1??:=1+????<1,???1??:=1????1+??1????,??1??:=?11+??1????=??1?????1+??1?????2,??2:=??1?1+?????Max0=??<8???*???-1<8,(3.43)and ??1,2 defined in Assertion 2.3 for the Adaptation Algorithm 1 and in Assertion 3.7 for the Adaptation Algorithm 2. The boundedness of the solution of the IBHE for finite initial conditions implies that of the additive disturbances via Assumption 3.6, that of the identification error via Assertions 2.33.7, and that of the tracking error since the solutions of the current and the reference one are both bounded. The estimated parameters and the estimated solution are bounded from Theorems 3.2 and 3.3. Those conclusions are now summarized in the result below.

Theorem 3.9 (closed-loop stability). If Assumption 3.6 holds with sufficiently small coefficient ??1, then all the signals in the closed-loop system remain bounded for all time so that the IBHE (and then the associate BHE) is globally stable if any of the Adaptation Algorithms 1 or 2 is used.

Remark 3.10. Theorem 3.9 becomes stronger than formerly stated in the absence of additive disturbances (i.e., ????=0) since asymptotic perfect tracking is furthermore achieved. Another result is that the identification error converges asymptotically to zero if the relative dead zone of the adaptation algorithms is removed, that is, the {????}80 is not zeroed at any sample as a result of a comparison procedure of the identification error and the additive disturbances. These results follow directly from Theorems 3.2 and 3.3 and (3.39). Note that the results may be directly generalized to the extension of using relative dead zones for adaptive algorithms by incorporating the estimation of the constants ??1,2 characterizing {????}80, provided they are unknown [16].

Remark 3.11. The closeness between the IBHE and that of its reference model has been considered for feasibility reasons since, in many ecological problems, where such models are commonly of interest, the environment characteristics, which directly influence the value of the environment carrying capacity, cannot be abruptly modified even in closed environments. The employed philosophy about relative dead-zone adaptation seems to be promising to be also applied to other ecological controlled problems involving, for instance, Kolmogorov-type ecological models or models for a biochemical aquariums [17, 18]. The obtained results seem to be promising for mathematical modeling and the use of adaptive control in other kinds of problems related to biology and ecology (see, e.g., [2022]) and to extend the research to governing extinction conditions [10].

4. Numerical Examples

Example 4.1. The objective is that the BHE solution {????}80 tracks a reference sequence {??*??}80 by synthesizing a control system in order that IBHE matches a reference model. Such a reference model is another IBHE defined by the appropriated sequences {??*??}80 and {??*??}80 to generate the sequence {??*??}80 associated to the desired {??*??}80 to be tracked by the BHE solution. The reference values are ??*??=1.55, and ??*??=170, for all ?????0, are used for the reference model and BHE, respectively, with the period ??=1 (one year). It is assumed that small local variations of the carrying capacity of the habitat are allowed. The evolution data of the population of cod in the North Atlantic Sea in the period 1952–1992 are given in [4], and exhibit a quasiperiodic behavior of the intrinsic growth rate. Based on these data, the intrinsic growth rate is modeled as a sinusoid defined by ????=2.3+0.15sin(????/100), which adjusts closely to those experimental data. The average value of ???? is 200, what is used as a fixed uncontrolled carrying capacity for the whole set of data. The initial population of the cod is ??0=500 being identical to the reference value. Figure 1 displays the population of cod if the BHE is controlled and if it is not controlled.
The population of the controlled BHE tracks perfectly the reference sequence as in the transient as in the stationary regime for all time since the true parameterization of the BHE is known, and then the control law is not adaptive. Note, however, the uncontrolled BHE does not track the reference. Figure 2 displays the controlled carrying-capacity sequence, which leads to a perfect tracking of the reference.

Figure 1: Evolution of the population of cod from ??0=500.
Figure 2: Evolution of the control sequence for the solution of Figure 1.

Example 4.2. A numerical adaptive example is now discussed. It is assumed that the environment carrying capacity of the reference model-environment carrying capacity is a constant value of 200, and the reference intrinsic growth rate is constant of value 2. The objective is to design the environment carrying-capacity sequence via estimation and feedback so that asymptotic tracking of the reference solution is achieved with positive solutions. It is assumed that no unmodeled dynamics or parametrical uncertainties are present so that a standard recursive least-squares algorithm is used for parametrical estimation of the intrinsic growth rate whose initial value is chosen and the initial value of the time-decreasing covariance gain is fixed to 1000. Figure 3 displays the reference BHE, the uncontrolled BHE, and the controlled one. It is seen that the reference BHE is asymptotically tracked, while the uncontrolled one differs from the reference. The carrying capacity of the unknown uncontrolled BHE varies from 180 to 220 according to the shape of Figure 4, and its intrinsic growth rate is constant, equal to 1.6. The carrying capacity exerting the control action is also displayed on Figure 4 with maximum and minimum admissible values being constrained by the adaptive controller to 280 and 200, respectively.

Figure 3: BHE Solutions: (a) reference, (b) uncontrolled, and (c) controlled for initial covariance gain 1000.
Figure 4: Environment carrying capacities: (a) uncontrolled BHE and (b) controlled.

5. Conclusions

This paper has considered the well-known BHE used in ecology within a control problem context where the control action on the linear IBHE is the inverse of the carrying capacity. An extended version of the standard BHE has been considered by incorporating additive disturbances. The overall control problem is firstly stated on the IBHE by taking advantage of its linear nature. This point of view is very feasible in close or semi-open environments where humidity, temperature, and other factors of the environment may be selected within certain margins. A reference is also defined, which describes the suited behavior for the system, the control action having an objective that the solution of the current BHE is able to perfectly track that of the reference (a). For feasibility purposes, the overall problem is stated in terms of a local variation of the inverse carrying capacity to perform the control action on the IBHE so that the reference (a) is fully tracked. This implicitly means that, in practice, the current parameters of the BHE and then those of the IBHE are locally deviated from those of its reference model. Then the method has been extended by incorporating an adaptive version for the case when the BHE parameterization is partly or fully unknown. The use of a relative adaptation dead zone freezes the adaptation when the identification error is sufficiently small according to an available absolute upper bound of the disturbances sequence so that the algorithm is proven to prevent potential instability caused by the presence of such additive disturbances. The second one incorporates to the estimation dead zone and estimates projection procedure by using a priori knowledge on the parameters of the BHE and the second algorithm, which have been presented. It has been proven that both adaptation algorithms and associate control law stabilize the current, provided that the absolute value of the additive disturbances grows not faster than linearly with the maximum of the solution with a sufficiently small slope. Some related numerical examples have been discussed to corroborate the theoretical results. The potential generalization to the general case, where the parameterizations of the current BHE and its reference one to be tracked are not close to each other, is direct although it may be not feasible in some practical cases, since the adaptation algorithms and the main stability results are formulated in a general way.


The authors are very grateful to the Spanish Ministry of Education for its partial support of this work through Project DPI 2006-00174. They are also grateful to the referees for their useful comments.


  1. N. J. Barrowman, R. A. Myers, R. Hilborn, D. G. Kehler, and C. A. Field, “The variability among populations of coho salmon in the maximum productive rate and depensation,” Ecological Applications, vol. 13, no. 3, pp. 784–793, 2003. View at Publisher · View at Google Scholar
  2. A. L. Jensen, “Harvest reference points for the Beverton and Holt dynamic pool model,” Fisheries Research, vol. 47, no. 1, pp. 93–96, 2000. View at Publisher · View at Google Scholar
  3. M. Holden, “Beverton and Holt revisited,” Fisheries Research, vol. 24, no. 1, pp. 3–8, 1995. View at Publisher · View at Google Scholar
  4. G. Stefansson, “Fish 480 (stockrec) spawning stock, recruitment and production,” November 2005, Course at the Department of Biology of the University of Iceland using data of the Marine Research Institute of Reykjavik, Iceland. View at Google Scholar
  5. R. J. H. Beverton and S. J. Holt, “On the dynamics of exploited fish populations,” Fisheries Investment, vol. 19, p. 1, 1957. View at Google Scholar
  6. C. Hui, “Carrying capacity, population equilibrium, and environment's maximal load,” Ecological Modelling, vol. 192, no. 1–2, pp. 317–320, 2006. View at Publisher · View at Google Scholar
  7. M. A. McCarthy, “The Allee effect, finding mates and theoretical models,” Ecological Modelling, vol. 103, no. 1, pp. 99–102, 1997. View at Publisher · View at Google Scholar
  8. M. De la Sen, “Parameter dependent Lyapunov functions for robust stability of time-varying linear systems under point delays,” Applied Mathematics and Computation, vol. 179, no. 2, pp. 612–621, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  9. S. Stević, “A short proof of the Cushing-Henson conjecture,” Discrete Dynamics in Nature and Society, vol. 2006, Article ID 37264, 5 pages, 2006. View at Publisher · View at Google Scholar · View at MathSciNet
  10. S. Alonso-Quesada and M. De la Sen, “Robust adaptive control of discrete nominally stabilizable plants,” Applied Mathematics and Computation, vol. 150, no. 2, pp. 555–583, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  11. M. De la Sen, “Robust stable pole-placement adaptive control of linear systems with multiestimation,” Applied Mathematics and Computation, vol. 172, no. 2, pp. 1145–1174, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  12. M. De La Sen and S. Alonso, “Adaptive control of time-invariant systems with discrete delays subject to multiestimation,” Discrete Dynamics in Nature and Society, vol. 2006, Article ID 41973, 27 pages, 2006. View at Publisher · View at Google Scholar · View at MathSciNet
  13. G. Feng, “Analysis of a new algorithm for continuous-time robust adaptive control,” IEEE Transactions on Automatic Control, vol. 44, no. 9, pp. 1764–1768, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  14. Z. Huang, S. Chen, and Y. Xia, “Incorporate intelligence into an ecological system: an adaptive fuzzy control approach,” Applied Mathematics and Computation, vol. 177, no. 1, pp. 243–250, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  15. R. Fazio and A. Jannelli, “Mathematical and numerical modeling for a bio-chemical aquarium,” Applied Mathematics and Computation, vol. 174, no. 2, pp. 1370–1383, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  16. F. Chen and C. Shi, “Global attractivity in an almost periodic multi-species nonlinear ecological model,” Applied Mathematics and Computation, vol. 180, no. 1, pp. 376–392, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  17. M. De la Sen, “On impulsive time-varying systems with unbounded time-varying point delays: stability and compactness of the relevant operators mapping the input space into the state and output spaces,” The Rocky Mountain Journal of Mathematics, vol. 37, no. 1, pp. 79–129, 2007. View at Publisher · View at Google Scholar · View at MathSciNet
  18. M. De la Sen and N. Luo, “On the uniform exponential stability of a wide class of linear time-delay systems,” Journal of Mathematical Analysis and Applications, vol. 289, no. 2, pp. 456–476, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  19. X. Ding and H. Su, “Dynamics of a discretization physiological control system,” Discrete Dynamics in Nature and Society, vol. 2007, Article ID 51406, 16 pages, 2007. View at Publisher · View at Google Scholar · View at MathSciNet
  20. N. Hritonenko, A. Rodkina, and Y. Yatsenko, “Stability analysis of stochastic Ricker population model,” Discrete Dynamics in Nature and Society, vol. 2006, Article ID 64590, 13 pages, 2006. View at Publisher · View at Google Scholar · View at MathSciNet
  21. N. Kamel, A. Mousavi, and A. Samraj, “Whitening of background brain activity via parametric modeling,” Discrete Dynamics in Nature and Society, vol. 2007, Article ID 48720, 11 pages, 2007. View at Publisher · View at Google Scholar
  22. M. De la Sen, “The environment carrying capacity is not independent of the intrinsic growth rate for subcritical spawning stock biomass in the Beverton-Holt equation,” Ecological Modelling, vol. 204, no. 1-2, pp. 272–273, 2007. View at Google Scholar