Файл: Кориков А.М. Математические методы планирования эксперимента учеб. пособие.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 27.06.2024

Просмотров: 173

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

- 9 5 -

3 . Все коэффициента регрессии незначимы.

В каждом варианте оптимум может быть близко, далеко

или о его положении нет информации (неопределенная ситуа­

ция).

Рассмотрим п е р в ы й вариант.

Если область оптимума близка, возможно 3 решения: окон­

чание исследования, переход к планам второго порядка и дви­

жение по градиенту.

Решение при неопределенной ситуации или удаленной об­

ласти оптимума одно и то же: движение по градиенту.

Известно, что движение по градиенту наиболее эффективно,

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

те выбираются решения, реализация которых приводит к полу­

чению значимых коэффициентов:

1)изменение интервалов варьирования;

2)перенос центра плана ;

3)отсеивание незначимых факторов }

4)увеличение числа параллельных опытов ;

5 ) достройка плана.

Достройка плана осуществляется несколькими способами.

1 . Методом "перевала" - у исходной реплики изменяют зна­

ки на обратные.

2 . Переходом к полному факторному эксперименту.

3 . Переходом к реплике меньшей дробности. . -

4 . Переходом к плану второго порядка (если область опти­ мума близка).

Реализация любого из этих решений требует больших аксяе-

- 9 6 -

риментальшх усилий. Поэтому иногда стоит пойти на некото­

рый риск и двигаться только

по

значимым факторам. '

Рассмотрим т р е т и й

вариант: все коэффициенты

регрессии незначимн (кроме

Ь0

) . Чаще всего это происходит

из-за большой ошибки эксперимента или из-за узких интерва­

лов варьирования. Поэтому возможные решения направлены на увеличение точности эксперимента путем улучшения методики

исследования или постановкой параллельных опытов, а тайке на расширение интервалов варьирования. Если область оптиму­ ма близка, то возможно также окончание исследования или

построение плана второго

порядка.

 

Рассмотрим принятие

решений в случае

н е а д е к в а т ­

н о й л и н е й н о й

п о д е л и . Если

линейная модель

неадекватна, значит не удается аппроксимировать поверхность отклика плоскостью. Формальные признаки (кроме величины

Б -критерия), по которым можно установить неадекватность

линейной

модели,

следующие;

 

 

 

 

Значимость

хотя бы одного эффекта

взаимодействия.

2 . Значимость суммы коэффициентов регрессии при квадра­

тичных членах Т

л . • . Оценкой этой

суммы служит разность

j

с 1

"

 

 

 

 

 

между 0 0

и значением

в

центра

плана

i j a . Ее ни

разность

превосходит ошибку опыта,

то

гипотеза о незначиыости

I I

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

Если область оптимума близка, то либо доследование з*—


- 97 -

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

Наконец, если поставлена задача построения интерполяцион­ ной фор*7лы, то на получении адекватной модели исследование

заканчивается, а в случае неадекватной модели принимается од-

*

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

них.

4 - 3 . 2 . П р и н я т и е р е ш е н и й п о с л е к р у ­ т о г о в о с х о ж д е н и я

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

- 98 -

учитывать адекватность (или неадекватность) линейной модели.

Об эффективности дваяеяш по градиенту судят по величи­

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

Если

к р у т о , е

в о с х о ж д е н и е

э ф ф е к т и в ­

н о й область оптимума

близка, то возмозность 2 решения:

окончание

исследования

и достройка линейного

плана

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

построению новой линейной модели.

 

Если

к р у т о е

в о с х о ж д е н и е

н е э ф ф е к ­

т и в н о

и область

оптимума близка, наиболее

типичные реше­

ния: окончание исследования (выбирается лучший опыт) или пост­ роение плана второго порядка для описания области оптимума.

Если линейная модель была неадекватна, то возможно и третье решение: выяснение причин неадекватности линейной модели.

Если область оптимума далека и линейная модель адекватна, но все же крутое восхождение оказалось неэффективным, то воз­ можное объяснение - в характере поверхности отклика (см.рис.

4 . 2 ) .


- 99 -

1 л

Рис. 4 . 2 .

I - исследованная область факторного простран­ ства в 1-ом цикле крутого восхождения ; П - исследованная область факторного простран­ ства во 2-ом цикле крутого восхождения.

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

Если область оптимума далека, линейная модель неадекват­ на я крутое восхождение неэффективно, то следует выяснить причины неадекватности модели.

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

ступить к построению плана второго порядка или кончать иссле­ дования.


- 1 0 0 -

§ 4 . 4 . Ч и с л о в о й

п р и м е р

с

и с п о л ь ­

з о в а н и е м

м е т о д а

 

к р у т о г о

в о с х о ж д е н и я

[32]

Оптимизацию процесса с применением метода крутого вое -

хождения начинают с получения линейного уравнения регрес­

сии. В этом случае исследователя интересуют в основном ли­

нейные члены, следовательно, целесообразно использовать

дробную реплику.

Рассмотрим процесс, оптимизируемый по четырем факторам.

Пусть

это будет

хл,

хг ,

х 5 ,

х 4 . Параметром оптимиза­

ции

будем считать

выход

готовой

продукции. Это - задача

на максимум с предельным значением параметра оптимизация,

равным 100%.

Поскльку число коэффициентов линейного уравнения при

Я = 4-равно пяти, можно использовать дробную реплику, со ­

держащую восемь точек. Для трехфакторного эксперимента это

будет полная реплика.

В связи с необходимостью получить несмешанные линейные

эффекты целесообразно использовать дробную реплику с опреде­

ляющим контрастом, равным I = х.,хгхьх^,

В соответствии с

этим матрица планирования будет иметь вид, представленный в

табл. 4 . 2 . Натуральные значения факторов X j приведены в

табл. 4 . 1 .

Результаты экспериментов приведены в табл .4.3. Для оцен­

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

по три параллельных опыта.


Таблица 4 . 1

Уровни

 

 

 

 

r

 

 

Х 1

! х г

 

s . i

* 4

§акторов_ I

 

+1

! 4,0 j 7,8

j 12,2!

5,3

0

j 3|2 • 6,3

 

Ю,Щ

4 , 1

!

2,4

j 4,8

 

8,2i

2,9

A i

!

0,8

i 1,5 i

2,0.

1,2

 

1

 

i

1

L

 

-•101 -

 

 

 

 

Таблица 4 П 2

*

i

т

-I

i * z j

ТОЧКИ

1

. } — —

I

j

+ I

! + l

i + i

2 j +1 ! - i ! +1 i - i

3 i - I i +1 j +1

i - I

4

! - I ! - I i +1

i + i

5

j +1, ! + I

i - i

 

 

6

j +1

I " 1

1 -

I

i - I

7

i - I i +1

- I

 

8

i - I

i - I

- J

 

 

Статистическую обработку результатов экспериментов начи­ наем с расчета дисперсии воспроизводимости. Вычисляем среднее арифметическое Ць по данным для параллельных измерений ; квадраты разностай между средним арифметическим и результатага! параллельных измерений ( у - Ц}г ; дисперсии воспроизводимоста з каздой точке Sэ 11 Значения * 1 занесены

табл. 4 . 3 .

 

 

 

 

 

 

 

 

 

Таблица 4.3

й

!

Параллельные

изме-Г

 

 

 

 

точ- т

 

 

 

 

4 i

I W ( f c - W 4 r ? v i

 

кя

Г

 

Г'

Чъ

5

 

1

¥i

I

I

 

 

 

1

 

I

i 68,15

 

66,50

65,90

66,85

1,690

!0,120

j 0,903

1,35

2

! 68,90 j 65,90

66,50

1 67,10

3,240

1 1 , 4 4 0

! 0,360

2,52

3

i 6

1 - 1 5

i

61,40

58,30

60,35

0,640

1,100

! 4,200 [2,93

4

i 62,12

 

61,50

58,60

60,74

1,900

0,578

! 4,580

3,03

5

j 72,00 ! 68,85

70,35

70,40

2,560

2,400

j0,003

2,48

6

i

71,10 ! 68,40

72,30

70,60

0,250

4,840

t 2 , 8 9 0

3,99

7

• 64,90 j 65,00

61,80

63,90

1,000

1*210 i 4^410 3,31 "

8

j 61,40 j58,80

61,90

60,70

0,490

? , 6 I 0

i 1 , 4 4 0

2,77