close

Вход

Забыли?

вход по аккаунту

?

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

код для вставкиСкачать
Кибернетика. Информационные системы и технологии
УДК 517.977.5
ЛИНЕЙНО-ВОГНУТАЯ ЗАДАЧА ОПТИМАЛЬНОГО УПРАВЛЕНИЯ С УПРАВЛЯЕМЫМИ
КОЭФФИЦИЕНТАМИ
О.Н.Кочеткова1, А.В.Бурдуковская2
1
Национальный исследовательский Иркутский государственный технический университет,
664074, г. Иркутск, ул. Лермонтова, 83.
2
Байкальский государственный университет экономики и права,
664003, г. Иркутск, ул. Ленина, 11.
Рассматривается задача оптимального управления линейной системой канонических гиперболических уравнений
с управляемыми коэффициентами. Исходная система линейна по состоянию с дополнительными условиями на
параметры задачи, при которых справедлив принцип максимума. Эта математическая модель удачно используется, например, для описания каталитических процессов в химии, и задача оптимизации в них моделирует процесс оптимизации работы химических реакторов с изменяющейся активностью катализатора. Предлагается алгоритм прямой точной релаксации. Алгоритм использует нетрадиционный вид формулы приращения целевого
функционала, что обеспечивает релаксацию без параметрической оптимизации.
Библиогр. 8 назв.
Ключевые слова: оптимальное управление; система канонических уравнений; допустимое управление; форм ула приращения; алгоритм; принцип максимума.
LINEARLY CONCAVE PROBLEM OF OPTIMAL CONTROL WITH CONTROLLABLE COEFFICIENTS
O.N. Kochetkova, A.V. Burdukovskaya
National Research Irkutsk State Technical University,
83, Lermontov St., Irkutsk, 664074
Baikal State University of Economics and Law,
11, Lenin St., Irkutsk, 664003.
The article examines the problem of optimal control of the linear system of canonical hyperbolic equations with controllable coefficients. The original system is linear in state with additional conditions on problem parameters, when the maximum principle is correct. This mathematical model is successfully used, for example, to describe catalytic processes in
chemistry, where the optimization problem models the optimization process of the operation of chemical reactors with
varying activity of the catalyst. An algorithm for direct faithful relaxation is proposed. The algorithm uses an unconventional type of the increments formula of a criterion functional that provides relaxation without parametric optimization.
8 sources.
Key words: optimal control; system of canonical equations; feasible control; increments formula; algorithm; maximum
principle.
Теория оптимального управления системами обыкновенных дифференциальных уравнений в настоящее
время достаточно высоко развита, созданы эффективные вычислительные методы поиска оптимального управления. В то же время имеется относительно малое число примеров их практического использования в управлении реальными естественно-физическими или инженерно-производственными и экономическими процессами,
т.к. описание сложных реальных процессов лишь в редких случаях может быть втиснуто в рамки систем обыкновенных дифференциальных уравнений. Как правило, состояние и управление реального процесса является
функцией не только времени, но и некоторых пространственных переменных. Отсюда наиболее естественно
описание этого процесса системой уравнений с частными производными, например [4] – [7]. Разработка здесь
общего анализа, универсальных методов крайне затруднена глубокими различиями типов уравнений, типов
начально-граничных условий, а также известной незавершенностью качественной теории при кусочнонепрерывных управлениях.
Постановка задачи
Пусть в заданной прямоугольной области P  S  T , S  s0 , s1 , T  t0 , t1 независимых переменных

s, t  P управляемый процесс определяется системой уравнений



___________________________
1
Кочеткова Ольга Николаевна, старший преподаватель кафедры математики, тел.: 89025669336.
Kochetkova Olga, Senior Lecturer of the Department of Mathematics, tel.: 89025669336.
2
Бурдуковская Анна Валерьевна, кандидат физико-математических наук, доцент, тел.: 89086412112, e-mail:
buran-baikal@mail.ru
Burdukovskaya Anna, Candidate of Physical and Mathematical sciences, Associate Professor, tel.: 89086412112, e-mail:
buran-baikal@mail.ru
8
ВЕСТНИК ИрГТУ №10 (57) 2011
Кибернетика. Информационные системы и технологии
1
u, s, t ,

2 
yt  A21u, s, t z  A22 u, s, t  y  f u, s, t 

zs  A11u, s, t z  A12 u, s, t  y  f
с начально-граничными условиями
zs0 , t   z 0 t , ys, t0   y0 s  ,
us, t U ,
здесь
(1)
(2)
(3)
U  E  компакт, не обязательно выпуклый.
n
1
2 
Вектор-функции f , f
непрерывны по своим аргументам вместе с частными производными по
удовлетворяют условиям Липшица по z, y .
z 0 t , y 0 s   заданные измеримые по t и s функции.
Под допустимыми управлениями будем понимать измеримые вектор-функции
мыми ограничениями вида (3).
Качество допустимого процесса оценивается функционалом
u  us, t  , стесненные пря-
J u    1 z s1, t , t dt   2  ys, t1 , s ds   F z, y, u, s, t dsdt  min .
T
S
z, y  и
(4)
P
z, y  для всех us, t U , s, t  P . В
задаче (1)  (4) исходная система уравнений с частными производными (1) линейна по состоянию x  z, y  с
управляемыми коэффициентами Aij  Aij u, s, t , i, j  1,2 .
Скалярные функции
1 ,  2 , F
дифференцируемы и вогнуты по
Формула приращения
Для задачи (1)  (4) выпишем формулу приращения
J u   J u~   J u    u~ H  , ~
x , u, s, t dsdt  u~ . (5)

Заметим, что u; x  x
стимые процессы и
s, t, u 
P
~  u  u, ~
x  x  x  xs, t , u~   базовый и варьируемый допуи u
H  , x, u, s, t    1 s, t , A11u, s, t z s, t   A21u, s, t  ys, t   f
+
 2  s, t , A12 u, s, t zs, t   A22 u, s, t ys, t   f
2 
u, s, t 
1
u, s, t 

 F  z , y , u , s, t  ,
F z, y, u, s, t  
,
z


T
1
T
 2  F  z , y, u , s, t  
  A21u, s, t   A22 u, s, t  
.

y

 s1   A11T u, s, t  1  A12T u, s, t  2  

В формуле приращения (5)
ций
( 2)
t
u~  0 , так как O  0; O  0,  OH  OF , OF  0
1
2
(6)
в силу вогнутости функ-
1 ,  2 , F .
Релаксационный алгоритм. Пример
Пусть на некотором k -том шаге,
k  0, 1, 2, ... , итерационного процесса u k s, t U  найденное, допу-
стимое кусочно-непрерывное управление, а




xk s, t   x s, t , u k ,  k s, t    s, t , u k  соответствующие ему
решения исходной (1)  (2) и сопряженной (5)  (6) задач. Для произвольного x построим функцию





Wk s, t , x, u~    u~ H  k , x k , u k , s, t  H  k , x k , u~, s, t  H  k , x k , u k , s, t

и в предположении о разрешимости условия максимума найдем
k
u x, s, t   arg max Wk s, t , x, u~  .
(7)
u~U
Очевидно, что
для любых
x, s, t  P .


(8)
ВЕСТНИК ИрГТУ №10 (57) 2011
9
Wk s, t , x, ux, s, t   0
Кибернетика. Информационные системы и технологии
Решение задачи (7) обусловлено обычным предположением о разрешимости условия максимума для функции Понтрягина при произвольных x [1].
Найденное из (7) управление подставим в исходную систему и решим задачу
k
k
1 k
zs  A11 u x, s, t , s, t  z  A21 u x, s, t , s, t  y  f  u x, s, t , s, t , 





 

k
k
2  k

yt  A21 u x, s, t , s, t  z  A22  u x, s, t , s, t  y  f  u x, s, t , s, t ,






0
0

z s0 , t   z t , y s, t0   y s .


Задача (9) представляет собой нелинейную по
z, y 
(9)
с возможными разрывами по состоянию правых ча-
стей, вызванных возможными разрывами по x функции u x, s, t  .
k
Предположим, что мы решили задачу (9). Обозначим ее решение 
решении построим управление


x s, t   z s, t , y s, t  . На этом
k
k
k
k
k
u~ k s, t   u  x s, t , s, t  .


k
k
~
~
Теперь формулу приращения (5) рассмотрим при u  u s, t  и u  u s, t  . Заметим, что в этом случае
k
~
x  x s, t  .
В силу неравенства
   
u~  0
k
имеем
k
J u~ k  J u k    Wk  s, t , x , u~ k dsdt ,


P
где в силу неравенства (8)
k
Wk  s, t , x , u~ k   0 .


Следовательно,
   
J u~ k  J u k .
k
Wk  s, t , x , u~ k   0 , s, t  Pk  P, mes Pk  0 , то u~ k  u k s, t , s, t  Pk , и, следова

k
~k
тельно, J u  J  u  .
 
Причем, если
 
Получена строгая релаксация без операции минимизации по скалярному параметру. Положим,
u~ k  u k 1 s, t , k  0,1,2,.., .

u~ k s, t   u k s, t , s, t  P ,
Если

k
Wk s, t , x k , u~ k  Wk  s, t , x k , u   0 .


k
k
Отсюда управление u  u s, t  удовлетворяет принципу максимума [3].
Пример
zs  u s, t , z 0, t   0,
yt  u s, t z , y s,0  0,
s, t  P  0,1 0,1,
u s, t   1.
1
1
0
0
J u    s1, t dt   ys,1ds  min .
Пример взят из [8].
10
ВЕСТНИК ИрГТУ №10 (57) 2011
то
x s, t   x k s, t ,
k
Кибернетика. Информационные системы и технологии
Информация о примере
H   1u   2 u  z  ps, t   u,
ps, t    1   2  z ,
 s1   2   u; 1 1, t   1,
 t2   0,  2  s,1  1;  2  s, t   1,

1
1
s
s, t   1   u s, t ds   u  , t d

0
0
s,1  z s, t .
0
Пусть начальным приближением является управление u s, t   0,75 .
u z , s, t   sign 
В этом случае
1
z 0 s, t   0,75  s, y 0 s, t   0,56  s  t.
 01 s, t   0,25  0,75s.
Решим задачу
u z, s, t   arg max  01 s, t   z s, t u
0
u 1
.
Очевидно, что
u z, s, t   sign 0,25  0,75s  z  .
0
Теперь решим уравнение
zs  sign 0,25  0,75s  z ,
z 0, t   0.
Предположим, что
zs, t   0,25  0,75s .
Тогда имеем уравнение zs  1 с начальным условием z 0, t   0 .
Его решение удовлетворяет нашему предположению. Если же
zs, t   0,25  0,75s , то решение
z s, t   s, zs  1, z 0, t   0 не удовлетворяет предположению. Следовательно,
0
0
z s, t   s, u  1, u~0 s, t   1 .
0
 
На этом управлении (не равном u s, t ) в силу алгоритма обеспечивается строгая релаксация без параметрической оптимизации, что и подтверждается расчетами
~0 s, t   1
Далее при u s, t   u
 
 
J u~0  1,5  J u 0  0,47 .
1
z1 s, t    s,
y1 s, t   s  t ,
 s1 s, t   s  2,


u z , s, t   sign  s1 s, t   z s, t  ,
1
u z , s, t   signs  2  z .
1
z s, t  : zs  u z, s, t   signs  2  z , z 0, t   0 .
Рассмотрим случаи, когда s  2  z  0 . Решением уравнения
1
1
zs  1; z0, t   0
ВЕСТНИК ИрГТУ №10 (57) 2011
11
Кибернетика. Информационные системы и технологии
является функция
z s, t   s ,
1
для которой справедливо
1
s  2  z  s  2  s  2s  2  0 (кроме s  1 ).
Предположим теперь, что 2  s  z  0 . Найдем решение уравнения
zs  1, z0, t   0 .
Для этого решения
z s, t   s
1
выполнено неравенство
s  2  s  2s  2  0 ,
которое противоречит предположению. Таким образом,
z s, t   s ,
1
определим
u~1 s, t   signs  2  s   1; u~1  u1 .
Библиографический список
1. Бурдуковская А.В., Васильев О.В. Оптимизация систем канонических гиперболических уравнений с гладкими ограниченными управлениями // Журнал вычислительной математики и математической физики. 2000. Т. 40, №1. С. 43–53.
2. Бурдуковская А.В., Васильев О.В. Об алгоритмах оптимизации в системах канонических гиперболических уравнений с
частными производными. Серия: Оптимизация и управление. Иркутск: Изд-во ИГУ, 1998. Вып.2. 55 с.
3. Кочеткова О.Н., Бурдуковская А.В. Алгоритм оптимизации системы канонических гиперболических уравнений с гладкими
ограниченными управлениями // Вестник ИрГТУ. 2010. № 7. С. 295–301.
4. Быков В.И., Яблонский Г.С., Слинько М.Г. Применение принципа максимума для оптимизации квазистационарных каталитических процессов с изменяющейся активностью катализатора // Proc. IFIP Techn. Conf. on Optimizat. Techn. Новосибирск,
1974. С. 11–16.
5. Васильев О.В. Принцип максимума Л.С. Понтрягина в теории оптимальных систем с распределенными параметрами //
Прикл. матем. Новосибирск: Наука, 1978. С. 109–138.
6. Васильев О.В., Срочко В.А., Терлецкий В.А. Методы оптимизации и их приложения (Оптимальное управление). Новосибирск: Наука, 1990.
7. Забелло Л.Е. Об условиях оптимальности в нелинейных инерционных управляемых системах с запаздыванием // Дифференц. ур-ния. 1990. Т. 26, № 8. С. 1309–1315.
8. Срочко В.А. Вариационный принцип максимума и методы линеаризации в задачах оптимального управления. Иркутск:
Изд-во Иркутского ун-та, 1989. 154 с.
УДК 681.5.003.23 (517.977.58)
МЕТОДИКА ИССЛЕДОВАНИЯ АЛГОРИТМОВ АВТОМАТИЧЕСКОЙ ПАРАМЕТРИЧЕСКОЙ
ОПТИМИЗАЦИИ СИСТЕМ С ИНТЕГРАЛЬНОЙ ШИРОТНО-ИМПУЛЬСНОЙ МОДУЛЯЦИЕЙ
Н.Н.Куцый1, Е.А.Осипова2
Национальный исследовательский Иркутский государственный технический университет,
664074, г. Иркутск, ул. Лермонтова, 83.
Рассматривается задача параметрической оптимизации систем с интегральной широтно-импульсной модуляцией. Приведены результаты исследования разработанного алгоритма параметрической оптимизации. Решается
вопрос о выборе числа настраиваемых параметров.
Ил. 4. Библиогр. 9 назв.
Ключевые слова: алгоритм автоматической параметрической оптимизации; беспоисковая самонастраивающаяся система; интегральная широтно-импульсная модуляция; чувствительность.
___________________________
1
Куцый Николай Николаевич, доктор технических наук, профессор кафедры автоматизированных систем, тел.: 89149178520,
e-mail: kucyinn@mail.ru
Kutsyi Nikolai, Doctor of technical sciences, Professor of the Department of Automated Systems, tel.: 89149178520, e-mail:
kucyinn@mail.ru
2
Осипова Елизавета Алексеевна, аспирант, старший преподаватель кафедры автоматизированных систем, тел.: 89501204839,
e-mail: olisa252@mail.ru
Osipova Elizaveta, Postgraduate, Senior Lecturer of the Department of Automated Systems, tel.: 89501204839, e-mail: olisa252@mail.ru
12
ВЕСТНИК ИрГТУ №10 (57) 2011
Документ
Категория
Без категории
Просмотров
6
Размер файла
3 072 Кб
Теги
оптимальное, линейной, управляемыми, управления, коэффициента, вогнутая, задачи
1/--страниц
Пожаловаться на содержимое документа