
книги из ГПНТБ / Хомяков Э.Н. Вопросы статистической теории оптимальных измерительных систем. Основание для расчета и проектирования
.pdfма имеем |
теперь |
уравнение |
Вольтерра. |
|
на^-ос.) по |
||||||
Заменяя |
нижний |
предел |
в соотношении (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.М) |
||
|
|
|
|
|
|
||
Таким |
образом, |
|
|
|
|
|
|
|
|
|
|
2б |
|
|
|
Если компоненты |
процесса |
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 . = |
1» |
о , . • • , о |
||
|
|
* , о , |
- • |
(,5.96) |
|
* n(t) |
— белый шум. |
о |
|||
|
|
|
|
||
Используя интегральное |
уравнение |
Винера — Хопфа, вид фи |
|||
зически |
реализуемого |
оператора |
фильтрации, а т а к ж е уравнение |
130