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

Лекции / Т4_ДГМ

.pdf
Скачиваний:
9
Добавлен:
16.05.2021
Размер:
400.66 Кб
Скачать

ˆ

 

ˆ

LB1)u(t) + (A21 LA11)Y(t) + LY(t)

(4-10)

W(t) = (A22

LA12 )W(t) + (B2

После не сложных преобразований структурную схему наблюдателя можно представить как

y

 

A -

11

L

 

21

LA

 

u

B2-LB1

 

 

1

W

 

 

 

 

 

 

 

p

 

 

 

 

 

A22 -LA12

 

На основании преобразованной структурной схемы динамика наблюдающего устройства этого типа будет определяться характеристическим уравнением вида

 

det( pI - A22 + LA12 ) = 0

 

Для

наблюдателя полного порядка, как ранее

было показано,

характеристическое уравнение имеет вид

 

 

det( pI - A + LC) = 0

 

4.2.3

Пример проектирования редуцированного

наблюдающего

устройства

Требуется спроектировать редуцированное наблюдающее устройство для объекта управления, заданного следующей структурной схемой. Известно, что доступной измерению является скорость исполнительного органа ω2 = x2 . Требуется спроектировать наблюдающее устройство для оценки переменных ω1 = x3 и Mу = x2 , т.е. определить структуру и

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

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

 

-1

 

 

 

 

x2

 

x1

M(p)

1

x3

 

1

1

 

Tм1p

 

 

Tcp

Tм2p

 

 

объект

 

 

-1

 

 

 

 

 

 

 

 

 

 

l11

 

y1

 

 

W1 = x2

 

 

1

 

 

1

 

 

 

 

-T

 

 

p

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

-

111

 

 

 

 

Tм2

 

 

 

 

 

 

1

 

 

 

 

 

 

 

Tc

 

 

 

 

 

 

-( 1

 

+

 

121 )

 

 

 

 

Tм1

 

Tм2

W2 = x3

 

1

 

 

1

 

 

 

Tм1

 

 

p

 

 

 

 

 

 

 

l21

 

 

Принимая во внимание уравнение редуцированного наблюдателя

ˆ

ˆ

LB1)u(t) + (A21 LA11)Y(t) + LY(t),

W(t) = (A22

LA12 )W(t) + (B2

определим для нашего объекта все матрицы, входящее в это уравнение. С учётом введённых обозначений переменных состояния векторно–матричное описание рассматриваемого объекта получит вид:

45

 

 

 

 

0

 

 

 

x1

 

 

 

 

 

 

 

 

 

 

 

 

1

x2

 

=

 

T

 

 

 

 

 

c

x3

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Tм2

 

x1

 

0

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

+

 

 

 

 

 

0

 

 

 

 

0

M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

1

 

 

 

 

 

1

 

 

Tc

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

3

 

T

 

 

 

 

 

 

 

 

 

 

 

м1

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

м1

 

 

 

 

 

 

 

 

 

 

 

Из этого описание несложно определить, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

1

 

A

 

= 0; A

 

 

=

 

1

0

; L =

l11

 

; B = 0 ; B

 

=

1

 

; A =

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Tc

 

11

 

 

 

12

 

 

 

 

 

1

2

 

 

 

21

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Tм2

 

 

l21

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Tм1

 

 

 

 

 

 

 

0

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

22

1

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

м1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычислим коэффициенты наблюдающего устройства, исходя из желания обеспечить его характеристическому уравнению наперёд заданное распределение корней, например стандартное распределение Баттерворта. С этой целью определим характеристическое уравнение наблюдателя и составим систему линейных уравнений, приравнивая его коэффициенты при одинаковых степенях оператора «р» к соответствующим коэффициентам стандартного распределения Баттерворта второго порядка:

 

 

 

 

 

 

) = p2

 

 

l

1

 

1

 

l

21

 

= p2

 

 

p + ω2

det( pI - A

22

− LA

12

+

 

11

p +

 

 

 

+

 

 

+ 1,4ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Tм2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Tc Tм1

 

Tм2

 

 

 

 

l11

=1,4ω

;

 

1

 

+

 

l21

 

= ω02н .

 

 

 

 

 

 

 

 

 

 

 

TcTм1

 

 

 

 

 

 

 

 

 

 

 

 

 

Tм2

 

 

 

TcTм2

 

 

 

 

 

 

 

 

 

 

 

Следовательно, коэффициенты наблюдающего устройства необходимо

 

 

 

 

 

 

 

ω2

 

1

 

 

выбрать равными: l

=1,4ω

T

; l

21

=

 

 

T T .

 

11

 

м2

 

 

 

 

 

c м2

 

 

 

 

 

 

 

 

 

TcTм1

 

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

46

 

 

 

 

 

 

 

l

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11

 

 

 

 

 

 

 

 

 

0

 

 

 

 

1

 

 

 

 

 

&

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

=

 

 

Tм2

 

 

Tc

1

 

+

1

 

 

 

+

 

 

y +

l11

 

y&

 

 

 

 

∆М

 

 

&

 

 

 

 

1

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

Tc

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

21

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

+

 

 

 

0

 

 

T

 

 

 

0

 

l21

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

T

 

 

 

 

 

м1

 

 

 

 

 

 

 

 

 

 

 

 

 

м1

 

 

м2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4.3 Проектирование линейных оптимальных систем для

управления объектами со случайными возмущениями типа «белого

шума».

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

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

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

 

 

 

 

 

 

 

r

;

 

 

Х

q

 

 

 

Y ZХ

qK

;

s,

 

(4-11)

 

 

Х .

 

 

 

 

где - выходная переменная; qr -векторный случайный процесс

типа «белого шума», возмущающий состояние X(t);

qK -белый шум

наблюдений, аддитивно наложенный на измеряемые сигналы.

Совместный

случайный

процесс

qr

и

qK являются

некоррелированным и характеризуется блочной матрицей интенсивности

 

 

 

t

0

 

 

 

 

 

t o 0

t

 

 

 

 

 

 

 

 

 

47

= Vvwxy z - wv

Доступным измерению является «r» переменных состояния из общего числа n. Задача построения линейной оптимальнойсистемы, включающей в себя наблюдатель и обратные связи по восстановленным переменным, состоит в определении такого управления u Y , при котором критерий:

(4-12)

достигает минимума при удержании выходной переменной z в нулевом заданном состоянии zз=0

Положительно определенные весовые коэффициентыты и определяют вклад в критерий качества средних значений квадрата отклонения управляемой переменной от предписанного (нулевого) значения и квадрата управления который при этом затрачивается. Сформулированная задача решается на основе принципа разделения. Применение данного принципа позволяет доказать, что матрица коэффициентов «*» оптимального регулятора, замкнутого по полному вектору состояния остается неизменной и в том случае, когда замыкание системы производится по оценкам вектора состояния \ , восстановленным наблюдающим устройством, параметры которого определены из условия минимума среднего квадрата ошибки восстановления. Таким образом, оптимальное управление, обеспечивающее

минимум критерия (4-12) формируется в виде:

* , где матрица коэффициентов обратных связей «*» как и при решении задачи

детерминированного оптимального линейного уравнения определяется.

выражением:

* + ,( ,

 

в котором значения матрицы «(» являются единственным неотрицательным

симметричным решением алгебраического уравнения Риккати:

 

0 , ( + ,( ,( (

48

При необходимости стабилизировать выходную переменную « » на уровнях, отличных от нулевого, требуемое управление может быть определено как:

\

 

0

 

* {

,

 

з+

з

где {з – передаточная функция замкнутой системы от входа до управляемой переменной, - задающее воздействие. Оценка полного вектора состояния объектаз \ выполняется установившемся оптимальным

наблюдателем, который реализуется в соответствии с выражением (4-8):

? ?| ] Y Z?|

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

uу

ωη

 

 

ωµ

-1

 

 

Wз(0)

1

x

 

y

B

C

p

 

 

 

A

 

 

 

-1

L

 

 

 

 

 

 

 

B

1

x

-C

y

p

 

 

 

A

 

 

 

K

 

 

 

 

1.

 

 

 

 

Рисунок 4.5 Установившийся наблюдатель является оптимальным в том смысле, что

при определенном выборе ( 4 ) матрицы коэффициентов « ] »

минимизируется среднее значение квадрата ошибки восстановления, } ? ?| , то есть:

49

Vvwx }, @} - wv,- ∞

где @ – наперед заданная положительная весовая матрица. Дифференциальное уравнение ошибки восстановления может быть

получено путем вычитания из первого уравнения (4-11) уравнения (4-8),

принимая при этом во внимание, что

Y Z

qK т.е

 

 

 

qr

 

 

\

\

\

 

 

K

 

]Z ]q

]Z

В результате, получим

 

 

qr

 

 

} ]Z }

b 1 ]

 

 

~qK

 

Решение задачи построения установившегося оптимального наблюдателя достигается выбором матрицы коэффициентов наблюдателя как

] /Z,t+ (4-13)

Если обозначить среднее значение вектора ошибок восстановления}€ Мy} z, то в выражении (4-13 ) –матрица

/ xy} }€ } }€ ,z

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

0 t /Z,t+ Z/ / /,

(4-14)

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

50

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

4.4 Проектирование линейных оптимальных систем для управления объектами со случайными возмущениями типа «цветного

шума».

Рассмотрим возможность построения линейной оптимальной системы для объектов управления с цветным шумом внешних воздействий и критерием качества (4-12). Пусть стационарный объект управления порядка « » задается уравнением:

 

 

&

 

u

 

 

 

 

 

r r

 

 

&

 

Y Z& & uK

 

(4-15)

 

 

& &

 

 

 

 

где ur и uK - случайные возмущающие воздействия с определенными статистическими характеристиками, полученными при фильтрации некорелированных белых шумов mr и mK через соответствующие фильтры первого порядка:

ur rur

rqr

 

 

KqK

(4-16)

uK KuK

 

 

 

(4-17)

Считается, что доступными измерению являются « 4 » переменных состояния объекта управления из общего числа « ».

Векторно-матричная структурная схема, соответствующая такому математическому описанию объекта может быть. представлена как показано на рисунке 4.6.

51

ωη

 

 

ω

 

 

 

 

 

µ

 

 

Gη

 

G

 

 

 

fη

 

µ

 

fµ

1

 

 

1

p

 

 

 

p

 

Aη

Bη

 

Aµ

y

u B0

 

1

x0

C

 

p

 

 

 

 

A0

 

 

 

Рисунок 4.6

При таком описании объекта задачу построения оптимальной системы управления стремятся привести к ранее рассмотренной, когда шум возбуждающий состояние объекта и шум измерений являются случайными процессами типа белого шума. С этой целью часть описания, относящуюся к моделированию цветного шума ur a , относят к самому объекту и в дальнейшем рассматривают математическое описание расширенного объекта, которое на основании уравнений (4-15), (4-16) может быть представлено как

 

 

 

 

 

 

 

&

 

&

0

 

 

 

 

 

 

 

 

 

 

r

 

S

~

• o

 

&

 

r o

r

 

 

ƒ

 

 

 

 

 

 

or q

(4-18)

 

&

 

 

r

u

 

 

 

 

 

 

 

 

r

 

0

 

 

0

 

 

 

 

 

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y Z&

 

0

&

uK

 

 

 

Q

 

 

 

our

 

 

 

Введём соответствующие обозначения блочных матриц в системе выражений (4-18), которые в общем случае, когда все переменные вектора состояния & возмущаются, будут иметь следующие размерности:

 

o 0 r 2 2 ;

&

 

 

 

A=

0

ƒ 2 1 ; our

& r

 

 

&

 

 

2 1 ;

52

Z Z& 0 4 2 ;

 

 

0

 

 

 

 

 

• or 2 .

 

С учётом принятых обозначений, уравнения расширенного

объекта можно переписать в виде:

 

•q

;

 

 

 

 

 

 

„ Z?

uK .

r

 

(4-19)

 

 

 

 

 

 

,

где uK - по-прежнему цветной шум.

Соответствующая (4-17) структурная схема показана на рисунке 4.7

 

 

ωµ

 

 

ωη

Gµ

1 fµ

 

 

 

 

 

p

 

G

 

Aµ

u B

1

x

y

p

C

A

Рисунок 4.7

Таким образом, включение описания формирующего фильтра в состав описания объекта управления изменяет первоначальную формулировку задачи. Теперь расширенный объект возбуждается воздействием типа «белого шума», поэтому коэффициенты жестких обратных связей при формировании оптимального закона управления могут быть как и прежде рассчитаны путем решения соответствующего алгебраического уравнения Рикатти (4-5).

Следует отметить, что расширенный вектор включают в себя не только переменные состояния объекта & , но и возмущения, действующие

53

Соседние файлы в папке Лекции