
Литература / Гришин Ю.П., Ипатов В.П., Казаринов Ю.М. Радиотехнические системы (1990)
.pdfРис. |
3.4 |
Uu^ = saJJ) = Ks(O) = E. |
(3.2) |
Белый шум на выходе |
СФ «окрашивается», и его |
корреляционная функция А’вых(т) по форме совпадает с Л^(т):
Квых(т)=0,5ад(т). |
|
|
(3.3) |
||
|
Из (3.3) следует, что дисперсия (мощность) шума на |
||||
выходе СФ |
|
|
|
|
|
^(О) = РШВЫХ = ^Ж |
. |
|
(3.4) |
||
|
Диаграммы |
рис. 3.3 отвечают случаю, когда |
порог |
||
z„ превышен и будет принято решение |
в |
пользу |
|||
гипотезы Н1. |
|
|
|
|
|
|
Рассчитаем вероятности ошибок дЛ1, рпс в оптимальном |
||||
обнаружителе |
детерминированного |
сигнала, |
пользуясь |
||
тем, что |
|
|
|
|
|
р„г = Р(Н,|Я0) = Р(с 2п|Я0) = J ir(z|770)dz; |
|
|
|||
|
|
zn |
|
|
|
рпс = Р(Я0|Я1) = Р(2<2п|Я1)= 7 W(z\HJdz, |
|
|
|||
|
|
— со |
|
|
( = 0,1. |
где |
lF(z|//,.)— ПВ корреляции z при |
гипотезе |
Hh |
||
на |
Графическая иллюстрация этих соотношений приведена |
||||
рис. 3.4, где площади заштрихованных областей |
равны |
||||
рлт |
(косая штриховка) и рпс (прямая |
штриховка). Так как |
z есть линейное преобразование нормального случайного
процесса (умножение |
на фиксированную функцию s(t) |
и интегрирование), то |
H/(z\Hi), где г = 0,1,— одномерные |
нормальные ПВ. Остается |
найти лишь их параметры: |
|||
среднее z |
и дисперсию D{z}. |
При |
отсутствии сигнала |
|
т__ |
т__ |
|
__ |
|
z = J y(0.r(0dz = f n(t)s(f)dt = O, |
так |
как |
/1(0 = 0. Появление |
|
о |
о |
|
|
|
сигнала на входе приводит к тому, что [см. (3.2)]
z = J y(t)s(t)dt - f [г(0 + /1(0 ]s(0df = fs2(0dz = E. oo о
40
Из физических соображений ясно, что дисперсия z, совпа дающая с дисперсией помехи на выходе СФ, не зависит от присутствия на входе сигнала и с учетом (3.4) D{z}=N0Е/2. Таким образом,
Рлт = 1 - Ф(й); Рпс = ф(л - я), |
|
(3.5) |
где Ф(х) = (1/^2л)'j exp(-/2/2)dz — интеграл |
вероятности; |
|
- со |
пороговый уровень; |
|
h = z„(yjN0E/2)~1 — нормированный |
||
q = y/2E/N0 — параметр обнаружения, |
равный |
отношению |
сйгнал/шум на выходе фильтра, согласованного с обна руживаемым сигналом s(t). График функции Ф(л) приведен на рис. 3.5. С учетом того, что Ф( — х) = 1 - Ф(л), выражение для рПС можно представить в виде
РПс=1-ф(<7-/!). |
(3.6) |
С помощью соотношений (3.5) и (3.6) осуществляется расчет обнаружителя в соответствии с принятым критерием оптимальности. Так, при использовании критерия Нейма на— Пирсона требуется минимизировать рпс при фиксиро ванном значении рлт. При этом из уравнения рлт=1-Ф(/|) следует найти нормированный порог Л = Ф-1(1 —рлт), где Ф-1()— функция, обратная Ф(х) (т. е. решение уравнения Ф(х)=у относительно х), и подставить полученное значение h в формулу для рПС (или рп0 = 1 -рпс). Зависимости
РПО = 1 -Рос = ф(<7 - h) = ф [<? - ф - (1* -рЛ1) ] |
(3.7) |
от q при фиксированных значениях вероятности ложной тревоги называют характеристиками обнаружения. Опи раясь на свойства интеграла вероятности Ф(х), легко установить, что зависимость рп0 от q является монотонно возрастающей, асимптотически стремящейся к единице при
41
<x>q*- . При q = 0 рпо = 1 -Ф[Ф-1(1 -рлт)]=рл1. Характерис тики обнаружения детерминированного сигнала приведены на рис. 3.6 (сплошные линии).
Часто бывает необходимо рассчитать минимальное значение параметра q, при котором достигается требуемая верность обнаружения, т. е. заданные значения рлт и рас. Это минимальное значение q = qMm определяет при задан ной спектральной плотности мощности шума NJ2 энергию сигнала Емяя, называемого иногда пороговым. Пользуясь
соотношениями (3.5) — (3.7), легко установить, |
что |
9ИИН = Ф’1(1-Рлт) + Ф’1(1-РПе)- |
(3.8) |
§ 3.2. ОБНАРУЖЕНИЕ СИГНАЛА |
|
СО СЛУЧАЙНОЙ НАЧАЛЬНОЙ ФАЗОЙ |
|
Рассмотрим сигнал, который не может считаться детерми нированным, так как содержит случайный параметр—фазу <р: ^(г;Э)=5'(г;<р). В общем виде модель такого сигнала можно записать как
s(r, <р) = S(r)cos [2nfot+у(0 + ф]=Re{5(0exp [/(27t/0t + ф) ]},(3.9)
где S(t) и y(t)—известные законы амплитудной и угловой модуляции; /0—-известная центральная частота; ф—слу
чайная начальная фаза с априорной ПВ И^/ф); S(t) = =5(г)еЛ(,)—комплексная огибающая сигнала s(f), являю
щегося |
реализацией |
s(t; ф) при ф=0: s(t)=s(t; 0). |
В |
соответствии |
с § 2.4 оптимальный обнаружитель |
должен формировать усредненное ОП (2.20) и сравнивать
42
его с порогом. Поскольку начальная фаза <р радиоимпульса
является неэнергетическим |
параметром, |
т. е. Е(<р) = Е = |
т |
(2.20) примет |
вид |
= (1/2) f S2(t)dt, выражение |
||
о |
|
|
п2z(<p) — Е
1= f |
exp |
lV0(<p)d<p, |
(З-Ю) |
|
— я |
т |
No |
|
|
где |
Пользуясь тем, |
что для любых |
||
z(cp)=f |
||||
|
о |
|
|
со со
функций u(t), v(t) f u(t)v(t)dt = f u±(z)y±(/)dz (равенство
Парсеваля для преобразования Гильберта), выражение для
z(<p) |
можно |
представить в виде |
|
|
||
z(<p) = Re |
1 г . . |
<P)dz |
|
|
||
-(/;*fy(/> |
|
|
||||
где |
|
. о |
s(t; ф)— аналитические |
|
отвечающие |
|
y(t) |
и |
сигналы, |
||||
y(j) |
и s(t; |
ф) (см. § |
1.3); * — знак комплексного сопряжения. |
|||
|
Так |
как |
y(t)=Y(t)exp(j2nfot), |
s(r, ф) = 5(/)ехр [j(2nfot + |
||
+ ср], где |
Y(t)— комплексная огибающая входной реализа |
|||||
ции |
y(f), |
то |
|
|
|
|
z(<p) = Re [zexp (—у'<р) ] = Z cos (ф - arg z). |
(3.11) |
|||||
|
В последнем |
равенстве |
|
|
z = H Y(f)S'(t)dt; Z=|z|.
(3.12)
2 о
Во многих задачах начальную фазу сигнала ф можно считать равномерно распределенной на интервале [ — л, л]: 1У0(<р) = ]/(2л), |ф|^л. При этом интеграл (3.10) с учетом
(3.11) имеет |
вид |
|
|
£ \ |
1 » |
2Z |
|
( |
v f exp |
No |
бф. |
/V Q) |
—соз(ф —argz) |
||
|
|
|
Воспользовавшись интегральным представлением модифи
цированной |
функции Бесселя нулевого |
порядка |
|||
|
1 |
" |
ехр[хсо5(ф-0)]бф, |
|
|
/о(х) = ^- f |
|
||||
|
2л |
|
|
|
|
окончательно |
получим |
|
|||
, |
/ |
E\Tf2Z' |
(3.13) |
||
/=еХр1 _ |
|
|
И |
||
|
\ |
/v0/ |
Vvo, |
|
43
|
|
Рис. |
3.7 |
Так как /0(х) ПРИ х^О монотонно зависит от своего |
|||
аргумента, |
то |
соотношение |
(3.13) позволяет решающее |
Н, |
1П |
записать как |
|
правило |
|
н„
Н.N
Z%Z„=^'[exp(E/N0)l„], |
(3.14) |
н„ 2 |
|
где |
1о 1(•)“Функция, |
обратная /0(х). |
|
||
|
Перепишем выражение для Z следующим образом: |
||||
Z=yjz2 + zl, |
|
|
|
||
где |
|
|
|
|
|
zi |
Re |
}*(Z)5*(/)d/ |
T |
cos |
|
|
(3.15) |
||||
|
Im |
= b'(W> . [2n/o(0 + Y(0]dr. |
|||
Z2 |
|
o |
Sln |
|
|
|
Таким образом, оптимальный обнаружитель сигнала |
||||
со |
случайной начальной фазой |
должен вычислять |
длину |
Z вектора с декартовыми составляющими Zj и z2. Как следует из (3.12), Z является абсолютным значением корреляции г комплексных огибающих принятого коле
бания У(г) и сигнала 5(z). При этом, согласно (3.15), zt и z2 есть корреляции принятой реализации у(0 с квадратур ными составляющими сигнала ^(z) = S(z)cos[2n/oz + y(z)] и j-±(/) = S(/)sin[2n/oz + Y(/)] — детерминированными колеба ниями, несущие которых сдвинуты по фазе на угол л/2
[s±(0 — преобразование Гильберта |
сигнала s(t) — s(c, |
0)]. |
Структура такого обнаружителя показана на рис. 3.7. |
па |
|
Отличие обнаружителя рис. 3.7 |
от приведенного |
рис. 3.1 состоит в наличии второго коррелятора и принятии решения по статистике Z, объединяющей выходные эффек ты обоих каналов. Если бы сигнал со случайной фазой
s(t; <р) обнаруживался |
как детерминированный, то |
при |
<р ~ + я/2 схема рис. 3.1 |
«нс замечала» бы s(t; ср) |
из-за |
44
слабой корреляции последнего с опорным сигналом кор релятора s(t). Благодаря тому что на рис. 3.7 опорные сигналы корреляторов находятся в квадратуре, статистика Z не зависит от <р, в результате чего устраняется вредное влияние случайности начальной фазы. Таким образом, инвариантность Z к начальной фазе сигнала s(t; <р) объяс няется тем, что значение <р влияет только на аргумент
корреляции (3.12) комплексных огибающих Y(t) и S(t), тогда как Z есть модуль z.
Иная реализация оптимального обнаружителя возмож на при использовании фильтра, у которого комплексная
огибающая импульсной характеристики H(t) = S*(T — t). Подобный фильтр согласован с сигналом s(t; <р), имеющим некоторое фиксированное значение <р, например <р = 0 [в этом случае фильтр согласован с первой из квадратурных составляющих сигнала, т. е. с $(0]. Огибающую на выходе
этого СФ |
Гвых(0 при воздействии y(t) на входе можно |
|
найти с помощью комплексного интеграла Дюамеля: |
||
Г.Ы.(0= I |
f Y(Q)H(t-Q)dG |
| f *(0y(e)s-(z-r»d0 |
|
- 00 |
- со |
При равенстве нулю сигнала за пределами интервала
наблюдения [О, |
Г], как следует из |
(3.12), |
yobIX(T) = Z. |
Таким образом, |
статистика Z может |
быть |
интерпрети |
рована как значение огибающей на выходе СФ в момент времени t = Т. Структурная схема обнаружителя на основе СФ приведена на рис. 3.8. Заметим, что вместо линейного детектора (ЛД) можно использовать любой, лишь бы его амплитудная характеристика была монотонной функцией огибающей входного процесса.
Для того чтобы рассчитать рлт, рлс в рассматриваемом случае, достаточно вспомнить, что отсчеты огибающей узкополосного нормального шума с дисперсией о2 распре
делены по закону Рэлея |
|
ip(Z|//0)=Rz^2<)QXp С—z2/C2cf2)L |
(3.16) |
Рис. 3.8
45
и подчиняются обобщенному |
закону Рэлея |
|
|
||||||||
|
( (Z/o2)exp |
/ |
- |
Z2+(/2 |
|
|
|
|
|||
|
|
|
+ |
м |
|
|
|
|
|||
IK(Z|H1) = ^ |
к |
|
|
2<у2 |
/ |
°\ d2 |
/’ |
|
(3.17) |
||
если к шуму добавляется сигнал с амплитудой С/м. |
|||||||||||
Как |
отмечалось ранее |
[см. выражения |
(3.2) |
и (3.4)], |
|||||||
на выходе СФ g2 = 7V0E/2, |
Uu-E. Поэтому |
|
|||||||||
рлт={ lK(Z|//0)dZ = |
— |
|
|
|
|
|
|
|
|||
z„ |
J /Ур |
|
|
2 |
|
|
2 |
2 |
|
|
|
|
z" |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
||
г, |
f |
|
Z |
|
f |
Z2 + E2\ |
f ZE\ , |
||||
A.c= f |
1----- exp------------ In I------- |
dZ. |
|||||||||
|
J |
^0E |
4 |
|
/ \N0EJ |
|
|||||
|
° |
|
2 |
|
|
|
2 |
|
2 |
|
|
Перейдя |
к нормированной |
переменной t = ZjyJN0Ej2, по |
|||||||||
лучим |
|
|
|
|
|
|
|
|
|
|
|
Рлт = ехР(-Л2/2), pnc = Q{h, q). |
|
|
|
|
(3.18) |
||||||
где h = ZJ,JN0EI2 — нормированный порог; |
q = yj2E!N0— |
||||||||||
параметр |
обнаружеш1я; |
Q(u, v) = |
|
|
|
o(«)dz — |
|||||
|
|
|
|
|
|
|
о |
|
|
|
|
табулированная (2-функция Маркума (интегральное распре деление Рэлея — Райса).
Для построения характеристик обнаружения необходи мо выразить нормированный порог h через заданную
вероятность |
ложной тревоги рят. |
Согласно (3.18), А = |
|
= |
21прлт. Подставив это в выражение для вероятности |
||
правильного |
обнаружения, придем |
к результату |
Рпо= 1 -Рпс= 1 -2(V~21прЛ1, ?)•
Характеристики обнаружения сигнала со случайной начальной фазой даны пунктиром на рис. 3.6.
Для определения порогового сигнала нужно |
решп п> |
||||
уравнение |
/>пс = Q(V ~ 21п/)Л1, </) |
относительно |
q: |
|
|
9м»я = Qi ’(х/-21прлт, рпс), |
|
|
(3.19) |
||
где £*(•,?£ |
■) — функция, |
обратная (2-функции |
по второмх |
||
аргументу. |
|
|
|
|
|
Cooi ношения (3.8) и |
(3.19) |
позволяют оцепить |
пог |
46
в пороговом сигнале, связанные со случайным характером фазы. Эти потери обычно характеризуют показателем
<?2мии(Рлт> Рис) 12 = Г g 2 W ~ 21Прдт, рпс) |
~12 |
|
41мни(Рлт>А,с)] |
|_Ф’1(1-/’лТ) + ф-1(1-Рис). |
|
где <71„нн(Рлт> Р„е) и |
92мим(Рлт> Рпс) —пороговые |
отношения |
сигнал/шум, необходимые для обнаружения с |
верностью |
рл1, р„с соответственно детерминированного сигнала и сиг нала со случайной начальной фазой. Величина £, показы
вает, |
во |
сколько раз |
следует увеличить энергию сигнала |
|||||
(т. е. |
его |
среднюю |
мощность рср |
при Т= const или |
||||
длительность |
Т при |
рср = const), |
|
|
||||
чтобы скомпенсировать |
снижение |
|
|
|||||
верности, |
обусловленное случай |
|
|
|||||
ностью начальной фазы. Обычно |
|
|
||||||
величину |
|
как и другие анало |
|
|
||||
гичные |
характеристики, |
выража |
|
|
||||
ют в децибелах: ^дБ= 101gЗави |
|
|
||||||
симость ^дБ=/(рпс) Для несколь |
|
|
||||||
ких значений р„ приведена на |
|
|
||||||
рис. 3.9 |
[5]. |
|
|
|
Рис. |
3.9 |
||
Как видно из приведенных кривых, значение потерь |
||||||||
зависит |
|
от |
заданных |
вероятностей |
ошибок, |
снижаясь |
с уменьшением значений рдт и рпс. Благодаря малым значениям Е, при малых рл1 и рпс ориентировочный расчет
порогового отношения сигнал/шум для модели |
сигнала |
со случайной фазой нередко проводят по более |
простой |
формуле (3.8), полученной для модели детерминированного сигнала.
§ 3.3. ОБНАРУЖЕНИЕ СИГНАЛА СО СЛУЧАЙНЫМИ АМПЛИТУДОЙ И НАЧАЛЬНОЙ ФАЗОЙ
Во многих практических задачах случайным пара метром обнаруживаемого сигнала наряду с началь ной фазой <р является и его амплитуда. В этом случае в модель сигнала (3.9) вводят дополнительный
случайный |
параметр — амплитуду А: |
||
л ЛЧ _ А С (Л |
ГЭ-г f 4 1 *,<4\ Xml- |
||
s V> Ф, zi; — |
v,|xit/ol -г [UM VJ — |
||
= Ие{Л5(/)ехр [/(2л/01 + <р)]}. |
|||
Когда амплитуда и начальная фаза независимы, сов |
|||
местная |
априорная |
ПВ Лиф fV0(A, ф)= Ил0(Л)Ил0(ф). |
|
Согласно |
(2.20), ОП |
|
47
00 |
я |
|
|
2z(A, <p)-E(A) |
(3.20) |
|
и/0(л)и/0(<р)с1ла<р, |
|
Q |
-ft |
|
где |
z(A, <р) = Ле(<р); E(A) = A2E; z(<p) и E имеют |
тот же |
смысл, что и в (3.10), т. е. представляют собой корреляцию принятого колебания y(f) с сигналом единичной амплитуды
s(t; <p) = Re{S(/)exp [/(2л/01 + <p) ]} в функции от фазы <р и эне ргию s(t; <р). Считая ПВ начальной фазы Й/0(ф) равно мерной и выполняя в (3.20) интегрирование по <р с учетом результатов § 3.2, получим
СО
0
где величина Z определена в (3.11), (3.14).
Так как подынтегральное выражение в (3.21) неотри
цательно, |
а /0(х)— монотонная функция своего |
аргумента |
|
при х>0, |
то |
/ также монотонная функция |
статистики |
Z и структура оптимального обнаружителя при любом |
|||
распределении |
1У0(Л) будет такой же, как и для сигнала |
со случайной начальной фазой и фиксированной амп литудой. Таким образом, структуры, приведенные на рис. 3.7, 3.8, сохраняют оптимальность и в том случае, когда у сигнала случайна не только фаза, но и амплитуда,
причем конкретный вид |
ПВ fV0(A) |
может влиять только |
|||
на |
пороговые значения |
Z„. |
|
|
|
|
Вероятность ложной тревоги определяется первым из |
||||
выражений (3.18). Для |
вычисления |
вероятности |
пропуска |
||
рпс |
необходимо задать вид распределения |
IVO(A), |
так как |
||
|
со |
|
|
|
|
Р„с = f pnc(A)W0(A)dA, |
|
|
|
(3.22) |
|
|
О |
|
|
|
|
где |
дпс(Л)— вероятность |
пропуска |
сигнала |
с конкретным |
значением А. Чаще других на практике встречается рэлеевс кая модель флуктуаций амплитуды
Л |
/ |
А2 \ |
|
( —ехР -т-у , Л>0; |
(3.23) |
||
И/о(^) = ) ° л |
\ |
2с V |
|
( 0, Л <0, |
|
|
хорошо согласующаяся, например, с реальными мерца ниями сигналов, отраженных от радиолокационных целей, федингом в линиях радиосвязи и т. и. Выбор параметра
48
ст л здесь произволен, так как истинная огибающая сигнала
в момент t есть Л|5(г)| = AS(f) и увеличение А во сколько-то раз при соответствующем уменьшении S(f) значения оги бающей не изменит. Обычно ст^ выбирают из условия
Л2=1. При этом средняя энергия сигнала Е(А) = А2Е равна
энергии сигнала Е при единичной амплитуде. |
Так |
как |
для рэлеевской случайной величины А2 = 2а2л, |
то |
для |
выполнения условия Л2=1 необходимо, чтобы СТл=1/2. Для рэлеевских флуктуаций А значение рпс можно вычислить, не прибегая к интегрированию в (3.22). Дейст вительно, мгновенные значения радиосигнала с начальной фазой, распределенной равномерно, и с амплитудой, флук туирующей по закону Рэлея, нормальны с нулевым
средним. Таким образом, при истинности |
гипотезы Hi |
на выходе СФ (см. рис. 3.8) наблюдается |
сумма незави |
симых нормальных процессов, один из которых (шум)
имеет дисперсию N0E/2 |
[см. (3.4)], |
а другой (сигнал) — |
с2лЕ2 = Е2/!. Поскольку |
эта сумма |
нормальна и имеет |
нулевое среднее, ПВ при гипотезе |
Нг описывается не |
(3.17), а (3.16), т. е. является рэлеевской (вместо обобщен ной рэлеевской), причем фигурирующий в (3.16) параметр
ст2 теперь |
равен N0Ej2 + Е2/2. Поэтому |
||
Zn |
Z |
|
|
f |
|
||
J^E/2 + E2/2 |
2(ад2 + Е2/2) |
||
При t-ZlJN0Е(1 + E/N0)/2 получим |
|||
Pnc=1~exP |
h2 |
(3.24) |
|
2(1 W/2) |
|||
|
|
где h имеет тот же смысл, что и в (3.18), а параметр обнаружения q соответствует сигналу с единичной ампли тудой, т. е. усредненному по А значению мощности сигнала при флуктуациях, описываемых (3.23). Для построения характеристик обнаружения служит равенство, следующее
из |
(3.24) |
и |
формулы для рлт в (3.18): |
п |
=1 —п |
=р ,+9/2 |
|
Г |
ПО 1 |
Г ПС |
гит |
Характеристики обнаружения сигнала со случайной фазой и рэлеевскими флуктуациями амплитуды на рис. 3.6 нанесены штрихпунктиром. Их особенность состоит в том,
49