1/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 Value11/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.7111/09/10 21:42:36 B -451.687 408041.624 1 3.115 7.7111/09/10 21:42:36 AB -414.701 11/09/10 21:42:36 C 198.007 78413.745 1 0.599 7.7111/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 411/09/10 21:42:36 Total 503676.981 711/09/10 21:42:36 Grand Mean -2908.26511/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 10ZADEL EQU 3OHE EQU 20KRIT VARIABLE (2#N$ABC1+2#N$ABC2)-2#N$OTH-0.1#X$PROST-1#N$OTKNORM1 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$ERL3ASSIGN 1,0GATE NU STAN1,BBBAAA1 SEIZE STAN1ADVANCE V$NORM1RELEASE STAN1UNLINK BUF1,AAA1,1TRANSFER .04,,DDDEEE LINK BUF2,FIFODDD ASSIGN 1+,1TEST L P1,2,OTHTRANSFER ,BBBOTH TERMINATEBBB GATE NU STAN2,CCCAAA2 SEIZE STAN2ADVANCE V$NORM2RELEASE STAN2UNLINK BUF1,AAA2,1TRANSFER .08,EEE,DDDCCC TEST L CH$BUF1,OHE,OTKLINK BUF1,FIFOOTK TERMINATEGENERATE 50000 savevalue kr,V$KRITTERMINATE 1GENERATE ,,,1FFF MARKTEST GE CH$BUF2,1SAVEVALUE PROST+,M1UNLINK BUF2,ABC1,1ADVANCE V$RAVNOMTRANSFER ,FFFABC1 TERMINATEGENERATE ,,,1KKK MARKTEST G CH$BUF2,ZADELSAVEVALUE PROST+,M1UNLINK BUF2,ABC2,1ADVANCE V$RAVNOMTRANSFER ,KKKABC2 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 B11/09/10 21:58:02 +9.01347 A B11/09/10 21:58:02 -81.1763 A^2 -324.705 B^211/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 = 15ZADEL = 6OHE = 20START TIME END TIME BLOCKS FACILITIES STORAGES 0.000 50000.000 41 2 0NAME VALUEAAA1 4.000AAA2 15.000ABC1 33.000ABC2 41.000BBB 14.000BUF1 10010.000BUF2 10008.000CCC 20.000DDD 10.000EEE 9.000ERL3 10006.000FFF 27.000KKK 35.000KR 10013.000KRIT 10003.000NORM1 10004.000NORM2 10005.000OHE 20.000OTH 13.000OTK 22.000PROST 10011.000RAVNOM 10007.000STAN1 10009.000STAN2 10012.000TAU 15.000ZADEL 6.000LABEL LOC BLOCK TYPE ENTRY COUNT CURRENT COUNT RETRY1 GENERATE 1094 0 02 ASSIGN 1094 0 03 GATE 1094 0 0AAA1 4 SEIZE 747 0 05 ADVANCE 747 0 06 RELEASE 747 0 07 UNLINK 747 0 08 TRANSFER 747 0 0EEE 9 LINK 1086 91 0DDD 10 ASSIGN 71 0 011 TEST 71 0 012 TRANSFER 63 0 0OTH 13 TERMINATE 8 0 0BBB 14 GATE 508 0 0AAA2 15 SEIZE 410 0 016 ADVANCE 410 0 017 RELEASE 410 0 018 UNLINK 410 0 019 TRANSFER 410 0 0CCC 20 TEST 192 0 021 LINK 192 0 0OTK 22 TERMINATE 0 0 023 GENERATE 1 0 024 SAVEVALUE 1 0 025 TERMINATE 1 0 026 GENERATE 1 0 0FFF 27 MARK 501 0 028 TEST 501 0 029 SAVEVALUE 501 0 030 UNLINK 501 0 031 ADVANCE 501 1 032 TRANSFER 500 0 0ABC1 33 TERMINATE 501 0 034 GENERATE 1 0 0KKK 35 MARK 494 0 036 TEST 494 0 037 SAVEVALUE 494 0 038 UNLINK 494 0 039 ADVANCE 494 1 040 TRANSFER 493 0 0ABC2 41 TERMINATE 494 0 0FACILITY ENTRIES UTIL. AVE. TIME AVAIL. OWNER PEND INTER RETRY DELAYSTAN1 747 0.597 39.945 1 0 0 0 0 0STAN2 410 0.496 60.502 1 0 0 0 0 0USER CHAIN SIZE RETRY AVE.CONT ENTRIES MAX AVE.TIMEBUF2 91 0 47.680 1086 91 2195.234BUF1 0 0 0.060 192 3 15.569SAVEVALUE RETRY VALUEPROST 0 983.181KR 0 1875.682FEC XN PRI BDT ASSEM CURRENT NEXT PARAMETER VALUE1098 0 50004.141 1098 0 13 0 50004.194 3 31 324 0 50024.576 4 39 401099 0 100000.000 1099 0 233. Таблица характеристик процесса обслуживанияМаксимальное значение критерия эффективности V будем искать, изменяя параметр TAU и величину задела [1]. V = С1Nобр - С2 Nотк - С3(Тпр3+ Тпр4) - С4Nотн , где: Nобр -количество обработок Nотк - количество отказов Тпр3 - время простоя первого станка Тпр4 - время простоя второго станка С1, С2 , С3,С4 - весовые коэффициенты Весовые коэффициенты: Таблица 1 Таблица характеристик процесса обслуживания |
tau | zadel | abc | otk | oth | prost | V | | 5 | 3 | 998 | 1347 | 7 | 301,861 | -2202,9305 | | | 5 | 997 | 1363 | 6 | 374,428 | -2288,214 | | | 8 | 1001 | 1369 | 5 | 509,681 | -2364,8405 | | 10 | 3 | 997 | 0 | 6 | 458,7 | 1758,65 | | | 5 | 995 | 0 | 6 | 532,1 | 1717,95 | | | 8 | 1000 | 0 | 7 | 618,315 | 1683,8425 | | 15 | 3 | 1001 | 0 | 5 | 624,631 | 1684,6845 | | | 5 | 994 | 0 | 8 | 940,56 | 1509,72 | | | 8 | 987 | 0 | 4 | 977,747 | 1481,1265 | | 17 | 3 | 969 | 0 | 7 | 3357,478 | 252,261 | | | 5 | 981 | 0 | 2 | 1926,388 | 996,806 | | | 8 | 982 | 0 | 7 | 1405,702 | 1254,149 | | 20 | 3 | 831 | 0 | 4 | 17676,305 | -7180,1525 | | | 5 | 853 | 0 | 3 | 14555,111 | -5574,5555 | | | 8 | 832 | 0 | 4 | 16842,063 | -6761,0315 | | 25 | 3 | 680 | 0 | 6 | 32682,004 | -14987,002 | | | 5 | 681 | 0 | 4 | 31300,383 | -14292,1915 | | | 8 | 663 | 0 | 2 | 33400,359 | -15376,1795 | | 30 | 3 | 565 | 0 | 2 | 40937,576 | -19340,788 | | | 5 | 561 | 0 | 7 | 43501,113 | -20635,5565 | | | 8 | 582 | 0 | 6 | 38063,421 | -17873,7105 | | |
Максимальное значение критерия эффективности (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 с.: ил.
Страницы: 1, 2, 3
|