CC BY
161
УДК 519.714.2 DOI 10.21661/Г-467707
С.С. Гусев
АНАЛИЗ МЕТОДОВ И ПОДХОДОВ РЕШЕНИЯ ЗАДАЧ МНОГОКРИТЕРИАЛЬНОГО ВЫБОРА В УСЛОВИЯХ НЕОПРЕДЕЛЕННОСТИ
Аннотация: в статье рассматривается анализ методов и подходов решения задач многокритериального выбора в условиях неопределенности. Приводится описание теории многомерной полезности, теории нечетких множеств Заде и метода анализа иерархий Саати о принятии решений в многомерных многокритериальных системах. Проводится анализ управления ядерными энергетическими установками на подкритическом уровне. В качестве методов управления ресурсами используется Парето-оптимизационный подход управления. Исследуется работа алгоритма оптимизации решений, состоящего из двух этапов.
Ключевые слова: идентификация, динамический объект, многомерные многокритериальные системы, ядерные энергетические установки.
S.S. Gusev
ANALYSIS OF METHODS AND APPROACHES OF SOLVING PROBLEMS WITH MULTI-CRITERIA CHOICE IN CONDITIONS OF UNCERTAINTY
Abstract: the article deals with analysis of methods and approaches solvingproblems with multi-criteria choice in conditions of uncertainty; describes the theory of multidimensional utility theory of fuzzy sets Zadeh and method of analysis of hierarchies Saaty decision-making in multiobjective multidimensional systems. The author has also provided the analysis of the management of nuclear power plants on the sub-critical level. Pareto-optimization approach of control has been used as methods of resource management. It explores the work of the optimization algorithm solutions consisting of two stages.
Keywords: identification, dynamic object, multidimensional multicriteria system, nuclear power plant.
1. Введение
Проектирование сложной технической системы обычно включает в себя этап параметрической оптимизации [1]. В настоящее время при этом чаще всего используется многокритериальная оптимизация, при которой качество функционирования объекта проектирования определяется некоторым набором критериев оптимальности.
Методы решения задачи многокритериальной оптимизации чрезвычайно разнообразны. Существует несколько способов классификации этих методов, например, классификация, основанная на содержании и форме использования дополнительной информации о предпочтениях лица, принимающего решения. В соответствии с этой классификацией выделяются следующие классы методов решения задачи многокритериальной оптимизации:
- методы, основанные на предварительном построении аппроксимации множества Парето;
- априорные методы;
- апостериорные методы;
- адаптивные методы.
Для пользователя наиболее удобны методы, основанные на непосредственном использовании множества Парето (а, тем самым, и фронта Парето). В этом случае лицо, принимающее решение, выбирает компромиссное решение на фронте Парето не формальными методами, исходя только из своих предпочтений, а также из ряда заданных ограничений, ряда факторов, ограничивающих выбор критериев лица, принимающего решения. Основным недостатком методов этого класса, сдерживающим их широкое применение, является их высокая вычислительная сложность. Известно значительное количество методов и алгоритмов приближенного построения множества Парето [2; 3]. Новым и высокоэффективным классом таких методов являются методы на основе генетических алгоритмов.
2 https://interactive-plus.ru
Содержимое доступно по лицензии Creative Commons Attribution 4.0 license (CC-BY 4.0)
Генетические алгоритмы принадлежат к классу эволюционных алгоритмов и обладают рядом характеристик, делающих их более предпочтительными, чем классические методы оптимизации. Так генетические алгоритмы применимы к задачам большой размерности и способны захватить Парето-оптимальные точки даже при однократном запуске алгоритма.
- метод VEGA (Vector Evaluated Genetic Algorithm);
- метод FFGA (Fonseca and Fleming's Multiobjective Genetic Algorithm);
- метод NPGA (Niched Pareto Genetic Algorithm);
- метод SPEA (Strength Pareto Evolutionary Algorithm).
В методе VEGA селекция производится по переключающимся частным критериям оптимальности. Для каждого из m частных критериев создается подпопу-ляция, содержащая Np/m индивидов, где Np - размер всей популяции. Индивиды в j-ю подпопуляцию отбираются с помощью пропорциональной селекции по критерию 0j(X). Далее подпопуляции смешиваются для получения популяции размера Np, а затем по общей схеме осуществляются скрещивание и мутация.
Метод FFGA использует процедуру ранжирования индивидов, основанную на Парето-доминировании. При этом ранг каждого из индивидов определяется количеством доминирующих его других индивидов данной популяции (так что чем ниже ранг, тем индивид ближе к множеству Парето). Пригодность индивида вычисляется на основе величины, обратной его рангу. Для отбора индивидов в следующее поколение используется процедура турнирной селекции.
В методе NPGA, в отличие от методов VEGA, FFGA, существует механизм поддержания разнообразия популяции. Метод основан на формировании попу-ляционных ниш.
Метод SPEA является самым сложным из числа рассматриваемых методов и, также, как метод FFGA, использует селекцию, основанную на Парето-доми-нировании. Для предотвращения преждевременной сходимости, метод использует популяционные ниши. Очень важным свойством метода SPEA является возможность априорного задания количества итоговых точек в искомой аппроксимации множества Парето.
В докладе исследуется работа специального алгоритма идентификации, учитывающего определенную информацию о параметрах объекта. Проводится анализ работы алгоритма принятия решений в условиях неопределенности на основе Парето-Нэш-УКУ-Шепли-комбинаций теории управления структурно-сложными системами. Анализируется связь точности идентификации и величины ошибки измерения.
В качестве методов оптимизации управления ресурсами рассмотрены динамические процессы принятия решений в условиях неопределенности с использованием стабильно-эффективных компромиссов (СТЭК) на основе Парето-Нэш-УКУ-Шепли-комбинаций теории управления структурно-сложными системами, дополненные и расширенные алгоритмами поиска стабильных решений на основе активных равновесий и векторного равновесия по Нэшу. К динамическим задачам принятия решений относятся задачи с неустойчивой во времени исходной информацией, то есть с течением времени изменяются состав и свойства альтернатив, набор критериев выбора, их относительная важность, а также предпочтительность альтернатив, критериев и исходов.
2. Методы принятия решений в задачах многокритериального выбора
Для решения задач многокритериального выбора в условиях неопределенности предложено множество математических методов, классификация которых приведена в [4]. Методы прикладной теории принятия решений различаются способом представления и обработки информации. Подход к проблеме выбора может основываться на отношениях порядка среди альтернатив - классической модели принятия решений, в которой каждой альтернативе ставится в соответствие некоторое число, либо на отношениях включения - поведенческой модели принятия решений, основанной на принадлежности альтернатив к некоторому множеству. Среди методов классического подхода наибольшей универсальностью и теоретической обоснованностью обладают методы теории полезности [5], методы теории нечетких множеств [6; 7; 9] и метод анализа иерархий [8].
Теория многомерной полезности позволяет для задач в условиях риска и неопределенности получить функцию многомерной полезности, максимальное
значение которой соответствует наиболее предпочтительному варианту. Многомерная функция полезности обычно получается как аддитивная или мультипликативная комбинация одномерных функций, которые строятся на основании опроса экспертов и позволяют провести ранжирование возможных исходов без взаимного сравнения альтернатив. При этом делается допущение о взаимной независимости критериев по полезности [4]. Процедура построения функции полезности требует привлечения значительных объемов информации и является достаточно трудоемкой. Достоинством этого подхода является возможность оценки любого количества альтернативных вариантов с использованием полученной функции. В случае неустойчивой исходной информации применение методов теории полезности становится малоэффективным.
Теория нечетких множеств, предложенная Заде, позволяет представить знания о предпочтительности альтернатив по различным критериям с помощью нечетких множеств. Формирование нечетких множеств является более простой и менее трудоемкой процедурой, чем построение функций полезности. Для выявления лучших вариантов по совокупности критериев необходимо иметь в распоряжении информацию о важности критериев и типах возможных отношений между ними. Теория нечетких множеств предоставляет различные средства для учета взаимных отношений критериев: использование весовых коэффициентов, нечеткие отношения предпочтения, нечеткий логический вывод на правилах определения лучшей альтернативы и т. д [6; 7; 9]. Широкие возможности представления знаний и простота вычислительных процедур делают эту теорию очень привлекательным инструментом для создания систем поддержки принятия решений. При этом требуется теоретическое и экспериментальное исследование получаемых системами результатов с целью проверки их адекватности, согласованности, надежности и т. д.
Метод анализа иерархий о принятии решений в многомерных многокритериальных системах (ММС), предложенный Саати, основанный на парных сравнениях альтернативных вариантов по различным критериям и последующим ранжированием набора альтернатив по всем критериям и целям.
Взаимоотношения между критериями учитываются путем построения иерархии критериев и применением парных сравнений для выявления важности критериев и подкритериев. Метод отличается простотой и дает хорошее соответствие интуитивным представлениям. Главным недостатком этого подхода является большое количество требуемой информации об объекте исследования, представляющей собой множество оценок предпочтительности, полученных в процессе парного сравнения альтернатив и критериев. Благодаря простоте метод хорошо подходит для решения динамических задач принятия решений.
3. Ядерная энергетическая установка и основные узлы ее сопряжения
Рассмотрим ядерную энергетическую установку (ЯЭУ), простейшая схема которой приведена на рис. 1 [4].
Рис. 1. Простейшая схема двухконтурной ЯЭУ с паротурбинным циклом
Основными элементами в ЯЭУ являются ядерный реактор (ЯР) и тепловое оборудование. Первый контур включает в себя ЯР, парогенератор (ПГ), циркуляционный насос Н1; трубопроводы горячего (от ЯР) и холодного (к ЯР) теплоносителя.
В теплосиловое оборудование второго контура ЯЭУ входят турбина (Т) с электрическим генератором (ЭГ), конденсатор отработанного пара (К), циркуляционный насос Н2 и т. д.
Рассмотрим подробнее основные элементы конструкции ЯР, изображенного на рис. 2. Основными элементами конструкции ЯР являются: 1 - управляющие
стержни, 2 - боковой отражатель, 3 - теплоноситель, 4 - биологическая защита, 5 - активная зона (АЗ), 6 - замедлитель, 7 - ядерное топливо.
Аппаратура систем контроля, управления и защиты - это комплекс механизмов, приборов, регулирующих устройств, предназначенный для обеспечения безаварийной эксплуатации ядерной установки, то есть исключения самопроизвольного разгона ЯР или отклонений технологических параметров установки от заданных значений.
Основную роль в обеспечении безопасности эксплуатации ЯР призвана сыграть система управления и защиты (СУЗ), на которую возлагаются функции по управлению цепной реакцией при пуске, переходе с одного уровня мощности на другой и остановке ЯР, а также быстрому прекращению реакции деления в случае возникновения аварийной ситуации.
1
Рис. 2. Основные элементы ЯР
Конструкция АЗ и других систем, определяющих условия ее работы, должны исключать превышение пределов повреждения тепловыделяющих элементов (ТВЭЛов), устанавливаемых для нормальной эксплуатации и нарушений нормальной эксплуатации, включая проектные аварии.
Не допускается превышение пределов безопасной эксплуатации ни при одном из следующих нарушений нормальной эксплуатации (с учетом действия защитных систем):
единичные отказы в системах управления реакторной установки;
- потеря энергоснабжения главных циркуляционных насосов (ГЦН);
- отключение турбогенераторов и потребителей тепла;
- потеря всех источников энергоснабжения нормальной эксплуатации;
- течи контура теплоносителя реактора, компенсируемые системами подпитки нормальной эксплуатации;
- непосадка предохранительных клапанов.
Конструкция и компоновка АЗ, реактора и другого оборудования первого контура (включая систему очистки теплоносителя) с учетом работы систем безопасности должны исключать при авариях, в том числе с разрушением или расплавлением топлива, образования объемов продуктов аварии, в которых может возникнуть самоподдерживающаяся цепная реакция.
Текущее состояние ЯР можно охарактеризовать эффективным коэффициентом размножения нейтронов k или реактивностью р, которые связаны следующим соотношением:
к -1 Р = Т.
Для этих величин характерны следующие значения:
- k > 1 - цепная реакция нарастает во времени, реактор находится в надкри-тичном состоянии, его реактивность р > 0;
k < 1 - реакция затухает, реактор подкритичен, р < 0; k = 1 - число делений ядер постоянно, реактор находится в стабильном критическом состоянии, р = 0.
Рабочие органы СУЗ, как наиболее эффективного органа аварийной защиты ЯР должны обладать достаточным быстродействием для перевода АЗ реактора в подкритическое состояние (кэфф < 1).
В случаях, если эффективность аварийной защиты недостаточна для длительного поддержания АЗ в подкритическом состоянии, когда кэфф не должен превышать 0,98, то есть кэфф < 0,98, должно быть предусмотрено автоматическое подключение другой системы остановки реактора, обладающей достаточной
8 https://interactive-plus.ru
Содержимое доступно по лицензии Creative Commons Attribution 4.0 license (CC-BY 4.0)
эффективностью для поддержания АЗ в подкритическом состоянии с учетом возможного высвобождения реактивности.
Рассмотрим СУЗ, как основную систему защиты для вывода ЯР из надкритического состояния (кэфф > 1) в подкритический уровень (кэфф < 1), используя математический аппарат с использованием обратной связи. В качестве образной связи будут рассмотрены компенсирующие стержни, используемые для аварийных средств защиты, воздействие которых приводит к замедлению цепной реакции деления нейтронов и, в крайнем случае, - остановке реактора.
4. Постановка задачи
Рассмотрим две противодействующие и в то же время сдерживающие цепную реакцию деления нейтронов системы ЯР: АЗ и СУЗ.
Активными средствами АЗ будут выступать ядерное топливо (стержни, постоянно находящиеся в АЗ для производства реакции) и теплоноситель, пассивными средствами АЗ будут выступать регулирующие стержни, замедлитель, имитаторы ТВЭЛов, и боковой отражатель - сплошная графитовая кладка толщиной 0,65м. Активными средствами СУЗ будут выступать компенсирующие стержни быстродействующей аварийной защиты (БАЗ), скорость введения которых в водный раствор АЗ в случае аварийной защиты составляет <0,9 секунды с давлением 16 МПа. Пассивными средствами СУЗ будут выступать борный раствор и компенсирующие стержни медленной аварийной защиты (МАЗ) время ввода которых в АЗ для замедления цепной реакции деления и, в крайнем случае, остановки реактора не столь быстро по сравнению с БАЗ и составляет ~18 секунд.
На рис. 3 приведена структура взаимодействия АЗ и СУЗ, где АС - совокупность активных средств коалиции, ПС - совокупность пассивных средств коалиции каждой из систем.
Рис. 3. Структура взаимодействия АЗ и СУЗ
Для упрощения решения задачи будем рассматривать только регулирующие (91), компенсирующие (9) и постоянно находящиеся (51) в АЗ стержни.
Система уравнений задается на основе динамики средних следующим образом:
(1)
Х1 = -Р31 • 4 • хз • К(Х1) х2 =-Рз2 • (1 - 42) • Х3 • К(Х2)
хз = -рз • 4 • х • К(X)
х4 = р4 • (1 - 41) • х1 • Х4 )
где р - эффективность воздействия одного объекта /-ого вида одной системы на один объекту-ого типа другой системы, 0 < р < 1; 4 - доля АС воздействия одной коалиции на АС противодействующей коалиции, 0 < 4 < 1; (1 - 4) - доля АС воздействия одной коалиции на АС противодействующей коалиции; х - текущая средняя численность объектов /-ого типа:
(X,, X, < 1
R( xt) =
v 1, xt > 1
Рассмотрим теоретико-игровую модель кооперативного взаимодействия двух взаимосвязанных и взаимозависимых систем АЗ - СУЗ в пошаговом варианте событий, где шаг - конечный интервал времени и число шагов конечно (к = 1, г). Каждое активное средство делает на шаге один ход. Шаг равен
Т
АТ = Т. г
В пошаговом варианте система (1) имеет следующий вид:
(к +1) = ^ (к) - Р31 • д2 • х3 (к) ■ Я(х ) х2(к +1) = х2(к) - РЪ2 ■ (1-д2) ■ хъ(к )■ Я( х2) х3 (к +1) = х3 (к) - Р13 ■ д ■ х1 (к) ■ Я(х3) х4 (к +1) = х4 (к) - Р14 ■ (1 - «71) ■ х (к) ■ Я(х4 )
где к = 1,2,..., 0 < д, < 1 (г = 1,2), 0 < р < 1 (г = 1,4; 7 = 1,2,4,4), х, > 0, х,(к) - численность к началу £-ого шага.
В качестве показателя терминальных потерь (!) выберем показатель, имеющий смысл суммарного перевеса по АС и ПС составляющим.
т
¿АЗ =аи ■ [х42 (Т) - х? (Т)]+ а12 ■ [х42 (Т) - х22 (Т)]+ а14 ■ |х4 ■ &
<0
T
JCY3 = а21 • [x2 (T) - x2 (T)]+ а22 • [x2 (T) - x2A (T)]+ a23 • J Xj ' d<
0
где JA3 - показатель потерь коалиции АЗ; Jcy3 - показатель потерь коалиции
СУЗ; ау - весовые коэффициенты, определяющие целевой приоритет каждой из
сторон в АС и ПС противоположных сторон (терминальная составляющая) или в увеличении интегральной скорости убывания АС противоположной стороны (интегральная составляющая) (0 <ау < 1, ал +аг2 +аг3 = 1, / = {1,2}), где значения
коэффициентов задаются в зависимости от тактики каждой из сторон.
Используя алгоритм общего вида, базирующийся на определении двух противодействующих и в то же время взаимосвязанных и взаимозависимых сторон -угроз и контругроз (УКУ), сформируем алгоритм получения сетевых приближений для задач данного класса: двухкоалиционных, двухкритериальных.
Сетевой алгоритм на первом этапе, когда формируется двумерная ортогональная равномерная сеть, позволяет получить начальное приближение для поиска точного решения на втором этапе оптимизации, когда формируется множество оптимальных сетевых решений, которые можно использовать в качестве начальных приближений для получения оптимального управления ММС. Даны исходные экспериментальные данные в виде таблицы 1.
Таблица 1
Временные замеры работы сетевого алгоритма
Глубина проникновения БАЗ в АЗ, мм 1350 — — — — —
Глубина проникновения МАЗ в АЗ, мм 270 540 810 1000 1200 1350
Время, с 1 4 8 12 16 18
Реализация сетевого алгоритма оптимизации осуществлена на алгоритмическом языке с использованием начальных условий.
5. Алгоритм идентификации Рассмотрим задачу взаимодействия двух противодействующих объектов двухкоалиционной динамической системы. Цель каждой из системы состоит в преодолении и получении приоритета над противодействующей системой с целью сохранения баланса всего комплекса в целом, то есть удержания цепной реакции деления нейтронов на заданном уровне, в данном случае подкритическом уровне ЯР при кэфф < 1.
Рассмотрим двухэтапный алгоритм получения УКУ-Шепли-оптимального управления прогнозом динамики конфликта [1] АЗ - СУЗ при следующих нулевых начальных условиях на первом такте. Исходные данные заданы условно, поэтому полученные оценки задачи не стоит рассматривать как результат решения исследовательской работы с практическим уклоном. Цель задачи состоит в определении Парето-УКУ-Шепли-оптимального вектора параметров при выведении элементов системы СУЗ из АЗ ЯР и ввода регулирующих стержней в АЗ ЯР для начала производства цепной реакции деления.
X (0) = 51 *2(0) = 91
X (0) = 0 ' Х4(0) = 0
P31 (0) = 0,99 P,2(0) = 0,99 P33 (0) = 0,99, Рзз(0) = 0,99
« i(0) = 0,9
«12(0) = 0,1
«21(0) = 0,2, «22(0) = 0,8
\ JА3 =«11 • k2 (T) - x2 (T)]+«!2 • X2 (T) - *22 (T)]
UСУ3 = «21 • к (T) - X2 (T)] + «12 • [x2 (T) - X4 (T)]
^ min ^ min
На первом такте прогноза происходит определение значений вектора дележа Шепли на основе Парето-Нэш-оптимизации и векторная оптимизация управления коалициями для минимизации отклонения от точки Шепли. Результатом оптимизации является точка сильной предпосылки игры (СТЭК при одно-тактовом прогнозе).
<
Результатом первого такта прогноза является определение Парето-УКУ-Шепли-оптимального вектора параметров:
^ = 0,496, дор1 = 0,573}.
Вектор состояния к концу первого такта:
хх = 51, х2 = 76, х3 = 2, х4 = 4,7 .
Начальные условия второго такта:
х1(1) = 51 x2(1) = 76 x3 (1) = 2 ' x4(1) = 4,7
P 1(1) = 0,99 P 2(1) = 0,99 P 3 (1) = 0,99, P 4(1) = 0,99
«„(1) = 0,9
«12(1) = 0,1
«21(1) = 0,2, «22(1) = 0,8
^ min
13З =«11 • [х2 (Т) - х2 (Т)] + «12 • [х? (Т) - х2 (Т)]
К СУЗ = «21 • [х12 (Т) - х32 (Т)] + «12 • [х22 (Т) - х42 (Т)] ^ т1П '
На втором такте прогноза происходит повторное определение значений вектора дележа Шепли на основе Парето-Нэш-оптимизации и векторная оптимизация управления коалициями для минимизации отклонения от точки Шепли. Результатом оптимизации является точка сильной предпосылки игры (СТЭК на двух тактах прогноза).
Парето-УКУ-Шепли-оптимальный вектор параметров на втором такте:
= 0,519, = 0,549}.
Вектор состояния после двух тактов прогноза:
х1 — 51, х^ — 59, х3 — 2, х4 — 7 .
На рис. 4 приведены изменения численностей и Парето-УКУ-Шепли-опти-мального параметризованного управления.
*i(0) = 51 xi 0) = 51 *i(2) = 51
= 91 = 76 = 59
*3(0) = 0 = 2 z3(2) = 2
*4<0) = 0 = 4,7 = 7
qf = 0,496, qf = 0,573 qf = 0,519, qf = 0,549
0 77 2 T
Рис. 4. Изменения численностей
<
<
На двух тактах формируется оптимальное управление активной и пассивной систем АЗ и СУЗ - , д°/.
На первом такте прогноза происходит изменение численностей АС и ПС СУЗ (х3 = 2, х4 = 4,7) и ПС АЗ (х2 = 76), тем самым, сохраняя баланс системы для
вывода регулирующих стержней из АЗ и сброса компенсирующих стержней СУЗ БАЗ и МАЗ. Оптимальный вектор параметров для обеих систем становится равным: = 0,496, д°р = 0,573 .
На втором такте численность регулирующих стержней АЗ и компенсирующих стержней СУЗ изменяются до значений: х2 = 59, х3 = 2, х4 = 7, тем самым, продолжая сохранять баланс системы на подкритическом уровне, а оптимальный вектор параметров для обеих систем становится равным: д°р = 0,519, д°°р = 0,549. Результаты в конце второго такта соответствуют точке сильной предпосылки игры.
На рис. 5 показано окончательное Парето-УКУ-Шепли-оптимальное параметризованное управление АЗ и СУЗ на двух тактах прогноза.
Рис. 5. Парето-УКУ-Шепли-оптимальное параметризованное управление АЗ
и СУЗ на двух тактах
Применение полученных программных управлений обеспечивает Парето-Нэш-УКУ-Шепли-компромисс с двухтактным прогнозом. В общем случае повторение данной процедуры на следующем временном интервале (т, 2Т) с измеренным вектором х(Т) позволяет получить программно-корректируемый закон управления (ПКЗУ).
6. Заключение
В докладе рассмотрен двухэтапный алгоритм получения УКУ-Шепли-опти-мального управления прогнозом динамики взаимодействия двух противодействующих систем: АЗ - СУЗ в условиях неопределенности на основе Парето-ком-бинаций теории управления структурно-сложными системами. При использовании рассмотренного алгоритма возможно точное определение неизвестных параметров динамического объекта при условиях, которые найдены в докладе. Практическое применение полученных программных модулей управления многомерными многокритериальными системами обеспечивает возможность получения Парето-компромисс с двухтактным прогнозом.
Список литературы
1. Норенков И.П. Основы автоматизированного проектирования: Учеб. для вузов. - М.: Изд-во МГТУ им. Н.Э. Баумана, 2006. - 448 с.
2. Соболь И.М. Выбор оптимальных параметров в задачах со многими критериями / И.М. Соболь, Р.Б. Статников. - М.: Дрофа, 2006. - 175 с.
3. Подиновский В.В. Парето-оптимальные решения многокритериальных задач / В.В. Подиновский, В.Д. Ногин. - М.: ФИЗМАТЛИТ, 2007. - 256 с.
4. Методы классической и современной теории автоматического управления: Учебник: В 5-и т. Т. 1: Теория оптимизации систем автоматического управления / Под ред. К.А. Пупкова и Н.Д. Егупова. - 2-е изд., перераб. и доп. - М.: Изд-во МГТУ им. Н.Э. Баумана, 2004. - 744 с.
5. Борисов А.Н. Диалоговые системы принятия решении на базе мини-ЭВМ. Информационное, математическое и программное обеспечение / А.Н. Борисов, Э.Р. Вилюмс, Л.Я. Сукур. - Рига: Зинатне, 1986.
6. Кипи Р.Л. Принятие решений при многих критериях: замещения и предпочтения / Р.Л. Кипи, X. Райфа. - М.: Радио и связь, 1981.
7. Нечеткие множества в моделях управления и искусственного интеллекта. - М.: Наука, 1986.
8. Алексеев А.В. Интеллектуальные системы принятия проектных решений / А.В. Алексеев, А.Н. Борисов, Э.Р. Вилюмс [и др.]. - Рига: Зинатне, 1997.
9. Саати Т.Л. Принятие решений. Метод анализа иерархий. - М.: Радио и связь, 1993.
References
1. Norenkov, I.P. (2006). Osnovy avtomatizirovannogo proektirovaniia. M.: Izd-vo MGTU im. N.E. Baumana.
2. Sobol', I.M., & Statnikov, R.B. (2006). Vybor optimal'nykh parametrov v zadachakh so mnogimi kriteriiami. M.: Drofa.
3. Podinovskii, V.V., & Nogin, V.D. (2007). Pareto-optimal'nye resheniia mnog-okriterial'nykh zadach. M.: FIZMATLIT.
4. (2004). Metody klassicheskoi i sovremennoi teorii avtomaticheskogo uprav-leniia. M.: Izd-vo MGTU im. N.E. Baumana.
5. Borisov, A.N., Viliums, E.R., & Sukur, L.Ia. (1986). Dialogovye sistemy prin-iatiia reshenii na baze mini-EVM. Informatsionnoe, matematicheskoe i programmnoe obespechenie. Riga: Zinatne.
6. Kipi, R.L., & Raifa, X. (1981). Priniatie reshenii pri mnogikh kriteriiakh: za-meshcheniia i predpochteniia. M.: Radio i sviaz'.
7. (1986). Nechetkie mnozhestva v modeliakh upravleniia i iskusstvennogo intel-lekta. M.: Nauka.
8. Alekseev, A.V., Borisov, A.N., & Viliums, E.R. (1997). Intellektual'nye sistemy priniatiia proektnykh reshenii. Riga: Zinatne.
9. Saati, T.L. (1993). Priniatie reshenii. Metod analiza ierarkhii. M.: Radio i sviaz'.
Гусев Сергей Сергеевич - соискатель ФГБУН «Институт проблем управления им. В. А. Трапезникова РАН», Россия, Москва.
Gusev Sergey Sergeevich - Ph.D. student at V.A. Trapeznikov Institute of Control Sciences of Russia Academia of Sciences, Russia, Moscow.
