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

ISSN 2079-3316 ПРОГРАММНЫЕ СИСТЕМЫ: ТЕОРИЯ И ПРИЛОЖЕНИЯ № 4(22), 2014, с. 243-256 УДК 517.97

Е. А. Трушкова, О. В. Фесько Модели управления параметризованной структуры

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

Ключевые слова и фразы: оптимальное управление, параметризация управления, дискретно-непрерывная модель, достаточные условия оптимальности.

Введение

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

Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (проект 14-31-50824-мол_нр).

© Е. А. Трушкова, О. В. Фесько, 2014

© Институт проблем управления имени В.А. Трапезникова РАН, 2014 © Восточно-КитАйский педагогический университет, 2014 © Программные системы: теория и приложения, 2014

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

В статье рассматриваются дискретизованные модели управления с постоянным и линейным законом изменения на шагах, которые приводят к двухуровневым дискретно-непрерывным моделям. Такие системы являются частным случаем более общего класса моделей, ориентированных на представление и исследование разнообразных систем неоднородной структуры [5-7].

1. Постановка задачи

Рассматриваемая модель управляемого процесса имеет вид:

(1) x(t) = f (t,x(t),u(t)), х (tr) = Xj, t e [tr,tF] ,

(2) u(t) e D„ (t) = [u(t) | u(t) < u(t) < u(t)} ,

(3) I = F (x (tF)) ^ inf e D.

Здесь x(t) — n-мерный вектор состояний, xi(t), i = 1,п — кусочно-гладкие, u(t) — р-мерный вектор кусочно-непрерывных управлений, ■u(t) и u(t) — заданные кусочно-непрерывные функции в R, t — время, D — множество допустимых элементов т = (х (t) ,и (t)), удовлетворяющих условиям (1)-(2).

Рассмотрим множество кусочных управлений U (Ь (t,w) ,т), где b(t, w) — непрерывная базовая вектор-функция, зависящая от векторных параметров w = (wi,...,w^) , а т e Z+ —число моментов переключений управления при разбиении временного отрезка ti = т0 < т\ <т2 < ... < тт+1 = tF. Размерность каждого параметра wi, г = 1, к зависит от то.

Выбор базовой функции в виде постоянной вектор-функции b (t, w) = (wi, W2,..., wv)T, к = p приводит к классу кусочно-постоянных управлений, а ее задание в виде линейной вектор-функции b (t,w) = (wit + W2,wst + W4,... ,w2p~it + W2P)T, к = 2p образует класс кусочно-линейных управлений.

Вместо исходной задачи оптимального управления на множестве допустимых управлений D„ (t) будем рассматривать суженную

задачу (1)—(3) на множестве кусочных управлений U (Ь (t, w), т), которая, очевидно, вкладывается в множество допустимых управлений исходной задачи.

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

Суженная исходная задача (1)—(3) на множестве кусочных управлений сводится к задаче условной конечномерной минимизации дифференцируемой функции многих переменных F (х (tp)) = G (w,t), где w — совокупность всех параметров управления, Ti,i = 1,т — моменты переключений управления.

2. Дискретизация непрерывной системы

В работе [8] описаны преобразования, приводящие суженную непрерывную задачу оптимального управления к ее дискретному аналогу.

Задачу (1)—(3) на множестве управлений U (Ь (t,w) ,т) можно представить в виде задачи оптимального управления для дискретно-непрерывной системы

z (I + 1) = д (l,z(l),w(l)), z(0) = z0 = xr, z € Mn,

(4) w(l) € W = {w(l) | w < w(l) < w] С Rfc, l

F (z(m + 1)) ^ min,

где m — число точек переключений управления при разбиении временного отрезка ti = т0 < т1 < т2 < ... < тт+1 = tp, а функция д (I, z(l), w(l)) — решение задачи Коши

(5) i(0 = f (€,х (О ,Ь (€,W (I))) , е € [п,п+1] , х (п) = z (I) ,

взятое в точке £ = r;+i. Непрерывная система (5) нижнего уровня в данном случае не содержит управлений, в то время как дискретная система верхнего уровня содержит к управлений. Заметим, что непрерывная задача (1)—(3) на множестве управлений U (Ь (t,w) ,т) и дискретная задача (4) эквивалентны.

По аналогии с достаточными условиями оптимальности Кротова вводится в рассмотрение функция ¡р(1, г) и строятся следующие конструкции:

К(1, г, ю) = <р(1 + 1, /(I, г, ю)) — <р(1, г), С(г) = ^(г) + <р(т + 1, г) — ^(0, г0).

т

Положим Ь = О (г (т +1)) — ^^ Я(1, г, ад).

1=0

Теорема 1. Пусть {z*,w*} —допустимая пара в задаче (4) и

z*,w*}

существует р>, такие что R(l, z*,w*) = max R(l,z,w) и G(z*)

wEW,z

minG(z). Тогда пара {z*,w*} —решение задачи (4).

Функцию ¡p(l, z) предлагается искать из условий (схема Кротова-Беллмана)

<р (I, z) = max <р (I + l,g (I, z, w)),

(6) wew

tp (m + 1,z) = -F (z), I e{0,...,m}.

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

3. Нелокальный метод улучшения на основе минимаксного принципа

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

Рассмотрим дискретную задачу оптимального управления

z(l + 1) = g (l, z(l), w(l)), l e L = {00,... pm - 1}, (7) z(o) = Z0, Z e R", w(l) e W с Rk, F (z(m)) ^ min .

Задача улучшения ставится следующим образом: пусть известен допустимый элемент задачи (7) (zl(1), wl(l)), требуется найти допустимый элемент задачи (7) (zll(l), wll(l)), такой что F (zll(m)) < F (zI(m)).

Для любой функции p(l, z) запишем конструкции, фигурирующие в достаточных условиях оптимальности:

т—1

L(z,w;p) = G(z(m)) — ^^ R(l,z,w),

1=0

где

G(z) = F(z) + <p(m, z) — p(0, zo),

R(I, z, w) = ip(l + 1, g(l, z, w)) — tp(I, z).

Заметим, что L(z,w;tp) = F(z(m)) на любой допустимой паре (z, w).

Одна итерация известного метода нелокального улучшения [9,10] заключается в следующем:

0) Имеем начальный допустимый процесс (zl(l),wl(l)).

1) Ищем p0(I, z) из соотношений:

(8) R( I, zI(l),wI (l))=minR( l,z,w1(l)), le{0,...,m — 1},

Z

(9) G( z1(m)) = maxG(z).

z

2) Строим функцию:

w(I, z) = arg max R(I, z, w; p0).

w£.~W

3) Решая систему

z(l + 1)=g (I, z(l), w(I, z(l))), z(0) = z0, находим улучшенный допустимый процесс (zll(l), w11 (I)).

Далее будем называть соотношения (8), (9) соотношениями метода нелокального улучшения, а соответствующую функцию <p0(l, z) — разрешающей функцией.

Здесь предлагается искать разрешающую функцию p°(l, z) из условия, чтобы функционал L(z(l),wl(l)) не зависел от z(l) (это специальный случай максимума по ( ) )

(10) R(l, z(l),wl(l)) = 0, le{0,...,m — 1},

(11) G(z) = —tp(0, zo).

Тогда соотношения (8), (9) перепишутся в более простом виде, а именно, в виде линейных рекуррентных соотношений

(12) p(l, z) = tp (I + 1, g(l,z,w1(l))) , le{0,...,m — 1},

(13) p(m, z) = —F(z).

Одна итерация построенного метода нелокального улучшения заключается в следующем:

0) Имеем начальный допустимый процесс (zl(l),wl(l)).

1) Ищем p0(I, z) из соотношений (12), (13).

2) Строим функцию:

(14) w(I, z) = arg max p0(I + 1, g(I, z, w)) — p0(I, z)).

3) Решая систему

z(l + 1)=g (I, z(l), ü(I, z(l))), z(0) = z0, находим улучшенный допустимый процесс (zll(l), w11 (I)).

Теорема 2. Вышеописанный метод нелокального улучшения гарантирует неухудшение функционала, т. е. выполнение неравенства F (zl(m)) — F (zu(m)) > 0. Если же для функции ip0, разрешающей соотношения (12), (13), условие

R(l, zI(l),wI(l); p0) = max R(I, zI(l),w; p0),

не выполняется хотя бы при одном значении I = т, то будет верно строгое неравенство

F (z\m)) —F (zll(m)) > 0.

Доказательство теоремы приведено в [8].

Заметим, что фактически табулировать и запоминать функцию w (I, z(l)) не требуется; она вычисляется при прямом счете; достаточно иметь только представление p0.

Соотношения рассмотренного метода (12), (13) для поиска разрешающей функции весьма удобны для программной реализации соответствующего алгоритма. При этом расчет по формулам (12)— (14) может производиться с использованием аппроксимации функции p0(l, z) при каждом фиксированном I (например, с помощью полиномов по методу наименьших квадратов), что даст возможность проводить анализ получаемой аппроксимации функции p0(l, z) в виде аналитических выражений. Соответствующий алгоритм состоит из следующих шагов.

(1) Разрешается система (7) при заданных wl(l) и -г1(0) = z0, определяется ml = (zl(l),wl(l)).

(2) Задаются сетки узлов в пространстве состояний и управлений: Sz(I) = {zu(l) — Az1, zu(l), zu(l) + Az1} x ...

... x {z1n(l) — Azn, z1n(l), z1n(l) + Azn) ,

Sw = {w— ,w— + h^, .. .,w+ — h^, w+} x ...

.. i к к I 7 к к ik fcl

... x |w—, w— + hw, ..., w+ — hw, w+],

где Az%, i = 1, n, h?w, j — 1,k, некоторые числа (регуляторы метода).

(3) На сетке узлов Sz (t) с помощью МНК ищется аппроксимация в виде полинома функции p0(I, z) из соотношений (12), (13).

(4) Разрешается при заданном z(0) = z0 система

z(l + 1) = д I I, z(l), arg max p0 (I + 1, g(l, z(l), w))

у wesw

определяется m11 = 11 ( ) , w11 ( ) .

(5) Если F(z11(m)) — F(z1(m)) < e, то в качестве текущего улучшенного управления выбирается w1 (1) = wu(l), и осуществляется переход на шаг 2, иначе итерации завершаются.

Существенно снизить время работы программы, реализующей описанный алгоритм, можно при решении задач управления некоторыми классами непрерывных систем за счет поиска разрешающей функции p0(I, z) из задачи Коши для уравнения в частных производных, являющейся аналогом соотношений (12), (13) [12,13].

4. Вычислительные алгоритмы

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

Первый алгоритм предназначен для поиска приближенно оптимального управления с помощью параметрических множеств управлений простой структуры и основан на сведении задачи (1)—(3) на множестве управлений U (Ь (t, w), m) к задаче многомерной оптимизации функции G (w,T, ß0, ßi, ß2), где моменты переключений т выступают тоже как параметры, ß0, ßi, ß2 > 0 — весовые коэффициенты для задачи с фазовыми ограничениями [14]. При этом для решения задачи Коши (1) использовался численный метод интегрирования систем обыкновенных дифференциальных уравнений, а сама задача оптимизации решалась комбинацией методов 1-го и 2-го порядков

,

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

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

Блок численного решения задачи Коши (РК/РКФ45)

Щ = 5 (х М)ДМ = О

Улучшение управления

Подсчет минимизируемого функционала

С(ш,Т,/30>)6ь/32)

Выбор штрафных коэффициентов

условия не выполнены

Комбинация градиентного спуска и метода Ньютона-Рафсона, проверка условий останова

условия выполнены

м> ,т

Рис. 1. Вычислительная схема улучшения управления

Во втором алгоритме осуществляется поиск начального приближения к решению во всей допустимой области значений ад на основе достаточных условий оптимальности за счет разрешения рекурсивной цепочки (6) (соотношения Кротова-Беллмана) относительно функции р(1, г) на сетке управлений |ад, ад + к,..., ад — к, ад |

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

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

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

Описание задачи, задание входных параметров

Рис. 2. Вычислительная схема поиска начального приближения и улучшения управления

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

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

5. Вычислительные эксперименты и приложения

Рассматривается задача динамической оптимизации подпитываемого биореактора для производства пенициллина путем анаэробной ферментации глюкозы [15-19]. Задача оптимального управления заключается в максимизации общего количества производимого пенициллина, используя в качестве управления интенсивность подачи субстрата (скорость подпитки). Необходимо максимизировать функционал

Р (х(Ьр)) = х2 (Ьр) х4 (Ър) ^ тах

в соответствии с динамикой системы, описываемой набором обыкновенных дифференциальных уравнений

и • Ж1 . . и • Х2

х 1(г) = п1х1 - —-—, х2^) = П2Х1 - 0.01Ж2 - ТТТТ;—, 500^4 500^4

. Ь,2Х\ 0.029x1^3 и / хз \

Х3(1' = - 0.47 - 1.2 - 0.0001 + х3 + щ V 500/ ,

где

х4(Ь) = и/500,

0.11ж3 0.0055ж3

Пл =_____, П2 =

0.006ж1 + 2 0.0001 +ж3(1 + 10ж3)'

Здесь хл,х2 из — концентрации биомассы, пенициллина и субстрата (в г/л) соответственно, а 14 — объем культуры (в л). Вектор начальных состояний: ж(0) = (1.5, 0, 0, 7)т.

Конечный продукт предназначен для потребления человеком. Поэтому для избежания вредных последствий его производство должно подчиняться некоторым условиям. По этой причине введены ограничения: 0 < хл < 40, 0 < х3 < 25, 0 < х4 < 10. Управление также ограничено: 0 < и < 50 (в л/ч); общее время протекания процесса фиксировано: £ р = 132 ч.

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

Ниже приведены результаты расчета, касающиеся той подобласти, в которой найдено оптимальное решение задачи. Начальное управление ад = (12.5,12.5, 37.5,12.5,12.5,12.5,12.5)т и моменты переключений интенсивности подпитки т = (5,10,15, 20, 25, 50)т.

В ходе работы программы удалось добиться улучшения управления ад = (14.54,11.24, 37.05,13.89,13.77,12.75, 9.41)т, сдвига моментов переключений т = (4.07,11.12,15.33,19.04, 24.84, 54.9)т и выполнения фазовых ограничений, доставив при этом значение функционалу Р (х (Ьр)) = 89.1, что является наилучшим результатом среди вышеперечисленных работ. В частности, в работе [18] было получено значение Р (х (Ьр)) = 88.013 с двадцатью точками переключений кусочно-постоянного управления. Значение траекторий в конце процесса: ж (Ьр) = (28.65, 8.85, 0.001,10.07)т.

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

6. Заключение

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

Рис. 3. Оптимальные состояния (слева) и управление (справа) для задачи по производству пенициллина

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

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

[1] Л. С. Понтрягин, В. Г. Болтянский, Р. В. Гамкрелидзе, Е. Ф. Мищенко. Математическая теория оптимальных процессов. М.: Физматгиз, 1961. ^ 243.

[2] Р. Беллман. Динамическое программирование. М.: ИЛ, 1960. ^ 243.

[3] В. Ф. Кротов. Методы решения вариационных задач на основе достаточных условий абсолютного минимума // Автомат. и телемех. 1962, № 12, с. 1571— 1583. | 243.

[4] В. Ф. Кротов. Достаточные условия оптимальности для дискретных управляемых систем // ДАН СССР, 1967. Т. 172, № 1, с. 18-21. | 243.

[5] В. И. Гурман, И. В. Расина. Сложные процессы.. — Новосибирск: Наука, 1990, с. 84—94. | 244.

[6] В. И. Гурман, И. В. Расина. Дискретно-непрерывные представления импульсных решений управляемых систем // Автомат. и телемех. 2012, № 8, с. 16-29. | 244.

[7] И. В. Расина. Дискретно-непрерывные .модели и оптимизация управляемых процессов // Программные системы: теория и приложения, 2011, № 5(9), с. 49-72. | 244.

[8] Е. А. Трушкова. Алгоритмы глобального поиска оптимального управления // Автомат. и телемех. 2011, № 6, с. 151-159. f 245, 248.

[9] В. Ф. Кротов, И. Н. Фельдман. Итерационный метод решения задач оптимального управления // Изв. АН СССР. Техн. киберн. 1983, № 2, с. 160-168. f 246, 247.

[10] V. F. Krotov. Global methods in optimal control theory. N.Y.: Marcel Dekker, 1996. f 246, 247.

[11] В. Ф. Кротов. Управление квантовыми системами и некоторые идеи теории оптимального управления // Автомат. и телемех. 2009, № 3, с. 1523. f 246.

[12] Е. А. Трушкова. Об одном классе задач оптимального управления для квантовых систем // Автомат. и телемех. 2013, № 1, с. 35-46. f 249.

[13] Е. А. Трушкова. Метод глобального улучшения для гамильтоновых систем с управляемыми коэффициентами // Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика, 2013. Т. 13, № 1(2), с. 95-99. f 249.

[14] В. Н. Квоков, Е. А. Трушкова, М. Ю. Ухин. Метод улучшения управления на имитационной модели объекта и его приложение к задаче оптимизации маневров нештатной посадки вертолета // Вестник Самарского государственного аэрокосмического университета имени академика С. П. Королева, 2009, № 1, с. 161-169. f 249.

[15] J. R. Banga, R. Irizarry-Rivera, W. D. Seider. ¡Stochastic optimization for optimal and model-predictive control // Computers and Chemical Engineering, 1998. Vol. 22, no. 4-5, p. 603-612. f 252.

[16] J. E. Cuthrell, L. T. Biegler. Simultaneous optimization and solution methods for batch reactor control profiles // Computers and Chemical Engineering, 1989. Vol. 13, no. 1-2, p. 49-62. f 252.

[17] T. Hirmajer, E. Balsa-Canto, J. R. Banga. DOTcvp: Dynamic Optimization Toolbox with Control Vector Parameterization approach for handling continuous and mixed-integer dynamic optimization problems, 2010. — 66 p. f 252.

[18] J. A. E. Larrosa. New heuristics for global optimization of complex bioprocesses, (2008), PhD thesis. 235 p. f 252, 253.

[19] R. Luus. Application of dynamic programming to differential-algebraic process systems // Computers and Chemical Engineering, 1993. Vol. 17, no. 4, p. 373377. f 252.

Рекомендовал к публикации д.т.н. В. И. Гурман

Об авторах:

Екатерина Александровна Трушкова

д.ф.-м.н., в.н.с. лаб. 45 Института проблем управления им. В.А. Трапезникова РАН; область научных интересов: теория оптимального управления, методы и алгоритмы приближенной оптимизации управления, параллельные вычисления.

e-mail: [email protected]

Олесь Владимирович Фесько к.т.н., н.с. ИЦСА Института программных систем им. А.К. Айламазяна РАН; East China Normal University; область научных интересов: приближенные методы оптимального управления, достаточные условия оптимальности, дискретно-непрерывные модели

e-mail: [email protected]

Образец ссылки на эту публикацию:

Е. А. Трушкова, О. В. Фесько. Модели управления параметризованной структуры // Программные системы: теория и приложения: электрон. научн. журн. 2014. Т. 5, № 4(22), с. 243-256.

URL

http://psta.psiras.ru/read/psta2014_4_243-256.pdf

Ekaterina Trushkova, Oles Fesko. Parameterized control models. Abstract. In this article discretized control systems with parametric structure are considered. They lead to descrete-continuous models. Parallel algorithms for finding an approximate optimal control are developed on the basis of these models. An example is given.

Key Words and Phrases: optimal control, control parameterization, descrete-continuous models.