- •3. Статистический анализ результатов эксперимента.
- •1. Найти вид математической модели
- •2. Определить параметры модели
- •3. Оценить точность приближения
- •4. Оценить погрешность в определении параметров
- •1. Найти вид математической модели
- •2. Определить параметры модели
- •4. Оценить погрешность в определении параметров Статистические характеристики модели и эксперимента.
- •Лекция 2
- •3. Оценить точность приближения
- •2.22. Анализ коэффициентов полинома
Лекция 2
3. Оценить точность приближения
4. Оценить погрешность в определении параметров Статистические характеристики модели и эксперимента.
Обычная последовательность реализации метода регрессионного анализа состоит в следующем:
1. Установление наличия зависимости между входами и выходами.
2. Выбор класса функций, которые могли бы описать связь между входами и выходами, если установлено, что она существует(например, полином (2.6).
3. Анализ числа степеней свободы. Первоначальное число степеней свободы модели определяется как разность между числом опытов и числом коэффициентов модели: f = n - m .Чем больше опытов, тем больше информации о модели, эта информация используется для того, чтобы определить ее коэффициенты. Отсюда ясно, что максимальное число коэффициентов модели, которое можно определить на основе проведенного эксперимента, равна числу опытов. Но в
5. Проверка адекватности. Процедура проверки адекватности требует точного знания структуры модели, в нашем случае - порядка полинома и значений его коэффициентов
6. Оценка значимости коэффициентов.
Для проведения статистической обработки эксперимента, к последнему предъявляются определенные требования:
1. Эксперимент должен включать опыты на воспроизводимость для определения дисперсии воспроизводимости. Без определения дисперсии воспроизводимости, строго говоря, нельзя принять никаких решений, входящих в комплекс статистической обработки, в частности нельзя проверить адекватность модели. При строгом выполнении этого требо
2. Выходные величины (
)
должны быть нормально распределенными
случайными величинами. Экспериментатор
практически не имеет возможности
проверить это требование, но обычно
считается, что оно соблюдается.
3. Ошибка в задании входных факторов (
)
должна быть пренебрежимо малой по
сравнению с ошибкой измерения выходных
величин (
).
Это требование обычно в эксперименте
соблюдается.
Адекватность модели эксперименту, число степеней свободы.
Расчет дисперсии воспроизводимости
Для расчета дисперсии воспроизводимости
необходимо провести каждый опыт несколько
раз, например,
раз. В общем случае число повторов в
каждом опыте может быть различным, но
чтобы не усложнять изложение, примем
одинаковое число повторений в каждом
опыте. Тогда таблица эксперимента будет
выглядеть следующим образом:
Таблица 2.3
№ опыта |
X0 |
X1 |
... |
Xm |
y |
y средн. |
|||
1 |
X01 |
X11 |
... |
Xm1 |
y11 |
y12 |
... |
y1l |
y1 |
2 |
X02 |
X12 |
... |
Xm2 |
y21 |
y22 |
... |
y2l |
y2 |
... |
... |
... |
... |
... |
... |
... |
... |
... |
... |
n |
X0n |
X1n |
... |
Xmn |
yn1 |
yn2 |
... |
ynl |
yn |
Сначала рассчитывается среднее значение выходной величины в каждых l повторенных опытах (потратили одну степень свободы):
(2.12)
Далее считается дисперсия воспроизводимости для каждого опыта на основе отклонений значений выхода в каждом опыте от среднего значения в этом опыте:
,
(2.13)
где
- число степеней свободы, с учетом того,
что одна степень свободы потрачена на
вычисление среднего значения.
Далее определяется среднее значение дисперсии воспроизводимости по всем опытам:
(2.14)
Этой дисперсии соответствует число
степеней свободы
-
общее число опытов в эксперименте за
выче
Расчет остаточной дисперсии
Далее определяется остаточная дисперсия на основе отклонений всех экспериментальных значений выхода от расчетных значений:
,
где
-
сумма квадратов отклонений экспериментальных
и расчетных значений,
-
число степеней свободы, соответствующее
этой дисперсии (общее число опытов за
вычетом числа коэффициентов модели).
Расчет дисперсии неадекватности
Имея в виду, что остаточная сумма квадратов отклонений складывается из суммы квадратов отклонений неадекватности и суммы квадратов отклонений воспроизводимости:
можно выделить сумму квадратов неадекватности:
Для определения дисперсии неадекватности рассчитывается число степеней свободы неадекватности по формуле:
Окончательно:
Проведя соответствующие подстановки, получим
(2.15)
Критерий Фишера
.Сопоставление дисперсий неадекватности и дисперсии воспроизводимости проводится по критерию Фишера:
(2.16)
Полученное расчетом значение сравнивается с табличным, определенным для соответствующих степеней свободы. Если F < Fтабл - уравнение адекватно.
