logo
Информатика и КГ_2014

21.3. Модель многомерного объекта

Предположим, что технологический процесс можно описать математической моделью вида

y = b0 + b1x1 +…+ bnxn + b12x1x2 +…+ bn–1,nxn–1xn,

где y – выходной параметр процесса; b0, b1,, bn1,n – искомые неизвестные коэффициенты процесса; x1, , xn – входные параметры процесса.

Соотношения такого вида называются уравнениями регрессии.

Например, для процесса (рис. 21.2), имеющего три входных параметра (фактора), математическая модель примет вид

y = b0 + b1  x1 + b2  x2 + b3  x3 + b12  x1  x2 + b13  x1  x3 + b23  x2  x3.

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

К

Рис. 21.2. Многомерный объект

оличество опытов в эксперименте рассчитывается по формулеN = 2n, где n – количество факторов. Входные воздействия принимают минимальные и максимальные значения. Для упрощения вычислений нужно перейти от физических переменных x1, , xn к кодированным по следующей формуле:

zi = , i = 1, 2, 3.

Здесь xi0 – значение фактора на базовом (нулевом) уровне, равное среднему значению между минимальным и максимальным значениями; ∆xi – интервал варьирования по данному фактору.

В случае трех входных параметров план проведения эксперимента имеет вид, представленный на рис. 21.3.

З

№ п/п

z1

z2

z3

уэ

1

–1

–1

–1

уэ1

2

+1

–1

–1

уэ2

3

–1

+1

–1

уэ3

4

+1

+1

–1

уэ4

5

–1

–1

+1

уэ5

6

+1

–1

+1

уэ6

7

–1

+1

+1

уэ7

8

+1

+1

+1

уэ8

начение –1 соответствует минимальному значению входного параметра, +1 – максимальному значению входного параметра.

В соответствии с методом наименьших квадратов производится вычисление коэффициентов:

Рис. 21.3. План эксперимента

, ,, , , , i=1,2,…,N

Коэффициент регрессии b (b0, b1,, bn1,n.) считается значимым, если выполняется условие

, ,

где Sb – среднеквадратичная ошибка в определении коэффициентов регрессии; tТ – табличное значение критерия Стьюдента, которое выбирается для числа степеней свободы f1 = m – 1.

Для расчета дисперсии воспроизводимости нужно выполнить дополнительно m опытов (m < N) по любой строчке плана, например, при значениях входных факторов на базовом уровне.

В результате получаются дополнительные значения экспериментальных данных yd1, yd2,, ydm.

Тогда

Sy2 = , yoc = , k = 1, …, m.

В табл. 21.1 приведены значения критерия Стьюдента.

Таблица 21.1