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

книги из ГПНТБ / Оптимизация процессов грузовой работы

..pdf
Скачиваний:
10
Добавлен:
27.10.2023
Размер:
26.66 Mб
Скачать

ПРИЛОЖЕНИЕ II

 

end;

fo r

/

: =

1 step

1

u n til

r

do C E

[ / ]

: =

C

[ / ]

C T

[ / ] ;

a

: =

0;

 

 

b : = 0; c : = 1; d : = I;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E 7

: A P

 

: =

10 f

18

X 0 ,9 ;

p

 

'■=

a\

i f

c =

p

 

th e n

go

to

 

F I

else

go

to

F 3 ;

 

F 3

: p

\ =

p

+

1;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F I

: l

 

: =

 

d\

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F 5 : / : = / +

1; Q ; = Q R [c ] + Q C [d ] ; i f Q R [ p ] < Q th e n go to F 2

 

else a

: =

a

+

1;

go

 

to

£ 3 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F 2

: QQ

: =

 

Q R \ p

-|- QC

[ / ] ] ;

 

i f

QQ —

Q <C

e

th e n

go to

FA

else

go

to

£ 5 ;

 

FA

: A B

 

: =

 

0;

fo r

i

: =

z

+

1 step

1

u n til

k do

C K

[ / ]

; =

C K

[j ] -{- W

[г,

/ ];

 

A B

: =

 

A B

 

+

(<

[У I

 

X

( C R

 

[ j ,

с ]

+

С С

[ j ,

d ]

 

+

С К

 

[у ])

X

( C R

[/,

p ]

+

 

C C [ j , l ] ) ) / C I J]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

end;

 

 

 

 

 

 

th e n go

to

F 6 else

 

 

: =

 

 

 

 

 

 

 

 

 

ß : =

/;

 

 

 

 

 

i f

A

B

^

 

A P

A P

 

A B ]

a

\ =

p \

 

 

 

 

 

c o m m e n t

О пределение

 

ка н д и д а та

в

о п ти м а л ь н у ю

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

 

FQ

: i f

р

ф

 

с

th e n go

to

F 3

else

/

: =

 

b;

i f

/ =

 

d

 

th e n

go to

F 7

else

go

to

E 2 \

£2

: / : = / +

1;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

£ 7

: p

; =

 

c;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

£ 4

: p

: =

p

-)-

1;

i f

Q C

[1 ]

<

Q

 

th e n

go

to

£ 1

 

else

b ; =

 

b

1;

go

to

£ 2 ;

£ 1

: Q Q : =

 

Q R

[/?]

 

 

Q C

[ / ] ;

i f

 

Q Q — Q <

E

th e n

go

 

to

£ 3

else

go

to

£ 4 ;

£ 3

: A B

: =

 

0;

 

fo r /

; =

 

1

step

1

u n til

r

do

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b e g in

C K

[ / ]

: =

0;

fo r

i ; =

г

+

 

1

step

 

1

u n til

k

do

 

C K

[ / ]

: = C K

[ / ]

+

 

«7

[i,

Л ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A B

: =

 

A B

 

 

(■

]

 

X

( C R

 

[у, с ]

-)- С С

[ j ,

d ]

 

-)-

C K

[ / ] )

X

 

( C R [ j ,

p \

-(-

 

C C ( j, l ] ) ) I C [ j ]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

end;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i f

l

=f= d

 

th e n

 

go

to

 

£ 5

else

A P

 

: =

A B ]

a

 

 

p ]

ß ; =

 

/;

 

 

 

 

 

 

 

c o m m e n t

 

О пределение

 

ка н д и д а та

в

о п ти м а л ь н у ю

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

 

£ 5

: i f

l = f =

d

th e n go

to £ 2

else

 

c : =

d;

a

; =

ß;

 

 

 

 

 

 

 

 

 

 

 

 

c o m m e n t

 

О пределение

 

члена

опти м а льно й

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

 

 

 

 

fo r

j

: =

 

1

step

1

u n til

 

r

do C P

 

[y ] : =

C C

[/,

d ]

+

C R [j,

c ];

 

 

 

 

 

fo r

/

: =

 

1

step

1

u n til

 

r

doi f

C P

[ j ] ^

C E

[ j ]

 

th e n

 

 

 

 

 

 

 

 

 

 

go

to

£ 6

 

else

Т Е

[c,

d ]

: =

0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

co m m e n t

 

В кл ю ч е н и е

 

в ы б р а н но го

 

вар иа нта

с т р у к т у р ы

 

в

о п ти м а л ь н у ю

после­

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

go

to

£ 7 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

£6 : fo r

i ; =

1

ste p

 

1

 

u n til

z

1do

fo r

 

h : =

 

1,

h - j-

1

w h ile

h^. M

 

do

TN [i,

h] : = N C

 

[t,

Л ];

v : =

 

0; m : =

0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L A : m : = m - \ - \ ] n

: = 0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

v - j- 1;

 

 

£ 3

: n

: =

 

n

-f-

 

1;

i f

Т Е

[ m ,

n ]

Ф

0

th e n

 

go

to

 

L I

else

 

 

 

fo r

 

i : =

 

1,

 

i +

1

w h ile

i= g : z —

1

do

N C

[i,

v] : =

TN [t, n ] ;

 

 

 

 

13*

195

ПРИЛОЖЕНИЯ

comment Оформление заголовка для следующей таблицы, т. е. занесение оптимальной последовательности в массив {У};

N C

[г, о] : = т

- j -

S

[z] —

1;

 

 

 

 

 

 

 

 

LI

: if

п =f= М

then

go

to L3

else

if

m ф

H

[ z]

then go

to

L4 else if z=j= k then

 

 

 

 

 

 

 

 

 

 

 

begin

M : =

v;

z

: =

z

 

l;

go

to

A

 

 

 

 

 

 

end;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

begin

array

P P

[/i];

 

 

 

 

 

 

 

 

 

 

 

 

В :

ft : =

ft -\-

1;

P P

[ft]

: =

1;

for

i

: =

1

step 1 until ft do

 

P P

[/i]

: = P P

[ft]

X

(I

? ( i ) t

N C

[i,

ft]);

if ft ф

v then

end

go

to В

else

print

( P P,

N C )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

end

end

 

Приложение

А Л Г О Л -П Р О Г Р А М М А В Ы Б О Р А

ІЙ

О П Т И М А Л Ь Н О Г О Ч И С Л А

 

К О Н Т Р О Л И Р У Е М Ы Х П А Р А М Е Т Р О В

 

begin integer

а2,

с,

і,

},

k,

l,

т ,

п,

р,

s,

w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

read

(п,

w);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

begin

real

a l,

q,

 

t;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

real array

A

[1

: ии],

В

[1

: tu],

F

[1

: w ] ,

 

Q

[1

: w ] ,

R

[1

: w ] ,

T

[I ; ш],

X [ 1 : и],

У [ 1 : пУ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

read

(t,

q.

В,

 

Q,

R ,

 

 

T)\

for

i: =

1

step

 

1

until

w

 

 

do

 

 

 

begin

F

[£]

: =

R

[г]

 

 

X

T

[t];a l

: =

F

[t];

 

 

 

 

 

 

 

 

 

 

 

 

A [t] : = 1/exp

(al);

A

[t]

: =

 

A

[£]

X

 

В

[t]

 

 

 

 

 

 

 

 

 

 

 

end;

print (A);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

begin

i

:=

1;

j

: =

1; X

[ j ]

: =

 

A

[£]; / : =

/

+

 

1; X

[i]

: =

T

[/];

 

 

 

j :=

І +

1;

X

 

( j )

 

 

 

:=

 

Q

[»];/

: =

/

+

I;

X

[/]

 

 

 

: = В [£];/ : =

/ +

 

X [/] := F [£];

m : = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

end;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

begin

H

: i

: =

 

i

+

1;

for

k : =

1

step

1

 

until

 

n

do

 

Y

[ k]

: =

0;

 

 

 

begin ft

:=

 

1; Y [ft]

: =

A

[£];

ft

: =

 

ft +

 

1;

 

Y

[ft]

: = T

[£];

 

 

 

 

 

ft : =

ft +

1;

у

[ft] : =

Q

[£];

ft : =

ft +

1;

У

[ft]:=

В

[/];

ft : =

ft +

1;

 

 

y [ f t ] : = f [ £ ]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 9 6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ПРИЛОЖЕНИЕ

III

end

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

end;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

begin j : = 0;

H I

: j

:=

/

+

5;

s : =

/;

ft : =

ft 4 -

5;

p

: =

ft;

 

 

 

 

Y

[ft]

; =

X

[/]

+

F

[t];

/

; =

/

1;

k

: =

 

k

-

1;

 

Y

[A]

: =

X

[ j ]

+

В [fl;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

: =

j

_

1;

k

: =

 

ft -

 

i;

Y

[A]

 

 

X

[/]

+

Q

[*];'

 

 

 

 

 

 

 

 

if

Y

 

[ k ] ^ >

q then go to H2 else go to ЯЗ;

 

 

 

 

 

 

 

 

 

 

 

 

 

H 2

: ft : =

ft

3;

go to

D l;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Я З

: j

:=

j

-

 

1;

ft ~

k

-

 

1;

Y

[ k]

: =

X

[/]

+

T [i];

 

 

 

 

 

 

if

Y

[ A ] >

t

then

go

to

Я 4

else go

to

Я5;

 

 

 

 

 

 

 

 

 

 

 

 

Я 4 : k : = k — 2; go to D l;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Я 5

: k : =

k —

1;

j

: =

 

p\

a ;==

К

[/];

К

[A]

: =

1/exp

(al);

 

 

 

 

j

: = ; -

1;

К

[A] ; =

К

[ k]

X

Y

[/];

k

: =

p

 

 

 

 

 

 

 

 

 

 

 

end;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D l : j : =

s; if

/ =

m

then

 

go

 

to

D 2

else

go

to

Я1;

 

 

 

 

 

 

 

 

 

D 2 : c :

=

k]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

begin A : = 0;

D 3 ;

ft : =

A -j-

1;

p

: =

k;

j

: =

0;

/ ; = / - ) -

1;

 

 

 

 

D 6

: if

К

[A] >

X

[/]

then go to D 4 else go to D5;

 

 

 

 

 

 

 

 

 

D 4

: j : =

j

-j- 5; if j

>

 

m then go to D7 else go to D6;

 

 

 

 

 

 

 

D 7

: s : =

/;

 

j

: =

j

— 4;

k

 

k

 

1;

if

К

[A] ^

X

 

[/]

then

 

 

 

go to D9 else go to D8;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D9 : / : =

/ +

 

1;

ft : =

 

ft +

 

I;

if

Y [ f t ] ^

X

 

[/]

then

 

go

to

DIO

 

 

else

go

to

D8;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

DIO : j : =

j

3;

ft : =

 

ft —

1;

if

Y

[ft]=s: X

[/]

then

 

go

to

D14

 

else

go

to

D13;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D14 : /' : =

/ -J- 1;

ft: =

 

ft +

 

1;

if

 

Y

[ft] ^

X

 

[/]

then

 

go

to

DIO

 

 

else

go

to

D15;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D 8

: j : =

s;

go to

D 16;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D ll

: j

: =

/

+

1;

go to

D16;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D13 : j : =

/ +

 

4;

go

to

D16;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

15 : / : =

/

+

3;

go

to

D16;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D16

: k : = p \

 

X

 

[ j ] : =

Y

 

[ft];

/

: =

/

+

1;

ft : =

ft +

1;

X

[/]

: =

Y

[ft];

j

: =

/

+

1;

ft : =

 

ft +

 

1;

X

[/]

: =

Y

[ft];

j

 

: =

j

+

1;

ft : =

ft +

1;

 

X

[/]

: =

Y

[ft];

j

: =

j

+

1;

ft : =

ft +

1;

X

[/]

: =

Y

 

[ft];

 

 

 

 

Я6

; m

 

/;

go

to

Я7;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D 5

: s : =

/;

j

: =

/ + 1 ;

ft ; =

ft +

1; if

F

[ f t ] ^

X

[/']

then

go

to

Я8

 

else

go

to

Я9;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Я 9

: / : =

/ — 2;

go

to

Я12;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Я 8

: j

: =

j

+

 

1; ft : =

ft -[-

1; if F

[ft ] ^

 

X

[/' ] then go to Я10

else go to Я 1 1

90

197

ПРИЛОЖЕНИЯ

Я10

: ft : = ft +

2;

go

to

Я 25;

Я П

: }

:=

j — 3;

 

 

 

 

 

 

 

 

 

 

Я12

: if

/ =

0

then

go

to Я13

else

go

to

Я14;

 

 

 

 

 

 

 

 

 

 

 

Я 13

: /

: =

у +

1;

go

to H20;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Я14

: у : =

j

3;

if

7

[ft]s £

X

 

[/]

 

then

go

to

Я16

else go

to

Я15;

 

 

Я15

: У:=

j

-j- 4; go to

Я20;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Я 16

: у : =

у

+

1;

ft :=

ft +

1; if

 

К[ft]s=CX

[/] then

go

 

to

Я11

 

 

 

else

go

to

Я 17;

Я17

: j

:=

j -(- 3;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Я20

: if у =

s then go

to

Я18 else go to Я19;

 

 

 

 

 

 

 

 

 

 

 

 

Я18

: У: =

ret -j- s;

 

/ : =

m;

m

: =

ret +

5;

a2

: =

s —

1;

 

 

 

 

 

 

 

Я21

: X

[у] :=

X

[ f ];

j

~

j

 

1;

 

/ : = / + l ;

 

if

/ =

o2

then

 

 

 

 

go to Я22 else go to Я21;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Я22

: у : =

s;

 

ft :=

p;

X

[/']

: =

 

У

[ft]; j

:=

j

+

1;

ft

 

 

:=

ft +

1;

 

X [/]

: =

К

[ft];

у : =

/

+

 

1;

ft : =

ft +

1;

X

[/}

: =

 

Y

 

 

[ft]; j

j

+ 1;

ft : =

 

ft +

1;

X

[y] :=

Y

[ft];

j

 

: =

/ +

1;

ft ; =

ft +

 

1;

 

 

X

[/']

:= К

[ft];

Я25

: if

k —

c

then

go

to

Я23

 

else

 

go

to

0 3 ;

 

 

 

 

 

 

 

 

 

 

 

Я23 : print (Я, X); if

j =

w

then

go

to

Я24

else

 

go

 

 

to

Я;

 

 

 

Я19

: ft : = p;

X

[ j]

:=

Y

[ft]; A

: =

ft +

1; у :=

j

+

1;

X

[ j ]

:= Y

[ft];

ft : =

 

ft +

1;

У: =

j

+

 

1;

X

[у]

; =

 

К

[ft];ft:=

ft +

 

1;

j : = j

+

1;

 

 

 

X []]

: =

Y [ft];

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ft : =

 

ft +

1;

У :=

j

+

1;

X

[y]

:=

Y

[ft];

j ~

 

j

+

1;

 

 

 

 

 

 

 

 

if у =

s then

go

to

Я25 else

go

to

Я26;

 

 

 

 

 

 

 

 

 

 

 

 

 

Я26

: I

: =

s;

X

[у]

:=

X

[/];

j

 

: =

у +

1;

s :=

s +

1;

 

 

 

 

 

 

 

 

if s >

m

then

go

to Я27

else

go

to

 

Я26;

 

 

 

 

 

 

 

 

 

 

 

 

 

Я27 : У :=

j — 1; go to

Я6;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Я 7 : у : =

у +

 

1;

-X

[/ ] : =

0;

if

 

у =

re

then go

to

Я25

 

 

 

 

 

 

 

 

else

go

to

Я7;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Я24:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

end

end

end

1 9 8

 

Приложение

П Р О Ц Е Д У Р А - Ф У Н К Ц И Я Д Л Я П О Л У Ч Е Н И Я

IV

С Л У Ч А Й Н Ы Х В Е Л И Ч И Н

 

С П Р О И З В О Л Ь Н Ы М И З А К О Н А М И

 

Р А С П Р Е Д Е Л Е Н И Я

 

Н А ЗН АЧЕН И Е

 

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

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

независимых случайных величин, равномерно

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

в интервале

[0, 1 ].

Пусть имеется случайная величина £, равномерно распределенная в интервале

[0, 1 ],

и непрерывная случайная величина т], функция

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

которой F (t) и плот­

ность распределения f ( t).

Тогда равенство

 

 

 

 

 

 

 

 

 

 

 

 

(П .1)

 

 

 

— СО

 

 

 

 

 

при / (0 >

0 задает взаимно-однозначное

соответствие между

переменными

| и т|.

Пользуясь равенством

(П .1), имеем:

 

 

 

 

 

1) для

равномерного

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

на

(а,

Ь)

 

 

 

т

=

О

при

t <

а, / > 6;

 

 

 

( b — а)-1

при

а;с; t ^ . Ь ,

 

 

 

 

 

 

 

тогда

 

 

 

 

 

 

 

 

 

 

 

г| = а+ 1

(6 — а);

 

(П.2)

2) для нормального распределения величин с дисперсией о2 и математическим ожиданием т

(П .З)

3) для показательного распределения

/ ( 0

0 при t < <mln;

=

 

Я, exp ['— Я, (t — Ь)] при / > J'nun,

тогда

(П .4)

4) для распределения Релея

t

ПРИЛОЖЕНИЯ

200

ПРИЛОЖЕНИЕ IV

Процедура-функция rav

Процедура-функция norm

Рис. П .1 . Блок-схема процедуры-функции

тогда

г) = (— 2а2 Іп | ) 1/2 = а V — 2 In | ;

(П .5)

5)для распределения Вейбулла

f (t) =

1 exp (—

 

= batа—1 exp (— Ы а),

тогда

 

 

 

 

 

 

 

 

*1

 

 

 

 

(П .6)

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

типа

 

\ - =

f

(t),

где 0 ^

t < оо,

 

11 =

г)!,

если

 

 

ТІ2 ,

если Т)а <

 

 

 

 

 

здесь

 

 

 

 

 

 

 

 

 

% =

- —

1п|;

 

 

 

 

 

 

 

 

-(П.7)

 

ПЯ=

а +

й ^ 2

|і — 6 ) ;

'201

ПРИЛОЖЕНИЯ

7)для распределения Пуассона в качестве значения случайной величины вы

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

 

 

 

 

 

 

 

m + 1

 

 

 

 

 

 

 

 

і < ;

 

 

- 2(=0

 

 

 

(П .8)

 

 

 

 

і=і

 

 

 

 

 

 

8) для

логарифмически-нормального

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

 

 

 

/ Ш

-

М

c x p f

I g < - l g m \

-

М

c::Df

l g ^ - l g a \

(П .9)

П

)

a t - 2 л

P V

2о3

)

~

2 я Ы

р \

2b2

) ’

 

/ 12

тогда т) = 10111, где г]., = а -)- b £ ь

9) для гамма-распределения

/(0

=

t b- 1

exp

(

(6 — 1) аь

 

 

 

тогда

(П .10)

Н а основании формул (П .2)— (П.10) составлена процедура-функция, с помощью которой могут быть получены случайные числа по любому закону с заданной плот­ ностью распределения. Блок-схема этой процедуры-функции приведена на рис. П. 1*.

О Б Р А Щ ЕН И Е К ПРО Ц ЕД УРЕ И Е Е О СО БЕН НО СТИ

Процедура & х (a, b, k) включает в себя: а, Ь — переменные типа real; k

переменную типа integer.

В данной процедуре а, Ь — параметры законов; k — номер закона.

При первом обращении к процедуре нужно положить k = 0, чтобы происходило присвоение величинам начальных значений. При последующих обращениях к за ­

дается числами от 1 до

10 в зависимости

от номера закона. Переменные а, Ь,

k

входные параметры. В

начале основной

программы необходимо описать мі

и м2

как переменные типа

real.

 

 

О ПИ САНИ Е ПРОЦЕДУРЫ

При первом обращении к процедуре-функции & х величинам «1 и «2 присваива­ ются начальные значения. При последующих обращениях по значению k опреде­

ляется вид закона распределения и реализуется одна из формул (П .2)— (П.10). Во всех операторах, реализующих формулы (П .2)— (П.10), участвуют вспомога­ тельные процедуры-функции гаѵ и norm, с помощью которых получаются псевдо­

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

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

по законам:

равномерному

в промежутке [0, 1 ]

и нормальному с

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

ожиданием

0 и дисперсией

1.

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

202

 

Приложение

П Р О Ц Е Д У Р А П О Л У Ч Е Н И Я

V

Г И С Т О Г Р А М М Ы И С Р Е Д Н Е Г О В Р Е М Е Н И

 

Р А Б О Т Ы С И С Т Е М Ы

 

Н АЗН АЧЕН И Е

Процедура & г служит для получения гистограммы непрерывной случайной величины,

 

 

 

 

 

N

N

суммарного времени работы

системы, а также для

получения сумм ^ t~[t

t 5lt

N

 

 

 

 

і= 1

/=і

 

 

 

 

 

 

Для построения гистограммы необходимо знать Т’зад и Л г“, поэтому до получения

гистограммы в процедуре определяются также эти величины.

 

Блок-схема процедуры представлена на рис. П .2.

 

О Б Р А Щ ЕН И Е К П РО Ц ЕД УРЕ

И ЕЕ О СО БЕН НО СТИ

 

 

Процедура

(fc,

т , N ,

At, е, k.2, k l , Т ,

t2, t3,

14, б) включает в себя: 1с, N ,

At, е, k2, Т , t2,

13, 14, б — переменные типа

real; kl — переменную типа integer;

пг — массив типа

real.

 

 

 

 

 

Пусть в основной

программе:

 

 

 

1С — время

работы системы;

 

 

 

m [0 : 13] — гистограмма,

массив чисел;

 

 

 

N — количество испытаний;

 

 

 

k2 — постоянная;

 

 

 

 

kl — постоянная для переключателя;

 

 

Т — суммарное время работы системы за N испытаний;

 

 

 

N

N

 

N

 

 

 

г=і

/=і

 

/=і

 

в— погрешность;

б— наибольшее число, представимое в ЦВМ.

Тогда обращение к процедуре может быть записано и без формальных параметров.

Процедура имеет следующие

особенности:

1. При первом обращении к

обязательно нужно положить к2 = 0. В этом

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

 

2. б — любое большое

число

такое,

что

если 12 V ö

V

б, то соответ­

ственно

12 : =

б V

Ö : =

б V

й

: =

б.

 

 

 

 

 

 

3. Процедура !РГ выдает на печать массив из трех чисел: At, число реализа­

ций

N ,

начальное время хО.

 

 

 

 

 

 

 

 

 

В программе, где будет использована процедура 2РТ, нужно обязательно

предусмотреть

переключатель с двумя

метками

 

 

 

 

 

 

 

 

switch L

: =

/VI,

N2,

 

 

где

Ml

: процедура

получения

tc;

 

go

to

L

[AI];

 

 

 

N 2

; процедура

получения

характеристик

надежности

невосстанавливаемых

 

 

систем (0>н.в).

 

 

 

 

 

 

 

 

 

203

ПРИЛОЖЕНИЯ

ПРИЛОЖЕНИЕ V

 

1)

x99

((fc —х О ) / Д t -|-

tz >

entier

.1

+

 

 

<N

 

 

<—

 

+

 

u

s _>

_^

+

 

 

 

 

<N

II

II

 

•*«*

 

II

F

 

 

 

 

 

>—•

 

 

 

s

 

 

 

«О

II

 

 

<N

 

 

t

 

 

 

 

CO

 

 

<—

Ю

 

о

О

+

A

 

 

CO

«N

 

-*«4

 

II

 

 

 

 

CO

«О

 

 

«Э

II

 

 

II

СО

 

 

 

f

 

 

t

 

 

<—

*

Ю

 

о

о

+

<o

A

 

A

 

 

CO

 

II

 

 

 

 

Рис. П . 2 . Блок-схема процедуры-функции

2 0 4

2 0 5

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