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

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

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

сдвига, выходы разрядов которого связаны с логически­ ми схемами, выделяющими четыре возможные комбина­ ции единиц и нулей: 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

р

 

Если входная последовательность описывается 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

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