
Литература / Гришин Ю.П., Ипатов В.П., Казаринов Ю.М. Радиотехнические системы (1990)
.pdfчто они пересекают аналогичные кривые для сигнала фиксированной амплитуды, соответствующие тем же зна чениям рЛ1. Объясняется это тем, _что эпизодические большие выбросы флуктуирующей амплитуды увеличивают вероятность обнаружения сигнала с малым значением q, в области же больших q провалы интенсивности флукту ирующего сигнала (замирания) резко замедляют рост раа как функции q.
§ 3.4. ОБНАРУЖЕНИЕ ПАКЕТОВ ИМПУЛЬСОВ
Пакетом или последовательностью № импульсов называют сигнал, образованный повторением с одинаковым интер валом (периодом повторения Та) N копий стандартного импульса s0(t; Эо). При этом копии отличаются друг от друга временем запаздывания и, быть может, значениями случайного векторного параметра Эо. Таким образом, общую модель пакета можно записать как
|
~N-l |
|
,s(/;9)=Re |
<iisl)(/-iT,;8i) , |
(3.25) |
|
i = 0 |
|
где ai— комплексные амплитуды, описывающие известный закон модуляции амплитуд и фаз импульсов внутри пакета (от импульса к импульсу), 9, — вектор неизвестных пара метров 1-го импульса пакета.
Когерентный пакет импульсон со случайной начальной фазой. Для пакета этого вида все N копий стандартного радиоимпульса имеют одну и ту же случайную начальную фазу и не содержат других случайных
параметров. |
Поэтому |
|
|
||
s(r; |
cn-i |
• |
) |
(3.26) |
|
(p) = Re< |
£ а(50(г-/Тп)exp[/'(2л/ог+ < )]>, |
||||
|
( |
i = 0 |
|
J. |
|
где |
<р — общая |
для всех |
импульсов случайная |
начальная |
фаза, подчиняющаяся равномерному распределению; S0(t) — известная комплексная огибающая одиночного импульса.
Очевидно, что при этом имеет место случай, рассмот ренный в § 3.2, так как сигнал (3.26) есть некая конкретная
•я-1 .
модификация (3.9) при S(f) = |
Однако струк- |
|
i = 0 |
тура обнаружителя может оказаться более удобной в ре ализации, если учесть специфику сигнала (3.26). Для
простоты ограничимся случаем, когда |
= — действитель |
ны *(arg« = O, л, 1 = 0, 1, ..., N— 1, т. е. |
at могут принимать |
50
как положительные, так и отрицательные значения). Под
ставив |
выражение (3.26) в (3.15), получим |
N-l |
I |
*!= Е |
W г2= £ a,z2(, |
i = O |
i = О |
где величины
Г1 т
Zi. = Re Y(t)S0(t-iTa)dt
о
т
= f y(J)S0(t - T* n)cos (2л/ог+yo(0)dz;
0
T
z2f = f y(‘)so(‘ -'TJsin (2л/01 + Yo(0)dz
0
(S0(l) и y0(t) — законы амплитудной и угловой модуляции стандартного радиоимпульса] могут быть сформированы как отсчеты на выходе фильтра, согласованного с одиноч ным импульсом (СФОИ), взятые в два момента времени, отстоящие друг от друга на четверть периода T0=l/f0 несущей радиоимпульса. Теперь структурную схему обна
ружителя |
пакета можно |
построить в |
соответствии |
с рис. 3.10, |
на котором |
ти — длительность |
одиночного |
импульса; Ен — накапливающий сумматор N выборочных значений (zu и z2i), опрашиваемый в момент окончания наблюдений (конца пакета) T=(N— 1)ТП+тя.
Показатели обнаружителя определяются соотношения ми (3.18) с учетом того, что в выражении для параметра
обнаружения |
q = s/2E/N0 должна фигурировать энергия |
||
|
JV- 1 |
Если амплитуды импульсов одина- |
|
всего пакета |
Е= Е |
||
|
i = 0 |
q = -J~Nq0, |
где q0— отношение сиг- |
ковы, то E=NE0 и |
|||
нал/шум на |
выходе |
СФОИ; |
Ео— энергия одиночного |
импульса. При проектировании обнаружителей часто необ
e(t-iTn-rll-Tll/4)
Рис. 3.10
51
ходимо знать минимальное число импульсов NM„„, обеспе чивающее при заданном q0 требуемые рл1, рпс. Для последовательности импульсов с одинаковыми амплитуда ми (|aj = const) с учетом (3.19)
v |
‘(^Z —21прдт, рпс)]2 |
^’мнн |
2 |
|
Qo |
Как отмечалось ранее, в случае малых рлт и рпс потери, связанные со случайным характером фазы, практически отсутствуют (^дВ~0) и потому вместо (3.19) можно вос
пользоваться |
выражением (3.8), так что |
л/ |
2 |
'’мин |
|
|
9о |
В последнем выражении равенство будет строгим для детерминированного пакета, в котором <р — известная вели чина и может без потери общности считаться равной
нулю. При этом |
схема |
рис. 3.10 |
упрестится— в ней |
|||
останется |
лишь один квадратурный |
канал |
и сравнению |
|||
с |
порогом |
будет |
подвергаться величина, |
накопившаяся |
||
в |
сумматоре. |
пакет. |
|
|
|
|
|
Некогерентный |
Рассмотрим пакет импульсов, |
у которого начальные фазы всех радиоимпульсов случайны и независимы друг от друга. Такой пакет называют
пекогерентным, его модель записывают в виде
•ф; |
|
Фо. |
Ф1> |
<Pn-i) = |
|
|
|
|
||
= Rej |
£ я,Л(г-г’^п)ехр[у(2л/ог + фЗ]1, |
|
|
|
(3.27) |
|||||
(1 = 0 |
|
|
|
|
J |
|
|
|
|
|
где |
и |
— действительные |
амплитуды |
импульсов |
пакета; |
|||||
S0(t) |
Тп имеют |
тот |
же |
смысл, |
что |
и |
в |
(3.25); |
||
ф;— случайные, |
независимые |
начальные фазы, |
подчиня |
|||||||
ющиеся |
равномерному |
распределению - |
И/0(<р;)= 1/(2 л), |
|ф(| < л.
Подставив соотношение (3.27) в (2.20) и повторив
выкладки, |
приведшие к |
(3.13), получим равенство |
||
( |
*EX- |
1 1 |
/ ехр |
2a-, Zj cos (<р; — arg z() |
/=ехр-----I I |
— |
|
dtp,, |
где величины Z, и argzf могут быть получены как отсчеты модуля и аргумента комплексной огибающей на выходе согласованного с одиночным сигналом фильтра в моменты времени /Т„ + т„ (т„ — длительность одиночного импульса),
52
е(1-1Тп-ти) |
e(t-T) |
in |
Рис. 3.11
N- 1
a E= £ E, — энергия пакета. После вычисления интегралов
i = O
и перехода к достаточной статистике In7 придем к правилу
где, как и в предыдущих аналогичных соотношениях, порог зависит от выбранного критерия и при исполь зовании критерия Неймана — Пирсона определяется задан ной вероятностью ложной тревоги рлт. В наиболее типич ном для практики случае прямоугольного пакета, в ко тором амплитуды а( одинаковы, структура оптимального обнаружителя имеет вид, показанный на рис. 3.11.
Характерной ее особенностью является критичность к виду амплитудной характеристики детектора Д, связанная с тем, что статистика, которая получилась бы на выходе сумматора при отклонении характеристики детектора от вида 1п70(-), не была бы взаимно однозначно связанной с £. Таким образом, в схеме рис. 3.11 оптимальный тип амплитудного детектора определен однозначно — это детек тор с характеристикой ln70(2Z//V0). Однако при слабых или сильных сигналах (<?0 1 или q0 » 1) возможны упрощения, основанные на приближениях функции ln70(x)s:x2/4 при х«1 и 1п70(х)«х при х»1. При отсутствии сигнала
(гипотеза Но), согласно (3.4), (3.16), Z2 = 2о2 — N0E. Поэтому
(2Z,/AO)2 =2qo и при q0 1 и q0 » 1 аргумент выражения 1п70(2ZJN^ с большой вероятностью будет соответственно мал или велик по сравнению с единицей. Аналогичные выводы нетрудно сделать и для случая истинности гипотезы Ну. Это позволяет, опираясь на приведенные ранее приближения для lnZ0( •), считать при q0 <к 1 оптимальным квадратичный детектор, а при q0 » 1 —линейный.
Расчет качественных показателей некогерентного обна ружения в общем случае является трудоемкой задачей, поэтому ограничимся лишь указанными ранее ситуациями,
ГИГПО Zf |
1 |
ул |
.1 |
1 |
1/0 |
А |
г* |
(/о |
А. |
53
N- 1
При ?0«l с порогом сравнивается сумма £ = £ ид[,
i = O где uai — отсчеты на выходе квадратичного детектора.
В этих условиях для получения представляющих практичес кий интерес значений рлг и рпс необходимо накапливать большое число импульсов N. Тогда на основании централь ной предельной теоремы величина £ будет распределена по нормальному закону при обеих гипотезах. Можно показать, что при q0 <е 1 появление сигнала приводит лишь к измене
нию среднего значения |
не влияя на дисперсию. Обозна |
||
чив среднее значение отсчетов при истинности Но и |
|||
через йдш и йдс, а |
дисперсию через Do и учтя, что при |
||
Т„ > т„ отсчеты ид>- |
независимы, |
получим |
|
Я^|Я0) = —L=exp |
2NDa |
J’ |
|
,/2яЛи>0 |
|
||
W^H^—=ехр |
2NDa |
(3.28) |
|
|
|
|
С учетом (3.28)
Введя нормированный порог й = |
окончательно |
получим |
ND0 |
|
|
Рлт=1-Ф(Л), Рпс=1-Ф(У^9д0-/!)> |
(3.29) |
где величину 9до = (йдс-йдш)/ч/^о можно считать отноше нием сигнал/шум на выходе детектора. Сравним эти результаты с полученными для когерентного пакета. Из выражений (3.29) и (3.5) следует, что для получения одинаковой верности обнаружения когерентного и некоге рентного пакетов должно выполняться условие
У^й?^до = >/^^о>
где N,* NHt— число импульсов, которое необходимо обра ботать в когерентном и некогерентном случаях. Если
54
учесть, что для слабых сигналов <7дО~<7о/2, то проигрыш во времени T=NT при обнаружении некогерентного пакета составит 4/ql раз. Так, например, если до = 0,1, то на обнаружение некогерентного пакета придется по тратить в 400 раз больше времени, чем когерентного. Эго означает, что некогерентная обработка слабых сиг налов практически лишена смысла; квалифицированный разработчик в условиях, когда режим слабого сигнала неизбежен (космическая связь, локационные и навигацион ные системы со сложными сигналами и др.), должен обеспечить возможность когерентного приема пакета.
В случае q0»1 ситуация в корне меняется. При неограниченном росте q0 требуемые рЛ1 и рпс можно обеспечить, обрабатывая лишь один импульс. При этом потери за счет незнания его фазы, как выяснено в § 3.2, невелики. Таким образом, некогерентная обработка силь ных сигналов почти столь же эффективна, как и когерент
ная.
Флуктуирующие пакеты. При изучении задачи обнару
жения флуктуирующих пакетов обычно ограничиваются двумя крайними случаями: 1) независимых флуктуаций,
описываемых моделью (3.27), в которой амплитуды импульсов полагают независимыми случайными вели
чинами, подчиняющимися распределению |
iy0(a;); |
2) друж |
ных флуктуаций, для которых в модели |
(3.26) |
а, = Ла0,-, |
где aOj — детерминированные множители, |
определяющие |
форму (закон амплитудной модуляции импульсов) пакета;
А — случайная |
величина, |
подчиняющаяся |
распределению |
W0(A). |
алгоритмы обнаружения |
флуктуирующих |
|
Рассмотрим |
|||
пакетов, начав |
с дружно |
флуктуирующего. Такой пакет |
не представляет ничего нового по сравнению с общей моделью сигнала со случайными амплитудой и начальной фазой. То, что этот сигнал является пакетом, находит отражение лишь в структуре комплексной огибающей S(t). Поэтому оптимальный обнаружитель имеет такую же структуру, как и для нефлуктуирующего когерентного пакета (см. рис. 3.10). Качественные показатели р„Т, р„с при рэлеевских флуктуациях А определяются первым из соотношений (3.18) и формулой (3.24), в которых следует считать
|
т |
N-1 |
. г |
Е—Ео £ |
£'0=- I So(t)dt. |
i = о |
2 J |
|
о |
55
Для некогерентного пакета независимо флуктуирующих импульсов с учетом независимости амплитуд и фаз
импульсов |
ОП |
|
|
|
|
|
|
||
N- 1 |
со |
|
fli2E0\ |
|
2aiZ,cos(<pl-argZ1) |
|
|||
■ 1 |
|
|
|
||||||
'=П |
2 л |
|
““vT/expL |
|
|
|
|||
1 |
= 0 |
|
|
|
|
||||
|
|
о |
|
|
|
|
|
|
(3.30) |
х Wz0(ai)da/d<pi, |
|
|
|
|
|
||||
где |
смысл |
обозначений |
Z,, |
Ео, argz; |
тот же, |
что |
|||
и в предыдущем пункте, а априорная ПВ |
для рэлеевских |
||||||||
флуктуаций |
|
|
|
|
|
|
|
||
W (д) = Ря‘ехР(-я-?)’ |
я‘ ^0; |
|
|
|
|||||
|
|
10, |
|
< 0. |
|
|
|
|
|
После |
вычисления интегралов |
в (3.30) |
|
|
|||||
|
N- 1 |
|
Z,2 |
|
|
|
|
|
|
1=С П ехР |
|
|
|
|
|
|
|||
N02(\+q02/l) |
|
|
|
|
|||||
|
i = 0 |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
где |
с = const. |
Перейдя |
к |
достаточной |
статистике |
In /, |
|||
получим правило |
|
|
|
|
|
||||
N- 1 |
И] |
|
|
|
|
|
|
|
|
О Е |
|
|
|
|
|
|
|
|
|
- = ° |
п0 |
|
|
|
|
|
|
|
реализуемое схемой рис. 3.11, в которой детектор с харак теристикой 1п/0(-) заменен квадратичным. Подчеркнем, что оптимальность квадратичного детектора, имеющая место для нефлуктуирующего некогерентного пакета лишь при <sc 1, для пакета независимо флуктуирующих им пульсов сохраняется при любых q0. Отметим также, что обнаружитель некогерентного независимо флуктуирующего пакета часто называют энергетическим приемником, так как вычисляемая в нем случайная величина £ в отсутствие шума пропорциональна энергии пакета.
§3.5. ОБНАРУЖЕНИЕ СЛУЧАЙНЫХ СИГНАЛОВ
Врадиоастрономии, радиоразведке, пассивной локации, биоэлектронике информация об источнике излучения либо каком-то явлении нередко связана с наличием или отсутст вием в наблюдаемом колебании y(t) реализации некоторого полезного ожидаемого случайного процесса. При этом решают задачи обнаружения случайного сигнала, описыва емого па языке и-мерных ПВ или функционала ПВ. Пусть в наблюдаемом колебании y(z) помимо белого шума n(z)
56
может еюдержаться реализация некоррелированного с n(t) нормального процесса s(l) с нулевым средним и корреля ционной функцией Ks(t, Л-т). Тогда при истинности гипотезы Я, процесс у(/) = ^(/) + и(/) как сумма некоррелиро ванных нормальных процессов будет также нормальным с корреляционной функцией Ky(t, Г-ьт), равной сумме
корреляционных |
функций |
з,(1) |
и n(t): |
|
|
Ky(t, t + i} = Ks(t, |
/ + т) + 0,5Яо5(т). |
(3.31) |
|||
Функционал |
ПВ |
у(/) |
при |
гипотезе Нх |
определится |
формулой (1.8): |
|
т т |
|
|
|
И/(у(1)|Я1) = сехр |
-1 |
V v |
|
G)y(G)dGdz2 ,(3.32) |
|
|
_ |
|
|
_ |
|
|
|
О О |
|
|
|
где X'71(zi> ?г)— обратная корреляционная |
функция у(1), |
получающаяся решением интегрального уравнения (1.7),
имеющего с |
учетом (3.31) вид |
|
f K,(G, |
V |
(3.33) |
z2)d/ = 5(/2-/1). |
||
О |
|
|
Разделив |
выражение (3.32) на функционал |
ПВ y(t) |
при гипотезе |
Но iy(y(/)|^0) [см. равенство (1.9)], после |
логарифмирования ОП и суммирования получившихся ин тегралов придем к правилу обнаружения случайного сигнала
ТHi
(3.34)
*«о
вкотором — порог, зависящий от избранного критерия, а
*>(')= Ь'(9) |
0) do. |
(3.35) |
О |
определяет некоторое |
линейное |
Соотношение (3.35) |
преобразование y(z), осуществимое, например, с помощью линейного фильтра (в общем случае с переменными параметрами). Поэтому обнаружитель, реализующий пра вило (3.34), можно построить по схеме, показанной на рис. 3.12 и повторяющей структуру корреляционного при
Рис. 3.12
57
емника (см. рис. 3.1) с той лишь разницей, что опорный сигнал sy(t) теперь формируется не автономно, а из самого наблюдаемого колебания ;>(/), пропускаемого через линей ный фильтр Ф.
Если процесс s{t) |
стационарен |
t + т) = /С5(т)] |
|
и достаточно широкополосен, т. е. его |
время |
корреляции |
|
т, подчиняется условию |
тк с Т, то |
можно, |
не внеся |
существенной погрешности, считать Ks(t — zj равной нулю за пределами интегрирования в (3.33) при любых . Тогда пределы интегрирования в (3.33) можно заменить бесконечными, превратив интеграл в обычную свертку. После этого переход к преобразованиям Фурье, согласно
теодэеме о свертке, приведет к |
равенству [^(/) + (Лго/2)] х |
в котором |
— спектральная плотность |
мощности случайного сигнала s(z), так что преобразо вание дФурье обратной корреляционной функции Ку 1 (/) = = l/[^s(/) + (Ar0/2)]. То же рассуждение позволяет записать (3.35) в виде свертки:
sy(t)= 1 y(e)h(t-e)de,
— оо
где /i(z) = 5(z) —(No/2) К~ r(z)— импульсная характеристика фильтра на рис. 3.12, имеющего в данном случае посто янные параметры. Коэффициент передачи этого фильтра Л(/) = 1 - (Л'о / 2) к;1 (/) = ^(/) / [£(/)+ (^ / 2)].
Будем полагать, что спектр мощности сигнала s(t) допускает прямоугольную аппроксимацию
£(/) = PC/(2F,), |ЖГ,; О, |/| <FB,
где Рс, F„ — средняя мощность и ширина спектра сигнала. Тогда фильтр Ф на рис. 3.12 оказывается идеальным
фильтром |
нижних |
частот с |
полосой |
пропускания |
Гв |
||||
и |
усилением в полосе пропускания РС/(РС + N0FB). |
Послед |
|||||||
нюю константу можно |
учесть |
непосредственно в |
пороге |
||||||
£п, |
считая |
коэффициент передачи |
ФНЧ равным |
|
|
||||
|
|
|
|
|
|
|
|
<3 36) |
|
£ |
Если |
вернуться |
к |
выражению |
(3.34) |
и в корреляции |
|||
колебания y(t) с реакцией sy(t) фильтра (3.36) на y(z) |
|||||||||
положить y(z) = 5 (z) + [y(z) —s^z)], |
нетрудно заметить, |
что |
|||||||
sy(f) и y(t) — sy(t) есть |
функции |
с неперекрывающимися |
|
e(t-T) |
in |
|
Рис. 3.13 |
|
|
|
[ —FB, FB], второй — за ее |
пределами [y(z) - sy(/)— то, |
что |
|
не пропускается ФНЧ]. |
Поэтому sy(t) и |
y(t) — sy(t) |
ор |
тогональны и их корреляция равна нулю. Таким образом,
(3.37)
о
и структура рис. 3.12 преобразуется в энергетический приемник (рис. 3.13), в котором решающей статистикой £ является энергия принятой реализации, пропущенной через ФНЧ.
Выражения для вероятностей ошибок рЛ1, рас энергети ческого приемника легче всего получить, воспользовавшись теоремой Парсеваля для разложения sy(t) в ряд Котельни кова. Так как sy(t) имеет отличный от нуля спектр только
в пределах полосы |
[ — |
Гв], то интеграл в (3.37) можно |
|
заменить суммой отсчетов sy(t) на интервале |
[О, Т], взятых |
||
через 1/(2Гв)с: |
|
|
|
2ГвГ |
|
|
(3.38) |
С= t Sy2[il(2F,)], |
|
|
|
i=i |
|
|
|
где множитель 1/(2FB) |
перед суммой опущен, так как он |
||
может быть учтен |
в |
пороговом значении. |
Значения р„, |
рпс могут быть найдены через интегральные %2-распределе- ния с 2FBT степенями свободы. Однако при соблюдении неравенства* 2FBT»1 величина £ как сумма мно гих независимых слагаемых может считаться нормальной
и для |
вычисления рлу, рпс достаточно определить сред |
||
ние и |
дисперсии £ при гипотезах |
Но и |
Нетрудно |
видеть, |
что ^=2F„TD {sy(t)}, где |
D {^(г)}—дисперсия |
(средняя мощность) на выходе ФНЧ. При гипотезе Но
С=(2ГВГ)РШ, а |
при гипотезе Нх t, = (2FBT)(Pm + P<.), где |
||
* Замена интеграла суммой (3.38) не приводит к потерям |
|||
информации только |
при |
этом условии, так как на поведение |
|
J,Q на отрезке |
[О, |
Г] |
влияют отсчеты sy(t) вне его и при |
небольшом значении произведения 2FBT пренебрегать их вкладом нельзя.
59