v0.5.7.final / Тема 2
.pdf51 Тема 02: Построение эмпирических статистических моделей ХТП |
|
|
|
|
|
|||||||||||
|
M y1 my 2 |
|
|
M y1 my y2 my |
|
|||||||||||
|
|
|
|
|
1 |
|
|
|
|
|
1 |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
M y |
|
|
|
|
|
|
||||
M y |
|
m |
yn |
y m |
|
m |
yn |
y |
|
m |
y2 |
|||||
|
|
n |
|
1 |
|
y1 |
|
n |
|
|
2 |
|
M y m y m |
||||
|
|
|
|
|
|
1 y1 n |
|
|
yn |
|
|
|
|
|
|
|
2 |
|
|
|
M yn myn |
|
В результате матрица дисперсий - ковариаций для экспериментальных значений
yi имеет вид:
|
|
|
|
y2 |
|
|
|
|
|
|
1 |
||
|
|
|
|
|||
|
|
|
COVy2 y1 |
|||
|
|
|||||
COVy |
||||||
|
|
|||||
|
|
|
||||
|
|
|
|
|||
|
|
|
|
|
||
|
|
|
COV |
yn y1 |
||
|
|
|
|
|
Если принять два допущения:
COVy y |
|
COVy y |
|
||
2 |
1 |
2 |
|
1 |
n |
|
|
COVy2 yn |
|||
y2 |
|
||||
|
|
|
|
|
|
|
|
||||
|
|
|
|
2 |
|
COV |
yn y2 |
|
|||
|
|
yn |
|
• |
о независимости измерений COVy y |
0 |
i j |
|
|
i |
j |
|
|
|
|
|
2 |
|
• |
об однородности дисперсии, т.е. несущественном отличии yi |
и их равенстве |
||
|
|
2 |
|
|
|
|
y |
|
|
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |
52 Тема 02: Построение эмпирических статистических моделей ХТП
то получается диагональная матрица дисперсий - ковариаций для измеряемых значений y с одинаковыми дисперсиями y2
COVy y2 E
4.3. Определение оценок дисперсий коэффициентов регрессии
Поскольку a - случайная величина, распределённая по нормальному закону,
ma M a .
По аналогии составим матрицу дисперсий-ковариаций для a:
COVa M a ma a ma T
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |
53 Тема 02: Построение эмпирических статистических моделей ХТП
|
|
2 |
|
|
|
|
|
|
|
|
a0 |
|
COVa0a1 |
COVan am |
|||
|
COVa1a0 |
2 |
COVa1am |
|||||
|
a1 |
|||||||
|
|
|
|
|
|
|
||
|
|
|||||||
|
COV |
|
COV |
|
2 |
|
||
|
|
a a |
0 |
a a |
|
a |
m |
|
|
|
m |
m 1 |
|
|
|
В соответствии с формулой для оценок коэффициентов уравнения регрессии
|
|
|
|
|
|
|
|
ˆ |
э |
m Lm |
|
||||
a Ly |
|
a |
y |
Для определения элементов матрицы дисперсий-ковариаций необходимо подставить два последних выражения в матричную формулу
Если в результате подстановки матрица дисперсий - ковариаций получится диагональной, то коэффициенты регрессии можно считать статистически независимыми
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |
54 Тема 02: Построение эмпирических статистических моделей ХТП
Выполним эту подстановку
|
|
|
|
M |
L |
y |
|
m |
|
|
y |
|
m |
T |
M |
L |
y m |
|
|
|
|
|
|
|
|
|
|||||||||||
COV |
|
|
L |
L |
L |
|
||||||||||||
|
|
|||||||||||||||||
|
a |
|
|
|
y |
|
|
|
|
y |
|
|
|
y |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
M L y my y my |
|
L |
|
||||||||
|
|
|
|
|
T |
|
|
T |
|
|||
|
|
|
|
|
|
|
|
|
||||
|
LM y my y my T L |
T |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|||||||
|
LCOV |
|
|
|||||||||
|
y |
|
L y my |
|
||
|
|
|
T |
|
|
|
|
LT
|
|
|
|
|
|
T |
|
|
|
T |
|
|
1 |
|
T |
|
|
|
T |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Lσ2 |
E L |
σ2 |
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
y |
|
|
|
|
y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
E
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |
55 Тема 02: Построение эмпирических статистических моделей ХТП
T
Поскольку COVy σ2y E , а матрица ( ) 1 симметрична,
|
|
|
|
|
T |
|
|
|
1 |
|
|
|
|
|
|
|
|||
COV σ2 |
|
|
|
|
|||||
|
|
||||||||
|
a y |
|
|
|
|
|
|
Назовём обратную матрицу ( T ) 1 корреляционной матрицей
|
|
|
|
|
|
|
|
|
|
C00 |
C01 |
C0m |
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
T |
|
|
|
|
|
C10 |
C11 |
C1m |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
C |
|
|
|
|
|
. |
||||||
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Cm0 |
Cm1 |
Cmm |
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |
56 Тема 02: Построение эмпирических статистических моделей ХТП
Тогда
|
|
|
|
|
|
C |
|
C |
|
C |
|
|
|
σ2 |
|
COV |
COV |
|
||
|
|
|
|
|
|
|
00 |
|
01 |
|
0m |
|
|
a0 |
|
2 a0a1 |
|
|
a0am |
|
|
|
|
|
σ2 |
|
σ2 C10 |
C11 |
C1m |
|
COVa1a0 |
σa1 |
COVa1am |
||||||||
|
|
|
|
|
||||||||||||||||
COV |
|
|
C |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
a y |
|
y |
|
|
|
||||||||||||||
|
|
|
|
|
|
C |
|
C |
m1 |
C |
|
|
COV |
|
COV |
|
σ2 |
|
|
|
|
|
|
|
|
|
m0 |
|
|
mm |
|
a a |
0 |
a a |
|
a |
m |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
m |
m 1 |
|
|
|
откуда следует:
• для дисперсий коэффициентов регрессии
σ2 |
|
σ2C |
jj |
, j 0, 1, ..., m |
a |
j |
y |
|
|
|
|
|
|
• для ковариаций коэффициентов регрессии
COV |
σ2C |
ji |
, j, i 0, 1, ..., m; |
i j |
a a |
y |
|
|
|
j i |
|
|
|
|
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |
57 Тема 02: Построение эмпирических статистических моделей ХТП
Таким образом, независимость коэффициентов определяется тем, будут ли недиагональные элементы в матрице корреляции C равны нулю.
Значения элементов этой матрицы определяются экспериментальными величинами x и видом функций x , т.е. зависят от того, как поставлен (спланирован) эксперимент.
Вслучае активного эксперимента (например, полного факторного эксперимента – ПФЭ и ортогонального центрального композиционного плана эксперимента - ОЦКП) его проводят так, чтобы матрица C стала диагональной, т.е. чтобы коэффициенты регрессии были статистически независимы.
Вслучае произвольного пассивного эксперимента матрица C оказывается недиагональной и поэтому коэффициенты будут статистически зависимы.
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |
58 Тема 02: Построение эмпирических статистических моделей ХТП
Матрица C называется корреляционной, т.к. с помощью её элементов в соответствии с можно рассчитать корреляции коэффициентов регрессии:
ra a |
|
|
|
C ji |
|
|
|
|
|
||
|
C jjCii |
|
|||
j |
i |
|
|
|
|
|
|
|
|
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |
59Тема 02: Построение эмпирических статистических моделей ХТП
4.4.Определение оценок дисперсии y2
Оценка y2 определяется из экспериментов.
Пусть выходная переменная y зависит от r входных переменных (независимых переменных x ).
Для оценки дисперсии проводятся два типа экспериментов:
•С изменением независимых переменных x1,...xr ;
•Параллельные опыты, когда независимые переменные не меняются.
4.4.1.Определение оценок дисперсий в экспериментах с изменением независимых переменных с различным числом параллельных опытов в каждой точке
а) Определение остаточной дисперсии SR2
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |
60 Тема 02: Построение эмпирических статистических моделей ХТП
S 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
R определяется из экспериментов с изменяющимися значениями |
|||||||||||||||||
(пассивный |
|||||||||||||||||
эксперимент) : |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x, y |
|
x1 |
|
|
|
|
|
|
|
|
|
|
xr |
|
yэ |
||
|
|
|
|
|
|
|
|
|
|
|
|||||||
n |
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
1 k1 |
|
x11 x11 |
|
|
|
|
|
|
|
|
|
x1r x1r |
|
y11э y1эk |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
1 kn |
|
xn1 xn1 |
|
|
|
|
|
|
|
|
|
xnr xnr |
|
ynэ1 ynkэ n |
|||
|
|
|
|
|
|
n |
k |
ˆ |
|
|
2 |
|
|
|
|
||
|
|
|
|
|
|
|
i |
|
э |
|
|
|
|
||||
|
|
|
|
|
yi |
|
yiu |
|
|
SSR |
|
|
|
||||
|
|
S |
2 |
|
|
i 1 u 1 |
|
|
|
|
|
|
|
||||
|
|
R |
|
|
|
n |
|
|
|
|
|
fR |
|
||||
|
|
|
|
|
|
|
ki p |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
i 1
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |