Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

v0.5.7.final / 02

.pdf
Скачиваний:
8
Добавлен:
13.02.2015
Размер:
510.81 Кб
Скачать

11Тема 02: Построение эмпирических статистических моделей ХТП

4.6.2.Качественное и количественное соответствие свойств ММ и

объекта моделирования

Качественное соответствие – это когда тенденции изменения переменных в реальном объекте и ММ совпадают.

При оценке количественного критерия соответствия следует использовать аппарат статистического (в нашем случае – регрессионного) анализа.

Получаемый в результате количественный критерий соответствия не должен компенсировать качественное несоответствие.

Строго говоря, при анализе количественного критерия соответствия должны сравниваться:

yij

- экспериментальные значения случайной величины yij , полученные в

j – ом параллельном опыте i – го эксперимента

с

myˆi - математическим ожиданием рассчитанной по модели величины

выходной переменной yˆi для значений входных переменных xi в i – ом эксперименте.

РХТУ им. Д.И. Менделеева

Кафедра информатики и компьютерного моделирования

12 Тема 02: Построение эмпирических статистических моделей ХТП

Если также ввести среднее значение в i – ом эксперименте

yi и рассчитанное

по ММ (уравнению регрессии) это же значение yˆi , полученное при тех же

величинах входных переменных

xi , что и в эксперименте, дающем

 

yi , то

будет справедливо:

 

 

 

 

 

 

 

 

 

эксп

 

myˆi

 

эксп

эксп

эксп

ˆ

 

ˆ

 

 

yij

 

 

yij

yi

 

yi

yi

 

yi

 

myˆi

 

 

 

 

 

 

 

 

 

 

S 2

 

 

Sa2d

 

 

 

S

2

 

 

 

 

e

 

 

 

 

 

 

 

y

Оценкой первых разностей будет дисперсия воспроизводимости

Se2 , которая

характеризует ошибку экспериментов.

 

 

 

 

 

 

 

Оценкой вторых разностей будет дисперсия адекватности Sad2 , которая характеризует ошибку уравнения (модели) в сравнении с экспериментальными величинами yi (если нет параллельных опытов в каждой экспериментальной точке – это не среднее значение, а просто измеренная величина).

РХТУ им. Д.И. Менделеева

Кафедра информатики и компьютерного моделирования

13 Тема 02: Построение эмпирических статистических моделей ХТП

Оценкой третьего слагаемого является дисперсия рассчитанного значения

выходной переменной (определяемая по аналогии с Se2 и SR2 ).

 

 

n

ˆ

 

myˆ

 

2

 

 

 

 

yi

 

i

 

 

SS y

 

 

i 1

 

 

 

 

 

 

 

S 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

p

 

 

 

 

p

 

 

 

 

 

 

 

 

где p - число значимых коэффициентов уравнения регрессии.

Аппарат дисперсионного анализа указанных выше трёх дисперсий Se2

SR2

S y2

позволяет решить две задачи:

 

 

• оценить адекватность уравнения регрессии ( Sa2d / Se2 ) с использованием критерия Фишера;

определить совместную доверительную область для истинных значений

коэффициентов регрессии

( S y2 / SR2 ).

РХТУ им. Д.И. Менделеева

Кафедра информатики и компьютерного моделирования

14Тема 02: Построение эмпирических статистических моделей ХТП

4.6.3.Оценка совместной доверительной области коэффициентов

регрессии

Отношение дисперсии рассчитанной величины выходной переменной y к остаточной дисперсии S y2 подчиняется распределению Фишера ( F ) с доверительной вероятностью , и условием их малой различимости является:

S 2

табл

 

y

 

 

F f y . , f R

S

2

R

 

 

 

В соответствии с логикой рассматриваемого анализа эти величины должны мало различаться, и граница области, где это условие выполняется, задаётся уравнением:

S y2

табл

 

SSy

 

p

табл

SR2

F f y . , f R

или

SSR

n p F f y . , f R

 

 

 

 

 

 

РХТУ им. Д.И. Менделеева

Кафедра информатики и компьютерного моделирования

15 Тема 02: Построение эмпирических статистических моделей ХТП

Величина

n

SSR yˆi yi 2

i1

-значение критерия, полученное при реализации программы минимизации.

Значение SS y

ˆ

SS y y

можно заменить матричным произведением:

T ˆ

 

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

myˆ y myˆ

a

ˆ

T

 

 

T

 

 

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

a

 

a

a a a

T

 

 

 

ˆ

 

 

 

 

a ,

 

 

 

так как в соответствии с эмпирическим и теоретическим уравнением регрессии для всех расчетных значений выходной переменной будут справедливы следующие матричные равенства:

РХТУ им. Д.И. Менделеева

Кафедра информатики и компьютерного моделирования

16 Тема 02: Построение эмпирических статистических моделей ХТП

 

 

 

 

 

 

 

 

 

 

 

 

 

a

ˆ ˆ

и

myˆ

y a

 

Подстановкой матричного произведения вместо SS y , получается квадратичная форма вида:

ˆ

T

 

 

T

 

ˆ

p табл.

 

 

 

 

 

 

 

 

 

 

a a

a a SSR

 

Fβf y . , fR

n p

Геометрической интерпретацией этой квадратичной формы является эллипсоид, оси которого пропорциональны собственным значениям матрицы

A T ,

определяемой из характеристического уравнения вида:

A E 0.

РХТУ им. Д.И. Менделеева

Кафедра информатики и компьютерного моделирования

17 Тема 02: Построение эмпирических статистических моделей ХТП

Для двух коэффициентов a0 и a1 , получается эллипс вида:

a1

aˆ1

ˆ

a

 

a0

0

 

 

Получена совместная доверительная область для коэффициентов (здесь a0 и a1) в линейной модели. Её можно сравнить с прямоугольником, образованным

оценками доверительных интервалов коэффициентов регрессии

РХТУ им. Д.И. Менделеева

Кафедра информатики и компьютерного моделирования

18 Тема 02: Построение эмпирических статистических моделей ХТП

Длинная, вытянутая доверительная область (собственные значения матрицы A существенно отличаются) указывает на то, что коэффициенты сильно коррелированны и значения коэффициентов плохо оценены.

Результатом высокой корреляции коэффициентов является то, что неправильно оцененное значение одного из коэффициентов можно сбалансировать при проведении подгонки компенсирующим исправленным значением другого параметра.

Подгонка даёт почти столь же хороший результат, как и результат, получаемый при использовании наилучших оценок. Поверхность суммы квадратов Crзадаётся уравнением:

 

 

 

p

табл.

 

 

 

 

 

Cr SSR SSy

SSR 1

 

 

Fβ f y , f R

n p

 

 

 

 

 

РХТУ им. Д.И. Менделеева

Кафедра информатики и компьютерного моделирования

Соседние файлы в папке v0.5.7.final