книги / Проектирование устройств фильтрации радиосигналов
..pdf
|
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 -К ;(у)Н , |
|
|
||
|
где KÎ (у) = 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.