Файл: Проверка выборочного распределения.docx

Добавлен: 30.10.2018

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

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

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

ПРИЛОЖЕНИЕ Б

Текст М-файла для проверки выборки на соответствие нормальному закону распределения по критерию Пирсона


qz=0.05; % выбрали уровень значимости

ResTabl=Tabl(1,1:6); % взяли первую строку

for k1=2:k, % берем остальные строки таблицы

if ResTabl(end,6)<5, % предыдущее npj<5 - будем суммировать

ResTabl(end,3)=Tabl(k1,3); % новая правая граница интервала

ResTabl(end,4:6)=ResTabl(end,4:6)+Tabl(k1,4:6); % суммируем

else % предыдущее npj>=5 - будем дописывать строку

ResTabl=[ResTabl;Tabl(k1,1:6)]; % дописываем строку

end

end

if ResTabl(end,6)<5, % последнее npj<5

ResTabl(end-1,3)=ResTabl(end,3); % новая правая граница

ResTabl(end-1,4:6)=ResTabl(end-1,4:6)+ResTabl(end,4:6);

ResTabl=ResTabl(1:end-1,:); % отбросили последнюю строку

end

kn=size(ResTabl,1); % число объединенных интервалов

ResTabl(:,1)=[1:kn]'; % новые номера интервалов

ResTabl(:,7)=(ResTabl(:,4)-ResTabl(:,6)).^2./ResTabl(:,6);

disp('Сгруппированная сводная таблица результатов')

fprintf(' j aj bj')

fprintf(' nj pj npj ')

fprintf(['(nj-npj)^2/npj\n'])

fprintf('%2.0f%12.5f%12.5f%6.0f%12.5f%12.5f%12.5f\n',ResTabl')

hi2=sum(ResTabl(:,7)); % сумма элементов последнего столбца

fprintf(['Статистика Пирсона chi2=%10.5f\n'],hi2)

fprintf('Задаем уровень значимости q=%5.4f\n',qz)

chi2qz=chi2inv(1-qz,k-2-1); % квантиль

fprintf(['Квантиль chi2-распределения Пирсона '...

'chi2(1-q)=%10.5f\n'],chi2qz)

if hi2<=chi2qz,

disp('Распределение подобрано верно, т.к. chi2<=chi2(1-q)')

else

disp('Распределение подобрано неверно, т.к. chi2>chi2(1-q)')

end



ПРИЛОЖЕНИЕ В

Текст М-файла для проверки выборки на соответствие нормальному закону распределения по критерию Колмогорова

pp=0.1; %уровень значимости

fprintf('---------------------------------------------------------\n')

fprintf(['Уровень значимости = %8.5f\n'],pp);

% критерий Колмогорова

[hkolm,pkolm,kskolm,cvkolm]=...

kstest(x,[x cdf('norm',x,Mx,Sx)],pp,0);

if (hkolm>0)

fprintf('Нулевая гиппотеза принимается\n')

else

fprintf('Нулевая гиппотеза отклоняется\n')

end

fprintf(['Статистика Колмогорова (лямбда) = %8.5f\n'],kskolm);% Статистика Колмогорова

fprintf(['Квантиль распределения Колмогорова = %8.5f\n'],cvkolm);% Квантиль распределения Колмогорова

fprintf(['Критический уровень значимости = %8.5f\n'],pkolm);% критический уровень значимости

figure

cdfplot(x); % эмпирическая функция распределения

xpl=linspace(xl,xr,500); % для графика F(x)

ypl=cdf('norm',xpl,Mx,Sx);

hold on % для рисования на этом же графике

plot(xpl,ypl,'r'); % дорисовали F(x)

hold off

set(get(gcf,'CurrentAxes'),...

'FontName','Times New Roman Cyr','FontSize',12)

title('\bfЭмпирическая и теоретическая функции распределения ')

xlabel('\itx') % метка оси x

ylabel('\itf\rm(\itx\rm)') % метка оси y


ПРИЛОЖЕНИЕ Г

Квантили распределения Колмогорова


p

1-p

p

1-p

p

1-p

0.99

0.44

0.50

0.83

0.15

1.14

0.90

0.57

0.40

0.89

0.10

1.22

0.80

0.64

0.30

0.97

0.05

1.36

0.70

0.71

0.25

1.02

0.02

1.52

0.60

0.77

0.20

1.07

0.01

1.63




ПРИЛОЖЕНИЕ Д

Значения функции Лапласа


х

Ф(х)

х

Ф(х)

х

Ф(х)

х

Ф(х)

0,00

0,0000

0,23

0,0910

0,46

0,1772

0,69

0,2549

0,01

0,0040

0,24

0,0948

0,47

0,1808

0,70

0,2580

0,02

0,0080

0,25

0,0987

0,48

0,1844

0,71

0,2611

0,03

0,0120

0,26

0,1026

0,49

0,1879

0,72

0,2642

0,04

0,0160

0,27

0,1064

0,50

0,1915

0,73

0,2673

0,05

0,0199

0,28

0,1103

0,51

0,1950

0,74

0,2703

0,06

0,0239

0,29

0,1141

0,52

0,1985

0,75

0,2734

0,07

0,0279

0,30

0,1179

0,53

0,2019

0,76

0,2764

0,08

0,0319

0,31

0,1217

0,54

0,2054

0,77

0,2794

0,09

0,0359

0,32

0,1255

0,55

0,2088

0,78

0,2823

0,10

0,0398

0,33

0,1293

0,56

0,2123

0,79

0,2852

0,11

0,0438

0,34

0,1331

0,57

0,2157

0,80

0,2881

0,12

0,0478

0,35

0,1368

0,58

0,2190

0,81

0,2910

0,13

0,0517

0,36

0,1406

0,59

0,2224

0,82

0,2939

0,14

0,0557

0,37

0,1443

0,60

0,2257

0,83

0,2967

0,15

0,0596

0,38

0,1480

0,61

0,2291

0,84

0,2995

0,16

0,0636

0,39

0,1517

0,62

0,2324

0,85

0,3023

0,17

0,0675

0,40

0,1554

0,63

0,2357

0,86

0,3051

0,18

0,0714

0,41

0,1591

0,64

0,2389

0,87

0,3078

0,19

0,0753

0,42

0,1628

0,65

0,2422

0,88

0,3106

0,20

0,0793

0,43

0,1664

0,66

0,2454

0,89

0,3133

0,21

0,0832

0,44

0,1700

0,67

0,2486

0,90

0,3159

0,22

0,0871

0,45

0,1736

0,68

0,2517

0,91

0,3186

0,92

0,3212

1,33

0,4082

1,74

0,4591

2,30

0,4893

0,93

0,3238

1,34

0,4099

1,75

0,4599

2,32

0,4898

0,94

0,3264

1,35

0,4115

1,76

0,4608

2,34

0,4904

0,95

0,3289

1,36

0,4131

1,77

0,4616

2,36

0,4909

0,96

0,3315

1,37

0,4147

1,78

0,4625

2,38

0,4913

0,97

0,3340

1,38

0,4162

1,79

0,4633

2,40

0,4918

0,98

0,3365

1,39

0,4177

1,80

0,4641

2,42

0,4922

0,99

0,3389

1,40

0,4192

1,81

0,4649

2,44

0,4927

1,00

0,3413

1,41

0,4207

1,82

0,4656

2,46

0,4931

1,01

0,3438

1,42

0,4222

1,83

0,4664

2,48

0,4634

1,02

0,3461

1,43

0,4236

1,84

0,4671

2,50

0,4938

1,03

0,3485

1,44

0,4251

1,85

0,4678

2,52

0,4941

1,04

0,3508

1,45

0,4265

1,86

0,4686

2,54

0,4945

1,05

0,3531

1,46

0,4279

1,87

0,4693

2,56

0,4948

1,06

0,3554

1,47

0,4292

1,88

0,4699

2,58

0,4951

1,07

0,3577

1,48

0,4306

1,89

0,4706

2,60

0,4953

1,08

0,3599

1,49

0,4319

1,90

0,4713

2,62

0,4956

1,09

0,3621

1,50

0,4332

1,91

0,4719

2,64

0,4959

1,10

0,3643

1,51

0,4345

1,92

0,4726

2,66

0,4961

1,11

0,3665

1,52

0,4357

1,93

0,4732

2,68

0,4963

1,12

0,3686

1,53

0,4370

1,94

0,4738

2,70

0,4965

1,13

0,3708

1,54

0,4382

1,95

0,4744

2,72

0,4967

х

Ф(х)

х

Ф(х)

х

Ф(х)

х

Ф(х)

1,14

0,3729

1,55

0,4394

1,96

0,4750

2,74

0,4969

1,15

0,3749

1,56

0,4406

1,97

0,4756

2,76

0,4971

1,16

0,3770

1,57

0,4418

1,98

0,4761

2,78

0,4973

1,17

0,3790

1,58

0,4429

1,99

0,4767

2,80

0,4974

1,18

0,3810

1,59

0,4441

2,00

0,4772

2,82

0,4976

1,19

0,3830

1,60

0,4452

2,02

0,4783

2,84

0,4977

1,20

0,3849

1,61

0,4463

2,04

0,4793

2,86

0,4979

1,21

0,3869

1,62

0,4474

2,06

0,4803

2,88

0,4980

1,22

0,3883

1,63

0,4484

2,08

0,4812

2,90

0,4881

1,23

0,3907

1,64

0,4495

2,10

0,4821

2,92

0,4982

1,24

0,3925

1,65

0,4505

2,12

0,4830

2,94

0,4984

1,25

0,3944

1,66

0,4515

2,14

0,4838

2,96

0,48846

1,26

0,3962

1,67

0,4525

2,16

0,4846

2,98

0,49856

1,27

0,3980

1,68

0,4535

2,18

0,4854

3,00

0,49865

1,28

0,3997

1,69

0,4545

2,20

0,4861

3,20

0,49931

1,29

0,4015

1,70

0,4554

2,22

0,4868

3,40

0,49966

1,30

0,4032

1,71

0,4564

2,24

0,4875

3,60

0,49984

1,31

0,4049

1,72

0,4573

2,26

0,4881

3,80

0,499928

1,32

0,4066

1,73

0,4582

2,28

0,4887

4,00

0,499968

 

 

 

 

 

 

5,00

0,499997




ПРИЛОЖЕНИЕ Е

Квантили распределения Пирсона


Число степеней свободы, f

Уровни значимости p


0,99


0,98


0,95


0,90


0,80


0,70


0,50


0,30

1

0,00016

0,0006

0,0039

0,016

0,064

0,148

0,455

1,07

2

0,02

0,040

0,103

0,211

0,446

0,713

1,386

2,41

3

0,115

0,185

0,352

0,584

1,005

1,424

2,366

3,66

4

0,30

0,43

0,71

1,06

1,65

2,19

3,36

4,9

5

0,55

0,75

1,14

1,61

2,34

3,00

4,35

6,1

6

0,87

1,13

1,63

2,20

3,07

3,83

5,35

7,2

7

1,24

1,56

2,17

2,83

3,82

4,67

6,35

8,4

8

1,65

2,03

2,73

3,49

4,59

5,53

7,34

9,5

9

2,09

2,53

3,32

4,17

5,38

6,39

8,34

10,7

10

2,56

3,06

3,94

4,86

6,18

7,27

9,34

11,8

11

3,1

3,6

4,6

5,6

7,0

8,1

10,3

12,9

12

3,6

4,2

5,2

6,3

7,8

9,0

11,3

14,0

13

4,1

4,8

5,9

7,0

8,6

9,9

12,3

15,1

14

4,7

5,4

6,6

7,8

9,5

10,8

13,3

16,2

15

5,2

6,0

7,3

8,5

10,3

11,7

14,3

17,3

16

5,8

6,6

8,0

9,3

11,2

12,6

15,3

18,4

17

6,4

7,3

8,7

10,1

12,0

13,5

16,3

19,5

18

7,0

7,9

9,4

10,9

12,9

14,4

17,3

20,6

19

7,6

8,6

10,1

11,7

13,7

15,4

18,3

21,7

20

8,3

9,2

10,9

12,4

14,6

16,3

19,3

22,8

21

8,9

9,9

11,6

13,2

15,4

17,2

20,3

23,9

22

9,5

10,6

12,3

14,0

16,3

18,1

21,3

24,9

23

10,2

11,3

13,1

14,8

17,2

19,0

22,3

26,0

24

10,9

12,0

13,8

15,7

18,1

19,9

23,3

27,1

25

11,5

12,7

14,6

16,5

18,9

20,9

24,3

28,2

26

12,2

13,4

15,4

17,3

19,8

21,8

25,3

29,3

27

12,9

14,1

16,2

18,1

20,7

22,7

26,3

30,3

28

13,6

14,8

16,9

18,9

21,6

23,6

27,3

31,4

29

14,3

15,6

17,7

19,8

22,4

24,6

28,3

32,5

30

15,0

16,3

18,5

20,6

23,4

25,5

29,3

33,5














Число степеней свободы, f


Уровни значимости р



0,20


0,10


0,05


0,02


0,01


0,005


0,002


0,001

1

1,64

2,7

3,8

5,4

6,6

7,9

9,5

10,8

2

3,22

4,6

6,0

7,8

9,2

10,6

12,4

13,8

3

4,64

6,3

7,8

9,8

11,3

12,8

14,8

16,3

4

6,0

7,8

9,5

11,7

13,3

14,9

16,9

18,5

5

7,3

9,2

11,1

13,4

15,1

16,3

18,9

20,5

6

8,6

10,6

12,6

15,0

16,8

18,6

20,7

22,5

7

9,8

12,0

14,1

16,6

18,5

20,3

22,6

24,3

8

11,0

13,4

15,5

18,2

20,1

21,9

24,3

26,1

9

12,2

14,7

16,9

19,7

21,7

23,6

26,1

27,9

10

13,4

16,0

18,3

21,2

23,2

25,2

27,7

29,6

11

14,6

17,3

19,7

22,6

24,7

26,8

29,4

31,3

12

15,8

18,5

21,0

24,1

26,2

28,3

31

32,9

13

17,0

19,8

22,4

25,5

27,7

29,8

32,5

34,5

14

18,2

21,1

23,7

26,9

29,1

31,3

34

36,1

15

19,3

22,3

25,0

28,3

30,6

32,8

35,5

37,7

16

20,5

23,5

26,3

29,6

32,0

34,3

37

39,2

17

21,6

24,8

27,6

31,0

33,4

35,7

38,5

40,8

18

22,8

26,0

28,9

32,3

34,8

37,2

40

42,3

19

23,9

27,2

30,1

33,7

36,2

38,6

41,5

43,8

20

25,0

28,4

31,4

35,0

37,6

40,0

43

45,3

21

26,2

29,6

32,7

36,3

38,9

41,4

44,5

46,8

22

27,3

30,8

33,9

37,7

40,3

42,8

46

48,3

23

28,4

32,0

35,2

39,0

41,6

44,2

47,5

49,7

24

29,6

33,2

36,4

40,3

43,0

45,6

48,5

51,2

25

30,7

34,4

37,7

41,6

44,3

46,9

50

52,6

26

31,8

35,6

38,9

42,9

45,6

48,3

51,5

54,1

27

32,9

36,7

40,1

44,1

47,0

49,6

53

55,5

28

34,0

37,9

41,3

45,4

48,3

51,0

54,5

56,9

29

35,1

39,1

42,6

46,7

49,6

52,3

56

58,3

30

36,3

40,3

43,8

48,0

50,9

53,7

57,5

59,7















ПРИЛОЖЕНИЕ Ж

Значение функции


u

0

1

2

3

4

5

6

7

8

9

0

0,3938

3989

3989

3988

3986

3984

3982

3980

3977

3973

0,1

3970

3965

3961

3956

3951

3945

3939

3932

3925

3918

0,2

3910

3902

3894

3885

3876

3867

3857

3847

3836

3825

0,3

3814

3802

3790

3778

3765

3752

3739

3726

3712

3697

0,4

3683

3668

3653

3637

3621

3605

3589

3572

3555

3538

0,5

3521

3503

3485

3467

3448

3429

3410

391

3372

3352

0,6

3332

3312

3292

3271

3251

3230

3209

3187

3166

3144

0,7

3123

3101

3079

3056

3034

3011

2989

2966

2943

2920

0,8

2897

2874

2850

2827

2803

2780

2756

2732

2709

265

0,9

2661

2637

3613

2589

2565

2541

2516

2492

2468

2444

1

2420

2396

2371

2347

2323

2299

2275

2251

2227

2203

1,1

2179

2155

2131

2107

2083

2059

2036

2012

1989

1965

1,2

1942

1919

1895

1872

1849

1826

1804

1781

1758

1736

1,3

1714

1691

1669

1647

1626

1604

1582

1561

1539

1518

1,4

1497

1476

1456

1435

1415

1394

1374

1354

1334

1315

1,5

1295

1276

1257

1238

1219

1200

1182

1163

1145

1127

1,6

1109

1092

1074

1057

100

1023

1006

0989

0973

0957

1,7

0940

0925

0909

0893

0878

0863

0848

0833

0818

0804

1,8

0790

0775

0761

0748

0734

0721

0707

0694

0681

0669

1,9

0656

0644

0632

0620

0608

0596

0584

0573

0562

0551

2

0540

0529

0519

0508

0498

0478

0478

0468

0459

0449

2,1

0440

0431

0422

0413

0404

0395

0387

0379

0371

0363

2,2

0355

0347

0339

0332

0325

0317

0310

0303

0297

0290

2,3

0283

0277

0270

0264

0258

0252

0246

0241

0235

0229

2,4

0224

0219

0213

0208

0203

0198

0194

0189

0184

0180

2,5

0175

0171

0167

0163

0158

0154

0151

0147

0143

0139

2,6

0136

0132

0129

0126

0122

0119

0116

0113

0110

0107

2,7

0104

0101

0099

0096

0093

0091

0088

0086

0084

0081

2,8

0079

0077

0075

0073

0071

0069

0067

0065

0063

0061

2,9

0060

0058

0056

0055

0053

0051

0050

0048

0047

0046

3

0,0044

0043

0042

0040

0039

0038

0037

0036

0035

0034