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

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

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

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

Добавлен: 02.07.2024

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

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

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

Глава VI

ДИАГНОСТИЧЕСКИЕ ПРОЦЕДУРЫ, ОПРЕДЕЛЯЮЩИЕ СОСТОЯНИЕ ЦВМ

25. СТАТИСТИЧЕСКИЕ ТЕСТЫ ДИАГНОСТИКИ

ЦВМ ПРИ РЕШЕНИИ ДИНАМИЧЕСКИХ ЗАДАЧ

Детерминированные методы программного контроля* Управ­ ляющая ЦВМ, входящая в контур управления автоматической системы, решет набор алгоритмов, связанных с управлением си­ стемой. Одним из основных требований, предъявляемых к алго­ ритмам, реализуемым с помощью ЦВМ, является точность выдачи результатов управляющей ЦВМ. Для проверки достоверности выданного результата, полученного в процессе работы системы управления, применяются программные методы системы контроля.

Метод двойного или кратного счета. Суть метода сводится к проверке правильности решения путем двойного или кратного просчета алгоритма. Недостатком метода является неполнота контроля, так как ЦВМ может работать неисправно, допуская

одну и ту же систематическую

ошибку.

Алгоритмический контроль,

использующий усеченный алго­

ритм решения.ч Отметим, что

лишь для отдельных алгоритмов

результаты счета основного алгоритма достаточно близки к усе­ ченному алгоритму. В этом состоит слабая сторона метода кон­ троля по усеченному алгоритму.

Метод, основанный на разработке аналитических условий,

позволяющих контролировать правильность вычислений. Эффек­ тивный контроль, основанный на проверке аналитических условий, возможен лишь для отдельных алгоритмов [15].

Метод

контрольного

суммирования.

Суть метода состоит

в том, что

подсчитывается

и запоминается

в ЦВМ контрольная

сумма всех чисел, входящих в контролируемый массив. Такого рода контроль не является полным, так как в случае наличия нескольких ошибок в контролируемом массиве контрольная сумма может совпасть, хотя и имеются ошибки в контролируемом мас­ сиве.

 

Указанные формы контроля

являются

детерминированными.

Они позволяют с

помощью специальных

программ, входящих

в

состав тестовых

задач,

проверить правильность

решения.

 

Постановка задачи статистической диагностики. Наряду с де­

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

методами

контроля

широко

используются

и

статистические

методы

[ 8 ] .

 

 

 

160


Рассмотрим следующую задачу. Пусть в управляющей ЦВМ реализуется алгоритм, основанный на вычислении операторного соотношения вида

F (t,Y

X (i), teT,

( V I . 1)

где F — известный оператор, зависящий от указанных аргумен­ тов; X (t) — известное входное воздействие; Y (t) — выходное воздействие.

Будем предполагать, что входное воздействие X (t) вводится в ЦВМ без ошибок. Такой случай может наблюдаться, например, когда устройство ввода работает исправно и его диагностика не проводится. Тогда в ЦВМ вместо операторного соотношения вида (VI.1) фактически реализуется оператор

 

 

 

Fk (t, Yk

(t)) =

X (f),

teT,

 

(VI.2)

где

индекс

k указывает

состояние

ЦВМ:

k = 0 —

состояние

исправности

ЦВМ,

k —'Л — состояние

неисправности;

F/t

t^>

Yk(t)\

оператор, фактически реализуемый в ЦВМ

вместо

оператора F

(t,

Y (t);

Yk (t)

— фактический

выход си­

стемы. .

 

 

 

 

 

 

 

 

Для

того

чтобы

провести стохастическую

диагностику ЦВМ,

с помощью которой следует установить, в каком состоянии нахо­ дится ЦВМ, подадим на вход оператора ( V I . 1) некоторое допол- . нительное случайное воздействие с известными статистическими характеристиками. Пусть дополнительное воздействие есть W (t);

тш (t),

Kw [t\,

t2)

— его

моменты

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

первых двух

порядков.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для

построения

теста

 

стохастической

диагностики

наряду

с оператором

( V I . 1)

рассмотрим

операторы

 

 

 

 

 

 

F

U,

Y, (t)}

-

X (t)

+

W (t),

t 6

T,

 

(VI.3)

если

F—линейный

 

оператор относительно Yx (t);

 

 

 

 

 

 

F(t,Y2(t)

 

=

W(l),

teT,

 

 

 

(VIA)

если

F—нелинейный

 

оператор

относительно

Y%

(t).

 

Здесь Y\

(/),

Y2

(i)

выходные

воздействия.

 

(VI.3),

Аналогично оператору.(VI.2), в ЦВМ вместо операторов

(VI.4),

реализуются

соответственно

операторы

 

 

 

 

 

Ft

U, Ykl

 

(t)]

-=

X (t)

+

W (t),

* (E T,

 

 

(VI.5)

 

 

 

Fk

It, Yk2(t)]

 

=

W{th

te

T.

 

 

(VI.6)

Задача стохастической диагностики ЦВМ состоит в том, чтобы по результатам выходных воздействий Ykl (t), Yk2 (t) построить диагностический тест, с помощью которого можно определить, в каком из двух состояний — исправном или неисправном — на­ ходится. ЦВМ, т. е. определить, какое значение 0 или 1 принимает индекс k.

11 Л. Т. Тарушкнна

161


Воздействие W (i) вырабатывается с помощью управляющей ЦВМ. Все необходимые данные о решениях, содержащих воздей­ ствие W (t), хранятся в памяти ЦВМ, а не передаются в контур управления автоматической системы.

Выработка дополнительного воздействия в ЦВМ. В качестве дополнительного воздействия W (f) возьмем нормально распре­ деленный процесс с известными моментами распределения. Рас­ смотрим два класса процессов:

1) стационарный процесс W (I), обладающий свойством эрго­ дичности;

2) процесс W (I), имеющий независимые приращения.

Для получения стационарного процесса W (t) в ЦВМ восполь­

зуемся

каноническим разложением

(1.37).

Тогда

 

 

 

 

 

"

 

 

 

 

 

 

' JIV

 

 

 

 

W(t)

xv cos cov

t -f- yx

sin ©у t,

wv =

^r-

 

(VI.7)

 

v = l

 

 

 

 

 

 

 

 

 

 

где xv,

yv некоррелированные

случайные величины, для

ко-

торых

 

 

 

 

 

 

 

 

 

 

 

 

 

М xvX

x v 2 = М yv\ yY2 = М xv yyi

= О,

 

 

 

 

 

 

Mxv

=

Myll=

 

О,

 

 

 

 

 

 

 

м 4 i

=

м

yli

=

i

 

 

 

 

 

 

 

K = f v 2 J

V, ( - 1 = 1 , л).

 

 

 

 

 

Выборочные значения для случайных величин xv,

г/д

получим

с помощью датчика

случайных

чисел,

учитывая

при этом,

что

xv< Ун—нормально распределенные

случайные

величины.

 

Для построения процесса W (/), имеющего,независимые прира­

щения, достаточно в каждый дискретный момент времени

опре­

деляемый последовательностью (11.34), положить

 

 

 

 

 

W(tt)

 

=

t

zv,

 

 

 

 

' (VL8)

 

 

 

 

 

v = l

 

 

 

 

 

 

 

где zv — независимые нормально распределенные случайные ве­

личины:

Mzv

= 0,

Mzl = 1.

Для

построения

диагностического теста из двух воздействий

(VI.7), (VI.8)

выберем то, с помощью которого получается наиболее

простое условие контроля с точки зрения его реализации на ЦВМ. Построение алгоритма диагностического контроля. Предполо­ жим, что для управления автоматической системой следует вы­ числить в ЦВМ оператор ( V I . 1) в некоторых точках tc 6 Т, обра­ зующих неубывающую последовательность (11.34). Для диагно­ стирования ЦВМ вычислим в ЦВМ либо оператор (VI.3), либо оператор (VI.4), в зависимости от того, линеен или нет оператор F. Вычисления проведем в некоторых точках t), взятых из последова­ тельности (П.34), при этом все необходимые данные для вычисле-

162


ния операторов (VI.3), (VI.4) в точке t) возьмем из оператора ( V I . 1), значения которого известны для всех моментов времени t. ^ t{. Таким образом, для некоторой последовательности моментов вре­ мени

 

 

 

 

0^t[<...<t'h^T,

 

 

 

 

 

(VI.9)

входящей в последовательность (11.34), имеем

решение

опера­

тора ( V I . 1)

и соответственно

либо оператора (VI.3), либо

опера­

тора (VI.4), т. е. для последовательности моментов времени

(VI.9)

имеем

последовательность

значений

выходных

координат

 

 

 

 

Уи (•'!),

Yk(t2)

 

 

Yk{t'h),

 

.

(VI. 10)

 

 

 

M '

i ) ,

Ykr{Q

 

 

Ykr(Q,

 

 

(

где

 

 

 

 

 

 

 

 

 

 

 

 

 

[ 1, если

F

линейный оператор,

 

 

 

' ~~ 1 2,

если

F — нелинейный

оператор.

 

 

 

Заметим,

что

значения

последовательности

( V I . 11)

являются

выборочными значениями некоторого случайного

процесса Ykr (t).

Так как оператор .F известен, а входное

воздействие W

(t)

имеет

известные моменты распределения, то можно предположить, что

известны и моменты распределения процесса

Ykr (t). Действи­

тельно,

данные

о значениях

моментов

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

процесса

Ykr (t)

можно получить до начала

работы системы управления,

используя методы определения моментов распределения, приме­

няемые на стадии проектирования

автоматических систем. Отсюда

следует, что условия диагностического контроля можно получить из условий, определяющих значения моментов распределения

процесса

Yk

(I).

 

 

 

 

 

 

 

 

 

1. Пусть Ykr(t)

стационарный

процесс,

обладающий

свой­

ством эргодичности.

Вычислим

в ЦВМ

интегралы

 

 

 

 

 

-

 

 

 

 

 

 

 

 

 

 

-р-

\ Ykr (0

dt

inykr,

 

 

 

 

 

 

 

hi о

 

 

 

 

 

 

 

 

 

4 -

h,

\

 

 

 

 

 

 

 

 

 

f [Ykr —tUykrYdl

=

оЦг-

 

 

 

 

 

 

hi о

 

 

 

 

 

 

 

Составим

функционал

 

 

 

 

 

 

V = [mukr~MYkr

 

( О ] 2 +

\<ttr-M

[Yl{r(t)-MYlir(t)]2}2.

 

(VI.12)

Функционал

( V I . 12)

будем рассматривать

как

условие

кон­

троля за

работой

ЦВМ.

А именно: если

V <

е, где

е > 0 — до­

пустимая точность, предъявляемая к работе ЦВМ, то ЦВМ ра­ ботает исправно, следовательно, k = 0; если V ^ 8, то ЦВМ ра­ ботает неисправно, следовательно, k — 1.

п *

163


Такие параметры, как величина отрезка [0, t'k\, количество точек в последовательности (VI.9), задание величины е, являются входными данными алгоритма статистического диагностирования

ЦВМ, их значения определяются

исходя из точности,

заложенной

в ЦВМ при ее проектировании,

и уточняются

в

процессе от­

ладки

ЦВМ.

 

 

 

2.

Пусть Ykr (t) — нормально

распределенный

процесс, имею­

щий независимые приращения. Рассмотрим задачу идентифика­

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

Ykr

(/) по выборочным зна­

чениям

реализации

( V I . 1.1) с

помощью

статистической

проверки

гипотез

(см. п. 8). Для

этого

составим случайные

величины

 

 

 

Zi =

Ykr{ix)

MYkr{t\),

 

 

=

Ykr (t,) -

Ykr

(/;_0 -

MYkr (t't)

+

MYkr (t't-i)

(t =

2, h).

Образуем статистики (11.45), которые определяют случайные

величины v, (t =

1,

/г).

 

 

 

 

 

 

Ставится гипотеза Я: значения v( являются выборочными зна­ чениями объема /г, полученными из генеральной совокупности, отвечающей случайной величине v, имеющей нормальный закон распределения с параметрами Mv = 0, Mv2 = Mz\.

Используя соотношения (11.48), (11.49), получим, что для ги­

потезы Я

критическая

область

уровня

а

равна

_

 

 

 

 

 

 

 

 

 

 

 

1

 

 

где

постоянная

с0

определяется

согласно

условию

(11.46);

для

гипотезы

Mv2 =

 

Mz\

область

принятия

гипотезы

 

 

 

 

*

<

- 5

г

| ) ( "

-

Т § '

< )

• < *

(V .. I4)

где

постояннные

 

с ъ

съ

определяются

согласно условию (11.47).

 

Условия ( V I . 13),

( V I . 14) являются

условиями

контроля,

по­

зволяющими производить диагностику ЦВМ. Действительно, если при заданном значении а выполняется неравенство, проти­

воположное

неравенству

( V I . 13),

а также выполняется

неравен­

ство ( V I . 14),

тогда гипотеза Я принимается, а это означает, что

моменты распределения

процесса

Ykr (t) соответствуют

моментам

распределения этого процесса, вычисленным в ЦВМ за отрезок времени, равный [0, ti\\. Таким образом, если гипотеза Я прини­ мается, то это означает, что ЦВМ работает исправно, т. е. к = 0; если гипотеза Я отклоняется, то это означает, что ЦВМ работает неисправно, т. е. к = 1.

3. Пусть Ykr (t) — мартингал имеющий нормальный закон распределения, приращения которого также нормально распреде-

164