Файл: Тарушкина Л.Т. Статистическая оценка параметров управляемых систем с помощью ЦВМ.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 02.07.2024

Просмотров: 108

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Здесь Х п р (0 — программное перемещение автоматической си­ стемы; X (/) — действительное перемещение.

Предположим, что процесс У (i) имеет нормальный закон распределения. В процессе перемещения автоматической системы в ЦВМ вычисляется вероятность

P\\Y(t + T)\^z},

где t— текущее время; т — величина упреждения; z— величина заданного допуска.

ио

x(t)

 

ЦВМ

Рис. 1. Управляющая ЦВМ по . выработке закона управления:

ИО — исполнительный орган

Коррекция траектории происходит дискретно в те моменты времени, начиная с которых указанная вероятность меньше заданной величины. Предположим, что

 

 

MY(f)

= m, Ку

(t,

t) =

cF (t),

 

где m,

с—неизвестные

постоянные;

F (t) — известная

функция.

Тогда

моменты

распределения процесса

У (f)

имеют

известную

 

 

4

 

 

 

 

 

 

Формиро­

Результат

Вычисление

Алгоритм

вание

запоми­

вероятности

коррекции

массивов

нается

 

Р

 

траектории

 

 

 

 

Алгоритм

Изменений

Алгоритм

 

 

обработки

нет

 

управления

 

 

Изменение

 

 

 

 

 

 

 

параметра

 

 

 

 

 

 

 

•я

 

 

 

 

 

 

 

Рис. 2. Решение алгоритмов

управления,

коррекции, статистической обработки

 

 

 

с помощью

ЦВМ

 

 

 

12


аналитическую структуру

и q = (т., с) — неизвестный параметр

системы управления.

 

 

В начальный момент при t — 0 зададим

априорное'значение

параметра q = q 0 . После

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

t — О сигнал Y (t)

измеряется и его значения записываются в память ЦВМ, обра­ зуя массив данных Rx — (У (^), . . ., Y Когда массив данных достигает требуемого объема, он поступает на статисти­ ческую обработку, в результате которой к моменту времени,

равному Т

уточняется значение параметра q. При выборе мо­

ментов коррекции на промежутке времени

[0, Тг] используется

априорное

значение параметра

q = q 0 ) на

промежутке времени

[7\, Т]

уточненное значение

параметра

q.

Блок-схеме алгоритмов системы управления с учетом после­ довательного уточнения значений параметра q на отдельных промежутках времени представлена на рис. 2.

2.КЛАССЫ ФУНКЦИЙ И АНАЛИТИЧЕСКАЯ СТРУКТУРА МОМЕНТОВ РАСПРЕДЕЛЕНИЯ

Определение класса. Класс случайных функций определяется вероятностными свойствами функций. Случайные функции отно­ сятся к одному и тому же классу, если они имеют общие вероят­ ностные свойства, определяемые первыми двумя моментами рас­ пределения.

При статистической обработке данных класс функций можно определить следующим образом.

1.Пусть в системе управления измеряется та координата, класс которой следует определить. В этом случае класс функций определяется по результатам обработки реализации измеряемой координаты.

2.Допустим, что в системе управления измеряется некоторая координата. Определяется класс функций неизмеряемой коорди­ наты, но связанный с измеряемой координатой. В этом случае задается математическая модель связи измеряемой и неизмеряе­ мой координат. Если модель такова, что по измеряемой коорди­ нате одназначно определяются моменты распределения неизме­ ряемой координаты, тогда класс функций измеряемой координаты однозначно определяет класс функций неизмеряемой координаты.

Впротивном случае в системе управления для неизмеряемой координаты можно указать лишь множество допустимых классов.

Отметим, что класс случайных функций в процессе работы системы управления может • изменяться, тогда для измеряемых координат задается конечное число допустимых классов функций. \ В процессе управления требуется по результатам обработки изме­ ряемых координат среди допустимого класса функций выбрать определенный, указав при этом значения моментов распределения данного класса.

13


Перейдем к рассмотрению отдельных классов случайных функ­ ций. Особое внимание будем уделять определению аналитической структуры моментов распределения, характеризующих данный класс.

Процессы непрерывные в среднем. Случайный процесс X (t) непрерывен в среднем, если для любого, значения tx 6 Т

l i m M [ X ( 0 — Х ( ^ ) ] 2

= 0, t£T.

<->/,

 

Пусть X (t) непрерывен в среднем и его математическое ожи­ дание равно нулю. Так как корреляционная функция является положительно определенным ядром, то по теореме Мерсера

W l t к)

(1.8)

\'=i 4

причем ряд сходится равномерно. Здесь cov и <pv (/) — соответ­ ственно собственные значения и собственные функции интеграль­ ного уравнения

т

q>v('i) = q>Jtf,('i. 4 ) ф ( 4 ) ^ 2 -

(1.9)

о

В силу того, что корреляционное ядро является положительно определенным, имеем cov > 0. Без ограничения общности можно считать, что система функций cpv (t), входящая в (1.8), является ортогональной.

Используя теорему Карунена [23, 24], получаем представ­ ление вида

со

Я Ю - Ё * ^ *

( 1 Л 0 )

 

v = i

\

 

 

где х\, ортонормированные

случайные

величины:

 

Mxv = 0,

Мх^

= 0

(v Ф (.1),

 

 

Mxl=l.

 

 

(1.11)

Разложение (1.10) является каноническим разложением. Раз­ личные представления канонических разложений даны в [19].

Ограничимся первыми п членами ряда в разложении (1.10). Дисперсия ошибки такого приближения за время Т равна в сред­ нем

т

п

т

^ (1.12)

= -L\Kx{t,

t)dt- £

-±-\<&{t)dt.

О

V = l

о

 

14


При известных значениях корреляционной функции Kx(ti, / 2 ) число п выбираем таким, чтобы величина ошибки (1.12) не пре­ восходила допустимого значения. В этом случае аналитическая

структура

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

функции определяется отрезком

ряда

 

 

 

 

п

 

 

 

 

КЛк,

к) =

2 - ^«М*х)<М* 8 ),

 

(J.13)

где cov — неизвестные

собственные

значения

корреляционного

ядра; cpv

(t) — известная

система

функций.

 

 

Если математическое ожидание tnx (f) не равно нулю и при­

надлежит

классу функций,

суммируемых с

квадратом,

тогда

 

 

 

 

со

 

 

 

 

 

mx{t)

= 2

bjv(t),

 

 

v = l

где Д, (t) — ортогональная система функций; bv — коэффициенты Фурье функции mx(t).

Вэтом случае аналитическая структура математического

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

 

 

 

т Л 0 = Е Ш * ) ,

(1-14)

где

fv (t) — известные

ортогональные

функции; Ьч

неизвест­

ные

параметры.

 

 

 

 

 

Ошибка такой аппроксимации равна за время Т

 

 

Т Г

k

Т2

т

k

 

 

mx (t) -

2

bvfv (t)

dt = J ml (t) Л J b\.

 

 

 

 

 

 

v = l

 

При известных значениях функции mx (t) число k выбирем таким, чтобы величина не превосходила допустимого значения.

Таким образом, для процесса непрерывного в среднем неиз­ вестные параметры q x = ъ . . ., bk), q 2 = (сох , . . ., а>„). входят в аналитическую структуру моментов распределения.

Процессы с независимыми приращениями. Рассмотрим возра­ стающую последовательность моментов времени

 

О < tx

<

t2 < . . . < ts <

Т.

(1.15)

Составим

разности

 

 

 

 

Z i

= X f t ) ,

zv =

X ( / v ) - X ( / v _ i )

(v = 2Ts).

(1.16)

Случайный процесс х- (t) называется процессом с независимыми приращениями, если случайные величины взаимно незави­ симы.

Имеем

X(tt)=t*v.

(1-17)

v = l

15


Обозначим через тг (tv) математическое ожидание величин zv, Mzv = тг (tv). Так как величины zv взаимно независимы, то

KAti,

£

M[zv-mz{Q}2.

(1.18)

 

v = l

 

Отсюда следует, что

 

 

 

Кх Vi, h)^Kx

(4,

к) < . . . ^ К х (ts, Q,

(1.19)

т. е. дисперсия процесса X (t) есть неубывающая функция вре­ мени.

Для определенности предположим, что

 

mm{tit

tt) = t,.-

(1.20)

Тогда

 

 

 

Kx(ti,ti)

= Kx(mm(titt,)t

min(tit

t,)) = Kx(ti> it)- (1-21)

Найдем каноническое разложение процесса с независимыми приращениями. Исходя из свойств функции ц (t) = Кх (t, t) определим в гильбертовом пространстве скалярное произведение

 

 

 

 

 

( Ф ( 0 ,

 

fW)=

J c p W f ( T ) d p . ( T )

 

 

и

рассмотрим

полную

ортонормированную

систему

функций

9 V

(^). Обозначим через

£ (^, ^) функцию, равную единице при

t ^

tt

и

нулю

при t

>

 

имеем

 

 

 

 

В

силу

равенства

Парсеваля

 

 

 

 

 

S

(Ф» (О,

.с с . У )

(ФУ (о.

с (*, */» =

(s (*;

*,), с (t,

щ,

тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

CO

 

 

 

' /

 

 

 

 

 

 

 

 

 

£

J Ф»(«) Ф

(s) J cP v (0 ф (0 = J d|i.(0 =

ft,

^)

(0, 0).

v = l

0

 

 

 

 

 

 

 

 

 

 

 

 

Используя

теорему

 

Карунена

[23,

24], получим следующее

каноническое представление

 

 

со

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X(t)

=

mx(t)-mx(0)

+

X(0)

+

£

xv\ ( p v ( s ) 4 Ф ) ,

(1,22)

 

 

 

 

 

 

 

 

 

v = i

о

 

 

где.ху —ортонормированные случайные величины, т. е. вели­ чины, удовлетворяющие условиям (1.11).

Ограничимся первыми п членами ряда в разложении (1.22).

16