
- •Министерство науки и образования украины
- •Инструктивно- методические указания
- •«Проверка качества последовательностей случайных чисел на равномерность»
- •Практическое занятие 11
- •План работы
- •1. Краткие теоретические сведения
- •2. Проверка последовательности случайных чисел на равномерность.
- •2.1 Проверка по гистограмме
- •2.2 Проверка с использованием критерия согласия Пирсона
- •2.3 Проверка по косвенным признакам
- •3. Проверка качества псч с помощью эвм в среде MatLab .
- •3. 1 Проверка по гистограмме
- •3.2 Проверка с использованием критерия согласия Пирсона
- •3.3 Проверка по косвенным признакам
2.2 Проверка с использованием критерия согласия Пирсона
Пусть х1, …, xn — выборка объемаn, представляющая собой результатnнезависимых наблюдений над случайной величинойX, относительно распределения которой выдвинута простая гипотезаН0 : FX(х) = F(х). (F(х) — теоретическая функция распределения, соответствующая гипотезеН0).
Наиболее распространенным критерием проверки этой гипотезы Н0 являетсякритерии 2 Пирсона.
Чтобы воспользоваться критерием 2 Пирсона, выборочные данныех1, …,xnследует предварительно сгруппировать, представив их в видеинтервального статистического ряда.
Пусть
—интервалы группировки;
1,…,N—частоты попадания выборочных значений в интервалыJ1,...,JNсоответственно(1 + ... + N = n).
Обозначим рk теоретическую (соответствующуюН0) вероятность попадания случайной величиныХ в интервал
.
Статистикой критерия 2является величина
.
Она характеризует отклонение эмпирической функции распределения
Fn*(х)(k/n—приращениеFn*(х)на интервалеJk)
от теоретической функции распределения F(х)(рk — приращениеF(х) на том же интервалеJk).
Поскольку
относительные частоты k/nсближаются с вероятностямирkприn, , то в случае справедливостиН0
значение величиныn2не должно существенно отличаться от
нуля.
Поэтому критическая область критерия 2задается в виде
К={tt},
где t = n2(х1, …,xn)— значение величиныn2, вычисленное для заданной выборки, а порогtопределяется по заданному уровню значимоститак, чтобыР{n2K/H0}=.
Нахождение tосновано на том факте (известном кактеорема Пирсона), что случайная величинаn2имеет приnпредельное распределениехи - квадратс (N—1) степенью свободы2(N—1).
На практике
предельное распределение 2(N—1)можно использовать с хорошим приближением
приn50иk5, .
При выполнении этих условий для заданного уровня значимости можно положитьt=21-,N-1, где21-,N-1— (1—)-квантиль распределения 2(N—1).
Таким образом, критерий согласия 2Пирсона состоит в следующем:
По заданному уровню значимости находится по таблице порог21-,N-1.
По заданной выборке х1, …,xnобъемаn 50определяется числоN интервалов группировки так, чтобыk5,
.
Вычисляется значение статистики 2n (х1, …,xn)=t.
Если t21-,N-1, тогипотезу Н0 отвергают.
Если t<21-,N-1, тогипотезу Н0 не отвергают.
Если случайная
величина X
дискретная,хk, — различные выборочные значения, аР{Х=хk}=рkв случае справедливостиН0,
то всегда можно определитьNинтервалов, содержащих ровно по одному
выборочному значению. Поэтому в данном
случае можно сразу считать, чтоk=mk,
,
гдеmk—частота
выборочного значенияхk.
На практике теоретическое распределение полностью определено редко. Чаще бывает известен предположительно только тип распределения, но неизвестны параметры его определяющие. В этом случае гипотеза о виде распределения, подлежащая проверке, имеет вид Н0: FX(х)=F(х;), =(1,...,r)и является сложной параметрической гипотезой.
Критерий согласия 2 Пирсонаприменим для проверки такой гипотезыН0 со следующими изменениями:
а) вероятности
рk=Р{XJk},вычисляют, заменяя неизвестные параметры
их оценками максимального правдоподобия:
б) число степеней свободыпредельного распределенияхи - квадратдолжно быть уменьшено на число неизвестных параметров и считаться равным (N—1—r).