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

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

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

 

А = Л - * - < V к , ,

 

 

>

 

 

Л . к. = А к +• X

к с - Ь 1(

.

'

(,2 . 90 )

 

''Fge_ О.'— матрица,

обратная

//..,.

 

 

 

 

 

Физически реализуемый алгоритм дискретной фильтрации бу-

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

(2.90) лишь

верхним

пределом

сум­

мирования. Вместо п здесь будет к. Дисперсии

ошибок оценок

па­

раметров совпадают с диагональными элементами матрицы Q

'Рассмотрим задачу амплкгуднон демодуляции, когда

 

Функциональные

производные

имеют

вид

 

 

 

'/Зри

помехе т и п а ' б е л о г о ' ш у м а

получим

 

 

 

k

N0

J

 

 

 

 

 

 

 

 

 

 

'

ik-OT.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

No

 

 

 

2 ke

- N v

 

* c '

-

^

• T„ -

 

 

 

Если ж е помеха

коррелированная, то оценка

величины Лу. зави­

с и т № з'НачеНИй

ла тех интервалах,

на

которых

 

 

П'р'й э1*Штиатрица

/ / 2

не диагональная

с

элементами

 

Я

= ( а г д - Г

'

(

coo t h

 

COS us0 t

d t

d<t.

(,2.46)

2

k t

'

 

 

 

 

 

 

 

 

Алгоритм

амплитудной

демодуляции

имеет

вид

выражения

. (2.90),

где

 

 

 

 

 

 

 

 

 

 

40

Если допустимо полагать, что на соседних интервалах времени наблюдения корреляция помех отсутствует, то матрица Я? будет диагональной с элементами

k e

" lk-i)T0 (k-<Ko '

а структура измерителя определяется выражением:

* ( U 0

[ J C Q S U - ^ - A " \ t i ^ ^ ^ ) - ^ ^ H > k ) c o s , w T ] d T , d T ( г . 9 9 )

 

\Ч-От„

 

 

 

 

 

 

 

 

Вариант

измерителя

следящего

типа

получим из

выражения

(2.99)

заменой

5 l k

на

\

.

 

»

 

 

Приравнивая

нулю

выражение

для

W(

найдем

алгоритм

оценки

А ь

в явном

виде:

 

 

 

^

 

 

\ = [~Г~

\ \

A ' \ t , * № « e t t - * ) d t d « c ]

'

 

 

 

( , k - t ) T 0

 

 

 

 

 

kTB

 

-\

 

 

 

 

0

 

 

 

* J U o

|j

cos

co^t

A ( 4 t , T r ) c j ( f ) d t - d < c - 1. .

(2 .100)

При

фазовой

модуляции,

когда

 

получим

. а

 

Л

- с

 

. ( u e c o s ^ e t ^ ' U t ^ ] d t .

°

 

«,2.105)

Если на отрезке времени 7'п

процесс

\ ( ( ) можно

полагать no­

 

 

 

i l

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

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

Функционал

(2.2)

имеет вид:

 

 

0 ч

Возвращаясь

к соотношениям (2.7), (2.8), найдем выражения

для функциональных

производных:

х \ А Ч з д { г ^ г ) - м ^ о * [ и в * + $ A \ { ) d i f ] ^ d S - d t

(.2.106)

H { t , t , A U ) , k ^ = f e \\ si n [ u s * j A t r ) d x ] *

t 't

0

*

A ( ^ l ( ) 5 L n [ u J + j

 

A U ) d E ]

d ? d j k .

(2.Ю7)

 

 

о

 

 

 

Положим для простоты,

чк1 аддпiинпын ш\м белый. При этом

вместо

cooiношений (2.10(>)

 

и (2.1ч7)

будем имен,

выражения:

^ , A № , i j t t ) } = ~ J

^ i m [ ^ J A \ x ) d x ) ^ ( c 3 ) -

 

0

t

,f

о

 

J sin [ щ в ? + j Л (X) d x ] d E

при t > * ,

.2 T

~~

j Sin [w o e + \ A I X ) d x ] dG"

п о и

-12

Алгоритм частотной демодуляции имеет вид в ы р а ж е н и я (2.90),

Н

 

 

 

( [ T - » n a x ( t , * ) l d t d * .

 

В частности, если T=2T0

(fc,/=l,2)

и

о ц е н и в а ю т с я - з н а ч е н и я

частоты j \

, и A 2

на

подынтервалах,

то

матрица

Нг имеет в и д :

 

 

 

2

- И

3

 

 

 

 

 

 

 

г* о

Т 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N .

2

 

3

 

 

 

 

 

 

 

 

 

 

Матрица

дисперсий

ошибок

будет

 

 

 

 

 

 

 

 

 

 

 

 

(.2. И

5)

 

 

 

 

 

1

г

 

 

 

Области

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

при

нахождений

изображены на

рис. 2.2. Величина

ц

является объемом

части

пирамиды,

пока-

Р и с . 2.2. Области имтег-рироиашш' при определении h.

• 'к

'43'

 

ланпоп на р и с 2.3 Jnpii ' 2 _ ~L*>

Н о

'

 

 

II

 

P l J i " .

'•' I. I ( . ' O M L ' T p i l ' l l Y K U H

till l"l'pMPl'l ПЦ11Я НСЛПЧ11НЫ

" j j f c f •

an.

ЛЛ1ШЧ1ТМЫ

ОПТИМАЛЬНОЙ о ц е н к и

 

 

ГАУССОВЫХ П Р О Ц Е С С О В

 

Л Tropinм оптимальной

оценки

гпуееонот процесса

получим,

приравнивая пулю первую функциональную производную от фучк

цппиала

 

 

 

 

тт

 

 

 

 

 

г

г

 

 

 

 

 

 

 

 

 

 

о

0

 

 

 

 

 

 

 

 

 

 

Р а з л а г а я

второе слагаемое

п

рял (2.'>) и

 

оставляя с л а ш е м ы е

ло

квадратичного

включительно,

получим:

 

 

 

 

ГJ[[Т

^"^4

о

 

 

 

 

 

 

 

 

 

a

 

>

,

 

 

*" I "

 

A

 

K - { t . ^ A v , t ) ^ ^ [ A ^ - M - t ) ] d t d r -

где

функции

h . , { t ,Л.1Д)^ и

n . ^ t , t , X U ) , M t ) \ совпадаю: г пер­

выми двумя

функциональными

производными

 

фупкн.иоиала(2.2).

 

Введем,

как и

прежде,

функцию

 

 

 

и

обозначим

 

 

 

 

 

 

 

 

 

•11

Используя выражение (2.115), (2.117), запишем уравнение правдоподобия:

т

**•*

Группируя второе и третье слагаемые

и у р а п ш и п , . (2.118) и

используя

функцию L(l,x); определенную

уравнением

получим

два варианта измерителя:

т

 

т

 

О

Ci

о

[ A ^ J - ^ W l d e - ^ d t

(.2.120

Физически реализуемы/! алгоритм оценки процесса получим из выражений (2.120) или (2.121) изменением пределов н т е г р н р о - иання

 

 

^

 

л,

 

 

 

^2

icl*^}

 

+

W * V t , < 5 - ) [ % ) - A ^ ] d e - \ d r .

 

 

 

Если

разложение

(2.115) осуществляется

в окрестности

средне­

го значения процессов, то оба варианта алгоритмов

д а ю т

 

 

 

- 4

t

) + f

U t / O - k \ * , Ч * ) \

d ^

 

 

{глт

 

 

 

 

о

 

 

 

 

 

Имеется

и еще

одна

возможность

для получения

разновидности

алгоритма оценки процесса. Если

в т о р о е . с л а г а е м о е

выражения

(2.1.14) разлагать в точке максимально правдоподобного

значения,

полагая

X(i)

Л (А) > то вместо

вы ра жени ft (2.120)

и

(2.121)

получим:

 

 

 

 

 

 

 

 

•15

Тт

ОQ

H*

T

T

t

A 0:) = <Kt) + \ Lv t,*) dt ^ H фр)

[I (S) - 0 (tf)] dor.

 

 

 

 

_o

 

 

0

 

 

 

 

 

 

 

В поля функции

 

 

 

1 1

"^П",®*)

 

 

 

 

•пишем

соотношение

(2.124) u

нидс

 

 

 

 

 

Л U) ~ A (X)

v ^ X 4 tч

-

 

dt ,

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

A ( t ) - 4 t ) * ] N J t f , * ) [ A ( * ) - H O ] cl t -

 

 

 

 

Алгоритм (4.121)

н а з ы в а ю т следящим,

если

п о л о ж и т ь

> (г) =

-•A* (<'*-) • " а

(>eiioiU"

Же

алгоритма

(4,120)

п р о щ е

всего

построить

т с п а Ц н о П п Ь ш

HiMepK тел ь, п о л о ж и в

Д ^ • ! {П , , (Н-)

 

 

 

 

 

 

 

 

\ ^ ) { \ [ 0 ' , v . , l o ]

*'

 

 

 

 

 

 

т

 

 

 

 

 

 

,

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

 

Л.Ч10|)цтм

123)

Наиболее прост

1-го называют

неслодйШНМ,

псскольку нрпипя части HP зависит oi

предшествующих

оценок.

Алгоритмы

(2.12-1)

т а к ж е

пеелсд и том»

'ii.n.'i, указывающие,

что оп­

тимальна*,!

оценка

может

быть подучена

л и н е й н о й ,

фильтрацией

оценок • максимального

правдоподобия.

 

 

 

 

 

Д л я

оценки

точности' измерения

процесса вычислим

функцию

 

 

j M t , T 0 - < , [ A \ t ) - A ^ y L A \ - i ! V - A ^ > .

 

 

Ксли

положить

 

^(т.)

-лч О .

то ич вь1раженпя

(2Л20)

непо-

средстьеппо

получим:

 

 

 

 

 

 

 

 

 

°т т

=

\ \ U ( . t , o - ) L az.pH

d a d j +

 

о о

 

т т

тт

 

так что •л i

 

 

 

e-J, (to

=

L U . t ) .

 

 

 

 

(2.12ч)

Такие ж е

потенциальные

возможности

имеет и

измеритель, постро

енный

на

основе aj

орнтма

(2.121).

Действительно,

 

 

т т

 

 

с о

 

т т .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D Q

Т

 

 

 

О Q

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

-

2 JJ

 

\\{J5,x)\f{K,*)

 

d e d x .

 

 

12.128)

Здесь

учтено,

что

 

 

 

 

 

 

 

 

 

Заметим, что из

уравнений (2.117)

и

( 2 . i H )

следует,

что

 

J U (Д,<0

H( . t, з-4)dLt-t J Lvt , * ) V f

\

* d t

= £ < 2 0 . V i . l i O )

 

 

 

 

 

• u

 

 

 

 

 

 

 

Умножая

это

равенство

на Yl (.<?,•()

и интегрируя

по о,

получим;

j J U ^ , t ) H C t , t J ) W ^ ^ ; d t d e T = Y ^ U )

- L ( . t , ) P

ч2ДМ)

о <з

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя

соотношение

(2.13!)

в

выражение

(2.128), имеем

 

 

 

 

т т

 

 

 

 

 

 

 

 

j b ^ ^ ^ U ^ + S J ^ ^ ^ L ^ ^ ^ H ^ ^ d e d ^

+

 

 

 

.

0 0

 

 

 

 

 

 

 

 

что и требовалось. Обозначив

т

HJ

о

 

 

равенства (2.131) получим интегральное уравнение

относитель­

н о

функции /. ( Л т ) :

 

 

 

т

 

 

 

J L(4 t,-c)a3-(,'r;^)d'r: + U( s t,Q) = ^ \Х,<3).

(2.<ЗЬ)

 

Рассмотрим неследяшип

измеритель, реализующий

алгоритм

(2.124). Учитывая, что

'

 

17

ч

[ х оо

А и % а л * ^ - out i]> ~ q u / о ;

 

-..Livt)-AU)IvXtVAa, )]> - о ,

(.8.1310

вй чем

н» п е р в о ю

com ношения (2.124)

 

 

 

г т

 

 

 

0 Q

 

 

T r

_

 

 

0 Q

 

 

 

 

-1

 

(J 0

H i ранено гни (2 130) умножением обеих е ю частей на Q(<3,{-)

иинпч рнропдппе.м по о получим:

ТТ

 

l\ИД

*)YJ"taWW/OD G "d * *

$

v.t,t)-L^Шйб)

 

•SO

 

п|.|рилоч1пе

 

 

.

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

.

 

 

 

По triaiviiiH

 

(2.136)

(2.135)» получим,

как и

прок,ip.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

К л о м у

же

результату

ирпвои>|

oiunni ичиг « ирсоОрнзоп.шпн,

выполненные

со

вторым

• ooiпопнч'Нем (1? 121),

Таким

образом,

потенциальная чочпоен,

оценки

процесса

и'Морнтелимн

(2.120),

(2 121) и (2.124) одинакова,

ч ю и

с\ч<лики*

о ж и д а т ь ,

поскольку

разложение функционала отношения нргл.'юиочобпа

осуществля­

лось в

«точке»,

близкой

к

пс'1 иппой

реяли о'цин

}Л),

л и б о

совпа­

дающей

с ней.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ч ю

касается

н ш е р и т е л я

(2.12.3),

то он

неоптималеп,

причем

^ t / t ) = w a ^ 5 J u i t > < j ^ ^

 

 

 

 

.

 

^ - 138 )

O I M C I I I M

 

 

<з о

 

 

 

 

 

 

 

 

 

 

 

 

одну

ночможноси,

упрощения

т м о р и т е д н

нормально­

го процесса. Пусть аддитивный шум белый,

а с т а л

связан

с про­

цессом

оператором

Нсмыцкого

(временной

 

аргумент

сигнала и

параметра совпадают) . При этом

 

 

 

 

 

 

 

 

U\txA\t)Alx)}--

 

 

 

 

 

 

 

 

 

.J-8v t-t) •

 

 

 

 

 

 

N 0

 

 

С1Л

 

 

 

a. A

 

 

l a . l i q )

Подставляя

выражение

 

(2.13°)

в

уравнение

(2.121)

при Х[Х)~

=- Л ' I)> по.г.-чим

 

 

 

 

 

 

 

 

 

 

 

 

 

с

> 1 L

J

(2.140 1

В соответствии'с уравнением (2.132),

Если сигнал высокочастотный и узкополоснын, то выражение

(2.141) содержит быстрооецплллрующее слагаемое,

которое Прн

последующем

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

в выражении

(2.133)

дает резуль­

тат,

.близкий

к

нулю .

 

 

 

В

частности,

если

 

 

 

то вместо выражении (2. МО)

и (2,133) будем

И Меть

 

 

 

 

т

 

^

 

Такие же уримненпя характеризуют измеритель процесса

при фа-

зоной

демодуляции

(2.101),

Различие

измерителен^

обусловлено

л и uiii

характером

формирования

функции

"k4

{t, ^WTV 1

Рассмотрим случпй

частотной

демодуляции,

когда

елравеллн-

оо соотношение (2.109). Используя

алгоритм

типа

(2.121), по*

лучим::

 

+ f

L ( Д Л ) { Ц L*. A U

 

+ ™ (li

 

A \ t ) « * <Д)

) l

Физически реализуемый алгоритм имеет вид:

- 4<?)]dG'} clT ,

С2ЛЦ5)

'19

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