Overall Objectives
Scientific Foundations
Application Domains
New Results
Contracts and Grants with Industry
Other Grants and Activities

Section: Scientific Foundations

Keywords : Computer Algebra, Control, Cryptography, Diagnosis, Dynamical Systems (dynamic system, dynamic system, dynamic system, dynamic system, dynamic system), Estimation, Fault-Tolerant Control (fault tolerance control, fault tolerance, control), Identification, Robotics, Signal Processing, Image Processing, Video Processing.

Fast parametric estimation and its applications

Parametric estimation may often be formalized as follows:

y= F( x, $ \upper_theta$) + n,(1)


Finding a “good” approximation of the components of $ \upper_theta$ has been the subject of a huge literature in various fields of applied mathematics. Most of those researches have been done in a probabilistic setting, which necessitates a good knowledge of the statistical properties of n . Our projectis devoted to a new standpoint which does not require this knowledge and which is based on the following tools, which are of algebraic flavor:

Linear identifiability

In most problems appearing in linear control as well as in signal processing, the unknown parameters are linearly identifiable : standard elimination procedures are yielding the following matrix equation

Im2 ${P\mfenced o=( c=) \mtable{...}=Q,}$(2)


How to deal with perturbations and noises?

With noisy measurements equation ( 2 ) becomes:

Im4 ${P\mfenced o=( c=) \mtable{...}=Q+R,}$(3)

where  R is a  r×1 column matrix, whose entries are finite linear combination of terms of the form  Im5 ${t^\#957 \mfrac {d^\#956 \#951 }{dt^\#956 },\#956 ,\#957 \#8805 0}$ , where  $ \eta$ is a perturbation or a noise.

Structured perturbations

A perturbation $ \pi$ is said to be structured if, and only if, it is annihilated by a linear differential operator of the form  Im6 ${\#8721 _\mtext finitea_k{(t)}\mfrac d^k{dt^k}}$ , where ak( t) is a rational function of  t , i.e.  Im7 ${\mfenced o=( c=) \#8721 _\mtext finitea_k{(t)}\mfrac d^k{dt^k}\#960 =0}$ . Note that many classical perturbations like a constant bias are annihilated by such an operator. An unstructured noise cannot be annihilated by a non-zero differential operator.

By well known properties of the non-commutative ring of differential operators, we can multiply both sides of equation ( 3 ) by a suitable differential operator $ \upper_delta$ such that equation ( 3 ) becomes:

Im8 ${\#916 P\mfenced o=( c=) \mtable{...}=\#916 Q+R^',}$(4)

where the entries of the  r×1 column matrix R' are unstructured noises.

Attenuating unstructured noises

Unstructured noises are usually dealt with stochastic processes like white Gaussian noises. They are considered here as highly fluctuating phenomena, which may therefore be attenuated via low pass filters. Note that no precise knowledge of the statistical properties of the noises is required.


Although the previous noise attenuation (It is reminiscent to what most practitioners in electronics are doing.)may be fully explained via formula ( 4 ), its theoretical comparison (Let us stress again that many computer simulations and several laboratory experiments have been already successfully achieved and can be quite favorably compared with the existing techniques.)with today's literature (Especially in signal processing.)has yet to be done. It will necessitate a complete resetting of the notions of noises and perturbations. Besides some connections with physics, it might lead to quite new “epistemological” issues  [73] .

Some hints on the calculations

The time derivatives of the input and output signals appearing in equations ( 2 ), ( 3 ), ( 4 ) can be suppressed in the two following ways which might be combined:

The numerical values of the unknown parameters  Im9 ${\#920 =(\#952 _1,\#8943 ,\#952 _r)}$ can be obtained by integrating both sides of the modified equation ( 4 ) during a very short time interval.

A first, very simple example

Let us illustrate on a very basic example, the grounding ideas of the ALIEN approach, based on algebra. For this, consider the first order, linear system:

Im10 ${\mover y\#729 {(t)}=ay{(t)}+u{(t)}+\#947 _0,}$(5)

where  a is an unknown parameter to be identified and  $ \gamma$0 is an unknown, constant perturbation. With the notations of operational calculus and  y0= y(0) , equation ( 5 ) reads:

Im11 ${s\mover y^{(s)}=a\mover y^{(s)}+\mover u^{(s)}+y_0+\mfrac \#947 _0s.}$(6)

In order to eliminate the term $ \gamma$0 , multiply first the two hand-sides of this equation by s and, then, take their derivatives with respect to s:

Im12 ${\mfrac d{ds}\mfenced o=[ c=] s\mfenced o={ c=} s\mover y^{(s)}=a\mover y^{(s)}+\mover u^{(s)}+y_0+\mfrac \#947 _0s}$(7)

Im13 ${\#8658 2s\mover y^{(s)}+s^2\mover y^^'{(s)}=a\mfenced o=( c=) s\mover y^^'{(s)}+\mover y^{(s)}+s\mover u^^'{(s)}+\mover u^{(s)}+y_0.}$(8)

Recall that Im14 ${\mover y^^'{(s)}\#8796 \mfrac {d\mover y^{(s)}}{ds}}$ corresponds to - ty( t) . Assume y0= 0 for simplicity's sake (If y0$ \ne$0 one has to take above derivatives of order 2 with respect to s , in order to eliminate the initial condition.). Then, for any $ \nu$>0 ,

Im15 ${s^{-\#957 }\mfenced o=[ c=] 2s\mover y^{(s)}+s^2\mover y^^'{(s)}=s^{-\#957 }\mfenced o=[ c=] a{(s\mover y^^'{(s)}+\mover y^{(s)})}+s\mover u^^'{(s)}+\mover u^{(s)}.}$(9)

For $ \nu$= 3 , we obtained the estimated value a :

Im16 ${a=\mfrac {2\#8747 _0^Td\#955 \#8747 _0^\#955 y{(t)}dt-\#8747 _0^Tty{(t)}dt+\#8747 _0^Td\#955 \#8747 ... d\#955 \#8747 _0^\#955 d\#963 \#8747 _0^\#963 y{(t)}dt-\#8747 _0^Td\#955 \#8747 _0^\#955 ty{(t)}dt}}$(10)

Since T>0 can be very small, estimation via ( 10 ) is very fast.

Note that equation ( 10 ) represents an on-line algorithm that only involves two kinds of operations on u and y: (1) multiplications by t , and (2) integrations over a pre-selected time interval.

If we now consider an additional noise, of zero mean, in ( 5 ), say:

Im17 ${\mover y\#729 {(t)}=ay{(t)}+u{(t)}+\#947 _0+n{(t)},}$(11)

it will be considered as fast fluctuating signal. The order $ \nu$ in ( 9 ) determines the order of iterations in the integrals (3 integrals in ( 10 )). Those iterated integrals are low-pass filters which are attenuating the fluctuations.

This example, even simple, clearly demonstrates how ALIEN's techniques proceed:

A second simple example, with delay

Consider the first order, linear system with constant input delay (This example is taken from  [60] . For further details, we suggest the reader to refer to it.):

Im18 ${\mover y\#729 {(t)}+ay{(t)}=y{(0)}\#948 +\#947 _0H+bu{(t-\#964 )}.}$(12)

Here we use a distributional-like notation where  $ \delta$ denotes the Dirac impulse and  H is its integral, i.e., the Heaviside function (unit step) (In this document, for the sake of simplicity, we make an abuse of the language since we merge in a single notation the Heaviside function  H and the integration operator. To be rigorous, the iterated integration ( k times) corresponds, in the operational domain, to a division by  sk , whereas the convolution with H ( k times) corresponds to a division by  sk/( k-1)! . For  k= 0 , there is no difference and H* y realizes the integration of y . More generally, since we will always apply these operations to complete equations (left- and right-hand sides), the factor  ( k-1)! makes no difference.). Still for simplicity, we suppose that the parameter a is known. The parameter to be identified is now the delay $ \tau$ . As previously, $ \gamma$0 is a constant perturbation, a , b , and $ \tau$ are constant parameters. Consider also a step input  u= u0H . A first order derivation yields:

Im19 ${\mover y¨+a\mover y\#729 =\#981 _0+\#947 _0~\#948 +b~u_0~\#948 _\#964 ,}$(13)

where  Im20 $\#948 _\#964 $ denotes the delayed Dirac impulse and  Im21 ${\#981 _0={(\mover y\#729 {(0)}+ay{(0)})}~\#948 +y{(0)}~\#948 ^{(1)}}$ , of order 1 and support  {0} , contains the contributions of the initial conditions. According to Schwartz theorem, multiplication by a function $ \alpha$ such that $ \alpha$(0) = $ \alpha$'(0) = 0 , $ \alpha$( $ \tau$) = 0 yields interesting simplifications. For instance, choosing  $ \alpha$( t) = t3- $ \tau$t2 leads to the following equalities (to be understood in the distributional framework):

Im22 $\mtable{...}$(14)

The delay $ \tau$ becomes available from k$ \ge$1 successive integrations (represented by the operator H ), as follows:

Im23 ${\#964 =\mfrac {H^k{(w_0+a~w_3)}}{H^k{(w_1+a~w_2)}},~t\gt \#964 ,}$(15)

where the wi are defined, using the notation zi= tiy , by:

Im24 $\mtable{...}$

These coefficients show that k$ \ge$2 integrations are avoiding any derivation in the delay identification.

Figure 1. Delay $ \tau$ identification from algorithm ( 15 )

Figure  1 gives a numerical simulation with  k= 2 integrations and  a= 2, b= 1, $ \tau$= 0.6 , y(0) = 0.3, $ \gamma$0= 2, u0= 1 . Due to the non identifiability over (0, $ \tau$) , the delay  $ \tau$ is set to zero until the numerator or the denominator in the right hand side of ( 15 ) reaches a significant nonzero value.

Again, note the realization algorithm ( 15 ) involves two kinds of operators: (1) integrations and (2) multiplications by  t .

It relies on the measurement of  y and on the knowledge of  a . If  a is also unknown, the same approach can be utilized for a simultaneous identification of a and $ \tau$ . The following relation is derived from ( 14 ):

$ \tau$( Hkw1) + a$ \tau$( Hkw2)- a( Hkw3) = Hkw0,(16)

and a linear system with unknown parameters ( $ \tau$, a$ \tau$, a) is obtained by using different integration orders:

Im25 ${\mfenced o=( c=) \mtable{...}\mfenced o=( c=) \mtable{...}=\mfenced o=( c=) \mtable{...}.}$

The resulting numerical simulations are shown in Figure  2 . For identifiability reasons, the obtained linear system may be not consistent for  t< $ \tau$ .

Figure 2. Simultaneous identification of a and $ \tau$ from algorithm ( 16 )


Logo Inria