
книги из ГПНТБ / Лихарев В.А. Цифровые методы и устройства в радиолокации
.pdfАналогичным образом определяются двумерные момен
ты
оо
ь {ЬТа) = 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т„ [у] — 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/и4а |
[У, Г д ] + |
/д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 [у] + |
2т4 ., [у, Гд ] |
|
[У, 2Г Д ] |
+ |
|
||||||
+ 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