Обчислюємо значення
. (50)
Значення для кожного x має бути більше від 0.01.
Рекомендована довжина послідовності: .
Примітка . На справді, стан x має приймати значення від до . Розробники керівництва НІСТ, вибрали іншу формулу, виграли у простоті, але програли (хоча зовсім не багато) у точності.
3.16. Перевірка випадкових відхилень-2 (Random Excursions Variant Test)
Даний тест є розширенням попереднього.
Нехай - двійкова послідовність довжини n. Перетворимо її у послідовність Х: (тобто 1®1, 0 ® -1). Cформуємо послідовність:
. (51)
Cформуємо послідовність:
(52)
Розіб’ємо її на блоки, що мають вигляд , тобто на послідовності, у яких перший та останній елементи дорівнюють нулю, а всі інші відмінні від нуля.
Будемо вважати, що стан будь якого не нульового елементу х блоку може приймати одно із значень : (-9), (-8), (-7),...8, 9, тобто .
Нехай J – загальна кількість блоків та - кількість появ станів х у всіх J блоках.
Обчислюємо значення для кожного х
. (53)
Значення для кожного x має бути більше від 0.01.
Рекомендована довжина послідовності: .
Примітка . Як і у попередньому тесті розробники керівництва НІСТ спростили значення x. На справді, стан x має приймати значення від до .
4. Математичне обґрунтування тестів
4.1.Частотний тест
Основною ідеєю тесту є використання нульової гіпотези: у послідовності незалежних випадкових величин з розподілом Бернуллі (Х чи ε, де Х = 2ε – 1, так що Sn = X1 + … + Xn = 2(ε1 + … + εn) - n) імовірність появи 1 складає ½. Відповідно до класичної теореми Лапласа, для досить великого числа іспитів розподіл біноміальної суми, нормалізованої по , близько наближений до стандартного нормального розподілу. Даний тест використовує апроксимацію оцінки близькості частки 1 до ½. Результати всіх наступних тестів спираються на результати проходження цього основного тесту.
Даний тест отриманий з центральної граничної теореми для випадкового блукання, Sn = X1 + … + Xn... Відповідно до центральної граничної теореми,
(54)
Цей класичний результат є основою для найпростішого тесту на випадковість. Це припускає, що для позитивного z
Відповідно до тесту, заснованому на статистиці , обчислюється розглянута величина |s(obs)| = |X1 + … + Xn| /, після чого рахується відповідна Р-величина, рівна . Тут erfc є (додатковою) функцією помилки
Приклад.
Вхід: e = 110010010000111111011010101000100010000101101000110
0001000110100110001001100011001100010100010111000,
n = 100.
Тест:
s100 = n1 – n0 = 42 – 58 = -16,
,
- тест пройдений.
4.2.Частотний тест в середині блоку
Метою тесту є виявлення відхилень частки «1» від ½ усередині під послідовностей, що не перекриваються, на які розбивається вхідна послідовність. Для оцінки відхилення використовується критерій .
Малі P - значення вказують на великі відхилення від рівної пропорції одиниць і нулів принаймні в одному з блоків. Рядок 0 і 1 (чи 1 і -1) розділений на множину непересічних блоків. Для кожного блоку обчислюється пропорція одиниць. Статистика хі-квадрат порівнює пропорції блоку з ідеалом ½. Статистика заснована на розподілі хі-квадрат зі ступенями свободи, рівними кількості блоків.
Параметрами цього тесту є М і N, так що n = МС, тобто, послідовність розбивається на N блоків, кожен довжиною М. Для кожного з цих блоків імовірність появи «1» оцінюється шляхом спостереження відповідних частот «1», , і = 1, ..., N. Сума
у припущенні гіпотези випадковості повинна мати розподіл з N ступенями свободи. Р-значення розраховується відповідно до вираження
Приклад.
Вхід: e = 110010010000111111011010101000100010000101101000110
0001000110100110001001100011001100010100010111000,
n=100,
M=10.
Тест:
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.