С учетом (8.27) нужные корреляционные функции имеют значения
kzz (s, о) = |
V ат V / (s~kT) f (a— kT— mT) + kuu (s— a), (9.102) |
|
m |
k |
|
|
|
|
|
kz№(S, iT) = |
f (s - k T ). |
(9.103) |
|
|
|
fe |
|
|
Здесь величины am — E [aka.h + m] |
определяют |
корреляцию |
последовательности |
амплитуд. |
Теперь, |
применяя преобразование |
Фурье, с учетом (9.68) и (9.69) находим |
|
|
Bzz (/, v) ^ 2 а ш £ |
F (/) е -/2я.г fF* (v) е/2я <*+«) rv + Кии {f) б (/ _ v)> |
m |
k |
|
|
|
(9.104) |
|
|
Сгю (/, iT) = Т (/) v а ._й е-/2«‘П. |
|
|
(9.105) |
|
|
|
k |
|
|
Рис. 9.14. Периодическая оценка амплитуд АИМ сигнала.
Как и в § 8.4, мы представим корреляцию амплитуд в частотной обла сти, введя периодическую функцию
= |
(9.106) |
Тогда (9.104) и (9.105) можно преобразовать к виду |
|
Вхг (/. v) = F (f) F* (v) М (v) ± 2 6 ( f - v - - L |
|
+ K „ „ ( / ) S ( f - v ) , |
(9.107) |
Cza(f,iT) = F(f)M(f)e-*”‘Tf, |
(9.108) |
где учтено, что
V |
p—j2nkTf |
г / v |
т j |
k |
|
ш,J V' |
|
|
(см. упражнение 4.4).
Подставив (9.107) и (9.108) в (9.67), получим уравнение для функ
ции передачи оптимального фильтра |
|
^ ( / ) |
( f — j ) н * [f — ~ ) + |
|
|
+ Kuu(f)H*(f) = F(f)M(f). |
(9.109) |
Замечая теперь, что М (f) периодическая, сделаем подстановку W(f) = = F* (f)Ii* (/). Тогда из (9.109) получается
|
Л р y ,w \^f___pi |
WJf) |
(9.110) |
|
R(f) |
|
|
|
|
где обозначено |
|
|
|
R { f ) = J r !i k - |
(9Л11) |
Разностное уравнение бесконечного порядка (9.110) допускает простое решение. Это связано с тем, что правая часть и первое слагае мое левой части есть периодические функции с периодом ИТ. Поэто-
H(f)
M f )
Рис. 9.15. Оптимальная фильтрация для синхронизированного АИМ сигнала.
Обозначены отклики на одиночный импульс единичной амплитуды.
му, функция W (f)/R (/) |
также |
периодична. |
Положим W (/) = |
= £?(/) Q (/), где Q (/) — периодическая |
функция. |
Тогда |
— |
1 |
( / — |
j = L(f)Q(f), |
|
т |
|
|
|
|
где мы обозначили |
|
|
|
|
|
L(f) =■ т М Н |
О - |
(9.112) |
|
Подставляя эти значения в (9.110), находим |
|
|
W(f) |
|
|
(9.113) |
|
|
\+L(f)M(f)' |
|
Окончательно, функция |
передачи фильтра имеет выражение |
Ж/) |
|
F* (/) |
M(f) |
(9.114) |
|
Kuu(f) |
Ll +L(f)M(f) |
|
|
|
Из этого выражения видно, |
что оптимальный фильтр |
представляет |
собой каскадное включение |
согласованного |
фильтра |
для |
сигнала |
/ (/), принимаемого в шуме |
со спектральной |
плотностью К ии (/), и |
фильтра, имеющего периодическую функцию передачи |
[L (/) |
и М (/) |
периодичны]. Если второй множитель в (9.114) аппроксимировать конечным рядом Фурье, то фильтр (с общим временем задержки пТ сек) может быть реализован на многоотводной линии задержки, как пока зано на рис. 9.15 [7—10].
Коэффициенты усиления в отводах есть коэффициенты ряда Фурье для D (/):
£>(/) = |
|
М(П |
■2 |
dk е/2 |
|
1 +L( f ) M( f ) |
|
|
2 |
dh d 2*ьт. |
dh = d. |
(9.115) |
|
k=—П |
|
|
|
Эффективность такого фильтра просто выражается через |
D (/). Если |
подставить (9.114) в (9.8), |
найдем |
|
|
|
|
|
со |
|
|
|
Лп*п = |
« о - 5 R (/) М (/) D (/) df. |
(9.116) |
Некоторое пояснение работы оптимального фильтра можно получить, рассматривая частный случай, когда отклик согласованного фильтра г (t) на одиночный импульс удовлетворяет условию Найквиста (6.123), в силу которого межсимвольные помехи устраняются. Из (9.112) мы имеем в этом случае L (/) = г (0). Поскольку межсимвольные помехи на выходе согласованного фильтра отсутствуют, можно ожидать, что часть фильтра с многоотводной линией задержки должна быть изъята, так как она снова создает такие помехи. Это верно в случае, если ам плитуды импульсов некоррелированы, т. е. М (/) = а 0 и D (/) = const. Но если корреляция в последовательности амплитуд имеет место, то оптимальный фильтр дает некоторый выигрыш за счет уменьшения влияния шума при небольших межсимвольных помехах. Этот эффект ощутим, если шум сравнительно велик.
В случае сильного сигнала мы имеем
|
|
|
|
|
|
|
D(f) |
M(f) |
1 |
(9.117) |
|
Н-г(0) |
M(f)^ г {0) |
|
|
|
И, если межсимвольные помехи в г (t) отсутствуют, фильтр на много отводной линии задержки должен быть исключен.
Упражнение 9.11. В задаче о периодической оценке амплитуд АИМ сигна ла предположим, что сигнал нестабилен во времени, как было изложено в § 8.5. В этом случае
|
|
|
х (t)=I>Hhf (t — h), |
|
|
|
|
|
k |
|
|
где |
случайный |
момент прихода определяется |
выражением |
= kT + 5^; |
k = Q, ± ] , ±2, |
. . . |
Пусть 8* — статистически |
независимые случайные вели |
чины с плотностью |
вероятности р (|). Получить условие оптимальности фильт |
ра |
Н (f) в форме, |
подобной (9.110). |
|
|
9.6. УСЛОВИЕ ФИЗИЧЕСКОЙ РЕАЛИЗУЕМОСТИ
При практической реализации рассмотренных нами оптимальных фильтров приходится учитывать различные ограничения, которые не были сформулированы в качестве предварительных условий в за даче оптимизации. Эти ограничения приводят к изменению характе ристик системы и уменьшают ожидаемый выигрыш. Такие условия, как, например, допустимая сложность или влияние неизбежных пара зитных связей, не являются принципиальными для проблемы филь трации. Они могут быть преодолены ценой увеличения затрат и неко торых конструкторских решений.
В отличие от этого, требование неопережающего отклика прин ципиально необходимо для физической реализации фильтра. В этом параграфе рассматривается физическая реализуемость в указанном смысле. Условие такой реализуемости можно сформулировать в виде ограничения на импульсную характеристику h (t\ s). Необходимо и достаточно, чтобы h (t; s) обращалась в нуль для всех t < s. Как и ограничение на длительность, встретившееся нам в задаче оптимиза ции сигнала (см. гл. 6), условие физической реализуемости фильтра обычно легче учесть прямым методом, а не с помощью множителя Лагранжа. Поэтому мы вводим дополнительное условие
|
h (t, |
s) = w (t — s) h (t, s), |
(9.118) |
где |
w (t) = |
1 |
для |
t ^ |
0; |
|
w (t) = 0 |
для |
t < |
0. |
Тогда в сформулированной в § 9.2 задаче о минимуме среднего квадрата ошибки изменяется функционал /. Мы получаем вместо (9.3)
со
/ = ^ w(t—s)w(t—<j)h(t, s)h(t, о) kzz (s, о) dsdo — |
|
—2 w(t— s)h(t,s)kza(s,t)ds + ke>a(t,t). |
(9.119) |
— 00
Вычислив градиент функционала (9.119) по переменной h (t, s) и приравняв его нулю, придем к необходимому условию минимума
со
^w(t ~s)w(t —a)h(t,a)kzz(s,a)da=w(t—s)kz<a(s,t). (9.120)
—оо
Учитывая значение (9.118) для w (/), перепишем (9.120) в виде
t |
|
^ h(t,o)kzz(s,o)do = kz(£>(s,t) при t ^ s . |
(9.121) |
—оо |
|
Это уравнение удается решить только в некоторых частных случаях. Если процессы z и о> совместно стационарны, так что kzz (s, о) и (s, t) можно представить соответственно в виде kzz (s — о) и
kzu, (s — T0 (9.121) превращается в уравнение Винера — Хопфа, достаточно хорошо исследованное. В этом случае
t
^ h(t — o)kzz(s— o)do = kza(s— t) при t ^ s
—оо
или, изменяя обозначения,
оо |
|
|
<\) h(x)kzz(t—x)dx = kaz(t) при t^O. |
(9.122) |
о |
|
|
Методы решения |
уравнения Винера — Хопфа, особенно для слу |
чаев, когда Kzz (/) и |
Кга (/) — рациональные функции от f, |
имеются |
во многих работах по теории систем связи и систем управления [11— 15]. Мы рассмотрим метод, предложенный Боде и Шенноном [16, 17], который не только снимает ограничения на рациональность функций, но дает физически прозрачную картину работы оптимального фильтра. Сначала заметим, что если принимаемый сигнал г (t) есть белый шум
с единичной дисперсией, так что |
kzz (s — а) = |
6 (s — а), то (9.122) |
имеет особенно простое решение: |
|
|
h{t) = K z(t) |
при t ^ O |
|
или, что то же, |
|
|
h(t)^w(t )kaz(t). |
(9.123) |
Следовательно, в этом частном случае оптимальный реализуемый фильтр имеет импульсную характеристику, такую же как без учета условия реализуемости, но усеченную при t — 0. Простое решение уравнения (9.122), предложенное Боде и Шенноном, основано на том, что в большинстве практически важных задач сигнал z можно пред варительно пропустить через физически реализуемый отбеливающий фильтр с постоянными параметрами, т. е. такой, что корреляция на его выходе будет соответствовать белому шуму единичной дисперсии. Тогда остается лишь включить последовательно с отбеливающим фи зически реализуемый фильтр, соответствующий рассмотренному слу чаю, как это показано на рис. 9.16.
В соответствии с (9.123) оптимальный фильтр, воздействующий на сигнал z, имеет импульсную характеристику
оо |
|
h2{t) = w{t)kaz{t) = w(t)(j h1(x)k<i>z(t-\-x)dx. |
(9.124) |
о |
|
Удобно обозначить через [X (/)]+ преобразование Фурье от |
функции |
w (i)x (0, т. е. от усеченной (произвольной) функции х (t). Тогда можно записать
я 2 (/) = [К<Гг (/)]+ = |
[Н\ (f) Kvz(/)]+. |
(9.125) |
Функцию передачи отбеливающего фильтра обычно находят с по |
мощью спектральной факторизации |
[11]. Если Kzz(f) |
можно пред- |