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

Т. Используя теорему Котельникова, выберем шаг дискрет­

ности 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),

[Т,

2Т],

составим функционал (IV.5) для

каждого

значения

k (k — 1, р).

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

Условия выбора корректирующего устройства Wk аналогичны условию выбора корректирующего устройства, основанного на определении неизвестных параметров матрицы A k (t), дающей минимум функционала (IV.29).

136