
- •Вычислительная система как объект моделирования.
- •Виды моделирования.
- •Методы вероятностного моделирования.
- •Основные этапы моделирования.
- •Моделирование случайных воздействий.
- •Случайные величины и их распределение
- •Мода случайной величины.
- •Медиана случайной величины.
- •Характеристики рассеивания случайной величины.
- •Корреляционный момент (корреляция).
- •Коэффициент корреляции:
- •Коэффициент автокоррелиации.
- •Правило трех сигма.
- •Экспоненциальное распределение.
- •Распределение Хи – квадрат.
- •Статистическое оценивание последовательности чисел.
- •Проверка качества последовательностей псевдослучайных чисел.
- •Соответствие по критерию согласия.
- •Критерий Хи – квадрат.
- •Критерий Колмогорова.
- •Проверка равномерности.
- •Тест пар.
- •Тест комбинаций (аналог теста пар).
- •Тест наибольшей из t.
- •Тест монотонности
- •Тест апериодичности, только для псевдослучайных чисел.
- •Моделирование случайных воздействий
- •Программный метод получения псевдослучайных чисел.
- •1 . Метод середины квадрата
- •2. Метод умножения.
- •3. Конгруэнтные процедуры.
- •Использование последовательностей равномерно распределенных случайных (псевдослучайных) чисел в задачах статистического моделирования.
- •Формирование последовательности случайных чисел имеющих неравномерный закон распределения.
- •Геометрическая интерпретация метода обратной функции
- •Метод обратной функции
- •Метод обратной функции для экспоненциального закона.
- •Метод Неймана
- •Метод ступенчатой аппроксимации
- •Специализированный метод для реализации последовательности псч по нормальному закону распределения.
- •Экспоненциальный закон распределения
- •Планирование машинного эксперимента
- •Тактическое планирование эксперимента
- •Тактическое планирование эксперимента
- •Технические и программные средства моделирования.
- •Технические средства моделирования.
- •2.Цифровые вычислительные машины
- •3.Гибритные вычислительные комплексы
- •4.Специализированные эвм
- •Программные средства моделирования.
- •Пакет прикладных программ “statistica”:
- •Лабораторная работа №1.
- •Вербальное описание на основе q-схем.
Соответствие по критерию согласия.
Критерии объединены общей методологией необходимо, вычислить величину (критерий), называемый “мерой расхождения”. И по данной мере расхождения определяем вероятность согласования функций теоретической и статистической.
F*(x) и F(x)
f*(x) и f(x)
Критерий Хи – квадрат.
Исходные данные – последовательность псевдослучайных чисел.
f(x) – теоретическая функция плоскости распределения.
Алгоритм.
Минимум 5 интервалов.
построение гистограммы.
Необходимо чтобы в каждом интервале попало не менее 10 значений, иначе критерий работать не будет.
вычисляем меру расхождения.
,
k – количество интервалов
гистограммы.
U – имеет закон распределение 2 с (k – 1) степенями свободы – количество связей положенных на последовательность.
при известных значениях U, k находим функцию распределения 2.
Критерий Колмогорова.
1) построение
гистограммы статистической функции
распределения
.
2)вычисляем меру расхождения
3)функция распределения Колмогорова. По заданной мере расхождения вычисляем вероятность согласования теоретического и эмпирического законов распределения.
Проверка равномерности.
Большинство генераторов равномерного распределения псевдослучайных чисел используются в качестве базовых при реализации более сложных законов распределения. Поэтому, промерив указанные последовательности по указанным критериям, выполняются дополнительные процедуры проверки равномерности.
Пусть каждое число представлено в двоичном виде.
Теоретически:
- количество чисел
- количество разрядов
Общий подход для тестов равномерности.
Известно
теоретическая плотность распределения
0 и 1 в разрядах числа.
можем вычислить количество 0 и 1.
,
Колмогорова.
Тест пар.
О равномерной последовательности судят по равномерности появления отдельных 0 и1 в отдельных разрядах числа.
Вычисляется величина Q – вероятность заполнения значений "1" разрядов числа.
Теоретическое
для каждого разряда значение вероятности
На практике:
- матожидание
доверительная
вероятность
Каждый разряд анализируется. В случае если число ½ попадается в данный интервал, то единица распределена равномерно.
Чем старше разряд числа, тем тщательно его поверяют.
Тест комбинаций (аналог теста пар).
За исключением того. Что проверяем закономерность мантисс числа (проверяем строки).
Алгоритм.
в качестве элементов последовательности
-число
единиц в мантиссе
строится ряд
,
l – число проверочных
разрядов l
k.
Rk – количество чисел имеющих k единиц в мантиссе.
теоретический ряд распределения (по биноминальному закону распределения)
p(k, l) – биноминальный закон
Эмпирическая функция распределения.
Хи – квадрат.
Тест многомерной равномерности.
(x1, … , xn)
…………………
Так как числа распределены по равномерному закону, они все заключены в куб с единичными сторонами.
Вероятность
попадания каждой из точек в гиперкубик
равна:
Выбор числа гиперкубиков.
Число гиперкубиков должно быть как можно больше, но при тестировании по критерию Хи – квадрат число должно быть равно 10 20 (20 30), поэтому существует ограничение на число гиперкубиков.
,
k – количество интервалов,
на который разбит гиперкубик.
Сторона гиперкуба.
из этой формулы следует, что количество отрезков разделяющих гиперкубик ограничено сверху.
,
t – число групп.