Файл: Тарушкина Л.Т. Статистическая оценка параметров управляемых систем с помощью ЦВМ.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 Когда массив данных достигает требуемого объема, он поступает на статисти ческую обработку, в результате которой к моменту времени,
равному Т1г |
уточняется значение параметра 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