Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Лабораторная работа 5.pdf
Скачиваний:
37
Добавлен:
26.03.2016
Размер:
138.67 Кб
Скачать

Лабораторная работа № 5.

Алгоритмы адаптивного управления с ЭМ

Цель работы: изучение свойств системы, синтезировнной алгоритмами адаптивного управления с эталонной моделью.

5.1. Теоретические сведения

Задачу синтеза адаптивной системы управления с эталонной моделью (ЭМ) содержательно можно сформулировать следующим образом. Заданы уравнения объекта и эталонной модели. Требуется синтезировать алгоритм адаптивного управления, т.е. алгоритм управления (основного контура) и алгоритм адаптации, при которых система глобально устойчива и ошибка слежения – разность между выходными сигналами основного контура и эталонной модели – сходится к нулю при стремлении времени в бесконечность.

Здесь предполагается, что эталонная модель задана, хотя она должна быть определена исходя из заданных требований к синтезируемой системе управления. Это связано с тем, что определение эталонной модели по заданным требованиям к системе управления является обычной задачей управления и не связано спецификой адаптивного управления.

Параметрическая сходимость. При адаптивном управлении с ЭМ основное целевое условие – это обеспечение сходимости к нулю ошибки слежения e(t) = y(t) – ym(t). Если параметры регулятора принимают идеальные значения, то, естественно, это условие будет выполнено. Однако из сходимости к нулю ошибки слежения не следует параметрическая сходимость – сходимость варьируемых параметров к идеальным значениям.

1

Параметрическая сходимость зависит от структуры («сложности») задающего воздействия. Если задающее воздействие простое, например константа, то по окончании процесса адаптации варьируемые параметры в зависимости от начальных условий могут принять различные значения. Однако когда задающее воздействие g(t) обладает таким свойством, что выполняется так называемое условие постоянного возбуждения, то сходимость к нулю ошибки слежения влечет за собой параметрическую сходимость.

Определение 5.1. Условие постоянного возбуждения n – векторного сигнала v(t) выполняется, если существуют положительные константы Т и α такие, что при любом t > 0

t +T

v(τ)vT (τ) αIn ,

t

где In – еденичая матрица порядка n.

Адаптивное управление по состоянию линейным объектом.

Постановка задачи. Пусть линейный объект описывается уравление

(n)

(n1)

 

a0 y + a1

y +... + an y = u ,

(5.1)

где у – выход, u – управление, ai (i = 0,1, ..., n) – неизвестные параметры; знак a0 известен. Эталонная модель задается уравнением

(n)

(n1)

 

y + a1

y

+... + an ym = β0 g(t).

(5.2)

m

m

 

 

Здесь ym – выход эталонной модели, ai (i = 1, 2, ..., n) и β0 – известные положительные постоянные, g(t) – задающее воздействие. Требуется

2

определить алгоритм адаптивного управления, при котором система глобально устойчива и ошибка слежения

e(t) = y(t) ym (t) 0 при t → ∞.

Ниже при записи решения используется (n × 1)-матрица

B = (0 0 ... 0 1)T

(5.3)

и (n × n)-матрица Р, которая является решением уравнения Ляпунова

PA + AT P = −Q ,

где Q – положительно определенная матрица, А – (n × n)-матрица

 

0

 

1

0

...

0

 

 

 

 

0

 

0

1

...

0

 

 

 

 

 

 

,

(5.4)

A = ...

 

...

...

... ...

 

 

0

 

0

0

...

1

 

 

 

 

 

 

 

 

a

n

a

n1

a

n2

...

a

 

 

 

 

 

 

 

 

1

 

 

 

в которой элементами последней строки являются коэффициенты уравнения эталонной модели.

Утверждение 5.1. Алгоритмом адаптивного управления с ЭМ (5.2) линейным объектом (5.1), обеспечивающим глобальную устойчивость и сходимость ошибки e(t) = y(t) ym (t) к нулю при t → ∞, является

 

g

 

(n1)

 

 

(5.5a)

u = k

0

(t) + k

y

+... + k

n

y = kT v ,

 

 

1

 

 

 

 

 

 

 

= −sign(a0 )ΓvBT Px ,

(5.5б)

 

 

k

3

где k = (k

0

k ...

k

n

)T

(n

+

l)-вeктop

варьируемых

параметров

 

1

 

 

 

 

 

 

 

 

 

 

регулятора,

v = (g

(n1)

(n2)

y)

– (n

+

1)-вектор

сигналов, Г-

 

y

 

y

...

произвольная

положительно

определенная (n

+

l) × (n

+

l)-матрица,

x = (e e ...

(n1)

 

 

 

 

 

 

 

 

 

 

 

e )T -вектор состояния.

 

 

 

 

 

Если в качестве Q принимается матрица qIn(q >0, In – единичная матрица n-го порядка), то, не нарушая общности, можно при записи уравнения Ляпунова принять q = l, т.е. рассмотреть уравнение

PA + AT P = −In ,

а значение q учесть при выборе матрицы Г.

Уравнение синтезируемой системы (5.1), (5.5а) совпадает с уравнением эталонной модели, когда

k

0

= k* = a

β

0

,

k

i

= k* = a a α

i

(5.6)

 

0

0

 

 

 

i

i

0

 

Коэффициенты ki* (i = 0, 1, .... n) по определению являются идеальными. И если для каких-либо коэффициентов алгоритма управления (5.5а) априорная информация позволяет вычислить их идеальные значения, то, естественно, нет необходимости для их определения использовать алгоритм адаптации

(5.5б).

Пример 5.1. Объект описывается уравнением

y +3y + a2 y = u ,

где a2 – неизвестный параметр; уравнение эталонной модели имеет вид

y + 2ym + ym = g(t) .

4

Определить алгоритм адаптивного управления, обеспечивающего ограниченность всех переменных и сходимость ошибки е = у уm к нулю при t → ∞.

Решение. В данном случае α1 = 2 и α2 = 1, матрицы А и В имеют вид

A =

0

1

 

0

 

 

 

, B =

 

.

 

 

1

 

 

1

2

 

 

Уравнение Ляпунова

p

p

0

1

0

1 p

p

 

1

0

11

12

 

 

+

11

12

 

= −

 

p21

p22 1

2

1

2 p21

p22

0

1

после перемножения матриц принимает вид

p12

p11 2 p12

 

 

 

p21

 

p22

 

1

0

p

22

p

21

2 p

22

 

+ p

 

2 p

p

 

2 p

 

= − 0

1

 

 

 

 

 

11

 

21 12

 

22

 

 

Это уравнение, учитывая равенство р12 = р21 можно записать в виде системы

2 p12 = −1,

p11 2 p12 p22 = 0,

2 p12 4 p22 = −1.

Эта система имеет следующее решение:

p = 1 ,

p

22

= 1

,

 

p = 3 .

12

2

 

2

 

 

11

2

 

 

 

 

 

 

Поэтому матрица Р имеет вид

 

 

 

 

 

 

 

 

 

P =

3/ 2

1/ 2

=

1

3

1

 

 

 

 

2

 

.

 

1/ 2

1/ 2

 

 

1

1

5

Так как известны коэффициенты a0 = 1, a1 = 3, то по формуле (5.6) находим

 

 

k0 = k0* =1, ki = ki* = 3 2 =1.

 

 

 

 

 

В данном случае v = (g

y

 

y)T

и x = (e

e)T

(е = у

уm). Приняв Γ =γI3 по

формуле с (5.5б) получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g

 

 

 

 

 

 

 

 

 

g

 

 

k0

 

 

 

 

3 1 e

 

 

 

 

 

 

 

= −

γ

 

y

 

 

= −γ(e +e)

 

y

 

.

k

 

 

(0 1)

 

 

 

 

 

1

 

 

2

 

 

1 1 e

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k2

 

 

y

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда и в соответствии с соотношением (5.5а) для адаптивного алгоритма управления находим

Адаптивное управление по выходу линейным объектом с единичным относительным порядком. Выше был рассмотрен случай,

когда в уравнение объекта не входили производные управления или, что то же, когда относительный порядок передаточной функции объекта был равен ее порядку, и все фазовые координаты были доступны измерению. Однако обычно не все фазовые координаты доступны измерению. И чтобы получить их, нужно дифференцировать выходную переменную, что не желательно изза помех, которые при этом возникают.

Ниже рассматривается адаптивное управление по выходу, т.е. такое управление, при котором в алгоритмах управления и адаптации используются только входной и выходной сигналы объекта и сигналы, получаемые путем их фильтрации.

Постановка задачи. Пусть задан объект с передаточной функцией

 

 

P

( p)

 

 

pn1 +b pn2

+... +b

 

W ( p) = k

 

0

 

= k

 

1

 

n1

(5.7)

 

R

( p)

 

pn + a pn1

+... + a

 

0

0

 

0

n

 

 

 

0

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

6

и выбрана эталонная модель с передаточной функцией

 

 

P

( p)

 

 

pn1 + β pn2

+... + β

n1

 

 

W ( p) = k

 

m

 

= k

 

1

 

 

;

(5.8)

m R

( p)

 

 

pn1

+... +α

 

m

 

m

pn +α

1

n

 

 

 

 

m

 

 

 

 

 

 

 

 

k0, bi (i = 1, 2, ..., n – 1), ak(k = 1, 2, ..., n) – неизвестные параметры объекта,

знак k0 известен; Рm(р), Rm(p) – устойчивые полиномы, и передаточная функция Wm(p) является строго вещественно-положительной, т. е. она устойчива и Re W() > 0 при все ω ≥ 0.

Требуется определить алгоритм адаптивного управления, при котором система глобально устойчива и ошибка e(t) = y(t) – ym(t) стремится к нулю при t → ∞. При этом в алгоритмах управления и адаптации должны быть использованы только доступные измерению сигналы (задающее воздействие, входной и выходной сигналы объекта) и сигналы, которые получаются путем их фильтрации, т.е. сигналы на выходе фильтров, на вход которых подаются указанные сигналы. Принимается, что уравнения фильтров в нормальной форме Коши имеют вид

 

 

 

 

 

 

v = Ev + Fu ,

 

 

 

(5.9a)

 

 

 

 

 

 

z = Ez + Fy .

 

 

 

(5.9б)

Здесь v = (v

v

2

... v

n1

)T

– (n – 1)-вектор переменных, получаемых путем

1

 

 

 

 

 

 

 

 

 

фильтрации входного сигнала (управления) объекта; z = (z

z

2

... z

n1

)T

 

 

 

 

 

 

1

 

 

 

(n – 1)-вектор переменных, получаемых путем фильтрации выходного сигнала объекта; Е — (n – 1) × (n – 1)-матрица, F — (n – 1) × 1 — матрица, и они имеют следующий вид:

7

 

0

1

0

...

0

 

 

0

 

 

0

0

1

...

0

 

 

0

 

 

 

 

 

 

...

...

... ...

 

F =

... .

E = ...

,

 

0

0

0

...

1

 

 

0

 

 

 

 

 

β

n1

β

n2

β

n3

...

β

 

 

1

 

 

 

 

 

 

1

 

 

 

 

Здесь в последней строке матрицы Е стоят коэффициенты полинома числителя передаточной функции эталонной модели.

Утверждение 5.2. Алгоритмом адаптивного управления с ЭМ (5.8) линейным объектом (5.7), обеспечивающие глобальную устойчивость и сходимость ошибки e(t) = y(t) – ym(t)к нулю при t → ∞, является

u = kTv v +kTv z + ky y + kg g ,

 

 

kv

= −sign(k0 )γve,

 

 

 

 

 

 

 

 

 

k z

= −sign(k0 )γze ,

 

 

 

 

 

 

 

 

 

ky

= −sign(k0 )γye,

 

 

 

 

 

 

 

 

 

kg

= −sign(k0 )γge ,

 

 

 

 

 

 

 

где kv = (kv1 kv2 ...

kvn1)T

и

k z = (kz1

kz2

... kzn1)T

векторы

варьируемых параметров регулятора, ky

и kg — скалярные варьируемые

параметры регулятора

v = (v

v

2

... v

n1

)T

и

z = (z

z

2

... z

n1

)T

 

1

 

 

 

 

1

 

 

 

 

выходы фильтров (5.9a) и (5.9б) соответственно.

Пример 5.2. Пусть объект и эталонная модель задаются передаточными функциями

8

W0 = k0

p +b1

 

, Wm = km

p +1

,

p2 + a p + a

2

p2 +3p + 2

1

 

 

 

где k0, b1, a1, a2 – неизвестные параметры, известен знак k0: k0 > 0. Требуется определить алгоритм адаптивного управления, обеспечивающего глобальную устойчивость системы и сходимость к нулю разности между выходами системы и эталонной модели.

Решение. Передаточная функция эталонной модели является строго вещественно-положительной, так как она устойчива, и вещественная часть частотной передаточной функции при любой частоте ω ≥ 0 положительна:

ReW ( jω) = k

 

 

2(1+ω2 )

> 0 .

m (2

ω2 )2

+9ω2

m

 

В данном случае n = 2, β1 = 1, и уравнения фильтров (5.9) принимают вид

v = −v +u ,

z = −z + y .

Из утверждения 5.2 для адаптивного алгоритма управления получаем

u = kvv + kz z + ky y + kg g ,

kv = −γve ,

kz = −γze ,

ky = −γye,

kg = −γge.

9