![](/user_photo/_userpic.png)
книги из ГПНТБ / Ширман, Я. Д. Разрешение и сжатие сигналов
.pdf(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. |