Файл: Тарушкина Л.Т. Статистическая оценка параметров управляемых систем с помощью ЦВМ.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 02.07.2024
Просмотров: 116
Скачиваний: 0
где Ak (t)— матрица-строка порядка п + 1» элементы которой являются непрерывными функциями.
Заметим, что динамическая система (IV.26) имеет вид системы (Ш . 2) .
Элементы Akj (t) матрицы Ak (t) под влиянием случайных воз действий могут изменить свои значения. Предполагается, что
функции |
Akj |
(t) имеют аналитическую |
структуру |
вида |
(IV.21), |
||||
причем |
Bkj(t), |
(pkjs |
(I) — |
известные |
функции |
t£2T; |
akjs |
— |
|
неизвестные |
коэффициенты; |
г — конечное число. |
|
|
|
|
|||
Требуется |
среди |
корректирующих |
устройств |
Wk |
(k |
= 1, |
р) |
||
выбрать такое, у которого дисперсия процесса Yk |
(t) |
минимальна |
за все время работы системы управления. Такое устройство об ладает максимальной нечувствительностью к изменениям внеш ней среды, изменениям динамических систем, описывающих ис полнительный орган и работу самого корректирующего устройства.
Определение минимального объема данных по априорным зна чениям моментов. В начальный момент работы системы управле
ния |
при t |
= |
0 зададим априорные значения для элементов ма |
||||
трицы |
АА |
(/). |
Априорные значения определяются теми значе |
||||
ниями, которые заложены при проектировании |
системы |
управ |
|||||
ления. |
|
|
|
|
|
|
|
Для процесса X (t) зададим полосу частот |
[Qlt й а ] , |
в |
кото |
||||
рой |
будем |
|
рассматривать спектральную плотность |
процесса |
|||
X (i), |
t£ |
Т. Используя теорему Котельникова, выберем шаг дискрет |
ности At перехода от непрерывного процесса X (t) к дискретному процессу. Предположим, что при проектировании системы управ ления определен тот минимальными объем статистических данных, соответствующих массиву Rx = (X . . ., X (tm)), который необходим для того, чтобы определить статистические характе ристики процесса X (0 с требуемой точностью, используя, на пример, методы, указанные в п. 6. Минимальный объем массива статистических данных R^ определяется с помощью ЦВМ сле дующим образом. Возьмем каноническое представление (1.37) процесса X (t), выберем число п1 таким, чтобы ошибка е л 1 , оп ределяемая формулой (1.40), не превосходила допустимой вели чины. Рассмотрим разложение (1.37) с конечным числом членов разложения, равным п х ;
X{t) = m t + S (xv cosfflv £ + yv sinсо.Д - v=l
c o v = - ^ - . (IV.27)
1
Будем предполагать, что процесс X (t) имеет нормальный за кон распределения. Тогда случайные величины xv, yv нормально распределены и являются независимыми ортонормированными величинами. Используя датчики случайных чисел, для величин xv, yv определим выборочные значения. Построим одну реали зацию для процесса, имеющего разложение (IV.27). Используя методы, изложенные в п. 6, определим время Tlt при котором
132
массив статистических данных Rx дает оценку для моментов рас пределения процесса X (t) с точностью-не ниже, чем допустимая точность. Выбирая время Тг минимальным, получим минималь ный объем массива статистических данных R^.
Из процессов Yk (t) выберем тот процесс, который имеет мини мальную дисперсию за время работы системы управления, рав ное Тг, т. е. такой процесс Yk (t), для которого
г, |
|
|
|
\ M[Yk(t)-MYk(t)fdt |
|
= |
|
о |
|
|
|
|
г, |
|
|
= rmn |
\M[Yt{t) |
— MYt(i)fdt. |
(IV.28) |
i ' = l, p |
о |
|
|
Процесс Yk (t) определяет априорное значение корректирую щего устройства Wk, которое будем использовать для управле ния системой на отрезке времени [0, Тг].
Уточнение параметров корректирующих устройств на' началь ном промежутке времени. Чтобы уточнить неизвестные пара
метры akjs |
функции |
Akj |
(t), |
входящей |
в корректирующее |
устрой |
||||||||
ство Wk |
(k = |
1, р), |
а также выяснить, правильно ли было выбрано |
|||||||||||
корректирующее устройство |
Wk |
на |
отрезке времени |
|
[0, |
Тг], |
||||||||
образуем |
в |
ЦВМ |
массив |
статистических |
данных |
Ry |
= |
(Y {tx), |
||||||
Y |
(t2) |
|
Y ( У ) , где |
Y (tt) |
= |
( Г х |
(*,), |
. . ., Yp |
(tt)) |
(i |
= |
7 ^ ) ; |
||
tm |
^ Тх. |
|
|
|
|
|
|
|
|
|
|
|
|
Для того чтобы по статистическим данным массива Ry произ вести уточнение априорных значений матрицы АА (t), "необхо димо, чтобы промежуток времени [О, Тх] был меньше, чем про межуток, равный [О, Т]. В противном случае априорные значе ния матрицы Ak (t), t£T не могут быть уточнены с требуемой точностью. Поэтому будем предполагать, что 7\ ==£ Т.
Пусть |
значения |
производных F<*> (t) (k = |
1, п) вычисляются |
|||
в ЦВМ с требуемой точностью. Представим значения |
интегралов |
|||||
|
|
|
т, |
|
|
|
|
|
i |
\Ak(t)LYk(t)dt, |
|
|
|
|
|
1 |
о |
|
|
|
4J-l |
оJ №k{t + |
x)LYk |
(t + т)-mj |
[Af t (t) |
LYk(t)-m,] |
dt |
в виде сумм интегралов, значения которых известны для данной реализации процесса. Yk (t), с неизвестными параметрами akis и их произведениями.
133
Аналогично функционалу (IV.3) для каждого корректирую щего устройства Wk составим функционал
|
|
|
|
|
Г, |
|
|
|
1 2 |
|
|
|
|
|
Vlk |
= |
•±- \ Ak(t)LYk(t)-mx |
|
|
+ |
|
|
|
||
|
|
|
N |
l |
T |
|
|
|
|
|
|
|
+ |
ir |
£ |
T " I [ A a ( * + |
T / ) L n ( |
f |
~ |
m J |
x |
|
|
||
|
|
1 = 1 |
{ |
0 |
|
I |
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
X |
[Af t |
(0 L K (0 — mx] dt — Ce~aIT< |
IJ |
|
(/e = |
1, p). |
(IV.29) |
|||||
Найдем минимальное значение функционала (IV.29) по всем |
||||||||||||
неизвестным параметрам akls, |
входящим в. матрицу Ak |
(t). |
|
|||||||||
Если априорные |
значения |
неизвестных |
параметров |
матрицы |
||||||||
Ak (t) в пределах |
допустимой |
точности, совпадают |
со значениями |
|||||||||
параметров |
ак/5, |
|
дающих минимум функционалу (IV.28) для |
всех |
||||||||
k (k — 1, р), то |
корректирующее устройство |
Wk |
выбрано |
пра |
вильно. Если хотя бы для одного значения /г априорные значе ния неизвестных параметров функции Ак (t) значимо отличаются от тех значений, которые доставляют минимальное значение функ ционалу (IV.29), тогда необходимо проверить, выполняется ли для выбранного корректирующего устройства соотношение (IV.28). Если для выбранного на промежутке времени [0, T - J корректи рующего устройства выполняется соотношение (IV.28), тогда на
промежутке, равном [0, 7\ + |
АТХ] |
остается то же корректирую |
щее устройство, которое было на |
промежутке времени [0, Тг]. |
|
Если же соотношение (IV.28) |
для |
корректирующего устройства, |
которое было выбрано на промежутке, равном [0, Тг], нарушается, то среди корректирующих устройств выбирается то корректирую
щее устройство, |
для |
которого |
имеет место соотношение |
||||
|
|
7-,+ДГ, |
|
|
|
|
|
|
|
\ |
M[Yk(t)-MYhWdt |
= |
|
||
|
= |
m i n |
j M.[Yt(t) |
— MYt{t)fdt. |
|
(IV.30) |
|
|
|
i = l , P |
T 2 |
определяется |
тем временем, ко |
||
Здесь величина времени АТг |
|||||||
торое необходимо для решения |
соотношения (IV.30) |
в управляю |
|||||
щей ЦВМ. |
, |
|
|
|
|
|
|
На промежутке времени, равном |
\ТХ + ATъ |
2Т], |
вычисляется |
соотношение, аналогичное (IV.30), с интегрированием по проме
жутку [Тг |
+ АТг, 2Т]. |
Определяется то корректирующее устрой |
||||
ство Wk, |
которое удовлетворяет условию (IV.30) на промежутке |
|||||
времени, равном [Тг |
+ |
A T 1 ( |
2Т\. Указанное устройство |
остается |
||
для всех |
моментов |
времени |
£<Е [Тг+ |
ATly 2Т], где |
Тг удов |
|
134 |
неравенству |
Тх <3 Т. |
|
|
||
летворяет |
|
|
Предположим теперь, что промежуток времени |
[О, |
Тх] |
|||||||||
больше, |
чем |
[О, |
Т]. |
Тогда |
на промежутке времени |
[0, |
Тг) |
бе |
|||
рутся |
априорные |
значения |
для |
неизвестных |
параметров |
af t j S |
|||||
функции |
Akj |
(f). |
Уточнение |
проводится на промежутке |
времени |
||||||
[ 7 , |
2Т]. |
|
|
|
|
|
|
|
|
|
|
Определение минимального объема данных при работе системы |
|||||||||||
управления. Согласно |
условию на |
промежутке |
(Т, |
2Т) |
процесс |
X (t) имеет независимые приращения, причем известны априор ные значения его моментов распределения.
Определим тот минимальный объем массива статистических данных Rx, который необходим для того, чтобы определить ста тистические характеристики процесса X (t), t£_ [Т, 2Т] с тре буемой точностью. Так как желательно для целей управления получить достаточно большой массив статистических данных
процесса X (t) за достаточно малый |
промежуток времени, то |
выбор шага дискретности процесса |
X (t) подчиним единствен |
ному условию: шаг должен быть минимальным, причем ошибки от достаточно мелкого шага дробления процесса X (t) не должны превосходить допустимой величины.
Предположим, что процесс X (f) и его приращения нормально распределены. Исходя из канонического разложения (1.22) про цесса с независимыми приращениями, определим величину e;l (t) ошибки от аппроксимации разложения (1.22) конечной суммой,
содержащей |
пх |
членов разложения: |
|
|
|
|
п, |
t |
|
|
|
X (t) = X(Г) + j |
xv \ cpv (s) dp (s), |
(IV.31) |
|
|
v = I |
T |
|
|
|
* e i 7 \ |
274, |
|
где p, (s) = |
К (s, |
s) определяется формулой (IV.25); xv |
— ортонор- |
мированные случайные величины, имеющие нормальный закон
распределения; |
cpv (s) — любая |
ортонормированная |
система. |
Ве |
|||
личину Е п 1 (t) |
определим из |
условия |
e; l l (t) < |
гп1 |
(2Т) |
< |
s0 , |
где s 0 — допустимое значение величины |
ошибки. |
|
|
|
|
||
Используя |
датчики случайных чисел, |
для величин |
xv |
опре |
делим выборочные значения. Построим одну реализацию про цесса, имеющего разложение (IV.31). Принимая во внимание ре
зультаты, приведенные в п. 8, |
определим время Т2, при |
котором |
|
массив статистических,данных |
RA. дает оценку для моментов рас |
||
пределения процесса X (/),. |
t£ [Т, 2Т] |
с точностью |
не ниже |
допустимой. Выбирая время Т2 |
минимальным, получим минималь |
||
ный объем массива статистических данных |
Rv . |
|
Указанную процедуру по определению минимального объема массива статистических данных лучше всего подготовить и про вести в период проектирования системы управления. Однако в некоторых случаях это можно сделать и в период работы си стемы управления на отрезке времени, равном [О, Т].
13.5
Если же использовать- в качестве критерия точности оцени вания не величину средней квадратической ошибки, а величину ошибок первого или второго рода при статистической проверке гипотез, то процедуру определения минимального объема мас сива статистических данных можно упростить. Для этого рас смотрим статистику
где величины z (^) определяются по формулам (11.37); tt — мо менты времени, удовлетворяющие условию
|
Т < t l |
< t2 <•• •< |
t m ^ 2Т. |
|
|
|
Статистика |
x 2 имеет |
распределение |
хи-квадрат с |
г—1 |
степе |
|
нями свободы. |
|
|
|
|
|
|
Используя |
реализацию (IV.31), вычислим |
значение |
(IV.32) |
|||
и обозначим его через |
%2. Тогда с вероятностью не |
ниже, чем |
||||
1 — р\ имеет |
место неравенство |
|
|
|
|
|
|
|
1<%1 |
|
|
|
(IV.33) |
где значение %1 определяется из таблиц распределения |
хи-квадрат. |
|||||
Зададим величину р\ равную, например, 5, |
1, 0 , 1 % . Фикси |
|||||
руем значение |
%1> определим то значение г — 1, соответствующее |
объему массива статистических данных, уменьшенному на одну
единицу, при котором имеет место |
неравенство (IV.33). |
|
|
Предположим, что время Т2 |
получения минимального |
объема |
|
статистических данных массива |
Rx |
= (X (^), . . ., X (tm)) |
удов |
летворяет неравенству Тг < 2Т. Тогда выбор корректирующих устройств сводится к следующему. На отрезке времени, равном
[0, Тг], |
принимаются априорные значения неизвестных пара |
|||
метров |
матрицы |
Ak (t) |
(k = 1, р). Выбор |
корректирующего уст |
ройства |
Wk на |
отрезке |
времени, равном |
[0, Т2], определяется |
исходя из формулы (IV.28) с интегрированием по отрезку вре
мени [0, Т2]. |
, |
|
|
|
|
|
|
Для уточнения априорных значений матрицы A k (t) с помощью |
|||||||
массива статистических данных |
R^ = |
(Y |
(tj, . . ., |
Y |
(/ т )), |
по |
|
лученного на промежутке, равном |
[Т, |
Т2], |
с шагом дискретности, |
||||
определяемым шагом дискретизации процесса X (t), |
t£ |
[Т, |
2Т], |
||||
составим функционал (IV.5) для |
каждого |
значения |
k (k — 1, р). |
Найдем минимальное значение функционала по всем неизвестным значениям параметров, входящих в матрицу Ak (г).
Условия выбора корректирующего устройства Wk аналогичны условию выбора корректирующего устройства, основанного на определении неизвестных параметров матрицы A k (t), дающей минимум функционала (IV.29).
136