
- •Общие вопросы моделирования
- •1.1. Прогнозирование расчетных ситуаций
- •1.2. Общие сведения о моделировании систем
- •1.2.1. Виды моделирования
- •1.2.2. Математическое моделирование Виды математических моделей
- •Структурные математические модели
- •Функциональные математические модели
- •1.3. Подготовка данных и обработка результатов моделирования систем
- •1.3.1. Случайные величины, законы их распределения
- •1.3.2. Основные виды теоретических распределений
- •1.3.3. Определение параметров эмпирических распределений
- •1.3.4. Доверительные интервалы. Надежность. Обеспеченность
- •1.3.6. Подбор подходящего теоретического распределения Критерии согласия
- •1.4. Планирование эксперимента
- •1.4.1. Полный факторный эксперимент
- •1.4.2. Дробные реплики
- •1.4.3. Общая схема планирования эксперимента Крутое восхождение
- •Этапы планирования эксперимента
- •1.4.4. Стратегическое и тактическое планирование
- •Стратегическое планирование эксперимента
- •Тактическое планирование
- •1.5. Обработка и анализ результатов моделирования систем
- •1.5.1. Общие положения
- •1.5.2. Метод наименьших квадратов
- •Основы метода
- •Последовательность обработки данных
- •1.5.3. Практические приемы подбора кривых
- •1.5.4. Подбор эмпирических формул по кривым
1.3.3. Определение параметров эмпирических распределений
Экспериментально полученные массивы чисел представляют собой случайную выборку, лишь приближенно характеризующую генеральную совокупность и тем более теоретическое распределение. Поэтому при обработке данных наблюдений или эксперимента следует говорить не о параметрах распределения, а об оценке этих параметров или, как иногда называют такие оценки, о статистических параметрах. В дальнейшем статистические параметры (оценки теоретических параметров) опытных величин будут обозначаться звездочкой. Так, если математическое ожидание, дисперсия, стандарт случайной величины х обозначаются тх, Dx и ах, то для соответствующих величин, полученных опытным путем, будем использовать обозначения тx*, Dx* и аx*.
Естественной оценкой математического ожидания т случайной величины х является среднее арифметическое ее наблюдаемых значений (или статистическое среднее):
(1.48)
Для оценки дисперсии, казалось бы, можно было пользоваться формулой (1.20), записав ее применительно к дискретному распределению, или формулой (1.21), но так можно поступать только при больших выборках. Для получения несмещенной оценки дисперсии при относительно малых выборках вводится поправочный коэффициент п /(п- 1), т. е. используется формула для исправленной дисперсии:
(1.49)
Оценка среднего квадратического отклонения (стандарта) производится по формуле
(1.50)
Заметим, что среднее
квадратическое отклонение среднего
арифметического п
одинаково распределенных
взаимно независимых случайных величин
в
раз меньше среднего
квадратического отклонения σ
каждой из величин,
поэтому для оценки стандартного
отклонения используется формула
(1.51)
1.3.4. Доверительные интервалы. Надежность. Обеспеченность
Работая с массивом случайных чисел, необходимо знать не только их среднее значение, но и реальные пределы, в которых это значение может изменяться. Среднее значение чисел данной выборки будет представлять собой оценку математического ожидания генеральной совокупности, которую нетрудно вычислить по формуле (1.48). Среднее квадратическое отклонение (1.51) указывает лишь на усредненный разброс значений вокруг математического ожидания, но не дает возможности установить реальные пределы этого разброса. Подобные прогнозы могут быть основаны на интервальных оценках параметров распределений.
Надежностью (доверительной вероятностью) называют вероятность, с которой исследуемая величина находится в пределах доверительного интервала. Площадь, выделенная на рис. 1.6, а, определяет доверительную вероятность того обстоятельства, что величина х заключена в пределах доверительного интервала [а, b]. Концы интервала а и b называют доверительными границами. Поскольку эти концы сами являются случайными величинами, то правильнее следует говорить не о вероятности попадания случайной величины в доверительный интервал, а о вероятности покрытия этой величины доверительным интервалом.
Иногда при выборе интервала за аргумент вместо доверительной вероятности р-а берут величину, дополняющую ее до 1, т. е. 1-р (вероятность выхода за интервал). Эту величину называют уровнем значимости
(1.52)
Обычно применяют 5; 1; 0,3; 0,1% уровни значимости, что соответствует интервалам, которые покрывают случайную величину соответственно с вероятностью 0,95; 0,99; 0,997; 0,999.
В тех случаях, когда интервалы отсчитывают от одного и того же значения (обычно 0 или - ∞), говорят об односторонних интервалах, если же обе границы изменяются, то говорят о двусторонних интервалах. В дальнейшем условимся доверительную вероятность, связанную с двусторонними интервалами, называть надежностью, а с односторонними - обеспеченностью.
Выбор количественной оценки надежности и обеспеченности представляет собой сложную задачу, требующую глубоких профессиональных знаний в данной предметной области. Эту задачу нельзя формализовать, поэтому приходится полагаться на интуицию и накопленный опыт. Понятно, что при проектировании постоянного моста нужно принимать обеспеченность расчетных параметров существенно более высокой, чем при проектировании временного моста. Приняв обеспеченность для паводка равной 0,98, вы вправе ожидать, что в процессе эксплуатации 100 мостов, построенных по данному проекту, в двух из них паводок когда-нибудь может превысить расчетное значение. Последствия такого превышения паводка для постоянного и для временного моста не сопоставимы друг с другом. При выборе обеспеченности полагаются на интуицию, но опираются на накопленный и проверенный практикой опыт.
Если задан массив случайных чисел xi и есть основания полагать, что эти числа распределены по тому или иному закону, то определить границы доверительного интервала можно по таблицам, составленным для соответствующих законов распределения. Обычно таблицы составляют для нормированных величин (1.26), что соответствует точности оценки (расстоянию до границы доверительного интервала), равной
(1.53)
Для обратного перехода можно использовать формулы (1.28)......(1.31).
Для определения границ, соответствующих заданной доверительной вероятности, удобно иметь таблицы функции распределения в «перевернутом» виде, т. е. они должны давать значения х в зависимости от выбранного значения F(x).