Design of Indirect MRAS-based Adaptive Control Systems
Design of Indirect MRAS-based Adaptive Control Systems
when
we
use
the
Adjustable model
The reference model, in this case referred to as the
adjustable model, will follow
ow the response of the process.
In the following discussions the terms adjustable model
and adaptive observer are used interchangeably. The goal
in process identification is to obtain a satisfactory model of a
real process by observing the process inpu
input-output behavior.
Identification of a dynamic process contains four basic steps
[1], [2]. The first step is structural identification, which
allows us to characterize the structure of the mathematical
model of the process to be identified. This can be done
don from
the phenomenological analysis of the process. Next, we
determine the inputs and outputs. Third step is parameter
identification. This step allows us to determine the
parameters of the mathematical model of the process.
Finally, the identified model is validated. When the
parameters of the identified model and the process are
supposed to be identical, the model states can be
considered as estimates of the process states. When the
states of the process are corrupted with noise, the structure
of the adaptive
daptive observer can be used to get
filtered estimates of the process states. When the input signal
itself is not very noisy, the model states will also be almost
free of noise. It is important to notice that in this case the
filtering is realized with minimum
imum phase lag [3]. However,
this adjustable is also able to deal with unknown or time
varying parameters [6].
In order to design indirect adaptive MRAS, two processes
will be followed: Firstly, determining the adaptive law for
variable parameters of
,
of the adjustable model.
Next, designing of PD
D adaptive controller based on
, .
Step1:Determining
Determining the process and adjustable model
constructions.
Any system can be described by either its transfer function
or its state space. In this case, the second order process is
given by the differential equation.
x1p x2p
x2p ap x2p bp uc
It would be rewritten in the state space form:
x1p
x1p
xp ; xp ; xp A p xp Bp uc
x2p
x2p
where:
0 1
0
Ap
;Bp b
0
a
p
where:
0 1
0
Am
;Bm b
0
m
m
Step2:Deriving the error equation.
e xp xm
de dxp dxm
dt
dt
dt
de
A p xp Bp uc A m xm Bm uc
dt
de
A p x p x m A p A m xm
dt
+ Bp Bm uc
de
A p e Axm + Buc
dt
where:
0
0
A A p Am
0 am ap
0
B Bp B m
bp bm
0 x1m
0
de
Ape
+
dt
0 am ap x2m
0
de
Ape
dt
x2m
0
b b uc
p m
0 am ap
uc bm bp
Therefore:
de
A p e o
dt
1
V (e) eT Pe ( o )T E ( o )
2
Where:
P is an arbitrary definite positive symmetrical matrix;
E is a diagonal matrix with positive elements which
determine the speed of adaptation.
p12
0
p
e
P 11
; E 11
p21 p22
0 e22
Step4:Determine
the
conditions
under
definite negative.
1
V (e) eT Pe ( o )T E ( o )
2
dV (e) 1 deT
de
Pe eT P
dt
2 dt
dt
1 d ( o )T
d ( o )
E ( o ) ( o )T E
2
dt
dt
1
eT Ap Pe eT PApT e
2
T
T
d
o T Pe o E
dt
dV (e) 1 T
e PApT Ap P e
dt
2
d
o T
T Pe E
dt
dV (e)
1
eT Qe
dt
2
The result ofthe adjustment law to be:
d1 dam
d
d dt dt
Pe E
0;
dt
dt d 2 dbm
dt dt
T
With:
0
x2m
0
uc
am 0 ap
, b
bm
p
T Pe E
E
which V (e) is
d
0
dt
d
d
T Pe
E 1 T Pe
dt
dt
e '
22
0
0 0 x2 m p11
e '11 0 uc p21
p12 e1
p22 e2
0 x2 m p21 x2 m p22 e1
e '
d
22
0
e
'11 uc p21 uc p22 e2
dt
d1 dam
T s
o 2
s 2 2o s o 2
Where:
o is a specific frequency.
is a damping factor.
PA
Ap P Q
0 1
q
Ap Am
; Q 11
q21
0 am
PA
From:
p11
p
21
Ap P Q
p12 0 0 0 1 p11
p22 1 am 0 am p21
p
12
p22
am p12 p21
am p22 am p21
p p21
12
p22 am p21
q12
q22
p12
Q
p22
p22
Q
am p22
am p12 p22
q
11
2am p22
q21
q12
q22
1 1
1
p21
q21 ; p22
am 2am
2am
J (t ) dt where R xm
0
a p bp .K d 2 o
After some calculations:
2o a p
2
K p o ; Kd
bp
bp
Setting
o 50; 0.7
Hence:
Kp
70 a p
2500
; Kd
bp
bp
70 a p
2500
; Kd
bp
bp