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

Текст из файла (страница 3)

11/09/10 21:09:09 Simulation in Progress.

11/09/10 21:09:09 The Simulation has ended. Clock is 50000.000000.

11/09/10 21:09:09 Simulation in Progress.

11/09/10 21:09:09 The Simulation has ended. Clock is 50000.000000.

11/09/10 21:09:09 Simulation in Progress.

11/09/10 21:09:09 The Simulation has ended. Clock is 50000.000000.

11/09/10 21:09:09 Simulation in Progress.

11/09/10 21:09:09 The Simulation has ended. Clock is 50000.000000.

11/09/10 21:09:09 Simulation in Progress.

11/09/10 21:09:09 The Simulation has ended. Clock is 50000.000000.

11/09/10 21:09:10 Simulation in Progress.

11/09/10 21:09:10 The Simulation has ended. Clock is 50000.000000.

11/09/10 21:09:32 SHOW anova (res,2,1)

11/09/10 21:09:32 ______________________________________________

11/09/10 21:09:32 ANOVA

11/09/10 21:09:32 _____________________________________________

11/09/10 21:09:32 Source of Sum of Degrees of Mean Square F Critical Value

11/09/10 21:09:32 Variance Squares Freedom of F (p=.05)

11/09/10 21:09:32 _____________________________________________

11/09/10 21:09:32 A 209051347.522 104525673.761 668.815 4.26

11/09/10 21:09:32 _____________________________________________

11/09/10 21:09:32 Error 1406564.869 9 156284.985

11/09/10 21:09:32 Total 210457912.392 11

11/09/10 21:09:32 ______________________________________________

11/09/10 21:09:32____________________________________________

11/09/10 21:09:32 Treatment

Level Count Mean Minimum Maximum 95% C.I. (SE)

11/09/10 21:09:32 A

11/09/10 21:09:32 ______________________________________________

11/09/10 21:09:32 1 4 1823.386 1786.401 1851.373 (1376.269, 2270.503)

11/09/10 21:09:32 2 4 1776.381 1714.503 1823.598 (1329.264, 2223.498)

11/09/10 21:09:32 3 4 -7054.076 -7908.097 -6476.600 (-7501.193,-6606.959)

11/09/10 21:09:32_____________________________________________

11/09/10 21:09:32 395.3289586

По результатам таблицы ANOVA видно, что параметр TAU является очень важным, так как его F-статистика намного превышает F_кр (668.815 > 4.26).

2.6 Отсеивающий эксперимент

Проведем отсеивающий эксперимент для определения наиболее важных факторов, влияющих на моделируемую систему. И выберем варьируемые величины исходя из результата опыта [1].

Текст программы GPSS:

TAU EQU 17

ZADEL EQU 3

OHE EQU 20

KRIT VARIABLE (2#N$ABC1+2#N$ABC2)-2#N$OTH-0.1#X$PROST-1#N$OTK

NORM1 VARIABLE (NORMAL(1,40,5))

NORM2 VARIABLE (NORMAL(1,60,10))

ERL3 VARIABLE (GAMMA(1,0,TAU, 3))

RAVNOM VARIABLE (UNIFORM(1,80,120))

GENERATE V$ERL3

ASSIGN 1,0

GATE NU STAN1,BBB

AAA1 SEIZE STAN1

ADVANCE V$NORM1

RELEASE STAN1

UNLINK BUF1,AAA1,1

TRANSFER .04,,DDD

EEE LINK BUF2,FIFO

DDD ASSIGN 1+,1

TEST L P1,2,OTH

TRANSFER ,BBB

OTH TERMINATE

BBB GATE NU STAN2,CCC

AAA2 SEIZE STAN2

ADVANCE V$NORM2

RELEASE STAN2

UNLINK BUF1,AAA2,1

TRANSFER .08,EEE,DDD

CCC TEST L CH$BUF1,OHE,OTK

LINK BUF1,FIFO

OTK TERMINATE

GENERATE 50000

TERMINATE 1

GENERATE ,,,1

FFF MARK

TEST GE CH$BUF2,1

SAVEVALUE PROST+,M1

UNLINK BUF2,ABC1,1

ADVANCE V$RAVNOM

TRANSFER ,FFF

ABC1 TERMINATE

GENERATE ,,,1

KKK MARK

TEST G CH$BUF2,ZADEL

SAVEVALUE PROST+,M1

UNLINK BUF2,ABC2,1

ADVANCE V$RAVNOM

TRANSFER ,KKK

ABC2 TERMINATE

Рис.1. Диалоговое окно генератора отсеивающего эксперимента

Рис.2. Процедура запуска модели

11/09/10 21:42:30 Model Translation Begun.

11/09/10 21:42:30 Ready.

11/09/10 21:42:30 SEM Procedure registered.

11/09/10 21:42:30 SEM_GETRESULT Procedure registered.

11/09/10 21:42:30 CHAST Procedure registered.

11/09/10 21:42:35 CONDUCT SEM()

11/09/10 21:42:35 **** Experiment in Progress. ****

11/09/10 21:42:35 Simulation in Progress.

11/09/10 21:42:35 A Simulation in an Experiment has ended. Clock is 50000.000000.

11/09/10 21:42:35 Simulation in Progress.

11/09/10 21:42:35 A Simulation in an Experiment has ended. Clock is 100000.000000.

11/09/10 21:42:35 "Run 1. Yield=1957.381969130512. TAU=10; ZADEL=3; OHE=5;"

11/09/10 21:42:35 Simulation in Progress.

11/09/10 21:42:35 A Simulation in an Experiment has ended. Clock is 50000.000000.

11/09/10 21:42:35 Simulation in Progress.

11/09/10 21:42:35 A Simulation in an Experiment has ended. Clock is 100000.000000.

11/09/10 21:42:35 "Run 2. Yield=1971.143408568981. TAU=10; ZADEL=3; OHE=20;"

11/09/10 21:42:35 Simulation in Progress.

11/09/10 21:42:35 A Simulation in an Experiment has ended. Clock is 50000.000000.

11/09/10 21:42:35 Simulation in Progress.

11/09/10 21:42:35 A Simulation in an Experiment has ended. Clock is 100000.000000.

11/09/10 21:42:35 "Run 3. Yield=1924.628865799248. TAU=10; ZADEL=8; OHE=5;"

11/09/10 21:42:36 Simulation in Progress.

11/09/10 21:42:36 A Simulation in an Experiment has ended. Clock is 50000.000000.

11/09/10 21:42:36 Simulation in Progress.

11/09/10 21:42:36 A Simulation in an Experiment has ended. Clock is 100000.000000.

11/09/10 21:42:36 "Run 4. Yield=1929.92535412587. TAU=10; ZADEL=8; OHE=20;"

11/09/10 21:42:36 Simulation in Progress.

11/09/10 21:42:36 A Simulation in an Experiment has ended. Clock is 50000.000000.

11/09/10 21:42:36 Simulation in Progress.

11/09/10 21:42:36 A Simulation in an Experiment has ended. Clock is 100000.000000.

11/09/10 21:42:36 "Run 5. Yield=-7357.309493947025. TAU=30; ZADEL=3; OHE=5;"

11/09/10 21:42:36 Simulation in Progress.

11/09/10 21:42:36 A Simulation in an Experiment has ended. Clock is 50000.000000.

11/09/10 21:42:36 Simulation in Progress.

11/09/10 21:42:36 A Simulation in an Experiment has ended. Clock is 100000.000000.

11/09/10 21:42:36 "Run 6. Yield=-7300.903184660247. TAU=30; ZADEL=3; OHE=20;"

11/09/10 21:42:36 Simulation in Progress.

11/09/10 21:42:36 A Simulation in an Experiment has ended. Clock is 50000.000000.

11/09/10 21:42:36 Simulation in Progress.

11/09/10 21:42:36 A Simulation in an Experiment has ended. Clock is 100000.000000.

11/09/10 21:42:36 "Run 7. Yield=-8553.776409936972. TAU=30; ZADEL=8; OHE=5;"

11/09/10 21:42:36 Simulation in Progress.

11/09/10 21:42:36 A Simulation in an Experiment has ended. Clock is 50000.000000.

11/09/10 21:42:36 Simulation in Progress.

11/09/10 21:42:36 A Simulation in an Experiment has ended. Clock is 100000.000000.

11/09/10 21:42:36 "Run 8. Yield=-7837.211632305294. TAU=30; ZADEL=8; OHE=20;"

11/09/10 21:42:36_____________________________________________

11/09/10 21:42:36 Alias Effect Sum of Degrees of F - for Only Critical Value

11/09/10 21:42:36 Group Squares Freedom Main Effects of F (p=.05)

11/09/10 21:42:36 _____________________________________________

11/09/10 21:42:36 A -9708.0708493249.342 1 1438.956 7.71

11/09/10 21:42:36 B -451.687 408041.624 1 3.115 7.71

11/09/10 21:42:36 AB -414.701

11/09/10 21:42:36 C 198.007 78413.745 1 0.599 7.71

11/09/10 21:42:36 AC 188.478

11/09/10 21:42:36 BC 162.923

11/09/10 21:42:36 ABC 167.156

11/09/10 21:42:36 _____________________________________________

11/09/10 21:42:36 Error 523972.270 4

11/09/10 21:42:36 Total 503676.981 7

11/09/10 21:42:36 Grand Mean -2908.265

11/09/10 21:42:36 ______________________________________________

11/09/10 21:42:36 Experiment ended.

Фактор А (TAU) является значимым, так как его F-статистика больше F_кр (1438,956>7,71). Факторы В (величина задела) и С (суммарная длина очереди) являются мало значимыми (3,115<7,71), (0,559<7,71).

2.7 Оптимизирующий эксперимент

Проведем оптимизирующий эксперимент, для того чтобы система нашла автоматически оптимальное решение и сравним его с уже полученным результатом [1].

TAU EQU 10

ZADEL EQU 3

OHE EQU 20

KRIT VARIABLE (2#N$ABC1+2#N$ABC2)-2#N$OTH-0.1#X$PROST-1#N$OTK

NORM1 VARIABLE (NORMAL(1,40,5))

NORM2 VARIABLE (NORMAL(1,60,10))

ERL3 VARIABLE (GAMMA(1,0,TAU, 3))

RAVNOM VARIABLE (UNIFORM(1,80,120))

GENERATE V$ERL3

ASSIGN 1,0

GATE NU STAN1,BBB

AAA1 SEIZE STAN1

ADVANCE V$NORM1

RELEASE STAN1

UNLINK BUF1,AAA1,1

TRANSFER .04,,DDD

EEE LINK BUF2,FIFO

DDD ASSIGN 1+,1

TEST L P1,2,OTH

TRANSFER ,BBB

OTH TERMINATE

BBB GATE NU STAN2,CCC

AAA2 SEIZE STAN2

ADVANCE V$NORM2

RELEASE STAN2

UNLINK BUF1,AAA2,1

TRANSFER .08,EEE,DDD

CCC TEST L CH$BUF1,OHE,OTK

LINK BUF1,FIFO

OTK TERMINATE

GENERATE 50000

savevalue kr,V$KRIT

TERMINATE 1

GENERATE ,,,1

FFF MARK

TEST GE CH$BUF2,1

SAVEVALUE PROST+,M1

UNLINK BUF2,ABC1,1

ADVANCE V$RAVNOM

TRANSFER ,FFF

ABC1 TERMINATE

GENERATE ,,,1

KKK MARK

TEST G CH$BUF2,ZADEL

SAVEVALUE PROST+,M1

UNLINK BUF2,ABC2,1

ADVANCE V$RAVNOM

TRANSFER ,KKK

ABC2 TERMINATE

Рис.3. Диалоговое окно генератора оптимизирующего эксперимента

Рис.4. Процедура запуска модели

11/09/10 21:57:57 Model Translation Begun.

11/09/10 21:57:57 Ready.

11/09/10 21:57:57 RSM Procedure registered.

11/09/10 21:57:57 RSM_MOVE Procedure registered.

11/09/10 21:57:57 RSM_FULLPLUSCENTER Procedure registered.

11/09/10 21:57:57 RSM_AUGMENTED Procedure registered.

11/09/10 21:57:57 RSM_CENTERTHEFACTORS Procedure registered.

11/09/10 21:57:57 RSM_GETRESULT Procedure registered.

11/09/10 21:57:57 CHAST Procedure registered.

11/09/10 21:58:02 CONDUCT RSM()

11/09/10 21:58:02 **** Experiment in Progress. ****

11/09/10 21:58:02 Simulation in Progress.

11/09/10 21:58:02 A Simulation in an Experiment has ended. Clock is 50000.000000.

11/09/10 21:58:02 Simulation in Progress.

11/09/10 21:58:02 A Simulation in an Experiment has ended. Clock is 100000.000000.

11/09/10 21:58:02 "Run 1. Yield=1961.381969130512. TAU=10; ZADEL=3;"

11/09/10 21:58:02 Simulation in Progress.

11/09/10 21:58:02 A Simulation in an Experiment has ended. Clock is 50000.000000.

11/09/10 21:58:02 Simulation in Progress.

11/09/10 21:58:02 A Simulation in an Experiment has ended. Clock is 100000.000000.

11/09/10 21:58:02 "Run 2. Yield=1938.857372592189. TAU=10; ZADEL=8;"

11/09/10 21:58:02 Simulation in Progress.

11/09/10 21:58:02 A Simulation in an Experiment has ended. Clock is 50000.000000.

11/09/10 21:58:02 Simulation in Progress.

11/09/10 21:58:02 A Simulation in an Experiment has ended. Clock is 100000.000000.

11/09/10 21:58:02 "Run 3. Yield=-1942.06772808578. TAU=20; ZADEL=3;"

11/09/10 21:58:02 Simulation in Progress.

11/09/10 21:58:02 A Simulation in an Experiment has ended. Clock is 50000.000000.

11/09/10 21:58:02 Simulation in Progress.

11/09/10 21:58:02 A Simulation in an Experiment has ended. Clock is 100000.000000.

11/09/10 21:58:02 "Run 4. Yield=-1513.918579316969. TAU=20; ZADEL=8;"

11/09/10 21:58:02 Simulation in Progress.

11/09/10 21:58:02 A Simulation in an Experiment has ended. Clock is 50000.000000.

11/09/10 21:58:02 Simulation in Progress.

11/09/10 21:58:02 A Simulation in an Experiment has ended. Clock is 100000.000000.

11/09/10 21:58:02 "Run 5. Yield=1908.306844669816. TAU=15; ZADEL=5.5;"

11/09/10 21:58:02 Simulation in Progress.

11/09/10 21:58:02 A Simulation in an Experiment has ended. Clock is 50000.000000.

11/09/10 21:58:02 Simulation in Progress.

11/09/10 21:58:02 A Simulation in an Experiment has ended. Clock is 100000.000000.

11/09/10 21:58:02 "Run 6. Yield=1950.7224334872. TAU=15; ZADEL=5.5;"

11/09/10 21:58:02 Simulation in Progress.

11/09/10 21:58:02 A Simulation in an Experiment has ended. Clock is 50000.000000.

11/09/10 21:58:02 Simulation in Progress.

11/09/10 21:58:02 A Simulation in an Experiment has ended. Clock is 100000.000000.

11/09/10 21:58:02 "Run 7. Yield=1913.700634287524. TAU=15; ZADEL=5.5;"

11/09/10 21:58:02 Using Model:

11/09/10 21:58:02 Y = -1374.78 -552.588 A +3477.12 B

11/09/10 21:58:02 +9.01347 A B

11/09/10 21:58:02 -81.1763 A^2 -324.705 B^2

11/09/10 21:58:02 Predicted optimum yield is 1812.81.

11/09/10 21:58:02 Optimum is in the local Experimental Region.

11/09/10 21:58:02 RSM_FitSurfaceToData() returns 4.

11/09/10 21:58:02 Experiment ended.

Оптимальное значение критерия V отражается в матрице LAB_BESTYIELDS, и оно равно 1812,810. Соответствующие этому результату: TAU равен 12,736, а ZADEL равен 5,531.

В результате проведения эксперимента получили следующую функцию регрессии, связывающую реакцию (критерий оптимизации V) с факторами TAU (фактор А) и ZADEL (фактор В):

Y = - 1374,78 – 552,588 A +3477,12 B + 9,01347 A B – 81,1763 A^2 – 324,705 B^2

Проведем эксперимент для полученных значений:

TAU = 15

ZADEL = 6

OHE = 20

START TIME END TIME BLOCKS FACILITIES STORAGES

0.000 50000.000 41 2 0

NAME VALUE

AAA1 4.000

AAA2 15.000

ABC1 33.000

ABC2 41.000

BBB 14.000

BUF1 10010.000

BUF2 10008.000

CCC 20.000

DDD 10.000

EEE 9.000

ERL3 10006.000

FFF 27.000

KKK 35.000

KR 10013.000

KRIT 10003.000

NORM1 10004.000

NORM2 10005.000

OHE 20.000

OTH 13.000

OTK 22.000

PROST 10011.000

RAVNOM 10007.000

STAN1 10009.000

STAN2 10012.000

TAU 15.000

ZADEL 6.000

LABEL LOC BLOCK TYPE ENTRY COUNT CURRENT COUNT RETRY

1 GENERATE 1094 0 0

2 ASSIGN 1094 0 0

3 GATE 1094 0 0

AAA1 4 SEIZE 747 0 0

5 ADVANCE 747 0 0

6 RELEASE 747 0 0

7 UNLINK 747 0 0

8 TRANSFER 747 0 0

EEE 9 LINK 1086 91 0

DDD 10 ASSIGN 71 0 0

11 TEST 71 0 0

12 TRANSFER 63 0 0

OTH 13 TERMINATE 8 0 0

BBB 14 GATE 508 0 0

AAA2 15 SEIZE 410 0 0

16 ADVANCE 410 0 0

17 RELEASE 410 0 0

18 UNLINK 410 0 0

19 TRANSFER 410 0 0

CCC 20 TEST 192 0 0

21 LINK 192 0 0

OTK 22 TERMINATE 0 0 0

23 GENERATE 1 0 0

24 SAVEVALUE 1 0 0

25 TERMINATE 1 0 0

26 GENERATE 1 0 0

FFF 27 MARK 501 0 0

28 TEST 501 0 0

29 SAVEVALUE 501 0 0

30 UNLINK 501 0 0

31 ADVANCE 501 1 0

32 TRANSFER 500 0 0

ABC1 33 TERMINATE 501 0 0

34 GENERATE 1 0 0

KKK 35 MARK 494 0 0

36 TEST 494 0 0

37 SAVEVALUE 494 0 0

38 UNLINK 494 0 0

39 ADVANCE 494 1 0

40 TRANSFER 493 0 0

ABC2 41 TERMINATE 494 0 0

FACILITY ENTRIES UTIL. AVE. TIME AVAIL. OWNER PEND INTER RETRY DELAY

STAN1 747 0.597 39.945 1 0 0 0 0 0

STAN2 410 0.496 60.502 1 0 0 0 0 0

USER CHAIN SIZE RETRY AVE.CONT ENTRIES MAX AVE.TIME

BUF2 91 0 47.680 1086 91 2195.234

BUF1 0 0 0.060 192 3 15.569

SAVEVALUE RETRY VALUE

PROST 0 983.181

KR 0 1875.682

FEC XN PRI BDT ASSEM CURRENT NEXT PARAMETER VALUE

1098 0 50004.141 1098 0 1

3 0 50004.194 3 31 32

4 0 50024.576 4 39 40

1099 0 100000.000 1099 0 23

3. Таблица характеристик процесса обслуживания

Максимальное значение критерия эффективности V будем искать, изменяя параметр TAU и величину задела [1].

V = С₁N_обр - С₂ N_отк - С₃(Т_пр3+ Т_пр4) - С₄N_{отн ,}

где:

N_обр –количество обработок

N_отк – количество отказов

Т_пр3 - время простоя первого станка

Т_пр4 - время простоя второго станка

С₁, С₂ , С₃,С₄ – весовые коэффициенты

Весовые коэффициенты:

Таблица 1 Таблица характеристик процесса обслуживания

Максимальное значение критерия эффективности (V=1712,65) было достигнуто при TAU = 16, а величина задела равна 3.

4. Графики зависимостей

Рис.5. Зависимость критерия эффективности от величины задела.

Рис.6. Зависимость критерия эффективности от интервала поступления.

Рис.7. Зависимость количества раз, когда станок занят от интервала поступления.

Рис.8. зависимость коэффициента использования станков от интервала поступления.

5. Анализ полученных результатов

В результате проведенного процесса моделирования были получены экспериментальные данные и построены графики зависимостей. Из графика зависимости критерия эффективности от интервала поступления сигналов (TAU*3), построенного по результатам значения таблицы характеристик процесса обслуживания видно, что максимальный критерий эффективности при данных коэффициентах был достигнут при интервале поступления сигнала, равным 16 (рис.6).

Решение найденное экспериментальным путем с помощью многочисленных прогонов (V=1812,81, TAU=15, ZADEL=6) оказалось близким со значением (V=1712,64, TAU=16, ZADEL=3) найденным программой после проведения оптимизирующего эксперимента, но не точное.

Проводя дисперсионный анализ видно, что параметр TAU является очень важным, так как его F-статистика намного превышает F_кр (668.815 > 4.26).

Также отсеивающий эксперимент показал, что параметр TAU является большим значением.

Заключение

Машинное моделирование – это эффективное средство решения сложных задач исследований, экспериментов и проектирования больших систем.

В данной работе основное внимание уделялось методам и этапам машинного моделирования в рамках общей методологии моделирования, изложенной в учебнике "Моделирование систем" [2].

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

Список использованной литературы

1. Моделирование систем: учебно-метод. Комплекс / А. И. Васильев; Дальневосточный государственный технический университет. – Владивосток: Изд-во ДВГТУ, 2008. - 172с.

2. Моделирование систем. Практикум: Учеб. пособие для вузов / Б. Я. Советов, С. А. Яковлев. – 3-е изд., стер. – М.: Высш. шк., 2005. – 295 с.: ил.

Характеристики

Список файлов курсовой работы