Тема: Изучение свойств случайных величин, планирование эксперимента и анализ данных

  • Вид работы:
    Курсовая работа (т)
  • Предмет:
    Математика
  • Язык:
    Русский
  • Формат файла:
    MS Word
  • Размер файла:
    496,44 Кб
Изучение свойств случайных величин, планирование эксперимента и анализ данных
Изучение свойств случайных величин, планирование эксперимента и анализ данных
Вы можете узнать стоимость помощи в написании студенческой работы.
Помощь в написании работы, которую точно примут!









Изучение свойств случайных величин, планирование эксперимента и анализ данных

Содержание

Введение

. Одномерные случайные величины

.1 Формирование выборки объёмом n=15

1.1.1 Вычисление среднего и дисперсии

1.1.2 Оценка нормальности выборки объёмом 15

.1.3 Определение доверительного интервала для математического ожидания

.1.4 Определение доверительного интервала для дисперсии

.2 Получение второй выборки объёмом 100

.2.1 Оценка нормальности выборки объёмом 100

.2.2 Вычисление среднего и дисперсии

.2.3 Проверка гипотезы о равенстве математических ожиданий этих двух выборок при условии равенства их генеральных дисперсий

.2.4 Оценка доверительного интервала для среднего первой выборки, используя данные второй выборки

. Двумерные случайные величины

.1 Выбор двух функций и построение корреляционного поля

.2 Изучение зависимости выбранного У от одного из факторов Х

.2.1 Вычисление условных средних У для фиксированных значений Х

.2.2 Вычисление условных дисперсий У для фиксированных значений Х

.3 Построение линии регрессии У по Х (эмпирической и приближённой)

. Дисперсионный анализ и планирование эксперимента

.1 Выбор факторов Х и функций отклика показателей качества У1 и У2, краткое описание эксперимента

.2 Составление плана эксперимента

.3 Составление матрицы эксперимента

.4 Дисперсионный анализ греко-латинского куба второго порядка

.5 Проверка условий применимости дисперсионного анализа, критерий Дункана для показателей качества Y1 и Y2

4. Регрессионный анализ

Заключение

Список литературы

математический ожидание дисперсия регрессия

Введение

Целью курсовой работы является изучение показателей качества (ПК), как случайных величин, и доказательство факта влияния на них нескольких факторов, действующих одновременно. По имитационной модели процесса необходимо получить значения двух функций отклика (ПК), выбрав несколько факторов и задавая им градации. Модель является таблицей EXCEL.

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

Одномерные случайные величины

.1 Формирование выборки объемом n=15

Используя модель переменных, выбираем функцию отклика Y3 и формируем выборку объемом 15. Выборка представлена в таблице 1.

Таблица 1 - Выборка объемом n=15

123,2130,0131,0132,8133,0135,8138,6141,0143,4143,4143,6148,4150,0150,4157,2

.1.1 Вычисление среднего и дисперсии

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

, (1)

где n - объем выборки;

yi - наблюдаемые значения выборки.

Определяем дисперсию:

, (2)

Для нашей выборки имеем:

Проверка наличия грубых погрешностей

Под грубой погрешностью измерения понимается погрешность, существенно превышающая ожидаемую при данных условиях. Она может быть сделана вследствие неправильного применения прибора, неверной записи показаний прибора, ошибочно прочитанного отсчета и т.п.

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

критерий "трех сигм" (надежен при числе измерений n>20);

критерий Романовского (применяется, если число измерений n<20);

критерий Шарлье (используется, если число наблюдений в ряду великоn>20);

вариационный критерий Диксона (мощный критерий с малыми вероятностями ошибок).

Для полученной выборки объема n=15 воспользуемся критерием Романовского:

- вычисляем отношение для каждого значения из выборки по формуле:


и сравниваем с табличным критерием.

по таблице 7.1 [3] на уровне значимости 0,05 для n=15 находим табличный критерий . Если окажется больше, то этот результат следует отбросить.


По результатам расчета, используя данные таблицы 1 делаем вывод о том что грубых погрешностей нет.

1.1.2 Оценка нормальности

Одним из способов проверки нормальности распределения является вычисление особых параметров выборочной совокупности результатов анализа, носящих название асимметрии А и эксцесса Е.

Воспользуемся описательной статистикой для нахождения значений асимметрии и эксцесса:

Столбец1Среднее140,127Стандартная ошибка2,383Медиана141,000Мода143,500Стандартное отклонение9,229Дисперсия выборки85,168Эксцесс-0,509Асимметричность0,023Интервал34,000Минимум123,200Максимум157,200Сумма2101,900Счет15,000Уровень надежности(95,0%)5,111

По величине асимметрии и эксцесса можно косвенно судить о нормальности распределения. Более достоверной является оценка с использованием дисперсий этих величин, которые являются функциями от кратности анализа:

(4)

(5)

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

и ,

то наблюдаемое распределение можно считать нормальным.

В нашем случае: и ,

и

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

Для оценки нормальности распределения можем воспользоваться не только критерием асимметрии и эксцесса, но и такими критериями, как:

критерий Пирсона (;

критерий Колмогорова-Смирнова;

критерий Шапиро-Вилка;

критерий Жака-Берра;

критерий Андерсена-Дарлинга;

критерий Дэвида-Хартли-Пирсона;

критерий Саркади;.

критерий Оя и др..

Воспользуемся критерием Дэвида-Хартли-Пирсона.

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

Статистика критерия имеет вид:

(6)

где R- размах, вычисляемый по формуле : R=Ymax-Ymin;

S - стандартное отклонение.

Гипотеза нормальности принимается, если .

Зададимся уровнем значимости , тогда

R=157,2-123,2=34,00

По таблице 75 [4] для n=15 и находим ,

,97<3.68<4.17,

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

1.1.3 Определение доверительных интервалов для математического ожидания

Доверительные интервалы для математического ожидания находим, используя критерий Стьюдента.

Рассмотрим случайную величину , которая согласно следствию из теоремы о распределении выборочных характеристик распределена по закону Стьюдента . При заданном значении , пользуясь таблицей П2 [1], вычислим значение из условия:

, (7)

где - надежность интервальной оценки.

α - генеральное среднее.

Из условия (7) получаем:

(8)

Таким образом, интервальная оценка надежности для неизвестной генеральной средней а имеет границы:

(9)

Выразим границы интервала через исправленную дисперсию . Так как =, то . Поэтому

. (10)

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

, (11)

По выборке объема 15 нормально распределенной найдено среднее значение 140,12. Построим доверительный интервал для математического ожидания с надежностью γ =0,95.

Пользуясь таблицей П2 [1] находим величину t(0,95;15)=2,15.

Тогда доверительные границы для математического ожидания с доверительной вероятностью 0,95:

;

Окончательно с надежностью 0,95 получаем, что параметр а заключен в интервале: .

1.1.4 Определение доверительных интервалов для дисперсии

По выборке объёма 15, имеющей нормальное распределение, найдено значение S=9.23. Найдем доверительный интервал для σ с надежностью γ=0,95

Доверительный интервал покрывающий с заданной надежностью находим по формуле:


По таблице приложения 4 [2] по данным γ=0,95 и n=15 находим q=0,46, подставляя значения в (12) получаем: 9,23(1-0,46)< 9.23(1+0.46)

Искомый доверительный интервал: 4,98<13.47

.2 Формирование выборки объемом более 60

.2.1 Вычисление среднего и дисперсии

С помощью модели эксперимента формируем вторую выборку (Y3) объемом 100

Таблица 2 - Значения функции Y3

№ опытаY3№ опытаY3№ опытаY3№ опытаY31113,8026131,8051140,2076148,202118,2027131,8052140,6077149,003119,8028133,0053141,278149,204120,8029133,2054141,3079149,205120,8030133,3055142,2080149,26121,2031133,6056142,4081149,407121,8032133,8057142,6082149,408122,8033134,2058142,8083150,009122,8034134,2059143,2084150,2010122,8035135,0060143,8085150,4011124,6036135,0061144,4086150,6012124,6037135,6062144,4087150,8013124,8038136,0063144,8088152,9014126,6039136,0064145,0089153,0015126,8040136,4065145,4090153,2016127,2041136,6066145,5091154,2017127,6042137,8067145,6092155,2018128,2043138,4068145,7093155,2019128,8044138,6069146,1094155,3020129,0045139,0070146,3095155,5021129,2046139,0071146,3096156,2022129,6047139,2072146,3097156,8023130,4048139,4073146,4098159,8024130,8049139,6074147,4099163,8025131,2050139,6075147,80100164,20

Определяем среднее результатов наблюдений по формуле (1)

139,6

Дисперсию вычисляем по формуле (2):

D=122,15

S=11.05

Проверим на наличие грубых погрешностей по критерию "трех сигм": Критерий "трех сигм" применяется для результатов измерений, распределенных по нормальному закону. По этому критерию считается, что результат, возникающий с вероятностью q £ 0,003, маловероятен и его можно считать промахом, если |xi - x| < 3σ, где σ - оценка СКО измерений. Данный критерий надежен при числе измерений n ≥ 20…50.

(13)

При проверке условия (13) грубых погрешностей не обнаружено.

1.2.2 Проверка "нормальности" полученной выборки

Воспользуемся критерием Дэвида-Хартли-Пирсона.

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

Статистика критерия имеет вид:


где R- размах, вычисляемый по формуле : R=Ymax-Ymin;

S - стандартное отклонение.

Гипотеза нормальности принимается, если

Зададимся уровнем значимости , тогда R=164,2-113,8=50,4

S=11,05

По таблице 75 [4] для n=100 и находим ;

,31<4,56<5,90. так как условие выполняется следовательно гипотезу о нормальности распределения принимаем.

Также для оценки нормальности данного распределения построим гистограмму и проведем описательную статистику.

Таблица 3 - Описательная статистика для выборки n=100

ПараметрЗначениеПараметрЗначениеСреднее139,627Асимметричность-0,113Стандартная ошибка1,105Интервал50,400Медиана139,900Минимум113,800Мода122,800Максимум164,200Стандартное отклонение11,052Сумма13 962,700Дисперсия выборки122,154Счет100,000Эксцесс-0,621Уровень надежности(95,0%)2,193

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

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

и ,

то наблюдаемое распределение можно считать нормальным.

В нашем случае: и ,

и ,

Рисунок 1 - Гистограмма для выборки n=100

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

1.2.3 Проверка гипотезы о равенстве двух математических ожиданий двух выборок в предположении равенства их генеральных дисперсий

Проверим нулевую гипотезу H0: M(Y1)=M(Y2) на уровне значимости 0,05 для двух выборок n=15 и m=100

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


Критерий Z - нормированная нормальная случайная величина, так как М(Z)=0, при справедливости нулевой гипотезы (Z)=1. Критическая область строится в зависимости от вида конкурирующей гипотезы.

Первый случай: нулевая гипотеза H0: M(Y1)=M(Y2). Конкурирующая гипотеза H1: M(Y1)≠M(Y2). В этом случае строим двустороннюю критическую область, исходя из требования, чтобы вероятность попадания критерия в эту область в предположении справедливости нулевой гипотезы была равна принятому уровню значимости 0,05

По таблице функции Лапласа (приложение 2 [2]) находим критическую точку по равенству

Если - нет оснований отвергать нулевую гипотезу, если - нулевую гипотезу отвергают.

Найдем для наших выборок:

По таблице 2 [2] находим


Так как следовательно наблюдаемое значение попало в область допустимых значений. Нулевую гипотезу о равенстве математических ожиданий принимаем на уровне значимости 0,05.

Второй случай: нулевая гипотеза H0: M(Y1)=M(Y2). Конкурирующая гипотеза H1: M(Y1)>M(Y2). В этом случае строим правостороннюю критическую область, исходя из требования, чтобы вероятность попадания критерия в эту область в предположении справедливости нулевой гипотезы была равна принятому уровню значимости 0,05

Найдем :

По таблице функции Лапласа (приложение 2 [2]) находим критическую точку по равенству



Так как попадает в область принятия гипотезы, нулевую гипотезу о равенстве математических ожиданий принимаем на уровне значимости 0,05.

Третий случай: нулевая гипотеза H0: M(Y1)=M(Y2). Конкурирующая гипотеза H1: M(Y1)<M(Y2). В этом случае строим левостороннюю критическую область, исходя из требования, чтобы вероятность попадания критерия в эту область в предположении справедливости нулевой гипотезы была равна принятому уровню значимости 0,05

По таблице функции Лапласа (приложение 2 [2]) находим критическую точку по равенству

Полагая, что находим:


Так как нулевую гипотезу о равенстве математических ожиданий принимаем на уровне значимости 0,05.

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

Генеральная совокупность (данные второй выборки) распределена по нормальному закону, причем параметр известен и равен 11,05. По выборке объема n=15 найдено 140,12 определим интервальную оценку для математического ожидания с надежностью .

Доверительный интервал имеет вид :


Находим t из соотношения 2Ф(t)=0.95 получаем Ф(t)=0.475. По таблице приложения 2[2] находим t=1,96.

Подставляем полученные значения в (15) получаем искомый интервал:


C надежностью 0,95 параметр а находится в интервале: ;

2. Двумерные случайные величины

Методы корреляционного и регрессионного анализов широко применяются для выявления и описания зависимостей между случайными величинами по экспериментальным данным. Для экспериментального изучения зависимости между величинами проведем 100 опытов. Результат каждого опыта дает пару значений (Y2i , Y3i).

.1 Построение корреляционного поля для двух функций (Y2 и Y3)

О наличии или отсутствии корреляции между величинами можно судить по виду поля корреляции. Корреляция на графике может быть положительной или отрицательной. Отрицательная корреляция( - корреляция, при которой увеличение одной переменной связано с уменьшением другой переменной, при этом коэффициент корреляции может быть отрицательным; положительная корреляция (- корреляция, при которой увеличение одной переменной связано с увеличением другой переменной, при этом коэффициент корреляции может быть положительным.

Рисунок 2 - Примеры корреляционных полей

а - отсутствие корреляции; б - линейная положительная и отрицательная связи; в -корреляционная положительная связь; г - корреляционная отрицательная связь;

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

Таблица 5 - Значения Y2 и Y3 при постоянных уровнях всех действующих факторов.

№ п/пY2Y3№ п/пY2Y3№ п/пY2Y3№ п/пY2Y31139,57153,0026112,17123,6051109,37131,8076138,37151,802134,57144,0027138,77155,2052115,17140,6077105,77123,203121,37144,8028137,37147,8053127,17137,6078108,77126,204129,77142,2029121,77137,2054151,77165,2079121,57137,005120,37135,8030120,17135,6055130,77138,2080121,17133,606112,37121,8031123,57133,0056133,37142,8081105,37119,807119,77132,2032148,77157,2057121,77141,2082135,97148,408115,37133,8033146,37159,8058129,37138,8083127,77148,209137,77155,2034114,77131,2059123,17144,6084136,37155,8010142,37156,8035120,77133,2060154,77161,2085108,77126,2011127,97141,4036131,97140,4061122,17137,6086132,37149,8012106,37120,8037119,97138,4062128,57142,0087138,77160,2013116,77134,2038125,77137,2063136,57143,0088105,37125,8014137,57155,0039134,37157,8064122,77134,2089138,77158,2015124,57139,0040106,37119,8065144,37153,8090117,77135,2016121,97137,4041115,17132,6066112,37136,8091133,37149,8017114,77132,2042131,97144,4067126,17138,6092121,17143,6018132,17143,6043123,57142,0068121,77137,2093113,77136,2019130,57146,0044125,97131,4069120,57141,0094137,97146,4020116,77136,2045130,57147,0070107,17122,6095121,77134,2021127,77145,2046139,77148,2071126,37142,8096135,97153,4022127,57142,0047113,37135,8072110,77135,2097123,37131,8023125,97143,4048107,37123,807398,37120,8098127,97139,4024134,37147,8049116,97138,4074123,57141,0099136,77153,2025141,37159,8050126,17146,6075109,17126,6010118,77137,20

Рисунок 3 - Корреляционное поле

Коэффициент корреляции вычисляется по формуле:


Рассчитаем коэффициент корреляции .

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

2.2 Изучение зависимости Y2 от одного из факторов -от X4

Значения Y2 сводим в таблицу, изменяя значение фактора X4=0, X4=1, X4=2, X4=3, X4=4.

Таблица 6 - Зависимость Y2 от X4

№ опытаX4=0X4=1X4=2X4=3X4=4

Продолжение таблицы 6

№ опытаX4=0X4=1X4=2X4=3X4=4150,3759,4790,37101,37114,57237,3759,6783,17113,87111,37344,3774,9489,97109,87124,17449,3767,8777,7790,37117,37540,3763,6794,17105,87139,57631,3752,3775,7795,67123,97749,3765,8777,57104,97131,97861,3767,7787,9787,67125,57953,3763,4799,57106,97141,771037,3782,9763,37134,17119,171128,3760,5768,3798,97125,971258,3756,3795,17108,87140,371346,3752,4765,3791,67116,171438,3772,9790,37103,87136,971549,3768,6786,7794,97123,371655,3764,6792,3792,27123,171725,3771,7791,5796,97114,571850,3756,5771,57112,17129,171938,3762,5779,7792,67129,372041,3765,6791,17117,87132,972155,3775,6777,77106,57125,372234,3774,9792,97114,87117,572352,3755,37100,57103,37142,372433,3774,8762,37114,57106,372550,3762,7790,9798,67127,97

2.2.1 Вычисление условных средних

X4=0)=44.49

X4=1)=65.36

X4=2)=83.83

X4=3)=103.96

X4=4)=125.65

2.2.2 Вычисление условных дисперсий Y2 для фиксированных значений X4

X4=0)=94.19

X4=1)=63.23

X4=2)=126.56

X4=3)=112.77

X4=4)=96.10

2.2.3 Построение линий регрессии

Приближенная линия регрессии:

Уравнение регрессии имеет следующий вид:

Вычислим коэффициенты и :

(18)

Уравнение регрессии примет вид: = 44.48 + 20.09x.

Рисунок 4 - Линия приближенной регрессии

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

Коэффициенты уравнения вычисляются по средним значениям функции по формулам (17) и(18).

Таблица 7 - Значения

X044,49165,36283,873103,974125,65

Уравнение регрессии примет вид:

Y = 44.48 + 20.09x.

Рисунок 4 - Линия эмпирической регрессии

Приближенная линия регрессии и эмпирическая совпадают.

3. Дисперсионный анализ и планирование эксперимента

Рассмотрим в качестве объекта исследования методы измерения прочности. Выберем ультразвуковой неразрушающий метод измерения прочности бетона.

Ультразвуковой метод основан на связи между скоростью распространения ультразвуковых колебаний в бетоне и его прочностью. Прочность бетона в конструкциях определяют по экспериментально установленным градуировочным зависимостям "скорость распространения ультразвука - прочность бетона" или "время распространения ультразвука - прочность бетона" в зависимости от способа прозвучивания.

Градуировочную зависимость "скорость - прочность" устанавливают при испытании конструкций способом сквозного прозвучивания.

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

Чем выше скорость распространения ультразвука, тем выше прочность бетона.

.1 Краткое описание эксперимента

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

Исследуемые показатели качества метода измерения прочности:

Y1 - точность;

Y2 - скорость ультразвука , 4500-4800 м/с.

Факторы, оказывающие влияние на показатели качества представлены в таблице 8.

Таблица 8 - Факторы и их уровни, изучаемые в эксперименте

ФакторыУровни факторовX1- квалификация персонала (A)- оператор с опытом работы более трех лет (A0); - оператор с опытом работы менее года(A1); - оператор без опыта работы (A2); - студент-практикант (A3); - (A4).X2- калибровка прибора (B)-прибор не прошедший поверку (B0); - калибровка произведена на влажных образцах (B1); - калибровка произведена на сухих образцах (B2); - калибровка прибора произведена не в срок (B3); - (B4).X3-производитель (C)- производитель 1 (C0); - производитель 2 (C1); - производитель 3 (C2); - производитель 4 (C3); - производитель 5 (C4).X4- условия проведения испытания (D)- температура (D0); - влажность (D1); - давление (D2); - загрязнение прибора (D3); - (D4).X5- условия хранения испытуемых образцов (E)- температура, влажность и давление соответствуют указанным в ГОСТ (E0); - температура не соответствует указанной в ГОСТ (E1); - влажность не соответствует ГОСТ (E2), - давление не соответствуют указанному в ГОСТ Е3 - условия не соответствуют (E4).X6 - вид материала, Х7 - содержание примесей; Х8 - возраст бетона; Х9 - плотность; Х10 - водоцементное отношение.Исследуются на постоянном уровне

3.2 Планирование эксперимента

Планирование эксперимента - комплекс мероприятий, направленных на эффективную постановку опытов. Основная цель планирования эксперимента - достижение максимальной точности измерений при минимальном количестве проведенных опытов и сохранении статистической достоверности результатов.

Полную систему ортогональных латинских квадратов для n=p можно построить, используя поля Галуа. Построим поле Галуа вычетов по модулю 5. В поле содержится пять различных элементов 0, 1, 2, 3, 4, 5. Составим таблицу сложения и умножения в этом поле:

СложениеУмножение0 1 2 3 4 1 2 3 4 0 2 3 4 0 1 3 4 0 1 2 4 0 1 2 3 1 2 3 4 2 4 1 3 3 1 4 2 4 3 2 1

Получим второй квадрат путем умножения (К=2)

2*0=0012342*1=2234012*2=4401232*3=1123402*4=334012

3*1=0012343*1=3341123*2=1123403*3=4401233*4=223401

Получим четвертый квадрат (К=4)

4*0=0012344*1=4401234*2=3340124*3=2234014*4=112340

При 5 уровнях в план можно ввести n+1 фактор. Число степеней свободы остаточной суммы будет равно нулю. Такие планы называют насыщенными. Построим насыщенный квадрат для n=5. Наложим для этого четыре полученных выше ортогональных латинских квадрата 55, составляющих полный ряд ортогональных латинских квадратов 55. Полученный план является насыщенным, так как число степеней свободы остаточной суммы, определяемое по формуле где k-число изучаемых факторов, равно нулю.

По этой модели проводится эксперимент из 25 опытов. Такая модель наиболее оптимальна для данного количества факторов, то есть позволяет учесть максимальное количество сочетаний факторов при минимальном количестве опытов.

Таблица 8 - Гипер-греко-латинский квадрат четвертого порядка

АВ012340C=0 D=0 E=0 F=0C=1 D=1 E=1 F=1C=2 D=2 E=2 F=2C=3 D=3 E=3 F=3C=4 D=4 E=4 F=41C=1 D=2 E=3 F=4C=2 D=3 E=4 F=0C=3 D=4 E=0 F=1C=4 D=0 E=1 F=2C=0 D=1 E=2 F=32C=2 D=4 E=1 F=3C=3 D=0 E=2 F=4C=4 D=1 E=3 F=0C=0 D=2 E=4 F=1C=1 D=3 E=0 F=23C=3 D=1 E=4 F=2C=4 D=2 E=0 F=3C=0 D=3 E=1 F=4C=1 D=4 E=2 F=0C=2 D=0 E=3 F=14C=4 D=3 E=2 F=1C=0 D=4 E=3 F=2C=1 D=0 E=4 F=3C=2 D=1 E=0 F=4C=3 D=2 E=1 F=0Таблица 9 - Матрица эксперимента n=5, N=25

Номер опытаАВСDEFX7X8X9X10100000011112011111111130222221111403333311115044444111161012341111711234011118123401111191340121111101401231111112024131111122130241111132241301111142302411111152413021111163031421111173142031111183203141111193314201111203420311111214043211111224104321111234210431111244321041111254432101111

3.3 Проведение модельного эксперимента

Проводим модельные эксперименты с выбранными значениями факторов, фиксируя остальные факторы. Заполняем матрицу эксперимента.

Таблица 10 - Матрица эксперимента n=5, N=25

Номер опытаАВСDEFX7X8X9X10Y1Y21000000111123,9726,372011111111131,6730,373022222111151,5748,374033333111187,0795,4750444441111209,57147,576101234111184,6763,3771123401111173,30143,678123401111195,9782,179134012111176,3362,3710140123111155,7763,4711202413111151,5744,1712213024111173,8350,37132241301111177,63109,97142302411111148,20126,2715241302111199,10117,37163031421111237,10133,3717314203111194,8335,8718320314111175,8378,87193314201111221,47199,37203420311111105,1085,37214043211111240,33122,87224104321111206,17178,37234210431111242,30134,1724432104111196,2781,47254432101111261,27192,47

3.4 Дисперсионный анализ

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

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

Метод дисперсионного анализа основывается на следующих предпосылках: распределение исходных случайных величин нормально; дисперсии экспериментальных данных одинаковы для всех условий эксперимента (т. е. для экспериментов, выполненных на различных уровнях изучаемого фактора).

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

3.4.1 Проверка применимости дисперсионного анализа

- проверим выборки на наличие грубых погрешностей по правилу "трех сигм" по формуле (13):

Для ПК Y1:

При проверке условия (13) для ПК Y1 грубых погрешностей не обнаружено.

Для ПК Y2:

При проверке условия (13) для ПК Y2 грубых погрешностей не обнаружено.

- проверка нормальности полученной выборки

Воспользуемся критерием Дэвида-Хартли-Пирсона.

По формуле (6) вычисляем значение критерия для ПК Y1:


R=Ymax-Ymin=261,26-23,96=237,00

Гипотеза нормальности принимается, если

По таблице 75 [4] для n=25 и находим ;

,34<3,43<4,71 следовательно условие выполняется, гипотезу о нормальности распределения принимаем на уровне значимости 0,05.

По формуле (6) вычисляем значение критерия для ПК Y2:


R=Ymax-Ymin=199,37-26,37=172,00

По таблице 75 [4] для n=25 и находим ;

,34<3,44<4,71 следовательно условие выполняется, гипотезу о нормальности распределения принимаем на уровне значимости 0,05. следовательно гипотезу о нормальности распределения принимаем.

- проверка гипотезы о неразличимости дисперсий

Условие однородности дисперсий - это равенство дисперсий. Его проверим по критерию Кохрена:

(18)

Оценка дисперсии определяется по формуле

(19)

Сформируем 5 выборок из ПК Y1 и Y2

Таблица 11 - Оценка дисперсии для ПК Y1

Выборки1234523,9731,6751,5787,07209,5784,67173,3095,9776,3355,7751,5773,83177,63148,2099,10237,1094,8375,83221,47105,10240,33206,17242,3096,27261,27S2у10765,815187,876278,833618,327335,37< G кр(5;5;0,05)=0,63

Так как условие Gp<G кр выполняется, то гипотезу об однородности дисперсий принимаем.

Таблица 12 - Оценка дисперсии для ПК Y2

Выборки1234526,3730,3748,3795,47147,5763,37143,6782,1762,3763,4744,1750,37109,97126,27117,37133,3735,8778,87198,3785,37122,87178,37134,1781,47192,47S2у2275,834680,071066,382832,192601,50

G кр(5;5;0,05)=0,63

Так как условие Gp<G кр выполняется, то гипотезу об однородности дисперсий для ПК Y2 принимаем.

Делаем вывод, что мы можем воспользоваться дисперсионным анализом.

3.4.2 Алгоритм дисперсионного анализа

1) определяем итоги для всех факторов на каждом уровне

Таблица 13 -Итоги по факторам А, B, C, D, E, F

Итоги фактора А 01234Y1403,8333486,0333550,3333734,33331046,333Y2637,63579,80643,30629,33730,80Итоги фактора В 01234Y1348,15415,05448,15531,85709,35Y2390,15438,65453,55563,95606,25Итоги фактора С 01234Y1509,93679,20477,8755,2333798,7Y2473,35543,65403,05553,85478,65Итоги фактора D 01234Y1521,53598,43640,5333675,6333784,7333Y2358,65418,65466,35558,25650,65Итоги фактора Е 01234Y1410,13496,6667642,9667660,63331010,467Y2343,25408,25483,45532,55685,05Итоги фактора F 01234Y1857,63621,2667670,2667531,5333540,1667Y2670,85447,05539,85373,15421,65

2) Сумма квадратов всех наблюдений:

, (20)

) Сумма квадратов итогов по каждому из факторов:

, (21)

, (22)

, (23)

, (24)

, (25)

, (26)

) Квадрат общего итога, деленный на число всех наблюдений (корректирующий член)

, (27)

) Суммы квадратов, обусловленные фактором:

SSA = SS2 - SS8, (28)B = SS3 - SS8, (29)

SSD = SS5 - SS8, (31)E = SS6 - SS8, (32)F = SS7 - SS8, (33)

6) Общая сумма квадратов:

SSобщ = SS1 - SS8, (34)

) Остаточная сумма квадратов:

SSост = SSобщ - ΣSSфакт, (35)

) Дисперсия фактора:

, (36)

) Дисперсия ошибки:

, (37)

10) Дисперсионное отношение:

, (38)

При расчете SSост исключим из расчета фактор В для Y1 и фактор С для Y2.

Таблица 14 - Результаты расчетов Y1и Y2

Y1Y2SS1550086,46301007,27SS2467247,09256071,04SS3417341,66247183,42SS4431587,31243584,10SS5422537,85251273,69SS6457154,51254223,11SS7429014,46251672,22SS8414959,28240600,06SSA52287,8015470,98SSB2382,376583,36SSC16628,032984,04SSD7578,5610673,63SSE42195,2313623,05SSF14055,1811072,16SS общ135127,1760407,21SSост2382,372984,04S²A13071,953867,74S²B595,591645,84S²C4157,01746,01S²D1894,642668,41S²E10548,803405,76S²F3513,792768,04S²ош595,59746,01

Значимость факторов проверяют по критерию Фишера:

Д.О. = > F 1- p (f1 ,f2) (39)

где р - уровень значимости;

f1 - число степеней свободы фактора, (f1 = 4);

f2 - число степеней свободы остатка, (f2 = 4).

Если дисперсионное отношение удовлетворяет неравенству, то принимаются нулевые гипотезы. Если дисперсионное отношение будет больше табличного, то соответствующая нулевая гипотеза отвергается, и влияние фактора считается значимым.

По таблице 7[2] значение критерия Фишера для всех факторов Y1, Y2:

F 0,05(4;4)=6.4

Таблица15- Итоговая таблица дисперсионного анализа для Y1

Источник дисперсииЧисло степеней свободы fСумма квадратовСредний квадратД.О.A452287,8013071,9521,95B42382,37595,591,00C416628,034157,016,98D47578,561894,643,18E442195,2310548,8017,71F414055,183513,795,90Остаток 42382,37595,59Общая сумма24135127,17

На ПК Y1 оказывают влияние факторы: А, С, E.

Таблица16- Итоговая таблица дисперсионного анализа для Y2

Источник дисперсииЧисло степеней свободы fСумма квадратовСредний квадратД.О.A415470,983867,745,18B46583,361645,842,21C42984,04746,011,00D410673,632668,413,58E413623,053405,764,57F411072,162768,043,71Остаток42984,04746,01Общая сумма2460407,21

На ПК Y2 факторы влияния не оказывают.

3.4.3 Вычисление множественного рангового критерия Дункана

Вычисление критерия Дункана будем производить по тем показателям качества, для которых обнаружено факторное влияние.

В качестве нулевой примем гипотезу Н0: отсутствие значимой разности.

При вычислении критерия Дункана необходимо:

проранжировать средние значения, расположив их в порядке возрастания

определить ошибку воспроизводимости результатов с соответствующим числом степеней свободы (таблица 15, 16);

определить ошибку для каждого среднего по формуле:

(40)

где S2ош - ошибка воспроизводимости результатов.

SY1(А) = SY1(С) = SY1(E)= 10,91

выписать из таблицы критерия Дункана значения рангов числом k-1 (где k - количество градаций фактора) с выбранным уровнем значимости 0,05 и числом степеней свободы f=4 и умножаем эти ранги на соответствующие ошибки средних.

Таблица 17 - Параметры влияющих факторов

ФакторыПараметрыЗначенияГрадации01234АСреднее значение80,7697,20110,06146,86209,26Ранги3,984,014,024,02РангSy43,4343,7643,8743,87СГрадации20134Среднее значение95,56101,98135,84151,04159,74Ранги3,984,014,024,02РангSy43,4343,7643,8743,87ЕГрадации01234Среднее значение82,02699,33128,59132,12202,09Ранги3,984,014,024,02РангSy43,43843,7643,8743,87

вычисляем разности средних по градациям и сравниваем полученные значения с соответствующими рангами для нахождения значимости разностей.

А4-A0642,50>43,87различие значимоеA4-A1560,30>43,87различие значимоеA4-A2496,00>43,76различие значимоеA4-A3312,00>43,43различие значимоеA3-A0330,50>43,87различие значимоеA3-A1248,30>43,76различие значимоеA3-A2184,00>43,43различие значимоеA2-A0146,50>43,76различие значимоеA2-A164,30>43,43различие значимоеA1-A082,20>43,43различие значимоеC4-C064,18>43,87различие значимоеC4-C1119,50>43,87различие значимоеC4-C2320,90>43,76различие значимоеC4-C343,47>43,43различие незначимоеC3-C0245,30>43,87различие значимоеC3-C176,03>43,76различие значимоеC3-C2277,43>43,43различие значимоеC2-C040,28>43,76различие незначимоеC2-C133,85>43,43различие незначимоеразличие значимоеC1-C06,43>43,43различие незначимоеE4-E0600,33>43,87различие значимоеE4-E1513,80>43,87различие значимоеE4-E2367,50>43,76различие значимоеE4-E3349,83>43,43различие значимоеE3-E0250,50>43,87различие значимоеE3-E1163,97>43,76различие значимоеE3-E217,67>43,43различие значимоеE2-E0232,83>43,76различие значимоеE2-E1146,30>43,43различие значимоеразличие значимоеE1-E086,53>43,43различие значимое

Анализ по критерию Дункана показал, что следующие градации влияющих факторов:

градации 0, 1, 2, 3, 4 фактора А для показателя качества Y1 являются значимыми

градации 0, 1, 2, 3, 4 фактора C для показателя качества Y1 являются значимыми

градации 0, 1, 2, 3, 4 фактора E для показателя качества Y1 являются значимыми

. Регрессионный анализ

Получим зависимость точности метода измерения прочности от факторов: А, С, E. Вычислим

z0j = (zmaxj + zminj)/2 (41)

∆zj = (zmaxj - zminj)/2 (42)

xj = (zj - z0j)/ ∆zj (43)

Составим матрицу планирования.

Таблица 18 - Матрица планирования

№z1z2z3z4x1x2x3x4Y212345678

Составим расширенную матрицу планирования

Таблица 19 - Расширенная матрица планирования

№x0x1x2x3x4x12x13x14x23x24x34x1234Y212345678

Уравнение регрессии с коэффициентами взаимодействия:

y = b0+b1∙x1+b2∙x2+b3∙x3+b4∙x4+b12∙x12+b13∙x13+b14∙x14+

+b23∙x23+b24∙x24+b34∙x34+b1234∙x1234 (44)

Список литературы

1. Воскобойников Ю. Е. Математическая статистика: учебное пособие/ Ю. Е. Воскобойников, Е. И. Тимошенко. - Новосибирск: НГАСУ, 2000.

. Гмурман В.Е. Теория вероятности и математическая статистика: Учеб. пособие для вузов /В.Е Гмурман. - М.: Высш. шк., 1997. - 480 с.

. Сергеев А. Г. Метрология: учебное пособие/А. Г. Сергеев, В. В, Крохин. - М.: Логос, 2000. - 408 с.

. Кобзарь А. И. Прикладная математическая статистика. Для инженеров и научных работников / А. И. Кобзарь - М.:ФИЗМАТЛИТ, 2006 - 816 с.

. ГОСТ Р 1.5. Государственная система стандартизации РФ. Стандарты. Общие требования к построению, изложению, оформлению, содержанию и обозначению. - М.: ИПК Изд-во стандартов, 2004

Похожие работы

 

Не нашел материал для курсовой или диплома?
Поможем написать качественную работу
Без плагиата!