Файл: Шаракшанэ, А. С. Испытания сложных систем учеб. пособие.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 21.10.2024

Просмотров: 130

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

лирования разбить на заранее выбранное число I интервалов. При выборе числа интервалов нужно руководствоваться общим положе­ нием, чтобы в каждоминтервале было не менее пяти значений случайной величины z. Если это условие не выполняется, то целе­ сообразно скорректировать априори выбранное число разбиений / путем слияния соседних интервалов так, чтобы число попавших точек в каждый интервал было больше или равно пяти [28];

2)рассчитать оценки математического ожидания mz* и диспе

сии аг*2

т

 

 

 

 

(7.4.36)

 

т

 

 

 

 

 

 

 

(7.4.37)

и с помощью их по таблицам распределения

функции

Лапласа

Ф(х)

найти вероятности pz появления значений случайной величи­

ны z

в каждом i-ш интервале (i =

l, 2, ...,

/). Если

гипотеза

Но: F(z) = N{mz, az2} верна, то оценки

(7.4.36)

и (7.4.37)

являются

асимптотически эффективными оценками;

 

 

3)

по данному уровню значимости р и найденному числу степе­

ней свободы k = l — 3 определить с помощью таблиц ^-распределе­ ния (табл. 7.4.9) границу %р критической области;

4) сравнить найденную границу %%со значением статистики,

рассчитанной по формуле (7.4.35). Если окажется,

что вычисленное

значение статистики %*2 меньше границы

, то проверяемую гипо­

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

 

гипотезу Но не

Вероятность ошибочно отвергнуть правильную

превосходит (3. Эту вероятность обычно называют

уровнем значи­

мости %2-критерия. Для %2-критерия характерна

некоторая неоп­

ределенность выбора числа и способа разбиения диапазона изме­ нения случайной величины z на интервалы. Поэтому при -малых объемах выборки к результатам проверки нужно относиться крити­ чески и рассматривать их с-учетом тех априорных сведений, кото­ рые обычно известны об изучаемых явлениях до моделирования.

В сомнительных случаях, когда %2-кр'Итерий принимает гипоте­ зу Н0 : F(z) =N {mz, о22} с невысоким уровнем значимости, а экспе­ риментатор убежден в обратном, целесообразно воспользоваться критериями, позволяющими оценить степень отклонения реального распределения от нормального. Эти критерии основаны на том, что

если случайная величина

z распределена

по нормальному зако­

ну N{mz, crz2}, то справедливы следующие соотношения [39]:

 

для нормированного среднего абсолютного отклонения

 

 

— =0,79788;

(7.4.383

,аг

«гг

"

 

V ПЛ

 

 

154


^

 

 

 

 

 

 

 

 

 

'

 

 

 

Таблица 7.4.9

#

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Критические значения хр

 

 

 

 

 

 

 

р

0,999

0 ,9 9

0 ,9 7 5

0 ,9 5

0,9 0

0,80

0 ,7 0

0,6 0

0,5 0

0 ,4 0

0 ,3 0

0,2 0

0,1 0

0,05

0 ,025

0,01

 

 

 

 

 

 

 

k

1

2

3

4

5

6

7

8 .

9

10

п

12

13

14

15

16

 

1

0,00000157

0,000157

0,000982

0,00393

0,0158 0,0642

0,148

0,275

0,455

0,708

1,07

1,64

2,71

3,84

5,02

6,63

2 0 ,0 0 2 0 0

0 , 0 2 0 1

0,0506

0,103

0 , 2 1 1

0,446

0,713

1 , 0 2

1,39

1,83

2,41

3,22

4,61

5,99

7,38

9,21

3

0,0243

0,115

0,216

0,352

0,584

1 , 0 0

1,42

1,87

2,37

2,95

.3,67

4,64

6,25

7,81

9,35

11,3

4

0| 0908

0,297

0,484

0,711

1,06

1,65

2,19

2,75

3,36

4,04

4,88

5,99

7,78

9,49

1 1 , 1

13,3

5

0 ' 2 1 0

0,334

0,831

1,15

1,61

2,34

3,00

3,66

4,35

5,13

6,06

7,20

9,24

1 1 , 1

1 2 , 8

15,1

6

0,381

0,872

1,24

1,64

2 , 2 0

3,07

3,83

4,57

5,35

6 , 2 1

7,23

8,56

1 0 , 6

1 2 , 6

14,4

16,8

7

0'598

1,24

1,69

2,17

2,83

3,82

4,67

5,49

6,35

7,28

8,38

9,80

1 2 , 0

14,0

16,0

18,5

8

0,857

1,65

2,18

2,73

3,49

4,59

5,53

6,42

7,34

8,35

9,52

1 1 , 0

13,4

15,5

17,5

2 0 , 1

9

1 ' 15

2 j09

2,70

3,33

4,17

5,38

6,39

7,36

8,34

9,41

10,7

1 2 , 2

14,7

16,9

19,0

21,7

10

1'48

2,56

3,25

3,94

4,87

6,18

7,27

8,30

9,34

10,5

1 1 , 8

13,4

16,0

18,3

10,5

23,2

12

2 , 2 1

3j57

4,40

5,23

6,30

7,81

9,03

1 0 , 2

11,3

1 2 , 6

14,0

15,8

18,5

2 1 , 0

23,3

26,2

14

3,04

4,66

5,63

6,57

7,79

9,47

1 0 , 8

1 2 , 1

13,3

14,7

16,2

18,2

2 1 , 1

23,7

26,1

29,1

16

3,94

5,81

6,91

7,96

9,31

1 1 , 2

1 2 , 6

14,0

15,3

16,8

18,4

20,5

23,5

26,3

28,8

32,0

18

4|90

7,01

8,23

9,39

10,9

12,9

14,4

15,9

17,3

18,9

2 0 , 6

2 2 , 8

26,0

28,9

31,5

34,8

20

5',92

8,26

9,59

10,9

12,4

14,6

16,3

17,8

19,3

2 1 , 0

2 2 , 8

25,0

28,4

31,4

34,2

37,6

25

8,65

11,5

13,1

14,6

16,5

18,9

20,9

2 2 , 6

24,3

26,1

28,2

30,7

34,4

37,7

40,6

44,3

30

11 ',6

15j0

1 6 18

18,5

2 0 , 6

23,Ф

25,5

27,4

29,3

31,3

33,5

36,3

40,3

43,8

47,0

50,9

35

14,7

18,5

2 0 , 6

22,5

24,8

27,8

30,2

32,3

34,3

36,5

38,9

41,8

46,1

49,8

53,2

57,3

40

17^9

2 2 ,2

24,4

26,5

29,1

32,2

39,9

37,1

39,3

41,6

44,2

47,3

51,8

55,8

59,3

63,7

45

21,3

25,9

28,4

30,6

33,4

36,9

39,6

42,0

44,3

46,8

49,5

52,7

57,5

61,7

65,4

70,0

50

24,7

29,7 •

32,4

34,8

37,7

41,4

44,3

46,9

49,3

51,9

54,7

58,2

63,2

67,5

71,4

76,2

60

31,7

37,5

40 ',5

43,2

46,5

50,6

53,8

56,6

59,3

62,1

65,2

69,0

74,4

79,1

83,3

88,4

70

39,0

45,4

48,8

51,7

55,3

59,9

63,3

66,4

69,3

72,4

75,7

79,7

85,5

90,5

95,0

100,4

80

46^5

53,5

57,2

60,4

64,3

69,2

72,9

76,2

79,3

82,6

36,1

90,4

96,6

101,9

106,6

112,3

90

54,2

61,0

65;6

69,1

73,3

78,6

82 /5

8 6 , 0

89,3

92,8

96,5

1 0 1 , 1

107,6

113,1

118,1

124,1

100

61,9

70,1

74,2

77,9

82,4

87,9

92,1

95,8

99,3

102,9

106,9

111,7

118,5

124,3

129,6

135,8


для коэффициента асимметрии

 

 

 

 

 

?! = —{1г—

т^ 3} = 0 ;

(7.4.39)

для коэффициента эксцесса

 

 

 

 

 

 

р2- 3

= £- ((г -

Отг)--1'----3 = 0 .

(7.4.40)

 

 

4

 

 

 

 

 

Чтобы узнать степень невыполнения равенств

(7.4.38) 4- (7.4.40),

нужно рассчитать статистики:

 

 

 

 

 

 

 

 

т

 

 

|

 

8 * —

1

V

I

 

-

 

 

ftl'Z J

 

 

*

 

 

 

т а

; = i

 

 

 

 

 

 

 

 

 

 

 

 

т

 

 

 

 

у * _

1

- V

i

 

-

* |3

 

 

f n Z 1 У

 

^

*3

 

 

7

а

' ~

1

 

 

m a z

 

1

 

 

 

 

 

 

т

 

 

 

 

 

В* —

1

_ V u

-

* и

 

U l z 1 у

 

^

*4

/

j

1

 

 

 

 

т * г

 

i

 

 

 

 

 

 

 

 

 

 

 

 

*

 

rn

 

 

 

 

 

 

 

 

 

 

 

 

mz~

m

 

 

 

 

 

 

 

/■=1

 

 

 

 

 

 

 

*2

1

m

 

 

*49

 

 

 

 

 

 

 

пг — 1

V

i

 

,

mz Г

 

 

Z

j

' ~

1

 

Когда статистики б*, Yi*. Рг* найдены, то с помощью таблиц, ко­ торые приведены в работе [39], нетрудно рассчитать количественные оценки, позволяющие определить, по каким нараметрам выбороч­ ное распределение наиболее сильно отличается от нормального.

Пример. Чтобы проиллюстрировать методику проверки нормальности распре­ деления результатов моделирования по некоторой выборке z u z2, , zm, примем в качестве ее элементов случайные значения х,-, которые приведены в табл. 7.2.4 (т. е. примем г ,= х i при t= 1,2, .. . , 50).

Следуя изложенной выше процедуре проверки:

1 ) определяем zmax=3,521; Zjnin= —2,526 и для 1=6 с учетом округления находим границы интервалов деления (— °о, —2], (—2 , —1], (—1 ,0], (0,1 J; (1, 2 ]; (2, + °э). Чтобы сумма вероятностей попадания в каждый интервал была равна единице, примем левой границей первого интервала — оо, а правой границей ше­ стого интервала +

2 ) находим частоты попадания случайной величины г,- в каждый интервал ( т.\=\\ т2'=Ь\ тз'=16; т 4'=18; т 5' = 8 ; пг6'= 3) Так как числа пц' и in6' мень­ ше пяти, то объединим первый интервал со вторым, а шестой интервал с пятым. Таким образом, границами уточненного разбиения интервала [zmax, zmin] будут

(—оо,—1], (—1,0], (0,1], (1, + оо), т. е. 1=4.

Для нового разбиения, получим т i=5; /7г2=16; т 3= 18; Ш4= 1 1 ;

156


3) по формулам (7.4.36) и (7.4.37) находим, что т г*= 0,21116; аг*=.1,160886

и с помощью таблицы распределения Ф(я) определяем вероятности попадания в каждый интервал (pi*=0,14917; р2* = 0,27941; р3*=0,323L7; рА* =0,24825);

4) вычислим выборочное значение %2— критерия по

формуле (7.4.35) : х*2=

= 1,47582;

 

 

 

 

2 для k—l—3=

5) сравнивая

значения %*2 с

критическими

значениями

= 4—3=1 и для

(3^0,2, получим,

что для всех

рассматриваемых

(3 справедливо

неравенство: х*2<Хр ■

 

 

 

 

Следовательно, анализируемую выборку результатов

моделирования можно

считать нормально распределенной, если априори выбираемый уровень значимо­ сти меньше 0,2.

Рассчитывая значения статистик б*= 0,63233675; у;* = 0,3797265; р2*=3,6929147

и сравнивая их с критическими значениями, приведенными в табл. 7.4.10, нетруд­ но убедиться в несущественности отклонения выборочного распределения от нор­ мального по крайней мере для общего уровня значимости, равного 0,05. Для ука­ занного уровня значимости критические значения соответственно равны бр=0.0б“

= 0,8644; Yi,р-=о,05 =0,533; Рг,р=о,05 =4,01.

Задача проверки гипотезы об однородности распределения двух выборок (критерий Уилкоксона). При решении задачи, связанной с проверкой гипотезы о принадлежности двух случайных выборок одной генеральной совокупности, можно воспользоваться критери­ ем Уилкоксона, который основан на применении теории инверсий. Суть этого критерия сводится к следующему.

Пусть выборка результатов хх, х2, ..., хп статистически независи­ ма с результатами z x, z2, ..., zm. Для этих условий проверка однород­ ности распределения двух выборок состоит в проверке адекватно­ сти законов распределения fi(x), F2(z), характеризующих их гене­ ральные совокупности. Законы распределения Fx(x), F2(z) априори неизвестны. При общей альтернативе Нх :FX=£F2 решение постав­ ленной задачи может быть найдено на основании критерия Смир­ нова, если для каждой выборки рассчитать эмпирические функции распределения. В частных случаях процедуры расчета эмпирических функций распределения можно избежать, воспользовавшись при проверке гипотезы Н0: Fх(х) =F2(z) порядковым критерием, каким и является критерий Уилкоксона. Состоятельность критерия Уил­ коксона доказана, когда в качестве альтернативы выступает одно­ сторонняя гипотеза о том, что величины х* стохастически больше или меньше величин Zj. Чтобы воспользоваться критерием Уилкаксойа, нужно:

1) построить общий вариационный ряд в порядке совместного возрастания случайных величин х*, z, (при i=l, 2, .... п; /= 1, 2,

. . . . т ) ;

2) определить число инверсий. Количество инверсий для данно­ го Xi определяется как число тех Zj, которые удовлетворяют усло­ вию Xi>Zj. Если для каждой пары наблюдений хг-, Zi определить функцию уц

(1, если x ^ z f ,

(7.4.41)

(О, если х х<

Г57


 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а

7.4.10

 

 

 

 

Критические значения

 

 

 

 

Критические

 

Критические значения р2,р

 

 

 

 

9

 

 

 

значения Т,

 

 

 

 

 

 

 

 

 

 

 

 

тп

 

 

 

 

 

9

 

 

 

тп

 

р

 

 

 

9

 

 

 

 

Г,0 '

0 ,0 5

 

0,10

0,90

 

0,95

0 ,9 9

 

0 ,0 5

0,01

 

0,01

0,05

0,9 5

0 ,9 9

 

 

 

 

 

0,9 5

0 ,9 9

 

 

 

 

 

 

 

 

 

 

 

 

 

11

0 , 9 3 6

0 , 9 0 7 -

0 , 8 9 0

0 , 7 4 1

0 , 7 1 5

0 , 6 6 8

2 5

0 , 7 1 1

1 ,0 6 1

50

4 , 9 2

4 , 0 1

2 , 1 3

1 , 9 5

16

0

, 9 1 4

0 , 3 8 8

0

, 8 7 3

0 , 7 4 5

0 , 7 2 4

0 , 6 8 3

3 0

0 , 6 6 1

0 , 9 8 2

100

4 , 4 0

3 , 7 7

2 , 3 5

2 , 1 8

21

• 0 , 9 0 0

0 , 8 7 7

0

, 8 6 3

0 , 7 5 0

0 , 7 3 0

0 , 6 9 5

3 5

0 , 6 2 1

0 , 9 2 1

150

4 , 1 4

3 , 6 6

2 , 4 5

2

130

2 6

0 , 8 9 0

0 , 8 6 9

0 , 8 5 7

0 , 7 5 3

0 , 7 3 6

0 , 7 0 4

4 0

0 , 5 8 7

0 , 8 6 9

2 0 0

3 ’ 98

3 , 5 7

2',51

2 , 3 7

31

0 , 8 8 3

0 , 8 6 3

0 , 8 5 1

0 , 7 5 6

0 , 7 4 0

0 , 7 1 1

4 5

0 , 5 5 8

0 , 8 2 5

2 5 0

3 , 8 7

3 , 5 1

2 ' 5 5

2 , 4 2

36

0 , 8 8 7

0 , 8 5 7

0 , 8 4 7

0 , 7 5 8

0 , 7 4 4

0 , 7 1 7

50

0 , 5 3 3

0 , 7 8 7

3 0 0

3 ’7 9

3 , 4 7

2 , 5 9

2 , 4 6

41

0 , 8 7 2

0 , 8 5 4

0 , 8 4 4

0 , 7 6 0

0 , 7 4 7

0 , 7 2 2

60

0 , 4 9 2

0 , 7 2 3

3 5 0

3 , 7 2

3 ^ 4 4

2 , 6 2

2 , 5 0

4 6

0 , 8 6 8

0 , 8 5 1

0 , 8 4 1

0 , 7 6 2

0 , 7 5 0

0 , 7 2 6

7 0

0 , 4 5 9

0 , 6 7 3

4 0 0

3 , 6 7

3 , 4 1

2 , 6 4

2 , 5 2

51

0 , 8 6 5

0 , 8 4 8

0

, 8 3 9

0 , 7 6 4

0

, 7 5 2

0 , 7 2 9

8 0

0 , 4 3 2

0 , 6 3 1

4 5 0

3',63

3 , 3 9

2 , 6 6

2

55

61

0 , 8 5 9

0 , 8 4 3

0 , 8 3 5

0 , 7 6 6

0 , 7 5 5

0 , 7 3 5

90

0 , 4 0 9

0 , 5 9 6

5 0 0

3^60

3 , 3 7

2 , 6 7

2 , 5 7

71

0 , 8 5 5

0 , 8 4 0

0

, 8 3 2

0 , 7 6 8

0

, 5 8

0 , 7 3 9

100

0 , 3 8 9

0 , 5 6 7

5 50

3 , 5 7

з;з5

2 , 6 9

2

5

81

0 , 8 5 2

0 , 8 3 8

0

, 8 3 0

0 , 7 7 0

0 , 7 6 1

0 , 7 4 3

125

0 , 3 5 0

0 , 5 0 8

60 0

3 j 5 4

3 , 3 4

2 , 7 0

2

6

91

0 , 8 4 8

0 , 8 3 5

0 , 8 2 8

0 , 7 7 1

0 , 7 6 3

0 , 7 4 6

150

0 , 3 2 1

0 , 4 6 4

6 5 0

3 ’5 2

з-,'зз

2 ,7 1

2

6

101

0 , 8 4 6

0 , 8 3 4

0

, 8 2 6

0 , 7 8 3

0 , 7 6 4

0 , 7 4 9

175

0 , 2 9 8

0 , 4 3 0

7 0 0

3 , 5 0

3 , 3 1

2 , 7 2

2

6

20 1

0 , 8 3 2

0 , 8 2 3

0 , 8 1 8

0 , 7 8 0

0 , 7 7 4

0 , 7 6 3

2 0 0

0 , 2 8 0

0 , 4 0 3

7 5 0

3 , 4 8

3 , 3 0

2 , 7 3

2 6

301

0 , 8 2 6

0 , 8 1 8

0 , 8 1 4

0 , 7 8 3

0 , 7 7 8

0 , 7 6 9

3 0 0

0 , 2 3 0

0 , 3 2 9

8 0 0

3 , 4 6

3 , 2 9

2,74:

2 , 6

401

0 , 8 2 2

0 , 7 1 5

. 0 , 8 1 2

0 , 7 8 5

0 , 7 8 1

0 , 7 7 3

4 0 0

0 , 2 0 0

0 , 2 8 5

8 5 0

3 , 4 5

3 , 2 8

2 , 7 4

2

6

501

0 , 8 2 0

0 , 8 1 4

0 , 8 1 0

0 , 7 8 6

0 , 7 8 3

0 , 7 7 6

5 0 0

■ 0 , 1 7 3

0 , 2 5 5

9 0 0

3 , 4 3

3 , 2 8

2 , 7 5

2

6

701

0 , 8 1 6

0 , 8 1 1

0 , 8 0 8

0 , 7 8 8

0 , 7 8 5

0 , 7 7 9

7 0 0

0 , 1 5 1

0 , 2 1 5

9 5 0

3 , 4 2

3 , 2 7

2 , 7 6

2

6

1001

0 , 8 1 3

0 , 8 0 9

0 , 8 0 7

0 , 7 8 9

0 , 7 8 7

0 , 7 8 2

1000

1 , 1 2 7

0 , 1 8 0

1000

3 , 4 1

3 , 2 6

2 , 7 6

2 , 6