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

книги из ГПНТБ / Хомяков Э.Н. Вопросы статистической теории оптимальных измерительных систем. Основание для расчета и проектирования

.pdf
Скачиваний:
15
Добавлен:
23.10.2023
Размер:
7.2 Mб
Скачать

причем коэффициенты A-v\

 

характеризуют степень

взаимной кор-

релнннонпон связи

белых

Ш У М О В :

 

 

 

 

 

<к..

 

 

 

, сч.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В к а ч е п н о

входных колебаний чаше всего выступает

аддитив­

ная

м о т о р н а я

смесь

 

 

 

 

 

 

 

 

О.шзкг)

нс1р;ч 11< > re я

и

более

сложные

ситуации,

когда

где

/ / ( О , 6'('О

 

некоторые

известные матричные

функции

 

Н ч.итном

случае,

входное

колебание

может

быть

скалярной

ф у н ы ш с й . При

'лом

матрицы

И(I)

и S(l)

имени

. л и ш ь

первую

строк\

племени»!),

отличных

от

Нуля.

 

 

 

 

Д л я

решения задач M I I O I

о.мерпой

фильтрации нормальных иро-

iifccon njin наблюдении на всей осп времени или па полуоси (>'<«• Эынаетсн, предпочтительным нахождение переходных функций из уравнения Випер-г -'Хопфм типа (2.100):

Умножая

.сирина обе части

(Гг. 15) на N*J 1 1 вводя функцию

 

 

-1

перепишем

уравнение (0Ц)..18;\Н)

в виде

t

 

 

либо

 

 

р

 

16 2\)

 

 

I по

 

 

Д л я оценки точности

многомерной фильтрации будем

опреде­

л я т ь матрицу

 

 

L

- ^ l.}Ui) »Je ..

СБ 22)

Матричная функции L{t т) удонлетворяет уравнению

ГИД

Оно

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

(2.133) при

 

 

. H i t , ю = N;4

$ ( д . - о = к ; ' ^ t - t y

 

Д л я

решения

уравнения (6.23) используется

метод

неопределен­

ных

коэффициентов.

 

 

 

 

 

 

Д л я

решения

задачи многомерной

линейной ' фильтрации и

экстраполяции

при наблюдении

па конечном

интервале использу­

ем метод Калмана . Наибольший интерес в

задачах

м.югомернон

фильтрации и

экстраполяции

взаимосвязанных

процессов

пред­

ставляют

вопросы

выявления

зависимости

качественных

показа­

телей

фильтрации

от учета "взаимосвязи

оцениваемых

процессов и

помех и вопросы эффективности многоканального наблюдения по

сравнению с автономным

случаем.

В рассматриваемых

задачах статистические характеристики

оцениваемых процессов и помех полагаются известными. Допуская,

что они могут быть известны с точностью

до некоторого конечного

числа неизвестных постоянных либо медленно, изменяющихся

пара­

метров,

можно построить параметрически

адаптирующиеся

фильт­

ры, если

дополнительно

оценивать и эти

параметры .

 

 

6,2.

СОВМЕСТНАЯ

О Ц Е Н К А В И Н Е Р О В С К И Х

П Р О Ц Е С С О В

 

 

 

П Е Р В О Г О

И В Т О Р О Ю П О Р Я Д К А

ПРИ Н А Б Л Ю Д Е Н И И

 

НА

ВСЕЙ ОСИ

В Р Е М Е Н И НА

Ф О Н Е

В З А И М Н О

Н Е З А В И С И М Ы Х

 

 

 

 

 

БЕЛЫХ Ш У М О В

 

 

 

В

этом

конкретном

случае

нмее-м

 

 

 

 

 

 

 

- 4

- 1

'

 

 

 

 

 

YKjU)=b

1

i

 

N .

 

11 >

 

 

 

 

 

 

 

 

 

 

 

- а"

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

No. «tf>

 

 

 

 

 

 

 

 

i

2 ^

 

 

 

 

 

 

 

 

 

I

 

161'

«-Рассмотрим решение уравнения (6.21) в клрссе физически не­ реализуемых 'фильтров, используя преобразование Фурье'

4"UuV) -

 

 

+

N CiuOl

'

(.6.26)

' Вычисляя выражение (6.26)

с

учетом

уравнения

(6.25), найдем

i

 

 

 

 

 

 

K i c o ) = - r

г ^

 

 

 

 

Va.ll)

<3г+ 9 < «

- j

^

n

 

 

 

 

 

 

где

а

 

'^ г

Обозначая виноj .>вские процессы

 

 

 

• u r , t t ) = X , l ^ = X \ t ) ,

 

 

 

з а п и ш е м в соответствии с выражениями

(2.187)

и (G.27) операторы

фильтрации «

спектральной

форме:

 

 

 

 

 

 

 

 

 

г

г

ц

 

 

(.5.10)

 

 

 

 

 

 

 

 

 

 

 

«За +

+

o i

 

 

 

 

^ г

U ^ )

 

 

 

 

 

 

 

 

Анализируя,

в частности,

алгоритм

оценки

процесса

и г с ( Л )

з а м е т и м ,

что если « а и а л

'производной

WV,*0

подавлен

шумом

( N o ^ M ^ ^ o ) , - ™

он отключается. Из''системы

обработки h.2 < ц<,>}.~о\,

'В противоположном случае ( . N 9 2 * M , ^2**"°)оценка процесса

иг2(т_)

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

только

интегрированием

производной.

 

/ ( л я

оценки 'точности

фильтрации

найдем

матрицу

 

 

 

 

 

 

 

 

 

 

• j u < J a 4 ,

1а'ким

образом

 

 

 

 

 

 

 

_ 2

 

П о л а г а я

 

 

 

 

 

 

 

 

 

так

что

 

 

 

 

 

 

 

 

 

 

р а з л а г а я

на отдельные

слагаемые

в ы р а ж е н и я

 

 

 

 

 

1

 

 

 

&

 

 

С

 

 

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

При

этом

 

А

=

С

" 6г

- а>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

6

 

 

 

 

 

 

 

 

 

•ur.

 

 

 

 

 

 

 

Если один из каналов

практически свободен от помех <ч^) ^г"'?0 )>

то

<У^

о.

П'рп опенке utTgOO

интегрированием

ior,U),

когда-

Ы 0 г ь . 5 . 1 | , а

о , ошибка

фильтрации

неограниченно

растет.

Если

подавлен

канал

процесса

u ^ c t )

(при

^ ' о ^ ъ . U<3,—0)>то автоном­

ная

оценка процесса " " ^ I t ) '

характеризуется дисперсией

 

 

 

 

 

 

 

 

 

 

 

' 1 .

 

(ДНО)

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим влияние на точность фильтрации взаимосвязи бе ­

лых

шумов.

 

 

 

 

 

 

 

 

 

ItycTt

 

 

Ы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Но

-

 

 

 

 

 

N

 

(6А г)

 

 

 

 

 

 

 

 

 

 

В

этом случае

имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

Но

=

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а с

учетом

выражения

(6.23)

 

 

 

 

 

163--.

• а к им образом,

 

 

\

В-

из

2 ^

uj

 

 

 

a

(SA4)

•UJ-,

ft

J

Л

 

 

 

 

C 0

+ с0

 

 

0

^

+ S (

 

что совпадает полностью о выражениями

(6.32)

н

(6.33).

М о ж н о

заключить, что взаимосвязь аддитивных

шумов

не

влияет

на по­

тенциальную точность оценки процессов, однако она н а к л а д ы в а е т своп отпечаток на структуру многомерного фильтра. Вместо выра­ жении (6.27) теперь получим

1(0-

jCO 05

Н и ж е синтезируются оптимальные фильтры, оперирующие с вход ными колебаниями на полуоси и на конечном интервале. Это по­ зволит выявить влияние длительности интервала времени наблю­ дения на качество фильтрации .

6.3 С Т А Т И С Т И Ч Е С К И Й С И Н И " » С Т А Ц И О Н А Р Н О Г О

Ф И З И Ч Е С К И Р Е А Л И З У Е М О Г О Ф И Л Ь Т Р А В И Н Е Р О Н С К И Х П Р О Ц Е С С О В П Е Р В О Г О И В Т О Р О Г О П О Р Я Д К А

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

ветствии с этим методом вычислим

9, 3*

сог

где

1ГН

в

В

 

(5.47)

i\\tvvpimy Ij^g-Cjc^) " находим

в виде

(со-jo.)

W > - j 6 )

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

Решение

системы

алгебраических уравнений для неопределен

ных коэффициентов

Ь

"' $

^kS = т \ 2 ) дает

 

 

 

 

' 2 т

дг

^ 7

 

 

 

 

 

 

' 2 2

 

 

 

 

В

=

 

 

 

 

~1Г

 

 

 

Таким образом,

J L

v&.Sb)

Где

Матрица результирующих ошибок фильтрации находится без обращения по Фурье выражения (6.50) с помощью соотношения

16.52)

со —

' 165

известного из теории интеграла Фурье:

2 -

 

 

 

 

16.5J)

 

 

i

,

 

 

 

 

 

При

взаимно независимых

белых

шумах

получим

 

Г ^

2

ь

i

 

 

3 , + <

 

 

 

C6.5V)

При

фильтрации вннеропских процессов

на

полуоси взаимо­

связь аддитивных шумов непосредственно влияет

на дисперсии

опенок:

 

 

 

 

 

о-

ш-

+ 2 ^

СБ.5Б)

'При а—О, сравнивая ' соотношения (6.55) н (6.56) с пыраже ниями (6.38) и (6.39) соответственно, получим

 

 

 

 

 

^

1 ;

 

 

 

е-.

{

4

>

1

 

 

 

 

2£ё_ =

 

 

 

 

Естественно, что дисперсия ошибки оценки процесса с помощью

физически реализуемого фильтра больше дисперсии

ошибки

оцен­

ки процесса

на всей

оси

времени.

Характерно

т а к ж е , что зависи­

мости ошмб' п фильтрации

от величины дг при

-~0 1.>\ог —«О оди­

накова дл я

операторов

с г л а ж и в а н и я

на всей

оси

времени

п им

полуоси.

 

 

 

 

 

 

 

 

 

Рассмотрим подробнее

матрицу

(6.53). I l p u K U - l

и О * - ^ , ^ 0 0

дисперсии не о б р а щ а ю т с я

в нуль и не равны

бесконечности

(но

могут приближаться

к

этим значениям, на границе

области

оире-

Uifi

 

 

 

 

 

 

 

 

 

Делении 5 1 ' %г • ^ г о о ч е В | 1 Л , | ° . если переписать матрицу (G.53)

ввиде

Г г— . I 2q7?

1

i

 

•и

^2

 

В предельном случае очень интенсивного

Шума в канале

процесса

ur\Cv) (,М0,(

Ь) дисперсии ошибок

еще будут конечными, но

если уровень шума

велик в канале w%(t)

( д г - * - о ) т о , как

отмеча­

лось выше, дисперсия измерения процесса t»a(0 растет неограни­ ченно. Это соответствует хорошо известному п технике факту неограниченного накопления ошибки при измерении ф а з ы путем ин­

тегрирования

частоты

 

При совместном измерении д в у х взаимосвязанных

процессов

имеет место

выигрыш » средней" квадратнче-скон ошибке:

 

N2,

 

в".

в , «

в-*

 

 

 

 

 

1

—• -• ,1., I

 

-•-

(6.60)

где

 

 

 

 

 

Если

T O

 

 

 

 

при

имеем •

Q ^ ~ i-

*,

^ 2 ft

 

График зависимостей

^ , ® г о т P и

« Д а "

на

рис. 6.1. Харак -

Т 6 7

 

 

 

 

-i

 

a

 

 

 

 

•0

'CC

PiIC.

G.I.

liiiiuiriiMociu iii-inrpiiiiiWi

n cpciiKM

i;n:i;ip;i i таткой

ошибке

от napaMOipon ft

чrf.мри

COISMITI uoii

фнлыр.шин

-

Riiiie|)OPri;nx ui»>U'"'cnR

испмп/и

и второго поритк;*.

терен рост показателен выигрыша с увеличением степени корре­ ляции белых шумов (<*.--- { ) .

6.4. С Т А Т И С Т И Ч Е С К И Й СИ И ТЕЗ Н Е С Т А Ц И О Н А Р Н О Г О Ф И З И Ч Е С К И Р Е А Л И З У Е М О Ю Ф И Л Ь Т Р А И Э К С Т Р А П О Л Я Т О Р А

В И Н Е Р О В С К И Х П Р О Ц Е С С О В П Е Р В О Г О И ВТОРОГО П О Р Я Д К А

Совокупность вип'еровских процессов первого н второго порядка характеризуется дифференциальным уравнением (Г>.9.4), п котором

П о л а г а я для простоты и 0, запишем сп-.тему дисперсионных уравнений (Г>.9!1)

d

t

 

S <

^

%

 

• d t

 

d t

 

H

^ 2 ' ' I

N

r S D ? 1 '

- ~

~

=2 p ( Д ) - —

- - -

1

(.6.62)

d t

 

^

 

к 0 1

rAQ 2

. У

Система

линейных

уравнений

(5.102)

"имеет вид

 

dx^t)

dt

s ; 4

^ t ) ;

dt

 

t

 

 

dt

 

 

Собственные

значении

матрицы коэффициентов системы (6.63)

равны:

 

 

 

 

16М)

Л г

Л,,

 

Получнн переволную матрицу и выполнив п р е о б р а з о в а н и я со­ гласно выражению (5.100), получим при Рц = 0 элементы матрицы дисперсий ошибок оценки векторного процесса:

1.6.6.5)

М (t)= в •

-

р;-

, — - — , :',

I №

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