Корреляционный анализ в системах управления. Учебное пособие
.pdf
X = 1 ∑m n x .
ni=1 i i
3.В том случае, когда статистические данные представлены
ввиде интервального вариационного ряда, при вычислении выборочного среднего в качестве значений вариант принимают середины интервалов.
Часто при решении практических задач необходимо определить характеристики, которыми являются мода и медиана.
В математической статистике мода Мо определяется по выборке, как варианта с наибольшей частотой.
Медианой называется варианта, расположенная в центре ранжированного ряда. Если ряд состоит из четного числа вариант, то медианой считают среднее арифметическое двух вариант, расположенных в центре ранжированного ряда.
Пример 1.4
Найти моду и медиану выборочной совокупности по данным, приведенным в примере 1.1.
Ответ: Мо = 24; Ме = 24.
Перечень основных числовых характеристик выборочной совокупности, используемых при экспериментальном исследовании данных
1.Размах вариационного ряда R = Xmax −Xmin . Этот показатель является наиболее простой характеристикой рассеяния и показывает диапазон варьирования величины. Этой характеристикой пользуются при работе с малыми выборками.
2.Выборочное среднее находится как взвешенное среднее арифметическое
|
|
|
|
m |
|
|
|
|
|
|
∑xini |
|
|
|
|
|
= |
i=1 |
, |
|
x |
||||||
|
|
|
||||
|
B |
|
n |
|
||
|
|
|
|
|
||
которое характеризует среднее значение признака X в пределах рассматриваемой выборки.
3. Выборочная дисперсия определяется по формуле
11
|
|
m |
|
|
|
|
|
∑(xi − |
|
B )2 ni |
|
|
x |
||||
D |
= |
i=1 |
|
|
|
|
|
|
|
||
B |
n |
||||
|
|
|
|||
и является мерой рассеяния возможных значений показателя X вокруг своего среднего значения; ее размерность совпадает с квадратом размерности варианты.
4. Выборочное среднее квадратическое отклонение
sB (X) =
D B описывает абсолютный разброс значений показателя X. Его размерность совпадает с размерностью варианты.
5. «Исправленная» дисперсия (вычисляемая при малых N,
|
|
|
|
m |
|
|
|
|
|
|
|
|
|
|
|
|
|
∑(xi − |
|
B )2 ni |
|
|
|
|
|||
|
n |
|
x |
|
|
|
|
||||||
n < 30) S2 (X) = |
D = |
i=1 |
|
|
|
|
|
|
и «исправленное» стан- |
||||
n −1 |
|
n −1 |
|||||||||||
|
B |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
||||||
дартное отклонение S(X) = |
|
S2 (X) . |
|
|
|
|
|||||||
6. Коэффициент вариации V = sB (X) |
100 %характеризует |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
B |
|
|
относительную изменчивость показателя X, т.е. относительный разброс вокруг его среднего значения xB . Коэффициент вариации является безразмерной величиной, поэтому он пригоден для сравнения рассеяния вариационных рядов, варианты которых имеют различную размерность.
Пример 1.5
Измерены два признака для одного элемента (например, длина и масса исследуемого образца). Получены следующие показатели: Х = 130 ед. дл., SХ = 5 ед, дл.; Y = 32 ед.. масс.; SY = 4 ед. масс. Какой показатель имеет большую вариативность?
Решение
Так как оба признака измеряются в разных единицах, то вариативность нельзя сравнить с помощью среднеквадратического отклонения (СКО). Необходимо вычислить относительный показатель вариации для каждого признака:
V |
= |
|
5 |
100 %=3,8 %; |
V |
|
|
= |
4 |
100 %=12,5 %. |
X |
|
130 |
|
Y |
|
32 |
|
|||
Таким образом, признак Y имеет большую вариативность, чем признак X.
12
Оценка параметров с помощью интервалов Оценка параметров с помощью интервалов заключает-
ся в определении интервалов, называемых доверительными, между границами которых с определенными доверительными вероятностями находятся истинные значения оцениваемых параметров. Интервальная оценка определяется двумя числами – концами интервала.
Доверительным интервалом называется интервал, который покрывает неизвестный параметр с заданной надежностью.
Пример 1.6
Путем измерений получены несколько блоков данных (N = 80), приведенных ниже:
2 4 2 4 1 1 1 2 0 6 |
1 2 1 2 2 4 1 1 5 1 |
0 2 4 1 2 2 1 1 1 1 |
|||
1 1 1 1 2 1 1 4 1 1 |
7 4 1 4 2 1 2 1 1 1 |
4 1 1 4 5 1 4 2 4 5 |
|||
1 6 4 1 1 2 4 1 1 1 |
0 0 4 6 4 7 4 1 1 5 |
|
|
||
|
|
|
|
||
2 4 2 4 1 1 1 2 0 6 |
1 2 1 2 2 4 1 1 5 1 |
0 2 4 1 2 2 1 1 1 1 |
1 1 1 1 2 1 1 4 1 1 |
||
7 4 1 4 2 1 2 1 1 1 |
4 1 1 4 5 1 4 2 4 5 |
1 6 4 1 1 2 4 1 1 1 |
0 0 4 6 4 7 4 1 1 5 |
||
Требуется выполнить следующие задания.
1.Получить дискретный вариационный ряд и статистическое распределение выборки.
2.Построить полигон частот.
3.Составить ряд распределения относительных частот.
4.Составить эмпирическую функцию распределения.
5.Построить график эмпирической функции распределения.
6.Найти основные числовые характеристики вариационного ряда (по возможности использовать упрощающие формулы):
а) выборочное среднее xB ;
б) выборочную дисперсию D(X);
в) выборочное среднее квадратическое отклонение s(X); г) коэффициент вариации V.
Решение
1.Для составления дискретного вариационного ряда отсортируем данные измерений по величине и расположим их в порядке возрастания:
0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 5 5 5 5 6 6 6 7 7.
13
Статистическое распределение выборки представлено в табл. 1.5, в которой первая строка – варианты (наблюдаемые значение), вторая строка – частоты появления этих вариант.
|
|
|
|
|
|
|
|
|
Таблица 1.5 |
|
|
|
|
Варианты и их частоты |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
Xi |
0 |
1 |
|
2 |
1 |
|
4 |
5 |
|
6 |
|
|
|
|
|
|
|
|
|
|
|
N i |
4 |
11 |
|
14 |
24 |
|
16 |
4 |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
2. Для построения полигона частот найдем относительные ча-
n
стоты wi = ni
ка X ( m ≤n ), и в данном примере M = 8), которые будем вычислять с одинаковой точностью.
Часто на практике для приближения теоретической функции распределения используют график накопленных частот и эмпирическую функцию распределения.
Графиком накопленных частот называется фигура, строящаяся аналогично гистограмме с той разницей, что для расчета высот прямоугольников берутся не простые, а накопленные от-
nн
носительные частоты, т.е. ni =wiн. Эти величины не убывают, и
таким образом график накопленных частот имеет вид ступенчатой «лестницы» (от 0 до 1).
Полигон частот – ломаная линия, соединяющая точки с координатами (xi,wi ) , представленная на рис. 1.3. Расчеты запишем в табл. 1.6.
|
|
|
Таблица 1.6 |
|
Относительные частоты и накопленные частоты |
||
|
|
|
|
Xi |
Ni |
Относительная частота w |
Накопленная частота wн |
|
|
i |
i |
0 |
4 |
0,050 |
0,050 |
|
|
|
|
1 |
11 |
0,161 |
0,211 |
|
|
|
|
2 |
14 |
0,175 |
0,188 |
|
|
|
|
1 |
24 |
0,100 |
0,688 |
|
|
|
|
4 |
16 |
0,200 |
0,888 |
|
|
|
|
14
|
|
|
Окончание табл. 1.6 |
|
|
|
|
Xi |
Ni |
Относительная частота w |
Накопленная частота wн |
|
|
i |
i |
5 |
4 |
0,050 |
0,918 |
|
|
|
|
6 |
1 |
0,018 |
0,975 |
|
|
|
|
7 |
2 |
0,025 |
1,000 |
|
|
|
|
Сумма |
80 |
1 |
– |
|
|
|
|
Рис. 1.3. Полигон частот вариационного ряда
3. Запишем ряд распределения (табл. 1.6) относительных частот в виде табл. 1.7, в которой первая строка – варианты (изучаемый признак), вторая строка – относительные частоты (частости).
Таблица 1.7
Распределение относительных частот появления признака
Xi |
0 |
1 |
2 |
1 |
4 |
5 |
6 |
7 |
|
|
|
|
|
|
|
|
|
Ni |
0,05 |
0,161 |
0,175 |
0,1 |
0,2 |
0,05 |
0,018 |
0,025 |
|
|
|
|
|
|
|
|
|
15
4. Эмпирическую функцию распределения найдем, используя накопленные частоты (табл. 1.7):
|
|
|
0, |
x ≤0; |
|
|
|
0,05, |
0 <x ≤1; |
|
|
|
||
|
|
|
0,213, |
1<x ≤2; |
|
|
|
0,388, |
2 <x ≤3; |
|
|
|
||
F |
* |
|
0,688, |
3 <x ≤4; |
|
(x) = |
|||
|
|
|
0,888, |
4 <x ≤5; |
|
|
|
|
5 <x ≤6; |
|
|
|
0,938, |
|
|
|
|
0,975, |
6 <x ≤7; |
|
|
|
||
|
|
|
1, |
x >7. |
|
|
|
5. Построим график эмпирической функции распределения (рис. 1.4), используя значения, полученные в п. 4).
Рис. 1.4. График эмпирической функции распределения
6. Для вычисления выборочного среднего xB и выборочной дисперсии DB с использованием приведенных выше формул целесообразно составить табл. 1.8.
Используя суммы, представленные в табл. 1.8, определяем искомые величины:
а) выборочную среднюю
|
|
|
|
m |
|
|
|
|
|
|
|
|
∑xini |
|
229 |
|
|
|
|
|
= |
i=1 |
= |
≈2,86; |
||
x |
||||||||
|
|
|
|
|||||
|
B |
|
n |
80 |
|
|||
|
|
|
|
|
||||
16
б) выборочную дисперсию
|
|
|
|
|
m |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∑(xi − |
|
B )2 ni |
|
|
|
|
|
|
|
|
|
|
|
|||||
|
x |
|
191,488 |
|
|
|
|
|
|||||||||||||||
|
D = |
i=1 |
|
|
|
|
|
|
= |
|
≈2,39; |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
B |
n |
|
|
|
80 |
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
в) выборочное среднее квадратическое отклонение |
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
sB (X) = |
D B |
= |
2,39 |
≈1,55; |
|
|
|
||||||||||||||
г) коэффициент вариации |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
V = sB (X) |
100 %= |
1,55 |
100 %≈54,2 %. |
|||||||||||||||||||
|
2,86 |
||||||||||||||||||||||
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
B |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 1.8 |
||
Расчетная таблица для вычисления выборочных величин |
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
Xi |
Ni |
|
Xi < N i |
|
|
|
|
|
(xi − |
|
B )2 |
|
(xi − |
|
B )2 Ni |
||||||||
|
|
|
|
|
|
|
|
x |
|
x |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
0 |
4 |
|
|
|
|
0 |
|
|
|
|
|
|
8,1796 |
|
|
12,7184 |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
1 |
11 |
|
|
|
|
11 |
|
|
|
|
|
|
1,4596 |
|
|
44,9748 |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
2 |
14 |
|
|
|
|
28 |
|
|
|
|
|
|
0,7196 |
|
|
10,1544 |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
1 |
24 |
|
|
|
|
72 |
|
|
|
|
|
|
0,0196 |
|
|
0,4704 |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
4 |
16 |
|
|
|
|
64 |
|
|
|
|
|
|
1,2996 |
|
|
20,7916 |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
5 |
4 |
|
|
|
|
20 |
|
|
|
|
|
|
4,5796 |
|
|
18,1184 |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
6 |
1 |
|
|
|
|
18 |
|
|
|
|
|
|
9,8596 |
|
|
29,5788 |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
7 |
2 |
|
|
|
|
14 |
|
|
|
|
|
|
17,1196 |
|
|
14,2792 |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Сумма |
80 |
|
|
|
|
229 |
|
|
|
|
|
|
|
– |
|
191,488 |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таким образом, можно сделать следующие выводы:
––величина xВ ≈2,86 характеризует среднее значение признака X;
––среднее квадратическое отклонение sB (X) описывает абсолютный разброс значений показателя X относительно среднего значения и в данном случае составляет sB (X) ≈1,55 ;
––коэффициент вариации V характеризует относительную изменчивость показателя X, т.е. относительный разброс вокруг его среднего значения x B , и в данном случае составляет V ≈54,2 %.
Ответ: x B ≈2,86 ; DB ≈2,39 ; sB (X) ≈1,55 ; V ≈54,2 %.
17
Пример 1.7
Построить вариационный ряд и график эмпирической функции распределения по результатам наблюдений (табл. 1.9).
|
|
|
|
|
|
|
Таблица 1.9 |
|
|
Результаты наблюдений для примера 1.7 |
|
|
|
||||
|
|
|
|
|
|
|
|
|
N i |
1 |
2 |
3 |
4 |
|
5 |
|
6 |
|
|
|
|
|
|
|
|
|
X i |
51 |
43 |
56 |
60 |
|
64 |
|
56 |
|
|
|
|
|
|
|
|
|
Решение
Построим вариационный ряд, упорядочив по возрастанию значения варианты, приведенные в табл. 1.10.
|
|
|
|
|
|
|
|
Таблица 1.10 |
|
|
|
Вариационный ряд |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
N i |
1 |
2 |
|
3 |
|
4 |
5 |
|
6 |
|
|
|
|
|
|
|
|
|
|
X i |
43 |
51 |
|
56 |
|
56 |
60 |
|
64 |
|
|
|
|
|
|
|
|
||
Искомая эмпирическая функция: |
|
|
|
||||||
|
|
|
0, |
x ≤43; |
|
|
|
||
|
|
|
0,16, |
43 <x ≤51; |
|
|
|
||
|
|
|
|
51<x ≤56; |
|
|
|
||
|
|
|
0,33, |
|
|
|
|||
|
|
Fб (x) = |
|
|
|
|
|
||
|
|
|
0,67, 56 <x ≤60; |
|
|
|
|||
|
|
|
0,84, 60 <x ≤64 |
|
|
|
|||
|
|
|
|
x >64. |
|
|
|
||
|
|
|
1, |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
Строим график эмпирической функции (рис. 1.5).
Рис. 1.5. График эмпирической функции распределения
18
Пример 1.8
Имеются результаты регистрации значений затухания сигнала xi на частоте 1000 Гц коммутируемого канала телефонной сети. Эти значения, измеренные в децибелах (дБ), в виде вариационного ряда представлены в табл. 1.11. Необходимо построить статистический ряд.
Таблица 1.11
Результаты наблюдений для примера 1.8
Ni |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
|
|
|
|
|
|
|
|
|
|
|
|
Xi |
25,79 |
25,98 |
25,98 |
26,12 |
26,13 |
26,49 |
26,52 |
26,60 |
26,66 |
26,69 |
26,74 |
|
|
|
|
|
|
|
|
|
|
|
|
Ni |
12 |
13 |
14 |
15 |
16 |
17 |
18 |
19 |
20 |
21 |
22 |
|
|
|
|
|
|
|
|
|
|
|
|
Xi |
26,85 |
26,90 |
26,91 |
26,96 |
27,02 |
27,11 |
27,19 |
27,21 |
27,28 |
27,30 |
27,38 |
|
|
|
|
|
|
|
|
|
|
|
|
Ni |
23 |
24 |
25 |
26 |
27 |
28 |
29 |
30 |
31 |
32 |
33 |
|
|
|
|
|
|
|
|
|
|
|
|
Xi |
27,40 |
27,49 |
27,64 |
27,66 |
27,71 |
27,78 |
27,89 |
27,89 |
28,01 |
28,10 |
28,11 |
|
|
|
|
|
|
|
|
|
|
|
|
Ni |
34 |
35 |
36 |
37 |
38 |
39 |
40 |
41 |
42 |
43 |
44 |
|
|
|
|
|
|
|
|
|
|
|
|
Xi |
28,37 |
28,38 |
28,50 |
28,63 |
28,67 |
28,90 |
28,99 |
28,99 |
29,03 |
29,12 |
29,28 |
|
|
|
|
|
|
|
|
|
|
|
|
Решение
Количество разрядов статистического ряда q следует выбрать минимальным, чтобы обеспечить достаточное количество попаданий в каждый из них; возьмем q = 6. Определим размер разряда:
|
h = |
(xmax −xmin ) |
= |
29,28 −25,79 |
=0,58 . |
|
||||||
|
|
q |
|
|
|
|||||||
|
|
|
|
|
|
6 |
|
|
|
|||
Сгруппируем наблюдения по разрядам (табл. 1.12). |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
Таблица 1.12 |
|
|
|
|
Сгруппированные наблюдения |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Ni |
1 |
|
|
2 |
|
3 |
|
4 |
|
5 |
6 |
|
|
|
|
|
|
|
|
|
|
|
|||
Xi |
25,79 |
|
26,37 |
|
26,95 |
|
27,53 |
28,12 |
28,70 |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
ni |
5 |
|
|
9 |
|
10 |
|
9 |
|
5 |
6 |
|
|
|
|
|
|
|
|
|
|
|
|||
s i =ni /n |
0,114 |
|
0,205 |
|
0,227 |
|
0,205 |
0,114 |
0,136 |
|||
n i=s i /h |
0,196 |
|
0,353 |
|
0,392 |
|
0,353 |
0,196 |
0,235 |
|||
19
На основе статистического ряда построим гистограмму (рис. 1.6), и график эмпирической функции распределения (рис. 1.7).
Рис. 1.6. Гистограмма, полученная по экспериментальным данным
Рис. 1.7. График эмпирической функции распределения полученный по экспериментальным данным
Рассмотренные примеры представления экспериментальных данных являются исходными для последующей обработки и вычисления различных параметров.
20
