Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Ширман, Я. Д. Разрешение и сжатие сигналов

.pdf
Скачиваний:
16
Добавлен:
22.10.2023
Размер:
13.3 Mб
Скачать

например, с помощью линий задержки. Для этого требуется обеспе­ чить изменение задержки на ± (//с) cos Ѳыако в диапазоне углов I 0 I < 0Макс- Функции линий задержки СВЧ могут выполнять моно­ кристаллы ЖИГ [82], элементы линзовых антенн и др. Эти элементы должны обеспечить возможность фокусировки антенны во всем рабочем диапазоне частот, не внося в оптической терминологии «хроматиче­ ской аберрации».

Рис. 2.1.13. Тело автокорреляционной функции дальность — угол широкопо­ лосного пространственно-временного сигнала (т = 0,05).

Широкополосность сигнала повышает в отдельных случаях также возможности и чисто у г л о в о г о р а з р е ш е н и я . Рассмотрим систему из двух ненаправленных приемных антенн при расстоянии / между ними. Под углом 0 к оси системы падает волна негармонических

колебаний (1). В антеннах наводятся э. д. с.

(t + I cos Ѳ/2с) с ком­

плексными амплитудами

 

—Я/2

Если система согласована для направления 0', то модуль коэффициента корреляции приходящих и ожидаемых колебаний по аналогии с (3) будет

со

 

 

 

J

[tfi (t,

0) Ui (t, ѳ') + ия (t,

0) U$ (t, 0')] dt

---- CO

 

 

( 10)

p ( M ' ) =

CO

 

 

J

0)l2+ l u2(t,

0) \4dt

 

— CO

 

 

§ 2. 1.9.

219

Подставляя (9) в (10) н замечая, что

 

$

еі2ли''~!'"и dt = 8 ( F ~ /■'),

 

—со

 

 

получаем

 

 

р(0,

0') = ро (0, 0') Рі (0, 0').

(П)

Здесь ро (Ѳ, 0') — модуль коэффициента корреляции при узкополос­ ном сигнале

J -*

1,0 ,

0,8

1-----

 

1----------

О )

 

0,6

----

 

 

 

 

~ -^ х

9 2

р0 (0, 0') = cos

л/

X

 

к,.

 

X (cos 0 —cos O') .

( 12)

a Pi (0, 0') - сомножитель,

связанный с широкополосностыо,

Рі(Ѳ, Ѳ') sin

; x

j

X(cos0 — cos O')

Я/, X

X (cos0 —cosO')/. (13)

 

Н

I 0.2

\

 

 

В

соотношении (13) ly —

 

 

 

 

 

 

— in/f0 — «эквивалентный

 

 

 

 

 

 

размер»

каждого

из облу­

 

 

1

 

 

 

чателей,

обусловленный

- в

- U

0

и

в

у

■ >

широкополосностыо сигна­

Рис. 2.1.14. Вертикальные

сечения

взанмо-

ла.

Расширение

полосы

корреляционной

функции

дальность— угол

сигнала приводит к уве­

широкополосного пространственно-временного

сигнала

при настройке антенны

на

сигнал с

личению этого «размера».

 

«нулевой» полосой.

 

 

В силу установленного ра­

 

 

 

 

 

 

нее

соответствия

между

коэффициентом корреляции и характеристикой направленности мож­ но сказать, что расширение полосы приводит к увеличению эквива­ лентной направленности каждого элемента системы. Оно как бы «заполняет» имеющийся разрыв между точками приема [41].

С увеличением разноса между элементами и с приближением к ним излучателей становится необходимым иногда учитывать кри­ визну фронта волны.

§ 2.1.10. РАЗРЕШЕНИЕ ПО НАПРАВЛЕНИЮ ПРИХОДА И КРИВИЗНЕ ФРОНТА ВОЛНЫ

Когда излучатели приближаются к элементам приемной антенной системы, появляется возможность разрешать не только по направлению прихода, но и по кривизне фронта волны, зависящей от расстояния до излучателя. Поясним это на примере линейной распределенной антенны.

220

§ 2.1.10.

Пусть на линейный раскрыв (см. § 2.1.4) падает волна от точечного источника И, находящегося на конечном расстоянии г от центра рас­ крыва (рис. 2.1.15), направление на источник составляет с осью рас­ крыва угол Ö. Расстояние/? = R (х, 0, г) до произвольной точки рас­ крыва будет

R (х, 0, /') = \А\л -(-х2-—2/-.VCOS0 .

(1)

Если источник И создает колебания ф0 (/), то в пренебрежении ам­

плитудными

различиями принимаются колебания

4’ {і, х) = Фо U — R {х,

0,

r)/cl.

 

Анализ

распределения

колеба­

 

ний

по

раскрыву

проведем

 

для

гармонического колебания

 

фо (t) =

e i 2nf°‘,

полагая,

что

 

кривизна

фронта еще

невелика

 

и в

разложении

R (х,

0,

г) по

 

степеням

отношения хіг можно

 

пренебречь

кубичными

члена­

 

ми, т. е.

 

 

 

 

 

 

R

(х,

0,

г) «

г X cos 0 +

Рис. 2.1.15. Пояснение геометрических

 

 

+

(л-2/2г) sin2 0.

 

(2)

 

 

 

соотношений.

Вводя комплексные амплитуды разрешаемых колебаний для двух различных положений излучателя 0', /•' и 0, г, находим нормирован­ ную автокорреляционную функцию

 

 

1/2

£ /2я [ « (.V, 0, r ) - R (л-. О', ОІЛ 0 dx

 

р(Ѳ. 0', г, /•')=

^

(3)

 

 

-1/2

 

 

Замечая,

что

 

 

 

 

j e

dt = C(u) + jS(u),

(4)

 

о

 

 

 

где функции С и S — косинус- и синус-интегралы Френеля, получаем

 

 

р(Ѳ, 0', г, г') =

 

 

= Ѵ Iе (“ і) — С(«а)]а -і-[S (иг)— S {и2)]2/\и1~ и 2\.

(5)

Здесь:

и1>2 = В/]/Л ±

Ѵ А ,

А = (IV2К0) (sin2 Ѳ/г — sin2 O'/r'),

В = (ZA0) (cos 0 — cos 0').

 

 

 

§ 2. 1. 10.

221

Когда кривизна фронта невелика, а различия в ней несуществен­ ны, пользуясь асимптотическим равенством (гг > 1)

 

2

)пи

 

приходим к выражению [(19), § 2.1.4] для

предельного случая плоской

волны.

 

(0 = 0' п В =

 

Когда отсутствуют угловые различия

0), то р ---

= С { у А ) Г\ 'ГА,

где при | Дг | = | г г' | <£ г величина

А ж (г —

— г ) В sin2 0/2.’Л,0

= уЛ/\

 

 

Кривая р = р (|уДг|) (рис. 2.1.16) характеризует разрешающую способность по дальности для согласованной р а с п р е д е л е н н о й линейной антенно-приемной системы. Если аналогичную разрешающую способность ввести и для дискретной системы, то расширение полосы

•частот сигнала приведет, как и в § 2.1.8, к приближению ее возмож­ ностей к возможностям распределенной системы.

§ 2.1.11. ОПТИМИЗАЦИЯ РАЗРЕШЕНИЯ ПРИ НАЛИЧИИ СЛУЧАЙНОГО НЕЭНЕРГЕТИЧЕСКОГО ПАРАМЕТРА МЕШАЮЩЕГО СИГНАЛА

Пусть наряду со случайной начальной фазой ф и релеевским ампли­

тудным множителем Ъ (Ь2 = 2) мешающий сигнал зависит от случай­ ного неэнергетического параметра а. Последний может описывать вре­ менное положение, длительность, частоту или производную частоты ко­ лебаний сигнала, направление прихода соответствующей ему волны,

ееполяризацию или кривизну фронта.

Доопытное распределение параметра сс считаем нормальным

р (а) = (2яа§)->/2 ехр [—(а — а 0)2/2о-2].

О )

222

§ 2.1.11.

Дисперсия al этого распределения характеризует априорную неопре­ деленность параметра, а величина а0— его математическое ожидание.' Послеопытное распределение а условимся также считать нормаль­

ным. Действительно послеопытное распределение параметра интен­

сивного п о л е з н о г о сигнала в шуме оказывается нормальным [15,

43, 70]. Обычно нас интересует случай и н т е н с и в н о г о

м е ш а ra­

ni, е г о сигнала.

В отсутствие полезного сигнала (а например, при

радиолокационном

обнаружении это наиболее вероятный

случай)

к параметру мешающего сигнала можно применить ту же закономер­ ность. В присутствии полезного сигнала нормальный закон послеопытного распределения можно с успехом применить для наиболее интересного варианта обнаружения, когда полезный сигнал значитель­ но слабее мешающего. Наконец, нормальное распределение, непремен­

но, имеет место, когда

достаточно велика априорная информация

о параметре и значение

мало: в этом случае послеопытное распре­

деление определяется доопытным. Из всего сказанного вытекает, что случай нормальной аппроксимации послеопытңого распределения па­ раметра а представляет достаточный интерес. Обсуждение приемле­ мости примятой гипотезы будем продолжать по ходу изложения.

Итак, пусть

принимаются

колебания

с комплексной амплитудой

 

U (і) = N(t) +

bUj_ (t , a) e ^ + aU2 (t),

(2)

где N (t), U1 (t,

a),

ff2 (t) — комплексные

амплитуды шума, мешаю­

щего и полезного

сигналов. Множитель

а, принимающий

значения

О или 1, подлежит статистическому оцениванию. Чтобы упростить анализ, случайный характер амплитуды и начальной фазы учиты­ вается лишь применительно к мешающему сигналу. Наличие этой случайности нельзя не учесть без искажения сущности задачи. Если бы все параметры мешающего сигнала были известны, его удалось бы сразу вычесть из принимаемого колебания.

Для выработки решения относительно а достаточно составить ло­

гарифм отношения правдоподобия

 

 

Л [U{t)\ = ln {Pl [ U (t) - U 2(t)]lPl [ff(*)]}.

(3)

Здесь

[ff (/)] — плотность вероятности шума и мешающего

сигна­

ла (имеется в виду обычный предельный переход от дискретной реа­ лизации к непрерывной, § 1.1.2)

со

Pi [ff(01 —

$ p1[U(t)\a]p(a)da,

(4)

 

—00

 

Py[U [t)\a] — условная плотность вероятности U {t) при

фиксиро­

ванном значении параметра

а.

 

Анализ (3) упрощается, если ввести логарифм отношения правдо­ подобия для трех вспомогательных случаев обнаружения:

1) полезного сигнала в отсутствие мешающего

 

Л0 [ff(01 = ln {Po [U{t) — U2’{t)\/p0[ff(Ol}-,

(5)

§ 2. 1. 11.

223

 

2)

мешающего сигнала

в

отсутствие полезного

(а = 0), когда

параметр а случаен,

 

 

 

 

 

 

 

 

 

 

л л а д -

in {Pl[U(t)\/p0

 

(6)

 

3)

мешающего сигнала в отсутствие полезного =

0), когда пара­

метр

а фиксирован,

 

 

 

 

 

 

 

 

 

 

Л, [U{t) I а] = In {рх [17(0 I a]lp0I*7(01}-

 

(7)

Здесь ро [U (г!)1 — условная плотность вероятности U (/) при а =

b =

=

0 (условная плотность вероятности одного шума).

 

 

из

В силу известных свойств логарифмов произведения п частного,

соотношений (3),

(5),

(6)

получим

 

 

 

 

 

А [£7(01 =

Лх [17(0 —

(01 —At [U{t)] + Л0 [U{t)].

(8)

 

В соответствии с (4),

(6), (7) и формулой полной вероятности

 

 

 

Лі [U(t)] = ln

^

exp (Л; [U(t) \а)}р(а) da.

(9)

 

 

 

 

—Л)

 

 

 

 

 

Введем следующие записи логарифмов отношения правдоподобия

для вспомогательных задач:

 

 

 

 

 

 

 

 

Л0 W (01 =

Re Z0 [U (01 -I- const,

 

(10)

 

 

Лх [U (0| a]

= Z\

[U (0 I a] + const.

 

(11)

Здесь

 

 

 

 

со

 

 

 

 

 

 

 

 

 

 

 

 

(12)

 

 

z 0[U(t)] =

- ±-

5

U(t)u:_(t)dt,

 

 

 

 

 

 

 

Л'о

 

 

 

 

 

 

Zi{U( 0 |a ] =

 

1

 

j

U(t) U[ (/, a) di

= 2, (a)

(13)

 

 

 

 

 

 

V2N0 (3,+No)

— нормированные корреляционные интегралы в задачах обнаружения сигналов с полностью известными параметрами и со случайной на­ чальной фазой, величина — средняя энергия мешающего сигнала.

Выражения (10)—(13) следуют из [(27)—(30), (32), § 1.1.31.

Из выражений (9), (11), (13) получим

со

 

Л; 1£7(01 = 1° 5 lexP Z?(a )] P(a) da-f-const.

(14)

Подынтегральное выражение (14) представляет собой с точностью до множителя послеопытную плотность вероятности параметра а (см. также § 3.1.3). В соответствии с принятым ранее условием исполь­ зуем для нее нормальную аппроксимацию

[exp Z\ (а)] р (а)

С

ехр

сім)2~

 

(15)

Ѵ2п стм

_ О

»

 

 

2ам

 

 

224

§ 2. 1. 11.

откуда, логарифмируя и подставляя (1), получаем

2? (а) — ( а —кц)2

(к а м)2 + 1 п С _

1п £ м -

(16)

2а?

2ам

ст0

 

Параметры аппроксимации С, аы, а м можно подобрать, сопоставляя значения левой и правой частей равенства (16), их первые и вторые производные в некоторой произвольной точке а = ах, в окрестности которой аппроксимация справедлива. Система уравнений, определяю­ щая параметры аппроксимации, имеет вид

In С = Z{ (ах) +

(а1—ССД,)2

(«!

ао)2+

Іп ам

(17)

 

2ом

2а02

во

 

« 1— «м _ « 1—«о___dZ\ (gt)

 

(18)

ам

а?

да

 

 

 

 

 

J _ _

J ___ д2 Z\ (a,)

 

 

(19)

Ом

во

ö a 2

 

 

 

 

 

Здесь уравнение (19) определяет величину

стм,

уравнение

(17) —

величину In С, а уравнение (18)— а м. Чем меньше значения сг0 и N 0, тем точнее определяется параметр а, меньше его послеопытная диспер­ сия о,?,, тем больше оснований для пренебрежения в (16) высшими членами разложения Zx (а) по степеням (а — а м)2, а значит, и для использования самой нормальной аппроксимации (15).

Используя эту аппроксимацию, выполним операции интегрирова­

ния

и возведения в степень,

предусмотренные

соотношением

(14).

В результате преобразований приведем (14) к виду

 

 

Ах [U (01 = ln С =

Z\

\U (t) I ax] + ДА [ü (t) | o j,

(20)

где

слагаемое

AA [U (t) | ссг]

назовем поправочным членом. В

соот­

ветствии с (14), (17)—(19) поправочный член будет

 

 

 

ДА [£/(/) К ] =

^

(а, — а 0)2д- Z\ (аг)

 

 

 

 

 

 

да2

 

 

 

 

— 2(ccj — а 0) dZ\ (ах)

dZ\(а 0

X

 

 

 

 

да

да

 

 

 

X

2 a* Z\ ( ах) -

-1

---- —1п ( 1

:d-Zl (Kl)

(21)

 

1— or;

 

 

<Эа2

 

 

öa2

 

 

 

 

величина его изменяется в зависимости от выбора ах. Поскольку левая часть равенств (20) от ссх не зависит, при этом меняется соотношение между слагаемыми в правой части (20).

Наряду с исследованным выражением А г [U (/)], в исходное (8)

входит выражение А х [U (і) U2 (/)]. Для

него справедливо соот­

ношение, аналогичное

(20),

 

Ах W ( t ) - U 2(t)} =

Z { { U { t) - U 2{i) |а о ]+

ДА [ U ( t ) - U 2 (0|а,1 .

 

 

( 22)

8 Зак. 1303

225

Значение а 2 может быть взято как равным, так и неравным а х. В за­ висимости от выбора значений а х и а 2 реализуются различные варианты оптимального или квазиоптпмального разрешения.

Пусть, например, значения а х и а 2

выбираются из различных

условий

 

ДАХ[U(t)\ а х] =

О,

ДЛа W (t) Uа ( 0 |а а1 = О,

т. е. полагаются в общем случае неодинаковыми. Тогда выражение

(8) сводится к алгебраической сумме

Л [U (01 « Z\ W

(t) - U2

(t) I ct2] -

—Z\ [£/ (0 I a x]

+ Re Z0

[U (t)].

Слагаемые этой суммы описывают выходные напряжения трех прием­ ников — двух квадратичных для оптимального приема мешающего сигнала (за вычетом и без вычета полезного) и одного линейного для оптимального приема полезного сигнала в отсутствие мешающего. В каждом из первых двух приемников производится своя о ц е н к а параметра а (с учетом доопытной информации, если она существенна), выходной эффект приемника выдается для соответствующего оценоч­ ного значения ссх или а 2. Любое из этих значений обращает в нуль выражение (21):

— первое из них а ь когда Z (ал ) определяется выражением (13) не­ посредственно;

— второе а„, когда U (t) в выражении (13) заменяется на U (t) —

-U2 (t).

Значения аъ а 2 можно выбрать также одинаковыми, приравнивая

нулю разность ДА [U {t) (J2 (t) | a] — ДА \U [t) | а]. Используя ко­ рень a 2 полученного уравнения

 

ДА IU (t) — U2 (t) I a i]

= ДА IU (0 (ail,

(23)

приходим

к единой обработке

 

 

 

 

A [U (01 = Z*[U (t) — U2 (t) I a x] -

 

 

—Z\ W (t) |aj]

+

R eZ 0 IU (t)].

 

С учетом

(12), (13) эта

обработка

преобразуется к виду

 

 

 

со

 

 

 

 

А [Щ*)1 =

— Re X U(t) R* (t, a j dt -f const.

(24)

 

 

No

 

 

 

Здесь R (t , a x) — решающая функция,

соответствующая обнаружению

когерентного сигнала при наличии одного мешающего сигнала со случайными амплитудой и начальной фазой для фиксированного зна­ чения его параметра а = аг.

Она определяется соотношением [(31), § 1.1.3] при і — 1, а в не­

сколько видоизмененной форме может быть представлена как

 

R V, а,) = Щ{і) -

РІ, а К) U, ах).

(25)

226

§ 2. 1. 11.

Здесь

ql = V 2 3 J N Ü и q%= }/23,/A 0

— отношения сигнал/шум

для

мешающего и полезного сигналов;

рх>2 (а) — нормированная

взаимокорреляционная функция

 

 

со

 

р 1 ’2 ^

= у г і ж

J

а) и '~^ du

(26)

где

 

 

 

 

I

1Ui iß, «) I2 dt,

3,

1 I U2{t) ]2

dt,

a выше было принято b\ = 2.

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

определения значения а ь

Так, при а0= сгм = 0 решением уравнения

(23) будет а*

= а 0, т. е.

в этом случае допустимо использовать д о-

о п ы т н у ю

о ц е н к у

п а р а м е т р а . При а0 > ам и интенсив­

ности мешающего сигнала,

значительно превосходящей интенсивность

полезного, в

качестве

можно принять

о ц е н к у

м а л ь и о г о

п р а в д о п о д о б и я

(см. также

определяется уравнением

dZj_ (а)

да

ам а к с и -

§3.1.4). Она

(27)

к которому сводится (23), когда интенсивность мешающего сигнала существенно превышает интенсивность полезного.

§ 2.1.12. КАЧЕСТВЕННЫЕ ПОКАЗАТЕЛИ КВАЗИОПТИМАЛЬНОГО ОБНАРУЖЕНИЯ ПОЛЕЗНОГО СИГНАЛА НА ФОНЕ МЕШАЮЩЕГО СО СЛУЧАЙНЫМ НЕЭНЕРГЕТИЧЕСКИМ ПАРАМЕТРОМ

Величина А = A [U it)], получаемая в зависимости от реализа­ ций шума и мешающего колебания после определения оценки согласно [(24)—(27), §2.1.11], является случайной гауссовой величиной. Дис­ персия и математические ожидания величины А в отсутствие и при наличии полезного сигнала характеризуют потенциальные возможности его обнаружения. Расчет проведем в предположении, что обработка [(24), §2.1.1] соответствует некоторому фиксированному оценочному

значению параметра

[136].

Математическое ожидание величины А (за вычетом константы) в отсутствие полезного сигнала равно нулю, а при наличии последнего в соответствии с [(24) §2.1.11] будет

 

_1_ Re

СО

 

 

 

 

Ао

 

U2(t)R*(t,

a1)dt = qly(q1, cti).

(1)

 

No

00

 

 

 

 

где

 

 

 

 

 

 

 

q2

 

 

 

 

 

 

 

 

(2)

 

ф(<7ъ ai) =

l — 9 *

I Pi.a («л.) |a.

 

 

 

<7? +

2

 

 

§ 2.1.12.

8*

227

Дисперсию величины Л при фиксированном аг как при наличии, так и п от­ сутствие мешающего сигнала представим в виде

Re j a j d t

1

J

U(t)R*(t, а у) dt = І

ШЯ V W

X* ('■ “ i) X (s, « i) X

Щ

 

2/Vо

—со

 

 

 

X p (а I a t) p (<|') p (b) dt ds da dip dö.

(3)

Здесь: U (t) — комплексная амплитуда входного напряжения [(2), §2.1.11] в от­

сутствие полезного сигнала (о 0); горизонтальная черта — знак

усреднения

по реализациям N (/), входящим в [(2), § 2.1.11 ]; р (Ь) =

b ë~b't2 (b »

0); р (і|>) =

= 1/2я, (I ф I <С л)і Р (alaі) — плотность вероятности

истинного значения па­

раметра мешающего сигнала UL{t, а) в [(2), §2.1.11]

для оценочного значения

параметра а г. Принимая нормальный закон распределения, полагаем

1

(tt—Ді)2

-сгЯ

p (a|c£ i)= y l s ,

(4)

 

Заметим, что в соотношении (3) мы пренебрегли зависимостью ошибки (а — о^) от значений Ь, ф, N (1), £А (О- Интегрируя (3) по t и используя [(25), (26), §2.1.11], получаем

 

со

 

 

 

 

 

 

ОО

 

 

 

 

 

С [N{t) + bUr { t , a ) ^ ] R * { t ,

a1)dt =

f

N(t) R* (t, a t) dt +

 

 

-----CO

 

 

 

 

 

 

-----CO

 

 

 

 

 

 

-ЬУгЗіЭгб е'ф Pi, 2 (a ) —

—У "

- p i , » ( a i) P i, i( a .

<*i)

 

(5)

 

 

 

 

 

Ql ~r -

 

 

 

 

 

 

где рЬ1 (а, ctj) — нормированная

автокорреляционная

функция

мешающего

сигнала

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P i . i ( a , a i ) = " T "

IUi (I ,

a)

U* (I, arf dt.

 

 

( 6 )

Аналогично проведем интегрирование по s. Замечая,

что N (і) =

0, N (t) N* (s) =

= 2N08 (t — s), выражение для

cr^ приведем к виду

 

 

 

 

 

ОО

 

 

 

 

ОО

 

 

 

 

 

 

 

Зл

J \R{t,

a 1)\2dl + ~ ^ - ^ [ ( a ) p ( a \ a 1)b,1p(b)dadb,

(7)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f( a) =

 

 

Ч

 

 

 

 

 

 

 

 

(8)

 

Р і,г (а )-

 

<?2 +2

Р і.2 (« і) Р і,і (<*>

« і)

 

 

Для гладкой функции / (а), описываемой

в окрестности | a — a x | <

3aMпервы­

ми тремя

членами разложения

в ряд

Тейлора, справедливо

асимптотическое

приближение

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

(а —(%і)2

da X

f (а{)

 

f" (ai).

(9)

 

1/2лГ стм

 

 

2а*

 

 

 

 

 

 

 

 

228

§ 2. 1. 12.

Соседние файлы в папке книги из ГПНТБ