Файл: Мастеров, В. А. Практика статистического планирования эксперимента в технологии биметаллов.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 21.10.2024
Просмотров: 50
Скачиваний: 0
Независимые оценки bi возможны только в том слу чае, когда система (6) распадается на (я+1) уравнений, каждое из которых содержит только один неизвестный параметр bi. Нетрудно убедиться, что это условие реали зуется, если в каждом уравнении суммы попарных произ ведений XiuXjи в левой части будут равны нулю, т. е.
£ X iu X iu = o;_i=hj. |
(8) |
U=l1
В этом случае матрицу планирования называют ортогональной. Соответствующая ей информационная матрица становится диагональной:
£ x l |
|
0 |
|
|
|
£ * ? „ |
|
|
|
|
|
и |
’ . |
|
|
|
|
|
О |
а « X W |
|
|
|
|
|
а |
|
|
|
|
|
|
' |
£ * |
L |
|
|
|
|
и |
|
а ее определитель равен: |
|
|
|
||
N |
N |
N |
|
N |
|
а = £ * 0V £ ; & • . . . ■ £ x L - . . . - £ * L ; |
|||||
0 = 1 |
0 = 1 |
0 = 1 |
|
0 = 1 |
|
соответственно определитель |
|
|
|
||
N |
N |
N |
|
N |
|
Дг = £ хЪи. £ |
Х\а-..,. £ |
х !иУи....' £ |
х пи2 , |
||
0=1 |
0=1 |
0=1 |
|
0=1 |
|
а неизвестные параметры |
|
|
|
||
|
|
N |
|
|
|
|
|
2 Xiuyu |
|
|
|
|
|
0 = 1 |
|
|
(9) |
|
|
N |
|
|
|
|
|
|
|
|
|
|
|
S |
* ? . |
|
|
|
|
о= 1 |
|
|
|
Из допущения нормальности |
распределения ошибок |
||||
измерения у и условия ортогональности |
(8) |
следует, что |
15
и оценки bi распределены нормально с центром распре деления (3;, а дисперсия оценки bi равна:
2 4,
и= 1
Если переменные кодированы, как указано выше (minХ{ = — 1, max^i== + l) и количество уровней выше и ниже основного уровня Х г= 0 одинаковы, то оценка b t еще больше упрощается:
N
( 10)
и=1
а дисперсия s2{6,} не зависит от номера i
= |
( 11) |
л
По закону накопления ошибок дисперсия s2{y} рас четных значений функции отклика (4)
з2 [у} = s* [Ь(} ■(1 + X 2 + Х\ + ■■•+ X I) = s- ^ i (1 + г \
где г2=Х\-\- Х\ +••• + Х2п— радиус гиперсферы, |
центр |
которой совпадает с центром эксперимента X t= 0 . |
Сле- |
л |
|
довательно ошибка прогноза у зависит только от длины радиуса г, а не от его направления в факторном прост ранстве. Это свойство матрицы планирования называют
ротатабельностыо.
Повторим необходимые |
и достаточные |
условия для |
применения расчетов модели с помощью |
формул (4), |
|
(10), (11): |
|
|
выполняются основные |
допущения регрессионного |
анализа, указанные на стр. 12; симметричность плана относительно центра экспери
мента |
|
£ Х 1и= 0 , |
(12) |
0=1 |
|
т. е. сумма элементов Хы любого столбца матрицы пла нирования равна нулю;
16
нормировка
I > X l = N, |
(13) |
и= 1
т. е. сумма квадратов элементов любого столбца равна числу опытов; ортогональность
% х1их 1и= о, i=hj,
и= 1
т. е сумма построчных парных произведений элементов Xiu, Xju любых двух столбцов матрицы планирования равна нулю.
СТАТИСТИЧЕСКИЙ АНАЛИЗ УРАВНЕНИЯ РЕГРЕССИИ
Анализ проводится после |
вычисления коэффициентов |
|
регрессии |
Возможны два |
варианта анализа: |
1)нет оценки для дисперсии воспроизводимости s2{y}.
2)есть оценка s2{y}.
Впервом варианте сначала вычисляется остаточная дис
персия
|
N |
|
|
|
2 |
\уи — уи? |
(14) |
|
2 _ u=1 |
||
|
*0 — |
N — 1 |
|
|
|
|
|
характеризующая разброс измеренных |
уи относительно |
||
средней арифметической у: |
|
|
|
У = |
(#i + |
У2 Н------- |
|
Затем подсчитывается дисперсия1
" |
л |
|
2 |
Iуи — уи\* |
|
м=1________ |
(15) |
|
si = |
N — d |
|
|
|
характеризующая разброс экспериментальных уи отно-
л
сительно уи, предсказанных по уравнению регрессии.
1 Дисперсия вида (15) называется дисперсией неадекватности мо дели; часто обозначается символом
ГОС . ЛVБ.П'
Число степеней свободы |
дисперсии |
s2 равно лц — |
— N— 1, а дисперсии s2 равно |
v z = N —d, |
где d — число |
членов уравнения регрессии |
(4). |
|
Далее составляется отношение дисперсий: |
|
|
F Q= |
s y s l |
(16) |
которое сравнивается с табличным значением критерия Фишера EVi.v.a для выбранного уровня значимости а.
Если критическое значение Е из приложения I не мень ше Е0, утверждается, что в исследованных интервалах переменных зависимость у от х* не описывается получен ным уравнением регрессии. Если Егс:Е0 и число точек N больше числа коэффициентов следует попытаться рас считать уравнение второго порядка относительно выб ранных факторов, вычислить для него дисперсию s^no
формуле (15), отношение F\=s\!s\ и проверить его зна
чимость по Е-критерию. Если окажется, что Е^Етабл, то в качестве модели процесса принимается предпослед няя модель.
Во втором варианте представляется возможным оце нить значимость отдельных коэффициентов регрессии и адекватность модели. Если матрица планирования ор тогональна, то доверительный интервал А6,- устанавлива ется для каждого коэффициента в отдельности:
= ± С | / |
•*{*} = |
(17) |
|
а~\ |
|
где tv.a — табличное значение ^-критерия |
Стыодента |
с v = N —1;
а — уровень значимости (см. приложение II).
Если |A&i|>|6j|, коэффициент bi признается незна чимым (вычеркивается из уравнения регрессии). Если матрица планирования не ортогональна, то доверитель ные границы i-того коэффициента, как и численные зна чения других коэффициентов, взаимосвязаны и необхо дим последовательный статистический анализ (подроб нее см. [3]).
Для проверки адекватности уравнения регрессии со поставляется дисперсия неадекватности s2A[формула
(15)] с дисперсией воспроизводимости э2{г/}:
= |
( 18) |
18
Гипотезе об адекватности уравнения регрессии соот ветствует условие:
F < F V*v,;a>
где Fv .v .a — критическое значение Е-критерия из при
ложения I; а — уровень значимости. Числа степеней свободы:
— N — d, v2 = N — 1.
На практике каждый из N опытов повторяется с ^ 1 раз или с > 1 раз повторяется опыт в центре эксперимента
X i= 0 . В первом случае получаем N раздельных |
оценок |
дисперсий воспроизводимости |
|
С |
|
2 \Уи1— Уи? |
|
$ 0 } = 1=1-------:----- , v = c - ;i . |
(19) |
С — X |
|
Гипотеза об однородности полученных дисперсий про веряется с помощью отношения наибольшей из диспер сий s,,{y} к их сумме по всем опытам:
, _ |
max Sg {t/} |
|
max — |
N |
(20) |
|
2 4 |
M |
|
u=\ |
|
Полученное значение Gm&x сравнивается с критичес |
||
ким значением критерия G |
а |
из приложения III для |
v i= c — 1, V2= N ( c — 1) и уровня значимости а.
Е сл и |
G m a x < G Vl; v.;a, ДИСПерСИИ |
СЧИТЭЮТСЯ ОДНОРОД |
|||
НЫМИ с |
н ад еж н о сть ю |
( 1 — а ) . |
Д л я |
д ал ьн ей ш и х |
р асч ето в |
п р и н и м ается оц ен к а ди сп ер си и в о сп р ои звод и м о сти |
|||||
|
N |
|
|
|
|
|
2 |
S2 { у и ) |
|
|
|
|
s2{«/} = fiL ^ -----, |
v = W(c — 1), |
(21) |
||
т. е. дисперсии усредняются. |
|
|
|
||
Если опыты повторялись |
только в центре экспери |
мента, то оценка дисперсии воспроизводимости
С
2 |yi — у\* |
|
*2 {У} = — ----- :-----, v = N ( c - 1). |
(22) |
с — 1 |
|
2* |
19 |