Файл: Пугачев, В. С. Основы статистической теории автоматических систем.pdf

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

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

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

Добавлен: 16.10.2024

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

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

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

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

 

b (I) =--

J

ср2 Ст) dx,

I

 

\

 

i

и

 

 

 

 

 

 

(0 = j Ф (Т) N (Т) dT-

}

 

^0

 

 

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

формулу

к

виду

 

 

t

 

 

 

Z (/) =

J ср (т) X (т) dx = Ub (l) + Zu, (/).

 

tо

 

 

 

( 10. 8)

(10.9)

Докажем, что преобразование входного сигнала, определяемое формулой (10.9), обеспечивает максимальное отношение сигнал/шум по сравнению с любой другой линейной системой при любом воз­ можном значении и случайной величины U. Если заменить вели­ чину U ее возможным значением и, то полезная часть сигнала

s (/) = ub (t).

Дисперсия шумовой части сигнала (10.9) на основании второй формулы (10.8) и известной формулы для дисперсии интеграла от белого шума определяется соотношением

L

 

Ош(/) = D {Zm(«1 = О j

<pa(T)dT = Gb (t),

(10.10)

tо

 

 

где G — интенсивность белого шума

N (t). Отношение сигнал/шум

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

квадратическому отклонению шумовой части этого сигнала.

На

осно­

вании предыдущих результатов это отношение определяется

фор­

мулой

 

 

I и I V ь (/)•

 

 

_

[s (0 I

_ Iи\ь (/)

( 10. 11)

7

КД7(Т)

УЩГ)

VG

 

 

Рассмотрим теперь фильтр, который выполняет те же операции, но умножает входной сигнал X (т) не на функцию ср (т), а на другую функцию ф (т). Полезная часть выходного сигнала такого фильтра и дисперсия шумовой части определяются формулами

i

Si (/) = и | ф (т) ср (т) dx, tо

t

Dm (t) = G j ф2(t) dx.


На основании известного неравенства Буияковского—Шварца

! 2 1 !

J 1|) (т) ср (т) dr ^ J ф2 (т) dr J ф2 (т) dr

имеем

 

 

 

 

 

 

 

/

I

 

 

 

I SI (О I ^ м У

Ь(/) [ ф2

т) d r .

 

 

 

 

^0

 

 

Отсюда

получаем следующую

оценку

отношения сигнал,'шум

на выходе

рассматриваемого

фильтра:

[ и I ■./

 

 

= I Sl I =

I* (01

b(i).

 

V D an(l)

 

 

V g

 

 

ф2 (т) dr

 

 

 

 

 

 

 

1

/ 4tо

 

 

Сравнив правую часть этого неравенства с выражением (10.11), видим, что у! ^ у. Таким образом, отношение сигнал/шум на выходе любого линейного фильтра не может превзойти такое же отношение на выходе фильтра, осуществляющего преобразование (10.9). Фильтр, осуществляющий это преобразование, обычно называют согласо­ ванным фильтром. Согласованный фильтр преобразует выходной сигнал таким образом, что в результате на выходе получается наи­ большее возможное отношение сигнал/шум.

Сравнивая выражение (10.9) для выходного сигнала согласован­ ного фильтра с выражением (10.7) входного сигнала системы, видим, что для воспроизведения полезного сигнала 5 (/) = Uф (/) без си­ стематической ошибки следует умножить выходной сигнал согласо­ ванного фильтра на ф (t)lb (t). Иными словами, следует подключить к выходу согласованного фильтра усилитель с переменным коэффи­ циентом усиления ф (t)lb (/). Выходной сигнал на выходе такой системы

Y(t) = Uy(t) + l ^ Z u, (().

( 10. 12)

Ошибка системы

Е (D = Y (/) - Ут (/) = У (0 - UФ (/) = Щ Zm(0.

Имея в виду, что математическое ожидание белого шума, а сле­ довательно, и математическое ожидание случайной функции Zm (t) равно нулю, а также принимая во внимание формулу (10.10), опре­ делим средний квадрат ошибки найденной оптимальной системы:

rf = М [Е2 (/)] = -УЩ- D [Zm(/)1 = G ^ ~ .

(10.13)

Легко понять, что найденная система обеспечивает минимум ди­ сперсии ошибки при нулевой систематической ошибке, или, что

17*

259


то же самое, условный минимум средней квадратической ошибки при нулевой систематической ошибке.

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

системы:

(10.14)

Y (/) = kZ (0 = kUb (0 + kZm {l).

Вычитая из этого выражения требуемый выходной сигнал YT (i) = = t/cp (/), найдем ошибку системы:

Е (0 = U [кЬ (/) — ср (/)] + kZm (I).

Отсюда, принимая во внимание формулу (10.10) для дисперсии ошибки шума на выходе согласованного фильтра, получаем следую­ щее выражение среднего квадрата ошибки системы:

11 = М [

Е2

(0 ] = уи ГП (0 — Ф (0 I2 + кЮЬ {().

(10.15)

Дифференцируя

эту формулу по к и приравнивая результат

нулю, получаем

уравнение для определения к:

 

-If- = 2у„ [kb (0 - ср (/)] Ь(0 + 2kGb (0 = 0.

Отсюда находим

Тиф (0

(10.16)

ЧиЪ (/) 4- G

Подставив найденное выражение к в формулу (10.15), находим минимальный средний квадрат ошибки оптимальной системы:

= ъМП + О = (0 Ф (0.

00-17)

где показана явно зависимость к от времени t. Сравнивая эту фор­ мулу с (10.13), видим, что найденная оптимальная система дает средний квадрат ошибки в

yub(() + G _

. .

G

УиЬ (/)

'

у„ь

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

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

260



Рис. 10.1. Оптимальная одно­

мерная система

 

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

соединение согласованного фильтра и усилителя

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

усиления X (t), определяемым формулой (10.16).

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

Изложенный метод дает возможность определять оптимальные системы как с бесконечной памятью,.так и с конечной памятью Т. В последнем случае следует во всех формулах заменить нижний предел интегрирования 10 на 1— Т.

В некоторых практических задачах форма полезного сигнала не­

известна,

но

известно

линейное

дифференциальное

уравнение,

определяющее полезный сигнал:

 

 

 

 

 

 

 

 

S

= a(1)S.

 

 

(10.18)

Интегрируя

это уравнение

при

случайном

начальном

условии

5 (/„) = U,

получаем

 

 

 

 

 

 

 

 

 

S (t) = U ехр | J

а (х) dx

 

 

 

 

 

 

 

 

I/о

 

 

 

 

Отсюда

видим, что

неизвестный

полезный

сигнал

ср (/)

связан

с коэффициентом а (t)

уравнения (10.18) формулами

 

 

 

 

Ф (0 =

ехр {J а (т) dTJ . а (*) =

 

 

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

запишем выражение для выходного сигнала

системы:

t

 

Y (t) = X (t) | ср (т) X (х) dx.

(10.19)

iО

 

Дифференцируя уравнение (10.19) по времени t и опуская для

краткости обозначение аргумента t, получим

 

t

 

У = к | ср (т) X (т) dx + А,ср (t) X (t).

(10.20)

Умножим это выражение на X и вычтем из него выражение (10.19), умноженное на X, в результате получим

X Y — XY = X2срХ, или У = - У + ХуХ.

(10.21)

261