Слайд 1Дисперсионный анализ
Однофакторный дисперсионный анализ
Слайд 2Дисперсионный анализ
Данный вид анализа применяют в тех случаях, когда необходимо
сопоставить не 2, а большее число результатов однотипных экспериментов.
Смысл дисперсионного
анализа заключается в следующем – из общей суммы квадратов дисперсии вычитают сумму квадратов отклонений по изучаемым факторам (межфакторная дисперсия). В результате чего получают остаточную сумму квадратов дисперсии, которая характеризует влияние различных факторов.
Слайд 3Таблица исходных данных
Для n разных уровней некоторого фактора проводят по
m измерений (для каждого уровня) величины y. Затем проверяют гипотезу
о том, что влияние фактора на средние значения для каждого уровня существенно.
Слайд 4Расчетные формулы
Вычисляют общее среднее и среднее для данного уровня фактора
.
Вычисляют остаточную дисперсию:
Вычисляют межфакторную дисперсию:
Слайд 5Дисперсия генеральной совокупности:
Определяют статистику:
Сравнивают ее с табличным значением
F(1,2), ((1=n-1), (2=nm-n)).
F*>F(1,2), фактор статистически значим
Слайд 6Дисперсионный анализ
Двухфакторный анализ
Слайд 7Исходные данные
Для двух факторного анализа необходимо задать m, как число
измерений величины у.
Пусть второй фактор k – принимает значения
от 1 до р, а первый фактор – от 1 до n.
Запишем у с тремя индексами k, i, j, где j – число повторений измерения ( ykij ).
Слайд 8Проверяются три гипотезы:
Влияние первого фактора – статистически значимое.
Влияние второго фактора
– статистически значимое.
Взаимодействие между факторами – статистически значимое.
Слайд 9Вычисления ведут по следующей схеме:
Определяют общее и частное среднее.
Дисперсия
генеральной совокупности
Межфакторные S1 и S2:
Слайд 10Дисперсия взаимодействия:
Остаточная дисперсия:
Проверка статистической значимости:
Слайд 11Пример
необходимо проанализировать зависимость пластичности полуфабриката от размера субзерна ( 1
фактор ) и объемной доли избыточных фаз ( 2 фактор
). В каждой точке испытывали по три образца
Слайд 12Расчеты
1,41; 0,081
; 632,9 ;
0,231 ; 0,095.
Поделив друг на друга, сравниваем с табличными значениями:
Для 1-го фактора: F0.05(3/40)=2,84
Для 2-го фактора: F0,05(4/40)=2,61
Для взаимодействия: F0,05(12/40)=2,00
После расчета выясняем, что размер субзерна не влияет на пластичность полуфабриката, на него влияет лишь доля избыточных фаз.
Слайд 13Дисперсионный анализ
Латинские квадраты
Слайд 14Латинский квадрат
– это квадратная таблица размером n х n элементов,
расположенных на поле квадрата таким образом, что каждый из них
встречается в каждом столбце и в каждой строке только по одному разу.
Если строкам, столбцам и элементам выписной таблицы поставить в соответствие уровни каких – либо факторов ( A, B, C ), то латинский квадрат можно рассматривать как план из эксперимента, позволяющий провести дисперсионный анализ с тремя факторами.
Слайд 15Оценка статистической значимости
– эффект взаимодействия.
ПФЭ=43=
Слайд 16Латинский квадрат 4х4
При проведении эксперимента порядок реализации опытов необходимо рандомизировать,
статистические свойства оценок при этом улучшаются.
Слайд 17Ортогональные планы
Латинские квадраты называются взаимно ортогональными, если при размещении их
элементов на поле общей таблицы каждая пара элементов двух квадратов
встречается только по одному разу.
Слайд 18Матрица плана эксперимента
На базе латинского квадрата
3х3
Уровни факторов С и D располагаются по полю таблицы в
виде ортогонального плана
Слайд 19Методика расчета
Вычисление сумм результатов по строкам, столбцам и одноименным
буквам: ∑ A,B,C,D ( по отдельности).
Вычисление вспомогательных расчетных сумм:
где
р – текущий индекс ячейки квадратов
Слайд 20Средние суммы квадратов по строкам, столбцам и латинским буквам:
;
;
;
Вычисление корректирующего члена
Вычисление сумм квадратов, характеризующих эффекты
строк, столбцов, латинских букв
.
Слайд 21Вычисление остаточной суммы квадратов
Вычисление оценок эффектов строк, столбцов и
латинских букв матрицы плана, определяемых как частные от деления соответствующих
сумм квадратов на числа степеней свободы, с которой они оцениваются
SSA/(n-1); SSB/(n-1); SSC/(n-1); SSD/(n-1).
Проверить по критерию Фишера значимость оценок эффектов строк, столбцов и буквенных оценок
Слайд 22Дисперсионный анализ результатов экспериментов
(без повторных опытов)
Слайд 23Дисперсионный анализ результатов экспериментов
(с повторными опытами)
Слайд 24Пример
Исследовали малоцикловую усталость (МЦУ) стали ВНЛ-3 в зависимости от чистоты
обработки поверхности.
Предварительные исследования показали, что, хотя с повышением чистоты
поверхности средние значения долговечности образцов несколько увеличиваются, для возрастающих классов чистоты поверхности значения МЦУ , как правило, попадают в пределы доверительных интервалов оценки средних значений для более низких классов чистоты.
Было отмечено также, что результаты испытания образцов из металла разных плавок, образцов, испытанных в разное время и на разных испытательных машинах, ложатся на кривые малоцикловой усталости с большим разбросом. Так, при max= 800МПа разброс значений долговечности отдельных образцов составлял (суммарно по всем имевшимся результатам испытаний):
Для МЦУ на базе тысячи циклов
(19-59) циклов для точения, 4;
(21-55) циклов для точения, 5;
(40-60) циклов для точения, 6;
(29-65) циклов для шлифования, 7.
Слайд 25Полученные результаты не могли быть объяснены только случайным рассеянием долговечности.
В этой связи была поставлена задача произвести дисперсионный анализ результатов
испытаний образцов, имеющих разную чистоту поверхности ( Х| ). В качестве характерных источников неоднородности условий испытаний были выбраны:
уровень напряжений цикла Х2;
плавка металла, Х3;
испытательная машина , Х4.
Х1 – чистота поверхности образца ( 4 5 6 7 ).
Х2 – уровень напряжения цикла ( нагрузка ) 100, 90, 80, 70 кгс/мм2
Х3 – плавка ( A, B, C, D ).
Х4 – разные испытания машины одного класса ( α, β, γ, δ).
Слайд 28Средние значения результатов экспериментов, тыс. циклов
В каждой ячейке плана было
испытано по 4 образца. Испытания проводились при рандомизации экспериментов.
Слайд 30Предварительные расчеты дисперсионного анализа
Слайд 31Вспомогательные расчеты
Суммы квадратов, характеризующие проверяемые эффекты, определены через вспомогательные суммы:
Слайд 32Дисперсионный анализ результатов
Средний квадрат
Слайд 33Выводы
Эффект от испытательных машин статистически не значим.
Шероховатость поверхности может
быть использована 4 или 5, т.к. эта сталь имеет невысокую
чувствительность к чистоте поверхности.
Уровень нагрузки статистически значим.
Эффект от различия в плавках может перекрывать эффект от влияния чистоты поверхности.