хибка
Оцінка похибки імітації, зумовленої наявністю в ІМ неідеальних генераторів випадкових чисел. Похибки, пов'язані з використанням у ІМ генераторів випадкових чисел, носять незначний характер, тобто результати моделювання не сильно залежать від вибору генератора випадкових чисел. З імовірністю 0.95 можна стверджувати, що точні оцінки компонент вектора відгуку потраплять в знайдений довірчий інтервал.
Визначення тривалості перехідного режиму
Як контрольованого параметра виберемо довжину черги ОСН3 надходить на пріборAA2. Проводити моделювання системи з кроком=1.
Tмод.12345678910L(OCH3)0.3350.4580.7190.9060.7560.8250.7550.7680.7720.729Tмод.11121314151617181920L(OCH2)0.7490.7420.7830.7770.7690.7720.7770.7750.7760.762
Видно, що довжина черги ОСН3 коливається приблизно біля значення 0.77, час пререходного режиму становить приблизно 5-6 одиниць. Наведемо графік цієї завісімотсі на рис. 3.
Рис. 3 - Залежність довжини черги ОСН3 від часу моделювання
Дослідження чутливості імітаційної моделі.
Розробка концептуальної моделі системи.
Час надходження заявок в систему моделювання позначимо як Тмод.
Визначимо вхідні і вихідні змінні для заданої системи моделювання.
Вектор параметрів Х: часи обробок заявок на пристроях АА1, AА2 - t (Aа1), t (AА2).
Хt (Aа1) t (AА2)
Вектор відгуку Y:
середній час обробки однієї заявки приладом АА1 - ave._time1.
середній час обробки однієї заявки приладом АА2 - ave._time2.
довжина черги L (OCH1).
довжина черги L (OCH2).
довжина черги L (OCH3).
YL (OCH1) L (OCH2) L (OCH3) ave._time1ave._time2
Оцінка чутливості:
Змінюючи значення вектора параметрів, певного концептуальною моделлю в пункті 3, отримаємо наступні значення вектора відгуків:
№ п / Пt (Aа1) t (АА2) L (Q1) L (Q2) L (Q3) ave. time1 (АА1) ave. time2(АА2)12.51.51.2400.2880.8012.4111.32222.11.50.2550.1060.4712.1301.56231.51.50.1850.1040.3541.5941.428421.80.2370.1130.4091.9451.569521.60.1950.1460.4852.0011.447621.10.7250.3360.2622.3691.043
(minXq, maxXq) - інтервал зміни q-ої компоненти вектора X.
(minY, maxY) - відгуки моделі, де minY і maxY означають відповідно вектори відгуку, отримані при мінімальному і максимальному значеннях q-ої компоненти вектора X.
Прирощення q-ої компоненти вектора параметрів моделі обчислюється за формулою:
,
яке і буде приростом вектора параметрів X при зміні тільки однієї компоненти q. Потім знаходиться прирощення n-й компоненти вектора відгуку:
.
Зміни вектора відгуку вибирається як. Отримані значення зафіксовані в таблиці:
11,7647119,44444-25,4826-17,0022-2,83831738,090185737,03704-115,217-78,838259,70549-16,84210532,449799
№ п / п,%,%,%,% 1 233301319252121312 31731228299324 51219251738255 63711579591732115
...
