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

книги из ГПНТБ / Лихарев В.А. Цифровые методы и устройства в радиолокации

.pdf
Скачиваний:
47
Добавлен:
25.10.2023
Размер:
14.41 Mб
Скачать

Аналогичным образом определяются двумерные момен­

ты

оо

ь {ЬТа) = j j " Щ (",. иа .&Гп ) g" (u.) g^u,) du,rfu2 =

—oo

= '«.»V +

^ i S A f ^ X

 

 

/=1

 

0 . 5 r - l

0 , 5 r - l

 

X S

( t - ) S ^ ( " ) -

( 2 Л 5 2 )

где £/^= sin ^/arc c o s

— полиномы Чебышева

второго

рода. Из соотношении (2.151) и (2.152) получим выра­ жения для среднего ти дисперсии а2 , функции корреля­ ции Bv(;kTu) и коэффициента корреляции iRv{kTn) после­ довательности u(kfn):

 

 

0 , 5 г — і

 

т1

= 0 , 5 , о2 =

0,25%•

У| i a r c s i n — ,

 

 

оо

Г0,5г—1

 

Bv ( k r a ) = « п ( /егц ) -

« ; = 8

J ] c/j

( i i )

 

Rv(kTa)=B»{k)

,

(2.153)

Теперь

определим

начальные

моменты v-ro порядка

по формуле (2.150а), когда wi(u)

—одномерное

нормаль­

ное распределение с нулевым средним и единичной дис­ персией. В соответствии с рис. 2.23 и формулой (2.150а) непосредственным интегрированием по каждой подобла­ сти переменной и получим

 

 

00

 

 

 

 

 

 

 

 

 

 

,/z =

J

g

w

(и) w

{u) йи=фо)

4

{F, [(0,5r

— 1) Ди] —

 

 

v

 

 

 

x

 

 

 

 

 

 

 

— 0 0

 

 

 

 

 

 

 

 

 

 

-

F,

[(0,5r

-

2) Дц]} - f {2bv)v {F, [(0,5/- - 2) Ди] -

 

-

F I [ ( 0 , 5 r - 3 ) A u ] } +

... +

 

(iA»)v {FI [(0,5г-і)Ди]

-

 

 

 

 

 

 

 

 

O.Sr—I

 

 

 

- F,

[ ( 0 . 5 r -

 

і -

1)Ди]} =

£

 

( І Д У ) Ч ^

[ ( 0 , 5 г - 1)Аы]

-

ISO

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

г—1

 

 

-

Fl [(0,5r - 1

- І ) Аи}} +

___

( / А у ) " { Л № - О . 5 ' +

 

 

 

f=0,5r

 

 

+

1)Ди] - F 1

{(і -

0,5г)Ди]},

(2.154)

 

 

 

 

X

 

где

Av = ljrl;

F1(x)=

у==

J е - ' 3 ' 2

— функция

—со

Лапласа.

Двумерные моменты найдем, используя разложение двумерной нормальной плотности в ряд '[46] при /Пі=0, а 2 = 1 :

ш 2 ( и і , « 2 ) , ) = | ] ^ Г , , ( " 1 ) ^ Г , ) ( " 2 )

 

оо

со

т^(!гТа)=^

g ^ F ^ ^ d u ,

j _f (ц_) F<1) (u2)du2 +

 

—СО

— 00

+ f l l { I g > i ) f ' ' ) > l ) r f U l

I ^ w ^ w J x

n = l — 00

 

 

III

Учитывая, что

со

со

и интегрируя по частям, получаем

со ґ Г г—1

X

n = i " - » = і

X

151

При изменении индекса

і от (0,5/'+1) до г—1

перемен­

ная

и

изменяется

от Аи до

—0,5г)

Ды. Следовательно,

 

 

 

 

со

|0,5г

 

 

 

 

тщ

(кТа) = т^+

vn £

Щ

(/At>)v-' F\n)[-

(0,5r -

і) Д« +

 

 

 

 

n—l

І = і

 

 

 

 

 

 

г—1

 

 

 

\

1 0,5г

 

 

 

+

J

((Au) v - , F} n >

[(г -

0,5г) Ц

Щ

( / Д ^ - ' X

 

 

 

 

 

 

г—1

 

 

 

 

 

XFln){-(0,5r-i)£ui\+

 

£

( і А о Г - ' Х

 

 

 

 

 

 

і = 0 , 5 ґ + 1

 

 

 

 

Х / ^ ' К і

- 0,5г) Д " ] ! - ^ ^ 1

( - 1 Г ^ ' .

(2.156)

 

На

основе

выражений для одномерных и двумерных

моментов распределений дискретных переменных на вы­

ходе аналого-цифрового преобразователя

(2.151)—

(2.154), (2.156) можно оценить эффективность

систем ко­

герентной фильтрации и режекции. Задача оценки эф­ фективности сводится к определению моментов распреде­ лений на входе порогового устройства при известных мо­ ментах на входе системы и заданном алгоритме ее функ­ ционирования. Наличие двух квадратурных каналов и функциональных преобразователей в схемах рис. 2.22 связано с тем, что разность фаз опорного колебания Ыоп(4) и принятого сигнала неизвестна и является слу­ чайной величиной. При известной разности фаз схема обнаружения полностью известного сигнала не-содержит нелинейных элементов (за исключением порогового устройства). Для упрощения анализа будем считать, что разность фаз известна. В этом случае мы придем к зада­ че обнаружения полностью известного сигнала. Теперь отпадает необходимость в учете нелинейности функцио­ нального преобразования (ФП) . Погрешность в опреде­ лении порогового сигнала при такой замене составляет величину около 1 дБ [57]. Далее примем во внимание тот факт, что ограничение и квантование (в особенности ма­ лоразрядное) расширяют спектр процесса на входе си­ стемы. Это обстоятельство облегчает нормализацию про­ цесса на выходе цифрового фильтра ( Ц С Ф ) .

152

При таких предположениях оценка эффективности си­ стемы рис. 2.22,а сводится -к определению средних и дисперсий нормальных распределений. Пусть ЦСФ задан импульсной переходной характеристикой •h(k), k = 0, 1, 2, ... , п—1. Тогда среднее и дисперсия на выходе ЦСФ вычисляются как

 

 

 

в—і

 

т,

[y\ =

ml[v\ __

h{k),

 

 

 

fe=0

 

 

 

1=0

1=0

 

где Rij — коэффициенты

корреляции отсчетов процесса

і/(-/г), определяемые

по

формулам

(2.153), (2.156). Рас­

пределения смеси сигнала с помехой и помехи на .выхо­ де функционального преобразователя (и детектора оги­ бающей) соответствуют при наличии постоянной состав­ ляющей закону Раиса (2.42), а при ее отсутствии — Релея (2.33). Однако без заметной'погрешности (порядка

0,5

дБ)

можно

использовать

нормальное

приближение.

 

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

обнаружения определяются по формулам

 

 

D=\—F,

 

 

 

где

индексы Sp

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

сигнала с по­

мехой

и одной

помехе; Fi(x)—функция

Лапласа.

 

Приведенная методика анализа справедлива в- том

случае,

если ширина спектра

процесса v(k) на выходе

преобразователя значительно шире полосы пропускания

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

A f > A F = \/пТл,

где п — число

выборок; Гд — период дискретизации.

 

Перейдем теперь к

анализу систем

режекции (см.

рис. 2.22,6). Пусть ЦОФ —система однократного череспериодного вычитания, а ФП — «.вадратор. В этом слу­ чае цифровой накопитель (ЦН) производит накопление величин

z(k) = [v{k)—v{k—\)f.

Вначале найдем моменты последовательности по форму­ лам

/я, [г] = £ (/г), тщ \г\ = zv (/г) z» (k~ 1).

153

Имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

/и, [г] =

2т„ [у] — и

[v, Г Д ] ,

 

 

 

 

 

 

 

пи [г] =

2m., [У] - f 6 л г 2 2 > . Г Д ] 8 т 1 3 [У,

Г Д

] ,

 

 

m 3

[г] =

2/л„ [У] 12//г15 [у, Г д ]

-

 

 

 

 

 

 

 

- 2 0 / я м

[ у , Г Д ] + 3 0 т , , [ о , Г Д ] ,

 

 

 

 

тл

[г] =

2 т 8 [У] 4-70//^ [У, Г д ] +

5 6 т 2 в [ у, Г Д

] -

 

 

 

-

16/га„ , Г д ] -

П2/я„[о,Гд],

 

 

 

п

[г] = 2 т 2 2 [У, Г д ]

/тг23 [у, 2Г Д ]

+

 

 

 

 

 

 

+

m 4 [y] — mlt[v,

Гд],

 

 

 

 

 

 

31

[г] =

/я, [о] +

8/и

[У, Г д ] +

12

[у, 2Г Д ]

-

 

 

 

 

- 6 / и . Л о .

Г д [ - 4 т 3 3 [ у ,

^д] + 0(А 2 1 ),

 

 

 

 

Д2 1

=

6 т . И 1 +

6 т 2 2 2 — 8 т 1 3 2

— 4от2 3 1 ^*0,

 

 

 

 

М

=

/и, [У] +

1 7 т 2 в

[У, Г д ] - ) - т 2

в [У, 2ГД ]

+

+

15/».,, [У, Г д ] -

 

2 6 т 5 3

[У, Г д ] -

8 m I 7 [у, Гд ] +

0(Д„),

 

Д3 1 =

10/лв1, +

1 0 т . ш + 30/я<и —30/яИ 1

— 2/я„а «*0,

т 2 2

(г) =

т 8 [у] +

4 ., [у, Гд ]

 

, Д ]

+

 

+ 12т 2 , [У, Гд ] -

 

8 т 3 5 [у, Гд ] - 8 т 1 7

[у, Гд ] +

0 ( Д И ) ,

 

Д2 2

=

48т.1 2 2

— 1 6 т 3

3 2 — 48/л1 5 2 -(- 16/гг1 И -{-

 

 

 

 

+

1 6 т 1 в 1

— 1 6 1 т 2 3 2 ^ 0 .

 

 

 

 

 

(2.157)

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

копителя с учетом лишь 'парной корреляции:

 

 

 

ті[г/]=лті{г],

 

 

 

 

 

 

 

 

 

m 2

[у]=пт2

[z] + 2(п—\)тц

[z],

 

 

 

 

 

 

ms

[у]=пт3

\z\ + 6 (п— 1) m 2 i [z],

 

 

 

 

 

 

m 4

[у]=nnii [z] + 6 (/г— 1) « 2 2 И + 8 (/г— 1) mla

[г].

(2.158)

Далее следует перейти от моментов (2.158) к куму­ лянтам по формулам (П.1.10) и ряду Эджворта (алго­ ритм 3 приложения 1), с помощью которого и рассчиты­ ваются характеристики обнаружения-

Такая методика применима к Схемам К-кратной ком­ пенсации, когда К не более двух, иначе соотношения бу­ дут очень громоздкими. На рис. 2.24,а, б приведены ха рактеристнкп обнаружения схем одно- и четырехкрат­ ной компенсации, рассчитанные по приведенной методике, а также полученные методом Монте-Карло (сплошные и пунктирные кривые). Как видно из рисунка, погреш­ ность расчетов не более 1 дБ. При этом было принято,

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

функция

сигнала экспоненциальная,

а помехи — гауссова.

 

 

 

При анализе амплитудно-фазовой

цифровой

системы

режекции (см. рис.

2.22,6')

удобно

использовать

метод

следа, рассмотренный выше. Зная корреляционную функ­

цию Bv(kTn)

процесса

v(k) на

входе ЦОФ Bv(kTa)

=

= mu(kTa)

тп], где

mv и т

(кТа) опр?дел ны соотно­

шениями '(2.154) и (2.156), по известной импульсной пе­

реходной характеристике h(k)

ЦОФ

определим корреля­

ционную функцию Вх(кТ„)

процесса

х(к)

[1]. При

дис­

кретном

времени

 

 

 

 

 

 

 

 

 

ВХ (АГП ) = 2

2 Bv

(kTa

- j - /Гц — рТа)

h (1Та) h

(рТп),

где h(l)

определена для

1 = 0,

1, 2, ... ,

К.

 

 

 

Для квадратичного функционального преобразователя

образуется

статистика

 

 

 

 

 

 

 

 

 

 

У = £

z(k)=t

 

 

[(x'(k)r+(x"(k)n

 

 

Кумулянты

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

wi(y)

помехи

и смеси

сигна­

ла с помехой определяются

по формулам (2.135),

(2.136),

(2.141),

(2Л42),

причем

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

матрица

BZn

определяется по корреляционной

функции Вх(кТп).

 

Рас­

чет характеристик обнаружения производится с помощью ряда Эджворта.

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

155

Рис. 2.24. Характеристики обнаружения фазовых ЦЧГЩ: а) одно­

кратная система (2.161а), б) четырехкратная система (2.1616).

По определению, коэффициент подавления помехи ра­ вен [43]

где Ррпхмощность помехи на входе системы подавле­ ния; Ррвых — мощность нескомпенсироваиного остатка помехи. 'В свою очередь, коэффициент подпомеховой ви­ димости определяется как

 

= /Ca/Cs,

(2.160)

где Ks=Ps

вх коэффициент прохождения

сигна­

ла через схему режекции. Рассмотрим цифровые системы режекции с однократным ЧІТВ и четырехкратным ЧПВ,

обеспечивающим равномерную скоростную

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

ку [41]. Соответствующие

алгоритмы

имеют вид

z, [А] = о ( А ) - о ( f t — 1 ) ,

(2.161а)

Zi[k] =

-^-v(k)-v(k-l)-\-b-^v(k-2)-

- v ( k ~ - 3

) - ~ v ( k -

4).

(2.1616)

Можно показать, что коэффициенты подавления вычис­ ляются в соответствии с выражениями [43, 70]

/Спі= 1 0,5[1 - і? г ) п )] - ',

К п 4 = [17,6—11,5 RP п ) 11,1 ЯР (27п) +3,5 Я Р (ЗГП ) +

 

+ 1,5ЯР(4Тп)]-\

(2.162)

где Rp(Tn)—коэффициент

междупериодной

корреляции

помехи на .входе системы

режекции.

 

Коэффициенты подпомеховой видимости для «опти­

мальной» скорости цели равны

 

 

 

К

X + R s

{ Т п )

(2.163)

 

П В 1 ~

1 - ^ ( 7 - п ) '

 

к

_ 1 7 , 6 + 1 1 , 5 ^ 5 ( 7 - п ) - 1 1 , 1 ^ ( 2 Г д ) -

 

1 7 , 6 - 1 1 , 5 / ? р ( Г п

) - 11,1^р (27'п )

+

 

- 3 , 5 ^ ( З Г п ) +

1 , 5 / ? я ( 4 Г п )

 

+ 3 , 5 / ? р ( 3 7 - п ) + 1,5/?р (47-п )

157

Ma рис. 2.25 приведена зависимость потерь в коэффи­ циенте подавления А/Сп однократной системы ЧПВ по сравнению с идеальной аналоговой однократной ЧПВ в зависимости от коэффициента .междупериодной корре­ ляции помехи ЯРа) на входе ограничителя при раз­ личной разрядности преобразования-m. Корреляционная функция помехи соответствовала гауссовой (2.100). На

рис. 2.26 приведены графики для

пв алгоритма

(2.1616)

в функции RP(T„)

при различном m фазовой

системы

режекции. Было

принято, что

выполняются

условия

(2.108), (2.109). Корреляционная функция сигнала соот-

0,2

0,U

0,6

0,8'йрп)

0,2

0,4

0,6

0,в Rp(Tn)

Рис.

2.25.

Потери в

коэффи-

Рис. 2.26.

Коэффициент

подпоме-

циенте

подавления.

ховоіі видимости

фазовой систе­

 

 

 

 

 

 

мы режекции.

 

ветствовала (2.99), а помехи— (2.100). Из рисунка вид­ но, что потери в /Спв растут с уменьшением разрядности преобразования и расширением спектра помехи.

Определяя Коэффициент ПОДПОМЄХОВОЙ ВИДИМОСТИ Кпъ

в фазовой системе, следует учитывать подавление сла­ бого сигнала сильной помехой в идеальном ограничите­ ле. В соответствии с работой [88] величина подавления в зависимости от ширины спектра помехи лежит в пре­ делах от 1 (широкополосная помеха) до 6 дБ (узкопо­ лосная помеха). На рис. 2.26 эти потери не учтены.

158

27

 

 

/77=5^'ЗГА

 

 

 

 

 

 

 

 

Аналог іовая У І

(

о

 

0,75 —o6pao~oi

/ / /

 

чгранича-,

 

 

 

з/А

 

 

 

 

ifh

 

іелем '

0,5

 

 

ill /7/

f /

 

 

Afsrn4f,0S

 

 

 

/ 1 І/

І

 

 

 

 

 

if

 

 

 

0,15

У І І/ 1 / /

 

 

n=tO

//

 

 

 

 

 

//

/

!/\

 

 

F=iO~3

 

 

 

 

K-2

-32

 

-8

0

8 o,35

 

 

-І6

Рис. 2.27. Характеристики обнаружения

систем

двукратной Ц Ч П К

Рассмотрим некоторые зависимости,

полученные с ис­

пользованием описанных выше методик. Влияние разряд­ ности преобразования на характеристики обнаружения системы двукратного череспериодного вычитания показа­ но на рис. 2.27. Указанная на рисунке разрядность соот­ ветствует квантованию нормальных процессов с нулевы­ ми средними и размахом ± З о > . Для перекрытия всего динамического диапазона приемника потребуется боль­ шая разрядность (2.146). В равной мере это заме-

'2

К=і __

 

 

 

-—"Т~

 

 

 

 

 

AfsTn- 0,05

 

 

 

 

п =го

 

 

 

 

 

27= <0,9

 

 

 

 

 

F=l о-5

 

Q.1

12

0.3

О,"

0.5

AfpTn

Рис. 2.28. Пороговые

сигнальї

амплитудно-фазовой

системы ЧПК.

159

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