Сутність методу Бокса-Уілсона.

Метод Бокса — Вілсона — метод оптимізації активного експерименту шляхом крутого сходження поверхнею відгуку (параметрів оптимізації) до оптимуму, суть якого полягає в наступному: рух у напрямі градієнта за наявності лінійного рівняння моделі здійснюється із центра експерименту послідовними кроками, які пропорційні добутку коефіцієнта регресії кожного фактора на значення його інтервалу зміни.

Застосовується, зокрема, для одержання моделей процесів збагачення корисних копалин та інших технологічних процесів, при гідродинамічних дослідженнях газліфтнихнафтових свердловин. Алгоритм Метод крутого сходження, або метод Бокса — Вілсона, поєднує істотні елементи методу Гаусса — Зейделя і градієнтного методу з методами ПФЕ або ДФЕ. Так, при використанні алгоритму крутого сходження кроковий рух з точки x_k здійснюється в напрямі найшвидшого зростання рівня виходу, тобто по grad , проте, на відміну від градієнтного методу, коректування напряму здійснюється не після кожного наступного кроку, а після досягнення в деякій точці  на даному напрямі часткового екстремуму цільової функції (рис. 1), аналогічно методу Гаусса — Зейделя. Важливою особливістю методу Бокса — Вілсона є також регулярне проведення статистичного аналізу проміжних результатів на шляху до оптимуму. Порядок виконання операцій при пошуку екстремуму за методом крутого сходження такий:

1) провадиться повний або дробовий факторний експеримент з центром у вихідній точці для визначення grad . Результати експерименту піддаються статистичному аналізу, який включає:

а) перевірку відтворюваності експерименту;

б) перевірку значущості оцінок коефіцієнтів b_i лінійної моделі об'єкта;

в) перевірку адекватності утвореної лінійної моделі

y = b₀ + b₁x₁ + ... + b_nx_n

досліджуваному об'єктові.

2) обчислюються добутки  де  — крок варіювання параметра x_i при проведенні ПФЕ, і фактор, для якого цей добуток максимальний, береться як базовий

;

3) для базового фактора вибирають крок варіювання при крутому сходженні ρ, залишаючи старий крок або впроваджуючи дрібніший;

4) визначаються розміри ρ_i за рештою змінних процесу . Оскільки під час руху градієнтом варійовані параметри повинні змінюватися пропорційно коефіцієнтам , які є компонентами вектора grad у(х), то відповідні ρ_i знаходяться за формулою

, де ρ і  завжди додатні, а коефіцієнт b_i береться зі своїм знаком;

5) проводяться уявні досліди, які полягають у завбаченні значень виходу у певних точках  факторного простору (див. рис. 1). Для цього незалежні змінні лінійної моделі об'єкта змінюються з урахуванням  таким чином, щоб зображуюча точка  виконувала кроковий рух у напрямі вектора grad , утвореного в п. 1, займаючи послідовно положення

.

Очевидно, що j-та координата k-ї точки визначається так:

x_k,j = x_1,j + k_ρ,j,j = 1,2,...,n.

Тоді

 робимо підстановку y_zav,k = ky_zav,1 − (k − 1)b₀,k = 1,2,...,m.

або ще зручніше y_zav,k = y_{zav,k − 1} + (y_zav,1 − b₀),k = 1,2,...,m. 6) уявні досліди продо жуються до тих пір, поки виконується нерівність

де y_max — максимально можливий вихід, який визначається з фізичних міркувань;

7) деякі з уявних дослідів (звичайно через кожні 2-3 уявних кроки) реалізуються на об'єкті для перевірки відповідності апроксимації об'єкта утвореним рівнянням (гіперплощиною). Спостережені значення y_exp порівнюються із завбаченими y_zav (див. рис. 1);

8) точка , де в реальному досліді утворено мaксимальне значення виходу, береться за нову початкову точку, і етап крутого сходження, описаний вище, повторюється;

9) оскільки кожен етап крутого сходження наближає зображуючу точку до області екстремуму , де крутість поверхні відклику менша, то для кожного наступного етапу ρбереться рівним або меншим попереднього;

10) пошук припиняється, коли всі коефіцієнти b_i(i = 1,2,...,n) лінійної моделі об'єкта виходять незначущими. Це свідчить про вихід в область екстремуму цільової функції.

⇐ Предыдущая45678910111213Следующая ⇒

Поделиться: