Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Коваленко И.Н. Полумарковские модели в задачах проектирования систем управления летательными аппаратами

.pdf
Скачиваний:
28
Добавлен:
23.10.2023
Размер:
7.56 Mб
Скачать

P -а{t) = X

(«)] du -1- V (/) P %(/) — TlPg (/)

 

о

 

(2.

14)

 

 

PAV) = 4 ! J Рг* [:I -

F («)] du + J

Л (*)) •

 

io

/-o

)

 

 

(^D

 

 

Система (2. 14) может быть решена численными мето­

дами, например, методом сеток. При v(0'=v = const ре­ шение получается в замкнутой аналитической форме, если применить преобразование Лапласа и известные теоремы операционного исчисления (о свертке и др.).

Опуская промежуточные выкладки, запишем решение системы (2. 14) в изображениях:

^ o ( s )=={5 + X + il + v [ l - g - ( s - f Х+ г р ] } - і ;

P1(S)=

P0(S). ѵ [ і —

g ( s + \ +

 

т|)] .

 

 

 

 

s +- X+ т]

 

 

 

 

 

Рг {8)= Р й{8)

S + V +

Ц

 

 

 

 

 

 

 

 

 

 

 

 

P3(s) = ^ü(5)

ѵХ

1—fffs-hX-b7))

I

1

(2. 15)

5 Н- 1

5 + X+ 7]

S

+ V - h 7]

 

^ М = лр0(5)(і 4

ф —ffO+X-H)]

T

=T~ +

 

s + X+ 1]

 

 

 

 

 

s+v+T]

 

,

ѵХ 1—g(s + X+ ij) .

 

1

 

 

5 + T)

S + X-j- Tj

s

+ V

+

1)

 

 

 

 

 

 

Для системы, способной функционировать (по край­ ней мере на отдельных участках траектории полета) с не­ исправным, но отключенным (во избежание ложной адап­ тации) контуром самонастройки, множество благоприят­ ных состояний

2 = Ы ,

і = 0, 1, 2, 3,

.(2.16)

Следовательно, вероятность безотказной работы системы в изображениях

Дл«)=2 о^ (5)-

(2Л7)

80

Математическое ожидание времени до первого вне­ запного отказа (полной поломки) системы можно опреде­ лить из выражения (2.9), ие прибегая к обратному пре­ образованию Лапласа, а именно:

ТоT>Q(s) |г=о— [X-j-T]-|-v(l —g-(X-l-Ti))] 1X

vX

1 — f f ( X - H )

 

1

X

X + T]

V

+ 7)

Ѵ+ -Г)

 

 

(2. 8)

Если безотказная работа контура адаптации необхо­ дима на всей траектории полета, то множество благо­ приятных состояний образуется из со0(/), coi(l'), u>2(t). В этом случае вероятность безотказной работы адаптив­ ной системы в изображениях

 

 

^cM=2

 

^(s)-

(2Л9)

 

 

 

1=0

 

 

Математическое

ожидание времени

до первого не-

скомпенсированного

параметрического

отказа двухкон­

турной адаптивной системы равно

 

Т„= ^

Pi (s)|j~o= {''•+

т1+

ѵ [l — g (x-f-1!)] I 1X

1=

0

V [ l

 

4 - T))]

 

 

X X

g (X

(2. 20)

 

 

 

 

 

V - f

t]

При известном законе распределения времени само­ настройки, применяя обратное преобразование Лапласа, можно получить характеристики безотказности адаптив­ ной системы в оригиналах.

В важном для приложений частном случае экспонен­

циального закона

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

времени адаптации

g(t) =це-!*':

 

 

 

 

 

 

g (s-f- X-(- т|)= Je- (■s+x+7i)

ldt -

5 +

• (2. 21)

 

 

 

 

 

“1"^ “h Щ

В этом случае:

 

 

 

 

 

Po(s)

 

s + p- +

X

1 }

 

(2.22)

(s +

v + X + TJ) (s

-I-

(J. -I-

X +

 

7)) — V(J.

81

P x (s) =

 

 

 

V

 

 

(s +

 

V -I- X+ 4 ) (s + |x + X+ V)) — Vjj.

 

 

 

 

 

 

 

X(s + fx -I- X-I- 4)

P*{s) =

 

(5 +

V-I- 4) [(5 4 - V

4 - X+ 4 ) (s + (X+ X+vj)-f-V(x]

P 3 (s) =

 

 

vX [2 (s"4-4) + у + X+ |x]______

(s 4

- 4) (s 4 - v 4- 4) [(s + V + X+ 4) X

3

 

 

 

 

 

X (s + fx +

1

 

 

 

 

X 4~ ij) — vjx]

P 4 (S)

 

V [4 (s + 4) (s + у + 1 ) + XT)(S + V 4 - 4 }] -

 

s (s 4- 4) (s + V4- 1)) [(s + V + X+ ij) X

 

 

(S+

4+ X+

Tj) [— vX-Г) - -

1] (S+ 4) ( s + V+ T))— Хч) ( s + Tj)]X

 

 

 

 

X (s + (X+

X+ T[) — V|x]

П р и м е н и в

о б р а т н о е

п р е о б р а з о в а н и е Л а п л а с а ,

чіш:

 

 

 

 

 

см

см

с м

п о л у -

Ро(0-

Jh— е-(х+т])/.

V4- [X

g —(V+ 1J.+X + 7])/ .

 

V 4- (X

 

 

 

 

 

РАО-

. g—(Х+і)Ц.

V 4- и

g—(V+IXX+TJ)/ .

 

V4fx

 

 

 

 

 

РА0--

X(X 4- (X— у)__ e_(v+7))i

 

 

 

X(X 4 - (X— v) — V(x

 

+

 

 

 

 

_____

g -u +чх.

v \

— e

 

00

 

-(ѵ-Цх+Х+Т))/.

j

 

 

(v — X) (v -

 

([X 4- X) ((X 4- v )

>

I см

 

 

X([X 4 - X— V )

-(v + i))/

 

 

 

X([X 4 - X— v) — V|x

 

 

 

 

v( v4 - | ) . - X)

c - g x 7|W

I

v^_

g—(V+IJ-+X4'7))/*

 

(v — X) (v 4- [x)

 

^Gx+XHv+tx)

 

 

 

P4 ( f ) = l - e - * .

 

 

 

 

 

 

Если время адаптации распределено по закону рав­

номерной плотности, т. е.

 

 

 

 

 

 

 

,

а < / < / > ;

 

 

 

g(t) = \ b-

a

 

 

(2.24)

 

 

О,

t > b ,

* < а , J

 

82

то

gIs + Ь+ 1D =

f e - 1 * + x+1» [Y [t - a)-

V (t- b)\

 

dt==

_ е -я($+Х +т])

_ е -(6-д) (j +).+ T))j ______ 1_______

(2. 25)

 

(s +

X + -г]) (b

a)

 

 

гдеу(г') — единичная функция.

Если при этом условия полета предусматривают не­ обходимость деятельности КС на івсей траектории, то ма­ тематическое ожидание времени до первого попадания

системы в поглощающее состояние соз(/)

равно

 

 

 

 

 

/“О

 

s=0

 

 

 

 

Н ^ + 'П+ ѵ

 

е- 0 ( х+ч)

 

_ е Сb о ) (Х + Tpj

 

X

 

(Х +

т)) (Ь — а)

 

 

 

 

 

 

 

 

X

1X++ Т) т

1 −: — а(Х+т|)

__ р— (й— я) (Х+7))

I

- ,

 

т

(X +

т)) (Л —

а)

 

V +

Т]

 

 

 

 

 

 

 

 

 

(2.26)

При экспоненциальном законе распределения време­

ни адаптации

 

 

 

 

 

 

 

 

 

_J__

______ X (X + ц —

ѵ)_______

 

 

 

Тп

 

[х (X + (J. — ѵ) —

Ѵ(х] (ѵ +

1))

 

 

 

X + ц

 

 

-L -__ —

_________ __________ S --------------- .

 

(2.27)

(v — X) (v + fj.) (X +

-rj) (p-4 - X ) ( [i+ v ) (ѵ + ц + Х + т,)

 

 

Из последней формулы нетрудно усмотреть, как вли­

яют на Тп параметры X, г), ѵ, р.

 

 

 

 

Тп от

- На рис. 2.5

приводятся графики зависимости

скорости

процесса

адаптации

для различных

значений

мгновенной интенсивности ѵ потока параметрических отказов. При бесконечно большой скорости адаптации (р—ѵоо) эффект влияния параметрических отказов сво­ дится к минимуму, кривая функции Г„=/(р) асимптоти-

-7-

„ =

1

 

Г,

Х

1

чески стремится к шах 1

------

1 + —-----

 

 

 

Х+

7) L

v-hi)J

 

83

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

(см. рис. 2.5) дают дополнительный материал для суж­ дения о необходимом

конечном быстродей­ ствии контура самона­ стройки . в целях обес­ печения достаточно вы­ соких субоптимальных значений Тп при _фик-

сированных К, т], ѵ. Реальные самона-

Рнс. 2.5. Графики зависи­ мости математического ожи­ дания времени до первого параметрического отказа от скорости процесса адаптации

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

А (0 = 2

(0^/(0.

(2-28)

/62

 

 

где <?,(/') — эффективность адаптивной системы в t-м со­ стоянии в момент t (это может быть, например, условная вероятность выполнения требований технического зада­ ния к точности стабилизации).

Таким образом удается отобразить процесс реального функционирования системы с учетом (внешних возмуще­ ний и изменения работоспособности элементов с течением времени і[29—34].

Учет частичной приспособляемости адаптивной систе­ мы стабилизации к внезапным отказам элементов ОК (к ^.-вариациям вектора внутренних параметров) приводит

84

к изменению некоторых коэффициентов системы интегродифференциальных уравнений.

Обозначим суммарную интенсивность внезапных отказов, компенсируемых за счет деятельности КС, бук­

вой

Тогда система (2.14)

принимает вид:

 

Р0(і)— (ѵ (0 + 4

- Tl)jD()W +

f P*(t, ti)g{u)du\

 

 

 

 

 

0

 

 

d p \ ( t , u )

d p t

(é, u )

 

 

 

— S

+

^

«

O

i

 

 

A W = ^ « W - ( ^ W + 4 )P SW;

о

 

I СЧ

 

t

 

 

 

 

CNJ

P-i W =x j дП 1 ~

F {n)] du~\-v{t)P2{t) — r\P3{t)\

 

Л Р ) =

(»l - ö

PAt) + \ p \ V ~ F m d u +

 

 

 

л -ч р л ъ + р м

-

 

Изменится и граничное условие

 

 

 

 

Рі& 0)=[ѵ(^ + ?]я0(^.

(2.30)

Алгебраическое суммирование интенсивностей отка­ зов вытекает из того, что пуассоновский поток является самопроизводящим по отношению к операции суперпо­ зиции.

Изложенный здесь подход к построению и исследова­ нию математической модели надежности многоконтур­ ных адаптивных систем может быть использован для аналитического изучения надежности СНС с иной струк­ турной организацией (например, СНС с двухкратной инвариантностью (39] и ультраустойчивых самоорганизу­ ющихся систем, основанных на методе избыточных пере­ менных [17]).

Отметим еще возможность сведёния немарковского процесса к марковскому путем расширения пространства состояний [21].

Пусть закон распределения случайного времени само­ настройки аппроксимируется суперпозицией законов Эрланга различного порядка

т

/ (о= 2

k'

(2-31)

/−1

 

 

85

где С) ■— постоянные положительные величины, в сумме равные единице;

л — плотность вероятности закона Эрланга і-го по­ рядка с масштабным параметром щ.

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

точно взять пары чисел і, г

І ^ г ^ т ,

—1)

п

особое состояние 0. В момент

достижения состояния

0

разыгрывается случайная ситуация с возможными исхо­ дами 1, 2,..,,т, с соответствующими вероятностями Р ь

Р2,-,Рт-

Если имеет место г'-й исход, то процесс попадает в состояние і, 0. Из состояния /, г( 0 ^ гі—2) возможен переход® состояние і, г+1, а из состояния і, ki—1 — воз­ можен переход лишь в состояние 0. Интенсивность пере­ хода в обоих случаях равна ц. (Пример применения это­ го способа к построению системы уравнений, описываю­ щих эволюцию состояний двухконтурной самонастраи­ вающейся системы с синусоидальным пробным сигналом дан в конце главы).

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

ность потока параметрических отказов

=

 

и

На временных интервалах, соответствующих і|)(^)>'1, темп возникновения параметрических отказов опережает скорость адаптации и наблюдается отличное от нуля рас­ согласование между действительным и желаемым пове­ дением адаптивной системы.

Таким образом, изучение переходных вероятностей дает возможность проследить эффект запаздывания адаптации относительно изменений среды, наблюдаемый

при работе адаптивных

систем с обратной связью по

критерию качества.

/

2.5. НАДЕЖНОСТЬ КОНТУРА АДАПТАЦИИ

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

86

щищен от собственных ошибок: в двухконтурі-юй системе они уже ничем не могут быть скомпенсированы. Это справедливо для любой иерархически упорядоченной системы, поскольку выход из строя управляющего уст­ ройства высшего ранга не восполняется деятельностью подчиненных регуляторов. Вопрос о предельно допусти­ мой относительной сложности КС при фиксированной сложности и надежности ОК всегда возникает при синте­ зе адаптивных систем и представляет далеко не академи­ ческий интерес.

Введем понятие о правильно организованном контуре самонастройки. Будем называть КС правильно организо­ ванным, если он удовлетворяет трем требованиям:

■1) реагирует на все отказы ОК;

2)его собственные неисправности не приводят к то­ му, что отказавший ОК признается исправным;

3)его отказы типа «ложный сигнал» или «короткое замыкание» не приводят к ложной адаптации.

Если КС наделен первыми двумя свойствами, то к устройству идентификации КС применима формула, по­ лученная в 1962 г. С. И. Злочевским, Р. А. Мирным, А. В. Розановым:

о.к (О +

Р д - 1 )

 

(2. 32)

Ру.At)

Р о А Г )

Р л

 

 

 

 

 

 

где Руи (/')— вероятность

безотказной работы

устройст­

ва идентификации за время t;

 

Ро.к(0— вероятность

безотказной

работы ОК за

время t;

 

 

 

 

 

Рд— коэффициент доверия

(условная

вероят­

ность того,

что отказ

ОК

действительно

имел место,

если

об

этом

сигнализирует

устройство идентификации КС).

 

По формуле (2. 32) рассчитана табл. 2. 1 с двумя вхо­ дами.

Из таблицы видно, что при достаточно высоких зна­ чениях коэффициента доверия (Рд^ 0,9) вероятность от­ каза устройства идентификации контура самонастройки должна быть, по крайней мере, на порядок ниже вероят­ ности отказа собственно основного контура (т. е. когда компенсаторные механизмы не действуют).

87

Таблица 2.1

Лек

 

 

р.

 

1

0,9

0,99

0,999

0,9999

 

0,9

0,987

0,998

0,9998

0,99999

1

0,99 -

0,999

0,9998

0,9999

0

1

0,999

0,9999

0,9999

1

1

1

0,9999

0,99999

1

1

1'

1

1

1

1

1

1

1

Если ограничиться рассмотрением объекта контроля

и устройства идентификации без структурной

избыточ­

ности, то из формулы

(2. 32)

и табл. 2. 1 можно сделать

при фиксированном Рл вывод об относительной сложнос­ ти Су.ц/Со.к, а именно:

С у м _

А У-><

(2. 33)

с о.к

0 , 1,

^ ° - к

 

где Лу ]І, Ло.„ — суммарная интенсивность

(внезапных от­

казов устройства идентификации и ОК соответственно. Уместно вспомнить, что в корректирующих кодах, например, в. коде Хэмминга с исправлением одиночных

ошибок,

отношение минимального числа

проверочных

разрядов /е к числу информационных

разрядов in{k?ü

j«log2

tri)

почти всюду (за исключением

коротких слов

с

3)

является правильной дробью.

С

увеличением

длины двоичного слова это отношение убывает, что хоро­ шо прослеживается на графике (рис. 2.6). В основе это­ го (кажущегося внешним) сходства лежат, быть может, глубокие причины. По-івидимому, любая разумная проце­ дура идентификации должна приводить к простым (отно­ сительно испытуемого объекта) физическим воплоще­ ниям.

Третье требование относится преимущественно к ис­ полнительным устройствам КС и связывается с представ­ лением о необходимости достижения их практической безотказности. Едва ли можно доверять коррекцию уст­ ройствам, собственная частота отказов которых сравни­ ма с этой же характеристикой корректируемого объекта. Суммарная интенсивность отказов корректирующих це­ пей контура самонастройки должна быть величиной, по

88

крайней мере, второго йорядка малости по сравнению с величиной г). (Справедливости ради заметим, что в не­ скольких известных нам практических примерах синтеза

самонастраивающихся систем это условие оказалось вы­ полненным) .

Рис. 2.6. Зависимость отношения минимального числа проверочных разрядов k к числу информационных разрядов т от длины двоич­ ного слова пі

Из этих рассуждений следует, что выход приведенной таблицы можно рассматривать как нижнюю грань для вероятности безотказной работы всего КС при заданных Ро.к и Рд. Вопрос о выборе значений Рд, к сожалению, до сих пор остается нерешенным. Интуитивное представле­ ние'связывает уровень коэффициента доверия Рл с тре­ буемой вероятностью достижения тех целей, для кото­ рых создается адаптивная система стабилизации или управления.

Впротивном случае действие КС не будет вне подоз­ рений, и возникнет надобность в коррекции самого кон­ тура адаптации, т. е. создания еще одного уровня управ­ ления.

Всвязи с изложенным, известный интерес представляют выска­

зывания нейрофизиологов Э. А. Асратяна и П. В. Симонова [4], каса­ ющиеся некоторых общих принципов компенсации нарушенных функ­ ций в высокоорганизованных биологических системах: «...Многоэтаж­ ная регуляция не предполагает примитивного дублирования. Каждый новый уровень приносит в регулируемую деятельность новые свойст­ ва, новое качество. Подобный принцип обладает, по меньшей мере, двумя следствиями. Подчиненные регуляторы способны функционн ровать в случае выхода из строя «командующих» механизмов.

89

Соседние файлы в папке книги из ГПНТБ