Файл: Батяев, Б. Г. Исследование надежности и точности линейных систем автоматического управления.pdf

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

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

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

Добавлен: 29.10.2024

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

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

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

32

 

Предположим,

 

что до момента

 

I

закончился,

по крайней мере,

один цикл,

Л гд а

вероятность

этого

события равна

 

 

 

 

P { T „ < i } = Р { ( т Л т * ) < в = W C O ,

(2 .1 5 )

где Тп - Тл' +

Та'

-

продолжительность,

л гго цикла,

W ( t ) -функция

распределения случайной величины

Тл .

Представим функцию распреде­

ления продолжительности цикла в виде свертки функций распределений

F(t)

и V (t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2 .1 6 )

Дифференцируя функцию W ( t )

по

t

,

получим плотность распределения

вероятностей ?& (t)

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

Тп

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

v s ( t ) = W '0 0 = l V ( t - ? ) / ( ? ) с / г f

(2 .1 7 )

 

 

 

 

 

 

 

 

 

О

 

 

 

 

 

 

где j О)

и v ( t )

 

-плотности вероятностей случайных

величин Тп и Тп’

 

Допустим, что

все

случайные

величины Т, , Тг

 

Т п имеют одну,

и ту же функцию распределения

W ( 0

со средним значением f-M{Tn) и

дисперсией &-£>(Тя), Тогда

эти числовые характеристики выразятся

через характеристики случайных величин

Та

и Тп

следующим образом

 

 

 

 

 

 

Т = М ( т ' + Т п ) = т \ т " ,

 

(2 .1 8 )

 

 

 

 

 

 

S)^SO(Tn+Tn) = «£>'+ so”.

(2 .1 9 )

 

5 .

 

 

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

ожидание числа циклов.

Зная функцию распре­

деления W СО и плотность вероятности

 

u K t ) одинаково распределен­

ных случайных величин

Т{ ,

найдем математическое

ожидание числа

циклов Л i t )

на промежутке

времени

(

О , t

)

 

 

 

 

 

 

Я СО = М [ 2(0]

=]Р лРпСО,

(2.20)

 

 

 

 

 

 

 

 

 

 

 

л - 1

 

 

 

 

где

2 ( 0

означает

число циклов,

происшедших на промежутке времени

( о

, t

) ,

 

/)t (7 )

- вероятность

завершения точно

п. циклов.

Учитывая равносильность

с о

б

ы т

и й и

[ t n < t }

,

найдем


33

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

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

/ £ ( 0 .

 

 

 

 

 

 

 

Р { * К О * л } = p [ t a н } = W n ( i ) .

'

(2 .2 1 )

Здесь

2

 

,

tc - О -

моменты включений

элементов,

a

Wa (t) -

функции распределения

t ^

,

определяемые как кратные свертки функ­

ции

W i t )

 

 

^

 

 

 

|

 

 

 

 

 

 

 

R i ( 0 = j W n . t ( ^ - ^ ) o ( W ( r ) f

щ -t) = W U ) . (2 .2 2 )

 

 

 

 

0

 

 

 

 

 

 

 

 

 

Исходя из формулы

( 2 .2 1 ), в результате

преобразований

определим

вероятность

завершения

п

циклов

 

 

 

 

 

 

 

= р {

=

гг} = р(*гСО

-р {1 Ь )* п + * }= М я11)-Ым № 2 .2 3 )

Подставляя

в выражение

(2 .2 0 )

значение

р п ( t)

из ( 2 .2 3 ),

получим

 

 

 

СЮ

 

СО

 

 

 

о°

 

ОО

 

 

 

 

н -1

 

- Т Ф ^ ( * Щ л ф Ж л К л(*)-2(*0ЬГлм (2 .2 4 )

 

 

 

П-\

 

 

 

п -I

 

п ч

 

 

откуда

 

 

 

 

 

 

СЮ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Я С О = И 2 W n ( i ) .

 

 

(2 .2 5 )

 

 

 

 

 

 

 

п--1

 

 

 

 

 

Преобразуем выражение ( 2 .2 5 ) ,

принимая во внимание форцулу

( 2 .2 2 ) .

Тогда подучим

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

СЮ

 

 

QO

 

 

 

 

 

 

 

 

 

 

j*

 

 

 

 

Л Ц = М [ » ( ф £ W „ ( t ) - I 2 j W h _ , ( i - T ) d W ( t r ) =

 

 

 

 

Дт/

 

 

л = /

о

 

 

 

 

 

 

i

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

( (< - г ) d w ( ? ) =

 

+

 

 

или

 

 

 

 

 

 

*

 

 

 

 

 

 

 

 

 

 

= W ( t ) + j f i ( t - ? ) d u / ( T ) .

 

(2 .2 6 )

 

 

 

 

 

 

 

О

 

 

 

 

 

 


34

Вьапе указанное интегральное уравнение Вольтерра называется уравнением регенерации.

Функция Л (t) является основной характеристикой процесса вос­ становления и равна средне^ числу циклов, происшедших на интерва­ ле времени ( О , t ),

Другой характеристикой данного процесса является производная этой функции 2 (t) , называемая средней частотой циклов в единицу времени

оооо

 

 

г а ) = Л ' СО = £

W a ( t ) =

 

 

 

 

(2 .2 7 )

Дифференцируя обе части уравнения

(2 .36) по

t

, получим уравнение

для средней частоты

циклов

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

гц) - iS(t) + j г(< - т)ьЯ(т) d z .

 

 

(2.28)

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

6 .

 

Функция готовности.

Для характеристики готовности элем

к выполнению предназначенных функций в любой момент времени

/

применяется

критерий, называемый функцией готовности.

 

 

 

 

Определение. Функция готовности F(t)

есть

вероятность

того,

что в момент Бремени

t

элемент находится

в исправном состоянии.

 

Критерий

Г ( t)

млеет больное

практическое'применение,

так как

позволяет оценивать готовность восстанавливаемого элемента

к работе

в произвольный

момент времени. Исходя из описания математической мо­

дели

элемента многократного действия, приведем формулу для

опреде­

ления

r ( t ) .

 

 

 

 

 

 

 

 

 

 

 

Пусть

п

- й цикл заканчивается до момента

/

на промежутке

времени ( Z ,Z+dT). Тогда вероятность этого события выражается

формулой

 

 

 

 

 

 

 

 

 

 

 

 

 

р ( Z < 1а

< T * d z ]

= d h ra( ? )*

lSn ( n

d Г, *-<■,W

#• (2.29)

Вероятность того, что до

момента

t

произошло точно

л

циклАз и

в момент t

элемент исправен

равна

 

 

('

 

 

 

р { ц ‘ < < Д '.,} = p { C < t ‘ f > T ; . , j =

12. Г )

С


„чктпьая несовместимость

и

I

событий £ t n < ( < ( „ , , } и определение

пункции готовности Г ( О

, получим

 

 

о©

 

 

 

 

Г(0

=

JTaCO •

 

 

(2.31)

 

 

 

 

а = 0

 

 

 

 

 

 

Подставляя

в

выражение (2 .3 1 )

вместо jfn ( О ее значение

из

 

(2 .3 0 ),

найдем

 

 

 

 

^

 

 

!

 

 

 

Г ( о - ] Т Г п С < ) = £ ] р (‘ -

 

=■

 

 

 

Л=0

П=0

о

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

= i-FCO + J [ i - F ( t - r ) ]

г(г)^г: .

 

 

(2.32)

 

 

 

(7

 

 

 

 

 

 

 

 

Вычислим стационарное значение функции готовности

К г

в ре­

жиме длительной работы элемента,

когда

t

—— »

 

 

 

 

 

 

Л" =

 

А’ О ; .

 

 

 

(2 .3 3 )

 

 

 

jf -Н-ОО

 

 

 

 

 

 

Для нахождения предела воспользуемся (|юрмулой Смита [35]

в

следующем

виде

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

°°

 

.

 

 

йт \ [ 1 - Р ^ - г ) } г Ю < 1 г

 

 

 

 

 

(2.34)

i

0 0

о

 

 

 

о

 

 

 

 

Подставляя

в

выражение (2 .34)

значение

Т

из ( 2 .1 8 ) . получим

&/n

/ f i - P C t - T ) ] a ( ? ) ^

=

+ Т

0

 

 

(2 *35)

t -

д

 

 

 

т

 

 

 

Принимая во

 

внимание выражения (2 .33)

и ( 2 . 3 5 ) , находим стационарное

значение

Г {

t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ОО

 

 

 

 

 

 

К = О т Г ( 0 = - г Д - г г f f l - f t s ) ]

 

 

 

 

г

 

 

r ' +

Т 0

 

 

 

 

 


36

или

Т '

п

( 2 .3 0

 

Выражение (2 .3 6 ) справедливо для произвольных функций распре­

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

большого промежутка времени.

Определим теперь функцию готовности восстанавливаемого элемен­ та на интервале ( t , t + © ) . Эта характеристика позволяет оценить готовность элемента к выполнению предназначенных функций на задан­

ном промежутке времени, функция готовности r ( t , 6 )

на интервале

( t , t + 6 ) равна

вероятности исправного

состояния элемента на этом

интервале времени.

В работе [22] функция

V { t 1в )

рассматривается в

качестве критерия надежности вычислительных устройств. Выведем фор­

мулу для

вычисления функции готовности Г ( t, $ )

применительно

к ма­

тематической модели восстанавливаемого элемента.

 

 

Пусть

п - й цикл заканчивается до момента t

на промежутке

времени

( ?

, Т + d t ) . Как известно, вероятность

этого

события

опре­

деляется

формулой (2 .2 9 ) . Вероятность

того,

что

до момента t

про­

изошло

П-

циклов и на интервале ( t

, t + 6

) элемент работал

ис­

правно

выражается формулой

 

 

 

 

t

 

 

 

 

j p { T „ . 1 > i + e - T ) w K( r ) d t : .

(2 .3 7 )

 

 

С

 

 

Тогда

вероятность

исправного состояния

элемента

на проме­

жутке

времени ( t

, t + 6 ) равна сумме

вероятностей

 

п -о

Подставляя в выражение (2 .3 8 ) вместо

it,®) ее значение из (2 .3 7 ),