CC BY
72
СРАВНИТЕЛЬНЫЙ АНАЛИЗ РЕШЕНИЙ ЗАДАЧИ ОПТИМАЛЬНОГО УПРАВЛЕНИЯ КОНКУРЕНТНЫМ ПОВЕДЕНИЕМ ПРЕДПРИЯТИЙ СОТОВОЙ СВЯЗИ
И.П. Болодурина, Т.А. Огурцова, Е.П. Маас
В работе предложена динамическая модель конкурентного поведения предприятий сотовой связи в виде системы дифференциальных уравнений с запаздыванием. Реализованы численные методы и алгоритмы для решения задачи оптимального управления поведением операторов сотовой связи в условиях конкурентной борьбы за потребителей услуг связи. Сформулирована и решена задача оптимального управления поведением предприятий сотовой связи с постоянным запаздыванием с использованием принципа максимума Понтрягина для систем с постоянным запаздыванием. Получены условия оптимальности для решения задачи управления динамикой развития абонентской базы фирмы № 1. Реализован метод линеаризации системы нелинейных дифференциальных уравнений с запаздыванием и проведен поиск оптимального решения задачи посредством операций улучшения управления, повторяемых в итерационной процедуре, на примере метода Шатровского.
Ключевые слова: математическая модель, запаздывание, оптимальное управление, принцип максимума Понтрягина, улучшение управления.
Введение
Телекоммуникационная отрасль относится к важнейшим секторам, обеспечивающим функционирование и согласованную работу всех отраслей экономики. Для описания конкуренции и взаимодействия нескольких фирм на рынке телекоммуникационных услуг, а также общей динамики развития абонентской базы двух конкурирующих экономических агентов в работе предложено использование логистической модели Лотки— Вольтерра с запаздыванием во времени. Решению задач, касающихся вопросов управления в компаниях сотовой связи посвящено значительное количество работ. В частности, в работе [4] построена модель, описывающая поведение предприятий сотовой связи на рынке в условиях неопределенности и предложен алгоритм управления поведением предприятия сотовой связи, реализованный по принципу обратной связи. В работе [9] исследованы модели рынка, описывающие поведение потребителей в отношении нескольких участников рынка, решена задача оптимального управления поведением предприятий сотовой связи линейно-квадратичного критерия оптимальности. Однако не решены задачи определения оптимального управления развитием предприятий сотовой связи для динамических моделей, имеющих запаздывание в фазовых переменных, что по большей части и определяет научную новизну работы.
Математические модели большинства экономических объектов и процессов являются нелинейными. При решении задач анализа и синтеза систем хорошо обоснованны только методы классической линейной теории, поэтому применяют различные методы линеаризации, при которых нелинейные зависимости заменяются эквивалентными линейными моделями. При решении многих оптимизационных задач для описания объектов и систем применяют модели на основе линейного приближения ряда Тейлора. Одним из направлений использования дифференциальной линеаризации являются числен-
ные методы, которые используются при синтезе систем управления и обработки информации, а также в идентификации параметров модели.
В первом разделе статьи рассматривается задача оптимального управления поведением предприятий сотовой связи с использованием принципа максимума Понтрягина для систем с постоянным запаздыванием. Во втором разделе предлагается метод линеаризации системы нелинейных дифференциальных уравнений с запаздыванием. В заключении сформулированы основные выводы и результаты работы.
1. Численное решение задачи оптимального управления поведением предприятия сотовой связи с учетом запаздывания
В настоящее время в России действуют три крупнейших оператора сотовой связи: МТС, Билайн, Мегафон и множество других менее развитых компаний, также предоставляющих услуги сотовой связи. Для описания конкурентного взаимодействия операторов сотовой связи на рынке телекоммуникационных услуг, а также общей динамики развития абонентской базы двух конкурирующих экономических агентов рассмотрим логистическую модель (1) Лотки-Вольтерра с запаздыванием во времени [3],
1Х1 () _ Х1 ()[^1 — У11Х1 ( — Т) — У12Х2 ( — Т)] — р11и1 () — р12и2 ()
\Х2 (^) _ Х2 ()[^2 — У21Х1 ( — Т) — У22Х2 ( — Т)] — р21и1 () — р22и2 () где х; () - число абонентов / -го оператора сотовой связи в момент времени ?;
£ 1, г = 1,2 - коэффициент прироста абонентской базы I-го оператора сотовой связи;
у1к, г, k = 1,2 - коэффициент взаимного влияния I-го и к-го предприятий, предоставляющих услуги сотовой связи;
р , i, k = 1,2 - коэффициент влияния средней стоимости минуты на прирост числа абонентов ¿-го оператора сотовой связи;
т - временной лаг (запаздывание), связанный с разницей во времени между изменениями в рыночной ситуации и моментом принятия управленческих решений с целью реагирования на эти изменения;
иг ) - средняя стоимость минуты пользования услугами связи /-го оператора в момент времени ?, удовлетворяющий ограничению (2)
а < и() < в, I е [0,Т] (2)
Число абонентов /-го оператора сотовой связи на начальном интервале
[-т, 0] задано функциями (t), / = 1,2
Хх () = Фх (X * е [-Т,0] (3)
Фазовые ограничения описаны неравенствами (4).
МI ^ X (^ П > 1 =12 (4)
Задача оптимального управления поведением предприятий сотовой связи состоит в оптимизации ценовой политики и1(Ь) первого экономического агента, с целью максимизации числа его абонентов.
T
J(u1) = I bx1 (t)dt ^ max (5)
0
Поставленная задача относится к классу задач оптимального управления с постоянным запаздыванием [1]. Поэтому для ее решения применим принцип максимума Понтрягина для систем с постоянным запаздыванием в предположении, что значение стоимости минуты связи фирм конкурентов U2=const. Для этого построим функцию Понтрягина (6), где yl = xt (t - т) .
H(t, x, y, u ,w, Л0) =-Ao bxx + Wi (^ Xi - Г11 Xj У1 - Г12 Xi У2 - A1U - P12 U2) +
+ W2 (^2 X2 - /21X2 У1 - Y22 X2 У 2 - P21U1 - P22 U2 ) ^
Пусть (x (t), Uj (t)) - локально оптимальный процесс в поставленной задаче.
Тогда с необходимостью существуют множитель > 0 и абсолютно непрерывная вектор-функция щ(t) = (t),щ2 (t)) , t е[0, T] такие, что выполняются следующие
условия:
- условие максимума H (t, x, y, uj, у, ) = maxH(t, x, y, uj, y, )
или Uj (t)( PjjVj (t)+puV 2(t)) = min u(t)(pxxW j(t) + p ij¥ 2(t));
a <uY < в
- сопряженные функции уг (t), i = 1,2, удовлетворяют системе дифференциальных уравнений
дН дН
¥i(t) = - — (t + т) - — (t) = ¥l(t + т )ynXl{t + т) + ¥2 (t + т )y2ix2 (t + т) - Л>ь - ¥1 (t)fo - Yiiyi(0 - УиУг (t')) дУ dxi
V2 (t) =-dH (t + T) -dH (t) = ¥1(t + T)Y12X1(t + T) + ¥2(t + T)Y22X2 (t + T) - ¥2 (t)(^2 - Г21У1(t) - Г22У2 (t))
дУ 2 dx2
на отрезке t0 < t < T-T, и ^(t) = y2 (t) = 0 на отрезке T-т < t < T;
- условия трансверсальности (T) = 0, (T) = 0;
- условия допустимости (1)-(4).
Введем функцию переключения: - р11щ 1 - р21щ2 = р(t). Для нахождения оптимального управления необходимо решить задачу
Щ (t)[-P11W1 - P21W2 ] ^ max, ГДе U = {u е R :а< щ < в, а > 0, в > 0}.
ueU
Из условия максимума функции Понтрягина находим оптимальное управление (7)
Р, р(() >
а, р(t) < ° (7)
е (а, в), р(t) = 0.
ui ) =
Наличие фазовых ограничений (4) учтем в функционале (5) путем введения штрафных слагаемых.
Для построения численной схемы приближённого решения поставленной задачи оптимального управления динамикой развития абонентской базы, проведем дискретную аппроксимацию непрерывной задачи с точностью первого порядка. Дискретная задача оптимального управления, аппроксимирующая непрерывную задачу (1)-(5), состоит в нахождении минимума функции (9) при ограничениях
(10)-(12).
q-1 _ д
1 (и1)=- X Ьх1 ^+А X
г=0 г=0
(х' Ц
Еmaxz {hJ | X |,0}
J=1
'{gj' (x')'0}
Е maX {g J |x ),0\
j =1
i+1
x.
i+1
i i A i i — V * i i — V i i i i А
= x1 + e1 x1 Дt - y11 x1 x1 Дt - y12x1 x2 Дt - p11u1 Дt - p12u2 Д
i i i i —V i i —V i i
= x2 + £2 x2 Д - 721x2 xi Дt - Y22 x2 x2 Д - p21U 2 Дt - P22 U 2 Д
i = 0, q -1
i (0)
= ^(iAt), i = - v,0, j = 1,2,
(9)
(10)
(11)
(12)
а < щ < в, а > 0, в > 0, I = 0, q -1.
где к,х) = п- х,- {() < £,• ((, х) = х,. - ^(() < 0 .
Параметры задачи подобраны с учетом анализа статистических данных абонентской базы, динамики предыдущей тарифной политики и тенденции развития рынка.
x
2. Поиск оптимального управления поведением
предприятий сотовой связи, посредством операций улучшения управления
Один из приближенных методов решения задач оптимального управления нелинейными объектами предложен Л.И. Шатровским [5]. Он относится к методам линеаризации, использующим первую вариацию, и основан на итерационной процедуре, в ходе которой при заданном в виде функции времени начальном приближении управления, на каждом шаге решается линеаризованная задача с использованием рекуррентных вычислений. В ходе итерационной процедуры произвели приближение функции управления щ (t) к оптимальному управлению щ (t). Достаточное количество итераций определили неравенством (13), где щ (t)
- функция управления, найденная на к-ом шаге, х - точность минимизации
функционала.
J(u” -1) - J(u к ) < X (13)
Начальное приближение управления выбрали равным и0 (t).
После аппроксимации дифференциальных уравнений отрезком ряда Тейлора с учетом вторых производных, задача улучшения управления щ (t) сводится к задаче поиска такой вариации vk (t), которая минимизировала бы приращение функционала (14)
J(u) = Уз(Т) ^min (14)
на решениях линеаризованной системы (15)
¿+1 ~ ~ ~ / \ А V
У1 = У1 + (а11 У1 + а12У2 + а13Уз - Рп^ )А?
* у2 = У2 + (а21 у1 + а22У2 + а23У3 — ^21^ )Д^ (15)
у3 = У3 + (а31 У1 + а32У2 + а33 У3 )Д^
С целью определения оптимальных динамических траекторий и оптимального управления в среде Borland Delphi 7.0 разработан программный комплекс, реализующий численное решение поставленной задачи оптимального управления конкурентным поведением предприятий сотовой связи, в основе которого лежит метод проекции градиента и метод улучшения управления (метод Шатровского).
По имеющимся статистическим данным абонентской базы, динамики предыдущей тарифной политики и тенденции развития рынка созданный программный комплекс позволяет находить оптимальные управленческие сценарии для различных предприятий, предоставляющих услуги связи. Значения функционала качества показаны в таблице.
Таблица
Значения функционала качества 3(щ), характеризующего прирост числа абонентов фирмы № 1 за период в 10 кварталов (млн. чел.)
Метод проекции градиента Метод Шатровского
Неуправляемая модель без запаздывания 91.23 91.03
Неуправляемая модель с запаздыванием 99.35 99.12
Управляемая модель без запаздывания 115.09 98.17
Управляемая модель с запаздыванием 127.002 108.79
Заключение
В работе решена задача оптимального управления поведением предприятий сотовой связи с использованием принципа максимума Понтрягина для систем с постоянным запаздыванием, а также методом последовательной линеаризации для поиска оптимального управления, посредством операций улучшения управления (методом Шатровского). В результате проведенного анализа решений установлено, что при использовании метода Шатровского не удалось улучшить функционал качества поставленной задачи по сравнению с методом проекции градиентов. Прирост числа абонентов фирмы № 1 при численном решении на основе принципа максимума Понтрягина на 16,71 % больше, чем при управлении, найденным по методу Шатровского. Таким образом, методы численного решения, базирующиеся на принципе максимума Понтрягина, являются более эффек-
тивным инструментом для наращивания абонентской базы предприятия в условиях конкурентной борьбы на рынке сотовой связи.
Работа выполнена при финансовой поддержке гранта РФФИ, проект 12-01-31325.
Литература
1. Андреева, Е.А. Управление системами с последействием / Е.А. Андреева, В.Б. Колмановский, Л.Е. Шайхет - М.: Наука. Гл. ред. физ.-мат. лит., 1992. - 336 с.
2. Болодурина, И.П. Дифференциальные уравнения с запаздывающим аргументом и их приложения: учебное пособие / И.П. Болодурина - Оренбург: Оренбург. гос. ун-т, 2006. - 101 с.
3. Болодурина, И.П. Управление ценой на предоставляемые услуги предприятий телекоммуникационной отрасли / И.П. Болодурина, Т.А. Огурцова // Проблемы управления. - 2011. - № 3. - С. 30-35.
4. Коблов, А.И. Модели и алгоритмы управления поведением предприятий сотовой связи в условиях неопределенности: дис. ... канд. техн. наук: 05.13.01: защищена 15.11.2007 / А.И. Коблов - ЮУрГУ, 2007. - 136 с.
5. Когут, А.Т. Модификация метода Шатровского решения нелинейных задач оптимального управления. / А.Т. Когут, А.А. Лаврухин // Омский научный Вестник. -Омск, 2005. - № 3. - С. 81-85.
6. Когут, А.Т. Применение алгоритмов линеаризации для идентификации и адаптивного управления в нелинейных динамических системах: монография / А.Т. Когут,
Н.А. Тихонова. Омский гос. ун-т путей сообщения. - Омск, 2008. - 126 с.
7. Прасолов, А.В. Динамические модели с запаздыванием и их приложения в экономике и инженерии / А.В. Прасолов - СПб.: Лань, 2010. - 192 с.
8. Шатровский, Л.И. Об одном численном методе решения задач оптимального управления / Л.И. Шатровский // Ж. вычисл. матем. и матем. физ., 2:3 (1962). -С. 488-491.
9. Модели прогнозирования развития региональных рынков и оптимальное управление поведением фирмы / В.И. Ширяев, И.Г. Гришин, А.И. Коблов // Вестник Пермского университета. Математика, механика, информатика. - 2007. - Выпуск 7(12). - С. 149-164.
Болодурина Ирина Павловна, д.т.н., профессор, зав. кафедрой прикладной математики, Оренбургский государственный университет (Оренбург, Российская Федерация), [email protected].
Огурцова Татьяна Александровна, ст. преподаватель кафедры прикладной математики, Оренбургский государственный университет (Оренбург, Российская Федерация), [email protected].
Маас Елена Петровна, студент кафедры прикладной математики, Оренбургский государственный университет (Оренбург, Российская Федерация), [email protected].
COMPARATIVE ANALYSIS OF NUMERICAL METHODS FOR SOLVING PROBLEMS OF OPTIMAL CONTROL THE BEHAVIOR OF ECONOMIC AGENTS IN THE CONDITIONS OF COMPETITIVENESS
Bolodurina I.P., Orenburg State University (Orenburg, Russian Federation) Ogurtsova T.A., Orenburg State University (Orenburg, Russian Federation)
Maas E.P., Orenburg State University (Orenburg, Russian Federation)
In work the dynamic model of competitive behavior of the enterprises of cellular communication in the form of system of the differential equations with delay is offered. Numerical methods and algorithms for the solution of a problem of optimum control by behavior of mobile operators in the conditions of competitive fight for consumers of communication services are realized. The problem of optimum control by behavior of the enterprises of cellular communication with continuous delay with use of the principle of a maximum of Pontryagin for systems with continuous delay is formulated and solved. Conditions of an optimality for the solution of a problem of management are received by dynamics of development of subscriber base of firm No. 1. The method of linearization of system of the nonlinear differential equations with delay is realized and search of the optimum solution of a task by means of operations of improvement of the management, repeated in iterative procedure, on the example of modification of a method Shatrovsky is carried out.
Keywords: mathematical model, delay, optimum control, management improvement.
References
1. Andreeva E.A. Upravleniye sistemami s posledeystviem [the management of systems with aftereffect]. Moscow, Nauka. The main editorial office of physico-mathematical literature, 1992. 336 p.
2. Bolodurina I.P. Differencialnie uravneniya s zapazdivausim argumentum i ih prilozheniya [Differential equations with retarded argument and their applications]. Orenburg. The state. un-t, 2006. 101 p.
3. Bolodurina I.P., Ogurtsova T.A. Upravlenie cenoy na predostavlyaemie uslugi predpri-yatiy telecommunicacionnoy otrasli [The management of the price for the provided services of the enterprises of the telecommunications industry]. Problemi Upravleniya [Problems of management]. 2011. No. 3. P. 30-35.
4. Koblov A.I. Modeli i algoritmi upravleniya povedeniem predpriyatiy sotovoy svyazi v usloviyah neopredelennosti [Models and algorithms for behavior management of the enterprises of cellular communication in the conditions of uncertainty]: the Dissertation of candidate of technical Sciences. SUSU, 2007. 136 p.
5. Kogut A.T. Modificaciya metoda Shatrovskogo resheniya nelineynih zadach optimal-nogo upravleniya [Modified method of Shatrovsky solution of nonlinear problems of optimal control]. Omskiy nauchniy Vestnik [Omsk scientific bullitin]. 2005. No. 3. P. 8185.
6. Kogut A.T. Primeneniye algoritmov lianerizacii dlya identifikacii i adaptivnogo uprav-leniya v nelineynih dinamicheskih sistemah [Application of algorithms of linearization for identification and adaptive control for nonlinear dynamic systems]. Omsk, The Omsk state. un-t of communication, 2008. 126 p.
7. Prasolov A.V. Dinamicheskiye modeli s zapazdivaniem i ih prilozeniya v economike i ingenerii [Dynamic models with delay and their applications in Economics and engineering]. St. Petersburg, Fallow deer, 2010. 192 p.
8. Shatrovsky L.I. Ob odnom chislennom metode resheniya zadach optimalnogo upravleniya [About one numerical method of the solution of problems of optimum control]. Zurnal vichislitelnoy matematiki i matematicheskoy fiziki [computational mathematics and mathematical physics]. 2:3 (1962), P. 488-491.
9. Shiryaev V.I., Grishin I.G., Koblov A.I. Jars B.M. Modeli prognozirovaniya razvitiya regionalnih rinkov I optimalnoe upravlenie povedeniem firmi [Models of forecasting the development of regional markets and optimal control of the behavior of the firm] Bulletin of the University of Perm. Mathematics, mechanics, information. 2007. Issue 7 (12). P. 149-164.
Поступила в редакцию 18 марта 2013 г.
