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

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

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

(hi — неопределенные множители Лангранжа) удовлетворяет диффе­ ренциальному уравнению Эйлера

S v - j j r s » - = 0 ,

( 3 )

а заданные граничные условия соблюдаются. Здесь Sy и Sy' — частные производные:

S y =dSldy, а S y = dSidy'.

Пользуясь условиями (1)—(3), рассмотрим три случая оптимизации спектральных распределений в соответствии со сформулированными критериями. Второй и третий случаи в несколько видоизмененной фор­ ме соответствуют работе Габора [5].

Случай 1. При заданной полосе сигнала П / 2 — Д (в области по­ ложительных частот) требуется найти амплитудно-частотную характе­ ристику сигнала \g(f) |, обеспечивающую м и н и м у м м о щ н о с т и п а с с и в н о й п о м е х и на выходе оптимального фильтра при за­ данной энергии сигнала. Движение отражателей и ограниченная про­ тяженность облака не учитываются.

Соответственно выражениям для мощности пассивной помехи на вы­ ходе оптимального фильтра

У = 2 c \ \ g i f ) \ * d f

и

и энергии импульса

3 = 2 ]\g(n\*df,

где С — постоянная, а множитель 2 учитывает область интегрирования

/ < 0, полагаем в (1)—(3):

У if) — \g if) \ 1 P = 2QA

Q = 2гД

S = 2 (Cif -

Ху*).

Уравнение (3) принимает вид

 

 

 

2Су3 — Ху — О,

 

 

откуда экстремальное распределение

 

 

\g if) Г = У2 if) =

W2С =

const.

(4)

Путем прямой проверки молено убедиться, что распределение (4) соот­ ветствует относительному минимуму функционала С/. Условия на гра­ ницах удовлетворяются только, когда они заданы одинаковыми.

В данном случае оптимальным оказывается прямоугольное спек­ тральное распределение \ gif)\ зондирующего сигнала (см. рис. 1.2.5, а). Оно соответствует наиболее полному использованию его полосы частот П для разрешения по дальности в присутствии распределенных отражателей. Спектральное распределение q (f) выходного сигнала так­ же прямоугольное. Огибающая выходного радиоимпульса W it) опре­ деляется выражением sin [лП it і0)]/л,П (t — ^0) и соответствует

§ 1.6.1.

109

рис. 1.2.5, б. Ширина радиоимпульса по первым нулям огибающей

составляет 2/Я.

 

значении полосы час-

Случай 2. При заданном а б с о л ю т н о м

стот Я

= / 2 — /і требуется найти амплитудно-частотный спектр сигна­

ла I g

(f) I, который

обеспечивает м и н и м у м

э ф ф е к т и в н о й

д л и т е л ь н о с т и

радиоимпульса на выходе оптимального филь­

тра.

 

 

 

Квадрат эффективной длительности определяется аналогично мо­ менту инерции массы, распределенной вдоль прямой,

СО

I со

 

т!=ю* 5

5 \ W ( t f d t .

(5)

Здесь I W (/) I — модуль выходного напряжения оптимального фильтра; к2 — коэффициент; і±— координата «центра тяжести» огибающей вы­ ходного импульса, определяемая из условия

I

а - к ) № ) \> < и = о .

(6)

— со

1

 

 

Не нарушая общности, далее считаем tx

= 0. Чтобы эффективная дли­

тельность т э прямоугольного импульса

была равна его длительности

ти, можно положить к =

J/12.

 

 

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

го радиоимпульса и ее производной

 

 

 

со

 

 

? (/)= 5

 

(7)

 

q'{f) = 5 (—/2л t)w(t)e-i2n^ dt,

(8)

 

— со

 

где

ш(0 = [Щ 0 е12я?°г + W *(*)e-/2*M]/2.

 

При условии оптимальной фильтрации и отсчета времени от центра тяжести выходного импульса | q if) | = q (f), I q' if) |2 = W (/)K Поэтому, используя теорему Парсеваля, получаем

$ |W (0 |2Ä = 2 jj q\{f)df,

со

со

 

 

4я2 J P \ W ( t ) f d t = 2 I

[q'{f)?df,

 

так что

 

 

 

(/c2/4h2) J

[q'(f)?df

5 Q4f)df.

( 9 )

 

 

—co

 

110

§ 1.6. 1.

Интегралы в выражении (9) достаточно вычислить для области частот Д < / < / 2. На краях интервала потребуем, чтобы q (fj) = q (/,) = 0.

Минимизация дроби (9) эквивалентна минимизации числителя при фиксированном знаменателе. Поэтому положим в (1)—(3):

У if) — q (/)> Р = ІУ')\ Q = У1, S = ( y ' f - X y \

Уравнение Эйлера (3) преобразуется к виду

У" + h y = 0.

Решение уравнения будет

y(f) = sin t y T t f - /„ ) ] .

Здесь X и /о определяются граничными условиями у (fi)=y {f2)=0, откуда

?(/) = '/(/) = sin {пт (/ — Д)//7],

(10)

где т = 1, 2, ... В дальнейшем выберем решение, соответствующее т — 1. При этом условии эффективная длительность выходного им­

пульса (9)

оказывается наименьшей

 

 

 

 

 

т э =

тк!2П = кІ2П,

 

(11)

а искомые функции равны

sin Ы {f /і)//7],

 

 

q if) =

 

 

U ( / ) I = K S T №

= 7

№ .

(12)

Выражение

огибающей

выходного отклика

найдем, подставляя

(12)

в выражение для w (t)

[(2),

введение к

гл.

1.6]. Тогда

 

W (t)/W (0) = 0,5 [ф + 0,5) +

ф (Ш — 0,5)]/ср (0,5);

(13)

здесь ф (z) = sin jtz/лг.

Кривые, рассчитанные по формулам (10), (12), (13), показаны на рис. 1.6.1, 1.6.2. Максимальная амплитуда бокового лепестка сжатого импульса составляет 7% от амплитуды основного.

Рис. 1.6.1. Амплитудно-частотные спект­

Рис. 1.6.2. Огибающая выходного

ры входного |g (f)|

и выходного

q (f)

импульса оптимального фильтра при

сигналов оптимального фильтра,

обес­

сигнале с амплитудно-частотным

печивающие минимум эффективной дли­

спектром I g (/) I,

соответствующим

тельности выходного импульса при за­

рис.

1.6.1.

данной абсолютной

ширине спектра П.

 

 

§ 1.6. 1.

111

Случай 3. Требуется найти амплитудно-частотный спектр зонди­

рующего сигнала, который обеспечивает м и н и м а л ь н у ю

э ф ­

ф е к т и в н у ю д л и т е л ь н о с т ь

радиоимпульса (9)

на

вы­

ходе согласованного фильтра. Заданы:

э ф ф е к т и в н а я

п о л о ­

с а

ч а с т о т

энергетического спектра

 

 

 

 

 

 

 

СО

 

 

 

с о

 

 

 

 

=

 

$ ( / - / о ) Ѵ ( М /

$

 

(14)

где /о — координата «центра тяжести» этого спектра

 

 

 

 

 

 

со

 

со

 

 

 

 

 

 

 

h =

I f

Ф (И dff 5

q2 ( f ) df ,

 

 

и

г р а н и ч н'ы е

у с л о в и я

? (/) - > 0

при / ->

±

оо.

 

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

 

 

 

 

 

со

 

 

 

 

 

 

 

 

 

ѵ =

\

W № d f

 

 

 

 

 

 

 

 

— со

 

 

 

 

при дополнительных условиях

 

 

 

 

 

 

 

 

 

 

00

 

 

 

 

 

 

 

 

 

 

1 Я~ ( f ) d f = const,

 

 

 

 

 

 

I

if — f q)2

q2 if) =

const.

 

 

Полагая в (1)—(3)

y(f)

= q(f),

имеем

P = (q')2,

Q0 = q2, Qi =

=

(/ — f0)2q. Уравнение Эйлера (3) приводится к виду

 

 

 

 

q" + \ h

+

h v - /о)2] Я = 0.

 

(15)

После замены переменной

 

 

 

 

 

 

 

 

 

 

/ = */Я2 + /о,

 

 

 

(16)

где

Я42 = —\

и замены

параметра Я0 =

Яа2Я придем

к одномерному

уравнению Шредингера

 

 

 

 

 

 

 

 

 

 

dqVdx2 +

— x2) q = 0.

 

(17)

Уравнение (17) имеет решения, конечные при х ->-±оо, если

Я = 2т + 1.

112

§ 1.6. 1.

• Эти решения имеют вид

Чт(х)*=?гхгі 2Н т{х),

(18)

где Н т (X) — полиномы Эрмита

H m{x) = { - \ y n ^ - - £ ; Z - x \

(19)

Функции qm (х), связанные с полиномами Эрмита, удовлетворяют интегральным соотношениям

со

 

 

со

 

 

^

[q!n{x)]2dx = (m + 0,5)

^

qll (x)dx,

(20)

— СО

 

 

0 0

 

 

со

x2q2m{x)dx = {tn-\-0,5)

со

qll (x)dx.

 

§

^

(21)

— СО

 

 

- - 0 0

 

 

Заменяя х в (20) и (21) по формуле (16),

подставляем (20) в (9), а (21)

в (14). Тогда получим

 

 

 

 

 

4я2 т! = (in

0,5) к2 Хг,

(22)

 

Я І = ( т + 0,5)/с2Д2.

(23)

Перемножая (22), (23) и извлекая квадратный корень, находим

 

 

2 я т э Л в =

( т + 0,5)/ѵ2.

(24)

Наименьшее значение т э при заданной П э получается Для т =

0, ког­

да т эЯ 8 = /с2/4зт. При этом в соответствии с (16)—(19) и (23)

 

q(f) = \g(n\2 = z - Kta- f°)1/4<

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

|g(f)| = e -* ! <f-fo)V8n!.

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

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

2.Наименьшая эффективная длительность тэ сигнала на выходе со­ гласованного фильтра при заданной абсолютной полосе частот П =

= / 2 — Д будет при амплитудно-частотном спектре вида

g ( f ) = V sin [я ( f i — f i W I I .

113

В связи со округленней амплитудно-частотного спектра уровень боко­ вых лепестков снижается по сравнению с предыдущим случаем с 21 до 7%. Проигрыш же по отношению к пассивной помехе невелик. Отноше­ ние мощности сигнала на выходе согласованного фильтра к мощности этой помехи уменьшается на

3. Если задана эффективная полоса частот сигнала на входе согла­ сованного фильтра П а, эффективная длительность выходного радио­ импульса т э минимальна при колокольной форме амплитудно-частот­ ных спектров на входе и выходе фильтра. Боковые выбросы выходно­ го сигнала при этом отсутствуют.

Наряду с амплитудной структурой выходного радиоимпульса пред­ ставляет интерес ф а з о в а я . Если комплексная амплитуда W (t) ра­ диоимпульса является вещественной функцией времени, последний не модулирован по фазе. Это желательно, например, при сочетании уко­ рочения с когерентной техникой. Амплитудно-частотный спектр вы­ ходного радиоимпульса

q{f) = 5 [ W W d W ^ j e - i ^ V d t

—со

симметричен при этом относительно несущей: q(fo + F) = q ( f o - А).

§1.6.2. СКРУГЛЕНИЕ СПЕКТРА И УРОВЕНЬ БОКОВЫХ ЛЕПЕСТКОВ

В§ 1.6.1 были отмечены известные преимущества сигналов с сим­ метричным «скругленным» амплитудно-частотным спектром. Однако специально не прослеживалось влияние скругления на форму огибаю­ щей выходного сигнала и уровень боковых лепестков. Для изучения этого влияния введем аппроксимацию спектра, которая при частных значениях параметров сводится к наиболее характерным: прямоуголь­ ной, колокольной и вида sin х/х.

Пусть спектральная плотность q (/) описывается выражением

 

?(/):

sin aF

0 (ß F + y) — <D(ß.F—у)

( 1)

 

a F

2Ф(у)

 

 

 

где F = f =p

a Ф (x) — интеграл вероятности

 

 

 

Ф(х) = - 4 = [ е - “2/2<іи.

(2)

 

 

 

1Т//2л; о

 

Если ß = 0 ,

тоСпектральная плотность (1) будет

 

 

 

q (f) — sin aF/aF.

(3 )

114

§ 1. 6. 2.

При а = 0 и у —>0

она описывается колокольной зависимостью

-

q(f) = e-V'F>i2)

(4)

 

 

 

 

а при а = 0, y/ß =

Я/2 и у -э-оо — прямоугольной:

 

I 1

для

I/7! <

Л /2,

 

= \0

для

IF I >

/7/2.

Аппроксимация (1) может быть использована как при согласованной, так и несогласованной фильтрации.

Огибающая W (7) напряжения на выходе фильтра (в отсутствие до­ бавочных фазочастотных множителей) сводится к обратному преобра­ зованию Фурье выражения (1). Полагая а = 0 и используя (2), полу­ чаем

W(t) =

"|/2яФ (y)

dF \ exp

- /2яА ^ dx.

(6)

 

—у

 

 

 

-----СО

 

 

Вводя взамен F новую переменную интегрирования

и — ßK + X — /23x//ß,

показатель экспоненты в выражении (6) приведем к виду

+ /2:n;7/ß)2/2 + /2я (и X + /2jt7/ß) 7/ß =

= — и212 — 2xt2/2/ß2 — ]'2 n tx /$ .

Замечая, что

^ e~“°/2 du

-УІТзх,

 

находим

 

 

■W( t) = W a e ~ 2ltJ

rsm (2яу7ф )

(7 )

 

2яy7/ß

 

 

где W0 = 2у/ß<P (у).

Найденное выражение (7) описывает огибающую W (7) как произ­ ведение двух функций времени, симметричных относительно начала отсчета 7 = 0: колокольной и вида sin xlx. В зависимости от соотно­ шения параметров ß и у уровень боковых лепестков в таком произве­ дении изменяется. Максимум первого бокового лепестка при у = 0 соответствует моменту времени 7Х, для которого 2ixy7x/ß « Зя/2, от­

куда 7Хл* 3ß/4y. Колокольный множитель уменьшает уровень

перво­

го бокового лепестка в некоторое число раз

 

р, = е2я*

(8)

откуда

 

у«* 3ß/47x = 3 я Y lg е/2 / 2 lg р,.

( 9 )

§ 1.6.2.

115

На рис. 1.6.3 и 1.6.4 показаны семейства кривых q ([) и

W (/) =

= 1 ^ ( 0 |,

рассчитанные по формулам (1) и (7). При построении кри­

вых было

зафиксировано отношение параметров y/ß = ПІ2

и связан­

ная с ним величина tx = 3ß/4y = 3/2Я. Здесь П — полоса частот для

предельного случая прямоугольной аппроксимации спектра (5), соот­

ветствующей у —>-оо и (.1 =

1. Боковые лепестки на рис.

1.6.4 при этом

не подавляются (ц = 1).

Они начинают подавляться

при конечных

значениях у по мере роста коэффициента ц.

Этому сопутствует скруг-

ление спектра на рис. 1.6.3.

для амплитудно-частотного спектра

Аналогично кривым рис. 1.6.3

выходного сигнала

q (/) на рис.

1.6.5 построены кривые амплитудно-

частотного спектра

|§•(/)] входного сигнала,

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

при согласованной

фильтрации,

когда

q (/)

= |£ ( /) |2. Соответствую­

щие друг другу кривые рис. 1.6.3, 1.6.4,

1.6.5 характеризуются одним

и тем же значением параметра ц. Сопоставляя кривые рис. 1.6.4 и 1.6.5, можно проследить снижение уровня боковых лепестков укороченного радиоимпульса по мере изменения формы амплитудно-частотного спек­ тра входного колебания от прямоугольной к колокольной. Огибаю­ щая амплитуды колебаний укороченного радиоимпульса изменяется при этом от sin х!х до колокольной*’.

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

недисперснонных

фильтров на более короткие

радиоимпульсы

без 4M?

 

 

 

Ответ ясен:

для

этого достаточно приблизить

выходной спектр

к sin х іх или

sin2 х/х2. Изменение выходного спектра в первом случае

можно проследить, пользуясь аппроксимацией спектра (1). При произ­

вольном а

спектру соответствует огибающая

 

 

 

 

 

W ( t ) =

t+а/2л

W 0(s)ds.

 

 

 

 

5

(10)

 

 

 

 

t—а/2л

 

 

 

 

Величина W0 (s) в (10)

определяется

соотношением (7).

Фильтровые

звенья

с

частотной

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

sin х /х рассматривались

в гл. 1.3**’.

 

 

 

 

 

 

*> Из-за округленна

спектра

в приемнике

наблюдаются

расширение

радиоимпульса и потери в пороговом сигнале. По данным [67], где спектр ап­

проксимировался прямоугольником

I F I <

/7/2,

а скругляющий множитель —

выпуклой

внутри него кривой 1 — В2 sin2 (nF/П),

расширение и потерн при

В2=0,92 соответственно будут в 1,5 раза и на 1,3 дБ. Расчетный уровень боко­

вых лепестков составит 0,16%.

Практически он несколько ’ повышается

из-за

неучтенных в [67] осцилляций спектра (см. § 1.6.4 и [102]).

[121].

**> Введение

корректирующего звена sin х/х

описано и в работе

Его импульсная

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

формируется по

многоотводному принципу

(см. рис. 1.3.5) на одной линии с основным звеном фильтра (см. рис. 1.8.9, а).

116

§ 1.6.2.

Рис. 1.6.3. Аппроксимации амплитудно-частотного спектра сигнала на выходе фильтра, промежуточные между прямоугольной и колокольной.

Рис. 1.6.4. Форма огибающей сжатого радиоимпульса при аппроксимациях амплитудно-частотного спектра входного сигнала (см. рис. 1.6.5).

Рис. 1.6.5. Аппроксимации амплитудно-частотного спектра на входе опти мального фильтра, соответствующие рис. 1.6.3.

§ 1.6.2.

117

§ 1.6.3. ВЛИЯНИЕ ОСЦИЛЛЯЦИЙ СПЕКТРА НА ЭФФЕКТ УКОРОЧЕНИЯ

На «скругленный» спектр q0(/) в

ряде случаев накладываются ос­

цилляции бq (/), которые существенно сказываются на

эффективную

длительность укороченного радиоимпульса т 3. Согласно

[(9), § 1.6.1]

последняя определяется отношением

интегралов от квадратов двух

функций частоты /: производной q'

(/)

спектральной плотности и самой

спектральной плотности q (/) = g0

(/) -ф 8q (f) выходного радиоимпуль­

са. При наличии осцилляций это отношение увеличивается.

Увеличение т9 обычно не связано с возрастанием длительности ос­ новного лепестка укороченного радиоимпульса или с изменением его

00

формы. Пусть по условию интеграл от осцилляции § бq (f) df равен ну-

 

 

 

•— со

 

лю, а протяженность ее ограничена областью частот | / — /і | ^

Д/7/2

<<£ П. Тогда приращение огибающей (для /0 =

0)

 

 

 

со

L

 

61*4/)=

J 8 q ( f ) e - M ‘ df

 

(1)

имеет нулевое значение при / = 0 и практически невелико в пределах

всего главного лепестка 11\

<

1/Л. В этом можно убедиться, заменяя

е—/2я// произведением e ~ j2ltlii

■e ~ i 2n( f ~ fi) первый сомножитель ко­

торого выносится за знак

интеграла, а второй близок к

единице,

так что по условию 6W (/)

 

0.

 

 

Поэтому увеличение эффективной длительности вследствие осцил­ ляций 8q (/) связано, в первую очередь, с ростом боковых лепестков укороченного радиоимпульса \ t \ > МП (остатков). Последнее проил­ люстрируем на нескольких простых примерах.

Результирующую спектральную плотность на выходе фильтра q (/) представим как наложение узкого колокольного спектра на более ши­ рокий

q (/о + F) = e -aF' + ce~bF\

где Ь > а (рис. 1.6.6, а).

Огибающая колебаний на выходе фильтра

W(t) = Y

-\гсУпІЬ& -{*іІЬ

имеет при этом протяженные остатки по обе стороны от основного ле_ пестка (рис. 1.6.6, б).

При наличии двух симметричных пиков в спектре q (/) (рис. 1.6.6, в) наблюдается наложение интерферирующих колебаний остатков с оги­ бающей, представленной на рис. 1.6.6, г. Вследствие изменения сдвига фаз колебаний результирующий остаток имеет вид биений

(рис. 1.6.6).

Пусть, наконец, спектр q if) состоит из М пиков одинаковой интен­ сивности. Если основной лепесток укороченного радиоимпульса име. ет интенсивность единица, то каждый пик спектрального распределе. ния дает в него вклад ММ. Обычно вне основного лепестка частные ос.

118

§ 1.6.3.

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