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

УДК 517.977

И.-Х.Д. Хишектуева

МЕТОД ПОИСКА НЕПОДВИЖНЫХ ТОЧЕК В ЗАДАЧЕ УЛУЧШЕНИЯ УПРАВЛЯЮЩИХ ПАРАМЕТРОВ

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

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

I.-Kh.D. Khishektueva

SEARCH METHOD OF FIXED POINT FOR IMPOVEMENT OF CONTROL PARAMETERS

The method of non-local improvement are considered to focus on a linear control parameters problem for optimization dynamical systems. Not-local improvement is provided with the solution of problem on fixed point for the continuous and unequivocal correspondence designed on the basis of projective operation on set of values of parameters. Article is devoted to algorithms of the solution of problem on fixed point for perturbation method at realization of not-local improvement control parameters.

Keywords: non-local improvement, projective methods, problem on fixed point.

1. Метод нелокального улучшения

Рассматривается линейная по управлению задача:

Ф(и) = p(x(t1)) + [ ((a(x(t), t), u) + d(x(t), t))dt ^ min , (1)

JT ' ' ueU

X(t) = A(x(t), t)u + b(x(t), t), x(t0) = x0, и eU , t e T, (2)

где матричная функция A(x, t), векторные функции a(x, t), b(x, t), функции p(x), d(x, t) являются

квадратичными по x и непрерывными по t на множестве Rn X T . Множество U с Rm выпукло и ком-

пактно.

Функция Понтрягина и стандартная сопряженная система имеют вид:

H(/, x,и,t) = H0 (/,x,t) + (H1(/,x,t),u), (3)

H0(/,x,t) = /,b(x,t))-d(x,t), H1(/,x,t) = AT(x,t)/-a(x,t)

/(t) = -Hx (/(t), x(t), w, t), t e T (4)

Введем модифицированную сопряженную систему:

p (t) = - Hx (p(t), x(t), w, t) - 2 Hxx (p(t), x(t), w, t) z (t) (5)

Для допустимых управлений u, v обозначим AvФ^) - приращение целевой функции; x(t, v), t eT -

решение системы (2) при u = v, x(t0, v) = x0; /(t, v), t e T - решение системы (4) при w = v , x = x(t, v),

/(t1,v) = -px(x(t1,v)) ; p(t,u,v), te T - решение системы (5) при w = u , x(t) = x(t,u), z = x(t,v)-x(t,u) с краевым условием

p(ti, u, v) = -Px (x(ti, u)) -1 Pxx (x(ti, u)) J (6)

где частные производные по x вычисляются при y = x(t1, v) - x(t1, u).

Поставим задачу об улучшении управления u0 e U : найти управление v eU с условием A v Ф^0) < 0.

Точная формула приращения целевой функции в рассматриваемой задаче принимает вид [1]:

A^(u) = -^^ H1(p(t,u,v),x(t,v),t)dt,v - uj, (7)

Для заданного a> 0 введем отображение Wa по формуле

Wa(u, v) = P,j (u + a\ H1(p(t, u, v), x(t, v), t)dt), u eU, v eU , (8)

J T

Необходимое условие оптимальности для управления и е и с помощью отображения (8) представляется в виде

и = Wа (и, и), — > 0. (9)

Формула (7) позволяет построить метод нелокального улучшения с помощью операции проектирования (8).

Первый проекционный метод нелокального улучшения: для заданных а > 0 и и0 е и сформируем отображение W—(у) = Wа(и0,у), уеи и найдем решение уравнения

V = W1а(v). (10)

Показано [1], что решение уа е и системы (10) обеспечивает улучшение. В силу известного свойства проекции имеем

^ Их( р(г, и0, у—), х(г, у—), г)йг, уа - и0^ >—| |у— - и 0|| .

Отсюда из формулы (7) следует уменьшение целевой функции с оценкой

1 и ц2

Ф(у—)-Ф(и0)<-------у—-и0 . (11)

а

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

неподвижной точки уа е и отображения W1а. В силу оценки (11) строгое улучшение гарантируется (в

том числе и для управления и0 е и , удовлетворяющего необходимому условию (9)), если уа Ф и0.

2. Численные эксперименты

В качестве иллюстрации работы метода рассмотрим следующий пример:

1 2

Ф(и) = -1 х2(г )йг

2 0

х(г) = и, х(0) = 1, и еи = [-1,1], г е Т = [0,2]

0

и0 = 0

Данный пример допускает аналитическое решение.

Функция Понтрягина и модифицированная сопряженная система принимают вид:

и 1 2

Н = ри — х

2

р(г) = х(г) +1 у(г), р(2) = 0.

Фазовая траектория, соответствующая управлению и0: х(г, и0) = 1, г еТ . Значение функционала для начального управления: Ф(и0) = 1.

Решение х(г,у) фазовой системы, соответствующее управлению у : х(г, у) = уг +1, г е Т.

Решение р(г, и0, у) модифицированной сопряженной системы:

р(г, и0, у) = у^ + г - (2 + у)

Сформируем отображение W1а(v) = Wа(и0,у):

г 4

W— (у) = Ри (и0 +—I Н1(р(г, и0, у), х(г, у), г)Ж) = Ри (—у - 2) =

2

-1, — у - 2 <-1

4 3

—у - 2, у е [-1; —), 34

4 2

—у - 2, 3

3

4

< 1 =

1, - 4 у - 2 > 1

3

3

-1, у е [-—;1]

4

Таким образом, задача о неподвижной точке имеет вид:

V - р (-—V - 2)

(12)

Уравнение допускает аналитическое решение

у = - 6 я -0,857142857.

7

при тф 0:

Численные расчеты применялись для исходного (12) и эквивалентного уравнения

4

у = у + т(у - Ри (- - у - 2)) (13)

Для решения обоих уравнений применялся метод простой итерации [2].

В ходе численных экспериментов установлено, что итерационный процесс для решения исходного уравнения (12) методом простой итерации сходится при значениях —<1. При стандартном значении а = 1 процесс не сходится по начальному приближению даже в минимально возможной машинной окрестности точного решения, при этом процесс итераций зацикливается с повторением значений -1 и -2/3. При уменьшении а до 0,7 наблюдается аналогичный тип расходимости. При дальнейшем уменьшении, начиная со значения а = 0,7 , итерационный процесс сходится, причем наилучшая сходимость получена при значениях в промежутке [0,01; 0,1]. Область сходимости по начальному приближению включает допустимое множество и =[-1;1]. При уменьшении параметра а вплоть до значений, сравнимых с машинным нулем, сходимость процесса итераций замедляется, и область сходимости по-прежнему включает все множество допустимых значений. Далее из-за ошибок округления он перестает сходиться и останавливается.

При увеличении а >1 процесс сходится, получаем точку у = -1.

Для модифицированного уравнения (13) при стандартном значении проекционного параметра а =1 варьировался параметр Т. При Т > 0 и его дальнейшем увеличении итерационный процесс расходится, наблюдается неограниченный рост последовательности приближений до +тс . При отрицательных значениях получаем сходимость на всей области допустимых значений управления, причем наилучшую по количеству итераций в диапазоне от -0,1 до -0,9. При уменьшении Т, начиная с -1, процесс итераций расходится, последовательность приближений неограниченно убывает до -тс .

Заключение

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

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

Литература

1. Булдаев А.С. Методы возмущений в задачах улучшения и оптимизации управляемых систем. - Улан-Удэ: Изд-во БГУ, 2008. - 260 с.

2. Самарский А.А., Гулин А.В. Численные методы: учеб. пособие для вузов. - М.: Наука. Гл. ред. физ-мат. лит., 1989. - 432 с.

3. Булдаев А.С., Хишектуева И.-Х.Д. Об одном методе улучшения управляющих параметров нелинейных систем // Инфокоммуникационные и вычислительные технологии и системы: матер. III Междунар. конф. (Россия, Бурятия, Улан-Удэ - оз. Байкал, 6-11 сентября 2010 г.). Улан-Удэ: Изд-во БГУ, 2010. - С. 79-82.

4. Хишектуева И.-Х.Д., Булдаев А.С. Оценка параметров динамических систем методом нелокальной оптимизации // Обобщенные постановки и решения задач управления: матер. V Междунар. симпозиума. Улан-Батор: Изд-во Монгольского государственного университета науки и технологии, 2010. - С. 187-189.

Хишектуева Ишин-Хорло Дамбадоржиевна, аспирант кафедры прикладной математики Бурятского государственного университета, тел. 89148363885, [email protected]

Khishektueva Ishin-Khorlo Dambadorzhievna, assistant at the chair of applied mathematics in the Buryat State University.