Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Бобцов Адаптивные системы.pdf
Скачиваний:
120
Добавлен:
14.04.2015
Размер:
3.2 Mб
Скачать

3.1. Компенсация возмущения для линейного объекта 77

Подставляя (3.60) в (3.7) с учетом леммы 3.1, получим

( ) =

[ ( )](~( ) − 0 0

 

( +

 

 

 

))

( + ) − =1

 

 

( )

 

1

 

 

1

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ 0( ) +

( ) + ( )

 

 

 

 

 

 

 

 

=1

 

 

 

 

 

 

 

 

 

 

 

 

= ( ) + 2( ),

 

[ ( ) ]

 

 

 

 

 

 

 

 

где функция ( )

=

~( ) ограничена и стремится к нулю,

 

 

 

 

 

 

( )

 

 

 

 

 

 

 

 

 

2( ) — экспоненциально затухающая функция. Теорема 3.1 доказана.

Таким образом, получен реализуемый закон управления (3.43)– (3.47), обеспечивающий основную цель (3.4) lim →∞ ( ) = 0.

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

3.1.4Числовой пример

Для демонстрации эффективности предлагаемого алгоритма управления представим результаты математического моделирования неминемально-фазового объекта (3.1)-(3.3) с относительной степенью − = 2 и параметрами

=

−35 0 1

0

, =

 

0 ,

=

0

, 0

=

1

,

 

 

−10

1

0

0

 

 

 

0

 

 

 

1

 

 

 

 

1

 

 

 

 

 

24 0 0

0

 

13

 

0

 

 

0

 

 

 

−50

0

0

1

 

 

−5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

10

, 1

=

4

=

6

. (3.61)

 

 

 

 

 

 

 

 

 

5

 

 

2

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

4

 

1

 

 

3

 

 

 

 

 

 

 

 

 

 

10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

78 Глава 3. Устойчивые объекты с запаздыванием

На рис. 3.1 представлены переходные процессы в замкнутой системе для частоты возмущающего воздействия = 1, 5, параметров

закона управления = 10, = 2 и запаздывания = 6. На рис. 3.2а,

3.2б представлены переходные процессы в замкнутой системе для запаздывания = 5, частоты возмущающего воздействия = 3 и

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

= 0, частоты возмущающего воздействия = 3 и параметров алгоритма управления = 10, = 3. На рис. 3.2в, 3.2г представле-

ны переходные процессы в замкнутой системе для возмущающего воздействия с нерегулярной составляющей 1.

Для иллюстрации эффективности разработанного алгоритма

компенсации возмущений рассмотрим случай с двумя гармониками. Выберем параметры возмущения (3.3) 2 = [1 2 3 − 1] и

2 = [2 4 − 1 5].

На рис. 3.3 представлены результаты моделирования замкнутой системы для частот 1 = 2 и 2 = 3, запаздывания = 4, пара-

метров закона управления 1 = 20, 2 = 20, = 2. На рис. 3.4

представлен результат моделирования мультигармонического возмущения, частоты которого изменились скачкообразно с течением времени.

Из рис. 3.1-3.4 видно, что цель управления (3.4) достигается при различных значениях запаздывания, даже при его значительном увеличении. Отсутствие временной задержки является частным случаем рассматриваемой задачи, где = 0. С увеличением

значений параметров закона управления и идентификация ча-

стоты и, как следствие, компенсация возмущения происходит быстрее. Показаны адаптивные и робастные свойства алгоритма компенсации возмущения.

1Нерегулярная составляющая моделировалась с помощью формирующего фильтра ( ) = 1

10 +1 , на вход которого подавалось белошумное воздействие (с ограниченной частотой) мощностью = 0, 1. Матрица входов для этого

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

3.1. Компенсация возмущения для линейного объекта 79

ωˆ

 

 

 

 

1

 

 

 

 

 

2,5

 

 

 

 

 

 

 

 

 

 

0,5

2

 

 

 

 

 

1,5

 

 

 

 

0

1

 

 

 

 

 

 

 

 

 

 

−0,5

0,5

 

 

 

 

 

 

 

 

 

 

t, c

0

10

20

30

40

−1

0

50

 

 

 

 

 

a

 

 

 

 

 

b

 

 

 

 

 

t, c

0

10

20

30

40

50

(а) Временная диаграмма оценки ча-

(б) Выходная переменная ( ) без

стоты ^( )

 

 

 

 

 

управления (a) и с управлением (b)

1

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

−0,5

 

 

 

 

 

 

 

 

 

ϕˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

ϕˆξ

0,8

 

 

 

 

 

ˆ

−1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Lξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−1,5

 

 

 

 

 

 

 

 

 

 

0,6

 

 

 

 

 

 

−2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,4

 

 

 

 

 

 

−2,5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−3

 

 

 

 

 

 

 

 

 

 

0,2

 

 

 

 

 

 

−3,5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t, c

−4

 

 

 

 

 

 

 

 

 

t, c

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10

 

20

30

40

50

−4,5

10

 

20

 

 

 

 

30

40

50

0

 

0

 

 

 

 

 

(в)

Временная

диаграмма

оценки

(г)

Временная

диаграмма

оценки

функций ^

 

и

^

 

 

функций ^( ) и ^

 

( )

 

 

 

( )

 

( )

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

25

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kp

 

 

 

 

 

 

 

20

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kd

0,5

 

 

 

 

 

 

15

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

−0,5

 

 

 

 

 

 

−5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t, c

−10

 

 

 

 

 

 

 

 

 

 

−1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t, c

10

 

20

30

40

50

−15

10

 

20

 

 

 

 

30

40

50

0

 

0

 

 

 

 

 

(д)

Временная

диаграмма

оценки

(е)

Временная

диаграмма

оценки

функции ^

( )

 

 

 

 

функций

 

( ) и

 

( )

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 3.1. Временные диаграммы оценки частоты возмущения =

1, 5, выходной переменной, параметров закона управления при за-

паздывании = 6 и коэффициентах идентификатора = 2, = 10

80

 

 

Глава 3. Устойчивые объекты с запаздыванием

ωˆ

 

 

 

 

 

 

2

 

 

 

 

 

4

 

 

 

k = 1, λ = 1

 

1

 

 

 

 

 

 

 

 

k = 10, λ = 5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

−1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

U = 0

 

 

 

 

 

 

 

−2

 

 

 

k = 1, λ = 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t, c

 

 

 

 

k = 10, λ = 5

t, c

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

10

15

20

25

30

−3

 

 

 

 

 

0

10

20

30

40

50

 

 

 

 

 

 

 

0

(а) Временная диаграмма оценки ча-

(б) Выходная переменная ( ) при

стоты возмущения = 3

 

 

запаздывании = 4

 

 

 

ωˆ

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

0,5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−0,5

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

h = 0

 

 

 

 

 

 

 

−1,5

 

 

 

h = 6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t, c

 

 

 

 

h = 15

t, c

0

 

 

 

 

 

−2

 

 

 

 

 

10

20

 

30

40

10

20

 

30

40

0

 

 

0

 

(в) Оценки частоты ^( ) при различ-

(г) Выходная переменная ( ) при

ном запаздывании

 

 

 

различном запаздывании

 

ωˆ

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

0,5

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

1,5

 

 

 

 

 

 

0

 

 

 

 

 

1

 

 

 

 

 

 

−0,5

 

 

 

 

 

0,5

 

 

 

 

 

 

−1

 

 

 

 

 

 

 

 

 

 

 

t, c

 

 

 

 

 

t, c

0

 

 

 

 

 

−1,5

 

 

 

 

 

 

20

40

60

80

100

20

40

60

80

100

0

 

0

(д) Временная диаграмма оценки ча-

(е) Выходная переменная ( ) без

стоты возмущения ^( ) с аддитивной

управления (a) и с управлением (b)

нерегулярной составляющей

при запаздывании = 4 и коэффици-

 

ентах идентификатора = 1 и = 10

Рис. 3.2. Временные диаграммы оценки частоты и выходной переменной ( )

3.1. Компенсация возмущения для линейного объекта 81

4

 

 

 

 

3,5

 

 

 

ωˆ1

 

 

 

ωˆ2

 

 

 

 

3

 

 

 

 

2,5

 

 

 

 

2

 

 

 

 

1,5

 

 

 

 

1

 

 

 

 

0,5

 

 

 

 

 

 

 

 

t, c

0

50

100

150

200

0

(а) Временная диаграмма оценки частот ^ ( )

1,2

 

L1

1

L2

 

Lξ1

0,8

Lξ2

 

0,6

 

0,4

0,2

t, c

0

0

50

100

150

200

(в) Временная диаграмма оценки

функций ^ ^

( ) и ( )

ˆ

ξ0

0,5

0

−0,5

t, c

−1

0

50

100

150

200

(д) Временная диаграмма оценки

функции ^

0( )

1,5

a

1

b

0,5

0

−0,5

−1

t, c

−1,5

0

50

100

150

200

(б) Выходная переменная ( ) без управления (a) и с управлением (b)

1

 

 

ϕˆ1

0

ϕˆ2

 

ϕˆξ1

−1

ϕˆξ2

 

−2

 

−3

−4

t, c

−5

0

50

100

150

200

(г) Временная диаграмма оценки функций ^ ( ) и ^ ( )

10

0

 

 

 

 

−10

 

 

 

 

−20

 

 

kp1

 

 

 

 

 

 

 

 

kp2

 

−30

 

 

kd1

 

 

 

 

 

 

 

 

kd2

t, c

−40

50

100

150

200

0

(е) Временная диаграмма оценки функций ( ) и ( )

Рис. 3.3. Временные диаграммы оценки частот возмущения 1 = 1,2 = 2, выходной переменной ( ) и параметров закона управления при запаздывании = 6 и коэффициентах идентификатора = 2,

= 10