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

Герасимов Адаптивное управление

.pdf
Скачиваний:
117
Добавлен:
14.04.2015
Размер:
364.62 Кб
Скачать

где =δ 2 / 2λ +θ2 / 4γ постоянная величина. Решая полученное дифференциальное неравенство, получаем:

V (ε(t)) eλtV (0) + λ λ eλt ,

откуда с учетом (2.5) следует, что

1

 

2

 

λt

 

 

 

 

 

 

λt

 

 

 

 

 

 

 

λ e

 

 

 

2

ε

 

e

V (0) + λ

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε(t)

 

λt

+

 

 

 

e

λt

(2.6)

2 e

V (0)

λ

λ

.

 

 

 

 

 

 

 

 

 

 

 

 

Из последнего неравенства следует экспоненциальная сходимость ошибки управления ε

к ограниченному

множеству с

границей

=

2

/ λ . При этом величину

можно

уменьшить путем

увеличения

коэффициентов

λ

и γ . Как следствие, величина θˆ

становится ограниченной.

Таким образом, алгоритм управления (1.6), (2.4) обеспечивает устойчивость в замкнутой системе и является робастным по отношению к внешнему возмущению. В то же время этот алгоритм имеет следующие недостатки, а именно:

при отсутсвии возмущения установившаяся ошибка ε(t) отлична от нуля, что видно из неравенства (2.6);

если подставить (2.4) в (1.6), то нетрудно видеть, что управление

пропорционально величине x2 . Следовательно при росте x амплитуда управления возрастает квадратично, в связи с чем практическая применимость такого закона (1.6) имеет существенные ограничения.

Рассмотрим решение, лишенное недостатков алгоритмов (1.6), (1.9) и (1.6), (2.4) за счет наделения нового алгоритма управления адаптивными и робастными свойствами.

Решение №3

Рассмотрим совместно с настраиваемым регулятором (1.6) алгоритм адаптации, параметрический дрейф в котором ограничивается линейной стационарной обратной связью:

&

(2.7)

θˆ = −σθˆ γxε ,

где σ постоянная положительная величина.

Проведем анализ устойчивости замкнутой системы, представленной объектом (2.1), регулятором (1.6) и алгоритмов адаптации (2.7) с помощью функции Ляпунова (1.8). Возьмем производную от функции и проведем ряд преобразований:

11

 

 

 

 

1

 

 

 

 

 

1 ~~&

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

1

 

~ &

 

 

 

 

 

 

 

 

2

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

θ~

 

 

 

ˆ

 

&

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V = 2

2εε

+

2γ

2θθ

= ε(λε

θ x δ )

γ θθ

 

= −λε

θ xε δε γ (

σθ γxε)=

 

 

 

 

 

 

 

 

&

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= −λε2 δε +

 

σ

θ~θˆ

= −λε2 δε

+

 

σ

θ~(−

θ~ +θ)= −λε2 δε σ

θ~2 +

σ

θ~θ =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ

 

 

 

 

γ

 

 

 

 

 

 

λ

2

 

λ

 

2

 

 

 

 

 

1

 

2

 

σ

~2

 

 

σ

 

~2

 

 

 

σ

 

~

 

 

 

 

 

 

 

σ

θ

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= −

2

ε

 

2 ε

 

 

 

δε ± 2λ

δ

 

2γ

θ

 

 

 

2γ θ

 

 

+

 

γ

θθ ±

 

2γ

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

= −

λ

ε

2

σ

 

 

~2

 

λ

ε +

 

 

1

 

 

 

 

 

2

 

σ

 

 

1 ~

+

 

 

 

1

θ

 

2

 

 

 

 

 

1

 

δ

 

2

+

σ

θ

2

 

 

 

 

 

 

 

θ

 

 

 

 

 

 

δ

 

 

 

 

 

 

 

 

 

 

θ

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2γ

 

 

 

 

 

 

2

 

 

 

2λ

 

 

 

 

 

 

 

 

γ

 

 

2

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

2λ

 

 

 

 

 

 

2γ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ

2

 

 

 

σ

~2

 

 

1

 

2

 

 

 

 

σ

 

 

2

 

 

 

λ

 

2

 

 

 

 

σ

~2

 

 

 

 

 

1

 

 

 

 

2

 

 

σ

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

≤ − 2 ε

 

 

 

 

θ

+

 

δ

 

+

 

 

 

 

 

θ

 

 

≤ −

 

2 ε

 

 

 

 

θ

 

+

 

 

δ

 

 

+

 

 

 

θ

 

.

 

 

 

 

 

 

 

 

 

2γ

2λ

 

 

2γ

 

 

 

 

2γ

 

 

 

2λ

 

 

 

2γ

 

 

Введем обозначение κ = min{λ, σ}. Тогда, считая λ, σ положительными, имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

2

 

 

 

 

 

1

 

 

~

2

 

 

 

 

1

 

 

 

 

 

2

 

 

 

 

 

 

σ

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V& ≤ −κ

 

ε

 

 

 

 

 

 

 

 

θ

 

+

 

 

 

δ

 

 

 

+

 

 

 

 

 

θ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

2γ

2λ

 

 

 

 

2γ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 / 2λ +σθ2 / 2γ

V& ≤ −κV +

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

постоянная

 

величина. Далее, применяя результат,

 

=δ

 

 

полученный в ходе второго решения задачи, окончательно получаем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V (ε(t)) eκtV (0) + κ κ eκt

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε(t)

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 eκtV (0)

 

κ

κ

eκt .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из последнего неравенства следует экспоненциальная сходимость ошибки

управления к ограниченному множеству с границей

 

=

 

2

 

/ κ .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Алгоритм

управления

 

(1.6), основанный

 

на

 

алгоритме

 

адаптации

(2.7), также

обеспечивает устойчивость замкнутой системы и является робастным по отношению к внешнему возмущению. В то же время алгоритм (1.6), (2.7) позволяет парировать недостатки робастного алгоритма управления (1.6), (2.4). Так, при отсутствии внешнего

возмущения верхняя граница может быть снижена до нуля за счет обнуления

коэффициента σ (т.н. гибридная σ -модификация). Кроме того, для уменьшения

нет

необходимости в значительном увеличении γ , которое влечет за собой рост амплитуды

управляющего воздействия. Снижение

можно обеспечить путем уменьшения

коэффициента σ .

 

Порядок выполнения работы

 

12

1. На основе данных, приведенных в таблице 1.1, провести моделирование адаптивной системы управления, полученной в лабораторной работе №1, в условиях действия на объект возмущения вида:

δ(t) = 0.25sin10t +0.5 .

По результатам моделирования построить три графика. На первом вывести x и xm , на втором u , на третьем θ~ =θ θˆ .

2. Заменить алгоритм адаптации (1.9) на нелинейную обратную связь (2.4) и повторить эксперимент для трех разных значений коэффициента γ . Определить влияние

этого коэффициента на свойства системы. По результатам моделирования для каждого γ построить два графика. На первом вывести x и xm , на втором u .

3. В п.1 заменить алгоритм адаптации (1.9) на робастную σ -модификацию (2.7). повторить эксперимент для трех разных значений коэффициента σ . Определить влияние этого коэффициента на свойства системы. По результатам моделирования для каждого σ построить два графика. На первом вывести x и xm , на втором u .

13

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

Цель работы: освоение принципов построения адаптивной системы управления многомерным объектом.

Методические рекомендации. До начала работы студенты должны ознакомится с принципом построения алгоритмов адаптации на основе стандартной модели ошибки с измеряемым состоянием [Ник Уш]. Лабораторная работа рассчитана на 2 часа.

Теоретические сведения. Рассмотрим пример типовой задачи адаптивного управления многомерным объектом с использованием эталонной модели. При этом воспользуемся принципами решения аналогичной задачи для объекта первого порядка (см. лабораторную работу №1).

Постановка задачи

Дан объект управления, модель которого описывается уравнениями вида

x = Ax +b u , x(0)

(3.1)

&

 

y = C x ,

(3.2)

где x Rn вектор состояния, y регулируемая переменная,

 

 

0

 

1

 

 

0

 

0

 

 

 

A =

 

M

 

M

 

 

 

 

0

 

0

 

a

0

a

 

 

 

1

0

 

L

0

 

 

0

 

 

 

1

 

L

0

 

 

 

0

 

 

 

 

 

 

 

 

, C = [1 0 L

0 0],

O

O M

 

, b =

 

M

 

 

 

 

0

 

L

1

 

 

 

0

 

 

 

a

2

L − a

n1

 

 

b

 

 

 

 

 

 

 

 

 

0

 

 

 

ai , i = 0, n 1 неизвестные параметры, b0 известный параметр.

Задача управления заключается в компенсации параметрической неопределенности объекта и обеспечении следующего целевого равенства:

lim(xM (t) x(t))= lim e(t) = 0 ,

(3.3)

t→∞

t→∞

 

где e = xM x вектор ошибки управления,

xM Rn вектор, генерируемый

 

эталонной моделью:

 

 

xM = AM xM +bM g ,

(3.4)

&

 

 

yM = CM xM ,

 

(3.5)

где g(t) задающее воздействие,

14

 

 

0

1

0

L

0

 

 

 

 

0

0

1

L

0

 

 

 

 

 

 

A

=

M

M

 

M

O

 

M

 

, b

M

 

0

0

0

L

1

 

M

 

 

 

 

 

a

M 0

a

M 1

a

M 2

L − a

M n1

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

0

 

 

 

 

, CM = [1 0 L 0 0],

=

M

 

 

0

 

 

 

 

 

a

 

 

 

 

M 0

 

Коэффициенты эталонной модели aM i ,i =1, n 1 строятся на основе метода

стандартных характеристических полиномов [Боб, Мир + Григ Лук].

Решение задачи

Предполагая параметры объекта известными, синтезируем регулятор, который обеспечит условие (3.3) с заданными динамическими показателями качества временем переходного процесса tп и перерегулированием σ .

Для синтеза регулятора сформируем ошибку слежения ε = xM x , возьмем ее производную и проведем ряд алгебраических преобразований:

 

 

 

 

0

1

0

 

 

 

 

0

0

1

e

= xM x =

 

 

M

M

 

M

&

&

&

 

 

 

 

 

0

0

0

 

 

 

 

 

 

 

a

M 0

a

M 1

a

M 2

 

 

 

 

 

 

 

L

0

 

x

M 1

 

L

0

 

 

 

 

xM 2

 

O

M

 

M

 

L

1

xM n1

 

 

 

 

xM n

 

 

 

 

 

L − aM n1

 

 

 

 

0

 

 

0

 

 

 

+

M

g

 

0

 

 

 

a

 

 

 

M 0

0

 

1

 

 

0 L 0

 

x

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0

 

 

1

 

L 0

 

 

 

 

1

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

 

M

 

 

M

 

O

M

 

 

 

M

 

M

u ±

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0

 

 

0 L 1

 

x

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a1

 

a2

 

 

 

 

x

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a0

 

L − an1

 

 

 

 

b0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

1

 

 

0

L 0

 

x

 

0

 

 

 

1

 

 

 

0

 

 

L 0

 

e

 

 

 

0

 

 

0

 

 

1

L 0

 

 

 

 

x

1

 

 

0

 

 

 

0

 

 

 

1

 

 

L 0

 

 

1

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

e

 

±

M

 

 

M

 

O

O

 

 

M

 

 

 

 

 

 

 

 

=

M

 

 

 

M

 

 

 

M

 

 

O

 

 

M

 

 

2

 

+

 

 

 

 

 

 

 

 

 

 

M

 

 

 

 

 

 

 

 

 

 

 

 

M

 

 

0

 

 

0

 

 

0

L 1

 

 

 

x

n1

 

 

0

 

 

 

0

 

 

 

0

 

 

L 1

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n1

 

 

 

 

aM 1 aM 2

 

 

 

 

 

 

 

 

x

n

 

 

 

aM 1

aM 2

 

 

 

 

 

 

e

 

 

aM 0

L − aM n1

 

 

 

 

aM 0

 

L − aM n1

 

n

 

 

 

 

0

 

 

 

 

0

 

0

 

 

 

 

 

L

 

 

 

0

 

 

 

 

x

 

0

0

 

 

 

 

 

 

 

0

 

 

 

 

0

 

0

 

 

 

 

 

L

 

 

 

0

 

 

 

 

 

1

 

 

0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

M

 

 

 

 

M

 

 

M

 

 

 

 

 

O

 

 

 

M

 

 

 

 

 

M

 

M

u +

M

 

g.

 

 

 

 

 

 

 

0

 

 

 

 

0

 

0

 

 

 

 

 

L

 

 

 

0

 

 

 

x

 

 

 

0

 

 

0

 

 

 

 

 

 

 

 

 

 

+ a a + a a + a

 

 

L − a

 

 

+ a

 

 

 

 

n1

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

xn

 

 

0

 

a

 

 

 

 

 

 

 

 

 

 

M 0

 

0

 

 

M 1

1

 

M 2

 

 

 

2

 

 

 

 

 

 

M n1

 

n1

 

 

 

 

 

 

 

 

M

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

Полученный результат запишем в более компактной форме:

e& = AM e + k[(aM 0 + a0 )x1 + (aM 1 + a1 )x2 +... + (aM n1 + an1 )xn + aM 0 g b0u], (3.6)

15

где k = [0 0

 

L 0 1]T . Введем обозначения:

 

 

 

θ1 = (aM 0 + a0 ),

θ2 = (aM 1 + a1 ), ..., θn = (aM n1 + an1 ),

после чего перепишем (3.6) в следующем виде:

 

 

 

 

 

 

&

AM e + k[θ

T

x + aM 0 g b0u],

(3.7)

 

 

 

 

 

e =

 

где θT = [θ

,θ

2

,...,θ

n

] вектор, содержащий параметры матрицы

A .

1

 

 

 

 

 

 

 

 

M

Нетрудно видеть, что, если

 

 

 

 

 

 

 

 

u =

 

1

[θT x + aM 0 g],

(3.8)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

0

 

 

 

 

то ошибка слежения подчиняется закону вида

 

 

 

 

 

 

e = AM e

 

 

(3.9)

 

 

 

 

 

&

 

 

 

 

 

 

и экспоненциально затухает со временем.

Однако в исходной постановке задачи параметры матрицы AM неизвестны. Следовательно закон (3.8) физически нереализуем. Заменим в 3.7 неизвестные параметры θ на оценки θˆ и получим настраиваемый закон управления:

u =

1

[θˆT x + aM 0 g].

(3.10)

b

 

0

 

 

 

Подставим последнее выражение в (3.7) и получим модель ошибок:

 

&

 

~T

x ,

(3.11)

e =

AM e + kθ

где θ~ =θ θˆ вектор параметрических ошибок.

Аналогично решению, преведенному в лабораторной работе №1, выбирается функция Ляпунова:

V = eT Pe + 21γ θ~Tθ~ ,

где P положительно определенная симметричная матрица, удовлетворяющая уравнению Ляпунова

AT

P + PA

= −Q ,

(3.12)

M

M

 

 

где Q произвольно выбранная симметричная положительно определенная матрица.

Далее берется производная от функции Ляпунова, проводятся простейшие алгебраические преобразования и выводится алгоритм адаптации аналогичный алгоритму

(1.9):

&

θˆ(0) = 0

(3.13)

θˆ = γxkT Pe ,

 

16

 

где γ > 0 алгоритм адаптации.

Таким образом, алгоритм адаптивного управления состоит из настраиваемого регулятора (3.10), алгоритма адаптации (3.13), в котором матрица P находится из (3.12).

Адаптивный регулятор (3.10), (3.13) обеспечивает:

1.Устойчивость замкнутой системы;

2.Асимптотическое стремление ошибки e к нулю;

3.Ограниченность сигнала θˆ . При этом вектор θˆ асимптотически стремится к θ ,

если вектор x удовлетворяет условию неисчезающего возбуждения:

t+Tx(τ)xT (τ)dτ >αI ,

(3.14)

t

 

где α > 0 , T > 0 постоянные величины.

Условие (3.14) эквивалентно условию наличия не менее k / 2 гармоник в векторе

x .

4. При увеличении коэффициента γ до определенного оптимального значения

увеличивается скорость сходимости параметрических ошибок θ~ к нулю. Дальнейшее увеличение γ приведет к снижению скорости сходимости ошибки и появлению

колебаний в системе.

Порядок выполнения работы

1. На основе заданных в таблице 3.1 значений времени переходного процесса tп и максимального перерегулирования σ сформировать эталонную модель в форме (3.4), (3.5). Построить график переходной функции модели, на котором показать tп и σ ;

2. На основе предположения, что параметры объекта известны построить и промоделировать систему управления с регулятором (3.8). Провести три эксперимента при следующих условиях:

использовать расчетные значения параметров объекта θ1 и θ2 ;

незначительно отклонить параметры объекта так, чтобы система не потеряла устойчивость;

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

По результатам каждого эксперимента построить траектории x(t) и xM (t) на одном графике и e(t) на другом.

3.Провести моделирование адаптивной системы управления с регулятором (3.10)

иалгоритмом адаптации (3.13). В ходе моделирования проиллюстрировать свойства 1-4 алгоритма управления. Для этого необходимо провести эксперименты при следующих условиях:

повторить три эксперимента п.п. 2 для фиксированного значения γ ;

используя расчетные значения параметров объекта, провести эксперимент с тремя различными значениями γ ;

провести один из пердыдущих экспериментов данного пункта при g(t) =1 .

17

По результатам каждого эксперимента построить траектории x(t) и xM (t) на

одном графике, e(t) на втором, θ~ =θ θˆ на третьем. 4. Сделать выводы по каждому пункту работы.

18

Таблица 3.1. Варианты заданий

 

 

 

Коэфф.

Время

Максимальное

 

 

Матрица

перерегули-

 

Вар.

передачи

переходного

Сигнал задания g(t)

A

рование

 

b0

процесса, tп

 

 

 

 

σ, %

 

 

 

 

 

 

 

1

0

1

1

0.16

0

2sign[cost] +3

 

 

 

1

1

 

 

 

 

2

0

1

2

0.3

0

cos t +3sin 2t +5

 

 

 

1

1

 

 

 

 

 

0

1

 

 

 

t

3

 

 

3

0.9

0

sign[cos 2τ]dτ

 

4

0

 

 

 

0

4

0

1

4

0.2

0

sign[sin 0.5t] + 2

 

 

 

1

4

 

 

 

 

5

0

1

5

0.6

0

10 cos 0.5t + 2sin t +12

 

 

 

2

5

 

 

 

 

 

0

1

 

 

 

t

6

6

0.3

0

sign[sin 0.5τ]dτ + 3

 

 

 

9

1

 

 

 

0

7

0

1

7

0.7

0

0.5sign[sin 0.7t] +1

 

 

 

2

3

 

 

 

 

8

0

1

8

0.1

0

5sin 0.5t + 4 cos 0.1t +8

 

 

 

2

2

 

 

 

 

9

0

1

9

0.9

0

sign[cost] +3

 

 

 

0

0

 

 

 

 

 

0

1

 

 

 

t

10

 

 

1

3.5

15

sign[sin 2τ]dτ

 

0

1

 

 

 

0

11

0

1

2

0.6

15

sin 5t + 0.5cos 0.2t + 2

 

 

 

1 6

 

 

 

 

12

0

1

3

0.9

15

2sign[sin 0.4t]+3

 

 

 

10

1

 

 

 

 

 

0

1

 

 

 

t

13

 

 

4

0.4

15

3sign[sin 0.5τ]dτ +3

 

1

5

 

 

 

0

14

0

1

5

0.2

15

2sin 0.2t +sin 0.1t +8

 

 

 

9

0

 

 

 

 

15

0

1

6

0.5

15

4sign[cos 2t]+5

 

 

 

9

1

 

 

 

 

 

0

1

 

 

 

t

16

 

 

7

0.5

15

3sign[cos 0.2τ]dτ +3

 

7

6

 

 

 

0

19

17

0

1

8

0.45

15

sign[sin 0.5t] + 2

 

 

 

1 7

 

 

 

 

18

0

1

9

0.15

15

0.8sin 2t + cos 0.8t + 2

 

 

 

10

6

 

 

 

 

 

0

1

 

 

 

t

19

1

0.7

15

sign[sin 0.3τ]dτ +1.5

 

 

 

2

2

 

 

 

0

20

0

1

2

1.2

0

7sign[cos 0.9t] +8

 

 

 

15

6

 

 

 

 

21

0

1

3

1.5

0

0.4sin 3t + cos 0.1t +1.5

 

 

 

11

1

 

 

 

 

 

0

1

 

 

 

t

22

4

0.8

0

6sign[sin 0.1τ]dτ +9

 

 

 

3

3

 

 

 

0

23

0

1

5

0.9

0

2sign[sin t] + 4

 

 

 

12

6

 

 

 

 

24

0

1

6

0.2

0

9 sin 0.2t +9 cos 0.1t +15

 

 

 

0

11

 

 

 

 

 

0

1

 

 

 

t

25

7

0.5

0

4sign[sin 6τ]dτ +5

 

 

 

13

2

 

 

 

0

26

0

1

8

0.9

0

4sign[cost] + 3

 

 

 

15

6

 

 

 

 

27

0

1

9

1.3

0

sin 0.1t + cos 5t + 2

 

 

 

13

1

 

 

 

 

 

0

1

 

 

 

t

28

1

1.6

0

9sign[sin 0.1τ]dτ +12

 

 

 

4

4

 

 

 

0

29

0

1

2

0.75

15

3sign[sin 4t]+8

 

 

 

9

7

 

 

 

 

30

0

1

3

0.65

15

7 sin 0.3t +8cos t + 20

 

 

 

1

9

 

 

 

 

20