
книги из ГПНТБ / Коваленко И.Н. Полумарковские модели в задачах проектирования систем управления летательными аппаратами
.pdfP -а{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