- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Discrete Dynamics in Nature and Society

Volume 2008 (2008), Article ID 793512, 21 pages

http://dx.doi.org/10.1155/2008/793512

## 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.

#### Abstract

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., [1–9]). 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 [1–9] with being the set of natural numbers, where , ( 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, and ,
for some ,
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 is sufficiently unfavorable to extinguish the
population since, from (1.1), a recruitment 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 resulting directly from (1.1) in the time-varying discrete linear system where , , and , 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], The IBHE (2.2) corresponds, through the change of variable , to an extended BHE which includes a sequence of additive disturbances: where Note that the sequence is an additive disturbance of the BHE which is identically zero if the additive disturbance sequence 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. * and , for all , and some .

*Assumption 2.2. *, for all , and is known.

Note that Assumption 2.1 implies that the BHE is stable since , and also controllable since so that , for all , so that in (2.2), and then the BHE solution 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 , then it always holds for some . 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 is uniformly bounded, then all solutions of
the BHE and the IBHE are uniformly bounded and positive, provided that .

*Proof. *.
Define .
Since ,
it follows from (2.2) that

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

*(1) A Reference*

BHE , for all ,
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 ,
with reference input ,
and parameter sequences and .

*(2) The Current*

BHE and *its corresponding* IBHE whose parameters are , are deviated from the suited ones and ,
which may include disturbance terms grouped in a sequence with control input obtained according to a feedback law.

*(3) The Control Law:*

provided
that is known, which is a relaxation of Assumption
2.2 which is parameterized by with .
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 and and a correcting sequence calculated from (2.8). The subsequent constraint 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 ,
then the control law (2.7)-(2.8) achieves an identically zero tracking error .*

*Proof. * It
follows, by complete induction, by assuming , for all ,
that

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 ,
provided that the disturbance sequence is upper bounded by a bounded sequence.*

*Proof. *Direct
calculations yield so that, from (2.7)-(2.8), where , for all , and

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, , and , for all . 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 .
Then the following properties hold. *

(i)* The control
parameter sequence is subject to the constraint , for all , which induces the parametrical constraints *(ii)* The tracking-error
sequence is subject to the constraint for all , which is trivially bounded, provided that the
disturbance sequence is bounded.*(iii)* A
positive carrying capacity is obtained from the controller synthesis if when is unknown but is known. The control parameter and the calculated from the control law are positive if is calculated satisfying with *

*Proof. *(i) Note that which
yields directly Property (i) since .

(ii)One gets
from (2.6)–(2.8) and Property (i) that 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., [17–19]). One takes the advantage that, since the parameters fulfill the constraint , for all , 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 ,
implying , are incorporated within
a disturbance sequence which can also include other unstructured
disturbance contributions in a sequence : and .
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.

*IBHE*

If ,
the resulting particular case of (3.2) is called the nominal IBHE, and its
corresponding inverse equation, whose solution is the sequence ,
is called the nominal BHE.

*Estimate of the IBHE*

where and are the estimates of and at the th sample.

*Identification Error*

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, with and , is the reference environment carrying capacity
inverse.

*Tracking Error*

*Adaptation Algorithm 1:*

is a recursive estimation of the
intrinsic growth rate, where and are
real sequences satisfying , 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 in the parameter-adaptive algorithm.

*Adaptive Control Law:*

with with the
sample-dependent controller parameters 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 , .
The above results are independent of the parameter-adaptive algorithm.*

*Proof. * Proceed, by induction, by assuming for any given .
Then one gets, from (3.8)-(3.9) into (3.3) by using , if and then using (3.6), 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 satisfying , for all ,
is known that the adaptation algorithm
(3.7) is updated by choosing ,
and **for
all , for some prefixed real constants , , and .
Then the following properties hold irrespective of the control law. *

(i)* The sequences , , and are bounded and have limits ,
and zero,
respectively.*(ii)*There exist
the subsequent limits *

*Proof. * Let be defined by .
Note, by subtracting in both sides of (3.7), that , where .
Since from
(3.4), for any absolute upper-bound ,
one
gets From (3.15), . Assume so that for any real , one gets, from (3.15), for all , provided
that and satisfies for any real .
Then ,
what implies that the sequences , and ,
and are
uniformly bounded and there exist limits , , , and . so that since either
or , for all , and Property (i) is proven. Since converges to a finite limit, Property (ii)
follows from 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 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 , ,
such that **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:*

being updated by
choosing ,
and for all ,
for some prefixed real constants , ,
and .*A Posteriori Estimation Via Projection on the Interval :*

Then Theorem 3.2(i)-(ii) holds. In addition, similar
properties hold by replacing

*Proof. * Define the sequence from (3.20)-(3.21) using a close
development to that in Theorem 3.2 so that it is monotonically decreasing for all since , for all , from the projection part (3.22) of the
estimation algorithm. Then is bounded and converges to a finite limit,
and is also bounded and converges to zero. Also for all , and . Since , for all ,
then there exists a finite limit , then converges to a finite limit and converges since converges (since converges) to finite limits as .
Thus Theorem 3.2(i)-(ii) follows. On the other hand, note that Since the normalized sequences and 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: The solution of the IBHE is equivalently calculated either as a function of the nominal parameters and or as a function of their estimates as follows, by considering the adaptive control law (3.8): 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 being sufficiently deviated from unity and 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. *, for all , and some prefixed real constants and . A
direct calculation with Assumption 3.4 yields while
similar expressions follow for replacements involving any two elements of the quadruplet as follows: The
above expressions (3.29)-(3.30) lead to the following result referring to explicit positive lower and
upper bounds of so that since .

Proposition 3.5. * Assume , for all .
Then the following constraints hold: *

*Proof. *It
follows directly from (3.29)-(3.30). The necessary condition is needed to keep in the right-hand-side upper bound of , for all 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 with known constants of the related linear upper bound and with sufficiently slow growth of with respect to . 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 such that

Then, from (3.28b), (3.31)-(3.32), and Assumptions 3.6–3.4, one gets 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 of the additive disturbances.

*Assertion 3.7. * If the Adaptation Algorithm 1 is used, then , for all , for some real constants and .

*Proof. *If ,
the assertion is true with and . If , then, from Theorem 3.2, (i) , for all .
Thus and
the assertion follows with since .

*Assertion 3.8. * If the Adaptation Algorithm 2 is used and Assumption 3.4 holds, then , for all , with

*Proof. *It follows from Assumption 3.4, (3.4) and (3.29), and the fact that is positively lower bounded, which yield
directly

From Assertions 2.3 and 3.7, for the Adaptation Algorithms 1 and 2, respectively, Assumption 3.6 and (3.28c), for any and some positive finite constants , and satisfying , , and . The last two inequalities of (3.39) follow directly by proceeding recursively with the preceding one by using since . It is obvious that, for each , it exists a unique , dependent in general on , such that such that and there is no other such that , where the same maximum is reached. A precise definition is with “” and “” being the conjunction and negation logic symbols. The substitution of (3.40) into (3.39) yields provided that , which holds if and only if (i.e., if Assumption 3.6 holds with sufficiently small ), where 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, , , , , and so on. The above result also holds for the Adaptation Algorithm 2 with specific constants redefined by and 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.3–3.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 ,
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., ) 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 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 characterizing ,
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., [20–22])
and to extend the research to governing extinction conditions [10].

#### 4. Numerical Examples

*Example 4.1. *The objective is that
the BHE solution tracks a reference sequence 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 and to generate the sequence associated to the desired to be tracked by the BHE solution. The
reference values are ,
and , for all , are used for the reference model and BHE,
respectively, with the period (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 ,
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 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.

*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.

#### 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.

#### Acknowledgments

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.

#### References

- 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 - 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 - M. Holden, “Beverton and Holt revisited,”
*Fisheries Research*, vol. 24, no. 1, pp. 3–8, 1995. View at Publisher · View at Google Scholar - 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
- 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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