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

книги из ГПНТБ / Цейтлин Я.М. Проектирование оптимальных линейных систем

.pdf
Скачиваний:
2
Добавлен:
24.10.2023
Размер:
9.35 Mб
Скачать

где

т (t) — неслучайный

сигнал, преобразование Лапласа

кото­

рого

есть

L [т (t)] =

М

(s);

п (t) — случайный

сигнал

со

спек­

тральной

плотностью Фпп (s);

 

 

 

 

 

 

 

 

z(f) =

H (р) т (О,

 

 

 

где Я (р) — заданный

оператор.

 

 

 

 

Ошибку є можно

представить в

виде

 

 

 

 

є = Ит — F (т + п) = (Я — F) т — Fn = гх +

S 2 ,

 

где 8j =

( Я — F) m и е2

= Fn — неслучайная

и случайная со­

ставляющие ошибки.

 

 

 

 

 

 

 

Будем

синтезировать

систему,

обеспечивающую

минимум

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

Иными словами, критериями оптимальности системы являются

 

I со

=

1 ®nnFFds = mm;

 

— ( С О

/ e i = 2 s r I

(H-F){H-F)MMds^c.

— І CO

 

Для решения такой задачи необходимо найти минимум нового

функционала

 

 

 

 

 

 

D

= D&2 +

ХІЄі,

где Я.— неопределенный

множитель

Лагранжа, также подлежа­

щий определению.

 

 

 

Функционал

D имеет вид

 

 

t'co

 

 

П=

I

[®n»FF +

HH-F)(H-F)MM\ds,

 

—і

со

 

 

т. е. совершенно аналогичен тем функционалам, которые были рас­ смотрены ранее.

Дифференцируя подынтегральное выражение по F, получим уравнение Винера — Хопфа в комплексной области

O n „ F — К (Я - F) ММ = £,

или

F (Ф„„+ ШМ) — %НММ = I .

Решением этого уравнения будет

F =

"

_

ммн

 

(Флп + ШИ)"

 

(<ьт

+

шм)+

+

 

 

 

 

Естественно, что найденная передаточная функция F зависит от множителя Лагранжа X, который находится в результате подста­ новки функции F в выражение

1'со

 

7 « « = 2 Й Г \

(H-F)(H-F)MMds^c

I со

ирешения полученного неравенства.

Аналогичным образом решаются и другие задачи с неслучай­ ными сигналами [1, 13].

22. УРАВНЕНИЕ ВИНЕРА—ХОПФА ВО ВРЕМЕННОЙ ОБЛАСТИ

Приведем вывод уравнения Винера—Хопфа и его решение во временной области. Вернемся к задаче выбора структуры системы, схематично изображенной на рис. 14.

Это система управления со случайным сигналом на входе х (t). Будем считать, что х (t) — случайная стационарная центрирован­ ная функция. Система управления, которая предполагается линей­ ной и стационарной, описывается функцией веса Р (t) или ее пре­ образованием Фурье—передаточной функцией F (s).

Через у (t) обозначен сигнал на выходе этой системы и через z (t) — идеальная функция на выходе. Разность этих двух сигна­

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

собой

ошибку

системы

 

 

 

Є

(0

= z ( 0 - # ( 0 .

(Ш.48)

Заметим,

что

систему

управления

всегда можно

привести

к виду, изображенному на рис. 14.

 

 

Нашей задачей

является такой выбор функции веса Р (f)

или передаточной

функции

F (s),

при

которых обеспечивается

минимизация

дисперсии

ошибки

е (/).

 

 

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

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

Искомая функция веса Р (t) находится как решение некоторого интегрального уравнения, которое получается следующим обра­ зом.

В выражении ошибки на выходе (Ш.48)

со

У(0=

J PVJxV-tj'db.

(Ш.49)

Возведя

в

квадрат

обе

части

(III.48),

найдем

 

 

 

 

 

 

 

Е 2

{ t

) =

гг

( 0

_

2

z

(t) у

(t)

+ у"

(t)

 

 

 

или с

учетом

(И 1.49)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

со

 

 

 

 

 

 

 

 

 

 

 

е2

(0 =

г2 (0 -

2z

 

J Р (tj) x(t-

tx)

dtx

+

 

 

 

 

 

 

 

 

 

 

 

 

 

—со

 

 

 

 

 

 

 

 

 

 

 

 

со

со

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

\

\

 

 

 

 

 

 

 

PitjPitJxit-tJxit-tJd^dtz.

(Ш.50)

Применив к обеим частям уравнения (III.50) операцию опреде­

ления

математического ожидания

и имея в виду, что

 

 

 

 

 

 

 

 

М

2

{t) \

=

КЕг(0)

=

DB,

 

 

 

 

 

так как

є (t)

— стационарный

центрированный

случайный

про­

цесс,

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

оо

 

 

 

 

 

 

 

 

 

 

 

 

DE = M

г2 (0 — 2 J x{t

— t1)z(t)P(t1)dt1

 

+

 

 

 

 

 

 

 

 

 

 

 

— со

 

 

 

 

 

 

 

 

 

 

 

 

 

 

со

со

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

J

\

х (t — tj) х (t -

t2)

P

(tj) P (t2)

dtx

dt

 

 

 

 

 

 

— CO — C O

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

CO

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

K * ( 0 ) - 2

J

 

 

 

KMPitJdt^

 

 

 

 

 

 

 

+

\

J Kxx(ti-ti)P(t1)P(t,)dt1dt2,

 

 

 

 

(Ш.51)

 

 

 

 

— CO — 0 0

 

 

 

 

 

 

 

 

 

 

 

 

 

поскольку

по

определению

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

[z2

(01

=

 

KB(0);

 

 

 

 

 

 

 

M

[x

(t-tj

z(t)]

 

=

Кхг

(

t -

t +

tx)

=

Кхг

(tj;

 

 

[M

[x (t—tj

x

(t-t2)}

 

=

Kxx

 

(t—tj—t+U)

=

Kxx

(tr-h)

 

=

 

 

 

 

 

 

 

 

~

Kxx (^1

^2)-

 

 

 

 

 

 

Определение

функции

веса

 

P

(t)

сводится,

таким

образом,

к минимизации

функционала

De.

Воспользовавшись известным

' методом

вариационного

исчисления,

представим функцию

веса

в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где Pm(t)

 

 

 

 

 

Р (0 = Pm(t)

+

ag

(t),

 

 

 

(ПІ.52)

— искомое

решение,

т. е. функция веса

оптимальной

системы, обеспечивающая

минимизацию функционала D e ;

g

(t) —

некоторая произвольная функция веса; а — произвольная кон­ станта.

Условие, которому должна удовлетворять искомая функция веса Pm(t), получится в результате приравнивания нулю производ­

ной -j^Dt

П Р И

а

= 0.

 

 

 

 

 

 

 

 

 

 

Подставив

(III.52)

в (III.51),

получим

 

 

 

 

 

 

 

 

 

 

 

 

со

 

 

 

 

 

 

 

 

 

 

А* =

Кв

(0) -

2

J Рт

(tj

Кхг (k)

dtx

-

 

 

- 2 а

 

J g

 

Кхг

(k)

dtx+

J

j

 

Pm(k)Pm(k) X

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

X

xx

(k

-

k) dk

 

+

«

J

J P

m

(k) g

(k)

X

 

 

K

 

 

 

 

 

 

 

 

 

 

 

 

— CO

— C O

 

 

 

 

 

X Kxx

(k -

k) dk

dk

+ « 2

j \

g

(k) g (k)

Kxx

(k

- k) dk dk.

(Ш.53)

Заметим, что пятое и шестое слагаемые в правой части (III.53) равны между собой. Действительно, переменив порядок интегри­ рования и индексы переменных интегрирования в пятом слагаемом

со со

J

J

Pm(k)g(k)Kxx(k-k)dkdk,

— со —со

 

получим

 

 

со

со

 

|

J

Pm(k)g(k)Kxx(k-k)dkdt2,

С О — ' С О

т.е. шестое слагаемое.

Всилу четности корреляционной функции

Кхх (к к) ~ Кхх {t% ti),

поэтому

 

со

со

 

 

 

 

 

 

 

 

 

 

 

 

— C O

с о

 

 

 

 

 

 

 

 

 

 

 

 

=

J

 

J Рт (t2)8

(к) Кхх

(tx

— t,) dtx

dt2.

 

 

 

— с о

— с о

 

 

 

 

 

 

 

 

 

 

Таким образом,

 

выражение

(III.53)

принимает

вид

 

 

 

 

 

 

 

со

 

 

 

 

 

 

 

 

 

DB

=

Кгг(0)

- 2

J

Рт

(tx)

Кхг

(tx)dtx+

 

 

 

 

 

 

 

 

 

— C O

 

 

 

 

 

 

 

 

 

со

 

со

 

 

 

 

 

 

 

 

 

 

 

 

— с о

— с о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

со

 

 

 

 

 

 

 

 

 

 

 

 

- 2 а

\ g(h)Kxz(ti)dtx

 

+

 

 

 

 

 

 

 

 

 

—со

 

 

 

 

 

 

 

 

 

 

 

со

 

с о

 

 

 

 

 

 

 

 

 

 

+

2 а

J

 

\ PM{U)g{t1)Kxx{t1-tjdt1dtt

 

 

 

+

 

 

 

— с о

— с о

 

 

 

 

 

 

 

 

 

 

 

 

 

со

со

 

 

 

 

 

 

 

 

 

 

 

+ а 2

— 0 0

—со

iti) § (*,) Кхх

(tx

- t 2

) dtx

dtt.

 

 

 

 

J

\ 8

 

 

 

Продифференцировав

обе части этого выражения

по а

и при­

равняв полученную

производную нулю

(при а

=

0), получим

 

со

со

 

 

 

 

 

 

 

 

 

 

 

|

 

J Рт

 

 

Кхх

(tx - t2) dtx dt2

-

 

 

 

— CO — с о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

со

 

 

 

*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

0,

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

со

со

 

 

 

 

 

 

 

 

dtx

=

0.

(III.53a)

jg(h)

\

 

 

 

 

Pm{tz)Kxx{ti-t%)dt2-Kxz{ti)

Функция g (ti) является произвольной функцией веса и, сле­ довательно, должна соответствовать некоторой физически осуще-

ствимой системе. Иными словами g (t)

= 0 при t

< 0 , а в осталь­

ном совершенно произвольна.

При

t ^> О для

 

удовлетворения

(III.53а) необходимо равенство

нулю

выражения,

заключенного

в скобки.

 

 

Таким образом, получается

условие

 

со

 

 

т{к)КхЛк-к)<Игхг{Ч)

= Ъ, t^O.

(ІП.54)

Этому условию и должна удовлетворять функция веса оптимальной (в смысле минимума дисперсии ошибки) системы.

Интегральное уравнение (III.54) носит название уравнения Винера—Хопфа (во временной области).

Заметим, что, вообще говоря, функция веса Рт (t), удовлетво­ ряющая (III.54), обеспечивает стационарное значение дисперсии ошибки, но не обязательно минимум. Однако из физических представлений очевидно, что решение, если оно существует, должно давать минимум, так как максимум дисперсии ошибки равен беско­ нечности.

То,

что функция

веса

Рт (t) обеспечивает именно минимум,

можно

показать

и

математически.

 

 

 

 

 

Продифференцировав

дважды

выражение

дисперсии ошибки

є (t),

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

со

со

 

 

 

 

 

 

 

 

 

^ D

s

=

2

j

J

KxAti-tJgiUgiWbdt,.

 

(HI.55)

 

 

 

 

 

 

— с о

— - c o

 

 

 

 

 

 

 

Правая часть (III.55) представляет собой

дисперсию сигнала ух (t)

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

на вход которой посту­

пает

сигнал х

(t):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С О

0 0

 

 

 

 

 

 

 

 

 

Dyi

 

=

2

j

j Кхх (U -t2)

g (ti) g (t2)

dti dt2.

 

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

Поэтому

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d*Dda? s

 

 

 

 

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

условие

=

0, определяющее

экстремум

дисперсии, соответствует

именно

 

минимуму.

 

 

Перейдем

к

решению

интегрального

уравнения

Винера—

Хопфа,

которое

согласно

(III.54) имеет

вид

 

 

 

 

со

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J Рт

(tx) Кхх (т -

h) dtx

К^х)

 

= О,

т ^ 0.

(Ш.55а)

>-« \

Если бы левая часть (III.55а) была равна нулю для всех значе­ ний т, получилось бы интегральное уравнение Фредгольма первого рода, которое можно решить с помощью преобразования Фурье при условии, что все функции Рт (t), Кхх (t) и Кхг (t) имеют пре­ образования Фурье.

Однако левая часть (III.55а) равна нулю только при положи­ тельном значении времени и, следовательно, ее преобразование дает функцию аргумента s, не равную нулю, в то время как в слу­ чае уравнения Фредгольма первого рода эта функция равнялась бы нулю. Поэтому уравнение (III.55а) желательно преобразовать в уравнение первого рода путем выполнения некоторых операций

во временной

области.

 

 

 

 

 

Перейдем

к

этим операциям.

 

 

 

 

 

Представим

четную функцию

Кхх

(т)

в

виде

свертки двух

функций КЇк

(т)

и

К7х (т):

 

 

 

 

 

 

 

j

К7Х (t2) Кїх (т -

Ъ) dt2

=

Кхх

(т).

(Ш.56)

Функции Kfx (т) и Кхх (т) определим следующим образом:

+

 

(т) =

0,

т < 0 ;

Кхх

 

 

(111.57)

Кхх (т) =

0,

т > 0 .

По определению Кхх

(т) может

иметь ненулевое значение как

при положительных, так и при отрицательных значениях време­

ни

T v

что характерно

для четной

функции.

 

 

 

 

Функцию Кхг (т)> вообще говоря,

не

являющуюся

четной,

также

представим

в виде свертки

двух

функций:

 

 

 

 

 

J Kxx{h)y{x-t2)dt2

=

Kxz{t),

 

 

(ІИ.58)

где

Кхх (т) определена

согласно

(II 1.57),

а у

(т) —

некоторая

новая

функция,

которая в отличие от

функции

КІх

(т)

может

иметь ненулевое значение при всех значениях т, как положи­

тельных,

так

и

отрицательных.

 

 

 

 

 

Подставив

(III.56)

 

и (III.58)

в

уравнение

(III.55а),

найдем

 

 

 

со

 

Г

со

 

 

 

 

 

 

 

 

 

\

P(h)

 

j

 

 

K7x{t2)Kt{x-k-t2)dt2-

 

 

 

— J

 

KJx{t2)y{x~t2)dt2

dh

= 0,

т ^ О .

 

Изменив

порядок

интегрирования,

получим

 

со

Г

со

 

 

 

 

 

 

 

 

 

 

j Кхх (k)

 

]

Р (h) Ktx

( x - h -

t2)

dti

-

V ( T -

h) dt2 =

0, T 5>0.

—oo

[_—со

 

 

 

 

 

 

 

 

 

126

Поскольку согласно (III.57) множитель Kl'x (t2) отличен от

нуля при t2

«S 0, то множитель, заключенный в квадратные скобки,

должен быть равен

нулю

при т ^

0

и

f2 <

0,

т. е.

 

со

 

 

 

 

 

 

 

 

 

 

 

\

P(h)Ktx(r

h — t2)dh

— y{x — t2)

=

0,

т ^ О ,

t2^0.

(111.59)

СО

 

 

 

 

 

 

 

 

 

 

 

Очевидно,

что

это

равенство справедливо

при

 

 

 

 

 

 

 

 

т — t2

 

О-

 

 

 

 

 

 

Обозначив

т —

t% = t, перепишем

(III.59)

в

виде

 

 

 

со

 

 

 

 

 

 

 

 

 

 

 

J

Pit^Ktxit-t^d^

Y(0 =

0,

t^O,

(111.60)

 

 

—со

 

 

 

 

 

 

 

 

 

По форме уравнение (III.60) напоминает исходное уравнение Винера—Хопфа (III.55а), однако имеет одно'существенное отли­

чие: первый член его левой части равен нулю при t ^

0. Действи­

тельно, Ktx (t)

=

0 при t

< 0 по определению и функция Р

(tj

=

= 0 при

при

t

< 0

по

условию

физической осуществимости

си­

стемы.

 

 

 

 

 

 

 

 

 

 

 

 

t < 0

 

 

 

Таким

образом,

левая

часть

(III.60)

при

может

быть

отлична от нуля только за счет члена у

(t), в то время как

левая

часть (III.55а) могла бы быть отличной от

нуля

при

т < 0 за

счет

обоих членов.

 

 

 

 

 

 

 

у

(t)

 

 

 

 

 

 

Представим

 

теперь

функцию

в

виде

 

 

 

 

причем

 

 

 

v

(0

=

у+

(0

+

v. (0,

 

 

(HI.61)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y+(t)

- 0 ,

t

< 0 ;

 

 

 

 

 

 

На основании (III.61) и (III.62) можно сделать вывод, что левая часть (III.60) отлична от нуля при t < 0 только за счет члена у (t) и, следовательно,

со

j P(tl)Ktx(t-tl)dh-y+(t)

= 0.

(111.63)

—со

Уравнение (III.63) справедливо при всех значениях времени t и является обычным уравнением Фредгольма первого рода. Такое уравнение, как было сказано, может быть решено в результате преобразования Фурье, т. е. путем перехода в комплексную область.

Напомним, что прямое преобразование Фурье функции вре­ мени / (t) определяется как

со

 

F(s) = j

f{t)t-«dt,

— со

а обратное преобразование Фурье — как

Из определений следует, что прямое преобразование Фурье функции времени, равной нулю при отрицательных значениях времени, имеет все полюсы, расположенные в левой полуплоскости. Прямое преобразование Фурье функции времени, равной нулю при положительных значениях времени, имеет все полюсы, рас­ положенные в правой полуплоскости. Прямое преобразование Фурье функции, отличной от нуля как при положительных, так и при отрицательных значениях времени, имеет полюсы, располо­

женные

как в левой, так и в правой

полуплоскости.

 

 

Перейдем к уравнению (III.63).

 

 

 

 

 

 

Преобразованием Фурье корреляционной функции Кхх

(t)

является

спектральная плотность Фхх

(s), преобразованием

Фурье

функции

веса Р

(t) — передаточная

функция

F

(s), преобразова­

нием Фурье функции у (t) — функция

Г (s).

Применяя преобра­

зование

Фурье

ко всем

членам уравнения

(III.63), получим

 

 

 

F (8)Ф+ ( s ) - r + ( s )

= 0,

 

 

(111.64)

где ФІх

(s) и Г +

(s) — преобразования

Фурье

функций

Ktx

и v+ (0

соответственно.

 

 

 

 

 

 

 

Из уравнения (II 1.64)

следует

 

 

 

 

 

 

 

 

f

( s ) = ^ P

 

 

 

 

( Ш - 6 5 )

что является окончательным выражением для передаточной функ­ ции искомой системы. Теперь остается установить способ опре­ деления функций Фtx (s) и Г + (s).

Из выражения

со

Kxx(t) = J

IQx{x)Ktx{t-%)dx

 

—оо

 

 

получаем

 

 

Фхх (s) =

Фtx (S) Ф7х (s).

(Ш.66)

Из определения функций Ktx (t) и К7х (t) свойств преобразо­ ваний Фурье и из того, что функции (s) и Ф7Х (s) являются преобразованиями Фурье функций Ktx (t) и К7Х (0. следует, что функция Ф ^ (s) имеет все полюсы в левой полуплоскости и не содержит полюсов в правой полуплоскости. Передаточная функ­ ция F (s) является преобразованием Фурье функции веса Р (t), равной нулю при t < 0. Поэтому все полюсы функции F (s) должны находиться в левой полуплоскости, что совпадает с требо-

ванием устойчивости системы. Поскольку в выражении (III.65), определяющем функцию F (s), функция Ofx (s) находится в зна­ менателе, то и все нули этой функции должны находиться в левой полуплоскости.

Таким образом,

функция

(s) представляет собой

множи­

тель, содержащийся

в Oxx(s)

и включающий все полюсы

и нули

функции Ф**^), расположенные в левой полуплоскости. Так как по формуле (III.66)

ф - (s)

то функция Ф~х~х (s) представляет собой второй множитель, содер­ жащийся в (s) и включающий все полюсы и нули функции Фхх (s)> расположенные в правой полуплоскости. Как уже ука­ зывалось, такое разложение четной функции на множители, обла­ дающие указанным свойством, называется факторизацией.

Из выражения

со

— со

вытекает, что

Ф*г(8) = Ф7*Г(8),

(III.67)

где ФХ2 (s) — преобразование Фурье взаимной корреляционной функции, т. е. взаимная спектральная плотность.

• По формуле (III.67)

Y(s)= Ф*г{5) .

Ф7*(*)

Эта функция является преобразованием Фурье функции у (t). Функция

 

 

 

Фхг

(S)

(111.68)

 

 

 

 

 

 

 

 

<-Фх~х (s) +

 

 

является преобразованием Фурье функции у

+ (t).

 

Из выражения

у (0

= у+

 

v_ (t)

 

 

 

(t) +

 

 

следует

 

 

 

 

 

 

T(s)

=

r + (s)

+ r_(s).

(111.69)

Эта операция носит название сепарации.

 

 

Так как функция Y+ (0 равна нулю при t < 0 , а функция v_

(t)

равна нулю при 4 >

0, то из (III.69)

следует,

что функция Г +

(s)

представляет собой слагаемое (компоненту, составляющую) функ­ ции Г (s), содержащее все полюсы в левой полуплоскости.

9 Я . М. Цейтлин

129

Соседние файлы в папке книги из ГПНТБ