Министерство науки и высшего образования РФ
Федеральное государственное бюджетное образовательное учреждение
высшего образования
«Сибирский государственный автомобильно-дорожный университет»
(СибАДИ)
Кафедра «Управление качеством и производственными системами»
КУРСОВАЯ РАБОТА ПО ДИСЦИПЛИНЕ
«Основы теории и планирования эксперимента»
На тему «Применение полного факторного эксперимента при проведении
исследований»
Выполнил(а):
Студентка группы УКб-15Э1
Граматчикова В. Е.
Проверила: к.э.н., доцент
Семенова Е.С.
Омск 2019
Министерство образования и науки РФ
Федеральное государственное бюджетное образовательное учреждение
высшего образования
«Сибирская государственная автомобильно-дорожный университет»
(СибАДИ)
Кафедра «Управление качеством и производственными системами»
ЗАДАНИЕ
на курсовую работу по теме:
«Применение полного факторного эксперимента при проведении исследований»
Студенту группы УКб-15Э1 Граматчиковой В.Е.
1. Срок сдачи выполненной курсовой работы до 15.01.2019 г.
2. Консультант к.э.н., доц. Семёнова Е.С.
3. Исходные данные к курсовой работе: учебная и научная литература российских и
зарубежных
авторов,
действующие
нормативные
документы,
касающиеся
рассматриваемой темы, ресурсы всемирной сети Internet.
4.Обязательный перечень содержания работы:
ВВЕДЕНИЕ
1 ОПРЕДЕЛЕНИЕ ЦЕЛЕЙ И ЗАДАЧ ЭКСПЕРИМЕНТА
1.1 Цели и задачи полного многофакторного эксперимента
1.2 Методика проведения эксперимента по определению устойчивости к
многократному изгибу материалов для средств защиты рук
2 ОПРЕДЕЛЕНИЕ ОТНОСИТЕЛЬНОЙ ПОГРЕШНОСТИ РЕЗУЛЬТАТОВ
ИЗМЕРЕНИЯ
2.1 Погрешности измерений
2.2 Определение относительной погрешности результатов измерений разрывной
нагрузки материалов для средств защиты рук
3 ПРОВЕДЕНИЕ КОРРЕЛЯЦИОННОГО АНАЛИЗА ЭКСПЕРИМЕНТАЛЬНЫХ
ДАННЫХ
3.1 Понятие корреляционного анализа
3.2 Корреляционная связь
3.3 Определение коэффициента корреляции
3.4 Проверка гипотезы о значимости коэффициента корреляции
4 ОБРАБОТКА РЕЗУЛЬТАТОВ МНОГОФАКТОРНОГО ЭКСПЕРИМЕНТА
4.1 Кодирование переменных
4.2 Построение матрицы планирования в кодированных переменных
4.3 Нахождение коэффициентов уравнения регрессии
4.4 Проверка коэффициентов на значимость
4.5 Проверка уравнения на адекватность
4.6 Интерпретация полученной модели
ЗАКЛЮЧЕНИЕ
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
Задание выдано
"12 " сентября
2018г.
(подпись)
Задание к исполнению принял
"12 " сентября
2018г.
(подпись)
3
СОДЕРЖАНИЕ
ВВЕДЕНИЕ...................................................................................................... 4
1 ОПРЕДЕЛЕНИЕ ЦЕЛЕЙ И ЗАДАЧ ЭКСПЕРИМЕНТА ....................... 6
1.1 Цели и задачи многфакторного эксперимента ...................................... 6
1.2 Методика проведения эксперимента по определению устойчивости к
многократному изгибу материалов для средств защиты рук ............................. 9
2 ОПРЕДЕЛЕНИЕ ОТНОСИТЕЛЬНОЙ ПОГРЕШНОСТИ
РЕЗУЛЬТАТОВ ИЗМЕРЕНИЯ ........................................................................... 12
2.1 Погрешности измерений ........................................................................ 12
2.2 Определение относительной погрешности результатов измерений
разрывной нагрузки материалов для средств защиты рук ................................ 16
3 ПРОВЕДЕНИЕ КОРРЕЛЯЦИОННОГО АНАЛИЗА
ЭКСПЕРИМЕНТАЛЬНЫХ ДАННЫХ ............................................................... 18
3.1 Понятие корреляционного анализа ...................................................... 18
3.2 Корреляционная связь ............................................................................ 20
3.3 Определение коэффициента корреляции ............................................. 24
3.4 Проверка гипотезы о значимости коэффициента корреляции ........... 26
4 ОБРАБОТКА РЕЗУЛЬТАТОВ МНОГОФАКТОРНОГО
ЭКСПЕРИМЕНТА ............................................................................................... 28
4.1 Кодирование переменных. ..................................................................... 28
4.2 Построение матрицы планирования в кодированных переменных .. 30
4.3 Нахождение коэффициентов уравнения регрессии............................. 32
4.4 Проверка коэффициентов на значимость ............................................ 33
4.5 Проверка уравнения на адекватность ................................................... 35
4.6 Интерпретация полученной модели...................................................... 37
ЗАКЛЮЧЕНИЕ ............................................................................................. 38
БИБЛИОГРАФИЧЕСКИЙ СПИСОК ......................................................... 39
4
ВВЕДЕНИЕ
Описание
поведения
и
прогнозирование
состояния
сложных
технических объектов затруднено отсутствием достаточного количества
математических моделей. В ряде случаев простым выходом из ситуации
является использование регрессионных закономерностей, получаемых в ходе
планирования эксперимента. К сожалению, регрессионные закономерности
не дают представлений о типе протекающих в системе процессов и не могут
быть использованы вне области планирования эксперимента. При этом их
достоинством является способность описания любого объекта при полном
соблюдении
алгоритма
активного
эксперимента.
Полный
факторный
эксперимент является наиболее легко реализуемым среди многочисленных
методов активного эксперимента.
Целью выполнения курсовой работы является изучение и применение
методов
планирования
полного
факторного
эксперимента
(ПФЭ)
применительно к технологическим задачам, освоение принципов составления
матрицы планирования ПФЭ, проведение расчета коэффициентов регрессии,
использование
статистических
критериев
для
оценки
однородности,
нормальности экспериментальных данных, значимости коэффициентов и
адекватности полученной математической модели.
Поставленная цель обусловила решение ряда задач:
1. Исследование теоретических основ и планирование эксперимента по
определению качественных характеристик материалов для средств защиты
рук, применительных в Северных регионах при решении технологических
задач.
2. Определение относительной погрешности результатов измерения.
3. Проведение корреляционного анализа экспериментальных данных.
4. Планирование ПФЭ, выбор числа и условий проведения опытов,
необходимых и достаточных для получения математической модели
процесса.
5
5. Интерпретация полученной математической модели при влиянии
факторов на устойчивость к многократному изгибу материалов для защиты
рук.
Объектом в курсовой работе выступают материалы для средств защиты
рук. Для защиты рук промышленностью выпускаются различные перчатки.
Средства защиты рук следует выбирать исходя из их назначения и свойств,
вида выполняемой работы и требуемой степени защиты. В данном случае
рассматривались средства защиты рук, рекомендуемых для работ на окрытом
воздухе в условиях Севера и Крайнего Севера (диапазон рабочих температур
от -60 С° до 100 С°) [11].
Предметом являются методы планирования ПФЭ.
Информационной базой для написания курсовой работы послужили
труды отечественных ученых, законодательная и справочная документация,
периодические издания и электронные ресурсы.
6
1 ОПРЕДЕЛЕНИЕ ЦЕЛЕЙ И ЗАДАЧ ЭКСПЕРИМЕНТА
1.1 Цели и задачи многофакторного эксперимента
Чаще всего эксперимент ставят для решения одной из двух основных
задач. Первую задачу называют экстремальной. Она заключается в
отыскании условий процесса, обеспечивающих получение оптимального
значения выбранного параметра. Признаком экстремальных задач является
требование поиска экстремума некоторой функции. Эксперименты, которые
ставят для решения задач оптимизации, называют экстремальными. Вторую
задачу
называют
интерполяционной
интерполяционной.
формулы
для
Она
состоит
предсказаний
в
значений
построении
изучаемого
параметра, зависящего от ряда факторов. Для решения экстремальной или
интерполяционной задачи необходимо иметь математическую модель
исследуемого объекта. Модель объекта получают, используя результаты
опытов. При исследовании многофакторного процесса постановка всех
возможных опытов для получения математической модели связана с
огромной трудоемкостью эксперимента, так как их число очень велико.
Задача планирования состоит в установлении минимально необходимого
числа экспериментов и условий их проведения, в выборе методов
математической обработки результатов и в принятии решений. Частным
случаем планирования эксперимента является планирование экстремального
эксперимента, т. е. процесс выбора их числа и условий проведения,
минимально необходимых для нахождения экстремальных экспериментов с
помощью
метода
Бокса
-
Уилсона,
называемого
методом
крутого
восхождения.
Метод Бокса - Уилсона предусматривает проведение экспериментов
небольшими сериями. В каждой серии одновременно варьируют все факторы
по определенным правилам. Эксперименты проводят так, чтобы после
математической обработки результатов предыдущей серии можно было
спланировать следующую серию.
7
При планировании экстремального эксперимента цель исследования
должна быть четко сформулирована и должна иметь количественную оценку.
Характеристику цели, заданную количественно, называют параметром
оптимизации. Параметр оптимизации является реакцией, или откликом, на
воздействие факторов, определяющих поведение процесса. Результаты
эксперимента
используют
для
получения
математической
модели
исследуемого процесса.
Математическая модель - система математических соотношений,
описывающих
изучаемый
процесс
или
явление.
При
планировании
эксперимента под математической моделью часто понимают уравнение,
связывающее параметр оптимизации с факторами. Такое уравнение
называют функцией отклика [3].
Планирование эксперимента — это процедура выбора числа и условий
проведения
опытов,
необходимых
и
достаточных
для
получения
математической модели процесса. При этом важно учитывать следующее:
стремление к минимизации числа опытов; одновременное варьирование всех
переменных, определяющих процесс; выбор четкой стратегии, позволяющей
принимать обоснованные решения после каждой серии экспериментов.
При использовании метода активного планирования весь эксперимент
обычно разбивается на несколько этапов. Информация, полученная после
каждого этапа, используется для планирования исследований на следующем
этапе. Планирование эксперимента позволяет варьировать ряд факторов и
получать одновременно количественные оценки всех проявляющихся
эффектов. При этом в отличие от классического регрессионного анализа,
избежать корреляции между коэффициентами уравнения регрессии. При
статистическом подходе математическая модель объекта или процесса
представляется в общем виде полиномом n-степени, т.е. отрезком ряда
Тейлора, в который разлагается неизвестная функция [4].
Под планированием эксперимента понимают процесс определения числа
и условий проведения опытов, необходимых и достаточных для решения
8
поставленной задачи с требуемой точностью. Эффективность использования
МПЭ при исследовании процессов в различных объектах объясняется тем,
что их наиболее важные характеристики являются случайными величинами,
распределение которых близко к нормальному.
Основная задача исследований — оптимизация, заключающаяся в
нахождении совокупности варьируемых параметров, при которых выбранная
целевая
функция
(параметр
оптимизации)
принимает
экстремальное
значение.
При решении задачи предполагается, что оптимизируется один
параметр, каждый из факторов управляем, результаты опытов должны
воспроизводиться. Для построения математической модели используется
полный
или
дробный
факторный
план
эксперимента,
обладающий
оптимальной матрицей планирования.
При планировании эксперимента цель исследования должна быть четко
сформулирована и должна иметь количественную оценку. Характеристику
цели, заданную количественно, называют параметром оптимизации.
Границы изменения факторов определяются так, чтобы обеспечить
условия
физической
управляемыми,
реализации
независимыми
процесса.
Факторы
(устанавливаться
на
должны
любом
быть
уровне
независимо друг от друга), однозначными (не являться функциями других).
На выбор интервала варьирования так же накладываются ограничения:
он не может быть меньше ошибки, с которой экспериментатор фиксирует
уровень фактора, и не может быть настолько большим, что верхний и
нижний уровень оказались за пределами области определения. Необходимо
отметить, что при выборе верхнего и нижнего уровней факторов необходимо
учитывать ограничения, связанные с свойствами объекта исследования.
Полный факторный эксперимент (2k) целесообразно проводить в том
случае, если он не продолжителен по времени и требует небольших затрат
(число варьируемых факторов k мало). Для каждого включенного в
эксперимент факоров устанавливают только два уровня: верхний и нижний.
9
Для полного факторного эксперимента, в котором реализуются все
возможные сочетания уровней факторов, число опытов определяется по
следующей формуле: N=2k, где N - число опытов; k - число факторов.
Эксперимент, в котором реализуются все возможные сочетания уровней
факторов, называется полным факторным экспериментом.
Вводят условное обозначение верхнего (+), нижнего (-) и основного
уровня (0).Затем строится план эксперимента (например, для трех факторов,
приведенный в таблице 1.1):
Таблица 1.1 – Матрица планирования
Номер точки
плана
1
2
3
4
5
6
7
8
Х1
+
+
+
+
-
Факторы
Х2
+
+
+
+
-
Значение параметра
оптимизации
Y1
Y2
Y3
Y4
Y5
Y6
Y7
Y8
Х3
+
+
+
+
-
Составленная матрица планирования должна соответствовать свойствам
полного факторного эксперимента:
1. Симметричность относительно центра эксперимента [5]:
,
(1)
где j - номер фактора; N - число опытов, j =1... k.
2. Условие нормировки:
.
(2)
3. Ортогональность матрицы:
,
(3)
j=u; j, u=0, 1, 2 … k.
1.2 Методика проведения эксперимента по определению устойчивости к
многократному изгибу материалов для средств защиты рук
Для определения номенклатуры показателей качества материалов для
10
средств защиты рук использовался ГОСТ 4.493-89 «Система показателей
качества продукции (СПКП). Материалы для средств защиты рук.
Номенклатура показателей».
Проводилось исследование некоторых факторов на устойчивость к
многократному изгибу материалов для средств защиты рук. На устойчивость
к многократному изгибу влияют такие факторы как [1]:
1. толщина, мм;
2. разрывная нагрузка, Н;
3. удлинение при разрыве, %.
Определение устойчивости к многократному изгибу материалов для
средств защиты рук осуществлялось по ГОСТ 8978-75 «Кожа искусственная
и
пленочные
материалы.
Методы
определения
устойчивости
к
многократному изгибу».
Устойчивость к многократному изгибу характеризуется устойчивостью
к образованию нарушений на лицевом покрытии или основе.
Существует три метода определения устойчивости к разрушению
искусственной кожи:
1. Метод определения устойчивости к разрушению искусственной кожи
при сжатии и изгибе образца цилиндрической формы вдоль оси цилиндра.
2. Метод определения устойчивости к разрушению искусственной кожи
при сжатии и изгибе образца по форме ромба.
3. Метод определения устойчивости к разрушению искусственной кожи
и пленочного материала при изгибе вокруг зажимов.
Первый метод предназначен
разрушению
искусственной
цилиндрической
формы
для
кожи
вдоль
определения
при
оси
распрямлением или растяжением.
Для проведения испытания применяют:
1. Прибор типа МИРЦ.
сжатии
цилиндра,
устойчивости
и
изгибе
к
образца
чередующихся
с
11
2. Раму для предварительного одноосного растяжения образцов с
зажимами, размещенными попарно в одной плоскости и параллельными
между собой.
3. Линейку металлическую по ГОСТ 427-75.
Образцы искусственной кожи растягивают на раме в направлении
длины образца. Зажимы рамы устанавливают на расстоянии (140±1) мм,
образец складывают втрое по ширине, закрепляют в зажимах и растягивают
на определенную величину. Значение этой величины, равное 75% от
удлинения при разрыве, вычисляют в миллиметрах по отношению к длине
образца между зажимами рамы.
Второй
метод
предназначен
для
определения
устойчивости
к
разрушению искусственной кожи при сжатии и изгибе образца по форме
ромба,
чередующихся
с
распрямлением
или
растяжением
образца.
Измеряют толщину образцов в миллиметрах в трех точках по средней
поперечной линии с погрешностью не более 0,01 мм по ГОСТ 17073-71. Для
режимов испытания с растяжением и без растяжения рассчитывают
минимальный зазор между зажимами, равный шести толщинам.
Третий
метод
предназначен
для
определения
устойчивости
к
разрушению искусственной кожи или пленочного материала при изгибе
вокруг зажимов с постоянным радиусом кривизны на заданный угол в
каждую сторону от вертикального положения образца (двойной перегиб),
находящегося при постоянно действующей растягивающей нагрузке.
Измеряют толщину образца в миллиметрах по месту изгиба на границе
между первым и вторым участками.
Результат каждого измерения устойчивости искусственной кожи и
пленочного материала к многократному изгибу выражают числом циклов
испытания в килоциклах с погрешностью не более 0,1% от измеренной
величины [2].
12
2 ОПРЕДЕЛЕНИЕ ОТНОСИТЕЛЬНОЙ ПОГРЕШНОСТИ РЕЗУЛЬТАТОВ
ИЗМЕРЕНИЯ
2.1 Погрешности измерений
Измерения не могут быть выполнены абсолютно точно. Всегда имеется
некоторая неопределенность в значении измеряемой величины. Эта
неопределенность
характеризуется
погрешностью.
Погрешность
измерения — отклонение измеренного значения величины от её истинного
(действительного)
значения.
Погрешность
измерения
является
характеристикой точности измерения.
Причины, приводящие к появлению погрешностей:
1. Ограниченная точность измерительных приборов.
2. Влияние на измерение неконтролируемых изменений внешних
условий (напряжения в электрической сети, температуры и т.д.).
3. Действия экспериментатора (включение секундомера с некоторым
запаздыванием, различное размещение глаз по отношению к шкале прибора
и т.п.).
4. Неполное соответствие измеряемого объекта той абстракции, которая
принята для измеряемой величины (например, при измерении объема
пластинка считается параллелепипедом, в то время как у нее могут быть
закругления на ребрах).
5. Нестрогость
законов,
которые
используются
для
нахождения
измеряемой величины или лежат в основе устройства прибора.
Выяснить с абсолютной точностью истинное значение измеряемой
величины, как правило, невозможно, поэтому невозможно и указать
величину отклонения измеренного значения от истинного. Это отклонение
принято называть ошибкой измерения. Возможно лишь оценить величину
этого отклонения, например, при помощи статистических методов. На
практике вместо истинного значения используют действительное значение
величины,
то
есть
значение
физической
величины,
полученное
13
экспериментальным путём и настолько близкое к истинному значению, что в
поставленной измерительной задаче может быть использовано вместо него.
Такое значение, обычно, вычисляется как среднестатистическое значение,
полученное при статистической обработке результатов серии измерений. Это
полученное значение не является точным, а лишь наиболее вероятным.
Поэтому в измерениях необходимо указывать, какова их точность. Для этого
вместе с полученным результатом указывается погрешность измерений.
Задачей эксперимента является определение числового значения
измеряемых величин с заданной точностью. Эта задача решается с
помощью прямых или косвенных измерений.
При прямом
физической
измерении осуществляется
величины
измерительных
с
количественное
соответствующим
приборов.
Отсчет
по
эталоном
шкале
при
прибора
сравнение
помощи
указывает
непосредственно измеряемое значение.
При косвенных
измерениях
физическая
величина находится
при
помощи математических операций над непосредственно измеренными
физическими величинами.
Точность
величиной
прямых
измерений некоторой
погрешности
или
ошибки,
величины X оценивается
измерений
относительно
действительного значения физической величины.
Действительное значение величины (
) – это значение физической
величины, полученное экспериментальным путем и настолько близкое к
истинному значению, что в поставленной измерительной задаче может быть
использовано вместо него.
По форме представления выделяют абсолютную (
) и относительную
(д) погрешности измерений.
Абсолютная
погрешность
измерения – это
погрешность
средства
измерений, выраженная в единицах измеряемой физической величины,
характеризующая
абсолютное
отклонение
измеряемой
действительного значения физической величины:
величины
от
14
–
.
(4)
Действительное значение находится по формуле:
,
(5)
где n – число измерений.
Относительная погрешность измерения – это погрешность измерения,
выраженная
отношением
абсолютной
погрешности
измерения
к
действительному значению измеряемой величины. Обычно относительную
погрешность выражают в процентах:
.
(6)
Если при измерениях получена относительная погрешность более 10%,
то говорят, что произведено не измерение, а лишь оценка измеряемой
величины.
В
лабораториях
физического
практикума
относительная
погрешность обычно составляет 1-10%. В научных же лабораториях
измерения некоторых физических величин, таких, например, как длина
световой волны, осуществляется с точностью порядка миллионной доли
процента.
Основная погрешность средства измерений — погрешность средства
измерений, применяемого в нормальных условиях.
Дополнительная погрешность средства измерений — составляющая
погрешности средства измерений, возникающая дополнительно к основной
погрешности вследствие отклонения какой-либо из влияющих величин от
нормального ее значения или вследствие ее выхода за пределы нормальной
области значений.
Стабильность средства измерений — качественная характеристика
средства
измерений,
метрологических
отражающая
характеристик
неизменность
(в
качестве
во
времени
его
количественной
оценки стабильности служит нестабильность средства измерений).
Нестабильность средства измерений — изменение метрологических
характеристик средства измерений за установленный интервал времени.
15
Класс точности средств измерений —
обобщенная характеристика
данного типа средств измерений, как правило, отражающая уровень их
точности, выражаемая пределами допускаемых основной и дополнительных
погрешностей, а также другими характеристиками, влияющими на точность.
Также, нужно учитывать погрешности по характеру проявления, такие
как:
Случайная погрешность — составляющая погрешности измерения,
изменяющаяся случайным образом в серии повторных измерений одной и
той же величины, проведенных в одних и тех же условиях. В появлении
таких погрешностей не наблюдается какой-либо закономерности, они
обнаруживаются при повторных измерениях одной и той же величины в виде
некоторого разброса получаемых результатов. Случайные погрешности
неизбежны, неустранимы и всегда присутствуют в результате измерения,
однако их влияние обычно можно устранить статистической обработкой.
Описание случайных погрешностей возможно только на основе теории
случайных процессов и математической статистики.
Систематическая погрешность — погрешность, изменяющаяся во
времени по определённому закону (частным случаем является постоянная
погрешность, не изменяющаяся с течением времени). Систематические
погрешности могут быть связаны с ошибками приборов (неправильная
шкала,
калибровка
и т. п.),
неучтёнными
экспериментатором.
Систематическую ошибку нельзя устранить повторными измерениями. Её
устраняют либо с помощью поправок, либо «улучшением» эксперимента.
Прогрессирующая (дрейфовая) погрешность —
погрешность,
медленно
меняющаяся
во
времени.
непредсказуемая
Обусловлена
она
нарушениями статистической устойчивости.
Грубая погрешность (промах) — погрешность, возникшая вследствие
недосмотра экспериментатора или неисправности аппаратуры (например,
если экспериментатор неправильно прочёл номер деления на шкале прибора
или если произошло замыкание в электрической цепи) [4].
16
2.2 Определение относительной погрешности результатов измерения
разрывной нагрузки материалов для средств защиты рук
Результаты прямых многократных измерений разрывной нагрузки
материалов тканей для средств защиты рук приведены в табл. 2.1.
Таблица 2.1 - Результаты прямых многократных измерений разрывной нагрузки и
устойчивости к многократному изгибу материалов для средств защиты рук
Относительная
Х
Абсолютная погрешность,
погрешность,
№
(разрывная
нагрузка, Н)
1
0,6
-0,56
-48,17
2
1
-0,16
-13,62
3
1,1
-0,06
-4,98
4
1,25
0,09
7,98
5
1,2
0,04
3,66
6
1,15
-0,01
-0,66
7
1,3
0,14
12,29
8
0,95
-0,21
-17,94
9
1,08
-0,08
-6,71
10
1,18
0,02
1,93
11
1,21
0,05
4,52
12
1,19
0,03
2,79
13
1,25
0,09
7,98
14
1,23
0,07
6,25
15
1,18
0,02
1,93
16
1,45
0,29
25,25
17
1,18
0,02
1,93
18
1,35
0,19
16,61
19
1,26
0,10
8,84
20
1,2
0,04
3,66
21
1,16
0,00
0,20
22
1,24
0,08
7,11
23
1,1
-0,06
-4,98
24
0,96
-0,20
-17,07
25
1,21
0,05
4,52
26
0,98
-0,18
-15,35
27
1,19
0,03
2,79
28
1,22
0,06
5,38
29
1,1
-0,06
-4,98
30
1,26
0,10
8,84
За действительный размер принимается среднее арифметическое при
многократных измерениях. Действительный размер рассчитывается по
формуле (5):
17
Н.
Для
каждого
ряда
наблюдений,
приведенных
в
табл.
2.1,
рассчитывается абсолютная погрешность с помощью формулы (4):
Н;
Н;
Н;
Н;
Н;
Н;
Н;
Н.
Аналогичным образом была рассчитана абсолютная погрешность для
наблюдений №5-30. Результаты представлены в табл. 2.1 (3 столбец).
Относительная погрешность рассчитывается по формуле (6):
;
;
;
;
;
;
;
.
Аналогичным образом была рассчитана относительная погрешность для
наблюдений 9-30. Результаты представлены в табл. 2.1 (4 столбец).
18
3 ПРОВЕДЕНИЕ КОРРЕЛЯЦИОННОГО АНАЛИЗА
ЭКСПЕРИМЕНТАЛЬНЫХ ДАННЫХ
3.1 Понятие корреляционного анализа
Зависимости между величинами (факторами, признаками) разделяют на
два вида: функциональную и статистическую.
Отличительной
особенностью
функциональной
зависимости
двух
величин считается тот факт, что значению одной из них всегда соответствует
одно или несколько точно определённых значений другой величины. Кроме
того, функциональная связь двух факторов возможна только при условии,
что вторая величина зависит только от первой и не зависит ни от каких
других величин. В случае зависимости величины от множества факторов,
функциональная связь возможна, если первая величина не зависит ни от
каких других факторов, кроме входящих в указанное множество.
В реальных ситуациях эти условия не выполнимы, поскольку всегда
существует огромное количество свойств самого объекта и окружающей его
среды, влияющих друг на друга, что приводит к невозможности учёта всех
взаимодействий (как качественно, так и количественно), поэтому такого рода
связи
не
существуют.
Естественно,
функциональная
связь
является
математической абстракцией и находит практическое применение только
тогда, когда определённая величина в основном зависит от соответствующих
факторов, а остальные связи пренебрежимо слабы.
При статистической зависимости изменение одной из величин влечёт
изменение распределения других величин, которые с определенными
вероятностями
принимают
некоторые
значения.
Функциональную
зависимость следует считать частным случаем статистической: значению
одного фактора соответствуют значения других факторов с вероятностью,
равной единице.
Значительно больший интерес представляет другой частный случай
статистической зависимости, когда существует взаимосвязь значений одних
19
случайных величин со средним значением других, при той особенности, что
в каждом отдельном случае любая из взаимосвязанных величин может
принимать различные значения.
Такого рода зависимость между переменными величинами называется
корреляционной, или корреляцией.
Корреляционный анализ — метод обработки статистических данных,
заключающийся в изучении коэффициентов (корреляции). Его применение
возможно в случае наличия достаточного количества (для конкретного вида
коэффициента корреляции) наблюдений из более чем одной переменной. При
этом сравниваются коэффициенты корреляции между одной парой или
множеством пар признаков, для установления между ними статистических
взаимосвязей.
Корреляционный анализ, как и другие статистические методы, основан
на
использовании
вероятностных
моделей,
описывающих
поведение
исследуемых признаков в некоторой генеральной совокупности, из которой
получены экспериментальные значения
и
.
Задача корреляционного анализа сводится к следующим более малым
задачам:
установлению формы (линейная, нелинейная);
направления (положительное или отрицательное) связи между
варьирующими признаками;
измерению тесноты связи;
и
проверке
уровня
значимости
полученных
коэффициентов
корреляции.
Популярность метода обусловлена двумя моментами: коэффициенты
корреляции относительно просты в подсчете, их применение не требует
специальной
математической подготовки. В сочетании
с простотой
интерпретации (принятие гипотезы о наличии корреляции означает что
изменение переменной А, произойдет одновременно с изменением значения
20
Б),
простота
применения
коэффициента
привела
к
его
широкому
распространению в сфере анализа статистических данных.
Когда исследуется корреляция между количественными признаками,
значения которых можно точно измерить в единицах метрических шкал, то
очень часто принимается модель двумерной нормально распределенной
генеральной совокупности. Такая модель отображает зависимость между
переменными величинами
и
графически в виде геометрического места
точек в системе прямоугольных координат. Изображенные на координатной
плоскости точки (xi, yi), где xi и yi – значения первой и второй переменных,
называются корреляционным полем (рис. 1а). Аналитическая функция,
аппроксимирующая
(приближенно
описывающая)
наблюдаемые
эмпирические значения, называется функцией регрессии (рис. 1б).
Рисунок 1 – Корреляционное поле (а) и функция регрессии (б)
Функция регрессии отражает тенденцию изменения одной величины
под действием другой и строится таким образом, чтобы эмпирические точки
корреляционного поля лежали как можно ближе к ней. Функция регрессии
может быть линейной, параболической, гиперболической, логарифмической
и др [6].
3.2 Корреляционная связь
По направлению корреляционная связь может быть положительной
(прямой) и отрицательной (обратной).
При положительной линейной корреляции более высоким значениям
21
одного признака соответствуют более высокие значения другого, а более
низким значениям одного признака - более низкие значения другого. При
отрицательной корреляции соотношения обратные.
Знак коэффициента корреляции зависит от направления корреляционной
связи: при положительной корреляции коэффициент корреляции имеет
положительный знак, при отрицательной корреляции - отрицательный знак.
При изучении корреляционной связи важным направлением анализа
является оценка степени тесноты связи. Понятие степени тесноты связи
между двумя признаками возникает вследствие того, что в действительности
на изменение результативного признака влияет множество факторов. При
этом влияние одного из факторов может выражаться более заметно и четко,
чем влияние других факторов. С изменением условий роль решающего
фактора может перейти к другому признаку.
При статистическом изучении взаимосвязей, как правило, учитываются
только основные факторы. Также с учетом степени тесноты связи
оценивается необходимость более подробного изучения конкретной данной
связи и значение практического ее использования.
В общем, знание количественной оценки тесноты корреляционной связи
позволяет решить следующую группу вопросов:
необходимость глубокого изучения данной связи между признаками и
целесообразность ее практического применения;
степень различий в проявлении связи в конкретных условиях
(сопоставление оценки тесноты связи для различных условий);
выявление главных и второстепенных факторов в данных конкретных
условиях путём последовательного рассмотрения и сравнения признака с
различными факторами.
Показатели тесноты связи должны удовлетворять ряду основных
требований:
22
величина показателя тесноты связи должна быть равна или близка к
нулю, если связь между изучаемыми признаками (процессами, явлениями)
отсутствует;
при наличии между изучаемыми признаками функциональной связи
величина показателя тесноты связи должна быть равна единице;
при наличии между признаками корреляционной связи абсолютное
значение показателя тесноты связи должно выражаться правильной дробью,
которая по величине тем больше, чем теснее связь между изучаемыми
признаками (стремится к единице).
Корреляционная зависимость определяется различными параметрами,
среди которых наибольшее распространение получили парные показатели,
характеризующие взаимосвязь двух случайных величин: коэффициент
ковариации (корреляционный момент) и линейный коэффициент корреляции
(коэффициент корреляции Пирсона).
Сила связи определяется абсолютным значением показателя тесноты
связи и не зависит от направления связи.
В зависимости от абсолютного значения коэффициента корреляции
корреляционные связи между признаками по силе делятся следующим
образом:
сильная, или тесная при коэффициенте корреляции
средняя (при 0,50<
<0,69);
умеренная (при 0,30<
слабая (при 0,20<
очень слабая (при
>0,70;
<0,49);
<0,29);
<0,19).
Простейшей, но достаточно информативной характеристикой тесноты
связи двух величин х и у является коэффициент корреляции. В теории
вероятностей
этот
показатель
определяется
с
помощью
других
вероятностных характеристик.
Наилучшим приближенным значением коэффициента корреляции (или
23
коэффициент Пирсона), который можно вычислить с помощью результатов
измерений, является величина rxy, определяемая по формуле:
,
(7)
где
;
(8)
;
(9)
.
(10)
Для расчета выборочного коэффициента корреляции нужно знать
выборочные средние для х и у, которые рассчитываются по следующим
формулам:
,
(11)
.
(12)
Значение коэффициента корреляции лежит в пределах от –1 до +1.
Если случайные величины x, y независимы, то коэффициент корреляции
обязательно равен нулю, обратное утверждение неверно.
Коэффициент корреляции характеризует значимость линейной связи
между параметрами:
при
при
при
значения
величины
величины
и полностью совпадают;
и
принимают противоположные значения;
и
практически не связаны друг с другом
линейным соотношением. Однако это не означает отсутствия каких-то
других (например, нелинейных) связей между параметрами;
при
однозначной линейной связи величин
и
нет, и
чем меньше абсолютная величина коэффициента корреляции, тем в меньшей
24
степени по значениям одного параметра можно предсказать значение
другого.
В первых двух случаях имеет место функциональная зависимость: зная
значение одного параметра, можно однозначно указать значение другого
параметра.
3.3 Определение коэффициента корреляции
Значение коэффициента корреляции не зависит от того, какую из двух
переменных назвать х, а какую – y. Оно также не зависит от выбора единиц
измерения этих переменных. Интерпретация величин коэффициентов
корреляции и оценка из значимости хорошо известны лишь для нормально
распределенных переменных. То есть, прежде чем делать строгие выводы о
корреляционной связи переменных, необходимо убедиться в том, что можно
принять гипотезу о нормальном распределении обеих переменных.
При исследовании статистической зависимости между разрывной
нагрузки (х) и устойчивостью к многократному изгибу (у) для материалов
средств защиты рук, получены промежуточные результаты вычисления
коэффициента корреляции (табл. 3.1).
Таблица 3.1 – Промежуточные результаты вычисления коэффициента корреляции
№
1
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
2
0,6
1
1,1
1,25
1,2
1,15
1,3
0,95
1,08
1,18
1,21
1,19
1,25
1,23
1,18
3
58
78
80
83,9
78
81
82,4
71
74
76,9
75,5
85
82
79
90
4
-0,56
-0,16
-0,06
0,09
0,04
-0,01
0,14
-0,21
-0,08
0,02
0,05
0,03
0,09
0,07
0,02
5
-19,72
0,28
2,28
6,18
0,28
3,28
4,68
-6,72
-3,72
-0,82
-2,22
7,28
4,28
1,28
12,28
6
0,31
0,02
0,00
0,01
0,00
0,00
0,02
0,04
0,01
0,00
0,00
0,00
0,01
0,01
0,00
7
388,75
0,08
5,21
38,23
0,08
10,78
21,93
45,11
13,81
0,67
4,91
53,05
18,35
1,65
150,88
8
11,00
-0,04
-0,13
0,57
0,01
-0,03
0,67
1,39
0,29
-0,02
-0,12
0,24
0,40
0,09
0,27
25
Окончание таблицы 3.1
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
1,45
1,18
1,35
1,26
1,2
1,16
1,24
1,1
0,96
1,21
0,98
1,19
1,22
1,1
1,26
86
75,5
80
79,6
71,4
77
89
79,3
69,2
71
70,9
72,6
85,3
74
76
0,29
0,02
0,19
0,10
0,04
0,00
0,08
-0,06
-0,20
0,05
-0,18
0,03
0,06
-0,06
0,10
8,28
-2,22
2,28
1,88
-6,32
-0,72
11,28
1,58
-8,52
-6,72
-6,82
-5,12
7,58
-3,72
-1,72
0,09
0,00
0,04
0,01
0,00
0,00
0,01
0,00
0,04
0,00
0,03
0,00
0,00
0,00
0,01
68,61
4,91
5,21
3,55
39,90
0,51
127,31
2,51
72,53
45,11
46,47
26,18
57,51
13,81
2,95
2,42
-0,05
0,44
0,19
-0,27
0,00
0,93
-0,09
1,68
-0,35
1,21
-0,17
0,47
0,21
-0,18
Для расчета коэффициента корреляции нужно рассчитать выборочные
средние для х и у по формулам (11), (12) соответственно:
Н;
циклы.
Выборочные дисперсии для х и у рассчитываются по формулам (8), (9):
Н;
циклы;
Нциклы.
Коэффициент корреляции вычисляется по формуле (7):
.
Так как коэффициент корреляции
>0,70, то зависимость между
исследуемыми величинами сильная или тесная. Следовательно, существует
некоторая положительная корреляция между разрывной нагрузкой и
устойчивостью к многократному изгибу. Насколько сильна связь между
26
исследуемыми переменными можно определить, используя метод проверки
статистических гипотез.
3.4 Проверка гипотезы о значимости коэффициента корреляции
Пусть двумерная генеральная совокупность (х, у) распределена
нормально. Из этой совокупности извлечена выборка объема n и по ней
найден
выборочный
коэффициент
корреляции
, который
оказался
отличным от нуля. Так как выборка отобрана случайно, то еще нельзя
заключить, что коэффициент корреляции генеральной совокупности также
отличен от нуля. Возникает необходимость при заданном уровне значимости
проверить нулевую гипотезу о равенстве нулю генерального коэффициента
корреляции.
Если нулевая гипотеза отвергается, то это означает, что выборочный
коэффициент корреляции значимо отличается от нуля (кратко говоря,
значим), а х и у коррелированы, то есть связаны линейной зависимостью.
Если нулевая гипотеза будет принята, то выборочный коэффициент
корреляции незначим, а х и у некоррелированы.
В качестве критерия проверки нулевой гипотезы принимают случайную
величину:
.
Величина
(13)
при справедливости нулевой гипотезы имеет распределение
Стьюдента с k=n-2 степенями свободы.
Обозначим значение критерия, вычисленное по данным наблюдений,
через
(по формуле (13)) и сформулируем правило проверки нулевой
гипотезы:
.
Для уровня значимости
и числа степеней свободы k=28
критическое значение критерия Стьюдента . Так как
, то нулевая
27
гипотеза отвергается, принимается альтернативная. Таким образом, можно
утверждать, что выборочный коэффициент корреляции получен по данным,
отобранным
из
генеральной
совокупности,
в
которой
исследуемые
переменные связаны между собой.
Доверительный интервал определяется по формуле:
.
(14)
Доверительный интервал определяется:
.
Распределение
уже при небольших n является приближенным
нормальным с математическим ожиданием:
.
(15)
Математическое ожидание определяется:
.
Доверительный интервал для :
.
По функции Лапласа
(16)
=1,96 доверительный интервал согласно
формуле (16) определяется:
.
Доверительный интервал для генерального коэффициента корреляции
равен:
.
Таким образом, существует некоторая положительная корреляция между
разрывной нагрузкой и устойчивостью к многократному изгибу.
28
4 ОБРАБОТКА РЕЗУЛЬТАТОВ МНОГОФАКТОРНОГО ЭКСПЕРИМЕНТА
4.1 Кодирование переменных
Для исследования влияния некоторых факторов на устойчивость к
многократному изгибу материалов для средств защиты рук были поставлены
эксперименты по плану ПФЭ 23. В качестве факторов, влияющих на
устойчивость к многократному изгибу y (циклы), были выбраны следующие:
1. z1 – толщина, мм:
2. z2 – разрывная нагрузка, Н:
3. z3 – удлинение при разрыве, %:
Требуется построить уравнение регрессии, учитывая все взаимодействия
факторов, проверить полученную модель на адекватность и произвести ее
интерпретацию.
План проведения экспериментов записывается в виде матрицы
планирования, в которой в определенном порядке перечисляются различные
комбинации факторов на двух уровнях. Знак «+» говорит о том, что во время
опыта значение фактора устанавливают на верхнем уровне, а знак «–»
показывает, что значение фактора устанавливают на нижнем уровне [7].
При проведении экспериментов получают значения исследуемой
величины y (табл.4.1) для каждого опыта (или серии опытов). Затем
переходят к построению математической модели.
Таблица 4.1 - Исходная матрица планирования ПФЭ 23
№
эксперимента
1
2
3
4
5
6
7
8
z1
+
+
+
+
-
Изучаемые факторы
z2
+
+
+
+
-
z3
+
+
+
+
-
y1
73
77
83,5
88
65,5
85,2
80
90,3
Результаты опытов
y2
y3
74
77,6
83,8
87,6
66
84,2
80,2
89,3
75
78,5
82
87,5
66,5
84
81,2
89,2
Под моделью понимается вид функции y = f(z1, z2,…, zk), которая
связывает изучаемый параметр со значениями факторов, лежащих в
29
интервале между верхним и нижним уровнями. Эта функция называется
уравнением регрессии..
Для обработки результатов проведенных экспериментов и дальнейшего
определения коэффициентов уравнения регрессии факторы приводят к
одному масштабу. Это достигается путем кодирования переменных.
Нижний уровень фактора zi обозначается через
через
, а верхний уровень –
.Тогда новые кодированные переменные xi будут определяться через
zi по формуле:
,
где
(17)
- центр плана;
– интервал варьирования.
Для каждого фактора находится центр, интервал варьирования и
зависимость кодированной переменной
от натуральной .
Центр плана находится по формуле:
.
Центр плана для
и
,
(18)
:
мм;
Н;
%.
Интервал варьирования находится по формуле:
.
Интервал варьирования для
мм;
Н;
,
и
:
(19)
30
%.
Данные расчетов представлены в табл. 4.2.
Таблица 4.2 – Кодирование факторов
Факторы
Верхний
уровень
Нижний
уровень
Центр
Интервал
варьирования
0,9
0,6
0,75
0,15
1,47
0,78
1,125
0,345
1,8
1,5
1,65
0,15
Зависимость кодированной
переменной от натуральной
При таком кодировании все новые переменные будут принимать
значения от –1 до +1.
4.2 Построение матрицы планирования в кодированных переменных
Линейное уравнение регрессии относительно новых переменных имеет
вид:
(20)
Для изучения влияние парных взаимодействий различных факторов на
исследуемый параметр, то уравнение регрессии записывают в виде
или
.
(21)
Если надо учесть другие взаимодействия, то число слагаемых
увеличивают.
Прежде, чем определять коэффициенты выбранной модели, матрицу
планирования записывают относительно новых переменных. Далее матрицу
дополняют (если это требует вид выбранного уравнения регрессии)
31
столбцами знаков «+» и «–», соответствующих уровням, на которых будут
находиться взаимодействия факторов (см. столбцы x1x2, x1x3, x2x3 в табл. 4.3).
Знаки этих столбцов получают с помощью исходной матрицы планирования.
Обозначим знак «+» или «–» в матрице планирования за xji, который
соответствует j-ому опыту (j=1,...,n) для i-го фактора (i=1,…,k). При этом
знак «+» показывает, что кодированная переменная принимает значение +1, а
знак «–» соответствует значению –1. Тогда знаки (или уровни варьирования)
для взаимодействия факторов xr и xp вычисляются простым перемножением:
xjr ∙ xjp , j=1,…,n [8].
В таблице 4.3 знаки для взаимодействия x1x2 получены таким образом:
для 1-го опыта (j=1)
x11x12= (+1)(+1) =+1,
для 2-го опыта (j=2)
x21x22= (-1)(+1) = -1,
для 3-го опыта (j=3)
x31x32=(+1)(-1) = -1, и т.д.
Обычно проводят несколько серий опытов для каждого эксперимента.
Это необходимо для проверки уравнения на адекватность.
Адекватность – это способность модели предсказывать результаты
эксперимента в некоторой области с требуемой точностью. Результаты
опытов в каждом j-ом эксперименте (j=1,…, n) записывают в правые столбцы
матрицы
планирования.
В
последнем
столбце
записывают
средние
выборочные значения полученных результатов для каждой серии опытов
(табл. 4.3). Если обозначить за yji значение результата, полученного в i-ом
опыте (i=1,…,m) для j-ого эксперимента (j=1,…,n), то выборочное среднее
для каждого эксперимента вычисляют по известной формуле:
(22)
Для каждого эксперимента по формуле (22) были рассчитаны средние
выборочные результатов:
циклы;
циклы;
циклы;
32
циклы;
циклы;
циклы;
циклы;
циклы.
В табл. 4.3 представлена матрица планирования с учетом всех
взаимодействий и средних значений отклика.
Таблица 4.3 – Матрица планирования для обработки результатов
№
эксперимента
1
2
3
4
5
6
7
8
Факторы
+
+
+
+
-
+
+
+
+
-
+
+
+
+
-
Взаимодействие
+
+
+
+
+
+
+
+
+
+
+
+
Результаты опытов
+
+
+
+
-
73
77
83,5
88
65,5
85,2
80
90,3
74
77,6
83,8
87,6
66
84,2
80,2
89,3
75
78,5
82
87,5
66,5
84
81,2
89,2
Среднее
результатов
74
77,7
83,1
87,7
66
84,46
80,46
89,6
4.3 Нахождение коэффициентов уравнения регрессии
Коэффициенты уравнения регрессии находят с помощью метода
наименьших квадратов.
Так как матрица планирования ПФЭ 2k должна удовлетворять
определенным требованиям (такие матрицы с заданными требованиями уже
построены),
то
формулы,
определяющие
коэффициенты
уравнения
регрессии, достаточно просты:
Расчет коэффициентов уравнения регрессии проводится по формулам:
(23)
,
(24)
33
(25)
и т.д., если учитываются другие взаимодействия.
Коэффициенты уравнения регрессии рассчитываются по формуле (23):
;
;
;
;
;
;
;
4.4 Проверка коэффициентов на значимость
Полученные коэффициенты необходимо проверить на значимость. Это
можно сделать с помощью критерия Стьюдента: если , то
значим , то
незначим и его полагают равным нулю в уравнении регрессии.
Критическую точку tкр. находят из таблиц распределения Стьюдента по
числу степеней свободы n(m - 1) и с заданным уровнем значимости б для
случая
двусторонней
критической
области.
Среднее
квадратическое
отклонение коэффициентов Sкоэф. зависит от дисперсии воспроизводимости
результатов по всем проведенным опытам
и вычисляется по формуле [9]:
.
Дисперсия
воспроизводимости
характеризует
(26)
ошибку
всего
эксперимента. В случае равномерного дублирования опытов (т.е. при
34
одинаковом числе наблюдений в каждом эксперименте) для расчета
используют формулу:
(27)
где n – число экспериментов (число строк в матрице ПФЭ);
m – число опытов (наблюдений) в каждом эксперименте;
yji – результат отдельного i-го наблюдения в j-ом эксперименте;
– среднее выборочное значение наблюдений для j-ого эксперимента.
Для облегчения расчетов нахождения дисперсии воспроизводимости
лучше использовать формулу (27) в ином виде:
(28)
В формуле (28) данные внутренние суммы
являются выборочными
дисперсиями результатов опытов для j-го эксперимента (j=1,…, n).
Для удобства расчеты заносятся в табл. 4.4.
Таблица 4.4 – Расчет выборочных дисперсий
№
Опыта
1
2
3
4
5
6
7
8
y1
y2
y3
у
73
77
83,5
88
65,5
85,2
80
90,3
74
77,6
83,8
87,6
66
84,2
80,2
89,3
75
78,5
82
87,5
66,5
84
81,2
89,2
74
77,7
83,1
87,7
66
84,46
80,46
89,6
у
у
1
0,49
0,16
0,09
0,25
0,53
0,21
0,49
у
у
0
0,01
0,49
0,01
0
0,07
0,07
0,09
у
у
1
0,64
1,21
0,04
0,25
0,21
0,53
0,16
S2
1
0,57
0,93
0,07
0,25
0,41
0,41
0,37
Далее находится сумма полученных значений дисперсий:
циклы.
Отсюда дисперсия воспроизводимости:
циклы.
Среднее квадратическое отклонение коэффициентов определяется по
формуле (26):
.
35
Из таблиц распределения Стьюдента по числу степеней свободы
n(m-1)=8·2=16 при уровне значимости α=0,05 находится . Следовательно,
Данное значение сравнивается с коэффициентами уравнения регрессии,
представленными в табл. 4.5.
Сравнивая полученное значение видно, что все коэффициенты кроме
,
больше
по
коэффициенты кроме
абсолютной
,
величине 0,3. Следовательно, все
значимы. Полагая, что
,
= 0, то
уравнение регрессии в кодированных переменных:
. (29)
4.5 Проверка уравнения на адекватность
Проверка на адекватность полученного уравнения регрессии со
значимыми коэффициентами осуществляется с помощью критерия Фишера:
если Fрасч < Fтабл , то уравнение адекватно, в противном случае – неадекватно.
Расчетное значение критерия Fрасч. определяется по формуле:
,
где
(30)
– дисперсия воспроизводимости;
– остаточная дисперсия (или дисперсия адекватности).
Остаточная дисперсия
вычисляется следующим образом:
,
где
(31)
– число экспериментов;
– число опытов в каждом эксперименте;
– число значимых коэффициентов в уравнении регрессии;
– значение изучаемого параметра, вычисленное по уравнению
регрессии со значимыми коэффициентами для j-ого эксперимента;
– среднее выборочное значение наблюдений для j-ого эксперимента.
Табличное значение критерия Fтабл. находят из таблиц критических точек
распределения Фишера по заданному уровню значимости
б и по
36
соответствующим степеням свободы k1= n - r и k2= n(m - 1). Степень
свободы k1 соответствует степени свободы числителя формулы (31) –
остаточной дисперсии
, а k2 – степень свободы знаменателя формулы
(22) – дисперсии воспроизводимости
.
Проверка полученного уравнения (29) на адекватность происходит по
критерию Фишера. Так как дисперсия воспроизводимости найдена в
предыдущем пункте, то для определения расчетного значения критерия Fрасч.
необходимо вычислить остаточную дисперсию
. Для этого нужно найти
значения изучаемого параметра по полученному уравнению регрессии ỹj
(j=1,…, 8), подставляя +1 или -1 вместо xi в соответствии с номером
j
эксперимента из табл. 4.1:
циклы;
циклы;
циклы;
циклы;
циклы;
циклы;
циклы;
циклы.
Далее рассчитывается остаточная дисперсия по формуле (31):
+
2+
2+
2+
2++
2
=1,56 циклы.
Расчетное значение критерия Фишера определяется по формуле (30):
.
По таблице критических точек распределения Фишера находится
табличное значение. Уровень значимости
1)=8*2=16. Таким образом, Fтабл=4,49.
k1=n-r=8-7=1; k2=n(m-
37
Так как, то уравнение регрессии (29) адекватно.
4.6 Интерпретация полученной модели
По уравнению (29) видно, что наиболее сильное влияние на
устойчивость к многократному изгибу оказывают фактор
нагрузка и
– разрывная
– удлинение при разрыве, так как они имеет наибольший по
абсолютной величине коэффициент. После этих факторов по силе влияния на
отклик
(устойчивость
к
взаимодействие факторов
многократному
изгибу)
идут:
двойное
- сочетание толщины и разрывной нагрузки.
Так как коэффициенты при отрицательны, то это означает, что с
уменьшением всех факторов и взаимодействий факторов значение отклика
будет возрастать, а с увеличением – убывать.
Запись
уравнение
регрессии
происходит таким образом:
(29)
в
натуральных
выражается через
переменных,
, которое
берется из последнего столбца табл. 4.3. Уравнение регрессии в натуральных
переменных:
.
Окончательный вид уравнения регрессии
натуральные переменные:
.
после преобразования в
38
ЗАКЛЮЧЕНИЕ
В
данной
курсовой
работе
проводилось
исследование
влияния
некоторых факторов на устойчивость к многократному изгибу материалов
для средств защиты рук, рекомендуемых для работ на открытом воздухе в
условиях Севера и Крайнего Севера, были поставлены эксперименты по
плану ПФЭ 23. В качестве факторов, влияющих на устойчивость к
многократному изгибу, были выбраны следующие: толщина; разрывная
нагрузка; удлинение при разрыве. Факторы определялись с помощью
номенклатуры показателей качества. Перед постановкой эксперимента по плану
ПФЭ, определялась относительная погрешность результатов измерений и
проводился
корреляционный
анализ
экспериментальных
данных.
Исследовалась статистическая зависимость между разрывной нагрузки и
устойчивостью к многократному изгибу для материалов средств защиты рук.
Так как коэффициент корреляции
и
устойчивостью
к
, то между разрывной нагрузкой
многократному
изгибу
существует
некоторая
положительная корреляция. Метод проверки статистических гипотез доказал
этот факт.
При
проведении
полного
факторного
эксперимента
полученное
уравнение регрессии показывает, что наиболее сильное влияние на
устойчивость к многократному изгибу средств защиты рук при низких
температурах, оказывают факторы разрывная нагрузка и удлинение при
разрыве, так как они имеет наибольший по абсолютной величине
коэффициент. После этих факторов по силе влияния на отклик (устойчивость
к многократному изгибу) идут: двойное взаимодействие факторов сочетание толщины и разрывной нагрузки. Так как все коэффициенты
отрицательны, то с уменьшением всех факторов и взаимодействий факторов
значение отклика будет возрастать, а с увеличением – убывать.
Полный факторный эксперимент относится к числу планов, которые
являются наиболее эффективными при построении линейных моделей.
39
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. ГОСТ 4.493-89 Система показателей качества продукции (СПКП).
Материалы для средств защиты рук. Номенклатура показателей. – Введ.
1989-28-12. – М. : Изд-во стандартов, 1990. – 17 с.
2. ГОСТ 8978-75 Кожа искусственная и пленочные материалы. Методы
определения устойчивости к многократному изгибу (с Изменениями N 1, 2)
– Введ. 1977-07-01. – М. : Изд-во стандартов, 1998. – 17 с.
3. Жуков В. К. Математическая статистика с элементами теории
планирования эксперимента: учеб. пособие: Л.В. Горская, В.Н. Пиунова, В.С.
Смирнова. Саратовский полит. институт, 1975, 103 с
4. Планирование
эксперимента
в
исследовании
технологических
процессов / Под ред. Э.К.Лецкого. - М.: Мир, 1977.
5. Бекряев В. И. Основы теории эксперимента. Учебное пособие. - СПб.:
Изд. РГГМУ, 2001 - 266 с.
6. Васильков Ю.В. Статистические методы в управлении предприятием
:учебное пособие / Ю.В. Васильков, Н. Иняц. – М. : Стандарты и качество ,
2008. – 279 с.
7. Реброва, Ирина Анатолиевна. Теория планирования эксперимента
[Электронный ресурс] : учебное пособие / И.А. Реброва. – Электрон. дан. −
Омск:
СибАДИ,
2016.
–
Режим
доступа:
http://bek.sibadi.org/fulltext/esd104.pdf, свободный после авторизации.
8. Соколовская
И.Ю.
Полный
факторный
эксперимент
/
И.Ю.
Соколовская // Методические указания для самостоятельной работы
студентов. – Новосибирск: НГАВТ, 2010. – 36 с.
9. Бирюков С.В., Чередов А.И. Метрология: Тексты лекций. – Омск:
Изд-во ОмГТУ, 2000, - 110 с.
10.
Сайт
«Ekonomstat».
–
Режим
доступа:
https://ekonomstat.ru/polnye-otvety-na-voprosy-k-ekzamenu-statistika/104ponyatie-korrelyacionnoj-svyazi.html (корреляция)
40
11.
Средства
индивидуальной
защиты
персонала
предприятий
атомной промышленности и энергетики: каталог-справочник. – М., 2015 –
257 с.
Отзывы:
Авторизуйтесь, чтобы оставить отзыв