Добавил:
Закончил бакалавриат по специальности 11.03.01 Радиотехника в МИЭТе. Могу помочь с выполнением курсовых и БДЗ по проектированию приемо-передающих устройств и проектированию печатных плат. Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Литература / Гришин Ю.П., Ипатов В.П., Казаринов Ю.М. Радиотехнические системы (1990)

.pdf
Скачиваний:
461
Добавлен:
10.09.2023
Размер:
16.52 Mб
Скачать

что они пересекают аналогичные кривые для сигнала фиксированной амплитуды, соответствующие тем же зна­ чениям рЛ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

'’мин

 

В последнем выражении равенство будет строгим для детерминированного пакета, в котором <р — известная вели­ чина и может без потери общности считаться равной

нулю. При этом

схема

рис. 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