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