2015-08-21
325
TEMPORARY Ran_Num; /*Визначення локальної змінної*/
/*Змінити початкове значення генератора випадкових чисел*/
Ran_Num = 248 # RN_Spost;
DoCommand("CLEAR OFF"); /*OFF викор. для збереження результата*/
/*Встановити початкове значення генератора випадкових чисел*/
DoCommand(Catenate("RMULT ",Ran_Num));
DoCommand("EP_Maz EQU ",RP_Maz);/*Встановити знач. 1-го фактора*/
DoCommand("EP_KamaZ EQU ",RP_KamaZ); /*Встанов.знач. 2-го фактора*/
DoCommand("START 100,NP"); /*Пропустити перехідний режим*/
DoCommand("RESET"); /*Почати збирати статистику*/
DoCommand("START 1,NP"); /*Запустити модель*/
END;
Модель розробленої системи наведена на рисунку 9.21.
Експеримент Eksp_Kor_TF розроблений для заповнення результатами моделювання глобальної матриці Rezult_TF і передачі її бібліотечній процедурі ANOVA. Процедура ANOVA призначена для виконаня багатофакторного (до шести факторів) дисперсійного аналізу, включаючи дво- і трифакторну взаємодію.
Процедура ANOVA має три аргументи. Перший – назва матриці результатів. Другий аргумент (необов’язковий) – розмірність матриці результатів, які можна використовувати для повторень. Кожен рівень цього виміру є спостереженням із різними початковими числами генератора випадкових чисел. Інформацію, пов’язану з числом повторень спостережень ANOVA використовує для оцінки стандартної помилки. Більша кількість додаткових спостережень (повторювань) підвищує інформативність оцінки.
|
|
|
Деякі плани експериментів не використовують розмірність для додаткових спостережень (наприклад, латинський квадрад). В цьому випадку другий аргумент процедури ANOVA вказується рівним нулю.
Рис. 9.21. Модель роботи транспортної фірми
Третій аргумент процедури ANOVA задає взаємодію факторів, які використовуватимуться у моделі. Якщо аргумент рівний одиниці, взаємодії не будуть враховуватись. Якщо аргумент рівний двом, в модель будуть включені лише двофакторні взаємодії. Якщо аргумент рівний трьом, в модель будуть включені дво- і трифакторні взаємодії.
Оператором INITIAL елементи матриці Rezult_TF ініціалізуються в стан UNSPECIFIED. Тоді процедурою ANOVA визначаються спостереження, які не були виконані.
Початкові значення факторів, кількість рівнів варіювання, значення інтервалів варіювання визначаються в процедурі експерименту.
Оскільки здійснюється двохфакторний дисперсійний аналіз (кількість автомобілів марки МаЗ – перший (А) фактор, кількість автомобілів марки КамаЗ – другий (В) фактор) і для кожного поєднання рівнів факторів виконуються три додаткових спостереження (N_Spost) для різних початкових значеннь (Ran_Num = 248 # N_Spost) генератора випадкових чисел, то процедура Plus-експерименту має три цикли WHILE / DO. Внутрішній цикл організований за числом додаткових спостережень, наступний внутрішній цикл – за кількістю рівнів 2-го фактора, а зовнішній цикл – за кількістю рівнів 1-го фактора.
|
|
|
У внутрішньому циклі за кількістю додаткових спостережень Plus-експерименту Eksp_Kor_TF перший оператор викликає процедуру запуску Eksp_Kor_Run. Ця процедура містить набір команд для ініціалізації і виконання одного спостереження. Кількість прогонів у ньому вказується процедурою DoCommand("START 100,NP"); /*Запустити модель*/.
Для передавання в процес моделювання команд CLEAR OFF, RESET і START процедура Eksp_Kor_Run може викликати бібліотечну процедуру DoCommand. За допомогою процедури DoCommand процесу моделювання також передаються команди для встановлення відповідних значень рівнів факторів.
Після завершення роботи процедури Eksp_Kor_Run експеримент Eksp_Kor_TF кожний раз зберігає результат (значення X$Prub) у відповідному елементі глобальної матриці Rezult_TF.
Процедура дисперсійного аналізу ANOVA розглядає кожний вимір матриці результатів як фактор. Розмірність матриці результатів повинна бути не меншою за кількість рівнів фактору або максимальну кількість додаткових спостережень (додаткові спостереження вказаються останнім індексом). Кількість рівнів кожного фактору не лімітуються, а лише обмежуються об’ємом віртуальної пам’яті комп’ютера. Нприклад:
