- •1. Случайные числа и случайные цифры. Таблицы случайных цифр.
- •3. Датчики случайных чисел.
- •4. Сравнение трех способов с практической точки зрения.
- •5. Метод псевдослучайных чисел.
- •6. Стандартный датчик псевдо-сдучайных чисел реализованный на эвм
- •7. Простые алгоритмы. Длина отрезка апериодичности
- •9. Алгоритм д. Неймана
- •10. Алгоритм д. Леммера
- •11. Тесты для проверки случайных цифр
- •12. Проверка псевдослучайных чисел.
- •13. Моделирование дискретных случайных величин.
- •14. Оптимизация метода интервалов
- •15. Моделирование случайных событий.
- •Моделирование случайных непрерывных величин.
- •Моделирование многомерной случайной точки.
- •Поправки к приближенным распределениям.
- •Разделение области моделирования случайной величины.
- •Общая характеристика методов.
- •Метод Неймана.
- •Модифицированный метод Неймана.
- •Метод Метрополиса
- •Моделирование усеченных распределений.
- •Выбор равномерно распределенных точек в сложных пространственных областях.
- •Простейший метод вычислений.
- •Геометрический метод.
- •Сравнение точности методов Монте Карло.
- •Сравнение трудоемкости различных алгоритмов Монте Карло.
- •31. Моделирование процесса переноса иМетод имитации для решения задач о прохождении излучения через слой.
10. Алгоритм д. Леммера
Наибольшее распространение получил алгоритм, предложенный Д. Л е м е р о м, который называют методом вычетов или методом сравнений.
В этом методе Ф (х) = Д (gx), т.е. , (6) где g- большое целое число.
Удовлетворительная последовательность псевдослучайных чисел получается, например, при g=, M=, =1 (L=P=);
11. Тесты для проверки случайных цифр
Для проверки таблицы случайных цифр ε1, ε2, …, εN М. Г. Кендалл и Б. Б. Смит предложили использовать четыре теста. В каждом из них цифры классифицируются по некоторому признаку и эмпирические частоты сравниваются с их математическими ожиданиями при помощи критерия χ2. Тесты эти: проверка частот — проверяется частота различных цифр в таблице; 2) проверка пар — проверяется частота различных двузначных чисел среди пар ε1ε2, ε2ε3, ε3ε4,…, εN-1εN; 3) проверка интервалов — проверяется частота различных интервалов между двумя последовательными нулями; 4) проверка комбинаций — проверяется частота различных типов четверок (abcd, aabc, aabb, aaab, aaaa) среди четверок ε1ε2ε3ε4, ε2ε3ε4ε5, . . .
В большинстве случаев используется эта система тестов, внося, однако, в нее некоторые изменения, среди которых отметим два. Вместо проверки интервалов обычно используют проверку серий: цифры ε k+1, ε k+2, ..., ε k+l образуют серию длины l, если ε k+1= ε k+2= ...= ε k+l, но ε k≠ ε k+1, ε k+l+1≠ε k+l. Вместо проверки пар ε1ε2, ε2ε3, ε3ε4,…проверяют независимые пары ε1ε2, ε3ε4, ε5ε6,…
С детерминистической точки зрения проверка частот и проверка независимых пар — важнейшие необходимые тесты. Проверка серий обобщает «критерий случайности», часто используемый в статистике. И только проверка комбинаций носит несколько искусственный характер.
Положительный результат любого теста означает. Только, что этот результат не противоречит гипотезе о случайности цифр ε1, . . ., εN, но, может быть, какой-нибудь другой тест эту гипотезу опровергнет .
12. Проверка псевдослучайных чисел.
В качестве основных тестов для проверки псевдослучайных чисел γ1, . . ., γN используют те же тесты, что и для проверки таблиц: проверяются первые десятичные цифры εk = Ц (10 γk).
Более детальную проверку распределения чисел γ1, . . ., γN можно осуществить с помощью критерия . Если расположить эти числа в вариационный ряд
То
Однако, при очень больших N построение вариационного ряда весьма трудоемко. ,.
Естественно, что при проверке псевдослучайных чисел всегда используют различные дополнительные тесты для проверки последующих десятичных цифр γk. Иногда десятичные (или двоичные) цифры чисел γk выделяют и проверяют независимо. На практике ограничиваются простыми проверками, и далекие цифры чисел γk обычно оказываются хуже проверенными.
В пользу указанной системы тестов можно привести аргумент практического характера: в литературе, пожалуй, нет примеров, когда числа, удовлетворяющие всем тестам, оказались бы непригодными для решения конкретной задачи (в которой не предъявлялись повышенные требования к точности решения); есть, однако, примеры неудачных расчетов с помощью чисел, которые не удовлетворяли одному из тестов.
Таким образом, строго говоря, разумность приведенной системы тестов — факт эмпирический. В действительности эти тесты не гарантируют универсальной пригодности чисел. Поэтому иногда целесообразно вводить дополнительные тесты, связанные с характером решаемых задач. Впрочем, успешное решение нужной задачи— самая лучшая проверка случайных чисел.