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

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

CC BY
248
17
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ЗАДАЧА ОПТИМАЛЬНОГО УПРАВЛЕНИЯ / ПРИНЦИП МАКСИМУМА / ДОСТАТОЧНЫЕ УСЛОВИЯ ОПТИМАЛЬНОСТИ / OPTIMAL CONTROL PROBLEM / MAXIMUM PRINCIPLE / SUFFICIENT OPTIMALITY CONDITIONS

Аннотация научной статьи по математике, автор научной работы — Срочко Владимир Андреевич

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

The Simplest Nonconvex Control Problem. The Maximum Principle and Sufficient Optimality Conditions

The optimal control problem with linear phase system and linear-quadratic functional is considered. The transition from the maximum principle to sufficient optimality conditions is fulfilled with the help of the notion of strongly extremal control. It means that in the problem of maximization of Pontryagin’s function phase or conjugate trajectory should be replaced with any admissible trajectory. Sufficient conditions give opportunity to obtain explicit expressions for extremal values of auxiliary problems contained in these conditions. Results are presented in the form of inequalities and equalities for a function with one variable with respect to time segment. A special situation is implemented in the analysis of the combined control with interior and boundary segments with respect to the constraint. At the point of connection of these segments there is a non-standard condition of maximum type. A positive factor is dual nature of obtained results: it is a pair of symmetrical relations, which operate independently. Their origin is connected with two types of strongly extremal controls with respect to phase or conjugate variables.

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

Онлайн-доступ к журналу: http://isu.ru/izvest.ia.

Серия «Математика»

2017. Т. 19. С. 184-194

УДК 517.97 MSG 49J15

DOI https://doi.org/10.26516/1997-7670.2017.19.184 Простейшая невыпуклая задача управления. Принцип максимума и достаточные условия оптимальности

В. А. Срочко

Иркутский государственный университет

Аннотация. Рассматривается задача, оптимального управления относительно линейной фазовой системы и линейно-квадратичного по паре «состояние - управление» функционала.. Переход от принципа, максимума, к достаточным условиям оптимальности проводится на. основе понятия сильно экстремального управления. Это значит, что в задаче на. максимум функции Понтрягина. необходимо заменить фазовую или сопряженную траекторию на. произвольную допустимую траекторию. Конструктивный характер достаточным условиям обеспечивает возможность получения явных выражений для экстремальных значений вспомогательных задач, фигурирующих в этих условиях. Результаты представляются в форме неравенств и равенств для функций одной переменной на. промежутке времени, что вполне соответствует принципу максимума.. Особая ситуация реализуется при анализе комбинированного управления с внутренними и граничными участками относительно ограничения. В точке сопряжения этих участков возникает нестандартное условие типа, равенства, как следствие состыковки двух неравенств.

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

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

В рамках теории достаточных условий оптимальности в задачах управления необходимо выделить, как минимум, основополагающие монографии [5; 7] и отметить, например, развивающие это направление в качественном и конструктивном плане работы [2; 4; 6]. Одна из проблем

Введение

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

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

1. Постановка задачи. Условия оптимальности

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

Ф(и) ->• шах, и € V, (1.1)

в которой

¿1

ф(и) = {с,х(и)) + J- ^

¿0

х = + &(£)«, ж(£о) = ж0, V = М-) € РС([*0,*1]) : К*)| < 1, t € [*0,*1] }•

Здесь € К - управление, ж(£) € Кп - фазовое состояние,

_ множество кусочно-непрерывных функций, с,х°, 7 > 0 -параметры задачи.

Предположим, что вектор-функции а(£),6(£) и матричная функция непрерывны на [¿0,^1]-

Приведем соотношения принципа максимума (ПМ) для задачи (1.1). Введем функцию Понтрягина

Н(1р, х, и, ¿) = (0, А({)х) + (0, &(£))« + (а(£), х)и — ^уи2 и сопряженную систему

ф = -Ат(г)гр - а{Ь)и, 1р(Ь) = с.

Определим стационарную точку функции Н по переменной и

ы(ф,х,1) = -((а(£),ж) + (Ь(1),ф) 7 ^

Тогда ^/-максимизирующее управление при ограничении |и| ^ 1 выражается по формуле

{ги(ф,хЛ), если \ъи(ф, х, ¿)| ^ 1, sign и)(ф,х,£), если \и)(ф, х, ¿)| ^ 1.

Пусть и(Ь) - допустимое управление (и € V), х(Ь,ь), ф(Ь,ь), £ € € [¿0)^1] - соответствующие траектории. Обозначим

Тогда принцип максимума для управления и(Ь) (ПМ(м)) представляется следующим соотношением: € [¿0,^1]

Гй(£), если|(а(*),ж(г,и))+ <7,

и\Ч = \ _

[ sign■u(t), если |(а(£),ж(£,«)) + (&(£), ф^, и)}\ ^7.

Следуя [1; 8], перейдем на уровень достаточных условий оптимальности. С этой целью достаточно в предыдущей формуле последовательно провести замены

ж(£, и) —> ж(£, у), V € V или ф(1,, и) —> ф(1,, у), V € У

(фиксированная траектория —>■ любая траектория, экстремальное управление —>■ сильно экстремальное управление).

В результате получаем достаточные условия оптимальности в начальном варианте.

Первое условие: \/у € V, £ € [¿0,^1]

С если|(а(*),ж(*,и))+ <7,

■и(£) = < _ (1.2)

[ гл(^), если |(а(£),ж(£,-и)) + (&(£),«))| ^ 7-

Второе условие: \/г> € V, £ € [¿0,^1]]

[«(*), если|(а(г),ж(г,и)) + (&(£),< 7, , ч и(Ь) = < _ (1.3)

[ sign■u(t), если |(а(£), ж(£, и)) + (&(£), ^ 7-

Представленные соотношения определены на множествах фазовых либо сопряженных траекторий, что существенно затрудняет их проверочную реализацию. В рамках данной задачи (1.1) имеется возможность конструктивно использовать экстремальные варианты полученных неравенств, что приводит к поточечным по времени £ € [¿0^1] условиям оптимальности, характерным для принципа максимума.

2. Преобразование условий

Найдем экстремальные значения скалярных произведений (a(t),x(t,v)),(b(t),ip(t,v)) для момента t € [íojíi] на множестве V.

Лемма 1. Пусть вектор-функция p(r,t), т € [to,t], t € [to,íi] является решением задачи Коши

pT(r,t) = —АТ(т)р(т, t), p{t,t) = a(t).

Тогда справедливо представление

t

(a(t),x(t,v)) = j(Ь(т),р(т, t))v{r)d,T + (p(t0,t),x°). to

Для доказательства необходимо найти производную

d

— (р(т, t),x(r,v))

с последующим интегрированием по г € [¿сь^]-

Следствие 1. Экстремальные значения для {a(t),x(t,v)} выражаются по формуле

t

min(max)(a(t),x(t,v)) = (p(to,t),x°) - (+) / \(Ь(т),р(т,Щс1т.

vëv vëv J

t0

Лемма 2. Пусть вектор-функция q(r,t), t € [t, t\], t € [to,íi] является решением задачи Коши

qr(T,t) = A(T)q(T,t), q(t,t)=b(t).

Тогда справедливо представление

ti

(,b(t),iP(t, v)) = J(а(т), q(r, t))v{r)dT + {q(h,t),c). t

Для доказательства необходимо найти производную

d

— (q(r, t),ip(r, v))

с последующим интегрированием по г € [t,t 1].

188 в. а. срочко

Следствие 2. Экстремальные значения для {b(t),ip(t,v)} выражаются по формуле

ti

min(max)(b(í),^(í,u)) = (q(ti,t),c) - (+) / | <a(r), q(r, t))\dr.

vëV v&V J

t

Приступим теперь к преобразованию достаточных условий оптимальности. Раскроем модульное неравенство в первом условии для некоторого момента t € [to,¿i]

-7 < (a(t),x(t, v)) + (b(t),tp(t, и)) ^7, v eV. (2.1)

Перейдем к эквивалентной экстремальной форме

таx(a(t),x(t,v)) + (b(t),ip(t,u)) ^ 7, v&v

min(a(t),x(t,v)) + (b(t),ip(t,u)) ^ -7.

vdv

Учитывая следствие 1, получаем

t

(b(t),ip(t,u)) + (p(to,t),x°) + J \ (Ь(т),р(т, t))\dr < 7, (2.2)

to t

(b(t),tp(t,u)) + (p(to,t),x°) - J \ (Ь(т),р(т, t))\dr > -r

to

Объединим эти неравенства одним выражением

t

I (b(t),ij(t,u)) + (p(t0,t),x°)\ + J \ {Ь(т),р(т, t))\dr < 7. (2.3)

to

Таким образом, полученное неравенство для функции одной переменной t эквивалентно соотношению (2.1).

Аналогичным образом обрабатывается модульное неравенство во втором условии с учетом следствия 2. В итоге приходим к дополнительному к (2.3) неравенству в момент t € [Í0)¿i]

ti

\(a(t),x(t,u)) + (q(t1,t),c)\+ J\(a(T),q(T,t))\dT ^ 7. (2.4)

t

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Сформулируем первый результат об оптимальности «внутреннего» управления u(t) = u(t), \u(t)\ ^ 1, t € [to,¿i] которое определяется условием стационарности функции Понтрягина

Hu(ip(t,u),x(t,u),u(t),t) = 0, t € [ío,il] ^

7 ^

Теорема 1. Пусть для управления и(Ь) = £ € [¿0,^1] выполнено неравенство (2.3) или (2-4) для всех £ € [¿0,^1]- Тогда это управление является оптимальным в задаче (1.1).

Пример 1.

1

Г г / \ 1 „ -

dt max,

Ф(и)=хг(1) + j (xi(t) +x2(t)^ju(t) - -~fu2(t) 0

x\ = X2, ¿2 = u, a;i(0) = —2, жг(0) = 1, |«(i)| ^ 1, i e [0,1]. В данном случае

1 2

H = Ф1Х2 + 02"U + X\и + X2и — -7и ,

01 = -и, 02 = -01 - и, 01 (1) = 1, 02 (1) = 0. Рассмотрим управление и = 0. Ему соответствуют траектории

X\(t, 0)=t-2, x2(t, 0) = 1, 0i (t, 0) = 1, 02(t, 0) = 1 - t.

Условие стационарности выполняется V7 > 0: Hu[t, 0] = 0, t € [0,1], т. е. управление и = 0 удовлетворяет принципу максимума.

Проверим условие (2.3). Вспомогательная вектор-функция p(r,t): Pi = 1, P2 = t — т + 1. Следовательно,

(6,0(i,O)) = l-t, (p(0,t),x°)=t-l, (b,p(T,t)) =t-T + 1.

Неравенство (2.3) принимает вид

t2

- + ¿^7, te [0,1].

Таким образом, для 7^3/2 неравенство выполнено, т. е. управление и = 0 является оптимальным.

К аналогичному результату приводит второе условие (2.4). Заключительное неравенство имеет вид

Ç-2î + |<7, ÎG[0,1] 7>f-

Рассмотрим далее вопрос об оптимальности граничных управлений и = ±1 на основе соотношений (1.2), (1.3). Для управления и = 1 достаточные условия представляются следующим образом

mm(a(t),x(t,v)) + (Ь(1),ф(1,1))

vdv

(а(1),х(1,1))+тт(Ь(1),ф(1,ь)) ^ 7,

1>€ V

В итоговой форме получаем неравенства для £ € [¿0,^1]

+ {р(1о,1),х0} - I Жт),р(т,тб,т>ъ (2.5) ¿0 ¿1

(а(1),х(1,1)) + ШъЬ),с) - I\(а(т),д(т,Щс1т> 7. (2.6)

г

Сформулируем результат.

Теорема 2. Пусть для управления = 1, £ € [¿0,^1] выполнено неравенство (2.5) или (2.6) для всех £ € [¿0,^1]- Тогда это управление является оптимальным в задаче (1.1).

Установим связь полученных условий с принципом максимума для управления и = 1 (ПМ(1)).

Предположим, что выполнено неравенство

(Ь(т),р(т,$) < 0, г € « € [io.ii]• (2.7)

Тогда на основании леммы 1 и следствия 1 заключаем, что

тт(й(4 ж(£, -и)) = (а(£), ж(£, 1)), £ € [¿о, ¿1]-г)€У

При этом условие (2.5) представляется в виде

Ш,х(1,1)) + (Ь(1),ф(1,1)) ^ 7, I € [to.ii].

Получили ПМ(1) как достаточное условие оптимальности в предположении (2.7).

Аналогично, неравенство

(.а(т),д(т,¿)) < 0, г € [¿,¿1], I € [io.ii] (2.8)

также гарантирует достаточность ПМ(1).

Замечание 1. В [1] показано, что (2.7), (2.8) эквивалентны.

Пример 2.

1

Ф(и) = I [(ж^) - -7и2(г)

о

„о ™ /п\ ™0

М —> тах,

¿1 = Ж2, ¿2 = и, Ж1(0) = Ж1, жг(0) = ж2, [^(¿)| ^ 1, ^ е [0,1].

Сопряженная система:

01 = ~и, 02 = -01 - U, 0l(l) = 0, 02(1) = 0. Рассмотрим управление и = 1. Соответствующие траектории:

xi(t, 1) = h2 + x°2t + х°1} x2(t, 1) = t +

0!(í, 1) = 1 - t, 02(í, 1) = \{t - l)2 - (t - 1). Вспомогательная вектор-функция p(r,t):

Pi{r,t) = 1, p2(T,t) = t-T + 1. В результате первое условие принимает вид

x4 + x°2(t + l) - 3Í + | > 7, Í € [0,1], что эквивалентно двум неравенствам

я? +^2 + ^ 7, Ж? + 2х°2 - ^ 7. (2.9)

Второе неравенство после несложных вычислений конкретизируется в виде

+ 1) + 3Í- | > ¿€[0,1], что приводит к неравенствам

Х°х + х°2 - ^ ъ х°х + 2х°2 + ^ 7. (2.10)

Подведем итог. Если начальное состояние фазовой системы (xi,x2) и параметр 7 > 0 удовлетворяют условиями(2.9) или (2.10), то управление и = 1 является оптимальным.

Замечание 2. Управление и = — 1 обслуживается аналогичным образом: неравенство (• • • ^ —7) и переход к операции тах.

v£V

3. Комбинированное управление

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

192 в. а. срочко

Зададимся, например, следующей структурой управления

u(t) = < (3.1)

\l, t<=[e,t i],

где

в € (io,il), ü(t) = -((a(t),x(t,u)) + (b(t),iP(t,u))y

Принцип максимума для управления u(t) обеспечивается неравенствами

\{a(t),x(t,u)} + {b(t),iP(t,u)}\ <7, te[t0,e\, (a(t),x(t,u)) + (b(t),ip(t,u)) ^7, te [в,h}.

Первое достаточное условие оптимальности получается после замены x(t,u) —> x(t,v), v € V и перехода к экстремальным значениям для {a(t),x(t,v)}. В результате приходим к неравенству (2.3) для t € [to,0] и к неравенству (2.5) для t € [0,t 1]. Отметим, что для t € [ö,t\] ip(t,u) = ф(1,1).

Особого внимания требует точка в сопряжения двух участков управления, в которой должны выполняться неравенства (2.2), (2.5) при t=9

т,ф(в,и)) + (Р(г0,в),хо) +! тт),Р(т,втт

{Ь(в),ф(в,и)) + {р(1 о,в),х°)~ I\(Ъ{т),р{т,в))\(1т> 7.

¿0

Отсюда получаем два условия типа равенства в

I\{Ъ(т),р(т,в))\<1т = 0 (Ь(т),р(т,9))=0, те [Ь,9],

т,ф(в,и)) + {Р(10,в),х°) = 7.

С учетом утверждения леммы 1 первое условие приводит к соотношениям

(а(в),х(в,у)) = (а(в),х(в,и)) = (р(*0,в),х°), V € V.

Следовательно, второе условие есть условие стационарности в точке в

Ни(ф(в,и),х(в,и),и(в),в) =0,

которое выполняется в силу ПМ(и).

Таким образом, два неравенства (^ 7, ^ 7) эквивалентны одному условию типа равенства

(Ь(Т),Р(Т,9))=0, ТЕ[Ь,9]. (3.2)

Сформулируем итоговое утверждение.

Теорема 3. Пусть для управления u(t) вида (3.1) выполнены следующие условия:

1) неравенство (2.3) для t € [to, в];

2) неравенство (2.5) для t € [0,ti];

3) равенство (3.2) для точки в € (to,t\).

Тогда u(t) - оптимальное управление в задаче (1.1).

Замечание 3. Условие типа равенства (условие ортогональности) появляется в точке в сопряжения двух участков управления u(t) и является условием состыковки двух неравенств, связанных с этой точкой. Следует отметить, что аналогичный результат имеет место для задачи с билинейным функционалом в точке переключения управления [3].

Замечание 4. Другие варианты комбинированного управления u(t) в рамках набора { ±1, Ti,(t) } обрабатываются по аналогичной схеме.

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

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

1. Аксенюшкина Е. В. Достаточные условия оптимальности для одного класса невыпуклых задач управления / Е. В. Аксенюшкина, В. А. Срочко // Журн. вычисл. матем. и матем. физ. — 2015. — Т. 55. — № 10. — С. 1070-1080.

2. Антипина Н. В. Линейные функции Ляпунова - Кротова и достаточные условия оптимальности в форме принципа максимума / Н. В. Антипина, В. А. Дыхта // Известия вузов. Математика. — 2002. — №12. — С. 11-22.

3. Антоник В. Г. Условия оптимальности типа принципа максимума в билинейных задачах управления / В. Г. Антоник, В. А. Срочко // Журн. вычисл. матем. и матем. физ. — 2016. — Т. 56. — № 12. — С. 2054-2064.

4. Батурин В. А. Приближенные методы оптимального управления, основанные на принципе расширения / В. А. Батурин, Д. Е. Урбанович. — Новосибирск : Наука, 1997. — 175 с.

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

6. Дыхта В. А. Неравенство Ляпунова - Кротова и достаточные условия в оптимальном управлении / В. А. Дыхта // Итоги науки и техники. Современная математика и ее приложения. — 2006. — Т. 110. — С. 76-108.

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

8. Срочко В. А. Достаточные условия оптимальности экстремальных управлений на основе формул приращения функционала / В. А. Срочко, В. Г. Антоник // Известия вузов. Математика. — 2014. — №8. — С. 96-102.

Владимир Андреевич Срочко, доктор физико-математических наук, профессор, Институт математики,экономики и информатики, Иркутский государственный университет, 664003, Иркутск, ул. К. Маркса, 1, тел.: +7 (3952) 521276 (e-mail: [email protected])

V. A. Srochko

The Simplest Nonconvex Control Problem. The Maximum Principle and Sufficient Optimality Conditions

Abstract. The optimal control problem with linear phase system and linear-quadratic functional is considered. The transition from the maximum principle to sufficient optimality conditions is fulfilled with the help of the notion of strongly extremal control. It means that in the problem of maximization of Pontryagin's function phase or conjugate trajectory should be replaced with any admissible trajectory. Sufficient conditions give opportunity to obtain explicit expressions for extremal values of auxiliary problems contained in these conditions. Results are presented in the form of inequalities and equalities for a function with one variable with respect to time segment. A special situation is implemented in the analysis of the combined control with interior and boundary segments with respect to the constraint. At the point of connection of these segments there is a non-standard condition of maximum type.

A positive factor is dual nature of obtained results: it is a pair of symmetrical relations, which operate independently. Their origin is connected with two types of strongly extremal controls with respect to phase or conjugate variables.

Keywords: optimal control problem; the maximum principle; sufficient optimality conditions.

References

1. Aksenyushkina E. V., Srochko V. A. Sufficient optimality conditions for one class of nonconvex control problems . Zhurnal vichislitelnoy matematiki i matematicheskoy fiziki, 2015, Vol. 55, №10, pp. 1070-1080. (in Russian)

2. Antipina N. V., Dychta V. A. Linear funtions of Lyapunov - Krotov and sufficient optimality conditions in the form of maximum principle. Izvestiya vuzov. Matematika, 2002, №12, pp. 11—22. (in Russian)

3. Antonik V. G., Srochko V. A. Optimality conditions of the maximum principle type in bilinear control problems. Zhurnal vichislitelnoy matematiki i matematicheskoy fiziki, 2016, Vol. 56, №12, pp. 2054-1064. (in Russian)

4. Baturin V. A., Urbanovich D. E. Approximate optimal control methods based on the extension principle. Novosibirsk, Nauka, 1997, 175 p. (in Russian)

5. Gurman V. I. The extension principle in control problems. Moscow, Nauka, 1985, 288 p. (in Russian)

6. Dychta V. A. The inequality of Lyapunov - Krotov and sufficient conditions in optimal control theory. Itogi nauki i tehniki. Sovremennaya matematika i eyo prilozheniya, 2006, Vol. 110, pp. 76-108.(in Russian)

7. Krotov V. F., Gurman V. I. Methods and problems of optimal control. Moscow, Nauka, 1988, 446 p. (in Russian)

8. Srochko V. A., Antonik V. G. Sufficient optimality conditions of extremal controls based on functional increment formulas. Izvestiya vuzov. Matematika, 2014, №8, pp. 96—102. (in Russian)

Vladimir A. Srochko, Irkutsk State University, 1, K. Marks St., Irkutsk, 664003, professor, Phone: +7(3952)521276 (e-mail: [email protected])

i Надоели баннеры? Вы всегда можете отключить рекламу.