
- •Математическое моделирование Учебное пособие
- •Донецк 2006
- •Содержание
- •Введение
- •1. Построение экспериментальных законов распределения
- •1.1. Общие положения
- •1.2. Статистические критерии согласия
- •Г) Критерий согласия Романовского
- •1.3. Построение закона Пуассона
- •1.4. Построение показательного закона
- •1.5. Построение нормального закона
- •2. Модели оптимизации
- •2.1. Принципы формирования моделей оптимизации
- •Задача производственного планирования
- •Задача оптимальной загрузки оборудования
- •Задача о смесях
- •Транспортная задача
- •2.2. Графический метод решения задачи линейного программирования
- •Алгоритм графического метода решения злп
- •2.3. Универсальный метод решения линейных задач оптимизации
- •Алгоритм симплекс-метода решения злп
- •Пример 2.3.1. Решить злп (2.2.1), (2.2.5) симплекс-методом.
- •Критерий оптимальности опорного плана
- •Переход к следующей симплекс-таблице осуществляют по правилам:
- •2.4. Двойственная задача линейного программирования
- •Свойства двойственных задач
- •2.5. Методы анализа конфликтных ситуаций с помощью матричных игр
- •Алгоритм принципа максимина (минимакса)
- •Решение. Этаматричная игра имеет размерность (3х4), т.Е. Игрок а имеет три стратегии, а игрок в – четыре. Запишем ее в нормальной форме.
- •Последовательность действий при решении игры
- •3. Регрессионный анализ
- •3.1. Однофакторные модели
- •3.1.1. Построение однофакторных моделей
- •3.1.2. Оценка качества моделей
- •Свойства коэффициента корреляции
- •Построение доверительного интервала для прогнозного значения
- •Пример 3.1.Исследовать зависимость объема прибыли от количества торговых точек. Сделать прогноз в предположении, что количество торговых точек будет увеличено до 25.
- •Вспомогательная расчетная таблица
- •Пример 3.2.Исследовать зависимость показателяуи факторахс помощью логарифмической, степенной и полиномиальной регрессий.
- •3.1.3. Модели рядов динамики
- •3.2. Автокорреляция данных и остатков
- •3.2.1. Автокорреляция данных
- •Пример 3.4. Исследовать на автокорреляцию динамический ряд:
- •Вспомогательная таблица для расчета коэффициента автокорреляции
- •3.2.2. Автокорреляция остатков
- •Причины возникновения автокорреляции
- •Вспомогательная таблица для расчета d-статистики
- •С помощью формулы (3.2.2) найдем d -статистику:
- •3.3. Мультиколлинеарность
- •Причины возникновения мультиколлинеарности:
- •Методы исследования мультиколлинеарности
- •Меры по устранению мультиколлинеарности:
- •3.4. Множественная линейная регрессия
- •3.4.1. Построение множественной линейной регрессии
- •Расчет элементов коэффициента
- •3.4.2. Матричный подход
- •Построение корреляционной матрицы
- •Построение модели множественной линейной регрессии и ее анализ
- •3.4.4. Нелинейные модели
- •3.4.5. Эластичность
- •4.Экспертные оценки и элементы теории графов
- •4.1. Ранговая корреляция
- •4.1.1. Экспертное оценивание
- •4.1.2. Этапы работ в системе экспертных оценок
- •4.1.3. Метод ранговой корреляции
- •Вспомогательные расчеты
- •Б) Случай многих экспертов
- •4.2. Элементы сетевого планирования
- •Основные элементы сетевого графика
- •Основные требования к сетевой модели
- •5. Индивидуальные задания для самостоятельной работы студентов по курсу “математическое моделирование”
- •5.1. Задания к разделу “Построение законов распределения”
- •5.2. Задания к разделу “Математическое программирование”
- •5.3. Задания к разделу “Регрессионный анализ”
- •Задание 2.
- •Задание 3.
- •5.4. Задания к разделу “Экспертные оценки и элементы теории графов” Задание 1.
- •Значение критерия Пирсона
- •Критерий Колмогорова
- •Критерий Колмогорова
- •Квантили распределения Стьюдента
- •Коефициентов автокорреляции
- •Литература
- •Пеніна Галина Геннадіївна, канд. Екон. Наук, доцент
1.2. Статистические критерии согласия
В результате наблюдений получено
эмпирическое распределение величины
Хс частотами.
В результате выбора закона определены
теоретические частоты
(k– число значений
или частичных интервалов).
а) Критерий согласия Пирсона
За меру расхождения частот принимают величину:
(1.2.1)
По уровню значимости
=
0,05 (0,1; 0,01) и числу степеней свободы
в таблице Пирсона (приложение А) находят
.
Если
,
то гипотеза о законе не отвергается, он
может быть использован. Если же
,
то гипотеза о рассматриваемом законе
распределения отвергается.
б) Критерий согласия Колмогорова
По вариационному ряду составляют
эмпирическую функцию распределения
.
Выбрав закон, формируют теоретическую
функцию
,
определяют максимальное отличие
и получают характеристику Колмогорова
.
(1.2.2)
Пользуясь специальными таблицами
(приложение Б), находят вероятность
.
Если эта вероятность меньше 0,01, то
гипотезу о выбранном законе распределения
отвергают. Если вероятность больше (или
равна) 0,01, то расхождения между эмпирической
и теоретической функциями признают
несущественными, а гипотезу о выбранном
законе распределения вполне согласованной
с экспериментом.
Замечание.Если0,29, то вероятность равна единице.
в) Критерий согласия Ястремского
Известный статистик Б.С.Ястремский доказал, что меру близости теоретического и фактического распределений можно характеризовать величиной
,
(1.2.3)
где
,
,
– теоретическая вероятность того, что
случайная величинаХ
примет значение
;k
– число групп;
при
.
Если
, то расхождение между теоретическим и фактическим распределениями несущественно.
Если
, то это расхождение существенно и его невозможно объяснить влиянием случайных факторов, поэтому теоретический закон распределения следует отклонить.
Г) Критерий согласия Романовского
Этот критерий используется для оценки степени приближения эмпирического распределения к теоретическому. Он тесным образом связан с критерием согласия Пирсона и состоит в том, что вычисляют величину
,
(1.2.3)
где
,q
– число степеней свободы,
.
Если
,
то это означает, что результаты испытаний
не противоречат гипотезе о законе
распределения. В противном случае
гипотеза отвергается.
Критерий В.И.Романовского является хорошим дополнением к критерию Пирсона.
1.3. Построение закона Пуассона
Этот закон используется во многих процессах, и в частности, в описании потока клиентов (заявок, требований) в системах массового обслуживания. Он описывает дискретные случайные величины, т.е. величины, которые могут принимать отдельные изолированные значения из некоторого интервала.
Предположим, что в результате наблюдений за случайной величиной получен дискретный вариационный ряд:
|
|
|
|
… |
|
|
|
|
|
… |
|
где
– значения случайной величины,
– частота соответствующих
,
.
Если предположить, что случайная величина, подчиняется закону Пуассона, то вероятности каждого значения должны определяться по формуле
,
(1.3.1)
где
– параметр закона, совпадающий с
математическим ожиданием:
.
Другими словами, чтобы сформировать закон, необходимо вычислить по опытным значениям математическое ожидание и подставить в формулу (1.3.1).
Определяя
для каждого значения
вероятность
,
находят соответствующие теоретические
частоты:
.
Пример 1.1. Было проведено 200 наблюдений (каждое длилось 2 мин.), в результате отмечалось следующее распределение покупателей:
-
0
1
2
3
4
5
6
7
41
62
45
22
16
8
4
2
Проверить, можно ли описать этот поток с помощью закона Пуассона (является ли он простейшим)?
Решение. Определим среднее число заявок (клиентов):
.
Сформируем функцию вероятностей Пуассона по формуле (1.3.1):
.
а) Проверим по критерию
Пирсона, подходит ли этот закон. Для
этого каждому значению
поставим в соответствие вероятность
,
определим
и рассчитаем
:
,
,
и т.д.
Все расчеты проведем в таблице 1.1.
Таблица 1.1
Проверка закона Пуассона по критерию Пирсона
|
|
|
|
|
|
0 |
41 |
0,1653 |
33 |
64 |
1,939 |
1 |
62 |
0,2975 |
60 |
4 |
0,067 |
2 |
45 |
0,2678 |
54 |
81 |
1,500 |
3 |
22 |
0,1607 |
32 |
100 |
3,125 |
4 |
16 |
0,0723 |
14 |
4 |
0,286 |
5 |
8 |
0,0260 |
5 |
9 |
1,800 |
6 |
4 |
0,0078 |
2 |
4 |
2,000 |
7 |
2 |
0,0026 |
0 |
4 |
2,000 |
|
200 |
1,0000 |
- |
- |
12,72 |
Замечание. |
Последнее
значение
|
Имеем,
,
.
По таблицам Пирсона (приложение А)
находим
.
В результате
.
Вывод: расчетное и критическое
значения достаточно близки, т.к. расчетное
значениевсе же больше критического, то возможность
применения закона сомнительна.
б) Проверим по критерию Колмогорова.
Предварительно найдем фактические ()
и теоретические вероятности (
),
затем сформируем эмпирическую (
)
и теоретическую (F)
функции распределения (таблица 1.2) и
найдем их разности.
Таблица 1.2
Проверка закона Пуассона по критерию Колмогорова
|
|
|
|
|
|
0 |
0,205 |
0 |
0,165 |
0 |
0 |
1 |
0,31 |
0,205 |
0,298 |
0,165 |
0,04 |
2 |
0,225 |
0,515 |
0,268 |
0,463 |
0,052 |
3 |
0,11 |
0,74 |
0,161 |
0,731 |
0,009 |
4 |
0,08 |
0,85 |
0,072 |
0,892 |
0,042 |
5 |
0,04 |
0,93 |
0,026 |
0,964 |
0,034 |
6 |
0,02 |
0,97 |
0,007 |
0,99 |
0,02 |
7 |
0,01 |
0,99 |
0,003 |
0,997 |
0,007 |
|
- |
1 |
- |
1 |
- |
По формуле (1.2.2) рассчитаем характеристику Колмогорова
.
Вывод:вероятность=1,
по данному критерию использование
закона вполне допустимо.
в) Оценим закон Пуассона по критерию Ястремского. Для нахождения величины Q, для чего построим вспомогательную расчетную таблицу 1.3.
Таблица 1.3
Вспомогательная расчетная таблица
|
|
|
|
|
41 |
33 |
0,165 |
0,835 |
2,323 |
62 |
60 |
0,298 |
0,702 |
0,095 |
45 |
54 |
0,268 |
0,732 |
2,049 |
22 |
32 |
0,161 |
0,839 |
3,725 |
16 |
14 |
0,072 |
0,928 |
0,308 |
8 |
5 |
0,026 |
0,974 |
1,848 |
4 |
2 |
0,007 |
0,993 |
2,014 |
2 |
0 |
0,003 |
0,997 |
2,006 |
|
|
|
|
14,368 |
Таким образом, Q=14,368. Посколькуk= 8,Z= 0,6, то по формуле (1.2.3) имеем:
А =.
Вывод:посколькуА< 3, то по этому критерию гипотеза о законе принимается.
г) Проверим пригодность закона Пуассона по критерию согласия Романовского.
Поскольку
= 12,72,q = 8 – 2 = 6, то
по формуле (1.2.3) имеем
.
Вывод:величинаменьше 3 – значит, в соответствии с этим
критерием гипотеза о законе Пуассона
не отвергается.
Замечание. |
В данном примере оказалось, что по трем критериям согласия из четырех закон Пуассона не отвергается. Следовательно, общий вывод состоит в возможности применения закона Пуассона к данному ряду. |