книги из ГПНТБ / Вайнштейн Л.А. Выделение сигналов на фоне случайных помех
.pdfСуммарный процесс (18.02) в высокочастотном |
канале |
||
можно также представить в |
виде: |
|
|
|
/ (0 = ef cos (uot—&f), |
(18.08) |
|
где |
___________________________ |
|
|
ef |
= V e2m + |
cos (&„г - &„) |
(18.09; |
есть огибающая суммарного процесса (сигнал-ф-помеха). 'Ге же формулы годятся и в канале промежуточной частоты а>0.
Детектирование заключается в нелинейном преобразо вании сигнала (например, -с помощью «линейного» или квадратичного детектора) и в пропускании через линей ный фильтр, отсеивающий высокочастотные составляющие. В дальнейшем под детектором мы будем понимать сово купность нелинейного устройства и фильтра. На выходе «линейного» детектора получается огибающая (18.09) входного процесса f(t), на выходе квадратичного детекто ра — ее квадрат
При отсутствии помехи эти детекторы выдают функции ет и ет' ПРИ отсутствии сигнала—функции еп и еп2 соот
ветственно. Однако при наличии как сигнала, так и по мехи результирующая функция на выходе детектора не сводится к сумме их „индивидуальных вкладов11. В самом деле, в выражении (18.10) наряду с суммой ет2 -\-е2п фигу
рирует член ^етеп cos (0 —& ), имеющий смешанный ха
рактер („интерференционный11, поскольку он определяется разностью фаз сигнала и помехи) и появляющийся вслед ствие нелинейности преобразования. Такие же члены
имеются и в выражении (18.09). В самом деле, будем счи тать помеху гораздо сильнее сигнала
Пользуясь формулой |
е |
<^е ■ |
|
|
(18.11) |
|
|
|
|
|
|
||
j/1 -|~х=1 |
х---- yq х2, |
(18.12) |
||||
мы можем выражение (18.09) представить в виде |
|
|||||
ef = en + em cos (&„-&„) + |
|
|
||||
+ |
е2 |
I |
|
1 |
, |
(18.13) |
|
1 ~ 4c°s2(’\-’\) |
|
||||
11Р
где второй член ет cos (&;г — &л), а также частично и по
следний член правой части, имеют интерференционный ха рактер, т. е. зависят от
При теоретическом анализе детектирования интерфе ренционные члены можно считать помехой. Тогда формулы (18.10) и (18.13) сразу позволяют рассчитывать интересный и важный эффект: подавление слабого сигнала сильной по мехой при детектировании их смеси. В самом деле, отно шение сигнала к помехе (по амплитуде) до детектирова
ния равно
г' = ^. |
(18.14) |
еп |
|
Если оно мало, то после квадратичного детектирования по формуле (18.10) будем иметь
|
е2 |
|
(18.15) |
|
г"=^- =(г')2 < г’, |
||
|
еп |
|
|
а после |
линейного детектора |
в силу формулы (18.13) |
|
|
е2 |
/ г> \з |
<18-16) |
|
|
<г'- |
|
Наоборот, слабая помеха |
при детектировании |
подав |
|
ляется |
сильным сигналом. В |
самом деле, если величина |
|
(18.14) велика, то на выходе квадратичного детектора бу дем иметь
|
|
е2 |
|
(18.17) |
г" = ^ = (г')2> г'. |
|
|||
При условии |
|
еп |
|
|
|
|
|
(18-18) |
|
|
|
|
|
|
вместо формулы (18.13) |
получаем |
|
|
|
ef = ^ + encos |
+ |
|
||
е2 |
Г |
1 |
. |
(18.19) |
+iT |
1-TTCOS 2 |
|||
tn |
L |
|
|
|
и на выходе линейного детектора |
|
|
||
|
|
4е2 |
|
(18.20) |
г"=—^ = (2Н)а>Н. |
|
|||
еп
120
При выводе формул (18.17) и (18.20) интерференционные члены не учитывались; предполагается, что их можно уни-
чтожить путем некоторого усреднения величин ef9 |
или ef |
|
на выходе детектора. |
|
|
Выражения (18.15) и (18.16) показывают, что при ма |
||
лых отношениях |
сигнал/помеха квадратичный |
детектор |
лучше линейного |
(он дает вчетверо большее значение г"). |
|
Наборот, при больших отношениях сигнал/помеха линей ный детектор лучше, поскольку он дает вчетверо большее значение г", чем квадратичный.
Выше мы считали, что детектор тем лучше, чем больше отношение сигнал/помеха на его выходе. Следует отметить,
что определение этого важного параметра при нелинейных преобразованиях суммы сигнала и помехи не вполне одно значно, поскольку интерференционные члены содержат какую-то информацию о сигнале и их можно причислять как к помехе, так и к сигналу. При линейных преобразова ниях неопределенности нет, поэтому отношение сигнал/по меха их характеризует более полно.
Полученный выше эффект подавления слабых сигналов показывает, что в видеоканале корреляционный метод вы деления сигналов (т. е. стробирование и накопление) дол жен быть менее эффективен, чем в канале высокой или промежуточной частоты. Задачу о выделении когерентной последовательности радиоимпульсов в этих каналах мы рассмотрим в § 20, а в следующем параграфе подготовим необходимый математический аппарат.
Заметим, что под коррелятором в литературе часто подразумевают прибор, образовывающий автокорреля
ционную функцию входного процесса /(0> |
т. е- величину |
00 |
|
^(Т)= p(Of(^-x)d/, |
(18.21) |
—00 |
|
где под f(t) следует понимать, например, |
сумму сигнала |
и помехи (18.02) за L периодов повторения. В видеокана ле радиолокационного приемника коррелятор (18.21) бес полезен, поскольку даже в отсутствие помехи n\t) он вы дает функцию /?.„(т), изображенную на рис. 19,а: она имеет максимум при т = 0 безотносительно к тому, с каким запозданием пришел полезный сигнал (потеря дальности цели) и имеется ли один или два полезных сигнала с раз личными временами запаздывания (потеря разрешающей
121
способности). Можно также показать, что при сильной по мехе коррелятор, работающий по формуле (18.2'1), приво дит к дополнительному подавлению слабого сигнала поме хой. Это подавление имеет тот же характер, что и при квадратичном детектировании.
Если в формуле (18.21) производить интегрирование по конечному промежутку времени и под f(t) понимать сум му высокочастотного сигнала и помехи, то получится «те кущая» автокорреляционная функция. Можно показать, что прибор, образовывающий такую функцию, лишь не
значительно отличается от квадратичного детектора, в ко торый он переходит при т = 0.
Выше мы говорили, что с помощью отношения сиг нал/помеха трудно оценить нелинейную обработку приня тых данных [например, по формулам (18.09), (18.10) или
(18.21)]. Поэтому оптимальные способы нелинейной обра ботки приходится искать, опираясь на более тонкую ста тистическую теорию, изложенную во второй части книги.
§ «9. ТЕОРЕМА КОТЕЛЬНИКОВА И СОПРЯЖЕННАЯ ТЕОРЕМА
Теорема Котельникова дает представление сигнала с ограниченным спектром. Возьмем интеграл Фурье
|
/'(0 = i ^emtF(w)du>, |
|
(19.01) |
||
где |
—ОО |
|
|
|
|
|
|
|
|
|
|
|
ОО |
f (/) di. |
|
|
|
|
F(co) = J |
|
(19.02) |
||
|
—00 |
|
|
|
|
Будем считать, что |
|
|
|
|
|
|
F (со) = О при со<—со |
и со |
со. |
(19.03) |
|
Это значит, |
что спектральная |
амплитуда |
F (со) отлична от |
||
нуля лишь в |
пределах—со <[ со <[со, |
и, |
следовательно, в |
||
этом интервале ее можно разложить в ряд Фурье, |
который |
||||
мы берем в |
комплексной форме |
|
|
|
|
|
|
. |
иаа> |
|
|
|
|
I--- |
|
|
|
|
|
Dае |
«> |
|
(19.04) |
1 22
Коэффициенты ряда Фурье будут равны
~. теасо
|
Da = ~ fe |
|
|
(19.05) |
|
2co ~ |
|
|
|
|
— CO |
|
|
|
Сравнивая выражение (19.05) |
с формулой (19.01), полу |
|||
чаем |
— -^Л |
|
(19.06) |
|
|
|
|||
|
ш \ |
'» / |
|
|
или, положив для сокращения |
|
|
|
|
|
Д^ = -^-, |
|
(19.07) |
|
|
<0 |
|
|
|
мы можем переписать формулы (19.06) и (19.04) в |
виде |
|||
|
Ра = Д^7(—аДО, |
|
|
|
F(v) = M V f (- аДО ег'ашА/ |
= |
} |
(19.08) |
|
|
|
|
||
|
Оо |
|
|
|
|
а=—со |
|
) |
|
Подставляя это выражение в формулу (19.01) и учиты |
||||
вая условие |
(19.03), мы получим |
|
|
|
|
2 |
i<D |
(t—a&t) |
|
|
J e |
du)== |
|
|
|
a——co |
~ |
|
|
|
— co |
|
|
|
|
oo |
~ |
|
|
|
= У )(аД0 4р- |
Sinr~ZvA<) |
(19.09) |
|
|
a=—OO |
|
|
|
ИЛИ |
|
|
|
|
|
co |
~ |
|
|
|
/(/) = у f (aAZ) |
(z ~ a^-L . |
(19.10) |
|
|
a=% |
|
|
|
Это есть |
теорема Котельникова. |
Она показывает, что |
||
если сигнал |
имеет ограниченный спектр, то его |
значение |
||
12.3
в любой момент |
времени определяется |
теми |
значениями, |
|||||
которые он принимает |
в моменты |
= аД£, |
отстоящие друг |
|||||
от доуга на интервал (19.07): |
|
суммы |
(19.10) в мо |
|||||
Следует отметить, что значение |
||||||||
менты /я = аД^ определяется только |
одним а-м слагаемым, |
|||||||
а все остальные слагаемые дают нули, |
поскольку выраже- |
|||||||
sinw(a — В)Д/ |
х |
|
|
/о |
|
|
||
ние ■ ~ |
■—-— обращается в нуль |
при а^р. |
||||||
ш (а — (J) Д/ |
|
|
|
|
|
|
|
|
Имеется и доугая теооема, которую |
мы будем назы |
|||||||
вать сопряженной теоремой Котельникова. |
Она относится |
|||||||
к сигналу конечной длительности, т. |
е. |
предполагает, что |
||||||
|
7(0 = 0 |
при |
— |
|
|
|
(19.11) |
|
Тогда функцию f(f) в |
интервале—можно пред |
|||||||
ставить с |
помощью ряда Фурье |
|
|
|
|
|
||
. тса/
оо1 ~
|
Ш = £ С„е |
‘ ■ |
(19.12) |
где |
а=—оо |
|
|
. nai |
|
|
|
~ |
|
|
|
t |
~1 и |
|
(19.13) |
С =4 fe |
|
|
|
2t |
2t \t |
J |
|
—t |
|
|
|
И |
|
|
|
|
A(d==-^-, |
|
(19.14) |
в виде |
t |
|
|
|
|
|
оо
7— —00
Подставляя это выэажэние в фэомулу (19.02) и учиты
вая условие (19.11), получим выэажение |
для комплексной |
||
спектральной амплитуды F (уь) |
сигнала |
конечной |
длитель |
ности |
|
|
|
00 |
|
~ |
|
F(a>) — V F (аДсо) |
sin( ----a.A2))A . |
(19.16) |
|
«=-00 |
(<о —ЯДо>)< |
|
|
124
Это и есть сопряженная теорема Котельникова. Она показывает, что для сигнала конечной длительности зна чение спектральной амплитуды при любой частоте одно значно определяется теми ее значениями, которые она при нимает при ш=аД<.о.
Отметим, что формулу (19.10) можно |
|
применить к |
||||||
функции |
|
|
|
|
|
|
|
|
|
Ш = |
+ |
|
|
|
|
(19.17) |
|
где t0 — произвольный момент времени, |
и |
мы получаем |
||||||
|
оо |
|
|
~ |
|
|
|
|
+ |
V /А + ^0 |
|
|
|
. |
(19.18) |
||
|
|
|
ш(/-аДП |
|
|
|||
Здесь fv как и |
f, есть |
любая функция с |
ограниченным |
|||||
спектром. Поэтому, обозначая |
через |
t |
и |
беря |
вместо |
|||
просто f, мы получаем соотношение |
|
|
|
|
|
|||
|
ОО |
|
sin со (t — t0—aht) |
|
||||
|
|
|
(19.19) |
|||||
|
|
|
co (t — ta — аДО |
|
||||
|
а=—оо |
|
|
|
||||
|
|
|
|
|
|
|
|
|
несколько обобщающее выражение (19.10). |
|
|
применима к |
|||||
Теорема Котельникова (19.10) и (19.19) |
||||||||
таким функциям |
времени f (t), |
которые |
разлагаются в ин |
|||||
теграл Фурье (19.01). Для случайных функций f (I) |
вместо |
|||||||
спектральной амплитуды |
F (со) |
вводится |
(ср. |
§ 3) |
спект |
|||
ральная интенсивность Sf(«). Поэтому для случайных функ
ций (например, для помех) теорема Котельникова, хотя и
очевидна с физической точки зрения, строго говоря, нуж дается в особом доказательстве, которое мы сейчас и при
ведем.
Пусть спектральная интенсивность Sf(co) удовлетворяет условию
S (со) = о при СО <) — со И со^>со, |
(19.20) |
аналогичному условию (19.03). Тогда корреляционная функ ция
ОО |
, |
|
1 л |
1<DZ |
(19.21) |
= 2тг |
е Sf dm’ |
: 25
связанная с Sf(w) интегральным преобразованием Фурье
Sf(<») = Je~^Rf^)dv, |
(19.22) |
—00
очевидно, может быть представлена в виде, подобном
(19.10),
|
|
00 |
|
~ |
|
|
|
|
Rf(r) = V Rf (аД0 |
si'~ ‘ |
, |
(19.23) |
|
|
|
а=-оо |
|
со^-аДО |
|
|
или, |
как в |
формуле (19.19), |
|
|
|
|
|
|
00 |
|
• ~ |
|
|
|
^(0 = V /?.(/o + 3AZ)sit^e--z°-^O , |
(19.24) |
||||
|
|
Р=—00 |
“ (’’■ |
0ДО |
|
|
где |
t0 — произвольный |
момент |
времени. |
Действительно, |
||
при |
выводе |
формул (19.10) и |
(19.19) были использованы |
|||
лишь свойства интегрального преобразования Фурье и
условие (19.03). |
|
для самой случайной |
|
Чтобы доказать формулу (19.10) |
|||
функции f (t), обозначим |
|
|
|
x(t) = V f (аДО -si”<0 |
~ azV) |
(19.25) |
|
a=-oo |
|
|
|
и докажем, что средний квадрат разности |
|
||
s(f) = x(O—/(О |
(19.26) |
||
равен нулю. Действительно |
|
|
|
еДО = х2 (0 — 2х"(0Т(0 + /ДО, |
(19.27) |
||
и, учитывая, что |
|
|
|
^(0 = /(0/(^-0, |
(19.28) |
||
мы будем иметь |
|
|
|
хД0 = У1 V Rf (рД^ — аДг1) |
si£M^-aA0 |
, |
|
~~ |
— |
Z(t — Д0 |
|
|
|
|
(19.29) |
х (0 / (О =У Rf(t — аДг1) ■ si”^(z-aA0 |
(! 9.30) |
||
126
f\t) = Rf(QY |
(19.31) |
По формуле (19,24) получаем (при t0 = — aAt, |
т— t0 — t) |
V Rf(— аДЦ-рдг) -si2C°{t=/?^—аД/), (19.32)
T |
ш(/- Д0 |
|
поэтому |
|
|
V/?z (^ — аДО sl^M(/~aA°=7(Fff(f). (19.33) |
||
« |
«(/—аД/) |
|
Полагая, кроме того, в формуле (19.24) |
t0 =— t, т = 0 и |
|
учитывая четность функции |
будем иметь |
|
но = /? (0) =у R (-1 + рдо |
=До7й, |
|
V
«о-рдп
(19.34)
откуда и получаем искомое соотношение
Й0 = 0, (19.35)
доказывающее формулу (19.10) для случайного процесса с ограниченным спектром.
§ 20. ГРЕБЕНЧАТЫЙ ФИЛЬТР И СОГЛАСОВАННЫЙ ФИЛЬТР (КОРРЕЛЯТОР) ДЛЯ ПОСЛЕДОВАТЕЛЬНОСТИ РАДИОИМПУЛЬСОВ
В данном параграфе мы рассмотрим последовательность
прямоугольных радиоимпульсов. |
Возьмем |
сначала |
один им |
|||
пульс, заданный следующим образом: |
|
|
||||
т (() = A cos <ogt |
при |
о |
|
|
||
2 |
|
|
||||
|
|
|
|
(20.01) |
||
т (/) = 0 |
, |
7'0 |
г0 |
|||
|
||||||
при t<Z — -у и |
Г>-у- • |
|
||||
Здесь <о0 — несущая частота, |
А — амплитуда импульса и |
|||||
То — его |
длительность. |
Спектральная амплитуда 34 (ш) оди |
||||
ночного |
импульса равна |
|
То |
|
|
|
|
оо |
|
|
|
||
|
|
2 |
|
(20.02) |
||
Л1(ш)= J |
m(t)dt=A $ е~м cos ^tdt |
|||||
|
—ос |
|
_ То_ |
|
|
|
127
или
|
(ю — а>о)Го |
si n |
(“ + "o)7'o |
|||
М (cd) = А |
si n------- |
п-------- |
|
9------- |
(20.03) |
|
__________-_________ I-------------------- |
CD —|— -Cl)-----------0 |
|||||
|
W ---- |
(l>0 |
‘ |
|
||
Графическое изображение |
этой |
функции |
дано на рис. 20 |
|||
для положительных частот со. |
При |
со < 0 |
получим такую |
|||
же кривую, поскольку |
|
|
|
|
(20.04) |
|
|
/И(— cd) = M(d). |
|
|
|||
Эти соотношения |
мы уже использовали в |
§ |
17. |
|||
Рис. 20. Спектральная амплитуда М (<и) прямоугольного
|
радиоимпульса. |
|
||
Импульс (20.01) |
симметричен |
относительно момента |
||
/ = 0. Беря прямоугольный импульс более общего вида |
||||
т (t) = A cos (cd t — 9) |
при т — ^-<2 t < т + |
1 |
||
т (0 = 0 |
22 (20.05) |
|||
при других t, |
I |
|||
мы вместо выражения (20.03) получим |
|
|||
|
|
. (“~ "о)7о |
|
|
|
sin------- 5-------- |
|
||
М (Со) = Де~‘'шх |
____________ |
|||
— со0 |
I |
|||
|
|
|||
|
((О + шо) Го |
|
||
|
Sin-------9-------- |
|
||
е- I (а.от-9)Z |
(20.06) |
|||
|
ш 4* |
“о |
||
|
|
|||
Если прямоугольные импульсы возникают беспорядочно, то мы получаем случайный процесс со спектральной интенсив ностью
S (cd) = д | М (cd) |2, |
(20.07) |
128
