
книги из ГПНТБ / Лихарев В.А. Цифровые методы и устройства в радиолокации
.pdfПередаточные функции, флуктуационные и динамичен ские ошибки фильтров с эффективной конечной памятью, рассмотренных в § 4.6, приведены в приложении 2. Во просы анализа цифровых фильтров рассмотрены в п. 4.2.3.
3.3.Ц И Ф Р О В Ы Е Э К В И В А Л Е Н Т Ы
АН А Л О Г О В Ы Х И З М Е Р И Т Е Л Е Й
3.3.1.Введение. Синтез цифрового эквивалента ана
логового измерителя состоит в следующем. Для непре-
—V —У
рывной плотности wn(u\X) синтезируется решающая
функция у (и). Она может соответствовать либо решению уравнений максимума правдоподобия (3.1), (3.2), либо уравнению оптимального дискриминатора (3.7), (3.8), либо оператору формирования апостериорной плотности
измеряемых параметров % на данном шаге п с учетом всей предшествующей информации [104]
|
|
00 |
оо |
|
_^ |
^ |
|
wn{l\u) |
= c |
J |
. . . j w |
|
.%i)w[u\Xlt |
Z |
Xn)X |
|
|
— 00 |
— 0 0 |
|
|
|
|
|
|
|
Х Д |
- Д і - , , |
(3.43a) |
||
с последующим |
формировании |
оценки |
|
|
|||
|
|
|
|
|
|
(3.436) |
|
|
|
|
|
Ф |
|
|
|
где до (А! |
|
А„) — априорная |
плотность |
вероятности |
из- |
||
меняющихся |
измеряемых |
параметров; w(u\ A l f . . . , An ) — |
|||||
функция |
правдоподобия; |
с — нормирующая постоянная. |
-*•
По аналоговому прототипу у (и) синтезируется циф ровой эквивалент, т. е. производится замена
- V
Т(ы)=>Т(«).
В устройствах обработки радиолокационной информа ции преимущественное применение находят измерители двух типов. В Р Л С обзора используются многоканаль ные неследящие измерители. Многоканальная структура обусловлена необходимостью измерения координат боль шого числа целей. В Р Л С слежения используются одноканальные следящие измерители.
240
Чаще всего многоканальный измеритель является устройством совместного обнаружения — измерения и строится в соответствии с уравнениями (3.9). Многока нальная система может быть построена различными спо собами (см. п. 2.2.3). При параллельном анализе всей
области Ф определения параметров X, разбитой на М каналов, возможно использование М измерителей или сочетание одного измерителя и адресного или ассоциа тивного запоминающего устройства. При последователь ном анализе области Ф используется один измеритель. Часто применяются комбинированные последовательнопараллельные измерители [98]. В дальнейшем мы будем рассматривать преимущественно один канал многока нального неследящего измерителя. Для следящих изме рителей мы ограничимся рассмотрением цифровых ди скриминаторов, поскольку цифровые сглаживающие це пи могут быть построены по оптимальным линейным фильтрам в соответствии с (3.43), а вопросы синтеза по : следних рассмотрены в § 4.5, 4.6.
Как отмечалось в п. 2.3.1, метод синтеза цифрового алгоритма по аналоговому прототипу имеет ряд недо статков. Тем не менее в двух случаях такие цифровые эквиваленты находят преимущественное применение: вопервых, когда производится обработка некогерентных сигналов; во-вторых, когда по условиям задачи требуется многоразрядное преобразование: при большом динамиче ском диапазоне входных процессов и значительной кор реляции между отсчетами процессов.
3.3.2. Цифровые многоканальные неследящие дально меры и дискриминаторы следящих дальномеров. Анало говыми прототипами цифрового многоканального даль номера в импульсной РЛ С обзора являются индикатор кругового обзора (ИКО) [63] и рециркулятор [50]. Прин ципы организации многоканальной структуры цифровых дальномеров рассмотрены в [98]. В каждом канале даль номера производится обнаружение бинарно-квантован ных сигналов. Обнаружитель каждого канала может быть любым: детектор движущегося окна, весовой би нарный, программный, накопитель позиций, либо весо вой вариант последовательной процедуры.
Остановимся на вопросе точности измерения дально сти такими дальномерами. Дальность цели определяется по номеру канала, в котором произошло обнаружение цели (выполнен критерий обнаружения). С учетом того,
16—1410 |
241 |
что распределение Истинной дальности цели относительно импульсов дальности равновероятно, среднеквадратнческая ошибка оценки времени запаздывания равна
а [ Г ] ^ 7 д / 2 ] / " 3 |
(при Г д ^ т и ) , |
(3.44) |
где Тп — длительность зондирующего импульса. |
||
В первом приближении |
о [Г] не зависит от отношения |
|
сигнал/шум q°-. Для пояснения этого факта |
зафиксируем |
q2 и будем уменьшать интервал дискретности по дально сти Гд. Тогда при Г д < т ц вследствие флуктуации времени запаздывания некоторые импульсы пачки будут попадать в смежные дальномерные каналы, что приведет к сниже нию вероятности правильного обнаружения цели. До тех пор пока можно пренебречь вероятностями ложных об наружений в смежных каналах, точность измерения вре мени запаздывания будет соответствовать (3.44). Нали чие ложных обнаружений данной цели в смежных кана лах увеличивает ошибку измерения времени запаздыва ния. Как правило, в подобных системах не предъявляет ся высоких требований к точности дальнометрии, по
скольку ошибка |
в определении местоположения |
цели |
в большей степени зависит от точности измерения |
ази |
|
мута >8ц и угла |
места 0Ц . |
|
Перейдем теперь к вопросу построения цифровых ди скриминаторов следящих дальномеров. Уравнение опти мального дискриминатора имеет вид (3.8). На практике
находят применение |
дискриминаторы типа |
|
|
||||
|
X = |
K-~L^\l) |
|
= Xu-F{s)^ |
|
(3.45) |
|
где Ь(и\Х) —информационный |
эквивалент |
условного ко |
|||||
эффициента |
правдоподобия А(и\Х) в общем |
случае не |
|||||
оптимальный |
[см. |
(3.9)]. |
Образование |
производной |
|||
в (3.45) заменяется формированием конечной |
разности: |
||||||
F[B) ^ |
L(u*\T-^)-L{u\T |
+ ^) _ |
|
( 3 |
|||
|
|
|
|
-> |
|
|
|
Чаще всего статистика Ь(и\Т) |
является |
корреляцион |
|||||
ным интегралом |
|
|
|
|
|
|
|
|
L(u\T)= |
j |
v(t—T)«оц(t)dt, |
|
(3.47) |
242
где uon(t)—опорное напряжение, согласованное по фор ме с излучаемым сигналом. Структурная схема дискри минатора, соответствующего уравнениям (3.46). (3.47), хорошо известна, она приведена на рис. 3.7. Функцию умножителей выполняют стробируемые усилители (СУ).
При построении цифрового эквивалента аналогового дискриминатора рис. 3.7 необходимо учитывать следую щее. Сигнальный вход дискриминатора нецелесообразно делать цифровым, так как это потребует высокоточного
аналого-цифрового |
преобразователя |
«временной |
интер |
||
в а л — код», не |
включенного в |
замкнутый контур |
следя |
||
щей системы. |
Если |
сигнальный |
вход |
будет аналоговым, |
то и вход опорного сигнала также должен быть |
аналого- |
|||
|
|
СУ |
ґ |
|
|
|
|
|
|
v(t-7) |
J |
|
+ |
eft) |
|
|
|
E |
|
СУ
" L
V, (t-T-1)
Рис. 3.7. Структурная схема дискриминатора задержки.
вым, что потребует включения цифро-аналогового преоб разователя «код — временной интервал» в цепь обратной связи. Перемножители и интеграторы в этом случае вы
полняются |
на цифровых элементах. |
|
|
|
На рис. 3.8 представлена структурная схема цифро |
||||
вого аналога дискриминатора рис. 3.7. |
Перемножение |
|||
принятого импульса v(t—Т) |
и двух опорных |
Vi(t—Т—та), |
||
vz(t—Т+Ти) |
производится |
в двух схемах |
И с |
заполнени |
ем временных интервалов импульсами масштабной ча стоты. Период Тл=Хц/п, где м ^ 8 . Функции интеграторов и вычитающего устройства выполняет реверсивный счет
чик, |
формирующий |
сигнал |
рассогласования. В |
на |
чале |
каждого зондирования |
счетчик устанавливается |
на |
|
0. Для уменьшения |
ошибок |
квантования в дискримина- |
16* |
243 |
торе следует уменьшать период следования масштабных импульсов 7"д, но это требует увеличения быстродействия реверсивного счетчика. Ширина дискриминационной ха рактеристики равна 2тц .
|
Более простои в реализации является схема, пред |
|||||||||
ставленная на рис. 3.9. Короткий |
импульс v\(t—Т—ти/2) |
|||||||||
длительностью т/=т„//г с выхода преобразователя |
«код-— |
|||||||||
временной |
интервал» подается на |
вход линии |
задержки |
|||||||
т п |
с .п |
отводами. Импульс |
цели |
с |
выхода |
приемника |
||||
v(t—Т) |
длительностью |
Тп поступает |
на |
входы |
схем |
|||||
И. В зависимости от величины временного |
рассогласова |
|||||||||
ния |
между |
отраженным |
и |
опорным |
импульсом |
будет |
||||
|
|
|
|
|
г |
|
|
|
|
|
|
|
|
|
|
_П_ |
|
|
|
|
|
|
|
t>(t-T) |
|
Рев. |
счетчик: |
Уст. 0 |
|
|||
|
|
|
|
|
|
|
|
|||
|
|
г |
ЧУ |
|
|
Efk) |
|
|
|
|
|
|
|
I |
v,(t) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 3.8. Структурная схема цифрового аналога дискриминатора задержки.
сформирован код, запоминаемый в регистре. В частно сти, при отсутствии рассогласования в регистре будет код (1, 1, . .. ., 1, 1). Код рассогласования, записанный в ре гистре, с помощью преобразователя преобразуется в двоичный код. В отличие от схемы рис. 3.8 здесь не тре буется быстродействующий реверсивный счетчик. Шири на дискриминационной характеристики равна ти -
3.3.3. Цифровые многоканальные неследящие частото меры и дискриминаторы следящих частотомеров. В им- пульсно-допплеровских Р Л С и непрерывных Р Л С для измерения допплеровской частоты используются набо ры (гребенки) полосовых фильтров [63], перекрывающих диапазон допплеровских частот /вшах—fomin- Число филь
тров M=(fDmax—fuminJ/Aif, где Af — полоса пропускания
244
фильтра. Обычно М составляет несколько сотен. Чаще всего применяются электромеханические фильтры, осно ванные на магнитострикционном эффекте .[115]. Построе ние цифрового эквивалента подобной гребенки в принци пе возможно на основе использования как рекурсивных, так и нерекурсивных цифровых фильтров. Разрядность
v,tt-r-z/2)_ |
1 |
и |
| |
|
oft-Г) |
1 |
1 |
|
|
|
1 * |
* * |
|
|
|
Регистр |
|
Уст О |
|
|
|
|
||
Преобра |
зователь |
е(к) |
||
|
кода |
|
|
Рис. 3.9. Структурная схема цифрового дискриминатора задержки.
аналого-цифрового преобразования должна быть не ме нее 8—12 разрядов.
Пусть каждый фильтр соответствует элементарному звену второго порядка последовательной канонической формы с передаточной функцией /^ - (z - 1 ) (4.21) рис. 4.4,6. Техническая реализация такого фильтра требует наличия
-j(k) і elf)
v,(k)
Рис. ЗЛО. Цифровой дискриминатор системы частотной автопод стройки.
четырех умножителей, двух регистров сдвига и двух трехвходовых сумматоров. Д а ж е при выполнении опера ций в последовательном коде объем оборудования зна чителен. Цифровые гребенки фильтров могут быть кон-
245
курентно-способными по отношению к современным ана логовым устройствам только при использовании больших интегральных схем.
Рассмотрим принципы построения цифрового дискри минатора системы частотной автоподстройки [116]. Его структурная схема приведена на рис. 3.10. Оба входа дискриминатора — сигнальный и опорного напряжения— цифровые. Дискриминатор содержит два перестраивае мых фильтра с передаточными функциями Ki(z) и /Сг(Х),
два квадратора, сумматор и устройство |
усреднения 2, |
с выхода которого и снимается сигнал |
рассогласования |
Рис. 3.11. Полюсы и нули частотного дискриминатора (а) и дискри минационная характеристика цифрового частотного дискримина тора (б) .
e(f) |
=Д1/(£). Передаточные функции фильтров таковы, что |
|||||||
Ki(z) |
Кг(г) — K(z), |
где |
K(z)—передаточная |
функция |
||||
восьмиполюсного |
цифрового |
фильтра |
Баттерворта |
[5]. |
||||
Оба |
фильтра |
настроены |
на |
близкие |
частоты так, |
что |
||
Ki(z) |
имеет |
резонансную |
частоту ниже |
опорной |
ю0 ш |
|||
a Кг(г)—выше |
|
опорной |
ш0п- Расположение полюсов и |
нулей на 2-плоскости изображено на рис. 3.11,а. Четыре нуля расположены так, чтобы сделать фильтр с большей
крутизной |
склонов частотной |
характеристики. |
||||
Передаточные функции обоих фильтров имеют вид |
||||||
К |
Ы |
|
|
K{\-2z-* |
+ z-*) |
|
4 |
, v |
' |
(1 — Л . г - і |
— B , z - 2 ) ( 1 — |
A2z-i—B2z~*y |
|
K |
Ы |
_ |
/ С ( 1 - 2 г - 2 + 2 - * ) |
|||
|
2 v |
' ~ ~ (1 — A3z-і |
— B 3 z - 2 ) (1 - |
Л ^ - 1 — B 4 z - 2 ) ' |
246
где |
|
|
|
|
|
|
|
i4rf = |
2exp(Pf 7"B )cos(Qi 7, „); |
|
|
|
|
||
В І = |
exp ( - 2Р,;7д); |
P, = - |
(Дш/2) sin (ir/8); |
|
|||
Р 0 = |
- ( Д ш / 2 ) |
cos (тс/8); Л = Р 4 , |
^ = |
Л ; |
|
||
Q , = c o o a + P 2 ; |
|
P 2 = (DO U + |
P I ; |
|
|
||
Q3 = |
t o ou — Л ; |
|
Q.i=«>0 u — pi\ |
|
|
||
До>— ширина |
дискриминационной |
характеристики; |
|||||
/ ( — постоянная, |
не зависящая |
от won (при шоп ^Дсо), ее |
|||||
величина выбирается из условия, чтобы |
K\(z)Kz(z) |
обес |
|||||
печивало усиление на частоте со0п, равное единице. |
|
||||||
Коэффициенты ВІ |
зависят |
только |
от Дсо — ширины |
дискриминационной характеристики и вычисляются за ранее. Центральную частоту дискриминационной харак
теристики |
можно |
подстраивать, изменяя коэффициенты |
|
А і, Аг, |
Л з, |
Л4 , для чего требуется вычислять четыре ко |
|
синуса |
и четыре |
произведения. |
Каждый четырехполюсный фильтр реализуется в виде последовательного соединения двух двухполюсных филь тров. Полная схема цифровой системы частотной авто подстройки содержит кроме дискриминатора входной аналого-цифровой преобразователь, цифровые сглажи вающие цепи и вычислитель, управляющий центральной частотой дискриминатор а.
З.ЗЛ. Цифровой фазовый дискриминатор. В аналого вых измерителях функции фазового дискриминатора обычно выполняет фазовый детектор. Для построения цифрового эквивалентна фазового детектора имеется, по крайней мере, две возможности. Во-первых, можно проквантовать выходное напряжение аналогового фазового детектора. Подобная схема используется в цифровых измерителях, когда частоты входного колебания и опор ного напряжения достаточно высоки. В этом случае ис пользуются цифровые сглаживающие цепи и в цепи об
ратной связи включается |
преобразователь код — часто |
|
т а — синтезатор |
частот. |
Во-вторых, можно построить |
цифровой аналог фазового дискриминатора, преобразо
вав в цифровую форму отсчеты |
входного колебания. |
||
Опорное колебание также должно |
быть |
представлено |
|
цифровым кодом. |
|
|
|
Дискриминационная характеристика |
фазового |
детек |
|
тора имеет вид |
|
|
|
Я д [ ' в ( 0 ї = ( 1 / 2 ) и ( * ) О і ( * ) c o s f o t a W - q p o n t f ) ] , |
(3.48) |
247
-где Овх, £>оп— амплитуды входного и опорного напряже ний; (рвх, фоп — фазы входного и опорного напряжений. Соответственно дискриминационная характеристика циф рового фазового дискриминатора
^д[е(/г)]= (1/2т с о з { ф в х ( й ) - ф о п ( / г ) ] . |
(3.49) |
Структурная схема, соответствующая уравнению (3.49), приведена на рис. 3.12. В качестве усредняющего устрой ства 2 удобно использовать простой в реализации опера тор экспоненциального сглаживания. Подобные дискри минаторы могут использоваться на не очень высоких частотах входного и опорного напряжений. Ограничения
vjt,y>sJt)] |
£(*) |
X
Рис. 3.12. Структурная схема цифрового фазового дискриминатора.
связаны |
с быстродействием преобразователя напряже |
н и е — код, умножителя и усредняющего устройства. |
|
3.3.5. |
Цифровые многоканальные неследящие угломе |
ры. В импульсных РЛС обзора измерение азимута про изводится цифровыми устройствами первичной обработ ки. Вся проблема углометрии может быть разделена на ряд задач: а) синтез алгоритма углометрии и определе
ние |
его параметров, б) оценка точности углометрии, |
в) |
техническая реализация алгоритма. |
В плане синтеза данная задача соответствует оценке параметра, постоянного на интервале наблюдения и она тесно связана с задачами обнаружения некогерентных сигналов (см. § 2.3). Статистические характеристики шу мов и эхо-сигналов приведены в п. 2.3.2.
В соответствии с (3.9) вместо максимизации по из меряемому параметру функции правдоподобия или ее информационного эквивалента можно максимизировать коэффициент правдоподобия или его информационный эквивалент. Поэтому для определения алгоритма опти мальной углометрии воспользуемся соотношением (2.58). С учетом (3.9) имеем
Ыр)= |
=Е'^(р)0л-,>с. (3.50) |
/=о |
1=0 |
248
За оценку принимается значение к, обеспечивающее max Ук($). В этом случае принимается р=&.
Фиксация момента времени, когда сумма (3.50) до стигает максимума максиморума, неудобна, поскольку
Считывание
Рис. 3.13. Преобразователь угла поворота в код.
неизвестна его величина. Удобнее определять момент прохождения через нуль выражения
M = |
2 £ G . ( R |
) ^ |
V k . ^ W ' ^ |
V k ^ Q |
(3.51) |
|
/ = 0 |
|
/=0 |
|
|
при |
выполнении |
условия |
(2.58). |
|
|
Входящая в формулу (3.51) функция <?л(|3) опреде |
|||||
ляется в соответствии с (2.41). Таким образом, |
алгоритм |
||||
оптимальной углометрии |
(определение |
центра пачки не- |
флуктуирующих сигналов) при малом отношении сиг нал/шум состоит в весовом скользящем суммировании выходных напряжений квадратичного детектора и опре делении момента времени, соответствующего максимуму суммы (3.50), либо в фиксации момента прохождения взвешенной суммой через нулевой уровень. Импульс от счета азимута подается в схему преобразователя угла поворота в цифровой код — рис. 3.13. Данный преобразо ватель основан на методе считывания И' работает сле дующим образом. На оси вращения антенны закреплен диск (кодовая маска), на котором фотоспособом нанесе на двоичная кодовая комбинация прозрачных и непро зрачных участков. Каждому дискретному углу поворота
?49