
книги из ГПНТБ / Лихарев В.А. Цифровые методы и устройства в радиолокации
.pdfсдвига, выходы разрядов которого связаны с логически ми схемами, выделяющими четыре возможные комбина ции единиц и нулей: 00, 01, 10, 11. По сигналу, соответ
ствующему выделению комбинаций ий (с»(3 = 0,1) вес W ар через блок ИЛИ записывается в накапливающий
сумматор. Накопленная сумма подается на вход цифро вого компаратора (ЦК) с двумя порогами: автозахват и сброс автозахвата. Согласно принятому условию о н'.
„ Сдвиг " из 5лока
селекции
Сброс'
Автозахвап*
Рис. 2.46. Структурная схема алгоритма (2.227).
чале автозахвата схема коммутируется таким образом, что процедура автозахвата начинается с появления двух единиц подряд и первым в накапливающий сумматор за
писывается вес Wn. При выполнении |
'Критерия |
сброса |
автомат возвращается в исходное состояние. |
|
|
Для нахождения вероятностей правильного и ложно |
||
го автозахвата за /г шагов (счет шагов |
ведется |
с появ |
ления двух единиц подряд) необходимо найти распреде ление суммы
п—1 |
|
5 = 2 " л |
(2.228) |
1=0 |
|
где переменные їїі на каждом шаге принимают одно из четырех значений весов: u0 = Wu, ui=Wi0, WiU їїг, іїз, . •.,
210
un-i = W<m, W0i, |
Ww, Wu. Если |
все четыре |
веса отличны |
|||
друг от друга, то переменные |
и* образуют |
одыосвязную |
||||
цепь Маркова с переходной |
матрицей |
|
||||
|
~ Л о |
Рог |
0 |
.° |
|
|
р |
0 |
0 |
Рю |
Ри |
|
|
|
Роо |
Рої |
0 |
0 ' |
|
|
|
О |
0 |
|
р10 |
р1у |
|
Если входная последовательность описывается v-связной цепью Маркова, то порядок матрицы Р определится по формуле
2V + I
Приводимый ниже анализ основан на методе харак теристических функций (см. п. 1.6.2). Условная характе ристическая функция суммы переменных (при условии начала процедуры с двух единиц подряд), образующих односвязную цепь Маркова, равна
< Ы ° ) = Е ^ и Р и Л . . . - л » „ _ . « » е х р { i v |
? 0 |
|||
|
|
|
|
(2.229) |
где S n |
означает |
сумму, |
взятую по всем |
значениям ин |
дексов |
СЙ, аз, • • •, а п , |
"принимающих |
значения 0 или 1; |
|
Pi — вероятность |
появления единицы в стробе автозахва |
та ,(2.225), (2.226); ри— вероятность появления едини
цы после единицы (2. 213) и т. д. |
|
||
Разностное |
уравнение |
(1.71), |
соответствующее |
(2.229), имеет вид |
|
|
|
фп+4+аіфп+з+агфті+г+а^фп+і+а 4фп=0, |
|||
где |
|
|
|
«1=—(<7oo+tyii), |
a2=iqmqн—"ЯІОЯОІ, |
||
|
йз = а4=0 . |
|
|
Окончательное |
выражение для условной |
характеристиче |
|
ской функции |
фп + 2 + аіфі г + і + а2ф7 г + з = 0. |
.(2.230) |
|
|
|||
Условные |
начальные моменты h-ro |
порядка суммы |
|
(2. 228) определяются через |
условную |
характеристиче- |
14* |
211 |
скую функцию в соответствии с (1. 73):
от<°> = « р п ( 0 ) = / > 1 р п , |
|
КН)=р*РпК- |
||||
< ° = АЛ, |
[ |
Л |
. |
+ |
Ри |
( 2 ^ , ) " ] . |
Последовательно |
вычисляя |
|
т ' Л ) , |
/ = |
1 / 2 , п , нахо |
дим искомые моменты tnh суммы (2.228). По формулам
(1.75) можно перейти от начальных моментов к куму лянтам и тогда функция распределения суммы записы вается в виде ряда Эджворта (см. приложение 1, алго ритм 3).
После обнаружения траектории она передается на со провождение, в процессе которого требуется подтверж дение того факта, что сопровождение ведется по отмет кам от цели. Дело в том, что при отсутствии отметок от цели в стробе (при потере цели или сопровождении «помеховой траектории») в алгоритмах экстраполяции ис пользуются координаты центра строба вместо координат отсутствующей отметки. Такое сопровождение может продолжаться неограниченно. Для своевременного пре кращения ложного сопровождения используется блок сброса сопровождения (рис. 2.43, 2.44). Если плотность единиц в стробах сопровождения на нескольких смеж ных циклах уменьшается, то этот факт свидетельствует о выходе цели за пределы строба и сопровождение сле дует прекратить. Простейший алгоритм подтверждения (сброса) сопровождения является детектором движуще гося окна.
Пусть критерием сброса сопровождения является на
личие трех нулей |
подряд. Тогда три разряда |
регистра |
||||
сдвига |
можно рассматривать |
как |
«движущееся окно», |
|||
причем |
в |
начальный момент со |
времени автозахвата |
|||
в первом |
разряде |
записана |
единица, а в момент сбро |
|||
са — во |
всех трех |
нули. Отсутствие |
импульса |
сброса на |
любом шаге после автозахвата означает подтверждение автозахвата (не сброс) на этом шаге. Число возможных комбинаций единиц и нулей (состояний автомата сбро са) в движущемся окне из" трех разрядов равно 2 3 = 8 . За начальное состояние примем комбинацию 100. Строго говоря, к возможным начальным комбинациям относят ся 101 и 111, но содержимое второго и третьего разря дов регистра в начальный момент не влияют на работу 212
автомата сброса, так как состояния последнего опреде ляются последующими комбинациями 'входной последо вательности. Кроме того, после осуществления автоза хвата можно осуществить сброс всех разрядов регистра, за исключением первого. Состояние, определяемое ком бинацией ООО, — поглощающее.
При наличии межобзорной корреляции, задаваемой переходной матрицей (2.213), переходная матрица цепи Маркова автомата сброса имеет вид
0 |
А о Л . О О О |
0 |
о - |
|
||||||
0 |
|
0 |
0 |
р00 |
рп |
о |
0 |
0 |
|
|
0 |
|
0 |
0 |
0 |
0 |
р10 |
A i |
0 |
|
|
р _ Рог |
О |
О |
О |
О |
О |
0 |
Poo |
(2.231) |
||
О |
Рго |
Ри |
О |
О |
О |
0 |
0 |
|||
|
||||||||||
О |
О |
0 |
р„ |
р01 |
о |
0 |
0 |
|
||
О |
0 |
0 |
0 |
0 |
pi0 |
Pn |
0 |
|
||
0 |
|
0 |
0 |
0 |
0 |
0 |
Q |
1 |
|
Для бернуллиевой входной последовательности в (2.231) следует положить poo = P i o = l — р , Poi = Pu = p, где р опре деляется из 1(2.225), (2.226). Аналогичным образом мо гут быть составлены переходные матрицы и для других критериев сброса.
Начальный вектор |
цепи |
есть вектор-строка |
р ^ = ( 1 , 0 , 0 , |
||
0, 0, 0, 0, 0). Вектор-строка |
состояний |
на /г-м шаге |
с мо |
||
мента автозахвата рк |
определяется |
по формуле |
(1.30). |
||
Элемент p f e | 8 вектора |
рк, |
соответствующий |
поглощаю |
щему (в нашем случае восьмому) состоянию, определяет вероятность сброса на k-u шаге. Для неоднородной це пи Маркова с вероятностями, зависящими от номера ша га (см. п. 1.3.2), методика решения сохраняется, только вместо возведения переходной матрицы в k-ю степень необходимо производить перемножение матриц, соответ
ствующих каждому шагу. |
|
|
_ |
||
ДЛЯ |
определения |
Среднего |
времени ДО Сброса Геб |
||
можно |
применить |
метод |
фундаментальных |
матриц |
|
(см. п. 1.5. 1) и воспользоваться |
формулой |
|
|||
|
|
3 ^ = |
2^. |
(2-232) |
213
где in — элементы первой строки фундаментальной мат рицы (1. 49).
Анализ процесса сброса сопровождения затрудняется, когда движущееся окно содержит более трех позиций, поскольку порядок переходной матрицы (2.231) равен 2". Поэтому рассмотрим другую методику получения пере ходной матрицы порядка я > 3 [81]. Под состоянием цепи Маркова будем понимать число единиц в я разря дах регистра сдвига. Такая цепь имеет ( я + 1 ) состояние. Состояния цепи статистически зависимы, поскольку на любых двух смежных шагах имеется (я—1) смежная по
зиция. На каждом шаге возможен переход лишь |
в три |
||||||||||||||
состояния: если |
на v-м шаге |
имелось |
і единиц, то |
на |
|||||||||||
(\>+1)-м шаге может быть: |
|
|
|
|
|
|
|
|
|
||||||
а) |
(і—1)—если |
|
при сдвиге |
одна единица |
вышла |
за |
|||||||||
пределы |
«движущегося |
окиа» |
и не добавилась |
новая |
|||||||||||
единица; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
б) |
і — если при сдвиге |
одна |
единица вышла |
и |
одна |
||||||||||
добавилась, |
либо не добавилось |
и не вышло |
ни одной; |
||||||||||||
в) |
( t + І ) — е с л и |
при сдвиге |
добавилась |
одна |
единица |
||||||||||
и одна не вышла за пределы движущегося |
окна. |
|
|
||||||||||||
Нетрудно |
видеть, что переход |
из |
( я + 1 ) |
состояний |
|||||||||||
в одно из трех состояний на (л>+1)-м |
шаге |
определяется |
|||||||||||||
состоянием |
на v-м |
шаге |
и входным |
воздействием |
на |
||||||||||
( V + 1 ) - M |
шаге. |
Такие события |
|
образуют |
односвязную |
||||||||||
цепь |
Маркова, |
переходная |
матрица |
которой |
является |
||||||||||
якобиевой |
(т. е. трехдиагональной). |
Вектор |
начальных |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
—• |
|
|
вероятностей |
цепи, |
как |
и ранее, |
имеет |
вид р=(\, |
О, |
О, ... , 0 ) т . Элементы переходной матрицы Р вычисляют ся так:
A f . * _ , = ( l |
—р)Рі,п(\. |
я). |
Pi,i = PPi.n(l- n ) + |
(] — Р)П |
—Pt.niUn)], |
Рг,г+1+Р[\-Рг,п(1,>г)}, |
(2.233) |
|
где / ? t , n ( l , n ) — у с л о в н а я |
вероятность наличия единицы |
в я-й позиции, покидающей при сдвиге движущееся окно,
при условии, что во всех я |
позициях |
получено і единиц. |
|||
|
Число |
комбинаций из і |
единиц в |
п позициях |
равно |
С'. |
Если |
же в я-й позиции имеется |
единица, то |
осталь- |
|
ные |
(і—1) |
единиц образуют Cl~J |
комбинацию. |
Тогда |
214
вероятность наличия единицы в п-й позиции при усло вии наличия /' единиц в п позициях равна
|
Рг,п(\,п) = |
С1-__\1Сп. |
(2.234) |
||
С учетом |
(2.234) соотношения |
|
(2.233) принимают |
вид |
|
P/.*-i = |
(l — Р)—— |
Ри=р— |
("(І — Р) |
j—ш |
|
|
К |
г ^ р ' |
^ |
^ 1 - |
(2.235) |
Дальнейший анализ процесса сброса сопровождения аналогичен рассмотренному выше. Ряд соотношений для цепей Маркова с якобиевой переходной матрицей при веден в (102]. В частности, компоненты вектора финаль ных вероятностей равны
|/.=.1 |
V |
S = l (=1 |
/ |
где r S = / e > l . |
|
|
|
Г Л А В А Т Р Е Т Ь Я
ОЦЕНКА ПАРАМЕТРОВ СИГНАЛОВ ЦИФРОВЫМИ УСТРОЙСТВАМИ
3.1.В В Е Д Е Н И Е
Данная глава посвящена вопросам синтеза структуры и анализа эффективности цифровых устройств оценки параметров сигналов, т. е. цифровых измерителей. При
этом |
оцениваемые |
параметры могут быть как постоян |
||
ными |
на интервале |
наблюдения, так и |
изменяющимися. |
|
В последнем случае |
возникает |
задача |
оценки процессов |
|
т. е. линейной или нелинейной |
фильтрации. |
Так же как и при обнаружении сигналов цифровыми устройствами (см. гл. 2), определение структуры цифро вых измерителей может осуществляться тремя методами: построением цифрового эквивалента аналогового алго ритма-прототипа, эвристическим и статистическим син тезом после аналого-цифрового преобразования с исполь зованием методов теории статистических решений.
Наряду со специфическими проблемами, присущими задачам синтеза аналоговых измерителей [7, 42, 43, 69, 104], при синтезе цифровых измерителей возникают до полнительные трудности. К их числу следует отнести сложность и неоднозначность выбора и оптимизации па раметров аналого-цифрового преобразования, а также синтез цифрового измерителя, который при заданной сложности его структуры обеспечивал бы минимум по терь эффективности. Кроме того, возникают дополни тельные трудности, связанные с учетом и уменьшением шумов квантования, ошибок округления, определением устойчивости следящих измерителей и входящих в их состав рекурсивных фильтров.
В эвристическом синтезе используются специфиче ские особенности процессов, преобразованных в цифро вую форму при малом числе разрядов преобразования. Определяя структуру цифрового измерителя, разработчи ки существенно используют накопленный опыт в части упомянутых особенностей. Подобные измерители просты в реализации и зачастую достаточно эффективны, однако вопрос об их оптимальности остается открытым.
216
При статистическом синтезе свойства смеси сигнала J помехой на выходе аналого-цифрового преобразователя рассчитываются по известным свойствам процесса на входе и параметрам АЦП. Последовательность дискрет ных переменных в общем случае может быть описана многосвязной цепью Маркова (см. § 1.7) с распределе нием вероятностей (1.20). Синтез измерителя параме тров, постоянных на интервале наблюдения, сводится к отысканию структуры устройств, определяющих иско мые параметры по максимуму соответствующей функции правдоподобия. В свою очередь, измеритель изменяю щихся 'параметров синтезируется в двух видах: либо в виде устройства, фиксирующего искомые параметры по максимуму (или центру тяжести) апостериорного распре деления измеряемых параметров с учетом предшествую щей информации, либо в виде устройства, усредняющего предшествующие единичные замеры.
Цифровые измерители, синтезированные статистиче ским методом, как будет ясно из дальнейшего, обладают двумя интересными и важными особенностями.
Во-первых, если последовательность дискретных пере менных является цепью Маркова, то структура таких измерителей инвариантна относительно статистических свойств процессов и параметров сигналов на входе ана лого-цифрового преобразователя. Последние влияют лишь на параметры измерителя. Однако при уменьше нии шага квантования Аиа—>-0, а также для слабых сигналов, структуры, полученные методом статистиче ского синтеза и до аналоговому прототипу, совпадают.
Во-вторых, сложность технической реализации, т. е. аппаратурные затраты, оцениваются сразу же, после определения параметров аналого-цифрового преобразо вания, т. е. разрядности т и периода временной дискре тизации Гд. Поскольку структура измерителя и ее пара метры определяются методами теории статистических решений, то при заданной сложности его статистическо го описания такой измеритель обеспечивает максимум эффективности, т. е. минимум потерь информации по сравнению с теоретически оптимальным аналоговым прототипом. Эффективность такого измерителя зависит от оптимизации параметров аналого-цифрового преобра зования. -
Следует отметить, что возможно построение цифро вых дискриминаторов, обладающих непараметрическими
217
свойствами, т. е. структурной инвариантностью и по стоянством параметров структуры относительно стати стических свойств помехи. Последняя должна обладать лишь определенными свойствами симметрии, например иметь симметричную одномерную плотность вероятности. Такие дискриминаторы малокритичны к виду помехи и обеспечивают возможность получения состоятельных оце нок при простой технической реализации измерителя. Все названные методы синтеза находят в настоящее время применение и рассматриваются в настоящей главе.
3.2.СИНТЕЗ И АНАЛИЗ О П Т И М А Л Ь Н Ы Х
ИК В А З И О П Т И М А Л Ь Н Ы Х И З М Е Р И Т Е Л Е Й
3.2.1.Синтез оптимальных и квазиоптимальных изме рителей. В л. 2.2.1 отмечалось, что в практике радиоло
кации всегда осуществляется сложное измерение, т. е.
|
S(t, |
-> |
—> |
|
|
|
|
обнаружение сигнала |
а, |
%) с оценкой его неизвест- |
|||||
ных параметров л.= (Лі, Ля, ... , |
fa)т, |
где |
а = (кц, |
... , ат) |
т— |
||
вектор неизмеряемых |
параметров. Параметры |
Л и а |
за- |
||||
|
|
-»• |
|
-> |
|
|
|
кодированы в сигнале |
S(t, |
a(i), |
K(t)), |
как правило, |
не |
линейно. При синтезе измерителей чаще всего предпола
гается, что принятая |
реализация |
соответствует смеси |
|
сигнала, шума n(t) и'помехи |
p(t): |
|
|
u(t) = S(t, |
a{t), i(t)) |
+ |
n{t)-\-p{t). |
Начнем рассмотрение с задач, когда вектор измеряе мых параметров X постоянен на интервале наблюдения. Подобные ситуации возникают, например, при измерении дальности и угловых координат в РЛС обзора, когда за время облучения цели Т о б л ее перемещением можно пре небречь, а также при оценке флуктуирующих параметров
Л, если время наблюдения Та |
меньше |
времени корреля- |
||||
ции т к |
каждой из |
компонент |
вектора |
Л=(Л,і, Ла, ... , Л()т . |
||
Для |
решения |
задачи |
синтеза должны быть |
заданы |
||
плотности вероятностен |
Wi(k) |
и wm(a), |
а также |
стати |
стические свойства помехи. Синтез оптимального измери теля состоит в формировании функции правдоподобия
218
измеряемых параметров |
|
|
|
|
|
|
wn (и | Я) = |
доп (и [ a, X)wm |
(a) |
da |
|
|
'А |
-> |
|
|
|
|
|
|
-> |
|
|
или. ее |
информационного |
эквивалента / (и | Я) |
на всей |
||
области |
Ф возможных значений вектора |
Я £ Ф |
и опреде |
ления вектора оценки Я, соответствующего максимуму
максиморуму функции правдоподобия. Неизмеряемые параметры а исключаются соответствующим усреднением.
Иными словами, устройство оценки должно формиро-
вать эффект Я, соответствующий уравнению максималь ного правдоподобия:
доп (ы | Я) = max {wn (и | Я)}. |
|
(3.1) |
|
Данный метод синтеза называется методом максиму |
|||
ма функции правдоподобия и |
его важным |
свойством |
|
является инвариантность точки |
максимума (или |
макси |
|
мума максиморума) функции |
праводоподбия |
к |
произ |
вольному взаимно однозначному преобразованию функ ции правдоподобия. Это обстоятелг г.во упрощает струк туры измерителей. Дело в том, і о непосредственная интерпретация соотношения (3.1) связана с построением
многоканальной системы, каждый канал |
которой наст |
роен на определенное значение вектора |
параметров %І, |
/ = 1 , 2 , . . . , М. |
|
->
Определение вектора Я, соответствующего уравнению
(3.1), связано |
со |
сравнением |
выходных |
эффектов |
/ (и | Яг-) от М |
каналов |
и выбором |
наибольшего |
выходного |
эффекта. Поскольку выходные эффекты каналов не нор мированы и изменяются в большом динамическом диапа зоне, то построение решающего устройства наталкивает ся на значительные трудности. Аналогичное затруднение возникает в измерителях времени запаздывания или углового положения, когда необходимо фиксировать мо мент времени Т, соответствующий максимуму функции правдоподобия.
219