State and Parameter Estimation Based on Filtered Transformation
for a Class of Second-Order Systems
Mehdi Tavan, Kamel Sabahi, and Saeid Hoseinzadeh


Abstract— This paper addresses the problem of state and
parameter estimation for a class of second-order systems with
single output. A new filtered transformation is proposed for the
system via dynamic vector and matrix. In this method, the
dynamics of the vector and matrix are derived by immersion
and invariance technique such that the state estimation
condition is guaranteed. Compared to the classical approaches
that persistency of excitation (PE) condition is required for
parameter convergence, the proposed method needs a weaker
one, so called non-square-integrability condition, in the
transformation via dynamic matrix. Simulation results are
concluded for a class of regressors which are not PE but satisfy
the new condition.

I. INTRODUCTION
The persistency of excitation (PE) condition on the
regressor is one of the main tools to ensure convergence in
parameter estimators for identification purposes. Recently, a
new procedure has introduced in [1] to design the parameter
estimators in which the PE condition is replaced with a
weaker one, called non-square-integrability condition. In this
procedure some dynamic regressors are designed to
guarantee the new condition for the determinant of an
extended regressor. Applications of this method for some
electrical examples show better transient and robust
performance properties [2, 3]. The use of dynamic regressors
can be found for adaptive observer design via filtered
transformation in [4] and the references cited therein. Besides
the PE condition for parameter convergence, an additional
disadvantage is overparameterisation in such a design. The
immersion and invariance (I&I) based filtered transformation
introduced in [5] reduces the order of the observer with
elimination of output injection error term in update law.
Application of this method for DC-DC boost converter is
evaluated via experiments in [6]. The filtered transformations
in [4, 5] are formed via a dynamic vector which yields to a
descent algorithm with singular descent rate that is at the
core of convergence problem in parameter estimation [7].
In this paper a filtered transformation is introduced via
dynamic matrix for a class of second-order systems with
single-output. The main contribution of this paper is to
propose an estimator based on the new transformation, which
asymptotically estimates unavailable state and unknown
parameters. The parameter convergence is ensured under the
non-square-integrability condition. Fulfilling this condition is
discussed for the estimator. Moreover, it is shown that the
estimator provides an arbitrary convergence speed. Also
M. Tavan and Saeid Hoseinzadeh are with the Department of Electrical
Engineering, Mahmudabad Branch, Islamic Azad University, Mahmudabad,
Iran (e-mail: m.tavan@srbiau.ac.ir).
K. Sabahi is with the Department of Electrical Engineering, Mamaghan
Branch, Islamic Azad University, Mamaghan, Iran.

simulation results are concluded for a class of regressors
which are not PE but satisfy the new condition.
II. PROBLEM STATEMENT AND PRELIMINARY
Consider a class of single-output nonlinear systems that
can be described in suitable co-ordinates by equations of the
form
̇

(

)

(
̇

)

(

)

(1)

(

)

(2)

where
is the measurable state,
is the
unmeasurable state,
is a constant vector of unknown
( ) are known
parameters,
( ) ( ) ( ) and
mappings. As mentioned in Introduction, our objective
consists of estimating and . As shown in [6], the average
model of some popular switch converters, such as DC-DC
boost, buck and buck-boost, can be represented in the form of
above equations.
For the system, the unmeasurable state can be assumed
to be generated by a parameter dependent change of
coordinate
[4]
[

(3)

]

where ( )
is an auxiliary dynamic vector, which
[
] is the
its dynamics are to be defined, and
unavailable vector. From (2) and (3) the dynamics of can
be written as
̇

[

]

(4)

with
(

[

̇)

]

where
and
are the zero matrix and zero
vector, respectively. Although using the vector filtered
transformation ( ) for such systems like in [5, 6] facilitate
the state estimator design, the limitation imposed by the PE
condition is still remained in the parameter estimation. To
overcome this problem the matrix filtered transformation is
introduced as
[
where
()

(5)

]

, in which
(
is identity matrix,

)
( )

,

,
is

an auxiliary dynamic matrix, which its dynamics are to be
[
] is the new unavailable vector. In
defined, and
turn, the dynamics of can be obtained as
̇

[ ]

(6)

̇

[

(

(
where
( )

)

̇)

(

]

)
.

and

(

)

As will become clear in the next section, the derivation of
the ( )-dynamics and the ( )-dynamics follows the idea
in [5, 6]. In the proposed approach, the advantages of I&I
technique with input-output filtered transformation is
combined to construct a proper immersion and an auxiliary
dynamic filter for the estimation purpose.

In order that the state be observable from the output ,
the following assumption is imposed.
Assumption 1: The system (1), (2) is forward complete,
i.e. trajectories exist for all
. The function ( ) and
its time derivative are bounded and | |
for all
and
.
A. Filtered Transformation via Dynamic Vector
This subsection is devoted to summarize the method
proposed in [5, 6] with some modifications. The change of
coordinate (3) and the dynamic equation (4) are considered
for design methodology in this subsection. Following the I&I
procedure [8] and [5, 6], define the estimation error as
(

)

(7)

where
is the estimator state, the dynamics of which
is to be defined, and the mapping
yet to
be specified. To obtain the dynamics of , differentiate (7),
yielding
̇

̇
̇

̇

̇

̇
Let

(

)[

]

(

)

(8)

]

]

)

(10)

[
], the
From (7) and the definition of
[
] with
estimation error can be decomposed into
and
. Also, with respect
to (3), the state estimation error can be expressed as ̅
. The following proposition establishes that there exists
the mapping and the dynamics of the filtered vector such
that the state estimation error ̅ converges to zero.
Proposition 1:
Assumption 1 with
(

)

Consider

( )

[

(

)

| |
̇

the

system

(10)

under

]

( )[

(11)

] ( )

(

(

)

)

(12)

where
and
are constant and
( )
positive definite matrices,
is a strictly
increasing function with bounded derivative. Then the system
(10) has a uniformly globally stable equilibrium at the origin,
and ̅
. If, in addition, ( ) is bounded, then
( ) is bounded and ̅ converges to zero.
Proof 1: Notice that the derivation of the signum
( ) is zero, due to the observation restriction
function
()
. Substituting (11), (12) into (10) yields the error
dynamics
̇

( ̇

[

)

]
(13)

| |

[

]

( ) has been used to get the last equality.
where | |
Consider now the following Lyapunov function candidate
( )

(

)

(14)

whose time-derivative along the trajectories of (13) satisfies
̇

where ̇ is given in (4), and the output dynamics ̇ can be
rewritten in term of , as

[

(9)

[

(

III. ESTIMATOR DESIGN
In this section, two classes of estimators are designed for
the system (1), (2). Although, the first one facilitates
estimator design using a filtered transformation via dynamic
vector, it requires the regressor persistency of excitation for
the parameter estimation. It is worth pointing out that the
second one ensures parameter convergence with PE
dismissal.

)

where ̇ is to be defined. Replacing (9) in (7), yields
̇

(

)(

̇

with
[

]

| |

(

| |

(

)
(

) )

(15)

for some constant
, where the positive definite
property of has been used to get the bound. The equations
(14) and (15) imply that the system (13) has a uniformly
globally stable equilibrium at the origin,
,
,
and then ̅
. Note that the above holds for
any . It is clear that, (12) is input-to-state stable dynamic,

then ( ) and its time derivative remains bounded for any
⁄
bounded . Consequently, , , ̇, and
being
̈
bounded, ̇ in (13) and then
in (15) are bounded.
Therefore, the convergence of ̅
to zero follows
directly from Lemma 1 in [9].
] converges to zero, an
Remark 1: If ̅ [
](
).
asymptotic estimate of is given by ̂ [
As is well known ([10], Lemma B.2.3), the parameter
estimation error
converges to zero if the regressor vector
( ) and its time derivative are uniformly bounded and the
regressor satisfies PE condition. Notice that, the bounded
condition is satisfied here for any bounded .
Remark 2: In comparison with the estimator proposed in
[5, 6], the gain matrix is added to filter dynamics in (12)
which can help to make the regressor exciting over a time
sequence arbitrary.
Remark 3: In comparison with the mapping in [5, 6],
the gain matrix is added in (11) and then in error dynamics
(13). The system (13) is in the form of damped nonlinear
oscillator in which the transient time can be improved
arbitrary by reducing the overshot via tuning.
B. Filtered Transformation via Dynamic Matrix
To dismissal the PE condition in parameter estimation,
the design methodology in this subsection is based on the
change of coordinate (5) and the dynamic equation (6). In
the light of the proposed design procedure in the previous
subsection, a general expression of ( ) can be expressed
by
(
(

)

( )[

)

( )

( (

)

( (

)[

(

| |

( )(

(

̇

| |

(

)
[

)

̇
(

(18)

))
.

](
[ ]

(21)

[
], the
From (16) and the definition of
estimation error can be decomposed into
and
. Also, with respect to (5), the state
estimation error of
can be stated as ̅
. A
choice of and
which achieves the design objective is
given in the following proposition.
Proposition 2: Consider the system (21) verifying
Assumption 1, with
( )
( )

( )

(22)

( )

(23)

where ,
and ( ) are chosen as in Proposition 1, in
addition, suppose that the eigenvalues of
are different.
Then the system (21) has a uniformly globally stable
equilibrium at the origin,
and ̅
. Moreover, if
( ) is bounded, then ( ) is bounded and ̅ converges to
zero. If in addition
( )
(resp.
( )
), then
converges to zero asymptotically (resp. exponentially).
Proof 2: Consider the following Lyapunov function
candidate
( )

(

)

(24)

whose time-derivative along the trajectories of (21), is
̇

| |

(17)

))
]

[

]

[

]

(25)

Replacing (23) and (24) in (25), we obtain

] ( )

( )
( )
where
choosing the update laws as
̇

[

(16)

where ( )
is to be defined, and ( )
( )
is a vector function in which, ( ) is defined
in Proposition 1. Notice that, the dynamics of the vector
function ( ) can be obtained in term of , as
̇

| |
̇

Now,

by

[ ]

] ( ))
)

(19)

̇

| |

(

)

(26)

which implies that the system (21) has a uniformly globally
stable equilibrium at the origin,
,
,
and then ̅
. The above holds for any
. Inserting
(23) and (24) into (20) yields
̇

| |

(

)

(

)

(27)

which is input-to-state stable with respect to perturbation
( ). Therefore, the dynamic equation (26) with
̇
induces
. On the other hand, for bounded ,
⁄ , boundedness of ̇ and ̈ can be obtained
, ̇, and
from (21)-(23) and (26) respectively. As a result, the
convergence of
and then ̅ to zero is guaranteed by
Lemma 1 in [9]. Note now that
(28)

(20)

and with respect to (6), (16)-(18), the resulting error
dynamics can be obtained as

for some
, where use has been made of positive
definite property of . For any ( )
, the expression
( ) ( ) is symmetric, and its eigenvalues are real and
nonnegative [11]. Also, if ( ) are the eigenvalues of ( ),
then ( ) are the eigenvalue of ( ) ( ), for all

. Let
, then

( ) be the minimal value of
( )

( )

( )

( ) for any

(29)

On the other hand, it is well known the determinant of a
square matrix is the product of all its eigenvalues. As a
( )
result, if
( )
then
. Now using (26)(29) and doing some basic bounding, yields
̇

| |

(

Remark 4: Replacing (23) and (24) in (21) yields to the
following error dynamics
[

IV. SIMULATION EXAMPLE

(30)

(
). It is clear that
for some
and ( )
exponential convergence of , and then can be directly
obtained from (30) due to the positiveness of ( ). Note
( ), if
now that for bounded
( )
then all
( ) ( )
eigenvalues ( )
[1]. As a result,
and from (30) asymptotic stability of
and then
is
guaranteed.

| |

Remark 8: Notice that, a sudden change in the system
parameters distorts the output signal in (1), (2). This can be
transferred to the
( )-dynamics in (27) by an output
⁄ . Then, the
feedback which appears in the term
dynamics are actuated to adapt ( ) with the changes, even
if the regressor ( ) is not dependent on the output signal,
i.e. ( ).

)

( )

̇

Generally, a necessary condition for any positive definite
matrix in (27) is that its eigenvalues be different.

]

in which the effect of overshoot induced by
be reduced arbitrarily by tuning .

Consider the system (1), (2) with
[
] and
( )]. Assume that
[
( )
is derived from the following equations

(

Remark 5: If
converges to zero, an
asymptotic estimate of is given by ̂
. Also, If
̅ [
] converges to zero, an asymptotic estimate of
). In turn, with respect to
is given by ̂ [
](
the definition of
, an estimate of unmeasurable state is
given by ̂
̂.
C. Discussion
The main question arises at this method is that, for a
( ) how to fulfill the condition
given regressor
( )
. The answer lies in the ( )-behavior tuning.
With respect to the first order differential equation (27), the
behaviour of ( ) consists of transient component and
steady state component. These components can be affected
by arbitrarily chosen of a possible output feedback, induced
by the term ⁄ , and the gain . In addition, the transient
component also depends on the initial value of ( ), i.e.
( ). Hence, some notes are remarked following.
Remark 6: In order to fulfil the condition
( )
, it
(
))
is reasonable to choose ( ) such that
(
.
However, the effect of this choice fades away over time.
Remark 7: The equation (27) becomes a simpler
differential equation if the positive definite matrix
becomes diagonal. In this case, (27) can be decomposed to
individual differential equations. Now, to satisfy the state
( )
, a necessary condition is that the entries on the
main diagonal of be different. This yields different steady
state components for each individual differential equation.

̇ ( )

( )

̇( )

(

)

(

( )

) ( )

(32)
(33)

where
are positive constant and
, and
( ) belongs to the following set of differentiable function
{

(31)
can

)

( )

| ̇( )
()
( )

()

( )
(34)

The set
is borrowed from [1] with some modifications.
( )
For any
, the condition ( ) ( )
insures
and then ( ) remains bounded. Also, the condition
()
( )
( ) and then ( ) is not PE
yields
because it converges to zero as
.
⁄
Let us assume
in (12) and (27). Then, the
steady state component of (12)-solution for
and
can be obtained as
[
[

(

)

(

)

( )]

(35)

( )

( )]

(36)

respectively. With respect to Remark 7, choosing
(
) for (27) yields to
[
where

]

(37)

is the steady state component of (27)-solution.

Remark 9: The regressors
because
( )
( )
(
)
, and
(
)
due to positive semi definite property of the
matrixes. Notice that, the matrix
, but
(
)
because
(
in which, by (34), ̇ ( )

)

̇( )
(

)(

)

.

A. Simulation Results
Consider the system introduced in this section with

(38)

( )

√

(39)
[3]

which
[1]. Simulations are run to evaluate the
performance of estimators introduced in Proposition 1 and 2.
All initial values correspond to the system and the estimators
are set to zero.
The time histories of the parameters estimation errors
̂ and the state estimation errors
̂ are shown for
estimator b1 (dashed line), estimator b2 (dashed-dot line),
and estimator B (solid-line) in Fig. 1. The estimator bi is
formed by replacing
in Proposition 1 and the
estimator B is formed by replacing
(
) in
Proposition 2, while the mapping ( )
and the gain
are chosen for all of them. The graphs in Fig. 1 show
)
the performance of the estimators for the set (
(
).

[4]

[5]

[6]

[7]
[8]

From the graphs in Fig. 1 it can be seen that, for
,
a bigger value of
has a better response in the parameter
estimation. In the all graphs, the estimator B has a faster
response in parameter and state estimation. The last graph of
Fig. 1 shows the absolute value of the ( )-determinant.
From (38) and (39) the determinant of ( ) is not square
integrable which is the core of convergence property of the
estimator B.

[9]
[10]

[11]

regressor extension and mixing," IFAC-PapersOnLine, vol. 49,
pp. 99-104, 2016.
A. Bobtsov, D. Bazylev, A. Pyrkin, S. Aranovskiy, and R.
Ortega, "A robust nonlinear position observer for synchronous
motors with relaxed excitation conditions," International
Journal of Control, vol. 90, pp. 813-824, 2017.
R. Marino and P. Tomei, "Global adaptive observers for
nonlinear systems via filtered transformations," IEEE
Transactions on Automatic Control, vol. 37, pp. 1239-1245,
1992.
M. Tavan, A. Khaki-Sedigh, M.-R. Arvan, and A.-R. Vali,
"Immersion and invariance adaptive velocity observer for a class
of Euler–Lagrange mechanical systems," Nonlinear Dynamics,
vol. 85, pp. 425-437, 2016.
M. Malekzadeh, A. Khosravi, and M. Tavan, "Immersion and
invariance Based Filtered Transformation with Application to
Estimator Design for a Class of DC-DC Converters,"
Transactions of the Institute of Measurement and Control, vol.
(Accepted), 2018.
R. Brockett, "The rate of descent for degenerate gradient flows,"
in Proceedings of the 2000 MTNS, 2000.
A. Astolfi, D. Karagiannis, and R. Ortega, Nonlinear and
adaptive control with applications: Springer Science & Business
Media, 2007.
G. Tao, "A simple alternative to the Barbalat lemma," IEEE
Transactions on Automatic Control, vol. 42, p. 698, 1997.
R. Marino and P. Tomei, Nonlinear control design: geometric,
adaptive and robust: Prentice Hall International (UK) Ltd.,
1996.
C. D. Meyer, Matrix analysis and applied linear algebra vol. 71:
Siam, 2000.

The next simulation is run to show the influence of the
gain
on the performance of the estimator B. All
conditions and gains are the same as previous simulation.
The time histories of the parameters and the state estimation
errors are shown for
(dashed line),
(dasheddot line), and
(solid-line) in Fig. 2. As it is
evident in Fig. 2, the overshot is decreased by increasing the
value of , and as a result, faster responses in the both
parameter and state estimation are obtained.
V. CONCLUSION
A new I&I based filtered transformation is introduced via
dynamic matrix for a class of second-order systems with
single-output. The main contribution of this paper is to
propose an estimator based on the new transformation,
which asymptotically estimates unavailable state and
unknown parameters. Parameter convergence is guaranteed
under non-square-integrability assumption for the matrix.
Moreover, it is shown that the estimator provides an
arbitrary convergence speed. Also simulation results are
concluded for a class of regressors which are not PE but
satisfy the new condition. Some practical results on
electrical systems are under way to show effectiveness of the
proposed design.

Figure 1. Transient performance of the estimator bi and estimator B.

REFERENCES
[1]

[2]

S. Aranovskiy, A. Bobtsov, R. Ortega, and A. Pyrkin,
"Performance enhancement of parameter estimators via dynamic
regressor extension and mixing," IEEE Transactions on
Automatic Control, vol. 62, pp. 3546-3550, 2017.
S. Aranovskiy, A. Bobtsov, R. Ortega, and A. Pyrkin, "Improved
transients in multiple frequencies estimation via dynamic
Figure 2. Transient responces of the estimator B for changs in .

