443.00Kb

advertisement
Параметрическое регулирование развития (эволюции) нелинейных динамических систем
Ашимов А.А., Сагадиев К.А., Боровский Ю.В., Искаков Н.А., Ашимов Ас.А.
В работе приводятся результаты исследования достаточных условий существования экстремали и точки
бифуркации экстремали задачи вариационного исчисления по выбору оптимальных законов параметрического
регулирования в среде заданного конечного набора алгоритмов при многопараметрическом возмущении.
Полученные результаты иллюстрируются на примере одной математической модели экономической системы
страны.
Ключевые слова: математическая модель, экономическая система, параметрическое
регулирование, экстремаль, функционал, ограничение, задача вариационного исчисления.
Как известно, поведение динамического объекта зависит как от значений внешних воздействий, так и от
значений внутренних параметров. Поэтому решение вопросов регулирования развития нелинейных
динамических систем, в том числе экономических систем стран, в отдельных случаях, можно свести к поиску
параметрических воздействий в виде экстремалей соответствующих задач вариационного исчисления. Так, в
последние годы в ряде работ [3-6] параметрические воздействия - экстремали соответствующих задач
вариационного исчисления (по выбору оптимальных законов параметрического регулирования в среде
заданного конечного набора алгоритмов) стали предлагаться для эффективного регулирования механизмов
рыночной экономики. В указанных задачах вариационного исчисления функционалы выражают некоторые
(глобальные, промежуточные или тактические) цели экономического развития, а фазовые ограничения и
ограничения в разрешенной форме представлены математическими моделями экономической системы из [7].
Математические модели из [7] содержат также ряд коэффициентов, изменение каждого из которых в
определенном интервале приводит к деформации (возмущению) рассмотренных задач вариационного
исчисления.
В известной литературе результаты исследования условий существования решений таких задач
вариационного исчисления по выбору оптимальных законов параметрического регулирования в среде заданного
конечного набора алгоритмов отсутствуют. В настоящее время широко исследуется параметрическое
возмущение задач вариационного исчисления. Так, в [1] параметрическое возмущение используется для
получения достаточных условий экстремума путем построения соответствующих S-функций и использования
принципа снятия ограничений. В [6] ставится вопрос об условиях устойчивости решений задач вариационного
исчисления (проблема Улама). Исследования этой проблемы сводится к нахождению условий регулярности, при
которых у функционала возмущенной задачи есть точка минимума близкая к точке минимума функционала
невозмущенной задачи. В [2] доказана теорема об условиях существования точки бифуркации для задачи
вариационного исчисления, функционал которой рассматривается на
пространстве Соболева
0
W pm () (2  p  ) и зависит от скалярного параметра  [0,1] . Таким образом, предлагаемый в работах [3-6]
подход параметрического регулирования развития нелинейных динамических систем требует исследования:
1) вопроса существования решения такой задачи вариационного исчисления в среде заданного
конечного набора алгоритмов;
2) влияния параметрического возмущения на результаты решения рассматриваемой задачи. Это в свою
очередь требует решения вопроса о существовании точек бифуркации экстремали задачи вариационного
исчисления в среде заданного конечного набора алгоритмов.
Ниже приводятся неохваченные известными работами результаты исследования достаточных условий
существования экстремали и точки бифуркации экстремали задачи вариационного исчисления по выбору
оптимальных законов параметрического регулирования в среде заданного конечного набора алгоритмов при
многопараметрическом возмущении.
1. . Условия существования решения задачи вариационного исчисления по выбору
оптимального набора законов параметрического регулирования на множестве сочетаний из p
параметров по r в среде заданного конечного набора алгоритмов
Пусть
R n - n-мерное евклидово пространство;
[t0 , t0  T ]
промежуток (времени); t [t 0 , t 0  T ] ; x  ( x , x ,..., x )  X  R
1
2
n
n
- фиксированный
- вектор состояния
1
2
l
l
системы; u  (u , u ,..., u )  W  R - вектор (регулирующих) параметрических воздействий;
1
W, Х – компактные множества c непустыми внутренностями соответственно;
и
Int (W )
Int ( X )
  (1 , 2 ,, m )    R m , 
- открытое связное множество;
f
f
f
отображения f ( x, u,  ) : X  W    R n
и
,
,
непрерывны в X W   .
u

x
Исследуемая задача вариационного исчисления рассматривается на базе следующей системы
обыкновенных дифференциальных уравнений, зависящей от параметров (u ,  )
dx
 f ( x, u ,  ) , x(t 0 )  x0
(1)
dt
для фиксированных начальных условий x0  Int ( X ) .
Пусть x (t ) - решение указанной задачи (1) на промежутке [t0 , t0  T ] при
постоянных значениях u W и
  .
Пусть x (t )  Int ( X ) . Решение задачи (1) для
1
2
l
выбранного u*  (u* , u* ,..., u* )  W обозначим через x* (t ) . Далее u* фиксировано.
Обозначим через  замкнутое множество в пространстве непрерывных вектор n l
функций C [t 0 , t 0  T ] , состоящее из всех непрерывных вектор-функций
удовлетворяющих следующим ограничениям.
x(t ), u(t )
x  X , u W , x j (t )  x*j (t )  x*j (t ) ,
(2)
t  [t0 , t0  T ] , j  1  n,   0 .
Пусть {F i ( x) : i  1  p} и G ( x)  0 - конечный набор непрерывных для x  X
F i
также непрерывны в X . Возможность
x j
выбора оптимального набора законов параметрического регулирования на множестве
сочетаний из p параметров по r и на промежутке времени [t0 , t0  T ] исследуется в среде
следующих алгоритмов (законов управления):
вещественнозначных функций. Все функции
U
ij

 kij F i ( x)  u*j , i  1  p, j  1  l .
(3)
Здесь, k ij  0 - настраиваемые коэффициенты. Использование набора r ( 1  r  l , здесь и
далее фиксировано) законов
U ij из (3) при фиксированных k ij в системе (1) означает
подстановку в правые части уравнений системы функций
{u js  Uis js } для r различных
значений индексов js (1  s  r , 1  js  l , 1  is  p ). При этом остальные
u j , где j не входит
j
в указанное множество значений j s , считаются постоянными и равными значениям u * .
Для решений системы (1) при использовании r законов управления вида
рассматривается следующий функционал (критерий):
K
{u js  Uis js }
t0 T
 G( x (t ))dt .
(4)
t0
2
Постановка задачи выбора набора законов параметрического регулирования на
множестве сочетаний из p параметров по r в среде заданного конечного набора алгоритмов
имеет следующий вид.
При фиксированном    найти набор из r законов
U  {U i
s js
, s  1  r}
из набора алгоритмов (3), который обеспечивает верхнюю грань значений критерия (4) -
K  sup
(5)
U
при выполнении условий (1,2) для заданного периода времени.
Следующая теорема доказывается в работе.
Теорема 1. При использовании любого выбранного набора законов
U  {U i
s js
, s  1  r } , где r  l из набора алгоритмов (3) при ограничениях (1) и
(2) существует решение задачи нахождения верхней грани критерия K:
t0 T
 G( x (t ))dt 
t0
При
этом
, , ki
если
sup
( ki1 j1 , ki2 j2 , , kir jr )
множество
.
возможных
(6)
значений
коэффициентов
законов рассматриваемой задачи ограничено, то указанная
верхняя грань для выбранного набора из r законов достигается. Для конечного
набора алгоритмов (3) задача (1)-(5) имеет решение.
(ki
1 j1
, ki
2 j2
r jr
)
2. Достаточные условия для существования точки бифуркации экстремалей задачи
вариационного исчисления по выбору оптимальных законов параметрического
регулирования в среде заданного конечного набора алгоритмов на базе математической
модели страны
Введем следующее определение, характеризующее такие значения параметра
котором возможна замена одного оптимального закона регулирования на другой.
Определение. Значение
 , при
*   называется точкой бифуркации экстремали задачи
   существуют как минимум два различных оптимальных набора из r

законов из (3), отличающихся хотя бы на один закон U ij , а в каждой окрестности точки 
(1)-(5), если при
найдется такое значение
   , для которого задача (1)-(5) имеет единственное решение.
Следующая теорема дает достаточные условия для существования точки бифуркации
экстремалей рассматриваемой вариационной задачи по выбору закона параметрического
регулирования в заданной конечной среде алгоритмов.
Теорема 1 (о существовании точки бифуркации). Пусть при значениях параметра
1 и 2 , ( 1  2 ,
1 , 2   ) задача (1)-(5) имеет соответствующие единственные
3
решения для двух различных оптимальных наборов из r законов из (3), отличающихся хотя
бы на один закон U ij . Тогда имеется хотя бы одна точка бифуркации
   .
Следующая теорема является непосредственным следствием теоремы 1.
Теорема 2. Пусть при значении
  1
регулирование с помощью некоторого
набора r законов из (3) дает решение задачи (1)-(5), а при   2 , ( 1  2 , 1 , 2   )
регулирование с помощью этого набора законов не дает решение задачи (1)-(5). Тогда
имеется хотя бы одна точка бифуркации
   .
Численный алгоритм нахождения бифуркационного значения параметра  .
Соединим точки 1 и 2 гладкой кривой T   . Разобьем эту кривую на n равновеликих
частей с достаточно малым шагом. Для полученных значений
 k  T , k  0  n,  0  1 ,  n  2
регулирования закона
U0
Uk
определяются оптимальные наборы r законов
и находится первое значение k, при котором этот закон отличается от
хотя бы по одному значению индекса. В этом случае точка бифуркации параметра
 * лежит на дуге
( k 1 ,  k ) .
Для найденного участка кривой алгоритм определения точки бифуркации с заданной
точностью  состоит в применении метода половинного деления. В результате находится
c  ( k 1 ,  k ) , с одной стороны от которой в пределах отклонения  от значения с
на кривой T оптимальным законом является U 0 , а с другой - в пределах отклонения  от
точка
значения с этот закон оптимальным не является. Из доказательства теоремы 1 следует, что
существует точка бифуркации  на указанной дуге. Эти алгоритмы были использованы при
нахождении множества точек бифуркации экстремали в следующем примере.

2. Пример нахождения точек бифуркации экстремалей для одной задачи
вариационного исчисления на базе математической модели экономической системы
Рассмотрим возможность нахождения точек бифуркации для одной вариационной
задачи [3] по выбору набора оптимальных законов параметрического регулирования
механизма рыночной экономики в среде конечного фиксированного набора алгоритмов на
базе математической модели [5] экономической системы страны. Фазовые ограничения и
ограничения в разрешенной форме исследуемой задачи вариационного исчисления по выбору
законов параметрического регулирования представлены следующими соотношениями.
dM  I

 M ,
dt
pb
dQ

 Mf  ,
dt
p
dLG
 rG LG   G  n p   nL sR L  nO (d P  d B ) ,
dt
(7)
(8)
(9)
4
dp
Q
 
p,
dt
M
(10)
 Rd  RS  L
ds s
d
S
 max 0,
, R  min{ R , R } ,
S
dt 
R


(11)
Lp 
(12)
1
LG ,

1
dp 
r2 LG ,

d B  r2 LG ,
(13)
(14)

   s 
x
1   
1     p 
1 


R  Mx ,
 0   0 pMf ,

,



(15)
d
(16)
(17)
 G  pMf ,
(18)
 1  
f  1  1 
x



 L  (1  n L ) sR d ,
I 
1
1
,
(19)
(20)
1

  (1   )n p
(1  n )
p
G
P
0
     ,
1
R S  P0A exp(  p t )
,
1  
0

 
(21)
 n0 (d  d )  n p   nL  (1  nL )n p sR  (   rG ) L ,
B
G
L
L
*
p
I
(22)
L
.

pP0 exp(  p t )
(23)
Начальные условия для этой системы для краткости не указаны.
Здесь обозначены: М – суммарная производственная мощность; Q – общий запас
товаров на рынке; LG – общий объем государственного долга; p– уровень цен; s – ставка
заработной платы; Lp – объем задолженности производства; dp и dB–соответственно
предпринимательские и банковские дивиденды; Rd и RS – соответственно спрос и
предложение рабочей силы; δ, v - параметры функции f(x), x – решение уравнения
f ( x)  s ; ФL и Ф0 – соответственно потребительские расходы трудящихся и
p
собственников; ФI – поток инвестиций; ФG – потребительские расходы государства; ξ норма резервирования; β – отношение средней нормы прибыли от коммерческой
деятельности к норме прибыли рантье; r2 – ставка процента по депозитам; rG – ставка
процента по облигациям государственных займов; η0 – коэффициент склонности
собственников к потреблению; π – доля потребительских расходов государства от
внутреннего валового продукта; np, n0, nL – соответственно ставки налогов на поток
платежей, дивиденды и доход трудящихся; b – норма фондоёмкости единицы мощности; μ –
коэффициент выбытия единицы мощности вследствие деградации; μ* - норма амортизации; α
– постоянная времени; Δ – постоянная времени, задающая характерный временной масштаб
процесса релаксации заработной платы; P0, P0A – соответственно начальные значения
численности трудящихся и общей численности трудоспособных; λp>0 – заданный темп
демографического роста; ω – душевое потребление в группе трудящихся.
5
Параметры модели и начальные условия для дифференциальных уравнений (7)-(11)
были получены на основе данных экономики Республики Казахстан за 1996-2000 годы [8]
(r2=0.12; rG=0.12; β=2; np=0.08; nL=0.12; s=0.1; n0=0.5; μ=μ*=0.012; Δ=1) или
оценены решением задачи параметрической идентификации (ξ=0.1136; π =0.1348; δ=0.3;
ν=34; η0=0.05; b=3.08; α=0.008; Q(0)=125000).
После последовательного исключения переменных стоящих в левых частях
соотношений (12)-(23), система (7)-(23) преобразуется в систему из пяти обыкновенных
дифференциальных уравнений, которая в силу ее громоздкости здесь не приводится.
Возможность выбора оптимального набора законов вида (3) параметрического
регулирования исследовалась:
-

на уровне одного из двух параметров  (j=1) ,
(j=2);
на промежутке времени [t0 , t0  T ] и в среде следующих алгоритмов.
1)U 1 j (t )  k1 j
M  M0
 const j ,
M0
2)U 2 j (t )  k 2 j
M  M0
 const j ,
M0
3)U 3 j (t )  k 3 j
p  p0
 const j ,
p0
4)U 4 j (t )  k 4 j
p  p0
 const j ,
p0
 M  M 0 p  p0 
  const j ,
5)U 5 j (t )  k 5 j 

p 0 
 M0
(24)
 M  M 0 p  p0 
  const j .
6)U 6 j (t )  k 6 j 

p 0 
 M0
Здесь M 0 , p 0 - начальные значения соответствующих переменных, а const j - оценка
соответствующего параметра, полученная по результатам параметрической идентификации
модели.
В рассматриваемой задаче использовался критерий (4) вида (среднее значение ВВП
за 1997-99 годы)
K
где
Y  Mf
1
T
t0 T
 Y (t )dt
t0
,
(25)
.
Замкнутое множество   C [t 0 , t 0  T ] в пространстве непрерывных вектор функций выходных переменных системы (1)-(23) и регулирующих параметрических
воздействий определяется следующими соотношениями
7
pij (t )  p ** (t )  0.09 p ** (t ),
( M (t ), Q (t ), LG (t ), p (t ), s (t ))  X ,
0  u j  a j , i  1,4, j  1,2, t  [t 0 , t 0  T ].
Здесь
aj
- наибольшее возможное значение j-го параметра,
pij (t )
(26)
– значения уровня цен
при U ij -ом законе регулирования; p (t ) - модельные (расчетные) значения уровня цен без
**
6
параметрического регулирования;
указанных параметров.
X
- компактное множество допустимых значений
В данной задаче вариационного исчисления рассматривалась ее зависимость от
двумерного коэффициента   (r2 , nO ) математической модели, возможные значения
которого принадлежат некоторой области (прямоугольнику)  на плоскости.
1
Случай неограниченных множеств H () мы не рассматриваем исходя из
экономических соображений. В этом случае, согласно лемме 1 верхняя грань критерия (25)
задачи (7)-(26) всегда достигается для некоторых значений коэффициентов набора из двух
законов.
В результате вычислительного эксперимента были получены графики зависимостей
оптимального значения критерия K от значений параметров (r2 , nO ) для каждого из 12
возможных законов U ij , i  1,6, j  1,2 . На рисунке 1 представлены указанные графики для
двух законов U 21 и U 41 , дающих наибольшее значение критерия в области  , линия
пересечения соответствующих поверхностей и проекция этой линии пересечения на
плоскость значений  , состоящая из точек бифуркации этого двумерного параметра. Эта
проекция делит прямоугольник  на две части, в одной из которых оптимальным является
закон управления U 21 , а в другой - U 41 , на самой проекции линии оба указанных закона
являются оптимальными.
7
Fig. 1.
3. CONCLUSIONS
Проведены исследования одной задачи параметрического регулирования развития нелинейных
динамических систем
Проведено исследование бифуркаций экстремалей одного класса задач вариационного
исчисления по выбору оптимального набора законов параметрического регулирования
динамических систем при многопараметрическом возмущении в среде конечного набора
алгоритмов содержащих по одному настраиваемому коэффициенту.
Доказан факт существования решения рассматриваемого класса задач вариационного
исчисления.
Предложено определение бифуркационной точки экстремалей; доказано утверждение
о достаточных условиях существования бифуркационной точки экстремалей и предложен
численный алгоритм обнаружения бифуркационного значения параметра для экстремалей
рассматриваемого класса задач вариационного исчисления.
Положения теоретических исследований проиллюстрированы на примере одной
задачи выбора оптимального набора законов параметрического регулирования механизмов
рыночной экономики в среде заданного конечного набора алгоритмов.
Полученные результаты могут быть рекомендованы для использования при разработке
и осуществлении эффективной государственной политики.
THE APPENDIX
Доказательство теоремы1 . Сопоставление набору значений коэффициентов
(ki
1 j1
, ki
2 j2
, , ki
r jr
) из набора законов U  {( U i j , k i j ), (U i j , k i j ),  , (U i j , k i j )}
11
11
2 2
2 2
r r
r r
соответствующих выходных функций и регулирующих параметрических воздействий
x (t ), u(t ) системы (1) при ее регулировании с помощью этого набора законов задает
непрерывное отображение H некоторого подмножества Rl  [0,) l в пространство
C n l [t0 , t0  T ] .
Полный прообраз H () множества  при отображении H замкнут согласно
теореме о замкнутости полного прообраза замкнутого множества при непрерывном
1
1
отображении. Множество H ij () не пусто, поскольку оно содержит начало координат Rl .
(При нулевых значениях коэффициентов функции x(t )  x* (t ), u (t )  u*  очевидно
удовлетворяет ограничениям (2)).

1
Сопоставление набору коэффициентов k  H () законов критерия (4) K для
решения системы (1) определяет непрерывную функцию
K : H 1 ()  [0, ) .
8
Следовательно, при выбранном наборе законов
задаче определения на замкнутом множестве
ограниченной функции
U
задача (1)-(5) равносильна
1
H () верхней грани непрерывной

y  K (k )
Эта функция является непрерывной в силу теоремы о непрерывной зависимости решения
системы обыкновенных дифференциальных уравнений от параметров, ограниченности этого
решения в силу включения x  X из (2) и непрерывной зависимости определенного
интеграла от параметра. Поэтому задача (1)-(5) для фиксированного набора законов U
всегда имеет решение, включающее конечное оптимальные значения критерия K * . Для
1
ограниченного множества H () это значение критерия достигается при некоторых

значениях коэффициента k (теорема о достижении наибольшего значения непрерывной
1
функции на компакте). Для неограниченного множества H () , может найтись

1
последовательность значений коэффициентов k из H () , соответствующие значения
критерия K для элементов которой стремятся к K * . Таким образом, доказан факт
существования решения задачи вариационного исчисления для случая одного закона
параметрического регулирования. Из конечности набора возможных законов регулирования
(3) следует справедливость леммы – факт существования решения задачи (1)-(5).
Доказательство теоремы 2. Соединим точки 1 и  2 гладкой кривой S лежащей в
области  : S  { (s), s [0, 1]},  (0)  1 ,  (1)  2 . Обозначим оптимальное значение
критерия K задачи (1)-(5) для выбранного набора законов регулирования
U  {U i
1 j1
,Ui
2 j2
, , U i
r jr
} и значения
 (s)
через K U (s ) . Функция
y  KU (s) является
непрерывной на [0,1] согласно теореме о непрерывной зависимости решения системы
обыкновенных дифференциальных уравнений, непрерывной зависимости определенного
*
интеграла и, в целом, в силу выше доказанной леммы. Функция y  max KU (s)  K (s) ,
U
дающая
решение
рассматриваемой
непрерывной на отрезке
задачи
(1)-(5),
следовательно,
также
является
[0,1] . Обозначим через (U )  [0,1] множество всех тех значений
параметра s, для которых KU (s)  K * (s) . Это множество замкнуто, как полный прообраз
{0} для непрерывной функции y  KU (s)  K * (s) . Множество (U )
может быть и пустым. В результате промежуток [0,1] представляется в виде следующего
замкнутого множества
конечного объединения, состоящего, как минимум, из двух замкнутых множеств (см. условия
теоремы)
[0,1]   (U ) .
U
9
0  (U ) для некоторого набора

законов U , соответствующего 1 , и 1 (U ) , то имеется граничная точка s множества
(U ) , находящаяся в промежутке (0,1) (будем считать, что s  - нижняя грань таких

граничных точек для множества (U ) ). Точка s также является граничной точкой

некоторого другого множества (U 1 ) и принадлежит эму. Для этого значения s точка
 ( s  ) является точкой бифуркации, поскольку при  ( s  ) имеется как минимум два
Следовательно, поскольку по условиям теоремы,
набора оптимальных законов, а при 0  s  s * один оптимальный закон - U . Теорема
доказана.
Литература
1.
2.
3.
4.
5.
6.
7.
8.
Ioffe A.D.., Tikhonov V.М. The theory of extreme tasks. M.: Science, 1974, (in Russian).
Bobylev N.A.,Emelianov S.V., Korovin S.К. Geometrical methods in variational tasks. M.: Magistr, 1998, (in
Russian).
Ashimov A.A., Borovsky Yu.V., Volobuyeva O.P., Ashimov As.A. On the choice of effective laws of
parametrical regulation of market economy mechanisms. “Automatics and telemechanics” journal, № 3, 2005,
105-112, (in Russian).
Ashimov A., Borovskiy Yu., Ashimov As. Parametrical Regulation of Market Economy Mechanisms // Proc.
18th International Conference on Systems Engineering ICSEng. 16-18 August, 2005. Las Vegas, Nevada, 2005,
189-193.
Kulekeev Zh., Borovskiy Yu., Ashimov A., Volobueva O. Methods of the parametrical regulation of market
economy mechanisms. Proceedings of the 15th international conference on systems science. 7-10 September
2004. V. 3. – Wroclaw: OWPW, 439-446.
Ashimov A., Borovskiy Yu.,Ashimov As. Parametrical Regulation Methods of the Market Economy
Mechanisms. Systems Science. Vol. 35, 2005. No. 1, 89-103.
Petrov А.А., Pospelov I.G., Shanain А.А. Experience of mathematical modeling of economy. M.:
Energatomizdat, 1996, (in Russian).
Ulam S. Unsolved mathematical tasks. M.: Science, 1964, (in Russian).
10
Download