А
нализ и синтез систем управления
УДК 681.5.015.63
ОПТИМАЛЬНЫЙ ПАРАМЕТРИЧЕСКИЙ СИНТЕЗ ПО КРИТЕРИЮ ЗАПАСА РАБОТОСПОСОБНОСТИ1
О.В. Абрамов, Я.В. Катуева, Д.А. Назаров
Институт автоматики и процессов управления ДВО РАН, г. Владивосток
Рассмотрена проблема синтеза аналоговых технических систем с учетом отклонений параметров от расчетных значений. Показано, что при отсутствии информации о закономерностях и количественных характеристиках этих отклонений наилучшее решение задачи параметрического синтеза может быть получено, если в качестве целевой функции принять «запас работоспособности». Обсуждены параллельные алгоритмы решения возникающих при этом задач многовариантного анализа и оптимизации.
ВВЕДЕНИЕ
Основные трудности, возникающие при решении задачи синтеза технических устройств и систем с учетом отклонений параметров от номинальных (расчетных) значений, связаны, прежде всего, с дефицитом информации о закономерностях случайных процессов вариации их параметров, а также высокой вычислительной трудоемкостью решения оптимизационных задач со стохастическим критерием.
Задача параметрического синтеза технических устройств и систем [1, 2] состоит в выборе номинальных значений параметров элементов исследуемого устройства хн = (х1н, ..., х ), обеспечивающих максимум вероятности выполнения условий работоспособности в течение заданного времени:
хн = а^тахР{Дхн, ?) е Dx, Уt е [0, Т]}, (1)
где Х(хн, 1) — случайный процесс изменения параметров; Dx — область работоспособности; Т — заданное время эксплуатации устройства.
Исходными данными в рассматриваемой задаче являются условия работоспособности, задаваемые обычно в виде допусков на выходные параметры:
Уп
у(0 < У,
(2)
1 Работа выполнена при поддержке РФФИ (проект № 0508-013998) и ДВО РАН (Программа № 15 ОЭММПУ РАН).
и математическая модель системы у(х(1)), определяющая зависимость выходных параметров
У(1) = {У/ (0}™= 1 от параметров элементов х(1) = = {х{( 1)}”= 1 в виде непрерывных зависимостей
У/0 = у/х(0), j = 1, т .
В большинстве случаев такая зависимость задается не в явной, а в алгоритмической форме, в частности, через численные решения систем уравнений (дифференциальных или алгебраических), описывающих процесс функционирования исследуемой системы.
Область пространства параметров элементов
х е Я ”, во всех точках которой выполняются условия (2), называется областью работоспособности Dx.
Л = {х е Я” : ут1п < у(х) < Утах}.
Условия работоспособности формируют в пространстве выходных параметров область допустимых значений Лу:
ЛУ = {У е Ят : Утш т У т Утах},
представляющую собой т-мерный ортогональный гиперпараллелепипед.
Решение задачи (1) базируется на анализе взаимосвязей выходных характеристик системы и параметров составляющих ее элементов (внутренних
параметров) с учетом технологического разброса, температурного и временного дрейфов этих параметров. Определенной характеристикой возможности системы выполнять заданные функции в условиях параметрических возмущений служит область работоспособности Dx, построенная в координатах параметров элементов системы. Совокупность значений внутренних параметров может быть представлена изображающей точкой в и-мер-ном пространстве этих параметров. Для обеспечения работоспособности системы эта точка должна находиться внутри области Лх. При этом расстояние от изображающей точки до границ Лх можно рассматривать как некоторый запас работоспособности системы.
Построение тем или иным методом области работоспособности Лх является, по сути, первым этапом решения общей задачи параметрического синтеза системы. Следующие этапы заключаются в назначении допусков и выборе номинальных значений параметров элементов проектируемой системы. Выполнение этих этапов связано с расчетами некоторых, как правило, статистических показателей, например, таких, как вероятность безотказной работы Р(х, 1) на интервале эксплуатации Т и т. п. Часто необходимая априорная информация о вероятностных свойствах отклонений параметров от расчетных значений отсутствует или является недостаточно полной. В таких случаях вместо статистических показателей могут быть использованы некоторые детерминированные критерии типа «запасов» (работоспособности, надежности и др.) [3, 4]. Расчет таких показателей не требует знания полных вероятностных характеристик случайных величин, фигурирующих в математической модели системы. В сравнении со статистическими они имеют и более ясную физическую интерпретацию, например, измеряемое определенным образом расстояние от номинальной точки до границы области Лх. Вместе с тем, расчет «запаса» работоспособности по каждому из внутренних параметров затрудняет оценку влияния размеров этих «запасов» на выполнение того или иного условия работоспособности системы в целом. Во многих случаях оказывается целесообразным рассматривать показатели типа «запасов» относительно не внутренних, а выходных параметров системы, ограничения на которые и составляют заданные условия работоспособности системы.
1. КРИТЕРИЙ ЗАПАСА РАБОТОСПОСОБНОСТИ
Как уже отмечалось, запас работоспособности можно рассматривать на уровне внутренних параметров (параметров элементов) или выходных параметров системы.
Запас работоспособности первого типа — на уровне внутренних параметров — позволяет оценить степень удаленности вектора внутренних параметров от границ области работоспособности, а, следовательно, пределы возможных вариаций параметров элементов, при которых не нарушаются условия работоспособности. Задача оптимального параметрического синтеза в этом случае сводится к нахождению таких точек внутри области работоспособности Лх (выбора таких номиналов параметров), которые находятся на максимальном в смысле выбранного критерия расстоянии от ее границ.
Область работоспособности Лх обычно неизвестна, поэтому выполнение условий работоспособности при выбранных внутренних параметрах проверяется в результате вычисления соответствующих выходных параметров и сравнения их с требованиями технического задания (областью допустимых значений выходных параметров Лу).
Можно говорить о запасе работоспособности второго типа, представляющем собой меру удаленности вектора выходных параметров у = (уДх), У2(х), ..., Ут(х)) от заданных требованиями технического задания границ области Лу.
Поскольку задача параметрического синтеза состоит в выборе номинальных значений внутренних параметров, будем называть выбор значений параметров по критерию запаса работоспособности первого типа прямой задачей, а выбор по критерию запаса работоспособности второго типа — обратной.
Любая комбинация внутренних параметров х = (х1, х2, ..., X”) однозначно определяет некоторую совокупность выходных величин у = (У1(х), У2(х), ..., Ут(х)) и, таким образом, некоторую точку у е Ят в т-мерном пространстве выходных параметров. При этом обратное отображение не всегда однозначно: одному и тому же набору значений выходных параметров могут соответствовать несколько различных векторов внутренних параметров.
Будем говорить, что совокупность внутренних параметров представляет допустимое решение (удовлетворяющее условиям работоспособности), если соответствующий им вектор выходных параметров лежит в т-мерном полиэдре Лу, задаваемом выходными ограничениями.
Чтобы иллюстрировать сказанное выше, рассмотрим случай с двумя внутренними параметрами (и = 2), значения которых обозначим х1 и х2, и двумя выходными параметрами у1 и у2. Область, образуемая ограничениями на выходные парамет-
ры (область работоспособности Ву) представляет собой в данном случае прямоугольник в плоскости выходных параметров (рис. 1).
Прямоугольник отображается на плоскость внутренних параметров в виде фигуры, выделенной жирными линиями (рис. 2). На рис. 2 также показано возможное решение задачи выбора номиналов х1н и х2н.
Рассмотрим более подробно задачу выбора оптимальных значений номиналов параметров по критерию запаса работоспособности на уровне выходных параметров.
Случай одной выходной и п внутренних переменных. Пусть условия работоспособности исследуемой системы заданы интервалом допустимых значений выходного параметра у(х): утіп т у(х) т утах.
При отсутствии какой-либо информации о закономерностях отклонений выходного параметра от номинального значения, оптимальным будет такое значение у, при котором обеспечивается максимальный запас работоспособности. Очевидно, что в рассматриваемом случае это
О _ у тах + у тіп Ун 2 ;
(3)
при котором запас работоспособности равен половине интервала допустимых значений.
Для нахождения номинальных значений внутренних параметров необходимо решить следующую оптимизационную задачу:
или
Хн = ащтіп| Ун - у(х)|
х н = ащтіп( ун - у(х))ь,
X
(4)
(5)
где Ь — положительное целое четное число.
Можно показать, что номинальное значение параметра, соответствующее середине поля допуска (3), будет оптимальным и в случае любого симметричного закона распределения вероятностей рассматриваемого параметра (например, нормального или равномерного). Оно гарантирует не только максимальный запас работоспособности, но и
-------------------------------------------------- максимальную вероятность нахождения параметра
Рис. 1. Область работоспособности Бу в пространстве выходных у(х) в пределах поля допуска [ут|п, УтахМ
Случай т выходных и п внутренних параметров. Если условия работоспособности заданы в виде системы независимых допусков (пределов допустимых вариаций), связанных с каждым из выходных параметров
параметров
у/тіп
т у(х)
т у/max,
j _ 1, т ,
то выбор вектора номиналов, максимизирующих запас работоспособности по каждому из выходных параметров, также не вызывает затруднений. Очевидно, что он будет определяться как центр пространства (гиперпараллелепипеда) допусков
о _ у і тах ун]
+ у і
2
і тіп
j _ 1, т .♦
Однако решение задачи нахождения опти-
0
мального вектора внутренних параметров хн =
Рис. 2. Область работоспособности В в пространстве
хп н) не всегда удается свести к реше-
внутренних параметров
нию т одномерных задач.
Необходимо ввести меру близости двух векторов у(х) и ун, которая будет служить в качестве критерия оптимальности (целевой функции). Во многих случаях следует учитывать различную степень важности отдельных выходных параметров.
Так, например, для схемы усилителя низкой частоты используют обычно следующие показатели (выходные характеристики): у1 — усиление по напряжению, у2 — амплитуда выходного импульса по низкой частоте, у3 — напряжение неискаженного сигнала на нагрузке, важность которых, вообще говоря, различна.
С учетом сказанного, в качестве целевой функции можно предложить следующий обобщенный показатель:
Y(x) =
j = 1
2(У,Сx) - (У/min + У/max)/2)
/max/
У,
/max y/min
(б)
где Уj — значение ]-й выходной характеристики, а. — положительный весовой коэффициент, Ь/. — положительное целое четное число. Значение выражения в квадратных скобках заключено в отрезке [—1, 1] для любой выходной характеристики у,, удовлетворяющей требованиям технического задания.
Значение целевой функции У(х) является мерой близости решения центру области возможных решений, т. е. подпространству пространства выходных параметров, в котором значения всех выходных переменных равномерно отстоят от своих верхних и нижних пределов.
Решение задачи оптимального выбора номинальных значений внутренних параметров по критерию запаса работоспособности сводится в данном случае к нахождению вектора внутренних параметров, доставляющих минимум обобщенному показателю (6).
Показатель (6) обладает свойством принимать малые значения в тех случаях, когда все выходные ограничения выполнены, и резко возрастать, если хотя бы одно ограничение нарушено. Использование такого показателя во многом напоминает метод построения поверхности отклика или метод штрафных функций, известный в теории оптимизации.
Решение задачи выбора номинальных значений внутренних параметров проектируемой системы, при которых обеспечивается максимальный запас работоспособности, может быть осуществлено методами поисковой оптимизации. При этом решение обратной задачи не требует построения области работоспособности в пространстве внутренних параметров Dx. Однако для достаточно сложных многопараметрических систем процесс решения
оптимизационной задачи (6) оказывается крайне трудоемким. Одним из путей снижения временных затрат на решение рассматриваемых оптимизационных задач может быть распараллеливание вычислительного процесса.
2. ПАРАЛЛЕЛЬНЫЙ АЛГОРИТМ ДИСКРЕТНОЙ ПОИСКОВОЙ ОПТИМИЗАЦИИ
При программной реализации параллельных алгоритмов представляется целесообразным использование возможностей как современных многопроцессорных вычислительных систем, так и распределенных многомашинных комплексов, связанных локальной сетью. Реализация подобных алгоритмов на многопроцессорных машинах, работающих под управлением операционных систем, поддерживающих многопоточность, не вызывает принципиальных затруднений. В то же время, в распределенных гетерогенных средах, необходимо самостоятельно реализовать механизмы загрузки данных, синхронизации процессов и балансировки вычислительной нагрузки между компонентами комплекса.
Главным критерием качества распараллеливания вычислений служит сокращение общего времени решения задачи. Возможности для распараллеливания вычислений ограничиваются не только числом имеющихся процессоров, но и особенностями вычислительного алгоритма, который может оказаться принципиально последовательным. К задаче параметрической оптимизации по критериям вида (6) применима модель передачи сообщений и, так называемая, master-slave-парадигма параллельного программирования. Один из процессоров назначается главным (master), он осуществляет динамическую балансировку загрузки, рассылает задания остальным подчиненным процессорам (slave) и формирует результаты. Распараллеливание здесь базируется на декомпозиции последовательного алгоритма вычислений, в качестве единицы параллелизма выступает задача однократного расчета модели системы (моделирования системы).
Универсальным и достаточно эффективным путем решения задачи (6) может стать применение параллельного аналога метода сканирования (слепого поиска).
Сущность метода заключается в том, что вся допустимая область пространства параметров разбивается на элементарные ячейки, в каждой из которых по определенному алгоритму выбирается точка: в центре ячейки, на ребрах или вершинах. В каждой ячейке осуществляется последовательный просмотр значений целевой функции и на-
m
хождение среди них экстремального значения. Точность метода, естественно, определяется тем, насколько плотно располагаются выбранные точки в области поиска.
Основное достоинство метода сканирования состоит в том, что при достаточно густом расположении точек всегда гарантируется отыскание глобального экстремума. Однако для этого требуется значительный объем вычислений, снизить который можно путем распараллеливания алгоритма. Наиболее простой алгоритм поиска экстремума методом сканирования (поиска на сетке переменных) заключается в том, что по каждой независимой переменной задаются приращения в соответствующем порядке, обеспечивающем заполнение всей исследуемой области равномерной и достаточно густой сеткой.
Поскольку номинальные значения параметров схемных элементов должны принадлежать ряду стандартных значений, регламентированных техническими условиями или ГОСТами, иногда предпочтительнее искать оптимальный вектор номиналов параметров на дискретном множестве номиналов Бн, соответствующем стандартным значениям и ограниченном областью допустимых значений Вх. В случае, когда эта область неизвестна, для каждого из параметров х. необходимо задать пределы их возможных изменений, например, построив описанный гиперпараллелепипед [5], или на основе известных допусков на каждый из внутренних параметров.
Пусть известны множества номиналов для каждого из п выбираемых параметров схемных элементов исследуемой системы:
12 x1 < x1
r1
< ... < x1
г 1 2 r 1 ,
nom1 = {x1, x1, ..., x1 },
1 2 rn 1 2 rn
nomn = {x« . xn , •••> xn h xn < xn < ... < xn .
Решение данной задачи предлагается провести как двухэтапную параллельную процедуру.
На первом этапе предлагается ограничить пространство поиска брусом допусков, проведя аппроксимацию области работоспособности описанным параллелепипедом. Используя изложенный в статье [5] параллельный алгоритм построения описанного бруса, построим брус В0 с Bd
B0 = {x е Rn| a0 m xt m b0, i = lTn},
0 ,0
где a t = min x.; b . = max x..
x e Dx x e Dx
Важно отметить, что при построении описанного параллелепипеда нет необходимости нахождения области работоспособности в пространстве
Рис. 3. Схема параллельвдго алгоритма дискретвдй oптимизации на мнoжестве нoминалoв парамет|юв
внутренних параметров Бх, что существенно уменьшает трудоемкость предлагаемого алгоритма [5].
На втором этапе сформируем дискретное множество номиналов внутренних параметров
тчвн ( вн 0 ^ 1 ^ г 0 . Ї Ї
ин _ {хн : а і т х, < Ь , , хі є пот,, і = 1, п}.
Ввн тлвн
каждой точке хн дискретного множества Вн необходимо найти значение целевой функции.
Искомый оптимальный вектор номиналов хнпт находим, решая задачу:
arg min Y( xнн).
(7)
Задача (9) эквивалентна задаче нахождения оптимального решения в смысле (4) или (5), а в многомерном случае доставляет минимум обобщенному показателю (6) на дискретном множестве номиналов параметров.
В простейшем случае нахождение решения задачи (7) сводится к полному перебору элементов
тлБН
множества Вн , для каждого из которых осуществляется расчет функции У( х Нн). Учитывая цикличность процедуры вычисления целевой функции, можно применить параллелизм по данным. Схема параллельного алгоритма дискретной оптимизации по критерию запаса работоспособности приведена на (рис. 3).
Множество DНН разбивается на непересекаю-
к
тлвн I I г тлвн .
щиеся подмножества Бн = . 11 {Бн.}, при этом каждому у-му процессору назначается свое под-
БИ БИ
х„ є D„
И
И
G8
CONTROL SCIENCES № О • 2007
множество DНН исходных данных. Таким образом, каждый у-й процессор рассчитывает целевую функцию для всех элементов множества DНН и находит оптимальный вектор номиналов параметров для своей подобласти. Результаты передаются главному процессору, который выбирает оптимальный вектор номиналов по всей области DНН. Такое разбиение всего множества поиска на непе-ресекающиеся подмножества составляет суть блока диспетчеризации параллельного распределенного процесса.
Для симметричного вычислительного кластера, состоящего из к равных по мощности вычислительных узлов, общее число точек разбивается на равные части для каждого из подчиненных процессов. В случае несимметричного кластера необходима предварительная процедура оценки трудоемкости типовой процедуры метода оптимизации, в качестве которой выступают однократное моделирование работы системы, проверка условий работоспособности и вычисление значений критерия запаса работоспособности. Вычислительная нагрузка делится между компонентами комплекса пропорционально их производительности [6].
По окончании работы программы диспетчеризации вычислительного процесса каждому вычислительному компоненту комплекса рассылаются
ГчВн
границы его подмножества Dнj• исходных данных. По окончании счета главный процессор получает результаты от подчиненных и формирует окончательные результаты дискретной оптимизации на
ГчВн
всем множестве Dн .
Данный параллельный подход обладает достаточно высокими характеристиками производительности, так как обмен между процессорами сводится к минимуму — назначению заданий и заключительной передаче результатов. Выполнение испытаний на процессорах не синхронизируется.
При одинаковых мощностях подобластей DНН и равных временных затратах на вычисление целевой функции на симметричных кластерах ускорение параллельного алгоритма сканирования практически достигает линейного. Наличие внутренних циклов порождает высокую масштабируемость алгоритма, главное условие эффективной реализации которого состоит в пропорциональной загрузке всех участвующих в вычислениях процессоров.
ЗАКЛЮЧЕНИЕ
Основные трудности, возникающие при решении задачи оптимального параметрического синтеза технических устройств и систем с учетом отклонений их параметров от расчетных значений, связаны с дефицитом информации о закономерностях этих отклонений и ее высокой вычислительной трудоемкостью.
Применение критерия запаса работоспособности позволяет получить решение задачи параметрического синтеза в случае отсутствия информации о стохастических закономерностях вариаций параметров, а применение технологии распределенных (параллельных) вычислений предложить эффективные вычислительные алгоритмы. Предложенный алгоритм, ориентированный на реализацию в распределенной вычислительной среде, позволяет выбрать номинальные значения параметров, обеспечивающих максимальный запас работоспособности исследуемой системы. Важный момент алгоритма состоит в том, что процессоры и станции локальной сети могут быть как однотипными, так и отличающимися друг от друга по своим вычислительным характеристикам. Такая схема повышает надежность вычислений, устраняется возможность блокировки системы.
ЛИТЕРАТУРА
1. Анализ и оптимальный синтез на ЭВМ систем управления / Под ред. А.А. Воронова и И.А. Орурка. — М.: Наука, 1984.
2. Абрамов О.В. Параметрический синтез стохастических систем с учетом требований надежности. — М.: Наука, 1992.
3. Антушев Г.С. Методы параметрического синтеза сложных технических систем. — М.: Наука, 1989.
4. Норенков И.П., Мулярчик СГ, Иванов С.Р. Экстремальные задачи при схемотехническом проектировании в электронике. — Минск: Изд-во БГУ им. В.И. Ленина, 1976.
5. Абрамов О.В., Катуева Я.В. Технология параллельных вычислений в задачах анализа и оптимизации // Проблемы управления. — 2003. — № 4. — С. 11—15.
6. Катуева Я.В. Балансировка загрузки несимметричного вычислительного комплекса при решении задачи статистического оценивания // Информатика и системы управления. — 2006. — № 2(12). — С. 88—93.
S (4232) 31-02-02, e-mail: abramov@iacp.dvo.ru,
gloria@iacp.dvo.ru, nazardim@iacp.dvo.ru
Статья представлена к публикации членом редколлегии
Б.Г. Воликом. □
ПРОБЛЕМЫ УПРАВЛЕНИЯ № О • 2007
G9