Файл: Мастеров, В. А. Практика статистического планирования эксперимента в технологии биметаллов.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