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

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

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

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

Добавлен: 09.04.2024

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

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

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

540

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

уравнение (14.1.1) приводится к линейному алгебраическому уравнению. Решая это уравнение и заменяя букву о буквой т, находим

(14.1.4)

Таким образом, мы имеем выражение функции g (t, т) в интере­ сующем нас интервале наблюдения t Т ^ т ^ t. Вне этого интервала ее следует принять равной нулю, так как значения вход­ ного сигнала, предшествующие моменту t Т, система не долж­ на использовать (или не может использовать, если момент t Т совпадает с моментом начала работы системы), а при т > t весовая

функция g (t,

т)

равна нулю вследствие условия физической воз­

можности

системы.

вместо /

(t, т)

по очереди функции

Подставляя

в

(14.1.4)

К ух (t,

т),

фі (т),

. . .,

(т), мы получим соответственно функции

(/,

т),

gw (t,

т), . . .,

gW (t, т).

После

этого весовая функ­

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

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

Пусть известная форма полезного сигнала характеризуется функцией Фі (<)• Случайную амплитуду полезного сигнала обозначим через Ui. Тогда задача сведется к решению одного интегрального уравнения (13.6.2). Поль­

зуясь формулой (14.1.4), находим весовую

функцию ga >;

' g4> (t,

т)

Фі (г)

 

к

'

 

 

После этого вычисляем по второй формуле (13.6.6) единственную в данном случае величину Ьці

t

t

* ii= j

ga Ht,x) Фі(т)йт = -і- j <pl(x)dx.

t - т

t-T

Решение единственного в данном случае уравнения (13.6.7) дает функцию

«):

фі (t)

Яі (0

*11 + C11

 

где сц = i/уц — величина, обратная по отношению ко второму начальному моменту случайной амплитуды сигнала. Наконец, формула (13.6.4) дает для весовой функции оптимальной системы следующее выражение:

g(t, x) = Xl (t) g<!> (t, T) = Ф1 (t) Фі (T)

(*11 + С ц )

Выясним подробнее смысл оптимальных операций над входной функцией в данном случае. На основании формулы (13.5.34) выходная переменная


§ 14.1. СЛ У ЧА Й БЕ Л О ГО ШУМА НА ВХ О Д Е

541

оптимальной системы выражается формулой

 

t

t

W*(t)=

|t / i j

q:f(T)dT+ j <Pi(t)X(T)dT

 

t - T

t - T

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

Простота решения уравнения (14.1.1) в случае, когда X (t) представляет собой белый шум, приводит к следующей идее реше­ ния уравнения (14.1.1) в общем

случае,

когда

X (t) не является

 

1

 

1

 

 

1

белым шумом. Если найти линей­

Х<V

\

V(t)

1

1

ную систему,

преобразующую

Г г (()

1

w~(t,x)

è.(t,T) І

 

1

данную

случайную функцию

 

1

 

1

 

1

 

 

X (і), т. е. сумму нерегулярной

 

части

полезного

сигнала и по­

Рис. 14.1.1.

мехи,

в белый

шум V (t), а

 

затем найти оптимальную систему для случая белого шума V (t)

на

входе,

то

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

этих двух

систем

и

будет

представлять собой

оптимальную

линейную

систему

для случайной

функции X (t)

на входе (обведена пунктиром

на рис. 14.1.1). Эта идея лежит по существу в основе всех известных методов решения интегрального уравнения (14.1.1).

Таким образом, решение интегрального уравнения (14.1.1) сводится к задаче канонического представления случайной функ­ ции X (t). Действительно, мы знаем, что проблема канонического представления случайной функции есть проблема ее выражения через белый шум — дискретный в случае канонического разложе­ ния и непрерывный в случае интегрального канонического пред­ ставления. Вследствие этого очевидно, что теоретической основой, на которой может быть построена общая теория оптимальных сис­ тем, является метод канонических представлений случайных функ­ ций. Далее мы рассмотрим применение обоих типов канонических представлений, когда случайная функция X (<) выражается через дискретный белый шум, т. е. через последовательность некоррели­ рованных импульсов, и когда она выражается через непрерывный белый шум. Применяя интегральное каноническое представление,


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

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

Заметим, что формулой (14.1.4) можно пользоваться как при­ ближенной и в том случае, когда случайная функция X (t) не является белым шумом, но близка к белому шуму. А именно фор­ мула (14.1.4) дает приближенное решение уравнения (14.1.1), если интервал корреляции т0 случайной функции X (t) (т. е. максимальное значение разности | т — а |, при котором кор­ реляционная функция К х (т, о) существенно отличается от нуля)

достаточно мал для того, чтобы функцию / (t,

т)/ G (т) можно было

считать

практически постоянной в любом

интервале изменения

т

длины

2т0. Для доказательства подставим выражение (14.1.4)

в

левую

часть уравнения (14.1.1), оставив

пока функцию G (т)

неопределенной. Тогда, принимая во внимание, что в данном случае при любом значении о в интервале t — Т + т0 ^ а

^ t — т0 подынтегральная функция практически равна нулю вне интервала — т0, о + т0), и вынося за знак интеграла значение функции / (t, т)IG (т) в середине а этого интервала, будем иметь

1

а+то

J к х (т, а) f(Gt’{^ ) dx »

j Кх { т, а ) Щ ^ - йхж

t - T

CT-То

 

-To

Отсюда видно, что если определить функцию G (а) формулой

G(o)= тоj Аж(о+ т), o)dr\,

-То

то формула (14.1.4) даст приближенное решение уравнения (14.1.1).

§ 14.2. Определение оптимальной линейной системы в случае стационарных сигнала и помехи

и бесконечного интервала наблюдения

Применим изложенную общую идею для нахождения оптималь­ ной линейной системы в частном случае, когда требуемый выходной сигнал W и входной сигнал Z содержат только нерегулярные части, представляющие собой стационарные и стационарно связанные


Рпс. 14.2.1.

§ 14.2. СЛ У ЧА Й С ТА Ц И О Н А РН Ы Х СИГН АЛА И ПОМ ЕХИ

543

случайные функции времени:

 

Z(t) = X(t), W {t) — Y (i),

(14.2.1)

а интервал наблюдения бесконечен: Т = оо. Для нахождения сис­ темы, преобразующей стационарную случайную функцию X (t)

вбелый шум, согласно изложенному в § 7.6, достаточно выразить

ееспектральную плотность в виде

S* (ю) = I Ф (ш) р,

(14.2.2)

где Ф (ію) — функция, все нули и полюсы которой лежат

в верх­

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

ей 1/Ф (s) и будет системой, пре­ образующей случайную функ­ цию X (t) в белый шум V(t) с еди­ ничной спектральной плотностью и интенсивностью, равной 2л.

Обозначим через | (і, т) не­ известную весовую функцию оп­

тимальной системы для случая белого шума V на входе. Тогда опти­ мальной системой для входной случайной функции X будет после­ довательное соединение системы с передаточной функцией 1/Ф (s) и системы с весовой функцией \ (t, т) (рис. 14.2.1).

Для определения весовой функции £ (і, т) оптимальной линей­ ной системы для случая белого шума V на входе имеем в данном случае одно уравнение (13.6.1) при Т = оо, правая часть кото­ рого представляет собой взаимную корреляционную функцию

куѵ (t

er) требуемого

выходного

сигнала Y (t) системы и белого

шума V (а). Применяя для решения этого уравнения формулу

(14.1.4),

находим

 

 

 

 

 

 

т) =

^

k

vv( t - x )

( —о о с т ^ і ) .

(14.2.3)

При т >

t функцию

I

(і,

т) следует положить равной нулю, так

как мы ищем физически возможную оптимальную систему. Из (14.2.3) видно, что оптимальная линейная система для белого шума V на входе стационарна. Ее передаточная функция Q (s) на осно­ вании (2.4.6) определяется формулой

ОО

 

&{s) = - ^ \ k yAl)e - *di .

(14.2.4)

о

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


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

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

¥ ( S) =

Q(s)

(14.2.5)

 

Ф(*)

 

Для полного решения задачи осталось теперь только найти взаим­ ную корреляционную функцию к,ІѴ(|) требуемого выходного сиг­ нала системы У и белого шума V. Так как белый шум V есть резуль­ тат преобразования случайной функции X устойчивой стационар­ ной линейной системой с передаточной функцией 1/Ф (s), то для определения взаимной спектральной плотности syv (со) случайной функции У и белого шума V можно применить формулу (7.4.21). В результате получим

 

s y x И

(14.2.6)

 

s y v (© ) — Ф (—ісо)

 

 

Зная взаимную спектральную плотность случайных

функций У

и V, можно найти их взаимную корреляционную функцию по

известной формуле теории вероятностей:

 

оо

оо

 

ку,( ! ) = j ѵ

( ^ і й Ф = j

<14-2-7>

—oo

—oo

 

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

Подставляя выражение (14.2.7) в (14.2.4) и после этого подстав­ ляя полученное выражение в (14.2.5), получим следующую форму­ лу для передаточной функции оптимальной линейной системы:

 

 

оо

оо

 

 

¥ (s) =

„ *

[ e-SdE

J

(

Syx W - é^d\i .

(14.2.8)

W

2n®(s)

J

 

ф (_ щ)

 

Эта формула является основной формулой теории экстраполяции, интерполяции и сглаживания стационарных случайных функций, разработанной Винером [12].

Двойной интеграл в формуле (14.2.8) во многих случаях легко вычисляется. Так, например, этот интеграл легко вычисляется, если отношение s,/x (р)/Ф (— ір) является дробно-рациональной функцией. В этом случае, полагая в (14.2.8) s = ісо и выполняя интегрирование, получим следующую формулу для частотной характеристики оптимальной линейной системы (см. [53], § 129):

¥ (ісо)

1

|~ Sy x (Ш) 1

(14.2.9)

Ф (ісо) і_ф( —іш) J+’

 

 

где квадратные скобки со знаком + внизу озпачают, что рацио­ нальная дробь, заключенная в эти скобки, должна быть разложе­