Приближенная оптимизация управления в параллельных вычислениях Текст научной статьи по специальности «Математика»

УДК 517.977:519.68 ББК 22.193:32.973.26-018.2

© В.И. Гурман, Е.А. Трушкова, А. О. Блинов

Институт программных систем им. А.К. Айламазяна РАН, Переславль-Залесский

E-mail: [email protected])

Приближенная оптимизация управления

1

в параллельных вычислениях

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

Ключевые слова: оптимальное управление, приближенные методы, улучшение управления, приближенно-оптимальный синтез, параллельные вычисления

© V.I. Gurman, Е.А. Trushkova, A.O. Blinov

A.K. Aylamazyan Institute of Program Systems, RAS, Pereslavl’-Zalesskiy E-mail: [email protected]...

Approximate optimization of control in parallel computations

The article is devoted to new approximate methods of optimal control on the base of sufficient optimality conditions and extension principle and their realization in parallel computations. The corresponding set of algorithms and software for the problems of improvement and approximate optimal synthesis of control are briefly described. The examples are given to demonstrate the efficiency of the methods proposed in parallel computations.

Key words: optimal control, approximate methods, approximate optimal synthesis, parallel computations.

Введение

Приближенные и вычислительные методы - обширная и ставшая самостоятельной область исследований и разработок в теории оптимального управления, нацеленных на эффективное решение практических задач. Исторически развитие этих методов началось одновременно с созданием современной теории оптимального управления. В числе основоположников отметим Д. Е. Охоцимского, В. А. Егорова, Т. М. Энеева, Л.И. Шатровского, Дж. Келли.

Одно из направлений в этой области базируется на достаточных условиях оптимальности Кротова и принципе расширения. Его плодотворность подтверждается многолетними разработками [1] - [4], отличающимися значительным многообразием подходов и результатов. В основном это различные итерационные процедуры улучшения управляемых процессов, как и в других школах. Спецификой является возможность априорного оценивания получаемых приближенных решений и использование характерного свойства вырожденно-сти прикладных задач и соответствующих специальных методов для поиска начальных приближений, что, как известно, является критическим моментом при использовании итерационных улучшающих алгоритмов. В последнее время интенсивно развиваются новые методы, ориентированные на параллельные вычисления на суперкомпьютерах кластерного

1 Работа выполнена при частичной финансовой поддержке Российского фонда фундаментальных исследований (проекты № 08-01-00274 и № 09-01-00170).

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

Цель данной статьи - представить собственный опыт работы в этом направлении, связанный с созданием, развитием и приложениями программно-алгоритмического комплекса улучшения и приближенно оптимального синтеза управления (ПК КСОМ) [5].

1. Комплекс алгоритмов улучшения и оптимизации управления

Алгоритмы ПК КСОМ формализованы исходя из задач управления непрерывной и дискретной системой в стандартной форме:

х = /(г,х,и), ге Т = [г,,гР], (1)

х(г +1) = / (г, х(г), и (г)), г е Т = г,, г! +1,..., гР, (2)

х е Х(г) с Rп, и е и(г, х) с Rр, х(гР) е Г, I = Р(х(гР)),

(г,, х(г,),гР) заданы, Р : Яп ^ Я . Реализуемые здесь методы построены или строятся на

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

В соответствии с теорией Кротова при помощи произвольной функции (р(г, х) строятся следующие конструкции:

Я(г, х, и) = (р(г +1, / (г, х, и)) -^(г, х), G( х) = Р (х) + (р(гР, х) - (р(г,, х(г,)).

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

Комплекс алгоритмов состоит из трех крупных блоков (АППРОКСИМАЦИЯ, СИНТЕЗ, УЛУЧШЕНИЕ) и предназначен для решения трех типов задач: (I) глобальный приближенно-оптимальный синтез управления; (II) улучшение и локально-оптимальный синтез управления и (III) управление в реальном времени. Задачи типа (I) в наибольшей степени отвечают указанной выше проблеме синтеза, однако требуют наибольших вычислительных ресурсов. Задачи типа (II) возникают в тех случаях, когда вычислительные трудности глобального синтеза оказываются непреодолимыми и заставляют решать проблему в два этапа. Вначале ищется оптимальная программа управления итерационными методами последовательного улучшения, и соответствующая траектория. Затем в некоторой достаточно малой окрестности последней строится приближенно-оптимальный синтез с требуемой точностью при предположении, что действующие возмущения не приводят к выходу результирующих траекторий за пределы этой окрестности. Локализация позволяет существенно уменьшить требования к вычислительным ресурсам за счет сужения расчетной области и возможного упрощения модели объекта с целью использования аналитических методов. Задачи типа (III) связаны с практической реализацией найденных законов управления, вычислительных экспериментов с целью адаптации их к реальным условиям функционирования управляемых систем. Здесь большое значение имеет принятая концепция управления как дискретного во времени. Она позволяет принципиально разделить этапы наблюдения-анализа текущей ситуации и принятия адекватных управленческих решений. Такое разделение делает полученные результаты не зависящими от принятого метода решения задачи наблюдения, следовательно, существенно расширяет круг возможных практических приложений в сравнении с подходами, где оба этапа не разделены.

Блок АППРОКСИМАЦИЯ предназначен для аппроксимации моделей объектов управления и функций Кротова при решении задач улучшения и синтеза. По существу он универсален и может использоваться для любых многомерных объектов. Аппроксимация про-

изводится в среднем по методу наименьших квадратов (МНК). Этот хорошо известный метод в данном случае дает возможность «развязать» аппроксимирующую конструкцию и сетку значений аппроксимируемой функции, которые при многомерной интерполяции жестко связаны. В частности, он позволяет генерировать регулярные («прямоугольные») сетки с постоянным или переменным шагом по каждой переменной, что очень важно для распараллеливания вычислений. Некоторые детали реализуемых в нем алгоритмов приводятся в нижеследующих описаниях основных блоков, где они интенсивно применяются.

2. Описание блока СИНТЕЗ

В этом блоке сосредоточены алгоритмы, реализующие метод приближенного синтеза оптимального управления для дискретных моделей управляемых систем на основе среднеквадратической аппроксимации (в частности, на интерполяции) соотношений динамического программирования при полиномиальном представлении функции Кротова-Беллмана [6]. При этом рассматривается постановка задачи «со свободным правым концом»: X(г) = Г = Ип. Фазовые ограничения, в том числе на правом конце, могут быть учтены с помощью известного метода штрафов.

Функция Кротова р выбирается в общем случае так, чтобы в окрестности некоторой опорной траектории приближенно выполнялись соотношения типа Беллмана:

Р(г, х) = р(г, х) - Н (г, х,р(г +1, х)) = 0, G( х) = р(гР, х) + (1 - а)Р (х) = 0,

Н (г, х,р(г +1, х) )= тах \ ррг +1, / (г, х, и)) - — а(и - и (г))2 |.

м<=и (г ,х)^ 2 )

Опорной может быть приближенно оптимальная, либо улучшаемая траектория, если синтезируется улучшающее управление, и (г) - улучшаемое управление, 0 < а < 1 - регулятор метода улучшения. При отсутствии опорной траектории принимается заведомо а = 0.

Аппроксимация функции Кротова производится некоторым многомерным полиномом

а

р(г, х) = !/т,(х), (3)

г=1

где g (х) = (^ 1 (х),..., gа (х))т - набор заданных базисных функций, а

/(г) = (/ (г),..., /а (г ))т - набор коэффициентов, подлежащих определению из условий

тп I (/ (г)g(хр (г)) - Н(г, хр (г), / (г + %(хр (г))))2, (4)

т1п I (/т (гР)g(хр (гР)) + (1 - а)Р(хр (гР )))2, (5)

/(гР)

в - номер узла сетки. Условия (4), (5) сводятся к решению системы линейных алгебраических уравнений А(г )[/а (г)] = В(г) относительно коэффициентов /(г) аппроксимирующего полинома.

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

\ т2 \ тп

рРг,х) = £ (х‘)к I (х2 У2 - I /и,.п (х" У

Л=0

к =0

" ) к

V к" =0 ))

допускающая интерполяцию на регулярной сетке с числом узлов т1 по переменной х{. Второй тип - регулярный полином порядка п

1

р(г, х) = у(г) + /((г) х + —хта(г) х + ., где / (г) - п -мерный вектор, о(г) - матрица п х п.

В вычислительном плане наиболее удобно строить постоянную во времени прямоугольную сетку Xp(t) = const, покрывающую область, где строится синтез. С этой целью (2) преобразуется к новым переменным - отклонениям от опорных значений.

Близость полученного приближенного синтеза оптимального управления ~(t, x(t)) к строгому оптимуму можно определить с помощью следующей верхней оценки Кротова:

A(t) = A(t +1) + max P(t, x) - min P(t, x)

x x

A(tF) = max G(x) - min G(x),

x x

где x пробегает заданную область. Найденное управление тем ближе к оптимальному, чем меньше эта оценка. Возможность вычисления оценки - это важное преимущество перед «чистым» методом Беллмана. Она позволяет организовать регулярную процедуру уточнения приближённого решения за счет увеличения числа узлов интерполяции и их расположения в фазовом пространстве, а также дает критерий ее остановки.

Схема соответствующего алгоритма представлена в общем виде на рис. 1 и более подробно на рис. 2, 3.

Важную роль здесь играют алгоритмы задания эффективных начальных приближений по методу восстановления функции цены (МФЦ) [6] в случае магистральных решений, характерных для многих прикладных задач оптимального управления. В этом случае аппроксимирующая конструкция для функции Кротова модифицируется следующим образом: 1=1 + р, где р - приближение по МФЦ, (р1 - полином (3).

Рис. 1

Рис. 2

Щ t,x ) = 'Гчг,(1 )gi( х >

1=1

Рис. 3 22

Достаточно легко реализовать параллельность описанного алгоритма при переборе значений в. Область формирования сетки разбивается на части, в каждой из которых строится часть матрицы А(г) и часть столбца свободных членов В(г) для системы уравнений с помощью исходных значений в узлах текущей части. Общая матрица получается в этом случае, как сумма всех построенных «частичных» матриц (это же справедливо и для свободных членов). На рис. 2, 3 этот вид потенциального параллелизма обозначен цифрой 1. Решение системы линейных алгебраических уравнений (обозначено цифрой 2 на рис. 2, 3) и поиск максимума функций многих переменных (обозначены цифрой 3 на рис. 2, 3) можно отнести к другим видам потенциального параллелизма. Как показывает практика, программная реализация этих видов параллелизма эффективна лишь при больших размерностях п , р исходной динамической системы (2).

3. Описание блока УЛУЧШЕНИЕ

В этом блоке сосредоточены алгоритмы, реализующие методы первого и второго порядка улучшения управления для дискретных моделей управляемых систем. Один из них основан на полиномиальном задании функции Кротова (3) и среднеквадратической аппроксимации соотношений Кротова-Беллмана (4), (5) в окрестности улучшаемой траектории и по существу описан в предыдущем разделе. Получаемый синтез при достаточно большом значении а обеспечивает улучшение.

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

Имеется элемент т = (т, л1 (г), и1 (г ))е О , требуется найти элемент ш11 = (т, л11 (г), и" (г))е О , такой, что I(ш11) < I(т1). Общие конструкции метода улучшения управления приведены в [4], где ш11 ищется как аппроксимация решения задачи: у(г +1) = в (г, у(г), v(г)), г е Т = {г,, г, +1,..., гР },

1 Р

у0(г +1) = у0(г) + -^v2(г), у(г,) = 0, у0(г,) = 0,

2 г=1

ра (у°(г^), У(гF))= ау°(гР) + (1 - а)Р(у(гР) + л1 (г)) ^ шш, где а - некоторое действительное число из полуинтервала (0,1] (регулятор метода), В (г, у(г), v(г)) = / (г, у(г) + л1 (г), v(г) + и1 (г))-/ (г, л1 (г), и1 (г)), л = у + л1, и = V + и1.

1

Функция Кротова ищется в виде (р(г, у0, у) =у(г) + /0(г) у0 + /т (г) у + — у1 о (г) у. Значения у(гР), /0 (гР), / (гР), о(гР) находятся из условий

у(гР) = -Ра(0), /0(гр) = -а, /(гр) = --^(Ра(01 уАл,п)-Ра(0)), I = 1П,

о« (гР) = -^ (Ра (012уАх,/;) - Ра (01 уАх,,;) + Ра (0)), I = Щ,

о, (гР ) = -(Ра (0 1 ^/ . ) - Ра (01 ) - Ра (01 ^/ . ) + Ра (0)) I, . = Щ,

где у - некоторое действительное число (регулятор метода), Ал;, Аи. - масштабы переменных л, и.. Здесь через Ра (01 уАл{/;) обозначено значение функции Ра в точке

(у “, Уl,., у^ у, УI■+l,., уп) = (0,0,■■■,0,1Аxi ,0,.,0).

Вводится в рассмотрение функция

На (г, у, V) = у(г +1) - а ][^2 + /т (г +1) в (г, у, V) +1 в (г, у, v)т о(г +1) в (г, у, V).

2 г=1 2

Эта функция разлагается в ряд в окрестности нуля:

1

На (г, у, V) = На (г ,0,0) + ут Ну (г ,0,0) + - ут Нщу (г,0,0) у +

Из условия

получается

+ ^ ^Нт (г ,0,0) + ^Нау, (г,0,0) у + 2 ^Наш (г ,0,0» | [хтНт (г,0,0) + ^На^ (г,0,0)у + 2^Нат (г,0,0)^ = 0, г е т \ {гР },

1

V" (г, у) = аге шах^ V Нт (г ,0,0) + V Нау, (г,0,0) у + 2 vTHavv (г ,0,0» J = = -Н-^ (г,0,0)Н (г ,0,0) + Н^ (г ,0,0) у ) у(г) = На (г,0,0), /°(г) = -а, /(г) = Ну (г,0,0), о(г) = Нщу (г,0,0), г е т \ {гР }.

Последние равенства, после замены производных их разностными аналогами, принимают вид:

N. (г) =

Vй(г,у) = -^(г)(Ма(г) + РТа(г)у), ге т\{гР},

^ Н

(ЗАщ)

На (0 I 28Аи./ п+.) - На (0 I 5Аи./ п+.) + На (0)) . = 1, р,

1 / п+1

^(г) = —^--(На(01 ЗАщ. . ,ЗАщ )-На(01 ЗАщ. .)

к1() З АщАи.1 а( к/п+к, "п+.) а( к/п+к)

к ]

На (0 I ЗАщ / п+. ) + На (0)) к, . = 1Р,

М' (г) = ЗАТ (На (01 ЗАи./п+1) - На (0)), . = 1,Р,

Р. (г) =---------------------1-(На (01 ?Ал; ,ЗАи. ) - На (01 ]Ал. )- На (0 I ЗАи; ) + На

уЗАл Аи 1 / 1 / п+. 1 / ' 1 / ”+ ' а

1 / п+1 '

г = 1, п, ] = 1, р,

у(г) = На (0), /0(г) = -а, / (г) = -2_ (На (01 уАл,,.) - На (0)), г = 1, п, уАл, о* (г) = 7^ (На (012уАлг/) - На (01 уАлг,) + На (0)), г = й,

(уАл г)

О;(г) =-

На (01 уАл1 /г ,уАл, /;) - На (01 уАлг/,.) - На (01 уАл,/ 1) + Н

(0)),

у1 Ал,Ас/ а'‘ ■ 1 /^

г, 1 = 1,п, ге Т \{гР},

где З - некоторое действительное число (регулятор метода). Здесь через На(01 уАл{п,8Аии ) обозначено значение функции На в точке

(г, Уl,., у^ у, , УI■+l,., уп, ^•■■, vl-l, V, vi+l,■■■, VР) = ^Д.-.Д^Ц ,0,...,0,0,...,0,ЗАи1Д...Д. Полагается

и" (г, л) = Vя (г, у + л1 (г)) + и1 (г), г е Т \ {гР }.

С помощью уравнений системы (2) находится соответствующая пара (и11 (г), л11 (г)). Полученный таким образом элемент ш11 = (т,л11 (г),и11 (г))е О и будет улучшением.

Замечание 1. В случае алгоритма первого порядка функция Кротова ищется в виде

р(г, у0, у) = у(г) + /0(г) у0 + /Т (г) у.

Следовательно, для этого случая все соотношения вышеописанного алгоритма второго порядка остаются в силе, если только положить о(г) = 0, г е т .

Схема описанного алгоритма представлена на рис. 4.

Вход

У1 . Ґ- . -т-у\+У1+р

н: р:іІ,ір:хІ є Д":/ : Я

Я”: Л: Л" -» Я:

*£«,>*£«>1 = =1 • • (0>* = ^ “1

Л. =0.1;/ :=0.5;г := 1.5; А _ := 0.5; 6 :=0;<5 :=0.5; а У тлі * / шах / шш шах ’ /і# :=0.5;Дг, := д£„ -дс^.і = 1„..,и;Дн; := = 1 />;е:=0.001

Алгоритм улучшения II порядка

Генерация сетки

/?:=1

Расчет конструкций метода

х?0),мЛ(і),ра(і),

и!(і)\?:=иП(і)\р

х!(і)\,:=хп(і)\!

©

Р“ := сїї^ шах(р(.х’(1Р ))- р{х1!( ) )р ))

иа(і):=ип(і)\

Выход

Рис. 4

Из блок-схемы алгоритма (рис. 4) видно, что расчет конструкций метода и прямой счет при нахождении управления в форме синтеза могут выполняться независимо для каждого набора параметров алгоритма (для каждого значения в). На рис. 4 этот вид потенциального параллелизма обозначен цифрой 1. При этом число итераций для различных наборов параметров заранее неизвестно и в общем случае различно, что затрудняет процесс равномерной загрузки вычислительных узлов при программной реализации данного вида параллелизма. К другим видам потенциального параллелизма можно отнести операции матричного исчисления (обозначены цифрой 2 на рис. 4) и поиск максимума функций многих переменных (обозначены цифрой 3 на рис. 4). Как показывает практика, программная реализация этих видов параллелизма эффективна лишь при больших размерностях п , р исходной динамической системы (2).

4. Параллельная реализация и вычислительные эксперименты

В целом анализ показывает, что данный алгоритмический комплекс обладает большим потенциальным параллелизмом. С учетом того, что число итераций в алгоритмах улучшения заранее неизвестно, динамическое распараллеливание предпочтительнее любых других типов. Эту концепцию автоматического динамического распараллеливания программ реализует Т-система — оригинальная российская разработка [7]. Она предоставляет язык программирования Т++ (очень простой параллельный диалект С++), который предназначен для эффективной реализации динамического распараллеливания.

Вышеописанные алгоритмы были реализованы в виде исполняемых модулей ПК КСОМ с учетом возможностей распараллеливания. Была проведена большая серия вычислительных экспериментов, с использованием методических примеров и прикладных задач, связанных с исследованием маневров безопасной нештатной посадки вертолета с определением границы безопасной зоны [8] и стратегий устойчивого развития регионов на их социо-эколого-экономических моделях.

Задача о посадке вертолета показательна тем, что потребовала использования различных аппроксимаций имитационной модели движения вертолета, не имеющей полного аналитического описания: силы и моменты, действующие на вертолет, были заданы лишь в форме компьютерных программ. Вначале была получена линейная аппроксимация модели объекта, на которой был проведен качественный анализ, в результате него сформировалось грубое начальное приближение глобально-оптимального управления. Далее это управление улучшалось на более сложной, нелинейной аппроксимации модели объекта с использованием рассмотренного выше алгоритма улучшения. Результаты для наиболее жесткого из рассмотренных сценариев нештатной ситуации уже после трех итераций позволили сделать вывод о повышении границы опасной зоны на 15% за счет оптимизации соответствующего посадочного маневра. Эффективность параллельной реализации программ при решении этой задачи демонстрируется в таблице 1.

Таблица 1

Программа аппроксимации модели динамической системы с помощью МНК

Число процессоров (ядер) 1 3 5 7 9 11

Время работы программы, с 3338.35 1237.14 729.92 632.00 631.21 588.02

Ускорение 1 2.70 4.57 5.28 5.29 5.68

П рограмма улучшения управления

Число процессоров (ядер) 1 3 5 7 9 11

Время работы программы, с 1029.85 351.99 218.83 159.60 130.71 110.29

Ускорение 1 2.93 4.71 6.45 7.88 9.34

Для исследования стратегий устойчивого развития использовалась одна из последних версий динамической социо-эколого-экономической модели региона, включающая блок активных инновационных процессов [9], которая не могла быть реализована в полном объ-

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

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

Описанный алгоритм оптимизации и некоторые другие алгоритмы (анализ чувствительности, сценарный анализ и т. п.) для социо-эколого-экономической модели реализованы в Т-системе. В вычислительных экспериментах по оценке эффективности распараллеливания программ были использованы наборы данных для Переславского [9] и Байкальского регионов.

Для Переславского региона суммарная размерность вектора состояния составила 54, а вектора управления - 56. Вычисления проводились на суперкомпьютере СКИФ МГУ «Че-бышёв» для 160 различных наборов параметров метода улучшения. Проводился запуск программы на различном числе узлов с замером времени работы. Полученные данные представлены в таблице 2.

Таблица 2

Программа улучшения управления

Число процессоров (ядер) 1 4 8 16

Время работы программы, с 10331 4184 2262 1160

Ускорение 1 2.469 4.567 8.906

Для Байкальского региона набор данных был сформирован совместно со специалистами БГУ, Бурятского и Иркутского научных центров СО РАН. Общая размерность вектора состояния, составила 3551, а вектора управления - 3588. Были проведены тестовые сценарные расчеты по анализу чувствительности целевого функционала к изменению элементов матриц прямых затрат и коэффициентов воздействия экономики на природный комплекс. Расчет в последовательном варианте занял около 6 часов. Вычислительные эксперименты с параллельной версией программы на суперкомпьютере СКИФ МГУ "Чебышёв" дали результаты, представленные в таблице 3.

Таблица 3

Программа анализа чувствительности

Число процессоров (ядер) 8 19 38 64

Время работы программы, с 3466 1483 785 520

Ускорение 6 15 28 42

Заключение

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

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

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

Литература

1. Кротов В. Ф., Гурман В. И. Методы и задачи оптимального управления. - М.: Наука, 1973.

2. Кротов В. Ф., Фельдман И. Н. Итерационные методы решения экстремальных задач // Моделирование технико-экономических процессов. М.: Изд-во МЭСИ, 1978.

3. Гурман В. И., Батурин В. А., Данилина Е. В. и др. Новые методы улучшения управляемых процессов. Новосибирск: Наука, 1987.

4. Гурман В. И. Принцип расширения в задачах управления. М: Наука, 1985.

5. Блинов А.О., Гурман В. И., Трушкова Е. А., Фраленко В. П., Программный комплекс оптимизации законов управления // Программные продукты и системы. 2009. № 2(86). С. 95-100.

6. Ухин М. Ю. Приближенный синтез оптимального управления. М.: Физматлит 2006.

7. Абрамов С. М., Есин Г. И., Загоровский И. М., Матвеев Г. А., Роганов В. А. Принципы организации отказоустойчивых параллельных вычислений для решения вычислительных задач и задач управления в Т-Системе с открытой архитектурой (ОрепТ8). // Программные системы: теория и приложения: материалы междунар. конф. Переславль-Залесский; М. Наука; Физматлит, 2006. Т. 1. С. 257-264.

8. Квоков В. Н., Трушкова Е. А., Ухин М. Ю., Метод улучшения управления на имитационной модели объекта и его приложение к задаче оптимизации маневров нештатной посадки вертолета // Вестник СГАУ. 2009. № 1 (17). С. 161-170.

9. Моделирование социо-эколого-экономической системы региона / под ред. В.И. Гурмана, Е.В. Рюминой. М.: Наука, 2001.