182430 (629465), страница 3
Текст из файла (страница 3)
;
По таблице критических точек критерия Стьюдента определим Ткр.= Т( q, f ) по уровню значимости q и числу степеней свободы f = N-2. Если |Тнабл|<Ткр, то гипотеза H0 – справедлива, т.е. коэффициент корреляции xy - незначим. В противном случае, нулевая гипотез H0 отвергается, т.е. случайные величины X и Y связаны линейной зависимостью (критическая область двусторонняя).
Р
ис.5. Критическая область критерия Стьюдента..
При использовании метода наименьших квадратов для вычисления коэффициента корреляции и построения уравнения регрессии предполагается, что X и Y имеют нормальное распределение.
8.3. Использование корреляционной таблицы для вычисления коэффициента корреляции
Если число экспериментов велико, то составляются корреляционные таблицы. Для этого среди результатов эксперимента выбираются xmin, xmax, ymin, ymax. Интервал [xmin, xma)] возможных значений X делим с шагом h1 на n частичных интервалов, Интервал [ymin,ymax] для Y делим с шагом h2 на m частичных интервалов. Границы интервалов по X записываются в 1-ый столбец, по Y - в 1-ую строку.
Для каждой пары (xi, yi) определяем в какую строку попало значение xi и в какой столбец yi. В клетку, расположенную на пересечении найденной строки и столбца, ставим палочку (или точку) . Операцию проводим для всех пар. Подчитываем число палочек (точек) в каждой клетке и записываем полученное число в клетку. Просуммируем числа, стоящие в 1- ой строке, получим частоту
- число пар (xi,yi), у которых первая координата попала в первый частичный интервал. Проведём суммирование по всем остальным строкам, полученные числа
заносим в последний столбец.
Таблица 7
| Y X, U | [y0, y1) y1*, v1 | [y1, y2) y2*, v2 | …… | [yj1,yj) yj*, vj C2 | …… | [ym-1, ym) ym*, vm | |
| [x0, x1) x1*, u1 | | | …… | | …… | | |
| [x1, x2) x2*, u2 | | | …… | | …… | | |
| ……… | ……… | ……… | …… | ……… | …… | ………… | ………… |
| [xi-1,xi) C1, xi*, ui | | | …… | | …… | | |
| [xn1,xn) xn*,un | | | …… | | …… | | |
| | | | …… | | …… | | N |
Просуммируем величины, которые стоят в первом столбце. Получим частоту
- число пар (xi, yi), у которых y попадает в первый интервал. Найдём суммы по всем столбцам. Полученное значение запишем в последнюю строку. Суммы полученных значений равны N:
По виду корреляционной таблице можно судить о виде корреляционной зависимости.
Вычислим середины частичных интервалов
;
i=1,…,n; j=1,…,m.
Внесем найденные значения в корреляционную таблицу. По таблице вычислим оценки математических ожиданий и дисперсий
;
;
;
;
;
;
.
Коэффициент линейной корреляции определяются по формуле:
.
Для простоты вычислений обычно используют замену переменных:
;
;
где С1 и С2 – значения xi* и yj* соответствующие максимальной частоте
. Желательно, чтобы клетка с данной частотой находилась в середине таблицы. Точку (С1,С2) называют ложным нулем. Переменные U и V – принимают значения: 0; 1; 2,…
,
,
,
;
.
При вычислениях используем, что
;
.
Коэффициент корреляции вычисляется по формуле:
.
Вернемся к исходным переменным:
;
;
;
.
Уравнения регрессии:
;
.
Графики функций пересекаются в точке
.
Пример:
Даны результаты 78 экспериментов:
| X | Y | X | Y | X | Y | X | Y |
| 73 | -291 | 57 | -219 | 61 | -241 | 68 | -264 |
| 69 | -270 | 71 | -281 | 62 | -243 | 62 | -240 |
| 72 | -279 | 66 | -262 | 63 | -245 | 70 | -277 |
| 72 | -282 | 76 | -302 | 71 | -282 | 70 | -279 |
| 65 | -254 | 70 | -275 | 65 | -252 | 65 | -253 |
| 67 | -264 | 68 | -267 | 70 | -276 | 70 | -275 |
| 56 | -216 | 74 | -290 | 70 | -276 | 63 | -248 |
| 70 | -276 | 68 | -266 | 63 | -246 | 63 | -243 |
| 63 | -248 | 71 | -283 | 73 | -284 | 67 | -264 |
| 64 | -253 | 60 | -237 | 68 | -271 | 68 | -267 |
| 70 | -276 | 56 | -222 | 59 | -227 | 55 | -213 |
| 67 | -262 | 71 | -281 | 64 | -256 | 56 | -218 |
| 60 | -234 | 68 | -269 | 79 | -309 | 58 | -223 |
| 80 | -313 | 66 | -257 | 77 | -300 | 70 | -278 |
| 71 | -278 | 60 | -235 | 78 | -310 | 59 | -236 |
| 74 | -292 | 70 | -275 | 66 | -255 | 68 | -263 |
| 68 | -271 | 69 | -276 | 63 | -252 | 69 | -268 |
| 65 | -256 | 72 | -282 | 69 | -274 | 63 | -243 |
| 73 | -291 | 70 | -277 | 74 | -291 | 70 | -271 |
| 63 | -243 | 69 | -270 |
Начало первого интервала x0 = 53, y0 = –321;
Длина интервала h1 = 5, h2 = 17.
-
Построить корреляционное поле для 4-ых столбцов X и Y и методом “натянутой нити” найти линейные функции регрессии.
-
Составить корреляционную таблицу. Вычислить коэффициент линейной корреляции, найти уравнения регрессий и построить их графики.
-
Проверить гипотезу о незначимости коэффициента корреляции.
Решение.
1. По последним столбцам X и Y находим:
xmin=55; ymin=-279;
xmax=70; ymax=-213;
На осях отображаем тот промежуток, где находятся значения X и Y. Представляя в виде точек пары чисел (x1; yj) строим корреляционное поле:
Используя метод “натянутой нити”, проведём прямую. На прямой выберем две точки (57, -220) и (69, -270), расположенные достаточно далеко друг от друга.. Подставляя значения в функцию y=ax+b, получим систему уравнений относительно a и b.
,
Получим решение a = - 4,17; b = 17,69. Уравнение линейной регрессии имеет вид: y = - 4,17 x + 17,69.
-
Найдём минимальные и максимальные значения X и Y среди результатов эксперимента:
xmin=55; ymin=-313; xmax=80; ymax=-213;
Составим корреляционную таблицу с шагом h1=5 по X и h2=17 по Y. Учитываем, что левая граница входит в интервал, а правая нет.
Клетка в шапке сверху содержит границы интервала по Y [yj, yj+1], значение середины интервала yj* и значение середины интервала для условной переменной V. Клетка в шапке слева содержит границы интервала по X [xi, xi+1], значение середины интервала xi* и значение середины интервала для условной переменной U.
Произвольная клетка таблицы содержит число результатов
, попавших в соответствующие интервалы. В нижней строке записываются суммы чисел в столбцах. В крайнем левом столбце – суммы чисел в строках.
Таблица 8.
| Y,V X,U | [321,-304) -312,5; -2 | [304,-287) -295,5; -1 | [287,-270) -278,5; 0 | [270,-253) -261,5; 1 | [253,-236) -244,5; 2 | [236,-219) -227, 5; 3 | [219,-202) -210,5 ; 4 | nx nu |
| [53,58) 55,5;-3 | . 1 | . . . . 4 | 5 | |||||
| [58,63) 60,5;-2 | . . . . 4 | . . . . 5 | 9 | |||||
| [63,68) 65,5;-1 | . . . . 9 | . . . . . 11 | 20 | |||||
| [68,73) 70,5;0 | .. 24 | . . . . 9 | 33 | |||||
| [73,78) 75,5;1 | . . . . 7 | . 1 | 8 | |||||
| [78,83) 80,5;2 | . . . 3 | 3 | ||||||
| ny, nv | 3 | 7 | 25 | 18 | 15 | 6 | 4 | 78 |
Переход к условным вариантам.
;
;
C1=70,5; С2=-278,5 – координаты клетки с максимальным числом результатов экспериментов.
,
,…,
;
,
,…,
;
Вычисляем средние:
.
Вычислим среднее квадратов:
.
Вычислим среднее квадратическое отклонение:
;
;
;
Коэффициент корреляции:
;
Находим статистические характеристики X, Y:
;
;
;
;
Уравнение регрессии:
;
;
;
(I);
;
(II)
Определим координаты двух точек для каждого графика:
| X | 60 | 75 | ||
| Y | -231,4 | -291 | ||
| Y | -300 | -220 | ||
| X | 76,46 | 58,06 | ||
Графики пересеклись в точке M(68; -263,2)
-
Проверим гипотезу о незначимости коэффициента корреляции. Наблюдаемое значение критерия:
; N=max{n, m};
n, m – число частичных интервалов по X и Y.
n = 6; m = 7; N = 7.
;
Tкр=T(0,05; N-2)=T(0,05; 5)=2,57 – по таблице распределения Стьюдента.
Так как |Tнабл|=6,65>2,57, то гипотеза отвергается, следовательно r xy значим.
Вывод
В курсовую работу вошли задачи, решаемые на стадии предварительного эксперимента. При решении этих задач использованы идеи и методы математической статистики, в частности ее разделы - оценивание параметров и проверка статистических гипотез. Используя эти методы, проверяются следующие гипотезы: о воспроизводимости результатов эксперимента, о виде распределения результатов эксперимента, о наличии корреляционных связей между факторами и переменной состояния и др.
Список литературы
1.Егоров А.Е., Азаров Г.Н., Коваль А.В. Исследование устройств и систем автоматики методом планирования эксперимента. – К.: Вища школа, 1986.
2.Бондарь А.Г., Статюха Г.А. Планирование эксперимента в химической технологии. – К.: Вища школа, 1978.
3.Кафаров В.В. Методы кибернетики в химии и химической технологии. – М.: Химия, 1971.
4.Колде Я.К. Практикум по теории вероятностей и математической статистике. – М.: Высшая школа, 1991.
5.Твердохлебов Г.Н., Бродский А.Л., Старобина Е.К., Кутакова Д.А. Методические указания по математическим методам анализа и планирования эксперимента для студентов всех химических специальностей. -Ворошиловград, 1985.
Приложение 1
(таблица значений функции Лапласа Ф(х))
(Таблица значений функции
| x | Ф(x) | x | Ф(x) | x | Ф(x) | x | Ф(x) | |||||||
| 0.00 | 0.0000 | 0.22 | 0.0871 | 0.44 | 0.1700 | 0.66 | 0.2454 | |||||||
| 0.01 | 0.0040 | 0.23 | 0.0910 | 0.45 | 0.1736 | 0.67 | 0.2486 | |||||||
| 0.02 | 0.0080 | 0.24 | 0.0948 | 0.46 | 0.1772 | 0.68 | 0.2517 | |||||||
| 0.03 | 0.0120 | 0.25 | 0.0987 | 0.47 | 0.1808 | 0.69 | 0.2549 | |||||||
| 0.04 | 0.0160 | 0.26 | 0.1026 | 0.48 | 0.1844 | 0.70 | 0.2580 | |||||||
| 0.05 | 0.0199 | 0.27 | 0.1064 | 0.49 | 0.1879 | 0.71 | 0.2611 | |||||||
| 0.06 | 0.0239 | 0.28 | 0.1103 | 0.50 | 0.1915 | 0.72 | 0.2642 | |||||||
| 0.07 | 0.0279 | 0.29 | 0.1141 | 0.51 | 0.1950 | 0.73 | 0.2673 | |||||||
| 0.08 | 0.0319 | 0.30 | 0.1179 | 0.52 | 0.1985 | 0.74 | 0.2703 | |||||||
| 0.09 | 0.0359 | 0.31 | 0.1217 | 0.53 | 0.2019 | 0.75 | 0.2734 | |||||||
| 0.10 | 0.0398 | 0.32 | 0.1255 | 0.54 | 0.2054 | 0.76 | 0.2764 | |||||||
| 0.11 | 0.0438 | 0.33 | 0.1293 | 0.55 | 0.2088 | 0.77 | 0.2794 | |||||||
| 0.12 | 0.0478 | 0.34 | 0.1331 | 0.56 | 0.2123 | 0.78 | 0.2823 | |||||||
| 0.13 | 0.0517 | 0.35 | 0.1368 | 0.57 | 0.2157 | 0.79 | 0.2852 | |||||||
| 0.14 | 0.0557 | 0.36 | 0.1406 | 0.58 | 0.2190 | 0.80 | 0.2881 | |||||||
| 0.15 | 0.0596 | 0.37 | 0.1443 | 0.59 | 0.2224 | 0.81 | 0.2910 | |||||||
| 0.16 | 0.0636 | 0.38 | 0.1480 | 0.60 | 0.2257 | 0.82 | 0.2939 | |||||||
| 0.17 | 0.0675 | 0.39 | 0.1517 | 0.61 | 0.2291 | 0.83 | 0.2967 | |||||||
| 0.18 | 0.0714 | 0.40 | 0.1554 | 0.62 | 0.2324 | 0.84 | 0.2995 | |||||||
| 0.19 | 0.0753 | 0.41 | 0.1591 | 0.63 | 0.2357 | 0.85 | 0.3023 | |||||||
| 0.20 | 0.0793 | 0.42 | 0.1628 | 0.64 | 0.2389 | 0.86 | 0.3051 | |||||||
| 0.88 | 0.3106 | 1.14 | 0.3729 | 1.40 | 0.4192 | 1.66 | 0.4515 | |||||||
| 0.89 | 0.3133 | 1.15 | 0.3749 | 1.41 | 0.4207 | 1.67 | 0.4525 | |||||||
| 0.90 | 0.3159 | 1.16 | 0.3770 | 1.42 | 0.4222 | 1.68 | 0.4535 | |||||||
| 0.91 | 0.3186 | 1.17 | 0.3790 | 1.43 | 0.4236 | 1.69 | 0.4545 | |||||||
| 0.92 | 0.3212 | 1.18 | 0.3810 | 1.44 | 0.4251 | 1.70 | 0.4554 | |||||||
| 0.93 | 0.3238 | 1.19 | 0.3830 | 1.45 | 0.4265 | 1.71 | 0.4564 | |||||||
| 0.94 | 0.3264 | 1.20 | 0.3849 | 1.46 | 0.4279 | 1.72 | 0.4573 | |||||||
| 0.95 | 0.3289 | 1.21 | 0.3869 | 1.47 | 0.4292 | 1.73 | 0.4582 | |||||||
| 0.96 | 0.3315 | 1.22 | 0.3883 | 1.48 | 0.4306 | 1.74 | 0.4591 | |||||||
| 0.97 | 0.3340 | 1.23 | 0.3907 | 1.49 | 0.4319 | 1.75 | 0.4599 | |||||||
| 0.98 | 0.3365 | 1.24 | 0.3925 | 1.50 | 0.4332 | 1.76 | 0.4608 | |||||||
| 0.99 | 0.3389 | 1.25 | 0.3944 | 1.51 | 0.4345 | 1.77 | 0.4616 | |||||||
| 1.00 | 0.3413 | 1.26 | 0.3962 | 1.52 | 0.4357 | 1.78 | 0.4625 | |||||||
| 1.01 | 0.3438 | 1.27 | 0.3980 | 1.53 | 0.4370 | 1.79 | 0.4633 | |||||||
| 1.02 | 0.3461 | 1.28 | 0.3997 | 1.54 | 0.4382 | 1.80 | 0.4641 | |||||||
| 1.03 | 0.3485 | 1.29 | 0.4015 | 1.55 | 0.4394 | 1.81 | 0.4649 | |||||||
| 1.04 | 0.3508 | 1.30 | 0.4032 | 1.56 | 0.4406 | 1.82 | 0.4656 | |||||||
| 1.05 | 0.3531 | 1.31 | 0.4049 | 1.57 | 0.4418 | 1.83 | 0.4664 | |||||||
| 1.06 | 0.3554 | 1.32 | 0.4066 | 1.58 | 0.4429 | 1.84 | 0.4671 | |||||||
| 1.07 | 0.3577 | 1.33 | 0.4082 | 1.59 | 0.4441 | 1.85 | 0.4678 | |||||||
| 1.08 | 0.3599 | 1.34 | 0.4099 | 1.60 | 0.4452 | 1.86 | 0.4686 | |||||||
| 1.09 | 0.3621 | 1.35 | 0.4115 | 1.61 | 0.4463 | 1.87 | 0.4693 | |||||||
| 1.10 | 0.3643 | 1.36 | 0.4131 | 1.62 | 0.4474 | 1.88 | 0.4699 | |||||||
| 1.11 | 0.3665 | 1.37 | 0.4147 | 1.63 | 0.4484 | 1.89 | 0.4706 | |||||||
| 1.12 | 0.3686 | 1.38 | 0.4162 | 1.64 | 0.4495 | 1.90 | 0.4713 | |||||||
| 1.13 | 0.3708 | 1.39 | 0.4177 | 1.65 | 0.4505 | 1.91 | 0.4719 | |||||||
,V 














