
Практическая №3 (Построение множественной регрессии) / Практическая_работа_№3_7b30415975695ee3fe394c5de8c9c5cc
.pdfПРАКТИЧЕСКАЯ РАБОТА №3
ПОСТРОЕНИЕ МНОЖЕСТВЕННОЙ РЕГРЕССИИ
Цель работы: получение умений и навыков в обработке результатов экспериментов методом наименьших квадратов
Теоретические основы
В научной и инженерной практике часто возникает необходимость изучить зависимость функции отклика — количественного изменения изучаемого явления — не от одной, а от нескольких факторов. При проведении экспериментов в такой множественной ситуации приходится записывать измеренные значения функции отклика и всех переменных, влияние которых изучается. Результатами наблюдений являются уже не два столбца значений y и x, как в случае проведения парного регрессионного анализа, а матрица результатов наблюдений. Поэтому применение множественного регрессионного анализа основано матричной записи систем уравнений и решении таких систем.
Прежде, чем приступить к выполнению работы №3 внимательно прочитать параграфе 3.3 раздела 3 лекционного материала. В праграфе 3.3 изложены все методики необходимые для выполнения отдельных частей работы №3. Поэтому необходимо внимательно читать параграф 3.3, не пропуская подпараграфы, и после освоения примера по практическому использованию изложенной информации, выполнить аналогичную процедуру и для исходных данных своего варианта задания.
Общие последовательность и содержание этапов выполнения практической работы №3 приведены ниже.
Порядок выполнения работы
1. Прочитать п. 3.3 раздела 3 лекционного материала.
2.Получить индивидуальный вариант задания у преподавателя и выпи-
сать исходные данные из табл. П3 (см. Приложение к работе 3) со-
гласно заданному варианту.
3.Найти средние значения отклика по результатам двух попыток про-
ведения опыта.
4.Записать систему уравнений подобную (3.16) раздела 3 лекционного материала, но для модели в виде гиперболоида (3.14) (см. план на рис. П3.1в Приложении).
5.Произвести соответствующие замены переменных (x1 → zi1; x2 → zi2; x1x2 → zi3, где i — номер опыта, т. е. соответственно по счѐту сверху вниз в системе (3.13). Записать систему уравнений в виде подобном
(3.17), но для своей модели.
6.Составить матричное уравнение (3.22) для определения 4-х неизвестных — параметров: а0; а1; а2; а3.
7.Определить из уравнения вектор столбец искомых параметров.
8.Построить график полученной функции отклика для рассматривае-
мой области варьирования факторов.
Требования к содержанию отчѐта по работе
Вотчѐте должна быть указана следующая информация:
1.Цель работы и исходные данные по варианту задания.
2.Система уравнений подобная (3.16), но для модели в виде гипербо-
лоида.
3.Система уравнений после замены переменных в виде подобном
(3.17), но для модели в виде гиперболоида.
4.Расчѐты матричного уравнения.
5.Аналитический вид полученного уравнения множественной регрес-
сии.
6.Графическое изображение функции отклика.
7. Выводы.
N.B. Все расчѐты представляются в виде равенства, в котором сна-
чала записывается расчѐтная формула, затем эта же формула, но с под-
ставленными численными значениями переменных, а в конце равенства записывается полученный результат.

ПРИЛОЖЕНИЕ 3
ИНДИВИДУАЛЬНЫЕ ВАРИАНТЫ ЗАДАНИЙ ДЛЯ РАБОТЫ №3
Найти уравнение множественной регрессии в виде в виде гипербо-
лоида:
y ao a1 x1 a2 x2 a3 x1 x2
согласно плану эксперименту, представленному на рис. П3.1.
Рис. П3.1. План эксперимента для моделей в виде плоскости и гиперболоида: курсивом обозначены номера опытов, обычным шрифтом внутри квадрата – значения кодированных факторов
Каждый опыт (1, 2, 3 и 4) проводили дважды, поэтому в качестве значе-
ния отклика y в каждом опыте надо брать среднее арифметическое ре-
зультатов двух попыток (см. табл. П.3).
|
|
|
|
|
|
|
Таблица П.3 |
|
|
|
Исходные данные по вариантам заданий |
||||
Ва- |
Номер |
Значение |
Результат опыта (y) |
Примечание |
|||
риант |
опыта |
факторов |
|
|
|
||
|
|
x1 |
|
x2 |
1-я по- |
2-я попыт- |
|
|
|
|
|
|
пытка |
ка (повтор) |
|
|
1 |
425 |
|
43 |
12,0 |
13,0 |
x1 —сила тока (А); |
|
2 |
550 |
|
43 |
13,0 |
14,0 |
x2 — расход аргона |
1 |
3 |
425 |
|
55 |
10,0 |
9,5 |
(л/мин); |
|
4 |
550 |
|
55 |
11,5 |
10,5 |
y — пористость плазмен- |
|
|
|
|
|
|
|
ного покрытия (%) |
|
1 |
300 |
|
32 |
11,0 |
12,0 |
x1 —сила тока (А); |
|
2 |
450 |
|
32 |
9,0 |
10,0 |
x2 — расход аргона |
2 |
3 |
300 |
|
44 |
10 |
9,0 |
(л/мин); |
|
4 |
450 |
|
44 |
8,0 |
7,0 |
y — пористость плазмен- |
|
|
|
|
|
|
|
ного покрытия (%) |

|
|
|
|
|
|
|
Продолжение табл. П.3 |
Ва- |
Номер |
Значение фак- |
Результат опыта (y) |
Примечание |
|||
риант |
опыта |
|
торов |
|
|
|
|
|
|
x1 |
|
x2 |
1-я по- |
2-я попыт- |
|
|
|
|
|
|
пытка |
ка (повтор) |
|
|
1 |
425 |
|
7 |
12,0 |
13,0 |
x1 —сила тока (А); |
|
2 |
550 |
|
7 |
13,0 |
14,0 |
x2 — расход водорода |
3 |
3 |
425 |
|
10 |
9,0 |
10,0 |
(л/мин); |
|
4 |
550 |
|
10 |
10,0 |
11,0 |
y — пористость покры- |
|
|
|
|
|
|
|
тия (%) |
|
1 |
300 |
|
4 |
11,0 |
12,0 |
x1 —сила тока (А); |
|
2 |
450 |
|
4 |
12,0 |
13,0 |
x2 — расход водорода |
4 |
3 |
300 |
|
7 |
8,0 |
9,0 |
(л/мин); |
|
4 |
450 |
|
7 |
9,0 |
10,0 |
y — пористость покры- |
|
|
|
|
|
|
|
тия (%) |
|
1 |
425 |
|
43 |
2850 |
2830 |
x1 —сила тока (А); |
|
2 |
550 |
|
43 |
3450 |
3410 |
x2 — расход аргона |
5 |
3 |
425 |
|
55 |
2340 |
2380 |
(л/мин); |
4 |
550 |
|
55 |
3020 |
2920 |
y — микротвѐрдость |
|
|
|
||||||
|
|
|
|
|
|
|
плазменного покрытия |
|
|
|
|
|
|
|
(МПа) |
|
1 |
300 |
|
31 |
3050 |
3100 |
x1 —сила тока (А); |
|
2 |
370 |
|
31 |
2990 |
2960 |
x2 — расход аргона |
6 |
3 |
300 |
|
43 |
2760 |
2770 |
(л/мин); |
4 |
370 |
|
43 |
2700 |
2680 |
y — микротвѐрдость |
|
|
|
||||||
|
|
|
|
|
|
|
плазменного покрытия |
|
|
|
|
|
|
|
(МПа) |
|
1 |
425 |
|
7 |
2900 |
2940 |
x1 —сила тока (А); |
|
2 |
550 |
|
7 |
3380 |
3360 |
x2 — расход водорода |
7 |
3 |
425 |
|
10 |
3120 |
3080 |
(л/мин); |
4 |
550 |
|
10 |
3440 |
3420 |
y — микротвѐрдость |
|
|
|
||||||
|
|
|
|
|
|
|
плазменного покрытия |
|
|
|
|
|
|
|
(МПа) |
|
1 |
300 |
|
4 |
2760 |
2780 |
x1 —сила тока (А); |
|
2 |
370 |
|
4 |
2660 |
2620 |
x2 — расход водорода |
8 |
3 |
300 |
|
7 |
2980 |
3080 |
(л/мин); |
4 |
370 |
|
7 |
2930 |
2910 |
y — микротвѐрдость |
|
|
|
||||||
|
|
|
|
|
|
|
плазменного покрытия |
|
|
|
|
|
|
|
(МПа) |
|
1 |
7 |
|
328 |
70 |
60 |
x1 — плотность тока |
|
2 |
15 |
|
328 |
400 |
380 |
(кА/м2); |
|
3 |
7 |
|
338 |
5 |
7 |
x2 — температура элек- |
9 |
4 |
15 |
|
338 |
150 |
130 |
тролита (К); |
|
|
|
|
|
|
|
y — кол-во трещин в |
|
|
|
|
|
|
|
электролитич. покрытии |
|
|
|
|
|
|
|
на ед. площади (мм-2) |

|
|
|
|
|
|
Продолжение табл. П.3 |
|
Ва- |
Номер |
Значение |
Результат опыта (y) |
Примечание |
|||
риант |
опыта |
факторов |
|
|
|
||
|
|
x1 |
x2 |
1-я по- |
2-я по- |
|
|
|
|
|
|
пытка |
пытка |
|
|
|
|
|
|
|
(повтор) |
|
|
|
1 |
8 |
330 |
56 |
46 |
x1 — плотность тока (кА/м2); |
|
|
2 |
14 |
330 |
320 |
360 |
x2 — температура электроли- |
|
10 |
3 |
8 |
340 |
10 |
12 |
та (К); |
|
4 |
14 |
340 |
180 |
160 |
y — кол-во трещин в элек- |
||
|
|||||||
|
|
|
|
|
|
тролитич. покрытии на ед. |
|
|
|
|
|
|
|
площади (мм-2) |
|
|
1 |
5 |
318 |
620 |
580 |
x1 — плотность тока (кА/м2); |
|
|
2 |
15 |
318 |
4 |
2 |
x2 — температура электроли- |
|
11 |
3 |
5 |
324 |
900 |
960 |
та (К); |
|
4 |
15 |
324 |
8 |
12 |
y — кол-во трещин в элек- |
||
|
|||||||
|
|
|
|
|
|
тролитич. покрытии на ед. |
|
|
|
|
|
|
|
площади (мм-2) |
|
|
1 |
6 |
316 |
640 |
660 |
x1 — плотность тока (кА/м2); |
|
|
2 |
14 |
316 |
14 |
12 |
x2 — температура электроли- |
|
12 |
3 |
6 |
326 |
980 |
1060 |
та (К); |
|
4 |
14 |
326 |
28 |
22 |
y — кол-во трещин в элек- |
||
|
|||||||
|
|
|
|
|
|
тролитич. покрытии на ед. |
|
|
|
|
|
|
|
площади (мм-2) |
|
|
1 |
15 |
67 |
165 |
155 |
x1 — межэлектродное рас- |
|
|
2 |
30 |
67 |
135 |
125 |
стояние (мм); |
|
|
3 |
15 |
317 |
228 |
232 |
x2 — расстояние от верх. |
|
13 |
4 |
30 |
317 |
200 |
210 |
торца вертикально располо- |
|
|
|
|
|
|
|
женной втулки (мм); |
|
|
|
|
|
|
|
y — толщина хромового по- |
|
|
|
|
|
|
|
крытия (мкм) |
|
|
1 |
15 |
17 |
153 |
143 |
x1 — межэлектродное рас- |
|
|
2 |
50 |
17 |
95 |
85 |
стояние (мм); |
|
|
3 |
15 |
367 |
235 |
240 |
x2 — расстояние от верх. |
|
14 |
4 |
50 |
367 |
105 |
110 |
торца вертикально располо- |
|
|
|
|
|
|
|
женной втулки (мм); |
|
|
|
|
|
|
|
y — толщина хромового по- |
|
|
|
|
|
|
|
крытия (мкм) |
|
|
1 |
30 |
117 |
130 |
124 |
x1 — межэлектродное рас- |
|
|
2 |
50 |
117 |
85 |
75 |
стояние (мм); |
|
|
3 |
30 |
267 |
160 |
168 |
x2 — расстояние от верх. |
|
15 |
4 |
50 |
267 |
90 |
98 |
торца вертикально располо- |
|
|
|
|
|
|
|
женной втулки (мм); |
|
|
|
|
|
|
|
y — толщина хромового по- |
|
|
|
|
|
|
|
крытия (мкм) |
|
|
|
|
|
|
Продолжение табл. П.3 |
Ва- |
Номер |
Значение |
Результат опыта (y) |
Примечание |
||
риант |
опыта |
факторов |
|
|
|
|
|
|
x1 |
x2 |
1-я по- |
2-я по- |
|
|
|
|
|
пытка |
пытка |
|
|
|
|
|
|
(повтор) |
|
|
1 |
0 |
20 |
154 |
156 |
x1 — расход электролита |
|
2 |
1,0 |
20 |
150 |
144 |
(л/с); |
|
3 |
0 |
360 |
240 |
238 |
x2 — расстояние от верх. |
16 |
4 |
1,0 |
360 |
170 |
176 |
торца вертикально располо- |
|
|
|
|
|
|
женной втулки (мм); |
|
|
|
|
|
|
y — толщина хромового по- |
|
|
|
|
|
|
крытия (мкм) |