
книги из ГПНТБ / Хомяков Э.Н. Вопросы статистической теории оптимальных измерительных систем. Основание для расчета и проектирования
.pdf
|
А = Л - * - < 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