Планирование и организация эксперимента

  Отняв от первого уравнения второе, получаем что 

  Из  чего следует, что коэффициент а₁=5,14

  Чтобы найти а₀ подставим найденный коэффициент в первое уравнение

6а₀+3,98*5,14= -38,73

6а₀= -20,465-38,73

6а₀= -59,195

а₀= -9,87

  Рассчитаем  экспериментальное значение для  каждого Y

Y₁= -9,87+5,14/1= -4,72

Y₂= -9,87+5,14/1,25= -5,75

Y₃= -9,87+5,14/1,5= -6,44

Y₄= -9,87+5,14/1,75= -6,93

Y₅= -9,87+5,14/2= -7,29

Y₆= -9,87+5,14/2,25= -7,58

  Построим  график

  Найдем  погрешность каждого эксперимента по формуле

  Е_i=|(y_т-у_р)/у_max|*100%

Ε₁=|(-4,76+4,72)/-4,72|*100%= 0,76%

Е₂=|(-5,59+5,75)/-5,59|*100%= 2,86%

Е₃=|(-6,76+6,44)/-6,44|*100%= 4,97%

Е₄=|(-6,4+6,93)/-6,4|*100%= 8,28%

Е₅=|(-7,12+7,29)/-7,12|*100%= 2,46%

Е₆=|(-8,1+7,58)/-7,58|*100%= 6,86%

Задание 4

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

Т а б л  и ц а 4

 

     Находим среднее значение в каждом опыте Ф-ла 

У_{1 ср} = (-1,18-1,28-1,28)/3= -1,247

У_{2 ср} = (-12,54-12,64-12,69)/3= -12,623

У_{3 ср} = (-15,24-15,22-15,24)/3= -15,233

У_{4 ср} = (1,49+1,5+1,38)/3= 1,457

     Находим построчную дисперсию каждой выходной величины Ф-ла 

S²₁= ((-1,18+1,247)²+(-1,28+1,247)²+(-1,28+1,247)²)/(3-1)=0,0033

S²₂= ((-12,54+12,623)²+(-12,64+12,623)²+(-12,69+12,623)²)/(3-1)=0,0058

S²₃= ((-15,24+15,233)²+(-15,22+15,233)²+(-15,24+15,233)²)/(3-1)=0,0001

S²₄= ((1,49-1,457)²+(1,50-1,457)²+(1,38-1,457)²)/(3-1)=0,0044

     Находим расчетный коэффициент  Кохрэна Ф-ла 

0,0033+0,0058+0,0001+0,0044=0,0137

G_p=0,0058/0,0137= 0,42

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

     Находим коэффициенты по ф-ле 

а₀= (-1,247-12,623-15,233+1,457)/4= -6,913

а₁=(1,247-12,623+15,233+1,457)/4= 1,328

а₂=(1,247+12,623-15,233+1,457)/4= 0,023

а₁₂=(-1,247+12,623+15,233+1,457)/4= 7,017

     Составим новый вид регрессионной  зависимости с учетом полученных  коэффициентов

У= -6,913+1,328Х₁+0,023Х₂+7,017Х₁Х₂

     Находим расчетное значение для  каждого эксперимента

У_{1 р}= -6,913-1,328-0,023+7,017= -1,248

У_{2 р}= -6,913+1,328-0,023-7,017= -12,624

У_{3 р}= -6,913-1,328+0,023-7,017= -15,234

У_{4 р} = -6,913+1,328+0,023+7,017= 1,456

Найдем  погрешность каждого эксперимента по формуле

  Е_i=|(y_т-у_р)/у_max|*100%

Ε₁=|(-1,247+1,248)/-1,247|*100%=0,08%

Е₂=|(-12,623+12,624)/-12,623|*100%=0,008%

Е₃=|(-15,233+15,234)/-15,233|*100%=0,007%

Е₄=|(1,457-1,456)/1,456 |*100%=0,07%

Необходимо  оценить найденные коэффициенты на статистическую значимость с использованием критерия Стъюдента по формуле

f=N*(m-1)

f=4*(3-1)=8

Если  расчетное значение коэффициента Стъюдента будет меньше табличного равного 2,3, то данный коэффициент удаляется из нашей регрессионной зависимости

t_i=|a_k|/S{} 

Найдем оценку дисперсии коэффициента по формуле

S²{}= S²_воспр/N*m          

     

 

S²{}=0,0034/3*4=0,0003

S{}=0,017

t₀=6,913/0,017=406,6

t₁=1,328/0,017=78,1

t₂=0,023/0,017=1,4

t₁₂=7,017/0,017=412,7

Так как  t₂ меньше табличного коэффициента, то он исключается из регрессионной зависимости

У= -6,913+1,328Х₁+7,017Х₁Х₂

Находим расчетное значение для каждого  эксперимента

У₁= -6,913-1,328+7,017=1,224

У₂= -6,913+1,328-7,017= -12,602

У₃= -6,913-1,328-7,017= -15,258

У₄= -6,913+1,328+7,017=1,432

Найдем  погрешность каждого эксперимента по формуле

  Е_i=|(y_т-у_р)/у_max|*100%

Ε₁=|(-1,247+1,224)/-1,224|*100%=1,88%

Е₂=|(-12,623+12,602)/-12,602|*100%=0,17%

Е₃=|(-15,233+15,258)/-15,233|*100%=0,23%

Е₄=|(1,457-1,432)/1,432 |*100%=1,74%

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

  Критерий  Фишера находим по формуле

  F_p=S²_ад/S²_воспр 
 

     

 

F_p=0,01/0,0034=1,94

    Так как расчетный критерий Фишера меньше табличного критерия равного 5,32, то модель адекватна экспериментальным данным и ее можно использовать для дальнейших исследований. 
 
 

Задание 5

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

Т а  б л и ц а

 

     Находим среднее значение в каждом опыте Ф-ла 

У_{1 ср} = (-6,94-6,97-6,97)/3= -6,96

У_{2 ср} = (-12,07-12,12-12,12)/3= -12,103

У_{3 ср} = (-13,26-13,33-13,27)/3= -13,287

У_{4 ср} = (-12,56-12,56-12,48)/3= -12,533

У_{5 ср} = (-9,58-9,51-9,42)/3= -9,503

У_{6 ср} = (-0,75-0,64-0,72)/3= -0,703

У_{7 ср} = (-9,14-9,02-9,07)/3= -9,077

У_{8 ср} = (5,73+5,69+5,62)/3= 5,68 

     Находим построчную дисперсию  каждой выходной величины Ф-ла 

S²₁= ((-6,94+6,96)²+(-6,97+6,96)²+(-6,97+6,96)²)/(3-1)=0,0003

S²₂= ((-12,07+12,103)²+(-12,12+12,103)²+(-12,12+12,103)²)/(3-1)=0,0008

S²₃= ((-13,26+13,287)²+(-13,33+13,287)²+(-13,27+13,287)²)/(3-1)=0,0014

S²₄= ((-12,56+12,533)²+(-12,56+12,533)²+(-12,48+12,533)²)/(3-1)=0,0021

S²₅= ((-9,58+9,503)²+(-9,51+9,503)²+(-9,42+9,503)²)/(3-1)=0,0064

S²₆= ((-0,75+0,703)²+(-0,64+0,703)²+(-0,72+0,703)²)/(3-1)=0,0032

S²₇= ((-9,14+9,077)²+(-9,02+9,077)²+(-9,07+9,077)²)/(3-1)=0,0036

S²₈= ((5,73-5,68)²+(5,69-5,68)²+(5,62-5,68)²)/(3-1)=0,0031 

     Находим расчетный коэффициент  Кохрэна Ф-ла 

(0,0003+0,0008+0,0014+0,0021+0,0064+0,0032+0,0036+0,0031)/8=0,0211

G_p=0,0064/0,0211= 0,3

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

     Находим коэффициенты по ф-ле 

а₀= (-6,96-12,103-13,287-12,533-9,503-0,703-9,077+5,68)/8= -7,311

а₁=(6,96-12,103+13,287-12,533+9,503-0,703+9,077+5,68)/8= 2,396

а₂=(6,96+12,103-13,287-12,533+9,503+0,703-9,077+5,68)/8= 0,007

а₃=(6,96+12,103+13,287+12,533-9,503-0,703-9,077+5,68)/8=3,910

а₁₂=(-6,96+12,103+13,287-12,533-9,503+0,703+9,077+5,68)/8= 1,482

а₁₃=(-6,96+12,103-13,287+12,533+9,503-0,703+9,077+5,68)/8= 3,493