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

u_lectures

.pdf
Скачиваний:
52
Добавлен:
22.03.2016
Размер:
2.16 Mб
Скачать

 

133

ˆ

- x&= (A - LC)(xˆ - x).

x&

Известно, что если собственные числа матрицы [A - LC], т.е. корни

характеристического уравнения системы ”объект-наблюдатель” лежат в левой полуплоскости, то система устойчива и при t 0 ошибка оценивания e= xˆ - x 0 , значит вектор xˆ x . Скорость сходимости оценок зависит от полюсов ”объекта и наблюдателя”.

Пример 1. Построить наблюдатель полного порядка для объекта второго порядка.

u

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

x1

 

 

 

a2

 

 

 

 

a

 

 

-

 

-

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

p

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 14.3. Структурная схема

Уравнения состояния, составленные для объекта на рис. 14.3, имеют матрицы А, В и С:

 

0

 

a

 

 

0

 

; C = [1 0].

A= -a

-a1

 

; B = a

 

 

 

2

2

 

2

 

 

Так как объект имеет второй порядок, то матрица L имеет вид

L

 

L = L11

.

 

21

 

В соответствии с методикой синтеза проведем вычисления:

0

 

a1

 

L11

 

0

 

a1

 

L11

0

-L11

a1

 

A - LC =

-a

2

-a

2

 

- L

×[1 0

]=

-a

2

-a

2

 

- L

0

 

=

-a

2

- L

-a

2

 

 

 

 

 

21

 

 

 

 

 

21

 

 

 

 

21

 

 

Найдем характеристическое уравнение системы:

p + L11

-a1

 

= (p + L11

)(p + a2

)+ a1(a2 + L21 )=

det[pI - (A - LC)]= a

 

+ L

p + a

 

 

2

21

2

 

 

 

 

=p2 +(a2 + L11 )p + L11a2 + a1a2 + a1L21 ;

Дж(p)= p2 + 2 ω0 p +ω0 2 .

134

Приравнивая коэффициенты при «р» в одинаковой степени найдем коэффициенты связей наблюдателя:

a2 + L11 = 2ω0

 

L11 = 2ω0 - a2 ;

 

 

 

;

 

1

 

 

 

 

 

L a

 

+ a a

 

+ a L = ω2

L21 =

ω02

- a1a2 - a2 (2ω0

- a2

)

;

2

2

 

a

11

1

1 21 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

Структурная схема объекта с наблюдателем приведена на рис. 14.4.

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2

 

 

 

 

 

x2

 

 

 

a

 

 

 

x1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L21

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xˆ

= x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xˆ2

 

 

 

 

a1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

3

 

 

 

-

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

xˆ 2= x4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 14.4. Структура наблюдателя с объектом

Связи наблюдателя необходимо вводить на входы интеграторов.

14.2. Редуцированный наблюдатель

Рассмотренный в предыдущей лекции наблюдатель называют наблюдателем полного порядка. Он оценивает весь вектор x, несмотря на то, что компоненты вектора “y”, входящие в состав “x”, могут быть измерены непосредственно. Для восстановления лишь тех переменных, которые не могут быть непосредственно измерены, наблюдатель может быть выполнен как редуцированный, или наблюдатель пониженного порядка. Все переменные состояния объекта, составляющие вектор “x”, можно разделить на измеряемые, которые образуют вектор “y”, размерности m, и неизмеряемые, образующие вектор w размерности (n-m), т.е. записать:

y

}m

x =

(14.8)

w

}n - m

135

Тогда уравнение x& = Ax + Bu может быть записано в виде:

y

A A

y B

&

= 11

12 ×

+ 1 ×u

 

w

A21

A22

w

B2

&

 

 

 

 

m n-m

или, что то же самое:

y& = A11 y + A12w+ B1u

w& = A21 y + A22w+ B2u

эти матрицы имеют размерности:

 

 

A11

-

m

×

m;

A22 - (n - m)

×

(n - m);

A12

-

m

×

(n - m);

A21 - (n - m)

×

m;

B1 - m × 1;

B2 - (n - m) × 1.

(14.9)

(14.10)

(14.11)

На основании второго уравнения (14.11) можно рассматривать часть системы с выходным вектором w, для которой входными воздействиями являются В2u и A21 y . Для этой части системы по принципам, изложенным

ранее, строится наблюдатель, на входе которого действуют векторы В2u и A21 y , а также вектор ошибки восстановления через некоторую матрицу L.

Матрица L в редуцированном наблюдателе играет ту же роль, что и в наблюдателе полного порядка. Вектор “w” неизмеряем. Однако он может быть измерен косвенно через вектор входного воздействия “u” и измеряемый вектор “y” в соответствии с первым уравнением системы

(14.10):

A12w = y& - A11 y - B1u.

Для получения в определенном масштабе вектора ошибки w~ надо умножить вектор восстановленных координат wˆ слева на матрицу – А12 и определить разность A12w - A12wˆ . После умножения на L ее следует ввести на

вход наблюдателя. Это поясняет структурная матричная схема (рис.4.4). Группируя входные каналы с матрицами В2 и 1 , A21 и LA11 , и перенося

сигнал “p y” со входа на выход наблюдающего устройства, что позволяет избежать операции дифференцирования, можно получить структурную схему на рис. 14.5, а затем, вынося точку суммирования за точку съема, схему на рис. 14.6, где обозначено:

136

G = A21 LA11 + ( A22 LA12 )L; F = A22 LA12 ;

Вектор “z” представляет собой вектор состояния наблюдателя. Он связан с вектором восстановленных неизмеряемых координат ωˆ и вектором “y” линейным преобразованием z = wˆ - L× y .

 

 

 

 

 

 

 

 

 

 

 

1

 

 

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

р

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A12 w

11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А21

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

A12 wˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

 

 

wˆ

 

 

 

 

 

В2

 

 

 

 

 

 

 

 

12

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А22

Рис. 14.5

y

 

 

 

 

 

 

А21-LA11

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

wˆ

 

B2-LB1

 

 

 

 

 

 

 

 

 

 

 

I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

A22-LA12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

G

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

wˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

I

 

 

 

 

 

 

 

 

 

TB

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

137

Учитывая, что в редуцированном наблюдателе полный вектор восстановленных координат формируется так

y xˆ = wˆ ,

можно записать, что z = Txˆ , где T = [ -L I ] . Тогда В2 - LB1 = TB .

Тем же образом, что и для наблюдателя полного порядка, решается вопрос выбора динамики или собственных значений матрицы A22 LA12 ,

через выбор матрицы L.

det [pI - (A22 LA12 )] = Дж( p).

Для схемы на рис. 14.4 уравнение наблюдателя:

ˆ

= (A22

- LA12 )w+(Aˆ 21 - LA11 )y + Ly&+(B2

- LB1 )u .

w&

Пример 2. Построить наблюдатель для восстановления переменной х2 в объекте, приведенном ранее.

x1 = 0 a1 × x1 + 0 u ;

y = [1 0]× x1 = x ;

&

 

 

 

 

 

 

 

&

 

-a2

 

 

x

1

x2

 

 

-a2 x2

a2

2

 

A22 - LA12 = [-a2 ] - [l11 ]×a1 = -(a2 +l11a1 );

 

 

det[pI - (A22 - LA12 )]= p + a2 +l11a1 = p + w0 = Дж(p).

 

ω = a +l a

l = ω0 - a2 ;

 

 

0

2

11

1

11

a1

 

 

 

 

 

 

 

 

 

A21 LA11 = −a2 l11 0 = −a2 ; B2 LB1 = a2 l11 0 = a2.

Схема наблюдателя приведена на рис. 14.8.

138

С помощью рассмотренных наблюдателей дополнительно можно решать и другие задачи, например.

1. Оценивать возмущения, приложенные к объекту (для этого необходимо в уравнения состояния объекта ввести возмущения как одну из координат).

2. Оценивать изменяющиеся параметры системы ki, Ti и на этой основе строить адаптивные системы.

u

 

 

 

 

 

x2

 

 

 

 

 

 

 

x1=y

 

 

a2

 

a1

 

 

 

 

-

-

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2

l11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

z

 

 

xˆ 2

a2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2 + l11a1

Рис. 14.8

139

Оценки координат можно использовать для построения модального регулятора.

140

15. ОСОБЕННОСТИ ДИНАМИКИ СИСТЕМЫ, ЗАМКНУТОЙ ЧЕРЕЗ НАБЛЮДАТЕЛЬ

Рассмотрим систему с модальным регулятором, где обратные связи замкнуты по переменным состояния, выработанным наблюдателем (рис. 15.1).

Система описывается уравнением:

x = Ax + Bu .

(15.1)

&

 

Наблюдателю соответствует уравнение:

ˆ

= (A - LC)xˆ + Bu + Ly .

x&

(15.2)

 

V

u

&

x

y

x I

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yˆ

 

 

 

 

&ˆ

 

 

I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

xˆ

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

наблюдатель

xˆ

Рис.15.1.Система с модальным регулятором и наблюдателем

141

Модальному регулятору соответствует закон управления:

u = V - Kxˆ .

(15.3)

Допустим, что V=0 (для простоты). Подставим (15.3) в (15.1) и (15.2) и объединим выражения (15.1) и (15.2) тогда получим:

x&= Ax - BKxˆ

;

(15.4)

ˆ

= (A - LC - BK)xˆ

x&

+ LCx

 

здесь учтено, что Ly=LCx.

x

A

142

 

 

 

 

-BK x

 

&

=

 

 

×

 

 

(15.5)

ˆ

LC

A - LC - BK

xˆ

&

 

 

 

 

x

 

 

 

 

Это уравнение системы с наблюдателем и модальным регулятором.

Перейдем от переменных (x,x)ˆ к переменным (x,е), где e = x - xˆ .

Первое уравнение получим путем добавления и вычитания слагаемого BKx к первому уравнению системы (15.4), второе уравнение получим, вычитая из

(15.1) уравнение (15.2).

x&= Ax - BKxˆ + BKx - BKx = (A - BK)x + BK(x - x)ˆ

 

ˆ

= Ax + Bu - (A - LC)xˆ - Bu - Ly = Ax - (A - LC)xˆ - LCx = (A - LC)(x - x)ˆ .

x&

- x&

Запишем в векторно-матричной форме:

x

A - BK

BK

x

 

&

 

=

 

×

(15.6)

 

0

e

 

 

A - LC

e

 

&

 

 

 

 

 

 

Матрица последней системы треугольная, поэтому характеристический полином замкнутой системы d(p) удовлетворяет уравнению:

d(p)= b(p) c(p),

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