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

книги из ГПНТБ / Приоритетные системы обслуживания

..pdf
Скачиваний:
14
Добавлен:
25.10.2023
Размер:
15.93 Mб
Скачать

Рассмотрим

интеграл

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

[ 1 1 Я ^ (2A j ( s j

 

 

J I V m ) e s t

d s '

 

 

 

2лі

 

( s

_ S

 

где мы

берем уже

 

ту часть

разреза у точки sm,

 

которая

входит в

круг с центром sm

и радиусом

а. После замены

s на sm

+ z

получим

 

&m

=

i£üL\\

+

П — f — Y

A, (sm) гіі-т &

dz.

 

(4.7)

 

 

 

2ЛІ

LJ

' J J

z

+

sm

/3=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Однако

на верхнем

берегу

разреза у точки s'm

 

г=хеія,

 

а

на

ниж­

нем — z = xe~in.

С учетом этого факта

при рассмотрении

каждого из

интегралов (4.7)

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

получаем

 

 

 

 

 

 

 

 

 

 

- ^ -

 

Г f У

 

Л, (s.) x'/vm

(e

V

m

) / e - "

 

dx —

 

 

 

 

 

 

a

'

;>0

 

 

 

 

 

 

 

 

 

 

 

 

 

-

f —

 

У

 

( s

« ) x i /

V m

(e~1

 

~тУ

 

e ~ x t

d

 

 

 

 

 

J

x

\

s m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

/5*0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

я

 

J ^ J x +

s„,

 

 

 

 

 

vO T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

Г

Г

Г

Аналогичным

образом получим

выражение

для

QJJ

=

I -f-

 

 

 

 

 

 

о

;>o

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сумма

+

ô7m

с учетом

sm

=

s0

-]- iVm

 

дает

 

 

 

 

 

 

 

 

 

,

-z,

 

_ 2 e € '

?•

(x -\ - s0 )

COS V m < -|- V m sin mt

 

 

 

 

 

 

 

 

л

J

 

 

(x + s 0 ) 2 + V ^

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

У

 

( s m )

s i n J2-

e~xt

dx.

 

 

 

 

 

(4.8)

 

 

 

 

 

/3=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

З а м е ч а н и е .

Если в

(4.8)

в интеграле заменить верхний предел

инте­

грирования a на оо, то

при

>-оо погрешность

от

произведенной

замены

будет

порядка

0 ( е ~ а

і ) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Далее, пусть имеем следующее

разложение:

 

 

 

 

 

 

 

 

-4- f0) cos Vmt -f Vm

sin Vmt

 

_

y

 

 

BT(t)xl.

 

 

(4.9)

181

Подставим (4.9) в (1.8). На основании сделанного замечания не­ трудно подсчитать при і->оо

2es°*

sin

\

ѵт

)

V m

t1

+ ( J / V „ ) - f 1

 

І>0

 

 

 

 

 

 

+ о (/<*•-*>o.

(4.10)

Вклад интервалов по разрезам у точки -s=So имеет несколько другой вид. Совершенно аналогично, находим

 

і2яі

J

J J

 

t J ^ x + s ,

v0

 

 

 

 

 

 

о / = 0

 

 

Ps„t

 

 

 

 

 

r ( l + - ! - + i )

 

j ] ( - s 0 ) - ^ ^ . ( S

; ) s i n

 

 

 

 

 

 

/ > 0

 

; > 0

v0

 

 

 

 

 

 

 

 

 

 

 

0 ( ^ 0 - « ) ' ) .

(4.11)

Так как a < s 0

— s*,

то

из

(4.4), (4.5),

(4.10) и (4.11) при t~+oo

получаем

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

соотношение

 

 

 

i - n f c ( * )

 

es,t

SS

 

 

 

 

 

 

 

 

x

sin - H

 

tl+Ulv0)+\

X

 

 

 

 

V q

 

X 2:SZ J

5BTf ( ^ . ( s j s h l

 

 

 

f 0(e&-e "),

 

 

 

 

 

 

f'+<//vmH-l

 

 

 

 

 

 

 

(4.12)

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

182

З а м е ч а н и е . На практике очень редкими являются случаи, когда на пря­ мой Re s = So кроме s0 у nh(s) есть и другие особые точки. Поэтому, вместо того

чтобы считать корни (3.5) и таким образом находить «о, расскажем об одном способе.

На прямой Im s= 0 функция nk (s) легко подсчитываете^ на ЭВМ. Тогда so

будет

ближайшим

к нулю

отрицательным

числом, в

котором

лк (s)

становится

равным бесконечности.

Разложение

Л&(А) в

окрестности

so

производится,

далее,

с помощью диаграммы

Ньютона. В некоторых случаях отсутствие других

особых

точек

ЯА(А) на прямой

Re s = s0 можно установить,

проверяя

устойчивость

мно­

гочлена R(,U),4)n)l(s—'So)

 

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

прямой.Re

s = s0

(см. [127]),

либо

про­

веряя

отсутствие

других корней уравнения

 

(4.5), лежащих

на

этой

же

прямой.

 

Иллюстрируя

сказанное

выше,

в этом

параграфе

мы

изучим

асимптотиче­

ское

поведение ф. р. Ri(t)

периода

занятости вызовов

приоритета

1.

Индекс 1

Пі(^)

в этом параграфе

будем

опускать.

 

 

 

 

 

 

 

 

 

 

Если длительность обслуживания имеет показательное распре­ деление с параметром Ь, то для я (s) получаем следующее квадрат­ ное уравнение (уравнение (3.4) в случае /г = 1):

аЬл2 (s) — (1 - L ab + bs) л (s) + 1 = 0.

(4.13)

Точки дискриминантного множества находятся как решения результантного уравнения, которое в свою очередь выводится ис­ ключением я (s) из (4.13) и из уравнения

2аЬп (s) (1 + ab + bs) = 0. Легко выводится результантное уравнение

R (Ф, Ф'я) = (1 + ab + bs)2 — 4ab = 0, откуда находим точки ветвления

(4.14)

о о

Обе точки ветвления являются точками ветвления второго порядка

для ri(s) и лежат

на отрицательной

полуоси

абсциссы,

причем

s0 <so. Интересующая нас ветвь яі (s)

легко находится

из

(4.13),

причем она должна

удовлетворять условию яі(0) = 1.

 

 

Как и раньше, обозначая

ni(s)=n(s) , имеем

 

 

 

/ ч

1 +

ab + bs V

(1 + ab + bsf — 4ab

 

 

я (s) = -

1

.

 

 

v ;

 

 

2ab

 

 

 

 

 

Выбираем контур интегрирования

Г так, как указано

на рис. 5,

где s 0 < s * < s 0 . Заметим, что на рисунке должно

быть р=0.

Функ­

цию я (s) в окрестности

точки ветвления s0 представим в виде

Я /S \ =

1

+ab + bs — Ь V(s

— ig) (s — î„)

 

 

^ '

 

 

2ab

 

 

 

 

183

Положим s0—70=---h. Выпишем разложение я (s) в ряд по степеням (s s0 ):

w

2a

l ab К

, . L

J _ ( s _ - } _

2ab

2a

 

^

"

/3=0

где

 

1 - f- ab -f- bs0

0

2ab

 

 

 

 

,1,

-

»'»

,

A,

- L

 

 

 

 

 

 

 

1

 

2a

 

2

2a

 

 

 

A,l

^

0,

 

A2i+i

=

 

^4=- •

~

r

^

- • - ( 4 - 1 5 )

 

 

 

 

 

 

2 a / я

П Г ( _ / +- r_ І Л

W

Теперь,

считая

по

формуле

(4.12),

находим

 

 

 

 

 

 

 

 

 

 

 

Г

/

2/ - f

1

1

 

 

 

So*

 

 

 

 

I / -1- ——

 

1 П (t)

:=

 

!>0

і^О S 0

 

'

т

 

г 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

2

 

 

 

 

 

 

 

~

 

 

 

 

3

 

 

 

 

О ( О = - ^ -

nikT (k + - | - ) Г к ~ + о ( О .

где обозначено

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H-s=ft

 

 

 

 

 

 

а числа

Л2/+1,

у -

0,

1, . . .

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

(4.15).

 

 

§5. Асимптотика времени ожидания

В§ 3 гл. 7 было отмечено, что особые точки пк-\ (s) являются особыми для преобразования Лапласа — Стилтьеса co?t(s) от ф. р. Wh(t) времени ожидания начала обслуживания вызовами приори­ тета k и что к особым точками функции coft(s) относятся и нули

уравнения

(3.7). Особые точки

ли (s)

будут

точками

ветвления

для (Oft(s),

а нули уравнения (3.7)

— полюсами. Так как

число

осо­

бых точек конечно, то теорема

§

2 позволяет

делать

выводы

об

их расположении относительно

начала

координат и мнимой

оси.

184

На рис. 6 указаны некоторые возможные конфигурации. Точки ветвления отмечены крестиками, полюсы — нуликами, причем оче­

видно, что все

особенности находятся в левой полуплоскости.

 

В случае

а) есть два полюса,

которые ближе к

мнимой

оси,

чем точки ветвления. Заметим, что

тогда обязательно

одна из

то-

 

IV

IV

IV

 

 

 

 

iR

X

ö

Рис.

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

Можно показать, что существуют соотношения между парамет­

рами a i , . . . , ük и

Ь\,...,

bh, при

которых имеет

место совпадение

полюса с точкой

ветвления.

 

 

Теперь, так же как и в § 4,

можно вывести

асимптотическую

формулу для 1 Wk(t)

при t-yoo

и любом k на основании в) рис. 2,

на котором указан и контур интегрирования Г в данном случае. Внутри контура Г могут оказаться полюса СОЙ(А) И, если нет сов­

падения точек ветвления на

Г с другими полюсами,

то для асимп­

тотической оценки 1 Wh(t)

будет

справедлива

формула

типа

(4.12) (конечно, с другими

коэффициентами), однако к правой ча­

сти такой формулы надо еще прибавить члены вида

 

 

Выч

<М*) est

1

 

 

где число q указывает, какой полюс

coft(s) из тех,

которые

оста­

ются внутри Г, берется. Совпадение полюсов и точек ветвления на контуре Г несколько изменяет результат, что будет показано на примере при k — 2.

185

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

co2 (s) = J l — а А ) [s +

a1_ — aln1(s)]{l-T-(s

+

a i • -fllJÏ! (s)) b2 }

(s- • a2 )

[1 + (s + ax — АІЛІ (s))

62] + a2

С помощью (4.13) и формул Виета можно исключить старшие сте­

пени lift (s),

если

знаменатель

o»2(s)

записать

в

виде

Л (s)яі (s) + ß ( s )

и умножить на A (s)m(s)

+B(s). Таким

образом,

мы вообще исключим

яі (s)

в знаменателе. Проделав все

эти

эле­

ментарные, но весьма трудоемкие выкладки,

получим

 

 

соа (s) =

Bs +

s^-f(Cs +

D)

/

s 2

4 - £ s + F

 

(5.1)

 

ßs -j-

s2

 

 

 

 

 

 

 

 

 

 

где обозначено

- flibi

A =

 

a,

1

 

1

2b1

 

2b,

2

 

 

2&1

 

 

_1

 

 

h

a,l"-a2

a,

 

 

 

 

 

 

 

В =

 

3

 

 

 

 

a, ôi

&1

 

 

 

 

 

 

 

 

C

= - î - ô .

^2

J

 

 

bi

 

 

 

 

2

 

02

B = T ( t - ' ) - t&i ( o A -

2 ( 1 + « A )

b, Oi&2

2 Ô!

«,2

3 1 - A . bi

l —

(5.2)

a A ) '

bi

 

 

 

 

 

(1

 

 

 

 

 

 

a

 

a 2 ( l

• albl — a2 b2 )

 

 

 

 

 

 

b i -

 

 

 

 

 

 

 

 

 

 

 

 

 

M 1 +

«1&2 —

+ 2 a 2 .

 

 

 

ß

 

 

b 2

(&i — bg)

 

 

 

 

 

 

 

Как

было

показано

в § 4, точки ветвления ( 0 2 ( 5 ) ,

являющиеся

корнями

дискриминанта

s2

+ Es + F=Q,

отрицательны.

Обозначим

их через

si и

s2,

причем

пусть

si>s 2 .

Корни знаменателя с\ и с2

дают нам полюса

юг (s),

где полюсом

может быть только тот ко-

186

рень, который не превращает в нуль числитель. Заметим, чго в за­ висимости от знака разности Ь\—Ь2 произведение с г С 2 = а может стать как положительным, так и отрицательным. Следовательно, оба корня знаменателя могут быть либо с одинаковыми (отрица­ тельными) знаками, либо с противоположными знаками, либо

А, //?'

д - — ш

s2

sf

СГС2

s,=c

 

4R

 

 

-iR

 

 

 

iR

 

A,

iR

 

 

O C,

 

Г

 

 

 

 

и s,

 

oC

src,=c2

s

s,

 

 

 

OC.

s"

 

 

 

OC2= CJ

 

 

 

oc,

 

 

 

 

 

 

-iR

-iR

 

 

 

 

 

Рис. 7

 

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

шениях между коэффициентами а ь

а%, Ьи

Ъ% все

эти случаи

воз­

можны. Все эти случаи приведены

на рис. 7, и здесь полюса с\,

с2

отмечаются нуликами, а точки ветвления Si, s2

крестиками.

 

Случай ж)

исключается, так

как он

противоречит

теореме

§ 2, гл. 7.

 

 

 

 

 

 

 

Случай а). Ближайшая особая

точка — полюс

Ci (рис.

7 а ) ) .

По формуле

(2.7)

 

 

 

 

 

 

l—W.(t)=l—Wt(0)

Г

s

e'1 ds.

(5.3)

 

 

2лі

J

 

 

 

s0too

187

Выбираем контур Г, как показано на рис. 7 а), и рассматриваем ин­ теграл

 

_ ! _ Г І ^ ? І £ ) e s ,

d s

- в ы ч ^ в "

 

 

и В ы ч - ^ - е *

 

 

 

2лі

J

s

 

 

 

 

 

s

 

 

s=0

 

 

 

 

s=0t

 

 

г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

to2(s)

определяется формулой (5.1). В таком случае

 

 

 

 

 

 

Выч

 

 

 

s=0

 

К

A +

DVF

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Л' Л -j-

B c r f с] +

(Ссх -j- D)V^

с? +

£ С і -f-

f

 

 

 

 

 

 

 

 

 

 

 

«1 (Ci — c 2 )

 

 

 

 

Аналогично

случаю,

рассмотренному

в § 4,

теперь можно

 

показать,

что

при ^ - > о о и при

t-+oo

справедливо

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

соотно­

шение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2ш'

S o + ' ° °

 

 

 

 

 

 

,

 

 

 

ft*

[2я(

 

 

 

J

s

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

М И (s)

es Ms

f

O(ss *0,

 

 

(5.5)

где

s 2 < s * < s r

Оценки

интегралов

l

h

]

по верхнему

 

и нижне­

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

му

берегам разреза у

точки ветвления

s1

производятся

так же, как и

в случае,

изложенном

 

в

§ 4.

Только

теперь

коэффициенты раз­

ложения

J£ïiiL

в р Я д будут

другими. Здесь

можно

записать

es, юAs)2

=

к

lA + Bs + 0] est

L

к

(Cs + D) / ( s -

S]) (s -

s2 ) g

s

t

 

s

 

 

s(a +

ßs +

s2 )

 

 

 

s (s — C!)(s c2 )

 

 

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

 

 

= 0.

 

В интеграле

 

 

& =

Г1

( C s + P ) V(s-sl)(s-si)

e S t d s

2лі И ]

s(s — c1)(s — c2)

 

188

сделаем замену s—Si=z

и положим Si — s 2 = li, после чего интеграл

I запишется в виде

 

1

 

2пі

(z + S l ) ( z — cî)(z — 4)

2яі

,1 fe>o

 

Нетрудно подсчитать коэффициенты ци, получающиеся при почленном перемножении рядов

Сг + Dt, где Dx ^ Cs - f D, ci = sl f

»'г (-»+т)

л > 0

Учитывая,

как

и в § 4

гл. 7,

что

на верхнем берегу разреза

z = xe

,

а на нижнем z=xe~ \

мы получаем

 

in

 

 

 

 

llг-ІЯ

 

 

 

 

 

^

^

Vhпкд*-кѵ.)

 

е-*'dx^j

+ О(e<s.-ô>< ),

 

 

 

 

о

 

 

 

 

где в

 

качестве ô берется

число

 

 

 

 

 

 

 

m i n ( s ^ s 2 , |Si — с ^ ,

Isx — c2 |,

— sx),

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

с7

Vh e s ' '

(5.6)

 

 

 

fc>0

 

Теперь из соотношений (4.3),

(4.5) и (4.6) находим

Wt(0)

189

fe>0

=

- К

у ^ * '

И*Г + T ) Г"~ "2 + О (e**'-0»).

 

 

2it

 

(5.7)

З а м е ч а н и е .

 

 

 

W 2 (0)

+

K

= 1 ,

(5.8)

 

 

 

 

 

a

 

 

так

как

V

A+DVF

равняется

полной вариации функции W%(t) на

величина д

 

интервале ( + 0, оо),

W2(0) есть

ее

скачок в

нуле, а полная вариация W2(t)

на

[—0,

оо)

равняется единице и дается суммой

(5.8).

 

В силу этого замечания член в (5.7), не содержащий экспоненты, равен нулю, что и должно было быть, так как 1—Wi{t)-*-0 при і->оо. Всюду в даль­ нейшем это замечание будем учитывать.

Случаи, представленные на рис. 7 в), д), з) и е), не очень отличаются от случая а). Сюда надо еще причислить и случай, когда оба полюса С\ и Сг отрицательны и расположены ближе к началу, чем точка ветвления st. Тогда к правой части (5.7) прибав­ ляется слагаемое qzec'i, а все остальное остается без изменения.

В случае в), так же как и в случае е), отсутствует слагаемое

В

случае

д)

в правой

части (5.7)

на месте qxeCit

фигурирует

сум­

ма

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-Всг

+ cï+(CcL

 

+ D)V

cî +

Ect

+F

 

 

q = K

 

 

 

•iV

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Положим c1 =

— и - f iV,

q = Qx -(-

iQr.

 

 

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i U = 2e-"'

(Qx cos

Vt - f Q2 sin W),

 

 

причем

и <

s r

Остальная

часть

(5.7) остается

неизменной.

месте

В

случае

совпадения

двух полюсов

с, =

(рис. 36), на

слагаемого qxe°J

в (5.7) стоит выражение

 

 

 

 

190

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