УДК 519.6

 

Карпенко А.П.
(д.ф.-м.наук, профессор МГТУ им.Н.Э.Баумана,
телефон: 263-65-26, E-mail: karpenko@nog.ru)

Федорук В.Г.
(к.т.н., доцент МГТУ им.Н.Э.Баумана,
телефон: 263-65-26, E-mail: fedoruk@comcor.ru)

 

 

Работа является четвертой в серии работ, посвященных адаптивным методам решения непрерывной конечномерной задачи многокритериальной оптимизации. В работах [1-3] рассмотрены методы аппроксимации функции предпочтений лица, принимающего решения (ЛПР), на основе планов первого и второго порядков, а также нейронных сетей и нечеткой логики. Данная работа посвящена рассмотрению адаптивных методов решения задачи многокритериальной оптимизации, построенных на базе указанных методов аппроксимации.

 

1.                Основные обозначения и общая схема методов

            Постановка задачи многокритериальной оптимизации приведена в работе [1]. Повторим основные обозначения, введенные там.

            Вектор  - n-мерный вектор параметров задачи (вектор варьируемых параметров), где

 -

замкнутое множество допустимых значений вектора параметров,  - ограничивающие функции.

             - векторный критерий оптимальности, определенный на множестве . ЛПР стремится минимизировать каждый из частных критериев .

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

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

,                                            (1)

где  - результат операции свертки (суперкритерий),

-

вектор весовых множителей,  - множество допустимых значений этого вектора.

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

.

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

.                                                              (2)

            Величина  полагается лингвистической переменной со значениями ”Очень-очень плохо”, ”Очень плохо”, …, ”Отлично”. Центры соответствующих нечетких множеств обозначаются  и имеют значения 1, 2, …, 9, соответственно.

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

.                                                          (3)

            Общая схема предлагаемых методов состоит из следующих основных этапов:

1)      задание ЛПР начальной точки , ;

2)      выполнение на основе планов первого порядка серии экстремальных экспериментов по максимизации функции предпочтений ;

3)      построение по результатам указанных экспериментов первой аппроксимирующей функции ;

4)      формирование области планирования для аппроксимации функции ;

5)      построение второй аппроксимирующей функции ;

6)      определение начальной точки  для следующей итерации.

            На этапе 2 используются симплекс-планы на основе правильных или неправильных симплексов, а также регрессионные планы первого порядка [1]; на этапе 3 – аппроксимация кусочно-линейной функцией [1], персептронной нейронной сетью и нейронной сетью с радиальными базисными функциями [3], а также аппроксимация на основе нечетких множеств [3]; на этапах 5, 6 – линейная и квадратичная аппроксимация [2].

 

2.                Этап 2 – экстремальный эксперимент на основе планов первого порядка

            В данном разделе рассматриваются экстремальные эксперименты по максимизации функции  на основе симплекс-планов, а также регрессионных планов первого порядка [1]. В подразделе 2.1 рассматривается экстремальный эксперимент на основе правильных симплексов, в подразделе 2.2 – на основе неправильных симплексов, в подразделе 2.3 – экстремальный эксперимент на основе регрессионных планов первого порядка и градиентного метода оптимизации с дроблением шага, в подразделе 2.4 – такой же эксперимент, но основе градиентного метода крутого восхождения [4].

 

            2.1 Правильные симплекс-планы (метод М2.1). Отметим прежде, что по количеству точек в спектре симплекс-планы являются насыщенными планам, т.е. требуют минимально возможного количества испытаний, однако имеют низкую помехозащищенность (сглаживание отсутствует) [4].

            Экстремальный эксперимент выполняется по следующей схеме (здесь и далее приводятся упрощенные схемы) [1].

1)                 Исходя из текущей начальной точки , СМКО (Система МногоКритериальной Оптимизации - Multi Criteria Decision Making System) строит начальный правильный симплекс , где . Если одна или несколько вершин симплекса оказываются вне области допустимых значений , то СМКО выдает ЛПР соответствующее сообщение и предлагает либо изменить начальную точку , либо уменьшить текущую длину ребра симплекса .

2)                 Для каждой из вершин ,  симплекса  СМКО выполняет следующие действия.

2.1)           Решает задачу (1) с весами .

2.2)           Полученные значения компонентов вектора параметров , а также соответствующие значения частных критериев оптимальности  предъявляет ЛПР.

2.3)           Запрашивает у ЛПР оценку предъявленного решения задачи многокритериальной оптимизации - требует ввести соответствующее значение лингвистической переменной .

3)                 СМКО находит минимальное и максимальное значения соответствующих четких значений лингвистической переменной

,      .

 

4)                 По приведенной ниже схеме СМКО отражает вершину  симплекса  относительно центра тяжести оставшихся вершин  - находит вершину  и новый симплекс .

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

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

4.3)           Если отражение всех указанных вершин не привело к успеху, то аналогично делается попытка отразить вершины, в которых величина  имеет следующее за  значение, и т.д. до вершин, в которых величина  имеет значение .

4.4)           Если имеет место вращение симплекса вокруг некоторой его вершины, то СМКО информирует об этом пользователя и предлагает (а) закончить вычисления, (б) уменьшить текущую длину ребра симплекса  и/или начать новый поиск из новой начальной точки.

4.5)           Если ЛПР выбирает вариант (а), то СМКО формирует для ЛПР соответствующее сообщение и завершает вычисления. Если выбирается вариант (б), то СМКО полагает  и выполняет переход к п. 1 (здесь  - свободный параметр метода).

5)                 Для вершины  СМКО выполняет действия, аналогичные действиям, указанным в п.2 (в результате становится известной величина ), полагает  и переходит к п.3.

            Заметим, что в этой схеме на всех итерациях, кроме первой, от ЛПР требуется оценка только одного решения.

 

            2.2 Неправильные симплекс-планы (метод М2.2). Эксперимент выполняется по упрощенной схеме метода Нелдера-Мида [5]. Упрощение заключается в исключении из числа операций над симплексом операции редукции (поскольку эта операция требует от ЛПР оценки значений функции предпочтений в  точках, в отличие от одной точки для операций отражения, сжатия и растяжения симплекса). Эксперимент выполняется по следующей схеме [1].

1)                 Исходя из текущей начальной точки ,  СМКО строит начальный правильный симплекс  (см. подраздел 2.1).

2)                 Для каждой из вершин  симплекса  СМКО решает задачу (1) и получает от ЛПР соответствующее значение лингвистической переменной  (см. подраздел 2.1).

3)                 СМКО находит минимальное, максимальное, а также первое, предшествующее максимальному, значения соответствующих четких значений

,     ,    .

 

4)                 СМКО отражает вершину  симплекса  относительно центра тяжести оставшихся вершин  - находит вершину  и новый симплекс . СМКО запрашивает у ЛПР оценку  решения, соответствующего точке  (см. подраздел 2.1). Если  и , то выполняется переход к п.5. Если , но , то выполняется переход к п.4, а если  - к п.6.

5)                 СМКО выполняет растяжение симплекса  в направлении  - находит новый симплекс  и запрашивает у ЛПР оценку . Если , то выполняется переход к операции отражения симплекса  (п.4), иначе – к операции отражения симплекса  (п.4).

6)                 СМКО выполняет сжатие симплекса  в направлении  - находит новый симплекс . МСКО запрашивает у ЛПР оценку  и выполняет переход к операции отражения симплекса .

 

            2.3. Регрессионные планы первого порядка и градиентный метод с дроблением шага (метод М2.3). Для исходных («натуральных») факторов  искомая линейная регрессионная модель имеет вид

,

где  - неизвестные коэффициенты регрессии [4]. Переход от исходных факторов  к нормированным (стандартным) факторам  выполняется по формуле , , , где  - основной (базовый) уровень фактора , а  - шаг варьирования этого фактора. Для нормированных факторов искомая регрессионная модель имеет вид

,

где , , .

            Экстремальный эксперимент в данном случае строится на основе слабо не насыщенных регрессионных планов первого порядка для количества критериев . Планы объединены в библиотеку планов , основные параметры которой представлены в Табл. 1. Здесь M – требуемое количество опытов,  - число степеней свободы, ПФЭ – полный факторный эксперимент, ДФЭ – дробный факторный эксперимент. Как это принято при планировании экспериментов, планы в библиотеки даны не для «натуральных», а для нормированных факторов [4].

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

 

Таблица 1. Основные параметры библиотечных планов первого порядка

 

            Схема эксперимента на основе регрессионных планов первого порядка и градиентного метода с дроблением шага имеет следующий вид.

1)                 Исходя из текущей начальной точки , СМКО строит область планирования в виде параллелепипеда , ). Здесь  - шаг варьирования фактора  на r-ой итерации.

2)                 СМКО выбирает из библиотеки планов соответствующий план первого порядка , , выполняет переход от стандартных факторов  к «натуральным» факторам ,  и проверяет принадлежность каждой из полученных точек области допустимых значений . Если одна или несколько точек плана ,  оказываются вне области допустимых значений , то СМКО выдает ЛПР соответствующее сообщение и предлагает либо изменить начальную точку , либо уменьшить величины .

3)                 Для каждой из точек плана ,  СМКО выполняет следующие действия.

3.1)           Решает задачу (1) с весами . Полученные значения компонентов вектора параметров , а также соответствующие значения частных критериев оптимальности предъявляет ЛПР.

3.2)                               Запрашивает у ЛПР оценку каждого из предъявленных решений - требует ввести соответствующее значение лингвистической переменной .

4)                 Для каждой из точек плана ,  СМКО формирует  значений базисных функций

,

в результате чего становится известной -матрица

.

5)                 СМКО решает систему линейных алгебраических уравнений (СЛАУ)

,                                                    (4)

где симметричная  матрица  есть информационная матрица (Фишера),  - вектор соответствующих четких значений лингвистической переменной ,  - вектор оценок коэффициентов регрессии для стандартных факторов. Поскольку для всех библиотечных планов матрица Фишера  является диагональной с диагональными элементами, равными количеству опытов , решение СЛАУ (4) не представляет труда.

6)                 СМКО выполняет следующие действия [4]:

6.1)                               Осуществляет проверку значимости оценок коэффициентов регрессии ;

6.2)                               Вычисляет коэффициент детерминации  и если , где , то выполняет переход к п. 7. Если , то система предъявляет ЛПР величину  и требует его разрешения продолжить итерации. При ответе «Да» выполняется переход к следующему пункту, а при ответе «Нет» - вычисления заканчиваются.

7)                 СМКО вычисляет градиент построенной функции регрессии и делает шаг в его направлении по следующей схеме:

7.1)                               Определяет компоненты  вектора  для натуральных факторов по формуле , .

7.2)                               Вычисляет текущую величину шага  в направлении градиента, где  - свободный параметр метода, а величина  выбирается из условия =.

7.3)                               Вычисляет координаты нового центра серии опытов

, где  -

нормированные компоненты вектора градиента для натуральных факторов.

7.4)                               Проверяет допустимость точки . Если точка недопустима, то а) уменьшает длину текущего шага , т.е. полагает , где  - свободный параметр алгоритма; б) выполняет переход к п. 7.3. Если точка допустима, то полагает  и переходит к п.1.

 

            2.4. Регрессионные планы первого порядка и градиентный метод крутого восхождения (метод М2.4). Схема эксперимента в этом случае аналогична схеме, рассмотренной в подразделе 2.3. Отличия имеются лишь в пункте 7, который рассматривается ниже.

7)                 СМКО вычисляет градиент построенной функции регрессии и делает шаг в его направлении по следующей схеме.

………….

7.5)                               Проверяет допустимость точки . Если точка недопустима, то уменьшает длину текущего шага  и выполняет переход к п. 7.3. Если точка допустима, то выполняет следующие действия.

7.5.1)     Решает задачу (1) с весами . Полученный в результате вектор параметров , а также соответствующие значения частных критериев оптимальности предъявляет ЛПР. Запрашивает у ЛПР оценку предъявленного решения - требует ввести соответствующее значение лингвистической переменной .

7.5.2)     Если , то на основе диалога с ЛПР либо завершает вычисления, либо уменьшает шаги варьирования факторов , , полагает  и переходит к п.1. Ситуация интерпретируется как достижение максимума функции предпочтений в направлении градиента.

7.5.3)     Если , то полагает, , и переходит к п.7.3. Ситуация интерпретируется следующим образом: в направлении градиента максимум функции предпочтений не достигнут и можно сделать еще один шаг в том же направлении.

 

3.                Этап 3 - построение вспомогательной аппроксимирующей функции

            В данном разделе рассматривается построение первой (вспомогательной) аппроксимирующей функции  на основе результатов серии экстремальных экспериментов, рассмотренных в предыдущем разделе. В подразделе 3.1 рассматривается кусочно-линейная аппроксимация на основе правильных симплекс-планов [1]; в подразделе 3.2- на основе неправильных симплекс-планов; в подразделе 3.3 - кусочно-линейная аппроксимация на основе регрессионных планов первого порядка [1]; в подразделе 3.4 - аппроксимация многослойной персептронной нейронной сетью (MLP-нейронной сетью) [3]; в подразделе 3.5 - аппроксимация нейронной сетью на основе радиальных базисных функций (RBF-нейронной сетью) [3]; в подразделе 3.6 - аппроксимация на основе теории нечетких множеств [3].

 

            3.1. Кусочно-линейная аппроксимация на основе правильных симплекс-планов (метод М3.1) [1]. Пусть по схеме, приведенной в подразделе 2.1, последовательно построены правильные симплексы , ,…,. Тогда построение кусочно-линейной аппроксимации  функции предпочтений ЛПР выполняется по следующей схеме (см. Рис. 1).

            1) СМКО строит линейные функции , ,…,, имеющие в вершинах соответствующих симплексов значения , ,…, , , а также гиперплоскости , ,…,, проходящие через общие точки смежных симплексов.

            2) СМКО строит кусочно-линейную аппроксимирующую функцию

                                                 (5)

где

,

,

…

.

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

            Рис. 1. К кусочно-линейной аппроксимации функции предпочтений ЛПР на симплексах: ; .

 

            Заметим, что приведенная схема аппроксимации учитывает тот факт, что с ростом количества итераций ЛПР обычно оценивает свою функцию предпочтений более точно.

            Поясним первый пункт приведенной схемы. Рассмотрим для примера построение линейной функции , где , ,  - искомые коэффициенты. Из условия прохождения функции  через точки , имеем следующую СЛАУ для определения этих коэффициентов:

                                 (6)

Таким образом, построение функций , ,…,   сводится к решению СЛАУ вида (6). Поскольку все точки ,  различны, матрица этой системы не вырождена и решение существует и единственно.

            Рассмотрим также для примера построение гиперплоскости . Положим, что симплекс  получен в результате отражения вершины  симплекса . Тогда общими для симплексов ,  являются m вершин , ,  и линейную функцию  можно записать в виде

,

где ,  - искомые коэффициенты. Здесь аналогично предыдущему коэффициенты  определяются из СЛАУ

Матрица этой системы не вырождена и ее решение существует и единственно.

 

            3.2. Кусочно-линейная аппроксимация на основе неправильных симплексов (метод М3.2) [1]. Аппроксимация на основе неправильных симплекс-планов выполняется по схеме п. 3. 1.

 

            3.3. Кусочно-линейная аппроксимация на основе регрессионных планов первого порядка (метод М3.3) [1]. Пусть по схеме, приведенной в подразделе 2.3, последовательно построены области планирования в виде параллелепипедов , ,…,. Построение кусочно-линейной аппроксимирующей функции  в данном случае выполняется по следующей схеме (см. Рис. 2).

            1) СМКО строит гиперплоскости , ,…, такие, что в гиперплоскости ,  лежит та из граней параллелепипеда, которую пересекает отрезок прямой .

            2) В качестве искомой кусочно-линейной функции  принимается функция (5), где , ,…,   - соответствующие линейные регрессионные модели,

,

,

…

.

            Как и в подразделе 3.1, приведенная схема аппроксимации учитывает тот факт, что с ростом количества итераций ЛПР обычно оценивает свою функцию предпочтений более точно.

            Поясним первый пункт приведенной схемы. Рассмотрим для примера построение гиперплоскости .     

            а) Пронумеруем вершины параллелепипеда =, ) с помощью бинарного отраженного кода Грея. В результате любые две соседние вершины параллелепипеда  получат номера, двоичное представление которых имеет только одну различающуюся позицию [6].

            б) Пусть  - ближайшая к точке  вершина параллелепипеда  (найти ее легко путем перебора). Используя отмеченное выше свойство кода Грея, найдем все  вершин , ,…,  параллелепипеда , которые являются соседями вершины .

 

            Рис. 2. К кусочно-линейной аппроксимации функции предпочтений ЛПР на параллелепипедах: ; .

 

            в) Вершины , , ,…,  объединим в  различных совокупностей ,  …., ,  по  вершин в каждой так, чтобы в каждую из этих совокупностей входила вершина .

            г) Из условия прохождения гиперплоскости  через все вершины совокупности ,  получим следующую СЛАУ для определения неизвестных коэффициентов , :

Матрица этой системы не вырождена и ее решение существует и единственно.

            д) Найдем уравнение прямой, проходящей через точки , :

.                                               (7)

            е) Из СЛАУ

найдем точку  пересечения прямой (7) с гиперплоскостью , .

            ж) Определим минимальное из расстояний от точки  до точки :

, .

            з) В качестве искомой гиперплоскости  примем гиперплоскость .

 

            3.4. Аппроксимация MLP-нейронной сетью (метод М3.4) [3]. Рассмотрим двухслойную MLP-сеть с сигмоидальными передаточными функциями, имеющую  входных нейронов,  нейронов в промежуточном слое и один выходной нейрон (см. Рис. 3). Обозначим выходы нейронов скрытого слоя . Тогда при заданных весах , , ,  и некоторых фиксированных  выход этой сети (значение искомой аппроксимирующей функции ) определяется как суперпозиция функций

,      , ,

где ,  - передаточные функции нейронов , , соответственно.

 

 

Рис. 3. Топология нейронной сети.

 

            Пусть , ,  - совокупность различных точек множества , в которых ЛПР определены значения его функции предпочтений . При использовании симплекс-планов – это те из вершин симплексов , ,…,, формирование которых рассмотрено в подразделах 2.1, 2.2 и на основе которых в подразделе 3.1 построена аппроксимирующая функция . При использовании регрессионных планов первого порядка – это совокупность точек всех использованных планов (см. подразделы 2.3, 2.4). Заметим, что в случае использования симплекс-планов требование  означает, что на этапе 2 должно быть построено минимум два симплекса , , т.е. .

            Построение аппроксимирующей функции  с помощью MLP-нейронной сети выполняется по следующей схеме.

1)                 СМКО строит двухслойную MLP-сеть с сигмоидальными передаточными функциями, имеющую  входных нейронов,  нейронов в промежуточном слое и один выходной нейрон (см. Рис. 3). Начальное количество нейронов  определяется из равенства

,

где величина  находится из условия

.

2)                 На обучающей выборке , СМКО выполняет параметрический синтез нейронной сети – определяет веса , , , . При этом в качестве оценки ошибки сети  при обработке вектора  используется величина , а в качестве оценки ошибки при обработке всех N векторов из обучающей выборки - величина . Здесь  - весовой множитель i-го примера из обучающей выборки, назначаемый по правилу экспоненциального возрастания весов опытов [3].

3)                 Если заданная точность аппроксимации не достигнута, то методом динамического наращивания узлов [3] СМКО увеличивает количество узлов в скрытом слое и возвращается к п. 2.

           

            3.5. Аппроксимация RBF-нейронной сетью (метод М3.5) [3]. Рассмотрим двухслойную RBF-сеть с сигмоидальными передаточными функциями в скрытом слое и линейной передаточной функцией  в выходном слое (см. Рис. 3). Веса , ,  нейронов скрытого слоя положим равными единице, так что . Тогда для некоторого входного вектора  при заданных весах ,  выход этой сети, т.е. значение искомой аппроксимирующей функции, равен

.

            Аналогично п.3.4, в качестве обучающей выборки используем выборку  , где .

            Построение аппроксимирующей функции  с помощью RBF-нейронной сети выполняется по следующей схеме.

1)                 Из выборки   СМКО случайным образом формирует обучающие выборки  . Здесь , ; , ;  - символ ближайшего целого большего.

2)                 СМКО строит двухслойную RBF-сеть, в которой каждому вектору  из первой обучающей выборки ставится в соответствие нейрон скрытого слоя с центром в точке , , . Таким образом, передаточная функция , , где  - параметр, определяющий ширину колокола функции .

3)                 СМКО повторяет некоторое фиксированное количество раз  следующие действия.

3.1)           Присваивает величинам ,  некоторые положительные случайные начальные значения.

3.2)           Путем решения СЛАУ

полученной на основе первой обучающей выборки, определяет веса  выходного нейрона.

3.3)           С использованием второй обучающей выборки решает многомерную, вообще говоря, многоэкстремальную, задачу условной оптимизации

,

где функция , а .

4)                 На основе результатов предыдущего пункта СМКО выбирает решение, обеспечивающее минимальное значения ошибки , а также соответствующие значения весов .

 

            3.6 Аппроксимация на основе нечетких множеств (метод М3.6) [3]. Пусть по схеме п. 3.5 построены обучающие выборки   , . Положим, что в рамках первой обучающей выборки  раз переменная  приняла значение ,  раз - значение  и т.д. до  и . Соответствующие входные векторы  обозначим , , . Термы  и термы  будем полагать нечеткими множествами, заданными на универсальном множестве :

, ;

, .

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

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

где  - весовые множители, определяемые модифицированным правилом экспоненциального возрастания весов опытов - . Здесь  - свободный параметр алгоритма, величина  находится из условия , I – номер опыта по порядку [3],  - малая положительная величина.

            Обозначим ,b полуширины функций принадлежности , , соответственно. Область допустимых значений параметров  определим как =; ; , где  - заданные константы.

            Построение аппроксимирующей функции  с помощью нечетких множеств выполняется по следующей схеме.

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

, .

2)                 СМКО  осуществляет параметрическую оптимизацию построенной нечеткой модели – повторяет некоторое фиксированное количество раз  следующие действия.

2.1) Присваивает величинам  некоторые допустимые случайные начальные значения.

2.2) Исходя из этих значений, решает на второй обучающей выборке  задачу многомерной, вообще говоря, многоэкстремальной условной оптимизации

,

где критерий оптимальности  имеет вид

, .

Здесь значения аппроксимирующей функции  определяются по следующей схеме.

2.2.1) Для данного вектора  определяется значение многомерной функции принадлежности .

2.2.2) Функцию принадлежности  «срезается» на уровне , .

2.2.3) Полученные нечеткие множества агрегируются – строится нечеткое множество

с функцией принадлежности .

2.2.4) Определяется четкое значение функции , соответствующее входному вектору  - т.е. выполняет операцию дефаззификации нечеткого множества , например, по методу центра тяжести.

3)                 На основе результатов предыдущего пункта СМКО выбирает значения параметров , обеспечивающие минимальное значение ошибки .

            После построения по рассмотренной схеме аппроксимирующей функции  ее значение в некоторой точке  вычисляется по схеме пп. 2.2.1 – 2.2.4.

 

4.     Этап 4 - формирование области планирования

            При использовании на этапе 2 симплекс-планов, симплексы , ,…, могут сформировать область, которая сильно вытянута и произвольно ориентированна в пространстве  (см. Рис. 4). Аналогично, при использовании регрессионных планов первого порядка, параллелепипеды , ,…, могут сформировать такую же область. Поскольку качество аппроксимации существенно зависит от выбора области планирования [4], возникает задача определения этой области так, чтобы она была в некотором смысле оптимальной.

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

            Заметим, что в общей постановке задача построения области планирования является задачей построения выпуклой оболочки , где ,  - вершины симплексов , ,…, или параллелепипедов , ,…, [7,8].

 

Рис. 4. К построению области планирования: ;  ;  - прямоугольная область планирования

 

            Наряду с системой координат  рассмотрим систему координат , оси которой параллельны ребрам искомого параллелепипеда . Матрицу перехода от первой системы координат ко второй системе обозначим , где  - косинус угла между ортами ,  указанных систем координат (здесь учтено, что множество допустимых значений  лежит в первом координатном угле системы координат ).

            В результате задача отыскания оптимального параллелепипеда  формализуется в виде многомерной задачи условной оптимизации