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

книги / Проектирование устройств фильтрации радиосигналов

..pdf
Скачиваний:
2
Добавлен:
12.11.2023
Размер:
9.63 Mб
Скачать

 

AI

 

P- (/T j = ф и+1р к / -

1) Tn + m t \ (ф “+l) r + 2

фй° (ф *)г

 

*=i

 

Использование такой

многоэтапной процедуры

позволяет

значительно упростить внутрипачечную обработку сигнала за

счет уменьшения размерности оценки и, кроме

того,

выбора

простого и удобного при реализации весового

коэффициента

Л Упрощение же процедуры обработки сигнала в

пределах

длительности пачки, когда загрузка вычислителя максимальна, позволяет уменьшить требования к его быстродействию. Заме­ тим, что упрощение внутрипачечной обработки сигнала достиг­ нуто ценой усложнения межпачечной обработки и поэтому мо­ жет быть рекомендовано при наличии достаточно длительных перерывов в принимаемом сигнале.

Предложенная процедура обработки пачечного сигнала до­

пускает некоторый произвол в выборе весового

коэффициента

Л

Потери в точности,

обусловленные выбором

неоптимально­

го

для внутрипачечной

процедуры (4.29) коэффициента Л,

компенсируются (в некоторых случаях) оптимальным межпа­ чечным сглаживанием.

Рассмотрим наиболее интересный для практики случай, ког­

да в пределах длительности сигнальной

пачки выполняется

условие Х(6+1) = Х(6), 6 = 1 ,2 , ..., N.

Такое условие часто

выполняется в радиолокационных и радионавигационных си­ стемах с направленным излучением при узкой диаграмме на­

правленности антенны и высоком темпе

обновления

данных.

В соответствии с (4.28) положим т = 1 .

Выберем весовой

коэффициент Л, оптимизирующий внутрипачечную

обработку

сигнала. Из выражений (4.25) и (4.27) получаем

 

 

Лопт(6) = (6 + о2/а2)-|,

/г= 1,

2,

N,

(4.36)

где

ап = С Р -(/7 ’П)Ст — средний

квадрат

ошибки экстраполя­

ции

оценки измеряемой компоненты Х(/)

к

началу

/-й пачки.

Найдем весовой коэффициент Q0nt(i), обеспечивающий опти­ мальное межпачечное сглаживание. Подставляя (4.36) в вы­

ражения для 5, В и R*, получим

 

 

5 =

[1 + о’ДЛ'а»)]-» =

C p-(/7’J Ст[СР-(у) Сг + Я * Г .

В =

[(Л/ + »2/ ^ ) - 1(W + °2hî)~l • • • (# ■+ «’/«S)"1],

 

 

 

 

 

R* = o2//V,

 

a из (4.35),

с учетом

F =

I, найдем

 

 

 

Оопт U) =

P" U T n) с т[CP- (/T J С1] - 1.

(4.37)

Заметим,

что Й0ПГ(У)

не зависит

от R и определяется

исклю­

чительно

соотношениями

между

элементами P“(yTJ:

ш<(/) =

= Рп (У)/Рц(у), / = 1 ,

2,

, л,

где ш, (У) — элементы весового

коэффициента £ï0IIT(y); PZ1(у) — элементы первого столбца ко­ вариационной матрицы P~(JT„).

Рекуррентную процедуру (4.29) можно упростить, если при внутрипачечной оценке не учитывать результаты обработки предыдущих пачек, положив

Л {/е) = 1//г, к = \,

2,

N.

 

(4.38)

В этом случае 5 = 1 , £ = 1 /JV, R* = O2/N

и оптимальный

коэф­

фициент межпачечного сглаживания

 

 

 

 

2о„т U) = P- UTn) с т[СР- (/Т„) Сг + R 'N ] - \

(4.39)

Упрощение внутрипачечной обработки

сигнала

вызывает,

как видно из сравнения (4.37) и (4.39), некоторое

увеличение

сложности межпачечной обработки.

Однако обе

процедуры,

эквивалентные одиночному оптимальному измерению при шуме наблюдения с дисперсией o2/N, также эквивалентны. Таким об­ разом, упрощение внутрипачечной обработки сигнала, ведущее к ее неоптимальности, компенсируется в данном случае опти­ мальным межпачечным сглаживанием, причем вся процедура обработки пачечного сигнала остается оптимальной.

Заметим, что рекуррентная процедура (4.29) при условии (4.38) эквивалентна операции вычисления в пределах длитель­ ности пачки среднего арифметического входного воздействия, которая требует меньших затрат времени, но большего объе­ ма памяти.

Выберем, наконец, наиболее простой, постоянный в преде­

лах длительности пачки весовой коэффициент:

 

A ( k ) = K \,

k — 1,2,

N.

(4.40)

Вычисляя

 

 

,JV

 

5 = 1

- ( 1

- К

 

ХУ .

 

В = [Ал(1 - A ',) 2(iV_1)/CÏ(l - / < 1)2('v-2)

Al],

 

Ki

1 4 - 0 - Kl)A

(4.41)

/?* = , 2 —/Ci i _ (1

-Kx) N »

найдем оптимальный коэффициент межпачечного сглаживания:

Вопт (J) =

Р~0'Тп) Сг [СР-(уГп) Сг +

 

 

Постоянный в

пачке коэффициент

/Ci 0)

можно найти из

уравнения

 

 

 

 

 

а» U) _______________г______________

 

 

~ N [2 - Кх {/Л* I1 ~

A, U l2'4'"1

 

 

1- К\ Ц)

|

К, (/)

 

(4.42)

 

N[ 2 - Ki (у))-

' 2^-K\U)

 

 

минимизирующего дисперсию ошибок измерения в момент окон-

чания пачки [17]. Это уравнение решается методом последо­ вательных приближений до прихода /-й пачки.

Процедура с постоянным в пачке коэффициентом Л отлича­ ется от первых двух рассмотренных процедур большей диспер­ сией эквивалентного шума наблюдения R* и, следовательно, меньшей точностью измерения. Однако потери в точности обычно незначительны. Зависимость R* от Ki в выражении (4.41) допускает следующее приближение (рис. 4.9):

R

N K i < \ ,

(4.43)

W , » l .

г — К\'

 

Поскольку межпачечная обработка сигнала

осуществляется

фильтром полной размерности, то ошибка экстраполяции неве­ лика, так что оптимальный ко­ эффициент К\, являющийся ре­ шением уравнения (4.42), оказы­ вается достаточно малым. Из (4.43) ясно, что при малом зна­

чении Ki, удовлетворяющем усло­ вию NKi<l, процедура с посто­ янным в пачке коэффициентом близка к оптимальной.

Таким образом, многоэтапная процедура фильтрации прерыви­ стых. сообщений с понижением порядка фильтра в пределах дли­ тельности сигнальной пачки по­ зволяет уменьшить требования к

быстродействию вычислителя без существенных потерь в точ­ ности фильтрации.

Та же задача может быть решена несколько другим спосо­ бом. Если в алгоритме (4.24) — (4.27) отказаться от вычисле­ ния весового коэффициента в пределах длительности сигналь­

ной пачки, то загрузка

вычислителя в пачке сведется Только

к реализации алгоритма

(4.24) с постоянным весовым коэффи­

циентом К. Вычисление

постоянного в пачке коэффициента К

осуществляется при этом в паузе между поступлением сигналь­ ных пачек, длительность которой в радионавигационных задачах обычно существенно превышает длительность пачки. В некото­ рых случаях такое упрощение вполне достаточно для реализа­ ции алгоритма фильтрации прерывистых сообщений с помощью вычислителя, имеющего невысокое быстродействие. При этом возникает задача оптимизации алгоритма (4.24) с постоян­ ным в пачке весовым коэффициентом К. Данная задача пред­ ставляет практический интерес также по следующей причине. Оптимизация процедуры обработки пачечного сигнала с помо­ щью соотношений (4.24)— (4.27) приводит к переменному (да­

же в стационарном режиме) оптимальному весовому коэффи­ циенту, так как с приходом каждого нового сигнального паке­ та импульсов в измерителе начинается новый переходный про­ цесс. Оптимизация процедуры (4.24) при постоянном в преде­ лах длительности пачки весовом коэффициенте позволяет опти­ мизировать широко распространенные следящие системы с по­ стоянными параметрами, работающие в пачечном режиме.

Решение задачи в общем случае приводит к громоздкой си­ стеме нелинейных алгебраических уравнений, требующих чис­ ленных методов решения [41]. Конечные результаты удается получить в частном, но практически важном случае, когда со­ стояние объекта в пределах длительности пачки существенно не меняется [39].

Пусть в пределах длительности пачки алгоритм оценивания

описывается уравнением

 

 

 

 

 

 

Х(6) = Х(6 — 1)

К [Х(6) НХ(/г — 1)],

6 =

1,2,

N.

 

 

 

 

 

 

 

(4.44)

Состояния измерителя в

 

А

л

 

Х(у')

/-й пач­

начале Х~(/) и в

конце

ки связаны менаду собой соотношениями

 

 

 

 

 

X” (у) = ФЛ1+1Х ( у — 1 ),

 

 

 

 

Х(У) =

[1-К(у)Н]'уХ-(у) +

 

 

(4.45)

N - 1

 

 

 

 

 

 

+ 2

[> -

к

(у) Н]"-'К (У) У ( N -

1),

 

 

1-0

 

 

 

 

 

 

 

где К (/) — постоянный

в

пределах длительности

j-й

сигналь­

ной пачки весовой коэффициент. Подставляя в (4.45) соотно­ шение У(/) = НХ_ (/) + К(1), где V ( l ) — белый гауссов шум на­ блюдения с дисперсией R = a 2, и полагая при этом, что шумом объекта за время приема пачки можно пренебречь, найдем ошибку измерения в конце приема j-й пачки:

е (У) =

X (у) -

X (у) = [I - К (у) Н |V

(у) -

 

N - 1

 

 

-

2 I1 -

К НИК (У) V (N - /).

(4.461

 

1=0

 

 

Из (4.46) нетрудно получить выражение для ковариацион­ ной матрицы ошибок измерения в конце j пачки Р(/7’п), которая связана с ковариационной матрицей ошибок экстра­ поляции к началу j пачки Р- (уТп) соотношением-

р (у т „ 1 = .[1 - к (У) H f p - (ут,,) {[I _ К ( У ) Н Н Ч

N

- 1

 

+ 2

[1 -К(У )Н ]гК(У)/?Кг (у){[1-К(У)Н]'}Г

(4.47)

1=0

Вычисляя

 

 

 

 

 

 

 

 

 

[1 - К (у )Н ]' =

1 -К ;(у)Н ,

 

 

 

где (у) = 1 ~

U )]lК U),

(4.48)

запишем (4.47)

в виде

 

 

 

 

р (уТп) =

[I -

К/v (У) н] Р- (/Тп) [I -

K:V(у) н ]г +

 

yv-l

 

 

 

 

 

 

 

+ 2

[i — К*(у) H] К (у) R K T(j) [ I - КГ(У)Н]7'.

 

/=О

 

 

 

 

 

 

 

 

 

где J& (/) =

K Î( y) Ly

 

(4.49)

Выполняя операцию суммирования, получим окончательно

P UTn) =

[I -

K;V (у) H] P- (уТп) [I - С

(У) Н]г +

 

 

 

 

+ к ;(у ) /? п к ^ (у ) ]7',

 

(4.50)

где R* определяется из выражения (4.41).

 

что обработка

Сопоставление (4.50) и (4.27) показывает,

"пачечного” сигнала

в соответствии с алгоритмом (4.44)

в рас­

сматриваемом случае

(Ф = 1 и Q = 0 в пределах пачки)

экви­

валентна одиночному измерению с весовым

коэффициентом

Kyv(f) при шуме

наблюдения с дисперсией R*. Это измерение

оптимально при

 

 

 

 

 

 

KN опт(У) =

Р-(уТ„) НГ[НР-(у’7'п) Нг+

/?*]- 1.

(4.51)

Коэффициент КоптО)» оптимизирующий

алгоритм (4.44), опре­

деляется методом последовательных приближений из выраже­ ний (4.48), (4.49) и (4.51).

Увеличение дисперсии эквивалентного шума наблюдения R* по сравнению с величиной a2/N приводит к потерям в точ­ ности оценивания, причем точно таким же, как и многоэтап­ ной процедуре с весовым коэффициентом (4.40).

5.ПРИМЕНЕНИЕ ФИЛЬТРОВ

ВРАДИОТЕХНИЧЕСКИХ СИСТЕМАХ

МЕСТООПРЕДЕЛЕНИЯ ОБЪЕКТОВ

Существенная трудность, возникающая при проектировании! радиотехнических систем местоопределения объектов [18, 27,! 43, 46, 49], связана с вычислительными затратами, необходи­ мыми для реализации оптимального фильтра. Размерность мо­ дели сообщения и размерность пространства, в котором реша­ ется задача местоопределения, обусловливают размерность

.вектора состояния объекта (число интеграторов оптимального фильтра). Количество каналов радиоприема, зависящее от ме­ тода местоопределения и числа радиомаяков, определяет раз­ мерность вектора измерений и влияет на структурную слож­ ность фильтра. Высокая размерность вектора состояния объек­ та и вектора измерений типична для задач местоопределения

.объектов и усугубляется нелинейным характером измеритель­ ных процессов в таких задачах. Возникающая при этом зада­ ча упрощения алгоритма обработки сигналов в радиотехни­ ческой системе местоопределения объектов может быть решена путем применения субоптимальных алгоритмов фильтрации, ■рассмотренных в предыдущих разделах. В данном разделе изла­ гаются некоторые способы упрощения алгоритмов фильтрации с учетом специфики радиотехнических систем местоопределе­ ния объектов. Показана целесообразность построения фильтра с предварительной компрессией данных. Разработан упрощен­ ный метод вычисления параметров фильтра с использованием вспомогательной системы координат. Задача проектирования радиотехнической системы местоопределения объектов рассмот­ рена для случаев использования наиболее распространенных методов: дальномерного и угломерно-дальномерного.

5.1. Проектирование и оценка характеристик устройств вторичной обработки информации в радиотехнических системах

Местоопределение объекта в пространстве осуществляется путем измерения координат объекта относительно радиомая-

ков, местоположение которых известно [25, 28, 43]. В зависи­ мости от выбора метода местоопределения и системы коорди­ нат, в которой оценивается вектор состояния объекта, уравне­ ние измерений, характеризующее измерительный процесс, мо­

жет быть линейным или нелинейным. Если построение оценки

ч

(X) вектора состояния объекта (X) осуществляется в такой системе координат, что измеряются непосредственно компонен­ ты вектора состояния объекта, то уравнение измерений оказы­ вается линейным. Например, угломерно-дальномерный метод местоопределения объекта на плоскости относительно одного радиомаяка с построением оценки вектора состояния объекта в полярной системе координат характеризуется линейным урав­ нением измерений

 

Y(*) =

H(A)X(*) +

V(ft),

 

 

(5.1)

где X(k ) = [р (k ) vp (k )

 

в (/г) v&(k) ... ] т — вектор

состояния

объекта; p(k)

и 0(Æ)— дальность

 

и азимут

объекта

относи­

тельно радиомаяка; vp{k)

и î>e(Æ)— скорости

изменения даль­

ности и азимута; Н(£) =

Н — матрица

измерений вида

 

'1

О

.

0

0

0

. 0

 

 

 

 

Н =

о

.

0

1

0

0 J :

 

 

 

о

 

 

V(ft) — вектор

ошибок

измерений

 

в дальномерном

и

угломер­

ном каналах радиотехнической системы.

Предполагается, что процесс формирования вектора изме­ рений Y происходит в радиотехнической системе достаточно быстро, и вектор состояния X за это время можно считать по­ стоянным. При этом дискретное время k в выражении (5.1) (для сокращения записи) можно опустить.

Если ошибки измерения в дальномерном и угломерном ка­ налах радиотехнической системы некоррелированы, то фильтр, осуществляющий построение оценки вектора состояния объек­ та, распадается на два самостоятельных. В одном из них фор­ мируются оценки дальности и ее производных, а во втором — оценки азимута и его производных.

В общем случае выбор линейного уравнения измерений при проектировании устройства вторичной обработки информации приводит к структурной схеме устройства, показанной на рис. 5.1. В этом устройстве осуществляется раздельная филь­ трация каждого измеряемого навигационного параметра Oj из совокупности НХ= \a\az ... а#]т с последующим преобразова­ нием координат, если это необходимо. Размерность вектора из­

мерений

Y = [ 6I&

••• ôw] r

равна числу измеряемых навигаци­

онных параметров. Векторошибок измерений V=[t>iU2

ол]т

характеризуется

нулевым

средним

£{V }=0 и

диагональной

ковариационной

матрицей

ошибок

измерений

R =£{V V r ).

Каждый

 

 

 

л

 

. . . ] г,

из N фильтров вырабатывает оценку A{=[diai

i= l,2 ,

N (а — производная а по времени).

 

 

В случае выбора декартовой системы координат уравнение из­ мерений имеет вид

bi = р, +

=

у

(л- -

х т )* +

- Урм,)2 + ®<.

(5.2>

где /= 1 ,

2, ....

N-,

х,

у — текущие

координаты

объекта;

хры. и у рм

— координаты t-ro радиомаяка.

 

Для полярной системы координат уравнение измерений име­

ет вид

 

 

 

 

 

 

 

 

b t — Y

Ррм,- +

P" — 2PPM,Pcos (<W/ — + v„

(5.3)

 

 

 

1 = 1 2 ,

N,

 

 

 

 

где p и 0 — текущие координаты объекта; ррм. и 0рм. — коорди­

наты i-го радиомаяка. Нелинейные уравнения (5.2)

и (5.3)'

могут быть представлены, в общем случае, в виде

 

Ь,. = М Н Х ) + УЬ

г = 1 ,2 ......... N,

(5.4)

где h — нелинейная функция;

Н — матрица, состоящая из еди­

ниц и нулей, с помощью которой из вектора состояния объек­

та X выделяется

вектор

координат. Например, в декартовой

системе координат

 

 

 

 

 

НХ =

'1

О

0

0

0

[xvx

■УК.

о

о

О

1

о

 

 

 

Путем линеаризации уравнения измерений (5.4) задача проек­ тирования устройства вторичной обработки информации в ра­ диотехнической системе сводится к проектированию расширен­ ного фильтра Калмана (см. раздел 2.3).

В зависимости от способа ввода информации различают па­ раллельный и последовательный фильтры. В первом случае на вход фильтра (рис. 5.2,а) результаты измерения всех навига­ ционных параметров поступают одновременно (в предположе­ нии о неподвижности объекта за время формирования вектора

измерений Y =

[b|bo... Ь^]т). Алгоритм построения

А

оценки X

по совокупности

измерений

Ь,-, (= 1 ,2 ,

...,

N, с некоррелиро­

ванными ошибками Vj описывается в

этом

случае

уравнени­

ем [66]

N

 

 

 

 

 

 

 

 

X = X - + 2

K/[b / - h i(HX-)]I

 

где Ki — весовой коэффициент (вектор,

размерность

которого1

равна размерности вектора

А

являющийся

решением-'

оценки X),

матричных ковариационных уравнений; Х~ — оценка вектора состояния объекта, построенная к моменту поступления вектора измерений Y.

Соседние файлы в папке книги