Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Адаптивные системы. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
886 Кб
Скачать

 

 

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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]