Адаптивные системы. Учебное пособие
.pdf
|
|
0 |
1 |
|
|
|
|
|
0 |
|
|
A a2,1 |
a2,2 |
|
, |
B b |
|
, |
|||||
|
|
|
|
|
|
|
|
|
2 |
|
|
неизвестные параметры модели составляют вектор
1, 2 , 3 T
a2,1, a2,2 ,b2 T .
Моделирование осуществляется с помощью дискретной модели x(k 1) Ax(k) Bu(k) Fq(k), x(0) x0 ,
|
|
|
|
|
|
|
|
|
|||||
где A I2 tA |
, B tB |
, |
F tF, |
t 0,1 и матрицы A и B |
|||||||||
линейно зависят от компонент вектора следующим образом |
|||||||||||||
|
|
|
1 |
|
|
t |
|
|
|
|
0 |
|
|
A( ) |
t |
1 t |
2 |
|
, B( ) |
t |
. |
||||||
|
|
|
1 |
|
|
|
|
|
3 |
||||
Математическая модель измерительного комплекса имеет вид y(k) Hx(k) r(k) ,
где H единичная матрица, что соответствует полному измерению вектора состояния, r(k) вектор последовательностей гаус-
совских величин с нулевым средним и дисперсией
0,32 |
0 |
|
|
R |
0 |
0,35 |
. |
|
|
||
Адаптивное управление, формируемое на основе минимизации математического ожидания классического квадратичного функционала
J (t0 ,T ) 12 M
имеет вид
u(k)
T |
T |
|
|
|
|
T |
|
|
|
(t)Cx(t) u |
|
, |
|||||
[x |
|
|
(t)Du(t)]dt |
|||||
t0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Dd |
B |
|
( (k))Sxˆ |
(k) , |
|
|||
|
|
1 |
|
T |
ˆ |
|
|
|
61
ˆ ˆ
где Dd tD , x(k), (k) оценки состояния и параметров, опре-
деляемые с помощью фильтров Калмана (5.12), (5.13). Матрица
(xˆ(k),u(k)) |
и вектор f (xˆ(k),u(k)) |
в (5.13) для заданного вектора |
|||||||||||||||||||||
параметров определяются следующим образом: |
|
||||||||||||||||||||||
|
|
|
|
|
|
(xˆ |
|
|
|
|
0 |
|
0 |
|
|
0 |
|
||||||
|
|
|
|
|
|
(k),u(k)) |
t(xˆ(k)) |
t(xˆ(k)) |
2 |
tu(k) , |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
f |
(xˆ |
(k))1 t(xˆ(k))2 |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(xˆ(k))2 |
. |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
Моделирование |
осуществлялось |
для |
модели, |
где матрицы |
|||||||||||||||||
|
|
, |
|
, |
|
и векторы x(0), |
имеют вид |
|
|
|
|
|
|
|
|||||||||
|
A |
B |
F |
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
0 |
1 |
|
|
|
0 |
|
|
|
|
0,51 |
0 |
|
|
|
|
|
|
|
|
A |
|
0,4 |
|
, |
B |
|
, F |
|
0 |
|
, |
||||||
|
|
|
|
|
|
|
|
|
|
1,4 |
|
|
1,3 |
|
|
|
|
|
0,55 |
||||
x(0) 101 , ( 0, 4; 1, 4;1,3)T
при следующих априорных данных
xˆ(0) |
|
7 |
|
|
Px0 |
1 |
0 |
|
|
|
|
|
1,3 |
|
, |
|
0 |
|
, |
|
|
||
|
|
|
|
|
|
1 |
|
|
|
||
|
|
0,5 |
|
|
|
1 |
0 |
0 |
|
|
|
ˆ |
|
1,7 |
|
, |
P 0 |
|
0 |
1 |
0 |
|
, |
(0) |
|
|
|
|
|||||||
|
|
1,5 |
|
|
|
|
0 |
0 |
1 |
|
|
|
|
|
|
|
|
|
|
||||
и
1,5 |
0 |
|
|
2,295 |
0,734 |
|
C |
1,7 |
|
, D 1, S |
0,772 |
. |
|
0 |
|
|
0,734 |
|
||
На рис. 5.2. |
приведены |
графики переходных процессов |
||||
x(k) x1(k), x2 (k) T , их оценок |
xˆ(k) xˆ1(k), xˆ2 (k) T |
и управлений |
||||
62
u(k) , причем жирная линия изображает управляющие воздействия u(k) , тонкая– состояния x(k) , апунктирная– оценкисостояния xˆ(k) .
0 |
|
|
|
|
|
5 |
|
|
|
|
|
0 |
|
|
|
|
|
0 |
20 |
40 |
60 |
80 |
100 |
Рис. 5.2. График переходных процессов, их оценок и управления
На рис. 5.3. приведены в порядке следования графики оценок векторапараметров (тонкаялиния) иихточныхзначений(жирнаялиния).
0.2 |
|
|
|
|
|
0.4 |
|
|
|
|
|
0.6 |
|
|
|
|
|
0.8 |
|
|
|
|
|
1 |
|
|
|
|
|
1.20 |
20 |
40 |
60 |
80 |
100 |
|
|
|
2 |
|
|
|
|
|
1.5 |
|
|
|
|
|
1 |
|
|
|
|
|
0.5 |
|
|
0 |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
3 |
0 |
20 |
40 |
60 |
80 |
100 |
00 20 40 60 80 100
Рис. 5.3. Графики оценок вектора параметров и ихточные значения
63
2. В модели объекта
x(t) Ax(t) Bu(t) Fq(t), x(t0 ) x0 ,
u(t) (t), u(t0 ) u0 ,
где
|
|
0 |
1 |
|
|
|
|
|
0 |
|
|
A a2,1 |
a2,2 |
|
, |
B b |
|
, |
|||||
|
|
|
|
|
|
|
|
|
2 |
|
|
неизвестным параметром является b2 .
Моделирование осуществляется с помощью дискретной модели
x(k 1) Ax(k) Bu(k) Fq(k), x(0) x0 , u(k 1) u(k) t(k), u(0) u0
где элементы матрицы B линейно зависят от следующим образом:
0 B( ) .
t
В математической модели измерительного комплекса H 1 0 , что соответствует измерению только первой компоненты вектора состояния, r(k) последовательность гауссовских величин с нулевым средним и дисперсией R 0,32 .
Адаптивное управление, формируемое на основе минимизации математического ожидания функционала обобщенной работы
t l t
J (tk ) M k p[xT (t)Cx(t) uT (t)D2u(t)
tk
vT (t)D1v(t) vonT (t)D1von (t)]dt ,
имеет вид
(k) D1 1W2 (k) ,
где W2 (k) является решением в обратном времени системы
64
ˆ |
gM (k lp ) xM (k lp ); |
gM ( j 1) 2gM ( j) AgM ( j) B( (k))u(k), |
T |
( j) tCgM ( j), |
W1 |
|
0 |
|
; |
W1 ( j 1) A W1 |
(k lp ) |
0 |
|
|||
|
|
|
|
|
|
ˆ |
T |
W1 |
( j) tD2u(k), W2 |
(k lp ) 0; |
W2 ( j 1) W2 ( j) B( (k)) |
|
|||
j k lp ,k lp 1,..., k 1, |
|
|
|
|
а xM (k lp ) в первом уравнении определяется с помощью прогнозирующей модели
xM ( j 1) AxM ( j) B(ˆ(k))u(k), xM ( j k) xˆ(k)
,
описывающей поведение объекта моделирования на интервале оптимизации [tk ,tk lp t].
Вектор xˆ(k) и значение |
ˆ |
являются оценками состояния и |
(k) |
параметра, которые определяются с помощью фильтров Калмана
(8.74),(8.75). Матрица (xˆ(k),u(k)) и вектор |
|
f (xˆ(k),u(k)) в (8.75) |
|||||||||||||||||||||
для заданного параметра определяются следующим образом: |
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
0 |
|
|
|
|
(xˆ(k))1 t(xˆ(k))2 |
|
|
||||||||
(xˆ(k),u(k)) |
|
|
|
, f |
|
|
ˆ |
|
|
ˆ |
|
ˆ |
. |
||||||||||
|
|
|
|
|
|
|
|
tu(k) |
ta2,1 (x(k))1 |
(x(k))2 ta2,2 (x(k))2 |
|
||||||||||||
|
|
|
Моделирование осуществлялось |
для |
модели, где |
матрицы |
|||||||||||||||||
|
|
, |
|
|
|
, вектор x(0) |
|
|
|
|
|
|
|||||||||||
|
A |
B |
, F |
и параметр имеют вид |
|
|
|
||||||||||||||||
|
|
|
|
|
|
0 |
1 |
|
|
|
0 |
|
|
|
|
0,51 |
0 |
|
10 |
|
|
||
A |
0,4 |
1,4 |
, B |
|
, F |
|
0,55 |
|
, x(0) |
, 1,3 |
|||||||||||||
|
|
|
|
|
|
|
1,3 |
|
|
0 |
|
1 |
|
||||||||||
при следующих априорных данных |
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
xˆ |
|
7 |
|
|
; Px0 |
1 |
0 |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
(0) |
1,3 |
|
|
1 |
|
; |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
ˆ |
|
|
|
0 |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1; |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
(0) 1,1; P 0 |
|
|
|
|
|
|
|||||||
65
и |
|
|
|
|
|
|
|
|
|
|
|
1,5 |
0 |
|
|
|
|
|
|
|
C |
1,7 |
, D1 1, D2 1 . |
|
|
|
|
|
|
0 |
|
|
|
|
|
На |
рис. |
5.4. |
приведены |
графики |
переходных |
процессов |
||
x(k) x1(k), x2 (k) T , их оценок |
xˆ(k) xˆ1(k), xˆ2 (k) T |
и управлений |
||||||
u(k) , причем жирная линия изображает управляющие воздействия |
||||||||
u(k) , тонкая – состояния |
x(k) , а пунктирная – оценки состояния |
|||||||
xˆ(k) . |
|
|
|
|
|
|
|
|
|
5 |
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
5 0 |
|
20 |
40 |
60 |
80 |
100 |
|
Рис. 5.4. Графики переходных процессов, их оценок и управления
На рис. 5.5. приведен график оценок параметра и его точное значение (жирная линия).
Заметим, что в данном случае параметр не доступен оцениванию.
66
1.4
1.3
1.2
1.1
10 20 40 60 80 100
Рис. 5.5. Графики оценки параметра и еготочное значение
3. В модели объекта
x(k 1) Ax(k) Bu(k) Fq(k), |
x(0) x0 , |
|||||||
где |
|
|
|
|
|
|
|
|
|
1 |
0,1 |
|
|
|
|
0 |
|
A |
0,1a2,1 |
1 0,1a2,2 |
|
, |
B |
|
0,1b |
|
|
|
|
|
|
|
|
2 |
|
неизвестные |
параметры |
составляют |
вектор |
||||
1 , 2 T a2,1 ,a2,2 T |
и матрица |
A линейно зависит от компо- |
|||||
нент вектора следующим образом |
|
|
|||||
|
|
|
|
1 |
t |
|
|
|
|
A( ) |
t |
1 t( )2 |
. |
|
|
|
|
|
|
1 |
|
|
|
В |
математической |
модели |
измерительного |
комплекса |
|||
H 0 |
1 , т.е. измеряется только вторая компонента вектора со- |
||||||
стояния, r(k) |
последовательность гауссовских величин с нуле- |
||||||
вым средним и дисперсией R 0,35 . |
|
|
|||||
Адаптивное управление, формируемое на основе минимизации математического ожиданиялокального квадратичногофункционала
67
J (k) 12 M xT (k 1)C1x(k) uT (k)D1u(k) ,
имеет вид
|
|
|
u(k) (B |
|
CB D) |
|
B |
|
CA( (k))xˆ |
(k) , |
|
|
|
|
T |
|
1 |
|
T |
ˆ |
|
где |
xˆ |
ˆ |
оценки состояния и параметров, определяемые с |
|||||||
(k), (k) |
||||||||||
помощью фильтров Калмана (5.12),(5.13). Матрица (xˆ(k),u(k)) и
вектор f (xˆ(k),u(k)) |
в (5.13) для заданного вектора параметров |
|||||||||||||||||||||||||
определяются следующим образом: |
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
0 |
|
|
|
|
|
(xˆ(k))1 t(xˆ(k))2 |
|
|||
|
|
(xˆ(k),u(k)) t(xˆ(k)) |
|
|
t(xˆ(k)) |
|
|
, f (xˆ |
(k)) |
2 |
tb u(k) |
. |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
2 |
|
|
|
|
2 |
|
|||
|
|
Моделирование осуществлялось |
|
для |
модели, |
|
где матрицы |
|||||||||||||||||||
|
|
, |
|
|
|
и векторы x(0), |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
A |
B |
, F |
|
имеют вид |
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
0 |
1 |
|
|
|
|
|
0 |
|
|
|
0,51 |
0 |
|
|
|||
|
|
|
|
|
|
A |
0,4 |
1,4 |
, B |
|
, F |
0 |
|
|
, |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
1,3 |
|
|
|
0,55 |
|
||||||||
|
|
|
|
|
|
|
|
|
x(0) |
10 |
|
|
|
|
|
|
T |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
, |
( 0,4; 1,4) |
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
||||
при следующих априорных данных |
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
xˆ |
(0) |
|
7 |
|
|
Px0 |
|
|
1 |
0 |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
1,3 |
, |
|
|
, |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
1 |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
ˆ |
|
|
|
0,5 |
P 0 |
|
1 |
0 |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
(0) |
|
1,7 |
, |
|
, |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
1 |
|
|
|
|
||||
и
C1 |
1,5 |
0 |
|
, D D1 |
|
23,969 |
7,976 |
|
|
1,7 |
|
1, C |
7,976 |
8,766 |
. |
||
|
0 |
|
|
|
|
|||
На рис. |
5.6. |
приведены |
графики |
переходных |
процессов |
|||
x(k) x1(k), x2 (k) T , их оценок |
xˆ(k) xˆ1(k), xˆ2 (k) T и управлений |
|||||||
68
u(k) , причем жирная линия изображает управляющие воздействия u(k) , тонкая– состояния x(k) , апунктирная– оценкисостояния xˆ(k) .
10 |
|
|
|
|
|
5 |
|
|
|
|
|
0 |
|
|
|
|
|
0 |
|
|
|
|
|
50 |
20 |
40 |
60 |
80 |
100 |
|
|
|
i |
|
|
Рис. 5.6. Графики переходных процессов, их оценок и управления
На рис. 5.7. приведены в порядке следования графики оценок вектора параметров и их точных значений (жирная линия).
|
|
|
|
|
|
|
2 |
|
|
|
|
|
0 |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
0.5 |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
0 |
20 |
40 |
60 |
80 |
100 |
20 |
20 |
40 |
60 |
80 |
100 |
Рис. 5.7. Графики оценок вектора параметров и ихточные значения
69
6. Ограничения по управлению и состоянию
При синтезе управлений на величину управляющего воздействия обычно накладываются ограничения. Достаточно часто эти ограничения накладываются в форме неравенств:
|
|
|
U1 |
ui (k) U2 |
i |
1,m |
, |
(6.1) |
|
|
|
maxi |
|
maxi |
|
||
где |
U1 |
,U2 |
– заданные величины (например, |
для автопилота в |
||||
|
maxi |
maxi |
|
|
|
|
|
|
самолете U1maxi ,U2maxi определеныпоусловиюбезопасности полета). Тогда, учитывая ограничения, будем полагать:
|
U |
1maxi |
, u |
(k) U |
, |
|
|
|
|
i |
|
1maxi |
|
|
|
|
|
(k), U1maxi |
ui (k) U2maxi , |
|
|||
ui (k) ui |
(6.2) |
||||||
|
|
|
, ui (k) U2max |
|
|
||
|
U |
2max |
. |
|
|||
|
|
|
i |
|
|
i |
|
Ограничения, задаваемые для вектора состояния, чаще всего, определяют качество функционирования управляемого объекта. Если все компоненты вектора состояния удовлетворяют заданным ограничениям, то это говорит о том, что управление обеспечивает хорошее качество функционирования объекта.
70
