Файл: Основы автоматического управления..pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 09.04.2024

Просмотров: 343

Скачиваний: 15

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

§ 14.8. СИСТЕМ Ы , О П И СЫ ВА ЕМ Ы Е

Д И Ф Ф Е Р . У Р А В Н Е Н И Я М И

591

весовой функции оптимальной системы:

 

 

Фі (t)фі (т)

( — оо < Т < 1 ) ,

(14.8.1)

е(*. т) = к(fcii+ eii)

где

(14.8.2)

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

Л

ФІ (*)

(14.8.3)

*п + сіі

 

(для удобства последующих выкладок мы пишем здесь просто г) вместо т)т1п). Логарифмируем формулу (14.8.1), а затем дифференцируем по t с учетом (14.8.2). В результате получим

g(t, т)= фі(0

фi(t)

g ( t , т) ф± (0

к { Ь п + С ц ) ’

где точкой, как и везде в дальнейшем, обозначено дифференцирование по t.

Приняв во внимание (14.8.3) и

положив

для краткости

 

.

Фі (О

 

(14.8.4)

 

Фі (0

 

 

приведем полученное уравнение к виду

 

 

g(t, т ) = ( а —

Т)'

(14.8.5)

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

мальной системы, достаточно заменить начальное

условие

g (т, т) =

= Ф'] (т)/к (6И -f сц) соответствующим слагаемым в

(14.8.5),

содержащим

б-функцию б (і — т) (см. §§ 4.4 и 4.5). Так как при любом т весовая функция

g (t, т) при переходе t через т изменяется скачком от 0 до ф|

)!к (6ц -j- cu ),

то б (t — т)

должна

войти

в правую часть (14.8.5) с

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

Фі (t)/k (Ьи +

di). Заменив в полученном таким путем уравнении б (f — т)

входным сигналом Z = Z (t),

а весовую функцию g (t, т) выходным сигналом

оптимальной системы

W* =

W* (t),

получим дифференциальное уравнение

оптимальной

системы:

 

 

 

 

 

 

1Р’* = ( а —

W* + ^ - Z .

(14.8.6)

Анализируем структуру оптимальной системы. Так как требуемый выходной сигнал w — W (t) в данном случае равен t/іфі (t), то из (14.8.4) вытекает следующее уравнение для требуемого выходного сигнала:

W = aW.

(14.8.7)

Сравнивая это уравнение с (14.8.6), видим, что оптимальная система полу­ чается из системы формирования требуемого выходного сигнала путем


592

ГЛ . 14. О П Р Е Д Е Л Е Н И Е

О П ТИ М А Л ЬН Ы Х Л И Н Е Й Н Ы Х

СИСТЕМ

замыкания £ ее ~отрицательной)

жесткой обратной связью

и включения на

входе усилителя с коэффициентом усиления г]/к. На рис. 14.8.1 показана оптимальная система. Пунктирной линией обведена система формирования требуемого выходного сигнала.

Величина rj, входящая в коэффициент усиления входного сигнала, может быть вычислена заранее и введена' в систему с помощью соответствующего программного устройства. Однако часто ее предпочитают вычислять непо­ средственно в системе в процессе ее работы (или в модели системы при ее моде­ лировании). При этом вычисление т) можно осуществить путем интегрирова­ ния дифференциальпого^уравпения, которому она удовлетворяет. Чтобы

получить это уравнение,^прологарифмируем формулу (14.8.3), а затем про­ дифференцируем с учетом (14.8.2) и (14.8.4). В результате получим

Д

. 2Ф1 (t)______ Фі (t)

^

Л

Д

Фі (і)

^(Ь ц + сц ) “

 

к

Отсюда вытекает дифференциальное уравнение для rj:

Ч = 2 а п - - у - .

(14.8.8)

Интегрируя совместно уравнения (14.8.6) и (14.8.8) (уже нелинейные), можно получить в каждый момент времени t одновременно выходной сигнал опти­ мальной системы W* и соответствующий средний квадрат ошибки т).

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

Рассмотрим общую задачу воспроизведения регулярного тре­ буемого выходного сигнала

W ( t ) = |] £ / ггМ<)

(14.8.9)

Г=1

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

IV-1

5 ( 0 = S ßkW"fc>(*).

(14.8.10)

fc= 0

 

а помеха N (t) = X (t) представляет собой белый шум интенсивно­ сти G (t). Подставив выражение (14.8.9) в (14.8.10) и изменив поря­ док суммирования, представим входной сигнал в привычном нам


§ 14.8. СИ СТЕМ Ы , О П И СЫ ВА ЕМ Ы Е Д И Ф Ф Е Р . У РА В Н Е Н И Я М И

593

виде (13.5.14):

 

 

 

N

 

 

(14.8.11)

 

 

z(t)= 2 Ur<?r(t)+x(t),

где

 

 

r = i

 

 

 

 

N- 1

 

 

 

 

 

 

 

 

 

(14.8.12)

 

Фг ( t ) = 2 Р » Г ( 0

( r =

l ,

 

 

fc=0

 

 

 

 

В данном случае

согласно

(13.6.4)

при

g(0)(i, т) = 0

и (14.1.4)

весовая функция

оптимальной системы

с бесконечной

памятью

(Т —t t0)

определяется

формулой

 

 

 

 

N

 

 

 

 

I z V’

=

г—1

 

( - о с < т <<),

(14.8.13)

 

 

 

 

 

 

где Xi (t), . . ., XN (t) — функции, определяемые системой урав­ нений (13.6.7), которая в данном случае, когда Ъро — 0 (р =

1, . . ., N), принимает вид

N

(Р = 1, ...,Л Т). (14.8.14)

2 Фрд cpq) ^g (0 — Фр (0

9=1

 

Формула (13.6.6), определяющая величины bpq, в данном случае принимает вид

Ъря= I ФР (1\1) (Т~ dx

(Р, 9 = 1 , ...,А Т).

(14.8.15)

*0

 

 

Легко убедиться непосредственной подстановкой, что требуе­ мый выходной сигнал W (t) удовлетворяет линейному дифферен­ циальному уравнению

W

Фі (t)

. Флг (t)

 

 

w f

ф; (*)

• • Фіѵ (t)

= 0.

(14.8.16)

 

 

 

w lN) Ф(Л (0

. .. Фw°(0

 

 

Совершенно так же, как в примере 14.8.1, дифференцируя формулу (14.8.13) N раз с учетом (14.8.14) и (14.8.15), можно убе­ диться в том, что весовая функция оптимальной системы в данном случае определяется линейным дифференциальным уравнением N-то порядка, достаточно просто связанным с уравнением (14.8.16). Однако необходимые для этого выкладки очень громозд­ ки, утомительны и мешают пониманию существа дела. Чтобы упростить задачу, целесообразно, во-первых, заменить ее задачей одновременного воспроизведения сигнала W и его производных

38 Под ред. В. С. Пугачева


594

Г Л . 14. О П Р Е Д Е Л Е Н И Е О П ТИ М А Л ЬН Ы Х Л И Н Е Й Н Ы Х СИСТЕМ

W , . . ., во-вторых, заменить уравнение (14.8.16) равно­ ценной системой уравнений первого порядка и в-третьих, вместо формул (14.8.13), (14.8.15) и уравнений (14.8.14) исходить непо­ средственно из уравнения (13.5.37) и формулы (13.6.25) и произ­ вести все выкладки в матричной форме. Заметим еще, что выклад­ ки почти не изменятся, если в правые части системы уравнений формирования требуемых выходных сигналов добавить произ­ вольные белые шумы. Это равносильно добавлению нерегулярных частей полезных сигналов, формируемых из белых шумов той же системой дифференциальных уравнений. Если нерегулярные части требуемых выходных сигналов формируются другой систе­ мой дифференциальных уравнений, то ее следует просто добавить

куравнениям, определяющим регулярные части этих сигналов.

Всоответствии со сказанньм рассмотрим задачу определения оптимальной многомерной линейной системы в случае,-когда век­ торный требуемый выходной сигнал W (t) определяется матричным дифференциальным уравнением

W = A W + V { t) ,

(14.8.17)

а векторный входной сигнал определяется формулой

Z (t) = B W (t) + Vi (t),

(14.8.18)

где А и В — матрицы коэффициентов (в общем случае перемен­ ных), V (t) и Уі (t) — некоррелированные векторные белые шумы, в общем случае нестационарные. На основании сказанного в § 13.5 можно без потери общности считать математические ожи­ дания белых шумов V (t) и Ѵі (t) равными нулю. Матрицу интен­ сивностей и взаимных интенсивностей составляющих белого шума V (t) обозначим G. Матрицу интенсивностей и взаимных интенсив­ ностей составляющих помехи Ѵ\ (t) обозначим Gt. Тогда матрицы корреляционных и взаимных корреляционных функций белых шумов V (t) и Ѵі (t) выразятся соответственно формулами

K v (t,

т) = М IV (t) V (т)т1 = G b ( t - т),

(14.8.19)

Кп (t,

г) = М [Ѵі (0 Ѵі (т)т] = Gfi (t — т).

(14.8.20)

Для краткости

мы не показываем явно зависимость матриц

А, В, G и Gi от времени t. Однако в случае, когда их аргумент обо­ значен другой буквой, мы будем указывать его соответствующим нижним индексом.

Для вывода дифференциальных уравнений оптимальной систе­ мы воспользуемся уравнением (13.5.37), определяющим матрицу весовых функций оптимальной системы, и формулой (13.6.25), определяющей матрицу вторых начальных моментов вектора ошиб-


§ 14.8. СИСТЕМ Ы , О П И С Ы ВА ЕМ Ы Е Д И Ф Ф Е Р . У РА В Н ЕН И Я М И

595

ки оптимальной системы *). Положив в (13.5.37) и (13.6.25) интер­ вал наблюдения Т равным полному времени работы системы, Т = = t tQ, перепишем эти равенства в виде

і

 

 

 

j g(t, т)Г2(т, a)dx=Twz (t, ст)

(?0< o < i) ,

(14.8.21)

<0

t

 

 

 

 

 

 

Н = ГШ((, t)— j g(t, x)Twz(t, t)tdx.

(14.8.22)

 

*0

 

 

Кроме этих соотношений нам понадобятся еще некоторые вспомо­

гательные соотношения.

Умнолшв уравнение (14.8.17) справа

на Z (ст)т и имея

в виду,

что требуемый выходной сигнал W (ст),

а следовательно,

и входной сигнал Z (ст) в момент ст < t не корре­

лирован со значением белого шума V (t) в последующий момент времени t, получим

М [W Z (а)»] = A M [W (t) Z (а)т1,

а <

t,

или

 

 

 

f wz(t, o) = ATwz(t, а),

о < t .

 

(14.8.23)

Из (14.8.18) ввиду некоррелированности Z(a) с Ft (t) при ст •< t следует

Гг (t, о) = M [ Z (t) Z (<т)т) = ВМ IW (г) Z (а)Н = BTwz(t, а).

Подставив сюда выражение Ги,2 (t, ст) из (14.8.21), получим

t

 

Tz (t, о) = В j g(t, т)Г2(т, ст)dx

(f0<CT<?). (14.8.24)

Аналогично из той же формулы (14.8.18) следует

Tzw(т, ст) = M I Z (т) W (ст)Н = В ТМ IW (т) W (ст)Н,

или

Tzw (т, ст) = B XTW(т, ст),

(14.8.25)

где, в соответствии с нашим условием, В х — значение матрицы В при t — X.

Теперь мы можем перейти к выводу дифференциального урав­ нения оптимальной системы. Дифференцируя уравнение (14.8.21) по t и имея в виду, что от t зависит и подынтегральная функция

*) Читатель, не владеющий основами матричной алгебры, может про­ пустить этот вывод и прочитать сразу первый абзац после уравнения (14.8.31),

вкотором формулируется окончательный результат, и после этого перейти

кчтению первого абзаца после формулы (14.8.37).

38*