
- •Часть 1. Имитационное моделирование. Факторный эксперимент
- •Часть 1. Имитационное моделирование. Факторный эксперимент
- •Предисловие
- •Введение
- •Глава I. Системное моделирование – инструмент управления в больших системах
- •1.1. Понятия системного подхода и большой системы
- •1.2. Эффективность больших систем
- •1.3. Управление в больших системах
- •1.4. Структура систем управления
- •1.5. Основные понятия системного моделирования
- •Глава II. Моделирование систем методом статистических испытаний
- •2.1. Принципы построения математических моделей
- •2.2. Требования, предъявляемые к математическим моделям
- •2.3. Моделирование больших систем методом статистических испытаний. Сущность метода статистических испытаний. Точность метода
- •Вопрос 1.
- •2.3.1. Формирование непрерывных случайных величин с заданным законом распределения
- •2.3.2. Приближенные методы формирования случайных величин с заданным законом распределения вероятностей
- •2.3.3. Моделирование системы массового обслуживания
- •1,07 1,09 1,14 Моделируемое время t
- •2.3.4. Получение наблюдений при моделировании
- •Прикладные задачи имитационного моделирования
- •Ориентированный процесс случайного блуждания как метод прогнозирования
- •2.4.2. Модифицированный имитационным моделированием метод экспоненциального сглаживания
- •Глава III. Оценка качества моделей и планирование статистических испытаний
- •3.1. Оценка качества моделей
- •3.1.1. Методы повышения качества оценок показателей эффективности
- •3.1.2. Пассивные методы повышения качества оценивания показателя эффективности
- •3.1.3. Активные методы повышения качества оценивания показателя эффективности
- •3.1.4. Косвенные методы повышения качества оценивания показателя эффективности
- •3.2. Планирование имитационных экспериментов
- •3.2.1. Общая схема испытаний
- •3.2.2. Полные факторные планы испытаний
- •3.2.3. Дробные факторные планы испытаний. Планирование испытаний
- •3.2.4. Анализ результатов испытаний
- •3.2.5. Оптимальные планы
- •Методы принятия решений по результатам испытаний
- •Общая процедура принятия решений
- •3.3.2. Проверка гипотез о параметрах
- •Принятие решений о стабильности условий испытаний
- •Заключение
- •Библиографический список
- •Предметный указатель
- •Оглавление
- •Глава I. Системное моделирование – инструмент управления в больших системах 6
- •Глава II. Моделирование систем методом статистических испытаний 43
- •Глава III. Оценка качества моделей и планирование статистических испытаний 147
- •191186, Санкт-Петербург, ул. Миллионная, 5
3.2.4. Анализ результатов испытаний
Завершающим этапом испытаний по полному или дробному факторному плану является статистический анализ полученных данных, который включает: оценку воспроизводимости результатов испытаний; оценку значимости коэффициентов регрессии и уточнение вида модели; проверку адекватности модели.
Целью проверки
воспроизводимости
является установление однородности
результатов испытаний, проводимых на
различных уровнях
.
В ходе опытно-конструкторских работ
при проведении лабораторных и стендовых
испытаний обычно используется освоенная
аппаратура, обеспечивающая стабильность
условий. Поэтому предварительное
заключение относительно воспроизводимости
ожидаемых результатов часто может быть
сделано до начала испытаний. При
необходимости проверки воспроизводимости
содержание задачи совпадает с задачей
проверки гипотезы о стабильности
условий испытаний. Пусть результаты
испытаний представлены, как показано
в табл. 2.8.
Проверке подлежит
гипотеза
,
где –
количество строк матрицы планирования.
Условием проверки гипотезы является
наличие параллельных опытов в каждой
из строк.
Тогда в каждой строке могут быть вычислены
;
,
,
где – количество повторных испытаний.
Таблица 2.8
Номер опыта |
План испытаний |
Результаты испытаний |
Вычисляемые оценки |
||||||||
|
|
|
|
|
|
|
|||||
1 |
+ |
– |
– |
. . . |
+ |
|
|
. . . |
|
|
|
2 |
+ |
– |
+ |
. . . |
+ |
|
|
. . . |
|
|
|
|
|
|
|
|
|
|
|
. . . |
|
|
|
|
+ |
– |
+ |
. . . |
– |
|
|
. . . |
|
|
|
Для проверки
гипотезы можно воспользоваться критериями
Кохрена или Бартлетта. Если
,
расчетное значение статистики
критерия Кохрена определяется по формуле
(1.49) и гипотеза
принимается, если
,
где
,
.
Если дисперсии
однородны (принята гипотеза
),
то дисперсия опыта (или, что то же самое,
дисперсия воспроизводимости)
подсчитывается по зависимости
,
(3.2.11)
где знаменатель
характеризует число степеней свободы
.
В общем случае,
подсчитывается как среднее взвешенное
значение
.
(3.2.12)
Проверка
значимости коэффициентов регрессии
позволяет лучше осмыслить математическое
описание процесса, а также уточнить
вид модели путем отсеивания факторов,
слабо влияющих на значение выходного
параметра. Проверка значимости
каждого из коэффициентов производится
независимо, с помощью проверки гипотезы
0 по
-критерию.
Расчетные значения статистики критерия
можно определить по соотношению
.
(3.2.13)
Если
,
то коэффициент
является значимым и соответствующий
фактор оказывает существенное влияние
на величину
.
Статистическая незначимость
может быть вызвана следующими причинами:
интервал варьирования
был выбран слишком малым;
уровень начального режима по фактору оказался близок к точке частного экстремума
;
велика ошибка опыта из-за влияния неуправляемых и неконтролируемых факторов;
данный фактор (совокупность факторов) не оказывают заметного влияния на величину выходного параметра.
Поскольку план ортогонален и коэффициенты оцениваются независимо друг от друга, оказавшиеся незначимыми коэффициенты могут быть отброшены без пересчета остальных.
Проверка адекватности заключается в подтверждении предположения, что полученная математическая модель достаточно верно описывает характер процесса. Формальное содержание гипотезы состоит в том, что предсказанные уравнением (расчетные) значения выходного параметра отклоняются от опытных на величину, не превышающую некоторый наперед заданный уровень, и модель пригодна для обоснования инженерных решений. Для проверки гипотезы оценивается дисперсия адекватности
,
.
(3.2.14)
Если дисперсия
адекватности не превышает дисперсии
опыта
,
то есть основание полагать, что модель
адекватно описывает процесс. Согласно
п. 1.3 для проверки гипотезы о дисперсиях
используется
-критерий.
Статистика критерия
.
(3.2.15)
Модель считается
адекватной процессу, если
,
где
,
.
Если
,то
для получения адекватного описания
необходимо увеличить порядок
аппроксимирующего полинома. Очевидно,
что проверка адекватности возможна
лишь в том случае, если
,
то есть число разных испытаний
превосходит количество включаемых в
модель факторов.
Если гипотеза
адекватности принимается, то модель
может использоваться для управления
процессом (при доводочных испытаниях
для определения значений факторов, при
которых достигается экстремальное
или заданное значение выходного
параметра). Незначимые коэффициенты
могут быть отброшены, однако при этом
необходима основанная на анализе
физической сущности явления
уверенность, что эти факторы не влияют
на выходной параметр. В противном
случае следует стремиться продолжить
испытания (расширив, если позволяют
технические возможности,
по незначимым факторам).
Если модель неадекватна, возможны следующие решения: перейти к модели более высокого порядка (когда можно предположить, что опыты проводились в области, близкой к оптимальной); продолжить испытания, уменьшив интервал варьирования, или перенося центр плана в другую точку.
Рассмотрим примеры использования ДФП при организации испытаний технических систем.
Планируются испытания двигателя на надежность при эксплуатации. В качестве основных эксплуатационных факторов выделены: механические воздействия при транспортировании, воспроизводимые на вибростенде температура и влажность воздуха , имитируемые попеременным термостатированием и выдержкой в термобарокамере; старение при хранении, достигаемое проведением ускоренных испытаний узлов двигателя ; излучение и прогрев элементов с помощью секционных термобарокамер с использованием электронагревательных приборов и холодильной установки, В качестве выходного параметра оценивалась величина давления в блоке цилиндров двигателя
. Для испытаний выделялось
двигателей. Поскольку для полного учета всех факторов (при числе уровней
) необходимо
образца, было принято решение применить ДФП , причем
,
. Тогда обобщенный определяющий контраст
. Получаемая в этом случае математическая
модель принимает вид
Поскольку реализация
плана
возможна при
,
оставшиеся четыре образца использовались
для проведения контрольной серии
испытаний. План испытаний показан в
табл. 2.9.
Планируются испытания двигателя, проводимые с целью эмпирического определения коэффициента усиления, величина которого в общем виде представляется как
, где
– изменение -того регулируемого выходного параметра двигателя,
– изменение регулирующего параметра по -тому каналу регулирования. Исполнительными органами систем регулирования являются регуляторы расхода, регуляторы давления, дроссели. Величина
обычно включается в ТЗ на разработку.
Таблица 2.9
Номер опыта |
|
|
|
|
|
Номер опыта |
|
|
|
|
|
1 |
– |
– |
– |
– |
+ |
7 |
– |
+ |
+ |
– |
+ |
2 |
+ |
– |
– |
+ |
+ |
8 |
+ |
+ |
+ |
+ |
+ |
3 |
– |
+ |
– |
+ |
– |
9 |
+ |
+ |
+ |
+ |
+ |
4 |
+ |
+ |
– |
– |
– |
10 |
0 |
0 |
0 |
0 |
0 |
5 |
– |
– |
+ |
+ |
– |
11 |
– |
– |
– |
– |
+ |
6 |
+ |
– |
+ |
– |
– |
12 |
– |
+ |
– |
+ |
– |
Из практики опытной отработки известно, что расчетные значения обычно не подтверждаются в ходе испытаний. Поэтому возникает потребность настройки регулирующих органов с последующей опытной проверкой. Для обеспечения эффективности организации работ по определению и настройке коэффициентов усиления по каналам регулирования двигателя находят применение полные и дробные факторные планы. Рассмотрим случай, в котором для проведения испытаний реализуется ПФП при числе испытаний (табл. 2.10). В каждом испытании значения факторов и ( и ) определяют крайние положения органов в выбранной области изменения этих факторов. Матрица планирования и полученные из опыта результаты представлены в табл. 2.10. В нижней строке таблицы приводятся вычисленные по формуле (3.2.8) значения коэффициентов
Таблица 2.10
Режим |
План |
Результаты испытаний |
|
|
||||||
|
|
|
|
|
|
|
|
|||
1 |
+ |
– |
– |
+ |
0,7166 |
0,6769 |
0,7166 |
0,5785 |
0,6721 |
0,6939 |
2 |
+ |
– |
+ |
– |
0,8100 |
0,8100 |
0,8700 |
0,7545 |
0,8111 |
0,7893 |
3 |
+ |
+ |
– |
– |
0,6615 |
0,6076 |
0,5250 |
0,6071 |
0,6003 |
0,6221 |
4 |
+ |
+ |
+ |
+ |
0,6071 |
0,5533 |
0,4941 |
0,5375 |
0,5480 |
0,5263 |
|
0,6578 |
– 0,0837 |
0,0216 |
– 0,0478 |
|
|
|
|
|
|
Рассмотрим порядок статистического анализа результатов испытаний. Для проверки условия воспроизводимости по формуле (3.2.11) определим
;
;
;
.
Затем вычислим расчетное значение статистики критерия Кохрена:
Из табл. П.2.7---------?
при уровне
значимости
и числе степеней свободы
,
находим
.
Поскольку
,
принимается гипотеза об однородности
данных (воспроизводимости результатов
испытаний). Следовательно, дисперсия
испытаний может быть определена по всем
испытаниям согласно зависимости
(3.2.12):
.
Из (3.2.9) видно, что погрешность оценивания
.
Для проверки значимости коэффициентов и уточнения вида модели вычислим расчетные значения статистики -критерия по формуле (3.2.13):
;
;
.
Из таблицы Приложения
при
и
получим
.
Следовательно, для
и
имеет место
и эти коэффициенты значимо отличаются
от нуля. Поскольку
коэффициент
оказался незначимым. Поэтому фактор
из дальнейшего рассмотрения исключаем.
Уточненная модель принимает вид
Для проверки
адекватности модели определим
предсказанные этой моделью значения
;
;
;
.
Согласно зависимости (3.2.14) мера
неадекватности модели оценивается
дисперсией
.
Тогда определяемое по (3.2.15) расчетное значение статистики критерия
.
Из табл.
П.2.6------------? при
,
,
находим
,
что позволяет принять гипотезу об
адекватности модели изучаемому процессу
и использовать ее в дальнейшем для
настройки двигателя.