Скачать

Исследование статистических характеристик случайной последовательности

Министерство Образования Республики Таджикистан

Таджикский Технический Университет

имени М.С. Осими


Кафедра «АСОИиУ»


Лабораторная работа №1


На тему: Исследование статистических характеристик

случайной последовательности.


Душанбе-2010


Лабораторная работа №1. Исследование статистических характеристик случайной последовательности

Цель работы:

1.Освоение методов оценки закона распределения и вероятностных характеристик случайной последовательности: математического ожидания, дисперсии, среднеквадратичного отклонения и автокорреляционной функции.

2.Освоение метода проверки гипотезы о законе распределения по критерию согласия хи- квадрат Пирсона.

3.Исследование свойств базовой псевдослучайной последовательности.

Теоретические сведения.

Оценка вероятностных характеристик.

Числовая последовательность Х1,Х2,...Хn, статистические характеристики которой требуется определить считается реализацией стационарной эргодической случайной последовательности Х1,Х2,...хn. Вероятностные характеристики случайной последовательности неизвестны и подлежат оценке с помощью соответствующих статистических характеристик числовой последовательности. При вероятностном моделировании последовательности Х1,Х2,...хn представляет собой совокупность результатов отдельных опытов. В данной лабораторной работе в качестве такой последовательности Х1,Х2,...хn рассматриваются псевдослучайные числа вырабатываемые генератором, построенным на М – последовательности (датчиком случайных чисел) в котором

M= g ⁿ -1 (1)

Где М – общее количество чисел, вырабатываемых генератором

g-основание системы исчисления

n- Количество разрядов в генераторе.

Генератор строится на базе регистра Хi(i=1,n),состоящего из ячеек , в которые записываются целые числа от 1 до g. Случайные числа М- последовательности снимаются с последнего элемента Хn. Числа записанные в ячейки Xm и Xn складываются по модулю g

R= Xm + Xn (2)

И приводится сдвиг чисел в регистре:

Xn-i= Xn-j-I (i=0,..n-2) (3)

В первую ячейку записывается содержимое сумматора Xi=R

Такая процедура повторяется М – раз, в процессе которой получается исследуемая базовая псевдослучайная последовательность Х1,Х2,...хn, где M=N. Для данной последовательности рассчитываются ее вероятностные характеристики.

Математическое ожидание M(Xi)=m оценивается по формуле:

m* =1/N ∑ Xi (5)

Дисперсия Dx оценивается по формуле:

Dx= 1/n-1∑(xi-mx)² (6)

Среднеквадратическое отклонение оценивается по формуле:

(м*=1/n∑xi) δ*= √D* (7)

i=1

Aавтокорреляционная функция (нормированная) представляет собой последовательность коэффициентов корреляции, зависящих от величины сдвига, как от аргумента.

K(r)=1/D · M((xi - m)(xi + r-m))

Ее оценка вычисляется:

K*(r)=1/D*(N-r-1)n-2i=1((xi-m*)(xi+r-m*))=1/D*(1/N-r-1)n-2i=1xixi+r-(N-r)/(N-r-1)m* (8)

Оценка закона распределения.

Одномерный закон распределения при большом объеме последовательности оценивается статистическим рядом, графическое изображение которого называется гистограммой. При малом объеме последовательности, когда N не превосходит несколько десятков, используется статистическая функция распределения, называемая также выборочной и эмпирической.

Для построения гистограммы диапазон возможных значений элементов последовательности разбивается на е участков точками U1,U2,Ue

Описание: Drawing1

Крайние точки Uo и Ue могут быть бесполезными. Длины участков ΔU могут быть необязательно одинаковыми. Если они различны, то чаще всего называются так, чтобы вероятности попадания на все участки были одинаково близки друг к другу. В связи с тем, что моделируемый генератор вырабатывает целые случайные числа от 1 до g , то участки выделяются точками U1=1;U2=2;Ue=g.

Статистический ряд- это совокупность чисел V1,V2,Ve , где Vj0- количество элементов последовательности, удовлетворяющее неравенствуUj-1 < Xi < Uj т.е попавших в j –участок. Графическое представление статистического ряда, т. Е гистограмму, удобно строить в относительных величинах. Поэтому производится нормировка:

e

∑ Vj / N=1 (9)

i=1

Статистическая (выборочная , эмпирическая) функция распределения F*(X) является оценкой для интегральной функции распределения и вычисляется по формуле :

0,если X

F*(X)= k/n, если Хk<Х<Хk+1

1, если X>Xn (k=1,2,..N-1) (10)

Где Xk-тый элемент вариационного ряда, т.е. последовательности, в которой элементы расположены в порядке возрастания числовых значений. Графическое представление функции распределения показано показано на рис.

Описание: Drawing2


Проверка гипотезы о законе распределения.

Гипотеза о законе распределения элементов последовательности задается названием закона и численным значением параметров. Она может быть задана плотностью вероятности в виде формулы или графика.

Иногда может быть задана интегральная функция распределения. Тогда знак F(x)можно всегда найти плотность вероятности как f(x)=p(x).

Для проверки гипотезы о законе распределения при большом объеме последовательности (n>100)пользуются критерием X² Пирсона. По построенному статистическому ряду (гистограмме) вычисляется статис. х² (Δ)

e

Δ= X² = ∑(Vj-NPj)² /(NPj) (11)

i=1

Где Pj- вероятность попадания элемента последовательности в j-ый участок

Vj- j-ый член стат.ряда, т.е. количество элементов последовательности попавших в j-ый участок

N – общее количество элементов последовательности.

Распределение Х² зависит от параметра r , называемого числом «степени свободы». Число степеней свободы r равно числу участков е минус число независимых условий, наложенных на частоты Pj =Vj / n (j=1,e). Примером такого условия может быть условие вида (9), которое накладывается при любом случае. Поэтому

R=e-1 (12)

Если для теоретического распределения задаются математическое ожидание, дисперсия и другие параметры, то число степеней свободы уменьшается на число таких параметров.

Для распределения Х² имеются специальные таблицы, по которым можно для каждого значения Х² и числа степеней свободы r найти вероятность P того, что величина, распределенная по закону Х² превзойдет его значение. Вероятность P, определенная по таблице, есть вероятность того, что за счет числа случайных причин мера расхождения теоретического и статистического распределения (11)будет не меньше, чем фактически наблюденное в данном серии опытов значения Х². если эта вероятность P весьма мала, то результат опыта следует считать противоречивым гипотезе о том, что закон распределения величины Х есть F(x). Поэтому эту гипотезу следует отбросить как неправдоподобную. Напротив, если вероятность P сравнительно велика, то можно признать расхождения между теоретическим и статистическим распределением вещественным. При этом гипотеза о том, что величина X распределена по закону F(x) можно считать правдоподобной или не противоречащей опытным данным.

На практике, если P оказывается меньше, чем 0,1, то рекомендуется проверить и по возможности повторить эксперимент. В случае, если опять появятся замеченные расхождения, то следует подобрать более подходящий для описания стат. Данных закона распределения.

Содержание исследования

В состав исследования, проводимого в данной лабораторной работе входит:

1. программная реализация базой псевдослучайной последовательности, вырабатываемой генератором случайных чисел при заданных преподавателем параметрах: g,n,m.

2.Определение оценок математич. ожидания, дисперсии, среднеквадратического отклонения и g коэффициентов корреляции (для r=1,g).

3. Построение гистограммы распределения.

4.Вычисление статист. по критерию Х² Пирсона.

5.Проверка гипотезы о равномерном распределении чисел от 0,1,2,..g, полученных генератором случайных чисел.

Таблица1

Критические точки распределения

Число степеней свободы rУровень значимости ά
0.010.0250.050.950.9750.99
16.65.03.80.00390.000890.00016
29.27,46,00,1030.0510.20
311.39,47,80,3520,2160,115
413,311,19,50,7110,4840,297
515,112,811,10,1150,8310,554
616,814,412,60,1641,240,872
718,516,014,10,2171,691,24
820,117,515,52,732,181,65
921,719,016,93,332,702,09
1023,220,518,33,943,212,56
1124,721,919,74,573,823,05
1226,223,320,05,234,403,57
1327,724,722,45,895,014,11
1429,126,123,76,575,634,66
1530,627,525,07,266,265,23
1632,028,826,37,966,915,81
1733,430,227,68,677,566,41
1834,831,528,99,398,237,01
1936,232,930,110,18,917,63
2037,634,231,410,99,598,26
2138,935,532,711,610,38,90
2240,336,833,912,311,09,54
2341,638,135,213,111,710,2
2443,039,436,413,812,410,9
2544,340,637,714,613,111,5
2645,641,938,915,413,812,2
2747,043,240,116,214,612,9
2848,344,541,316,915,313,6
2949,645,742,617,716,0214,3
3050,947,043,818,516,815,0

Таблица 2

Варианты заданий к лабораторной работе

qmnx1x2x3x4x5x6x7
133501211--
22371111111
32571111000
432511111--
52471110011
633500011--
7336012120-
8346111222-
9356111002-
102371010101
112671010101
122270101020
133351111000
142261111111
153451210122

Интерфейс программы


Введённые значения: q-2, M-4, N-7

X1-1, X2-1, X3-1, X4-1, X5-0, X6-0, X7-1

Листингпрограммы

Public m, n, q, r, xe As Integer

Private Sub Command2_Click()

End

End Sub

Private Sub ok_Click()

list.Clear

If tm.Text = "" Or tn.Text = "" Or tq.Text = "" Or tx(0).Text = "" Or tx(1).Text = "" Or tx(2).Text = "" Or tx(3).Text = "" Or tx(4).Text = "" Then

MsgBox ("Сначала введите все значения")

Exit Sub

Else

Call fun

End If

End Sub

Private Sub fun()

Dim reg, regt As String

Dim xi(6), yi(100000), p(2), p0, p1, p2 As Integer

Dim xe As Double

n = CInt(tn.Text)

m = CInt(tm.Text)

q = CInt(tq.Text)

For i = 0 To (n - 1)

reg = reg & CStr(tx(i).Text)

Next

mg = q ^ n - 1

list.AddItem ("M=q^n-1 = " & mg)

list.AddItem ("")

For i = 1 To mg

yi(i) = Right(reg, 1)

r = CInt(Mid(reg, m, 1)) + CInt(Right(reg, 1))

If r >= q Then

r = r - q

Else

End If

If r = 0 Then p0 = p0 + 1

If r = 1 Then p1 = p1 + 1

If r = 2 Then p2 = p2 + 1

regt = reg

reg = CStr(r) & Mid(regt, 1, (n - 1))

list.AddItem ("Генератор=" & reg)

list.AddItem ("Число=" & yi(i))

list.AddItem ("R=" & r)

Next

list.AddItem ("")

p(0) = p0 / mg

p(1) = p1 / mg

p(2) = p2 / mg

For w = 0 To q - 1

list.AddItem ("p" & w & "=" & p(w))

Next w

For j = 0 To q - 1

xe = xe + ((1 / q) - p(j)) ^ 2

Next

xe = (1 / q) * xe

list.AddItem ("")

list.AddItem ("Критерий Пирсона=" & Round(xe, 12))

Dim mx, dx, kr, k As Double

mx = 0

For j = 1 To mg

mx = mx + CInt(yi(j))

Next

mx = mx / mg

list.AddItem ("")

list.AddItem ("Математическое ожидание=" & mx)

For j = 1 To mg

dx = (CInt(yi(i)) - mx) ^ 2

Next

dx = dx / (n - 1)

list.AddItem ("")

list.AddItem ("Дисперсия=" & dx)

list.AddItem ("")

list.AddItem ("Среднеквадратичное отклонение=" & Sqr(dx))

list.AddItem ("")

For i = 1 To q

For j = 1 To (mg - i)

k = k + ((CInt(yi(j)) - mx) * (CInt(yi(j + i)) - mx))

Next

kr = k / (mg - i - 1)

list.AddItem ("Автокорреляционная ф-ия (" & i & ")= " & Round(kr, 12))

Next

End Sub


Список используемой литературы

1. Ли И. Т., Лабораторный практикум. Имитационное моделирование экономических процессов, Душанбе 2008 год