Материал: 3647

Лабораторная работа №3 Отыскание оптимальных условий функционирования системы. Крутое восхождение по поверхности отклика

1. Цель работы

Цель данной работы - приобретение практических навыков нахождения оптимальные условия функционирования системы методом крутого восхождения по поверхности отклика.

2. Общие сведения

Метод отыскания оптимальных условий функционирования процессов был предложен Боксом и Уильсоном в 1951 г. [1, 2]. Они предложили использовать последовательный метод изучения поверхности отклика, напоминающий итерационный метод решения задач вычислительной математики.

Исследователь вначале ставит небольшую серию опытов для локального описания небольшой поверхности отклика полиномом первой степени. Далее он движется по поверхности отклика в направлении градиента линейного приближения. Если одного линейного приближения оказывается недостаточным, исследователь повторно ставит небольшую серию опытов и находит новое направление для движения по поверхности отклика. Такой шаговый процесс движения по поверхности отклика продолжается до тех пор, пока исследователь не попадет в «почти стационарную область», где линейное приближение окажется недостаточным. Здесь ставится большая серия опытов для описания поверхности отклика полиномом второго порядка.

Опишем процедуру поиска экстремума функции отклика. Ставится ПФЭ или, чаще всего, ДФЭ с шагом варьирования Δх_i , i=1,2,..к. Получают линейное описание поверхности отклика:

, (1)

где коэффициенты регрессии представляют собой частные производные:

(2)

Градиент функции отклика y задается выражением:

(3)

где i, j, ….k – единичные векторы (орты) в направлении координатных осей.

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

X_i^’ = X_i⁰ + , i =1,....k, (4)

где с – константа, шаг изменения факторов при крутом восхождении.

При постановке эксперимента всегда приходится переходить к натуральным переменным х_i. В натуральных переменных величина шага принимается равной:

С_i^’ = с* Δх_i. (5)

Надо помнить, что направление движения по градиенту не инвариантно к изменению интервала варьирования независимых переменных Δх_i. Инвариантными остаются только знаки составляющих градиента.

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

Пример

Исследуется система, модель которого описывается уравнением:

y= 200-4*(X₁-5)²-4*(X₂-5)²-4*(X₃-5)²+e.

Ошибка измерений «е» имеет нормальный закон распределений с дисперсией, равной =1. Необходимо найти соотношение факторов, обеспечивающее наибольшее значение показателя «у».

Реализован ПФЭ 2³, в первых строках табл. 1 записаны интервалы варьирования и указаны кодовые обозначения переменных. В строках 7-14 дана матрица планирования.

Таблица 1 – Планирование эксперимента при выборе оптимального

(максимального) значения зависимой переменной

При постановке эксперимента переходили к натуральным значениям переменных Х_i. Переход от натуральных значений переменных к кодированным значениям осуществляется по формуле:

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

Для движения по поверхности отклика в направлении крутого восхождения нужно независимые переменные (факторные переменные) изменять пропорционально величине соответствующих коэффициентов регрессии при натуральных переменных с учетом их знака:

где с – константа, выбирается из соотношения, чтобы шаг изменения факторов при крутом восхождении не превышал величины 1,0.

Пользуясь этими величинами, выполнена серия мысленных опытов 1 и 2 (строки 19-20) для крутого восхождения по поверхности отклика.

Опыт 3, соответствующий строке 21, был реализован и дал хороший результат. Далее был реализован опыт 4 (строка 22), который дал худший результат, чем в предыдущем опыте. Дальнейшее продвижение по градиенту не целесообразно, т.к. приводит к уменьшению зависимой переменной.

В результате крутого восхождения зависимая переменная увеличилось до величины 200,33. При традиционных методах исследования решение этой задачи потребовало бы длительных и дорогостоящих усилий [1, 2].

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

3. Порядок выполнения лабораторной работы

1) Ознакомиться с методикой нахождения оптимальных условий функционирования системы и методом крутого восхождения по поверхности отклика.

2) Получить у преподавателя вариант задания.

3) Составить матрицу планирования полного факторного эксперимента. Рассчитать коэффициенты линейного уравнения регрессии.

4) Определить градиент функции отклика, задав его выражением.

5) Выбрать величину шага для движения по направлению градиента.

6) Выполнить два мысленных шага, убедиться в правильности движения по направлению градиента (значение зависимой переменной «у» должно возрастать).

6) Реализовать модельный опыт, убедиться в возрастании значения зависимой переменной «у».

7) Реализовать следующий модельный опыт, убедиться в том, что попали в «почти» стационарную область поверхности отклика, при котором не стало дальнейшее улучшение зависимой переменной «y».

8) Определить соотношение факторов х₁, х₂, х₃, обеспечивающее наибольшее значение показателя «у».

9) Сделать выводы по результатам нахождения оптимальных условий функционирования системы.

10) Оформить отчет по выполненной работе.

4. Содержание отчета по лабораторной работе

1) Название, цель работы и номер варианта задания.

2) Математическая модель объекта и интервалы варьирования переменных.

3) Матрица планирования эксперимента с наложением ошибок измерений.

4) Вычисления коэффициентов уравнения регрессии в кодированных переменных. Оценка значимости линейных коэффициентов регрессии и адекватности модели.

5) Запись уравнения регрессии в натуральных значениях факторных переменных.

6) Таблица с планированием эксперимента и результатами крутого восхождения по поверхности отклика.

7) Значения факторов х₁, х₂, х₃, обеспечивающих наибольшее значение показателя «у».

8) Выводы по выполненной лабораторной работе.

5. Контрольные вопросы

1) В чем сущность метода крутого восхождения по поверхности отклика?

2) До каких пор продолжается шаговый процесс движения по поверхности отклика?

3) Опишите процедуру поиска экстремума функции отклика.

4) Запишите выражение для градиента функции отклика.

5) Чему равны частные производные функции поверхности отклика, если она локально может быть описана линейным уравнением?

6) Как нужно изменять независимые переменные для движения по поверхности отклика в направлении крутого восхождения?

7) При движении по градиенту, что будет свидетельствовать о достижении почти стационарной области?

6. Список литературы

1. В.В. Налимов, Н.А. Чернова. Статистические методы планирования экстремальных экспериментов. М.: Наука, 1965. -340с.

2. Адлер Ю.П., Маркова Е.В., Грановский Ю.В. Планирование эксперимента при поиске оптимальных условий. – М.: Наука, 1976. – 280 с.;

7. Варианты индивидуальных заданий

Освоить методику и получить навыки нахождения оптимальных условия функционирования системы методом крутого восхождения по поверхности отклика на модели объекта.