Файл: Медич Д. Статистически оптимальные линейные оценки и управление.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 11.04.2024
Просмотров: 246
Скачиваний: 1
Показать, что оптимальная оценка х = |
А°у |
вектора х |
имеет |
||||||
вид: |
|
|
|
|
|
|
|
|
|
|
х = |
Е(ху')[Е(уу')]-' |
|
|
у. |
|
|
||
Иными словами, показать, |
что |
|
|
|
|
|
|
||
А° = Е{ху')[Е{уу')]-і |
|
= |
|
РхуР~\у |
|
||||
в предположении, что матрица Руу |
несингулярна. |
|
|
||||||
Показать также, |
что |
корреляционная |
матрица ошибки |
оценки |
|||||
X = X — X имеет вид |
|
|
|
|
|
|
|
|
|
р |
р |
__ р |
р— 1 |
р |
|
|
|||
* |
— ~ — ' X X |
' |
х У г |
ци |
1 |
УХ' |
|
|
|
|
XX |
|
|
|
а |
а |
|
|
|
3-14. Какой физический смысл имеют результаты решения задач 3-11 и 3-13, если функция плотности распределения f(x, у) —гауссов- ская?
Г л а в а ч е т в е р т а я
ЭЛЕМЕНТЫ ТЕОРИИ СЛУЧАЙНЫХ ПРОЦЕССОВ
И ПОСТРОЕНИЕ МОДЕЛЕЙ СИСТЕМ
Целью настоящей главы является построение мо делей систем, для которых впоследствии будут сформу лированы и решены задачи оценки и управления. Основ ной метод главы заключается в объединении результа тов гл. 2 и 3. Здесь используется вероятностное описа ние векторов начальных условий, возмущений системы и ошибок измерения в дискретной и непрерывной моде лях линейных систем из гл. 2. Такое описание можно было ввести непосредственно в начале главы. Однако предварительное изложение некоторых идей теории слу чайных процессов способствует лучшему пониманию ко нечных результатов.
Теория случайных или стохастических процессов изу
чает явления, управляемые вероятностными |
законами, |
в которых время или какая-либо другая |
величина |
является параметром. Строгое и законченное |
изложение |
теории выходит за рамки данной книги, поэтому здесь формально рассмотрены только аспекты теории, сущест венные для наших целей. Читатель, интересующийся более подробным изложением теории случайных процес сов примерно на том же уровне, на каком она рассма тривается в этой главе, найдет более полезными учебни ки Парзена [Л. 4-1] и Папулиса [Л. 4-2]. Более тради-
120
ционное с инженерной точки зрения изложение |
теории |
|
можно найти, например, в |
книгах Лэннинга и |
Бэттина |
[Л. 4-3], Давенпорта и Рута |
[Л. 4-4], Ли [Л. 4-5], |
Джейм |
са [Л. 4-6] и др. |
|
|
После изложения в § 4-1 теории случайных процес сов в § 4-2 и 4-3 вводятся две модели систем. Эти моде ли, так называемые гауссовская марковская последова тельность для дискретных линейных систем и гауссовский марковский процесс для непрерывных линейных систем, с успехом используются в инженерной практике моделирования систем, подверженных влиянию случай ных начальных условий, возмущений и ошибок измере ния. Причина заключается в том, что эти модели оказа лись достаточно простыми для аналитических и вычис лительных целей и в то же время настолько содержа тельными, что с их помощью удалось получить сущест венные результаты.
4-1. ЭЛЕМЕНТЫ ТЕОРИИ СЛУЧАЙНЫХ ПРОЦЕССОВ
Случайный процесс
Рассмотрим явления, развитие которых во времени подчиняется вероятностным законам, например движе ние заряженной частицы в случайно флуктуирующем магнитном поле или отклонение межпланетного космиче ского аппарата от расчетной траектории в результате случайных ошибок в его двигательной и управляющей системах. В обоих указанных случаях шестимерный век тор состояния включает три координаты положения и три координаты скорости в соответственно выбранной системе координат. Множество вариантов развития во времени любого такого явления называется случайным процессом. Короче говоря, случайный процесс представ ляет собой множество случайных векторов, упорядочен ных по времени. Для большей строгости введем следую щее определение.
Случайным |
процессом |
называется |
семейство |
случай |
||
ных векторов |
{x(t), |
tŒ.1), |
зависящих |
от параметра |
t, все |
|
значения |
которого |
принадлежат некоторому множеству |
||||
индексов |
I. |
|
|
|
|
|
Определение не исключает случая, когда х является скаляром. Однако в общем случае будет рассматривать ся векторный случайный процесс.
121
Множество индексов / может быть абстрактным. Однако выше уже говорилось о времени как о параме тре случайного процесса. В дальнейшем будут исполь зоваться только два множества индексов времени. Пер вое является множеством дискретных моментов времени l={tk: k = 0, 1, . . . } , где th<tk+u и моменты времени не обязательно равно удалены друг от друга. Для упро
щения обозначений заменим tu на k. Тогда |
множество |
|||||||
индексов |
примет |
вид / = {6:& = 0, 1, . . . } . |
В |
этом |
случае |
|||
процесс |
называется случайным |
процессом |
с |
дискретным |
||||
временем. |
|
|
|
|
|
|
|
|
Другое используемое здесь множество индексов вре |
||||||||
мени |
включает |
интервалы |
на |
оси времени, |
такие как |
|||
/ = {/ : 0 ^ / ^ : Т) или I = {t : t^to). |
Тогда |
процесс |
назы |
|||||
вается |
случайным |
процессом |
с непрерывным |
временем. |
Понятие «семейство случайных векторов» в определе нии случайного процесса особенно важно. Оно означает, что случайный процесс состоит из множества или из ансамбля случайных векторов, определенных на множе стве индексов. Ансамбль может включать в себя как счетное, так и несчетное множество элементов. Для на ших целей представляет интерес только последняя си туация.
Комбинация понятий времени и ансамбля подразу мевает, что случайный процесс фактически является
функцией двух переменных. Это можно записать с по |
|||||||||
мощью |
обозначения |
[x(t, |
ю), IŒI, |
C O Œ Q } , |
где Q |
— выбо |
|||
рочное |
пространство |
рассматриваемого |
эксперимента. |
||||||
При постоянном значении t x(t, |
•) |
является |
вектор- |
||||||
функцией, определенной на выборочном пространстве Q, |
|||||||||
т. е. случайным вектором. Если |
же |
ш постоянно, |
то |
||||||
x(-, со)—вектор-функция времени, |
представляющая |
со |
|||||||
бой одну из |
возможных |
реализаций |
или |
выборочных |
|||||
функций |
случайного процесса. |
|
|
|
|
|
|||
Если выборочное пространство Й дискретно, то часто |
|||||||||
вместо |
слова |
процесс |
используется |
термин |
цепь. |
Поэто |
му можно говорить о случайной цепи с дискретным вре менем и о случайной цепи с непрерывным временем. Подчеркнем, что о цепи говорят, если ансамбль случай ного процесса содержит счетное множество элементов. Поскольку случайные цепи в дальнейшем не встречают ся, их теория здесь не рассматривается.
Теперь приведем несколько примеров случайных про цессов.
122
Пример 4-1. Рассмотрим скалярный случайный процесс [x(t), t^3?0}, где x(t) = Vo sin t и V0—непрерывная случайная величина, рав номерно распределенная на интервале [—1, lj, т. е.
|
|
|
|
|
|
|
|
|
I |
0 |
в остальных точках. |
|
|
|
||||
|
Ансамбль |
здесь |
состоит |
из |
несчетного |
множества |
элементов, |
|||||||||||
a |
{x(t), |
t^tts], |
очевидно, |
является случайным |
процессом |
с |
непре |
|||||||||||
рывным |
временем. |
Ансамбль |
|
,2 |
|
|
|
|
|
|||||||||
процесса |
можно |
представить |
+ 7 |
|
|
|
|
|
||||||||||
в виде рис. 4-1. |
|
|
постоянно |
|
|
|
|
|
t |
|||||||||
|
Для |
некоторого |
О |
|
I |
р Ъ ^ - ' |
i |
|
||||||||||
го значения |
? = |
/ i ^ 0 , |
x(tt) |
бу |
|
3sz |
|
|||||||||||
-1 |
|
|
||||||||||||||||
дет, |
очевидно, |
случайной |
вели |
•л jI |
|
|
|
|||||||||||
чиной. Любая |
синусоида явля |
+ J |
|
|
t |
|||||||||||||
ется |
выборочной |
функцией |
про |
- ч |
I |
|
|
i |
||||||||||
цесса. |
|
|
|
|
|
|
|
про |
О |
\ |
I |
/ |
|
|
|
|||
|
Если дискретизовать |
-1 |
- |
i |
|
|
|
|
||||||||||
цесс, |
рассматривая |
только |
мно |
|
|
|
|
|||||||||||
жество |
моментов |
|
времени |
|
|
|
|
|
|
|
||||||||
{tk |
:/„ = Ая/4, |
*=0, |
1, . . . } , то |
|
|
|
|
|
|
|
||||||||
в |
результате получится |
случай |
|
|
|
|
|
|
|
|||||||||
ный процесс с дискретным |
вре |
|
|
|
|
|
|
|
||||||||||
менем. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
Пример 4-2. В качестве вто |
|
|
|
|
|
|
|
||||||||||
рого |
примера |
рассмотрим |
ска |
|
|
|
|
|
|
|
||||||||
лярный |
случайный |
процесс |
с |
|
|
|
|
|
|
|
||||||||
дискретным |
временем |
{у |
(к), |
|
|
|
|
|
|
|
||||||||
k=0, |
1, |
. . . } , где для |
каждого |
|
|
|
|
|
|
|
||||||||
k |
предполагается, |
что |
у(к) — |
Рис. 4-1. Представление случайно |
||||||||||||||
гауссовская |
случайная |
величи |
||||||||||||||||
на с нулевым средним. Возмож |
го |
процесса |
в |
виде ансамбля. |
|
|||||||||||||
ную |
зависимость |
дисперсии |
у (k) |
от |
|
|
|
|
|
|
||||||||
k учтем, записывая плотность рас |
|
|
|
|
|
|
||||||||||||
пределения |
y{k) |
в |
виде |
|
|
|
|
|
|
|
|
|
1
Î(У, |
к): |
(k) |
Ѵ2л |
|
|
|
|
||
где |
a2(k) |
=Е[у*(к)]. |
Предполагается, |
|
что |
значение y(k) |
для заданного k |
статистически не зависит от его зна чения в любой другой точке множе ства индексов.
Две типичные выборочные функ ции процесса показаны на рис. 4-2.
Этот процесс называется чисто случайной гауссовской последова тельностью. Прилагательное гауссовский используется с целью подчерк нуть тот факт, что при известном k распределение значений y{k) по ан самблю гауссовское, а определение
?
3 "
о
Рис. 4-2. Выборочные функ ции дискретного случайного процесса.
123
«чисто случайный» указывает на статистическую независимость во времени элементов во всех выборочных функциях. Термин «последо вательность» возникает в силу того, что любая выборочная функция представляет собой просто последовательность случайных перемен ных, определенных в дискретные моменты времени.
В дальнейшем процесс этого типа подробно исследуется при опи сании возмущений системы и ошибок измерения в дискретных линей ных системах.
Описание случайного процесса
Из приведенных рассуждений ясно, что описание случайного процесса включает истинное (или вообра жаемое) изображение ансамбля и вероятностный закон, управляющий ансамблем. Здесь исследуются возможные способы описания такого вероятностного закона.
Рассмотрим случайный процесс с дискретным време нем {x(k), kŒ\}, где / = {&:&=1, ..., N}, a N— некото рое постоянное положительное целое число. В этом слу чае естественным описанием процесса будет совместная
функция |
распределения |
вероятностей |
N |
случайных |
|||||||
/г-векторов |
х(\), |
|
|
x(N), |
определенная |
для любых |
|||||
п-векторов |
хх, ..., |
|
xN |
|
|
|
|
|
|
|
|
|
F(x\ |
|
х»)=Р{х(\)<^х\ |
|
x(N)^x»l |
(4-1) |
|||||
|
Эквивалентное описание можно ввести с помощью |
||||||||||
совместной |
плотности |
распределения |
вероятностей |
|
|||||||
|
f(xl |
|
xN) |
= |
• - j |
г-^-я |
|
ж ' |
(4 -2) |
||
|
, |
ѵ |
|
; |
|
дх\...дх1п...дх^ |
...dxNn |
ѵ |
' |
||
где |
х'{і=\, |
|
п\ |
/ = 1 |
|
N) является |
і-й компонен |
||||
той |
n-вектора xj |
в |
предположении, |
что |
указанные |
nN |
частных производных существуют, или через совместную характеристическую функцию
|
|
|
|
|
N |
x' |
(t) |
si |
(4-3) |
|
9(s\ .... |
s")-- |
|
ехр / £ |
|||||
для любых и-векторов s'(t =l , ..., |
N), |
где |
/=]/"—1- |
|
|||||
Если множество индексов является дискретным, но |
|||||||||
бесконечным, т. е. I = {k:k=l, |
2 . . . } , |
или же непрерыв |
|||||||
ным, ограниченным |
или |
неограниченным, |
т. е. / = {^:0^ |
||||||
ssZt^zT} |
ИЛИ I={t |
: t^sO}, то |
подобное |
|
описание |
не бу |
|||
дет столь естественным. Однако с практической |
точки |
||||||||
зрения |
представляется |
разумным, |
что |
случайный |
про- |
124