
- •1.1 Краткий исторический обзор
- •1.2 Математические методы и моделирование экономических процессов
- •1.3 Этапы математического моделирования
- •1.4 Классификация математических моделей
- •Вопросы для самопроверки
- •Тема 2. Модели производства
- •2.1 Производственные функции
- •2.1.1 Понятие производственной функции одной переменной
- •2.1.3 Формальные свойства производственных функций
- •2.1.4 Характеристики производственной функции
- •2.2 Задача производителя
- •2.3 Учет налогов
- •2.4 Функции спроса на ресурсы
- •2.5 Модели ценообразования
- •Вопросы для самопроверки
- •Тема 3. Функция полезности
- •3.1. Множество благ
- •3.2. Функция полезности и ее свойства
- •3.3. Предельная полезность и предельная норма замещения благ
- •3.4. Оптимальный выбор благ потребителем
- •3.4.1. Модель задачи оптимального выбора
- •3.5. Взаимная задача к задаче оптимального выбора благ потребителем
- •Вопросы для самопроверки
- •Тема 4. Балансовые модели
- •4.2 Экономико-математическая модель межотраслевого баланса
- •4.3 Коэффициенты прямых и полных материальных затрат
- •4.4 Агрегирование показателей межотраслевого баланса
- •4.5 Анализ экономических показателей
- •4.5.1 Модель затрат труда
- •4.5.2 Модель фондоемкости продукции
- •4.6. Динамическая модель межотраслевого баланса
- •Вопросы для самопроверки
- •Тема 5. Моделирование финансовых операций
- •5.1. Наращение и дисконтирование
- •5.1.1 Проценты и процентные ставки
- •5.1.2 Наращение по простым процентам
- •5.1.3. Сложные проценты
- •5.1.4. Номинальная и эффективная ставки процентов
- •5.1.6. Учет инфляции при наращении процентов
- •5.1.7. Эквивалентность простых и сложных процентных ставок
- •5.1.8. Дисконтирование и наращение по учетной ставке
- •5.1.9. Наращение по учетной ставке
- •5.1.10. Сравнение методов наращения
- •5.1.11. Сравнение методов дисконтирования
- •5.2. Потоки платежей, ренты
- •5.2.1. Основные определения
- •Ренты бывают постоянные и переменные.
- •5.3. Наращенная сумма потока платежей
- •5.3.1. Наращенная сумма годовой ренты
- •5.3.2.Наращенная сумма годовой ренты с начислением процентов m раз в год
- •5.4. Современная величина потока платежей
- •5.4.1. Современная величина годовой ренты
- •5.4.2. Современная величина годовой ренты с начислением процентов m раз в год
- •5.4.3. Современная величина p – срочной ренты ( m = 1)
- •5.4.4. Современная величина p – срочной ренты при начислении процентов m раз в год
- •5.5 Доходность финансовой операции
- •5.5.1. Различные виды доходности операций
- •5.5.2. Учет налогов и инфляции
- •5.5.3. Поток платежей и его доходность
- •5.5.4. Мгновенная доходность
- •5.6. Кредитные расчеты
- •5.6.1. Показатель полной доходности финансово-кредитной операции
- •5.6.2. Баланс финансово-кредитной операции
- •5.6.3. Определение полной доходности ссудных операций с удержанием комиссионных
- •5.6.4. Методы сравнения и анализа коммерческих контрактов
- •5.6.5. Планирование погашения долгосрочной задолженности
- •Вопросы для самопроверки
- •Тема 6. Математическое и компьютерное моделирование
- •6.1. Классификация видов моделирования
- •6.2. Достоинства и недостатки имитационного моделирования
- •6.3. Типовые задачи имитационного моделирования
- •6.4. Социально-экономические процессы как объекты моделирования
- •6.5. Примеры задач имитационного моделирования
- •Вопросы для самопроверки
- •Тема 7. Сущность метода имитационного моделирования
- •7.1. Метод имитационного моделирования и его особенности
- •7.2. Процесс имитации
- •7.3. Формулирование модели
- •7.4. Оценка адекватности модели
- •7.5. Экспериментирование с использованием имитационной модели
- •7.6. Понятие о модельном времени. Механизм продвижения модельного времени
- •7.7. Интерпретация и реализация результатов моделирования
- •Вопросы для самопроверки
- •Тема 8. Имитационная модель глобальной системы
- •8.1. Основные компоненты динамической мировой модели
- •8.2. Концепция «петля обратной связи»
- •8.3. Основные петли «обратных связей» в мировой модели
- •8.4. Основные переменные в мировой модели
- •8.5. Структура модели мировой системы
- •8.6. Основные результаты экспериментов на модели мировой системы
- •Вопросы для самопроверки
- •Тема 9. Метод Монте-Карло и проверка статистических гипотез
- •Тема 10. Моделирование случайных событий
- •10.1. Моделирование простого события
- •10.2 Моделирование полной группы несовместных событий
- •10.3 Моделирование дискретной случайной величины
- •10.4 Моделирование непрерывных случайных величин
- •10.4.1. Метод обратной функции
- •10.4.2. Моделирование случайных величин с показательным распределением
- •10.4.3. Моделирование случайных величин с равномерным распределением на произвольном интервале (a, b)
- •10.4.4 Моделирование случайных величин с нормальным распределением
- •10.4.5. Моделирование случайных величин с усеченным нормальным распределением
- •10.4.6 Моделирование случайных величин с произвольным распределением
- •Вопросы для самопроверки
- •Тема 11. Системы массового обслуживания
- •11.1. Основные понятия. Классификация СМО
- •11.2 Понятие марковского случайного процесса
- •11.3 Потоки событий
- •11.4. Уравнения Колмогорова. Предельные вероятности состояний
- •11.5. Процесс гибели и размножения
- •11.6. CMО с отказами
- •11.7. СМО с ожиданием (очередью)
- •Вопросы для самопроверки
- •Тема 12. Модели управления запасами
- •12.1. Основные понятия
- •12.2. Статическая детерминированная модель без дефицита
- •12.3. Статическая детерминированная модель с дефицитом
- •12.4. Стохастические модели управления запасами
- •12.5. Стохастические модели управления запасами с фиксированным временем задержки поставок
- •Вопросы для самопроверки
- •ЛИТЕРАТУРА
140
моделировании дискретных случайных величин фактически используется та же процедура, что и при моделировании ПГНС.
10.4 Моделирование непрерывных случайных величин
10.4.1. Метод обратной функции |
|
|
|
|
Пусть имеется некоторая непрерывная случайная величина |
x , заданная функцией |
|||
распределения |
F(x) . Можно доказать, что значения этой |
функции |
равномерно |
|
распределены в интервале (0,1). Поэтому между случайной величиной |
z , |
равномерно |
||
распределенной в том же интервале, и функцией распределения случайной величины x |
||||
существует взаимно однозначное соответствие, т. е. |
|
|
|
|
|
z F(x) . |
|
|
(10.1) |
Отсюда следует, |
|
|
|
|
|
x F 1(z) . |
|
|
(10.2) |
Следовательно, если уравнение (2.1) имеет аналитическое решение, то для
моделирования случайной величины |
x |
можно использовать |
датчик случайных чисел, |
||||||||
генерирующий величину z , и затем осуществить расчет по формуле (2.2). |
|||||||||||
10.4.2. Моделирование случайных величин с показательным распределением |
|||||||||||
Пусть имеется случайная величина |
|
x |
с |
|
показательным |
распределением. Функция |
|||||
распределения имеет вид: |
|
|
|
|
|
|
|
|
|
|
|
|
F ( x) |
1 e |
x |
, |
|
|
|||||
|
|
|
|
|
|||||||
где – параметр распределения. |
|
|
|
|
|
|
|
|
|
|
|
Применив метод обратной функции, получим: |
|
|
|
|
|
||||||
z F (x) 1 e |
x |
, |
|
||||||||
|
|
|
|||||||||
откуда |
|
|
|
|
|
|
|
|
|
|
|
x |
1 |
ln(1 z) . |
|
|
|
(10.3) |
|||||
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
Учитывая, что случайная величина |
(1 z) имеет также равномерное распределение в |
||||||||||
интервале (0,1), соотношение (10.3) можно заменить соотношением |
|||||||||||
|
x |
1 |
ln(z) . |
|
|
|
|
||||
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
Экспоненциальный закон распределения применяют для моделирования следующих явлений:
времени поступления заказа на предприятие;
посещения покупателями магазина супер-маркета;
времени телефонных разговоров;
срока службы деталей и узлов в компьютере.
На первый взгляд, это распределение является надуманным. Поэтому рассмотрим предельную теорему о суперпозиции потоков. Предположим, что можно наблюдать k независимых потоков событий (рис. 10.4)

141
Рис. 10.4
В каждом таком потоке можно наблюдать m j элементарных событий ( j 1,...,k
Интервалы между tkj событиями — это независимые случайные величины,
).
распределенные по по неизвестному закону с математическим ожиданием j . Спроектируем моменты всех событий на общую ось времени и рассмотрим случайный
интервал времени
t T (k)
между двумя событиями полученного суммарного потока,
|
k |
|
j |
|
|
состоящего из n событий , где n |
|
m |
. |
|
|
|
|
|
|||
|
j 1 |
|
|
|
|
Теорема. Если сделать предельный переход и устремить |
n , то распределение |
||||
случайной величины интервала t T (k) |
в суммарном потоке событий, состоящем из |
элементарных потоков, устремится к экспоненциальному закону с математическим ожиданием
k
M (t)
k j
j 1
.
Следствие. Поток заявок, интервал поступления которых в некую систему имеет экспоненциальное распределение, является простейшим.
Пример 10.1. Допустим, что имеется некая крупная фирма. Клиенты фирмы - это физические и юридические лица. Каждый из них может иметь набор планов и расписанных дел на значительном интервале времени. Однако если рассмотреть суммарный поток обращений этих клиентов к служащим фирмы по разным вопросам, то интервал времени между двумя последовательными обращениями в соответствии с рассмотренной теоремой является случайной величиной, распределенной по экспоненциальному закону.
10.4.3. Моделирование случайных величин с равномерным распределением на произвольном интервале (a, b)

142
Датчик случайных чисел генерирует случайные величины с равномерным распределением в интервале (0,1). Если же нужно моделировать случайные величины с равномерным распределением в интервале (a,b) то можно воспользоваться методом обратной функции.
Для рассматриваемого случая выражение (2.1) примет вид:
z F (x) x a , b a
Откуда
x
a
z(b a)
.
На практике применяется и другой способ Вместо границ интервала задаются среднее
задания |
равномерного распределения. |
||
значение |
случайной величины |
xcp |
и |
величина интервала x . Тогда определение возможного значения случайной величины с равномерным распределением может быть произведено по формуле
x xcp
x(z
0,5)
.
Равномерное распределение можно использовать при расчетах по сетевым графикам работ, в военном деле (времени выдвижения воинской части или ее подразделения на исходный рубеж, времени марша)
10.4.4 Моделирование случайных величин с нормальным распределением
Метод обратной функции для нормального распределения неприменим, так как после подстановки соответствующей функции распределения выражение (2.2) не имеет аналитического решения. Поэтому в данном случае применяется другой метод.
Согласно центральной предельной теореме теории вероятностей при сложении достаточно большого числа одинаково распределенных независимых случайных чисел получается случайная величина, имеющая нормальное распределение.
Как показали исследования, уже при сложении более десяти случайных величин с равномерным распределением в интервале (0,1) получается случайная величина, которая с точностью, достаточной для большинства практических задач, может считаться распределенной нормально.
Процедура розыгрыша нормально распределенной случайной величины заключается в следующем.
1.Сложим 12 случайных величин с равномерным распределением в интервале (0,1), т. е. составим сумму
v
12
zi
i 1
.
Использовав известные теоремы о сумме математических ожиданий и дисперсий независимых случайных величин, можно установить, что в данном случае случайная величина v имеет следующие характеристики:
математическое ожидание:
12 |
|
|
1 |
|
|
M (v) M (zi ) 12 |
6 ; |
||||
|
|||||
i 1 |
2 |
|
|||
|
|
|
|
||
дисперсия: |
|
|
|
|
|
12 |
|
1 |
|
||
D(v) D(zi ) 12 |
1 . |
||||
12 |
|||||
i 1 |
|
||||
|
|
|
|

|
143 |
|
среднее квадратическое отклонение: |
(v) |
D(v) 1. |
2. Нормируем и центрируем случайную величину v , т. е. перейдем к величине
|
v M (v) |
v 6. |
|
(v) |
|||
|
|
2.От нормированной и центрированной величины перейдем к случайной
|
величине y с заданными параметрами M ( y) и ( y) по формуле |
|
|
|
|
y M ( y) ( y) , |
|
|
|
где |
M ( y) – известное математическое ожидание случайной величины |
y |
; ( y) |
– |
известное среднее квадратическое отклонение случайной величины |
y . |
|
|
Любые сложные работы на объектах экономики (ввод информации из документа в компьютер, проведение переговоров, ремонт оборудования и др.) состоят из многих коротких последовательных элементарных составляющих работ. Причем количество этих составляющих настолько велико, что условие выполнимости центральной предельной теоремы не вызывает сомнений. Поэтому при оценках трудозатрат всегда справедливо предположение о том, их продолжительность — это случайная величина, которая распределена по нормальному закону.
10.4.5. Моделирование случайных величин с усеченным нормальным распределением
|
Усеченное |
нормальное распределение |
случайной |
величины x |
задается четырьмя |
||||||||||||||||||||||
параметрами: математическим ожиданием M (x) , |
средним квадратическим отклонением |
||||||||||||||||||||||||||
(x) , а также минимальным и максимальным значениями |
x1 и |
x2 |
(точками усечения). |
||||||||||||||||||||||||
|
Функция распределения случайной величины |
x |
определяется равенством |
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0, |
|
|
|
|
|
|
|
|
|
x x |
; |
|
|
|
|
|
|
||||
|
|
|
|
F (x) |
|
|
(t) |
|
(t |
)] A, |
x |
|
x |
x |
|
; |
|
|
|||||||||
|
|
|
|
|
[ |
0 |
0 |
|
2 |
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
1 |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x x |
, |
|
|
|
|
|
|
||||
|
|
|
|
|
1, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
x M (x) |
|
|
|
|
x |
M (x) |
|
|
|
|
x M (x) |
|
|||||||||
где |
A |
|
|
; t |
|
|
|
|
|
; |
t |
|
1 |
|
|
|
|
|
; t |
|
|
|
|
2 |
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
0 (t2 ) 0 (t1) |
|
|
(x) |
|
|
1 |
|
|
(x) |
|
|
2 |
|
|
|
(x) |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Существуют также формулы для расчета математического ожидания, дисперсии и среднего квадратического отклонения случайной величины x . Однако с достаточной для практики точностью при моделировании случайной величины с усеченным нормальным распределением можно обойтись без расчетов по формулам.
Для определения возможных значений случайной величины с этим распределением можно использовать алгоритм, схема которого приведена на рис. 10.5.