Файл: Френкс, Л. Теория сигналов.pdf

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

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

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

Добавлен: 16.10.2024

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

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

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

В (9.62) мы использовали тот факт, что Рх (0) — 1, Рх (/) — преоб­ разование Фурье от р х (£).

При большом отношении сигнал/шум мы можем пренебречь Кии (/) в (9-58), и характеристика оптимального фильтра принимает вид

(9,63)

l°<n[! 1- ~ [ р,(я +6, (-т

Рис. 9.9. Характеристика оптимального фильтра, компенсирующего искаже­ ния за счет считывающего окна случайной ширины.

Для ширины окна, распределенной равномерно в интервале 0,5 т0 < 1,5 т0, оптимальная характеристика приведена на рис. 9.9» Здесь же для сравнения показаны характеристики компенсаторов

(инверсных фильтров) для усредненного канала G (/) и для канала GTo (f), соответствующего среднему значению ширины окна. Отличия, связанные с флюктуациями, особенно ощутимы при высоких частотах.

Упражнение 9.5. В задаче о случайном канале (см. рис. 9.7) принять, что канал имеет характеристику идеального фильтра низких частот с частотой

среза

W, так что

 

 

 

 

 

 

G (/) = - ~

[1 +sign (W —1/|)],

но величина W случайна и равномерно распределена в интервале 0 < Д, < W <

<

Полагая

K u u ( t )

= 0,01,

найти и нарисовать оптимальную харак-

***(/)

теристику Н (f) для оценки формы стационарного в широком смысле сигнала х. Упражнение 9.6. Найти оптимальный фильтр для оценки формы сигнала

в канале с мультипликативной помехой, аддитивным шумом и случайной пере­ даточной функцией (инвариантной во времени), т. е.

Т0 : to (0 = х (t — T),

СО

Tx:z(t)= j g(t — т) w(x)x(T)dT + u(0,

—00

где w (t), g (t), u (t) и x (t) — статистически независимые величины.

241


9.4.ОЦЕНКА ИМПУЛЬСНЫХ АМПЛИТУД

Всистемах типа радиолокации или АИМ, когда полезная инфор­ мация имеется только в амплитуде импульсов, применять устройства, оптимальные в смысле непрерывной оценки формы сигнала, может ка­ заться нецелесообразно. Действительно, такой путь приводит к квазиоптимальным решениям, поскольку строго оптимальным является устройство, минимизирующее средний квадрат ошибки оценки ам­

плитуды.

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

x(t)=a.f(t)

 

z (t)

M(f)

 

G(f)

0

-

 

Огпсчетный

Входной

 

 

импульсfit)

 

и. ft)

 

 

момент

со случайной

 

( f )

t - t n

амплитудой а

 

К

2]

 

ии

J

/

 

 

 

 

^

1 | у 0 ) - < 4 J

Рис. 9.10. Оценка амплитуды импульса.

предполагается непараметрическим, т. е. с постоянными парамет­ рами. Это допущение не снижает общности, поскольку представляет интерес отклик фильтра только в момент t = t0. Структура выходного сигнала фильтра в другие моменты не имеет значения, так что нет не­ обходимости рассматривать фильтр с переменными параметрами. Итак,

Г0:(о(г0) = а

(9.64)

и условие оптимальности (9.4) для h (t) принимает вид

 

оо

 

J kzz(s, о) h(t0o)do = kz<iJ (s, t0).

(9.65)

—oo

Хотя в данном случае входной процесс не стационарен, как это имело место в предыдущих задачах, в решении (9.65) все же удобно использовать частотные представления. Эквивалентное (9.65) условие в частотной записи можно получить, если воспользоваться подста­ новкой

ОО

 

h(t0 — о)= J H *(v)e-i2m{t° - a) dv,

(9.66)

---оо

изатем применить преобразование Фурье по переменной s. Эго дает

оо

 

j Я* (v)e~,2nvt°Вzz (/, v)dv — Cza(f, t0),

(9.67)

242


где обозначено

со

B z z ( f , V ) =

j I k z A

s - c r ) Q —

i 2 z i ( f s — v a ) d s d e ,

(9.68)

 

—oo

 

 

 

 

CO

 

 

с гм (/,

t0) = j

fc2(B(s,

(0) e-W*l* ds.

(9.69)

Согласованный фильтр

Рассмотрим классическую задачу оценки амплитуды импульса известной формы при известном времени прихода и при наличии аддитивного, стационарного в широком смысле шума с нулевым сред­ ним. Из рис. 9.10 видно, что

оо

: 2 (/) = a ^f{t — a)g(a)do + u(t) = ar(t) + u(t).

(9.70)

— со

 

 

Вычислив корреляционные функции

 

 

kzz (s, о) — Е г (s) ar (а)] + kuu (s—а) =

aV (s) г (а) -\-kuu (s~o),

kz®(s, t0) = E [aar (s)] =

aV (s)

(9.71)

и использовав (9.68), (9.69), найдем

Bzz(f, v) = t f R (/) R* (v) + Kuu (/) 6 (f —v), Сги(Д t0) = a2R(f).

Условие (9.67) для функции передачи оптимального фильтра при­ нимает вид

оо

a?R(f) J R*(\) ff*(v)e~2nv(»dv-l-H*(/)K'uu(/)e~ /'2lTB|>==a2£((/). (9.72)

Вводя неизвестный пока коэффициент пропорциональности а, полу­ чаем решение в виде

• W(f)

gfl* (В e~l2rlfl°

(9.73)

Каи (/)

 

 

причем подстановка (9.73) обратно в (9.72) дает

а —

аа

(9.74)

I R (/) I2 df

1+а2

 

 

Каи (/)

 

Мы получили согласованный фильтр для сигнала r(t), принимаемого в шуме со спектральной плотностью Кии (/) [4—6]. Этот фильтр есть обобщение согласованного фильтра, полученного в примере 6.3. Там

243


мы нашли фильтр, максимизирующий выходной отклик в момент /= ^ 0, при входном сигнале единичной энергии. Можно было предполагать, что такой фильтр окажется оптимальным для оценки амплитуды сиг­ нала. Для белого шума это действительно так, но, как следует из (9.73), если спектр шума неравномерен по частоте, можно реализовать дополнительный выигрыш за счет большего усиления на тех частотах, где плотность шума мала. Реализация такой возможности для одного частного случая рассмотрена в следующем примере.

b(t)

Рис. 9.11. Реализация согласованного фильтра для треугольного импульса в шуме с экспоненциальной автокорреляцией.

Пример 9.6. Треугольный импульс в шуме с экспоненциальной корреляционной функцией. Пусть г (t) — треугольный импульс, по­ казанный на рис. 9.11, а аддитивный шум имеет спектральную плот­ ность вида

Kuu(f)

__N0__

(9.75)

1 +

(f//o)3

 

 

Согласно (9.73) функция передачи согласованного фильтра имеет вид

sin2 nTf c- j 2nft0

(9.76)

(я/)2

244

Положив tQТ и представив sin nTf через экспоненциальные функ­ ции, можно получить

Я ( / ) = ~

/2л/

2я/0

(1—2е—/2л?7

+ е' - /2jt(27')f ). (9.77)

Л^о

 

 

Отсюда непосредственно

следует реализация фильтра, показанная

на рис. 9.11. Согласно (9.8) получаемый средний квадрат ошибки оцен­

ки

амплитуды

имеет значение

 

 

 

 

 

 

 

оо

 

Imin = km>(t0, /„) —

Г kza(s, t0) h(t0~s) ds= а2 —а2 а Г

^ df

 

 

 

 

d

 

J

(/)

 

 

------------ .

(9.78)

 

 

 

 

—oo*««,(/)

 

 

 

 

 

Для сравнения

рассмотрим

 

 

случай, когда для фильтрации

 

 

этого

сигнала

используется

 

 

фильтр,

согласованный при бе­

 

 

лом шуме. Такой фильтр

полу­

 

 

чается из рис.

9.11

удалением

 

 

блока (1/2 я /0)2

и подстановкой

 

 

оптимального значения а.

Если

Рис. 9.12. Сравнение оптимального и ква-

это

проделать,

средний

квад­

зиоптимального согласованных фильтров

рат ошибки (см. упражнение

для больших отношений сигнал/шум.

9.9)

принимает значение

 

 

 

 

 

 

 

Iw = '

 

а-

(9.79)

 

 

 

 

J| R(f) 12d/]2

 

 

 

 

 

 

 

 

 

 

 

1+ J IЯ tf) I* *«.(/) df

 

Для

сравнительно больших отношений сигнал/шум можно получить

 

 

 

 

1R (/) I2

 

 

 

 

'W

Каи (П■df

[j I R (f) I2 Kau (f) df]

(9.80)

 

 

 

 

 

 

 

[jltf(f) I2 df]s

Это отношение показано на рис. 9.12 как функция параметра с = 2nf0T. Ясно, что выигрыш за счет строго оптимального согласованного филь­ тра ощутим лишь в случаях, когда полоса шума равна или меньше ширины полосы сигнала.

Упражнение 9.7. Выяснитть структуру согласованного фильтра (рис. 9.10) при дополнительном ограничении на площадь усиления

Д = \H(i)\*dfJ .

—ОО

245


Упражнение 9.8. Показать, что фильтр (см. рис. 9.10), который максими­ зирует пик импульсного отклика 1г — Е \у (/„)], при фиксированной мощности шума на выходе

/2== f KuU( f)\HWdf

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

Упражнение 9.9. Показать, что если согласованный фильтр для импульса г (I) в белом шуме

H(f)=aR* (/) е~>'2п1/°

используется с шумом иной плотности, то минимальный средний квадрат ошибки оценки амплитуды определяется выражением

W'

^fj |Д(/л»<*/]2

+ -

 

i Кии (/)! R (Dl2 df

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

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

Если форма принимаемых импульсов точно не известна, то сог­ ласованный фильтр будет, строго говоря, квазиоптимальным. Будем рассматривать г (t) как реализацию некоторого случайного процесса и найдем фильтр, минимизирующий усредненный по реализациям сред­ ний квадрат ошибки оценки амплитуды. Можно свести эту задачу к фильтрации сигнала / (t) известной формы при наличии случайной дисперсионное™ канала G (/), т. е. к задаче, рассмотренной ранее в связи с оценкой формы сигнала. Получается та же задача, что на рис. 9.10, но G (/) следует считать случайной. Итак:

То:® (О = а-

 

оо

 

Tt :z(t) = a. j / ( / - i ) g ( S ) d | + u(*).

(9.81)

 

— оо

 

Нужные корреляционные функции имеют вид

 

kzz(s, сг) = a2 §f ( s — l)f(o— i\)E[g(l)g(i\)]dtdn + kuu(s— o),

 

_

«

(9-82)

/Gc„(s, /„)~ a2

Z)Elg(t)]dZ.

 

— ОО

246