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

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

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

ма имеем

теперь

уравнение

Вольтерра.

 

на^-ос.) по­

Заменяя

нижний

предел

в соотношении (5.19)

ручим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

М

* * —

 

 

<•"»

При этом функция

Q ^ i o - ) , как и

Нд'ца)

имеет полюсы только в

верхней

полуплоскости ш.

 

 

 

 

 

Таким

образом,

 

соотношения (5.20), (5.17), (5.15),. (5.12)

дают

решение

задачи синтеза

физически реализуемой системы

оценки

процесса

 

при т =-э :

 

 

 

 

 

 

 

 

г

-

 

Т

п

 

 

'•

 

 

Оценка

процесса

k(t)

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

дисперсией

коррелиро­

ванной

составляющей

 

 

 

 

 

 

^ »

 

j

— = ~

t i m .

 

5

 

V.J*C^

AW

 

Нг 0

 

 

 

\\\

 

 

и спектральной плотностью мощности остаточного белого шума:

Н

В случае единственной реализации

j

имеем тривиальный результат;

W ) = ^ o o ;

.

Г / \ a t - t ) = K V 4 i - t ) + r 4 , . - B ^ - t \ i * * ;

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

Введем коэффициенты эффективности многоканальной обра­ ботки соотношениями

130'

Set)

Рассмотрим некоторые

частные случаи.

6.8. ОПТИМАЛЬНЫЕ

АЛГОРИТМЫ ОЦЕНКИ ПРОЦЕССА

Д Л Я НЕКОТОРЫХ СЛУЧАЕВ КОНКРЕТИЗАЦИИ

КОРРЕЛЯЦИОННЫХ

ФУНКЦИИ АДДИТИВНЫХ ПОМЕХ

Рассмотрим вначале часто встречающийся на практике случай, когда шумы n(t) некоррелнрованы во времени. При этом

Обработка входных колебаний сводится к взвешенному усред­ нению:

SKI**,-

- ^

"

°

(

Если шумы в каждом канале взаимно некоррелнрованы, то матриц» V днагоналъна;

^ 0

так что

п.

- г

.

И

в-.'

 

kz.\

k

 

131

 

 

 

 

(5.33)

•При

^ -

- • • • = <з\г =

<r

 

имеем

среднее

арифметическое;

 

 

 

 

.

S e t - * } •

CS.'lO

 

 

h.

 

 

Д л я

коэффициента эффективности в случае (5.31) имеем вы­

ражение

 

 

 

Пусть теперь коррелированная составляющая

шумов

представ­

л я е т

собой n-мерный стационарный процесс с

полной

.корреля­

цией,

когда

 

 

К СО) =

-

 

 

 

 

(5.56)

 

 

 

 

 

 

 

 

V

^

"

г- 2

 

 

 

Используя

соотношения (5.12) и

(5.13),

получим

 

 

 

 

г

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

(5.Ь8)

 

к,е=ч

w.e

 

 

 

. Функция Н\Ц—т), необходимая

при решении

уравнения (5.19).

имеет следующее изображение

по

Фурье:

 

 

 

Н1 t i w ) =

 

2 < С 0

£

 

 

 

С1+ a") ( a *

j c o f 0

)

причем 1

 

 

 

 

 

 

 

(5АО)

кд,)Л=г *1

и

В соответствии

с

(5.20) получим

 

 

 

 

 

Не

 

С5.М)

 

 

 

 

 

 

Таким

образом,

 

 

 

 

 

 

 

 

 

 

 

 

Если компоненты

процесса

n(t)

некоррелнрованы м е ж д у

cow

Сой и

 

 

 

 

.

 

'

то

 

 

 

 

 

ч5ЛЭ)

 

 

 

 

 

 

 

 

 

 

14 u»t,vт

1 ;

1.5 Л V)

 

 

 

 

 

 

Прй~этом

 

 

 

Ко

 

 

 

 

 

 

 

 

Рассмотрим также случай,

когда "v/О

многомерный

вине-

ровский

процесс

первого порядки

с-полной

корреляцией

 

Такой процесс получается при идеальноминтегрировании вектор­ ного белого шума с корреляционной мйтринеб вида

> | 1 Л ; ...

3

Ч5А7).

- - • »

Соотношения метода неопределенны* коэффициентов приво* дят к следующему результату:

15Л&)

J . 1 3 3

Д а л е е имеем

d

 

 

-4

 

d = S

К

> Г а Т к / ( о

15.50

d

Если положить, что белые шумы, формирующие винеровский процесс, независимы и имеют одинаковую спектральную плотность

мощности, то при

К

= ~=— §

получим

- \ l - 7 ~ , d = — - , <3\

Величина Я?,

зависи . от времени, поскольку

дисперсия

вине-

ровского процесс? первого порядка, определяющая

погрешность

оценки процесса

при

единственной реализации,

равна

В • t.

 

В рассмотренных

з а д а ч а х при групповой оценке процесса

основ­

ную роль играют операции взвешенного усреднения по а н с а м б л ю

реализаций . Л и н е й н а я . фильтрация

учитывает

корреляционные,

свойства помете. Представляет интерес

обсудить

вопросы, связан ­

ные с использованием априорной информации относительно оцени­ ваемого процесса.

6.3. ИСПОЛЬЗОВАНИЕ АПРИОРНОЙ ИНФОРМАЦИИ _ПРИ_ХРуППОВОЙ ОЦЕНКЕ ПРОЦЕССА

П р и конечных интервалах времени наблюдения," как правило, допустима аппроксимация оцениваемого процесса полиномом

ы' •

1 - 0

134

Л о г а р и ф м

функции правдоподобия векторного

параметра

^ +

• • • •> л w "\ ' при наблюдении аддитивной смеси

(5,0 имеет

вид

 

Дифференцируя по параметрам A-v , получим систему уравнений максимального правдоподобия:

N

где

<*.. - U £ b ' . U ^ t V dt dt ;

J

1

0

"l

I,к.енPH

 

6

 

 

 

т т

rt

a a

Вводя матрицу Ф, состоящую из элементов ^ t ] , н ей о б р а т н у о , допишем соотношение для вектора искомых оценок:

 

 

X = ф " 1

• V .

1.5.58)

В простейшем случае, когда

допустимо полагать

>,(Д)v x „ струк­

тура обработки

определится

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

 

 

\ \

£' к

U , t ) 1 W d t d t

 

 

Т Т

п-

 

 

15.S9)

0

 

 

 

Дисперсия

оценки

среднегоравна

 

В частных случаях (5.28), (5.31) и (5.34) имеем соответственно

ч

J t 4

« H ^ d t

 

«

 

^ о = - " 1 - ^

Г 5 — ^

J

\

- Т

S N ;

^

 

 

у ' 1

 

 

• цен

 

"кС

*

I

1 (

г

При едит'ственнон оеализлции

дисперсия оценки среднего значения равна

 

<*f * -

И ^;; Ct,Tr)dt cjt

 

(5Л5)

 

-J

 

i t

 

u .

 

 

 

 

 

 

Эффективность групповой

обработки

при этом характеризуется

выражением

 

.

 

 

 

 

 

 

 

 

т т

 

\

 

 

 

 

 

 

~

JJ £ . Л ( t ^ d t d t

 

($.6Ь)

 

 

а 0

у

н

, ц

 

 

 

 

 

3 j

 

П \ - . 4 t , x ) d i d r

 

 

:B

простейшем

частном

V o

iJ

(5.63) Я з " * 1

 

 

случае

процесса. Ис­

 

Рассмотрим

теперь

задачу оценки нормального

пользуя выражения

(2.4) и (5.1), запишем слагаемые

логарифма

функционала отношения правдоподобия,

зависящие от

\(t):

 

 

 

 

тт

 

 

 

 

 

 

 

 

т т

 

 

 

 

 

 

 

 

 

- ; ^ ( . r ) ] d t d t - \ \ \ - b t & f t

 

*&)\ ^

d't .

 

 

0 0

 

 

 

 

 

 

 

 

 

Структура алгоритма оптимальной оценки процесса

определит­

ся

выражением

 

,

 

 

 

 

 

 

 

= j № ^ ) d t

 

+ 5

X

„ (t/t) u (*) dt

,

(5.68)

 

• о

 

-

o

.w,««V

< C

C

 

 

T

i';36

а функция M ' l j t ) ' определена уравнением

Для средней"квадратической ошибки Оценки процесса, как прежде, имеет место соотношение

< /t< Л ) - U t , °^ - .

Для нахождения функции Ц(,т) используется интегралЬНОР уравнение (2.133) с учетом соотношения (2.132), в котором еле дует полагать

Рассмотрим для" иллюстрации случай некоррелированных флук­ туации, когда

 

• H < A , * v i i ' ^ c H t - * ) .

 

w >

При этом соотношения

(2J32) и (2.133) примут следующий

вид:

 

^ ^ • O - ' ^ N U t t o

•,

^

 

t5.1V)

 

т

 

 

 

 

 

 

 

 

2 • J u (Л (&Л (5,®) d? + L (,i,<*) - V СЛ.»),

( д а

где"

*

'

Л

 

 

 

 

 

При

единственной

реализации

 

 

 

 

 

У ..It)- W

+ П . Д Ц , -

 

(5Л7)

«огда

 

 

 

 

 

 

 

 

 

< a , <Д) п. д*)

> -

\ .

8. U

- О

,

CS.7&)

ошибка

оценки нормального

tfpeaecca равна

 

 

 

a

t

 

 

 

 

 

 

 

<*.(Д) - V C t Л " ) " \ ^ , ^ Ж ^ Л )

d < c

>

^.79)

 

 

0

_

_ -

 

 

 

 

' Де

^3- ( . t . ' t ' )

является

решением

уравнения

 

J

/ Г37.

J U C t , ^ 4 ^ O d ^ + V n M i l . t » = M U , f ) ,

15.80)

о

Эффективность групповой оценки нормального процесса опре­

деляется

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

 

 

 

 

 

 

 

 

 

г Э

1 Д 1 Д )

 

 

 

П о л о ж и м

Д Л Я

примера,

что

вннеровский

процесс пер­

вого

порядка

с корреляционной функцией

 

 

 

 

 

Ж * Л ) - В • m i n C . t . ' c )

 

 

Д Ш )

Интересуясь

стационарным

решением уравнения (5.75),

найдем

частотную

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

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

фильтра

 

Т.1К

Ч Т О

 

 

 

 

 

 

 

 

И

случае

фильтрации К(() по

единственной реализации

ч.(Л)

. л ь .

 

 

 

 

о

 

 

 

^

имеем

 

 

 

 

 

 

 

 

 

Исполь.«уя урапнения (5.81) и (5.85),

получим

 

 

11с л и

 

 

 

 

 

 

 

 

то

 

 

 

 

 

 

 

 

 

 

 

 

 

J

1-1

N t

 

 

 

 

 

 

 

Ut-'l)

 

 

 

 

При

 

 

 

 

 

 

 

 

 

Y: = 4

= const

имеем

q

= >Га .

45.88)

 

 

I ГШ |

Для определения оптимального совместного измерителя про­ цесса Ц1) можно использовать метод Кэдм'ана (I03J, если извест­ но дифференциальное уравнение

А НУ)

 

 

d t

 

 

 

 

 

где f f f )

некоторая непрерывная

функцн

 

 

m(t)

6елы»и шум, д л я . которого

 

 

 

 

 

 

 

 

 

CS.90)

Уравнение

(5.89)

определяет диффузионный

марковский про

цесс; если_к_тому ж е белый шум — гауссовып,

то процесс

будет

нормальным

и м а р к о в с к и м .

 

 

 

 

Д л я вннёровского

процесса

первого

порядка

с корреляционной

пункцией

вида (5.82)

уравнение

(5.89)

выглядит

наиболее

просто;

Учитывая многоканальность наблюдения, введем векторный rt-мерный процесс:

1 \ Ъ

= { W , 0 , . . - , o \ •

1542)

Запишем для него дифференциально - матричное

уравнение

d t

 

 

 

В данном примере (5.91)

 

 

 

f C t ) - 0 , < i r \ C t ) f r \ ' f ^ O > = Q C - t ) S ( > t - ' c ) - B l - § ( , t - t ) ,

где / — матрица, и м е ю щ а я

элементами

нули,

за исключением

единицы в левом диагональном углу.

 

 

Ансамбль н а б л ю д а е м ы х

реализаций

в соответствии с уравне ­

нием (5.92) перепишем в виде

 

 

С5.95)

где

' < , 0 , . . • , 0 . "

 

№ = * H . =

о , . • • , о

 

 

* , о ,

- •

(,5.96)

* n(t)

— белый шум.

о

 

 

 

 

Используя интегральное

уравнение

Винера — Хопфа, вид фи­

зически

реализуемого

оператора

фильтрации, а т а к ж е уравнение

130

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