Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
1627.pdf
Скачиваний:
43
Добавлен:
07.01.2021
Размер:
1.43 Mб
Скачать

где j – номер столбца в матрице планирования; i – номер строки. Дисперсии коэффициентов уравнения регрессии следующие:

 

Sbj2 Sвосп2

n

Sbjj/ 2

Sвосп2

n

 

/ji2 ,

 

 

x2ji ,

x

 

 

 

 

 

i 1

 

 

 

 

i 1

 

 

 

2

2

n

 

2

2

 

2

 

k

2

2

(4.67)

 

 

 

 

 

 

Sbju

Sвосп

xjixui ,

Sb0

Sb/ 0

xjjSb/ jj .

 

 

i 1

 

 

 

 

 

 

j 1

 

 

 

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

сии, постоянна для всех точек, находящихсяИна равном расстоянии от центра эксперимента. Экспериментатору заранее неизвестно, где на-

4.5.2. Ротатабельные планы второго порядка

Ротатабельным называют планирование, для которого диспер-

ходится та часть поверхности откликаД, которая представляет для него особый интерес, поэтому следует стремиться к тому, чтобы количест-

сия параметра оптимизации ~y , предсказанного уравнением регрес-

во для всех равноотстоящихботАцентра эксперимента точек. Действительно, удаление от центра точек 5-8 в 2 раз меньше, чем удаление точек 1-4 и, следовательно, коэффициенты уравнения регрессии оп-

во информации, содержащееся в уравнении регрессии, было одинако-

ределяются с разл чнойид сперсией. Бокс и Хантер предложили ротатабельные планы 2-го порядка. Для того чтобы композиционный план

был ротатабельным, величину звёздного плеча выбирают из условия

С

 

k 1

 

 

 

k

 

 

 

 

2

4

при k 5

и 2

4

 

при k 5.

(4.68)

 

 

 

Значения звёздных плеч и числа точек в центре ротатабельных планов приведены в табл. 4.22.

Матрица планирования для ротатабельного плана представлена табл. 4.23.

88

Таблица 4.22

Значения звёздных плеч в ротатабельных планах второго порядка

Параметр

 

 

Значения параметров при числе независимых факторов

плана

 

 

2

 

 

3

 

4

 

5

5

 

6

 

 

6

 

7

7

Ядро плана

22

 

 

23

 

24

 

25

25-1

 

26

 

 

26-1

 

27

27-1

Звёздное

1,414

 

1,682

2,000

 

2,378

2,000

2,828

 

2,378

 

3,333

2,828

плечо

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Число

точек

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в

центре

5

 

 

6

 

7

 

10

6

 

15

 

 

9

 

21

14

плана n0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ротатабельный план второго порядка

 

 

 

 

Таблица 4.23

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Номер опыта

 

 

 

 

 

 

 

 

Факторы

 

 

 

 

 

 

 

 

Результат

 

 

 

 

 

 

X0

 

X1

 

X2

X1X2

 

X12

 

 

X22

 

 

yi

Ядро плана

 

1

 

 

+1

 

-1

-1

+1

 

+1

 

 

 

+1

 

 

y1

 

 

 

2

 

 

+1

 

+1

-1

-1

 

+1

 

 

 

+1

 

 

y2

 

 

 

3

 

 

+1

 

-1

+1

-1

 

+1

 

 

 

+1

 

 

y3

 

 

 

 

 

 

 

 

 

 

И

 

 

 

 

 

 

 

 

4

 

 

+1

 

+1

+1

+1 +1

 

+1

 

 

y4

Звёздные

 

5

 

 

+1

+1,414

0

0

 

+2

 

 

 

0

 

 

 

y5

точки

 

 

6

 

 

+1

-1,414

0

0

 

+2

 

 

 

0

 

 

 

y6

 

 

 

 

 

 

 

 

 

 

 

 

Д

 

 

 

 

 

 

 

 

 

 

 

 

7

 

 

+1

 

0

+1,414

0

 

0

 

 

 

+2

 

 

y7

 

 

 

8

 

 

+1

 

0

-1,414

0

 

0

 

 

 

+2

 

 

y8

 

 

 

 

 

 

 

 

 

А

 

 

 

 

 

 

 

 

 

 

 

Центр

 

 

9

 

 

+1

 

0

0

0

 

0

 

 

 

0

 

 

 

y9

плана

 

 

10

 

 

+1

 

0

0

0

 

0

 

 

 

0

 

 

 

y10

 

 

 

 

 

 

 

 

б

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11

 

 

+1

 

0 0

0

 

0

 

 

 

0

 

 

 

y11

 

 

 

12

 

 

+1

 

0

0

0

 

0

 

 

 

0

 

 

 

y12

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

13

 

 

+1 0

0

0

 

0

 

 

 

0

 

 

 

y13

УчитываяСспецифический характер ротатабельного плана в общем виде, можно получить формулы для расчёта коэффициентов уравнения регрессии и их дисперсий:

 

 

 

 

 

 

 

 

A

 

 

k

 

 

 

 

 

 

 

b

 

 

 

2 2 k

2 0y 2 c jjy ;

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

n

 

 

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

bj c/n jy ;

 

 

 

 

A

 

k 2 k jjy c

k

 

b

jj

 

 

 

c2

2 1 jjy 2 c 0y ;

 

 

 

 

n

 

 

 

 

 

 

 

j 1

 

 

 

 

 

 

 

 

 

bju

c2

juy ;

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.69)

(4.70)

(4.71)

(4.72)

89

 

 

 

 

 

 

 

 

2

 

 

2A 2

k 2 2

 

 

 

 

 

 

 

 

 

 

 

Sb0

 

 

 

 

 

 

 

 

 

Sвосп

;

 

(4.73)

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

2

 

 

A k 1 k 1 c2

2

 

 

 

 

 

 

 

 

Sbjj

 

 

 

 

 

 

 

 

 

 

 

 

 

Sвосп ;

(4.74)

 

 

 

 

 

 

 

 

 

n

 

 

c2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Sbju

 

 

 

 

Sвосп ,

 

 

(4.75)

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

n

 

 

 

 

 

 

 

n

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

jy

 

где 0y X0i yi ;

 

juy X ji Xui yi ;

 

 

X ji yi

;

 

i 1

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

i 1

 

 

n

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

jjy X 2 y

 

;

 

c

 

 

;

 

 

A

 

 

;

i

 

n

 

 

 

2 k 2 k

 

ji

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i 1

 

 

 

 

 

X 2ji

 

 

 

 

 

 

 

 

nk

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

n

n n .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k 2 n1

 

 

 

1

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Матрица ротатабельного планирования оказывается неортого-

нальной, так как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

X 2

0;

 

 

n

 

X

 

 

0;

j u.

(4.76)

 

 

 

 

X

0i

 

 

X 2

 

2

 

 

 

 

 

i 1

ui

 

 

 

 

 

ji

 

 

ui

 

 

И

 

 

 

 

 

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно, если какой-либо из квадратичных эффектов ока-

 

 

 

 

 

 

 

 

 

 

 

 

 

Д

 

 

 

зался незначимым, то после его исключения коэффициенты уравне-

ния регрессии необходимо пересчитать заново.

 

 

 

 

При использован

 

ротатаАельных планов второго порядка дис-

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

­ находят остаточную сумму квадратов

2

 

 

n

~

2

 

 

(4.77)

С S1

 

yi yi

 

 

 

 

 

i 1

 

k 2 k 1

 

 

с числом степеней свободы d.f

1

n

;

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

­ по опытам в центре плана определяют дисперсию воспроиз-

водимости

 

n0

 

 

 

 

 

 

2

 

~

2

 

(4.78)

S2

 

y0i y0i

 

i1

счислом степеней свободы d.f2 n0 1;

­находят сумму квадратов, характеризующих неадекватность

90

S2

S2

S2

(4.79)

3

1

2

 

счислом степеней свободы d.f3 d.f1 d.f2 ;

­проверяют адекватность по F-критерию:

S2

/d.f

3

 

 

F

3

 

,

(4.80)

 

/d.f2

S22

 

 

уравнение адекватно, если F FТ .

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

1.Как зависит число опытов отДвида Ипринимаемой математической модели?

2.Чем можно объяснитьАширокое распространение полиномиальных моделей?

3.Дайте определениебполного факторного эксперимента.

4.Что характеризуют -коэффициенты?

5.Перечислитеиэтапы планирования и реализации полного факторного эксперимента.

6.Что называютСкод рованием факторов? Зачем его проводят?

7.Геометрическое представление планов типа 2k . Постройте

9.Перечислите свойства матрицы планирования полного факторного эксперимента.

10.Что называют рандомизацией опытов? Зачем ее проводят?

11.Какие опыты называют параллельными?

12.Как и для чего проводится проверка однородности дисперсии параллельных опытов?

13.Что означает понятие воспроизводимости эксперимента?

14.Как оценить ошибку эксперимента?

15.Какой метод применяется при расчете коэффициентов уравнения регрессии? Запишите формулу расчета b-коэффициентов.

91

16. Что называют взаимодействием факторов и как оно учитывается при планировании полного факторного эксперимента?

17. Что называют взаимодействием первого, второго, третьего и т.д. порядка? Как определяется число возможных взаимодействий

факторов?

 

 

18.

Какиесуществуютспособыпроверкизначимостиb-коэффициентов?

19.

Чем может быть обусловлена незначимость коэффициентов

уравнения регрессии?

 

 

20.

Как и для чего проводится проверка адекватности уравнения

регрессии?

 

 

21.

Что называют дробным факторным экспериментом?

22.

Дайте определение дробной реплики полного факторного

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

 

 

23.

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

24.

Какие планы называют насыщенными?

25.

Явление смешивания оценок -коэффициентов в дробном

26.

Что называют генерирующим соотношением и определяю-

щим контрастом?

 

И

 

 

27.

Ортогональные и ротатабельные планы второго порядка.

28.

 

 

Д

Определение звёздных плеч и количества опытов в центре

планов второго порядка.

 

 

29.

Определен е коэфф циентовАуравнения регрессии и провер-

ка их значимости в ортогональных и ротатабельных планах.

30.

 

б

 

Оценка адекватности модели, построенной с помощью пла-

нов второго порядка.

 

 

 

и

 

5. ПЛАНИРОВАНИЕ ЭКСПЕРИМЕНТОВ ПРИ ПОИСКЕ

 

ОПТИМАЛЬНЫХ УСЛОВИЙ

 

С

 

 

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

92

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]