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

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

.pdf
Скачиваний:
16
Добавлен:
22.10.2023
Размер:
13.3 Mб
Скачать
Рис. 1.9.4. Осциллограммы коле­ баний, огибающие которых опре­ деляют амплитудно-частотные спектры радиоимпульсов без вну­ триимпульсной модуляции дли­ тельностью 2 мкс (а) и 8 мкс (б).

тров не требуется, отпадает необходимость повторного гетеродиниро­ вания, а схема (рис. 1.9.3) переходит в рассмотренное ранее видоизме­ нение схемы (рис. 1.9.1, а),’ состоящее в замене полосового фильтра фильтром сжатия. Структура схем не нарушается, если гетеродинное

колебание имеет вид eJ <й/2+ 2яѴ ), а не е‘ьг~, как было принято для сокращения выкладок.

По заданным полосе частот П и длительности анализа Т, исходя

из вида функции

е'6

можно оценить минимальную длитель­

ность

импульсной

характеристики

фильтра Тф =

Т (1

+ Яд/Д/г), кото­

рая соответствует

 

полученным ранее

результатам;

имеется

в виду,

что

тг яй Т.

 

 

 

 

 

 

На

рис. 1.9.4

показаны в качестве

примера выходные

колебания спект­

роанализатора,

построенного

на

основе ультразвуковой линии с есте­ ственной дисперсией, при воздействии на него прямоугольных радиоимпуль­ сов без внутриимпульсной модуля­ ции длительностью 2 и 8 мкс.

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

Еще одной важной особенностью рассматриваемого метода является возможность моделирования не толь­ ко прямого, но и обратного преобра­ зования Фурье. Двойные Фурье-пре- образования с использованием сжа­

тия

открывают интересные возможности обработки сигналов за

счет

дополнительных операций над однократно преобразованны­

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

Некоторые элементы схемы двойного Фурье-преобразования про­ водят сходные между собой операции. В зависимости от выбора знаков производных частот гетеродинных колебаний (одинаковые, противо­ положные) эти элементы можно объединять или аннулировать.

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

§ 1.9.1.

179

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

§1.9.2. ПРИЛОЖЕНИЕ СПЕКТРАЛЬНОГО АНАЛИЗА СО СЖАТИЕМ

КОБНАРУЖЕНИЮ СИГНАЛОВ С НЕИЗВЕСТНОЙ ЧАСТОТОЙ КОЛЕБАНИЙ

Пусть ожидаемое колебание

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

Uq (f) е'2лх/, где %— неизвестная

расстройка частоты относительно

несущей /0. Если U (і) — комплексная амплитуда принимаемых коле­ баний сигнала и шума или только шума, то выражение для модульного значения корреляционного интеграла имеет вид:

1_

Z

2

СО

—со

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

должен

быть

многоканальным.

 

 

 

Во избежание многоканальности можно использовать технику

сжатия.

Повторяя ход вывода

[(4), § 1.9.1],

получаем

 

 

 

Т/2

 

 

 

 

 

Z

I Ui(t) е/м®-о**й

,

(1)

 

 

2

 

 

 

 

где Ux (t) =

U (f)U* (t) — комплексная амплитуда

преобразованных

в смесителе

колебаний; Ur (t) =

U0 (t)e,bt'—комплексная амплитуда

гетеродинного колебания;

Ѳ =

пх/ 6 — временная

координата. Со­

гласно (1) модульное значение корреляционного интеграла вырабаты­ вается схемой, состоящей из смесителя, дисперсионного фильтра и детектора; в зависимости от частоты колебаний меняется положение сжатого радиоимпульса. Схема обеспечивает быстрый поиск по ча­ стоте за время длительности сигнала. В обычном своем виде, без филь­ тра сжатия она не может быть оптимальной*).

Степень неоптимальности выясним путем расчета, пренебрегая расстройка­ ми по времени и частоте.

Пусть принимается колокольный' импульсный сигнал без модуляции (или со скомпенсированной модуляцией) несущей, имеющий длительность тс на уров­

не е“ л/4 яг 0,46, [197]

—я/2/Тс

/2л/с t

(2)

«с (() = е

е

Его спектральная плотность будет

 

 

—л (/—!г)гTg

(3)

Ё с (/) = т о е

 

*> Она неоптимальна и без гетеродина (с растягивающим дисперсионным фильтром [37]).

180

§ 1.9.2.

Сигнал (2) гетеродинируется. Опорными являются колебания

 

 

М 0 = ч Г Л' 2/Т‘-е /Л(2^ + М=)

(4)

с

постоянной

производной частоты

рг. Результирующие колебания

и0 (t) =

=

Uc (0 и* (t)

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

 

 

 

«о {І) =

г

яГ-Щэ /2я/о t

(5)

где /о = /с — /г — промежуточная

частота, а

 

 

 

1/^ э= 1/-^-ИМт;

(6)

 

 

1/т* =

І/т* + І/т*.

(7)

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

Ы /) = т0Эе_Л(?- Ы’т8э.

(8)

Пренебрегая несущественной постоянной задержкой в фильтре, его импульсную V (t) и частотную К (/) характеристики задаем в аналогичной форме

 

 

ѵ ( і ) =

Г

я і Ч ^

е

1 3 я ^

\

 

(9)

 

 

K(f) = xфЭе_ л й

 

 

 

 

(10)

 

 

1/Тфэ =

1/Тф +

/(I*.

 

 

(И)

 

Спектральная плотность выходного импульса gB (f) = g0 (f)K (f)

при /ф =

— /о = /в имеет вид

 

 

 

 

 

 

 

 

 

где

 

£в(/) = Лтвэе — я U —fr,)1тв

 

( 12)

 

 

2

—- гг2

 

т2

 

 

(13)

 

 

 

„ і

 

 

 

 

 

^В Э

 

С0Э

*

l (j)3j

 

 

 

 

 

 

Л— т-оэ Ч-фэ/Т-ЕЭ>

 

 

(14)

а сам импульс будет

 

 

 

 

 

 

 

 

 

 

 

“ в (t) = A e -п іг/ЧЭ

(

 

(15)

 

Величина отношения сигнал/шум по мощности пропорциональна | А |2 и об­

ратно пропорциональна

 

 

 

 

 

 

 

 

 

с о

СО

 

 

 

 

 

 

 

 

 

I

j и а —ѳ) «у (в) de

d t =

1 / У

R e (l/T |3^ l/ T 2) =

(l/ T |-H / T 2) - 1/2.

— с о

— с о

 

 

 

 

 

 

 

 

 

В целом она пропорциональна

произведению | Л |а (1/т|, + 1/т2) 1^2,

которое

в соответствии с (6), (II),

(13), (14)

преобразуем к виду

 

 

 

(І/тФ+ 1/тг )1/_ [ 0 /тс +

1Ітг +

1/тф)2_;К і1ф +

21-1/2

(16)

 

I1г)2]

§ 1.9.2.

181

Выражение (16) приобретает наибольшее значение при согласованной обработке (1/Тф + 1/т£ = 1/т£, (.іф -|- рг = 0), равное 0,5тс. Отношение (16) к величине

0,5то определяет коэффициент использования энергии k. При тг -> оо имеем

k —2/Ѵ (тф/тс^ т с/Тф)--г|-([.іф -ф- |іг)2ТфТс •

(17)

Значение k равно единице только для согласованной обработки. При | {.ігх® |

> 1 ее реализация требует дисперсионного фильтра (рф = —рг 0). Когда дисперсия фильтра отсутствует (рф = 0, а рг Ф 0), подбирая его полосу 1/тф, можно обеспечить максимально возможное значение /гопт < 1. Из условия dkldТф = 0 находим

 

Тф опт =

т с/

-ф- р г Тс

 

 

и

 

 

 

 

 

 

—^

2/(і +

1^"р.г Тс + 1 ) .

 

 

Если при медленном поиске без сжатия (рг Т с Сі ) значение /г0пт ~

1.

то

при быстром (рг Тс 2>

і) значение

копт^І /Ѵ 0,5рг т^ .

 

 

Таким образом,

быстрый поиск со сжатием, для которого игт с2

)ф>

1,

обеспечивает значительно более высокую чувствительность, чем бы­ стрый поиск без сжатия. Преимущества сжатия в отношении разрешаю­ щей способности по частоте [197] были рассмотрены качественно

в§ 1.9.1.

§1.9.3. ОСОБЕННОСТИ АНАЛИЗА ПОСЛЕДОВАТЕЛЬНОСТЕЙ ПЕРИОДИЧЕСКИ СЛЕДУЮЩИХ РАДИОИМПУЛЬСОВ

Пусть на спектроанализатор со сжатием подается последователь­ ность из М импульсов периода Т длительностью ти без модуляции не­ сущей. Гетеродинное напряжение анализатора за время М Т изменяет частоту на Д/г = \/т и (ѵ > 1), его фильтр не лимитирует разрешающую способность.

Длительность «сжатой гармоники» в этом случае 1/ДД, = т„/ѵ менее длительности т„ подаваемого импульса. За время тш соответст­ вующее элементу временного разрешения при согласованной обра­ ботке, можно наблюдать ѵ элементов разрешения частоты по 1ІМТ каждый, просматривая в целом интервал частот ѵ/МТ.

Спектр периодической импульсной последовательности имеет ча­ стотный период 1/7\ Его можно просмотреть за время длительности импульса, когда ѵ ^ М. При этом изменение частоты колебаний гете­ родина за время длительности импульса тн не велико: ѵ « ( 1 к 2)М. Поэтому можно, не нарушая качества обработки, перейти от пилооб­ разной модуляции частоты к ступенчатой, с изменением частоты за период посылки Д/г/М = ѵШт„ ^ 1/ги. Такую же структуру можно придать и импульсной характеристике фильтра, используя линию за­ держки на время М Т со съемом на контуры (см. § 1.4.4).

Устройство спектроанализатора еще более упрощается, если в м е- с т о л и н и и з а д е р ж к и н а в р е м я М Т использовать р е ц и р к у л я т о р с л и н и е й з а д е р ж к и н а в р е м я

182

§ 1.9.3.

Т, Допускающий М эффективных циркуляций одиночного радио­ импульса. Рецнркулятор (рис. 1.9.5) в отличие от обычного совмещает функции накопления и гетеродинирования при ступенчатом законе

развертки

Рис. 1.9.5. Спектроаналпзатор-рециркулятор.

изменения частоты опорных колебаний. Для этого в цепь обратной связи включен однополосный преобразователь частоты с двумя гете­ родинами. Разность частот монохроматических колебаний гетероди­ нов соответствует интервалу между ступенями ѵШг„ ^ 1/тп. На эту

Рис. 1.9.6. Огибающая входного радиоимпульса спектроанализатора-рецирку­ лятора (а), опорное колебание разностной частоты (б) и сжатый сигнал (в).

величину изменяется частота колебаний радиоимпульса при каждой циркуляции. Результаты спектрального анализа воспроизводятся в те­ чение периода колебаний разностной частоты гетеродинов Мтп/ѵ (рис. 1.9.6), длительность сжатого радиоимпульса менее этого периода

§ 1.9.3.

183

в М раз, его положение относительно начала периода (момента прохож­ дения синусоиды через нуль при нарастании) характеризует расстрой­ ку частоты входных колебаний б/. Как видно из рис. 1.9.6, одновре­ менно с частотой можно приближенно измерять время запаздывания радиоимпульса.

По аналогии с рис. 1.4.21 когерентное сложение импульсов раз­ ных частот достигается при определенных фазовых соотношениях. Они обеспечиваются, если разность частот гетеродинов рециркулято­ ра кратна частоте циркуляции 1 /Т. Для этого требуется цепь син­ хронизации (на рис. 1.9.5 не показана).

Обеспечивая когерентное накопление радиоимпульсов с точно неизвестной частотой колебаний, схема рис. 1.9.5 помогает их обна­ ружению на фоне шума. Если ѵ )>> 714, когерентное накопление, охва­ тывая совокупность участков исходных импульсов, каждый из которых составляет небольшую часть полной длительности, должно дополняться накоплением сжатых импульсов, когерентным или некогерентным.

Поэтому слишком большие

значения ѵ ^

1,5 М выбирать нецелесо­

образно.

 

 

 

 

 

 

 

 

В свою очередь, при ѵ = Л4

сжатые

импульсы могут

оказаться

далеко от вершины огибающей

(пунктир

на рис. 1.9.6, в), что

при­

водит к энергетическим потерям.

 

 

 

 

По

принципу

действия

спектроаиализатор-рециркулятор

(рис.

1.9.5)

должен

быть более

чем в М раз

широкополоснее

обычного

рециркулятора,

без

частотного

сдвига в цепи

обратной

связи.

Это

предъявляет повышенные требования к линии

задержки.

 

 

Достоинством рециркуляторного метода является его сравнитель­ ная простота. Существенным его ограничением является малый интер­ вал однозначного анализа частот, который определяется величиной, обратной периоду циркуляции.

§ 1.9.3.

ЧАСТЬ II

ОСНОВЫ ПРОСТРАНСТВЕННО-ВРЕМЕННОГО РАЗРЕШЕНИЯ КОЛЕБАНИЙ

С РАЗЛИЧНОЙ СТЕПЕНЬЮ КОГЕРЕНТНОСТИ

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

1. В развитие анализа временного разрешения когерентных сиг­

налов (гл. 1.1,

1.2) рассматривается

п р о с т р а н с т в е н но-

в р е м е н и 6 е

разрешение (гл. 2.1): угловое, поляризационное, ком­

бинированное (временное и угловое) и т. д. Наряду со случайными амплитудой и начальной фазой учитывается наличие других неизвест­ ных параметров когерентных колебаний, например направления при­

хода [56, 58, 69, 70, 70а, 90, 104,

107,

109, 116, 135а, 136, 143—145,

147, 150, 151, 159, 163, 172, 173,

177,

191 и др].

2. Анализируется разрешение пространственно-временных сигна­ лов с р а з л и ч н о й с т е п е н ь ю к о г е р е н т н о с т и , в том числе обнаружение и разрешение временных сигналов (гл. 2.2). Рас­ смотрение ведется только в приближении гауссовой статистики*). Учитывается изменение неизвестных параметров во времени, требую­ щее самонастройки оптимальных систем обработки [34, 45, 54, 70, 79, 86, 97, 104, 107, 116, 140, 143, 147, 150, 151, 153, 154, 182, 188, 195

и др.].

3. Как развитие разрешения по координате рассматривается р а з ­ р е ш е н и е п о к о о р д и н а т а м и и х п р о и з в о д и ы м (гл. 2.3). Развиваются вопросы обнаружения радиолокационной цели в конечном облаке пассивных отражателей (§ 2.3.5) с учетом скорост­ ных различий [43, 70, 137, 149]. Анализируется влияние пространст­ венной нестационарности облака [137].

Г л а в а 2.1

ЗАКОНОМЕРНОСТИ ПРОСТРАНСТВЕННО-ВРЕМЕННОГО РАЗРЕШЕНИЯ КОГЕРЕНТНЫХ СИГНАЛОВ

Принципы оптимизации совместного приема колебаний в несколь­ ких точках пространства § 2.1.1 — 2.1.2 аналогичны принципам прие­ ма для одной точки, хотя и имеют некоторую специфику. В §2.1.3—

*> С состоянием теории решении для иегауссовой и вообще непараметричес­ кой статистики можно ознакомиться по работам [141, 142, 148, 150, 164, 176, 196, 198, 199, 200].

185

2.1.6 эти вопросы конкретизируются для линейных, поверхностных п объемных антенн. Рассматриваются случаи, когда шумы, на фоне которых разрешаются когерентные колебания: 1) имея внутриприемное происхождение, дельта-коррелироваиы вдоль антенны п 2) когда они порождены внешним изотропным тепловым полем. В §2.1.7—2.1.10 обсуждается разрешение по поляризации, влияние дискретной струк­ туры антенной решетки на угловое разрешение, особенности углового разрешения широкополосных пространственно-временных сигналов и совместного разрешения по направлению прихода и кривизне фронта волны. В § 2.1.11—2.1.12 обсуждается случай, когда мешающий сигнал имеет дополнительный к амплитуде и начальной фазе случайный пара­ метр .

Приступая к чтению §2.1.1, желательно повторно просмотреть

§1.1.2— 1.1.3.

§2.1.1. ОПТИМИЗАЦИЯ СОВМЕСТНОГО ПРИЕМА В НЕСКОЛЬКИХ

ТОЧКАХ ПРОСТРАНСТВА

Пусть в нескольких точках пространства р. = 1,2, ..., р0 располо­ жены приемные антенны (рис. 2.1.1). Для каждой из них задана эф­ фективная площадь приема как функция координат источника излу­ чения и спектральная плотность мощности шума N 0 (последняя счи­ тается одинаковой для всех антенн).

Реализации шумового напряжения uß (t) в каналах будем ха­ рактеризовать совокупностью их дискретных значений по Котельни­

кову для моментов времени t = Д

(^ = 1, 2,

....

Я0).

Совокупность

скаляров р-го канала

образует

м н о г о м е - р н ы й

в е к т о р

ии =

\uß Ді),

Uy. (Д),

...,

(До)) размерности

Д,,

совокупность же

Ро таких векторов (р0А0

скаляров)

образует

вектор

и размерности

р070 =

V. По аналогии с [(1), § 1.1.2] введем с к а л я р н о е

п р о и з ­

в е д е н и е

двух ѵ-мерных

векторов

 

 

 

 

 

 

 

 

|1о Ä-0

 

 

 

 

 

 

 

UV = 2

 

2 (Д) (Д) At

 

 

( 1)

 

 

 

 

ц = 1Х=1

 

 

 

 

 

или в предельном случае At

 

0

 

 

 

 

 

 

 

 

 

 

ДМ-о Г

 

 

 

 

(2)

 

 

 

UV

=

J uß (t)vß (t)dt.

 

 

 

 

 

 

д= 1

 

 

 

 

 

Рассмотрим с о в м е с т н у ю

п л о т н о с т ь

в е р о я т н о ­

с т и

д и с к р е т о в

рс

(ии иг,

..., мѵ) шумовых внутриприемных

случайных процессов (t). Ее можно считать плотностью вероятности многомерного вектора и, все р0Я0 = ѵ дискрет которого независимы и имеют одинаковую дисперсию. По аналогии с [(19), § 1.1.2] получим

рс (и) = (яЛД )-0,5ѵехр( — и2/ЛД),

(3)

где операция и2 = ии определяется соотношениями (1), (2). Пусть наряду с белым гауссовым шумом С действуют т мешающих сигналов

186

§ 2.1.1-

} — 1, 2, .... lit со случайными релеевскнми амплитудами й равнове­ роятными начальными фазами (процесс В) и на их фоне обнаруживает­ ся + 1)-й, полезный сигнал (А). Отношение правдоподобия описы­ вается тогда выражениями [(22), (23), (24), § 1.1.3], векторы Г; и r,-j_

определяются [(15), § 1.1.3].

Во всех указанных выражениях скалярные произведения опреде­ ляются согласно (1), (2). После перехода к комплексным аплитудам

соотношение (2) преобразуется к

виду

 

и ѵ = Т Ë

Re

§ V*(f)VZ(t)dt.

(4)

[i= 1

—СО

 

’>

1

 

’>

______ Устройство

Решение

----------------- >

 

обработки

 

Рис. 2.1.1. Пояснение сущности обработки колебаний, поступающих из несколь­ ких точек пространства.

Если ввести совокупность комплексных

амплитуд

(7), р =

1 ,2 ,...

..., р0, подобную [(33), §1.1.2],

 

 

 

 

Wn (/) = -і- (0 2

S

V* (s) П (S) ds,

 

(5)

X = 2

1 — со

 

 

 

то ей будет соответствовать многомерный вектор

 

 

w = v (uv)-j-vj_(uvj.).

 

(6)

Вектору (6) соответствует, в свою очередь, совокупность комплекс­ ных амплитуд (5).

РІспользуя (4)—(6), от выражений [(15), § 1.1.3] для решающих

векторов г;

и ггі_ перейдем далее к соотношениям для решающих функ­

ций Riii {t).

Аналогично

[(31),

§ 1.1.3],

последовательно увеличивая

число і мешающих сигналов до т,

приходим к рекуррентным соотноше­

ниям =

2, ..., т +

1)

 

 

 

 

Ri»(t) = UiVt{t) + *2

K j R i M

 

 

 

 

/= і

 

 

 

и»

»

 

(7)

 

 

S

J

UlM(s)RfM(s)ds

 

 

 

JLlo

©О

 

 

 

 

 

I UjK(s)RJK(s)ds

 

 

X

1

— CO

 

§ 2. 1. 1.

187

При i = m + 1 они дают окончательные значения решающих функ­ ций, определяющие опорные напряжения каналов при многоканальной корреляционной обработке

R y i ( 0 —

R ( m - j - I) д (0

(р —1, 2,..., р0).

 

О т н о ш е н и я п р а в д о п о д о б и я

[(22)—(24), § 1.1.3]

яв­

ляются, как и ранее,

м о н о т о и и ы м п

ф у н к ц и я м и

в ы-

р а ж е и и и

 

 

 

 

 

г (u) =

urm+1,

 

(8)

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

z (ц) = / ( игш+1)2 + (иг(ш+І)і)2,

(9)

если она равновероятна. Выражения (8) и (9) соответствуют реальной части и модулю к о м п л е к с н о г о к о р р е л я ц и о н к о г о и н т е г р а л а

 

 

 

а д / ? ; ц о

dt —

 

 

 

 

 

Цо

 

«>

 

=

2

Т

J Uß (l)R£(t)dt.

(10)

\1

1

-----ОС

 

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

Значения решающих функций R^ (t) определяются также из системы ин­ тегральных уравнений

Ко

со

qVxtf-

s)RVi(s)ds = Uß (t),

 

у .= 1

j

( И )

со

 

 

 

 

где [і= 1 , 2 ..........р0, а

 

 

 

 

 

Фд* it. s) =

1

т

(12)

2

UfVL (О V]„ (s) .

 

 

0

/ = і

 

Система (11) аналогична [(46), § 1.1.3] и получается из [(45), § 1.1.3], (4) — (6).

Функция (12) — взаимокорреляционная функция комплексных амплитуд т мешающих колебаний в антеннах с номерами р, х для моментов времени t, s, нор­ мированная по отношению к спектральной плотности мощности шума N0. Урав­ нение (11) можно записать в виде

§

Г фд *(^’ s)R.Js)ds = N0Uß (t).

(13)

У . =

1 ---- СО

 

188

§ 2. 1. 1.

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