ВЕСТНИК ПЕРМСКОГО УНИВЕРСИТЕТА
2012
ЭКОНОМИКА
Вып. 3(14)
РАЗДЕЛ II. ЭКОНОМИКО-МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ
УДК 517.929 + 330.4
КРАЕВЫЕ ЗАДА ЧИ ЭКОНОМИЧЕСКОЙ ДИНАМИКИ С ПРИБЛИЖЕННЫМ ВЫПОЛНЕНИЕМ КРАЕВЫХ УСЛОВИЙ. КОНСТРУКТИВНОЕ ИССЛЕДОВАНИЕ *
В.П. Максимов, д. физ.-мат. наук, проф. кафедры информационных систем и математических методов в экономике
Электронный адрес: [email protected]
А.Л. Чадов, ст. преп. кафедры информационных систем и математических методов в экономике
Электронный адрес: alchadov@yandex. т
Пермский государственный национальный исследовательский университет, 614990, г. Пермь, ул. Букирева, 15
Рассматриваются линейные краевые задачи для систем функционально'-дифференциальных уравнений с числом краевых условий, превышающим размерность системы. В экономической динамике краевые задачи связаны с исследованием достижимости заданных показателей функционирования экономической системы. Исследуется разрешимость таких задач в случае, когда допускается приближенное выполнение краевых условий. Предлагаемый подход использует теоремы, условия которых могут быть проверены с использованием современных средств вычислений.
Ключевые слова: функционально-дифференциальные уравнения; краевые задачи; конструктивные методы; доказательные вычисления.
Введение
Краевые задачи экономической динамики тесно связаны с реальными задачами экономики, в которых требуется дать ответ на вопрос о возможности достижения заданных значений некоторых целевых показателей развития экономической системы. При этом нередко число целевых условий, сформулированных на основе плановых или желаемых показателей, оказывается больше, чем это допускается математической теорией, требующей строгого согласования числа условий с размерностью соответствующей динамической модели. Обычно такая модель имеет вид системы обыкновенных дифференциальных уравнений, или, в более общем случае, системы функционально -дифференциальных уравнений, позволяющей учитывать при моделировании реальных процессов эффекты запаздывания и наличие внешних импульсных (шоковых) воздействий на систему [9, 10, 12 13, 15]. Краевая задача, т.е. система уравнений динамики и краевых условий, с числом краевых условий,
превышающим размерность пространства состояний системы без ограничений, называется переопределенной и требует для ее исследования специальных подходов. Один из таких подходов основан на расширении пространства состояний, позволяющем согласовать число краевых условий с размерностью расширенного пространства. Этот подход применительно к задачам экономической динамики подробно описан в работе [12]. Другой подход (см. [13]) использует идею ослабления ограничений, при котором происходит последовательное расширение множества допустимых состояний. В настоящей работе проблема переопределенности решается путем перехода к постановке краевой задачи, при которой допускается приближенное выполенение всех или некоторых краевых условий, заданных в виде равенств. При этом допустимый уровень погрешности может задаваться индивидуально для каждого краевого условия. Содержательный смысл приближенного выполнения краевых (целевых) условий представляется весьма естественным
*Работа выполнена при поддержке РФФИ, проект 10-01, и компании «Прогноз», г. Пермь. © Максимов В.П., Чадов А.Л., 2012
для реальных задач экономической динамики. Основной математический результат работы в этой части (теорема 1) опубликован без доказательства в кратком сообщении [14].
Исследованию краевых задач для систем дифференциальных уравнений и их обобщений посвящена обширная литература (см., например, [1, 2, 3, 19] и приводимую там библиографию). В этой работе речь идет о направлении исследований, связанном с теоретическим обоснованием и практической реализацией компьютерного (computer-assisted) исследования линейных краевых задач. Целью такого исследования является установление факта разрешимости краевой задачи и построение гарантированных оценок погрешности приближенных решений. Основу излагаемого здесь подхода составляют приемы приближенного описания множества решений линейного функционально-дифференциального уравнения с гарантированной оценкой погрешности в сочетании со специальными теоремами, условия которых могут быть проверены в результате доказательного вычислительного эксперимента, использующего современные компьютерные технологии и системы (Maple, Mathematica и др.). В случаях, когда известные достаточные признаки разрешимости краевой задачи оказываются неприменимыми, обсуждаемый подход может оказаться единственно возможным для получения результата [20]. Различные варианты этого подхода применительно к обыкновенным дифференциальным, интегральным уравнениям и уравнениям с частными производными занимают заметное место в современной литературе, начиная с основополагающей монографии Каучера и Миранкера [21]. Из числа недавних отметим работы [22-24].
1. Предварительные сведения
Мы ограничиваемся здесь краевыми задачами
Lx = f, lx = р (1)
с линейным ограниченным оператором L : ACn [0;T] ^ Ln[0;T] и линейным
ограниченным вектор-функционалом
l : ACn[0;T] ^ Rm . Здесь ACn[0;T] -
пространство абсолютно непрерывных функций
x : [0;T] ^ Rn , суммируемых z : [0,T] ^ Rn ,
по
= x(
где p|-max р,-1
Ln [0,7 ] - пространство Лебегу функций
... INL» = |0гКг)И’
для a = col(al,...,an)eR".
+ x
всюду предполагается, что главная часть оператора l , - оператор q = l v , где
t
{Vz) (t) = jz(s)ds , имеет представление
0
{Qz)(t) = z(t)--к(^Ф(^А te[0,T]. (2)
Здесь элементы k (t,s) ядра к(t,s) измеримы на множестве {(t,s):0 „ s „ t „ T} и таковы, что на этом множестве
K(t\ i,j = \,...,n,
где функция к суммируема на [0,t]. В этом случае функционально-дифференциальная система lж = / охватывает дифференциальные уравнения с сосредоточенным и/или распределенным запаздыванием и интегро-дифференциальные системы Вольтерра. Пространство всех решений однородной системы Lx = 0 имеет размерность n. Пусть базис в этом пространстве. Матрица
А' = (х]..xj называется фундаментальной
матрицей (для определенности будем считать, что ^(0) = E ,- единичная n х n -матрица). Задача Коши
Lx = f, x(0) = а однозначно разрешима при любых f е Ln [0,T] и а е R” и ее решение представимо в виде
л t
x(t) = X (t )a + jC(t, s) f (s)ds, (3)
где C(t, s) - матрица Коши.
Вектор-функционал l в задаче (1) имеет представление
rT
/x = xPx(0)+j Ф(s)x(s)ds,
(4)
Систематическое изложение теории краевых задач (1) дается в монографиях [1, 2, 8]. Ниже
где элементы т х п -матрицы ф измеримы и ограничены в существенном на [0,г ], т -постоянная т х п -матрица. Мы будем считать, ЧТО компоненты / г = \ , векТОр-
функционала / = со/(/1,...,/и) образуют линейно
независимую систему. Вектор-функционалы вида (4) исчерпывают класс линейных ограниченных вектор-функционалов,
определенных на АСп [0,т]. Краевые условия 1Х = р охватывают многочисленные классы конкретных краевых условий, встречающихся в приложениях, в том числе двух- и многоточечные, интегральные, нагруженные интегральные и др.
В случае т = п необходимое и достаточное условие однозначной
разрешимости задачи (1) — обратимость
матрицы IX = (/х,,..., 1хп) . Теория и реализуемая схема доказательного
x
AC
вычислительного эксперимента, ориентированного на исследование задачи (1), изложена в [1, 2, 8, 16], там же приведены конкретные
примеры эффективного исследования конкретных краевых задач. В основе реализованного подхода лежит известная теорема об обратном операторе [17, с.141], из которой следует, что если удастся найти такую обратимую п х п -матрицу г, для которой
IIх-г <|^> (5)
то матрица гх тоже обратима и, следовательно, задача (1) однозначно разрешима. При этом однозначно разрешимыми оказываются все краевые задачи с г и х, попадающими в указанную выше окрестность матрицы г. Этот факт играет принципиальное значение для установления факта однозначной разрешимости краевых задач с неточно заданными параметрами.
Оказалось, что матрицу г целесообразно искать в виде г = Тх“, где 1 : АСп[0,Т] ^ Яп - вектор-функционал,
близкий кг, а матрица ха со столбцами из АСп [0,Т] и свойством ха (0) = Е дает для
оператора Ь : АСи[0,Т] ^ !И[0,Т], близкого к ь,
см _
достаточно малую «невязку» А = ЬХа и, таким образом, служит приближением для фундаментальной матрицы х. Степень близости Т иг, ь и ь и степень малости а , гарантирующие однозначную разрешимость задачи (1), определяется специальными теоремами [1, с. 211-256], [2, с. 226-249]. При
этом фактическое построение матрицы гха и достоверная проверка неравенства (5) стали возможными с развитием современных компьютерных технологий. Эти технологии предъявляют определенные требования к операторам ь и Т — они должны принадлежать специальным классам
«вычислимых» операторов [1, с. 217].
2. Переопределенные краевые задачи и их е -приближенная разрешимость
В случае т > п задача (1) не может быть корректно разрешимой (т.е. всюду и однозначно разрешимой), необходимое и достаточное условие разрешимости такой задачи может быть записано как условие
ортогональности правой части {/
пространству решений однородной
сопряженной задачи ([1, с. 24]; [19, с. 11]). Таким образом, свойство существования точного решения переопределенной задачи является «тонким» (не грубым) свойством,
которое не может быть установлено в результате вычислительного эксперимента, оперирующего с приближенными данными и/или использующего вычисления с конечной точностью. Кроме того, в прикладных задачах, где краевая задача (1) возникает как модель реальных изучаемых процессов, упомянутое тонкое свойство либо указывает на неадекватность модели, либо приводит к необходимости изменить постановку задачи. Подход к преодолению проблемы переопределенности, связанный с расширением основного пространства и обобщением понятия решения, был предложен в [5], его систематическое изложение можно найти в [1,
2, 8]. Конструктивная реализация этого подхода подробно описана в [12]. Здесь мы используем другой подход.
С учетом того, что в любом случае на практике доступно для построения лишь приближенное решение (т.е. функция, дающая достаточно малую невязку при подстановке в уравнение и краевые условия), естественной представляется следующая постановка переопределенной краевой задачи (1).
Зафиксируем е = со1 ,..., ет }, ег. ... 0,
1 = 1 ■ Будем называть е -приближенным
решением краевой задачи (1) такое решение Х уравнения ь ж = / , что
т.е. краевые условия гХ = р выполняются приближенно с точностью, определяемой заданным вектором е .
Наша цель — сформулировать условия е -приближенной разрешимости в форме, позволяющей производить их проверку с помощью вычислительного эксперимента в ситуации, когда параметры краевых условий (матрицы т и ф ) могут быть заданы неточно, с известными оценками погрешностей. Для формулировки таких условий введем следующие обозначения.
Каждой матрице в, элементы которой могут принимать значения из заданных отрезков, поставим в соответствие пару матриц В и В, элементами которых являются соответственно левые и правые концы упомянутых отрезков. Через в обозначим множество матриц, элементами которых являются всевозможные сочетания левых и правых концов соответствующих отрезков. Значки . 5 - и ' при необходимости будем использовать применительно к элементам матрицы в. Так, например, Ь е[^ Ь ],
Ъ1]=\Ь1] Через ^(г) обозначим пхп-
матрицу п . Я(1) = лЧОехр||г(г)(/г|.
Пусть, далее, ха (/) - приближенная
фундаментальная матрица, ха (0) = Е ; А(г) - ее невязка: (ьха) (г) = А (г), г е [0,Т ];
г
Л(г) = |_А(г)] + Я (г)| |_А(г)]^г . Здесь и ниже для
0
матрицы В = Ь} через |_В ] обозначена
матрица {]/>, |}. Обозначим г (/) = ха(/) - \(/>,
¥{1)=х\1)+т, у со = {>-,(/)},
у (г)„ Уу(г)„ Уу(г), (е [0,Г];
Т г
© = |ф(5)7(s)ds , g(г) = |С(/, 5)/(5^5 ,
0 0
ь=Р-к.
Определим матрицу А равенством
А=Т+0.
Теорема 1. Пусть гэп^ = г для всех А = {а^} > а,, „ а ” а» . Пусть, далее, найдутся
последовательности индексов {/^, , /г} и
{Л > • • • > Л } м такая последовательность нулей и единиц {^...,1/.}; |/4е{0Д}, & = 1,...,г, чдао
г . ••• а . Ъ. +(-1)Г1е.
'1Л \}г '1 1 ; '1
| max
и
а ■ Ь. + (-1) 'е.
‘гіг ’ '•
(6)
| min
а . ■■■ а .
\3\ \)г
‘гА
а .
'гЛ
Тогда краевая задача (1) е -приближенно разрешима.
Доказательство. Представление
решения (3) сводит задачу об е -приближенной разрешимости к задаче о разрешимости системы линейных алгебраических неравенств. Действительно, применяя вектор-функционал г к обеим частям представления (3), получаем 1х = 1Ха +
и, таким образом, разрешимость краевой задачи разрешимость системы неравенств
е -приближенная — это
а. а -Ъ
^41
1=1
(7)
где А = Ц} = IX , Ъ, = Д - (%),..
Критерий разрешимости (7) дает теорема С.Н. Черникова [18, с. 66-70]: необходимым и достаточным условием совместности системы (7) ранга г >0 (г - ранг матрицы, составленной из коэффициентов при ее неизвестных) является существование в ее матрице такого отличного от нуля минора г -го порядка
Л, =
11
'1Л
и такого решения
системы
\и —Ъ. \ — є. , к = 1,...,г,
| гАг А- I 5 5 5
ЧТО ДЛЯ всех ; = 17 соотношение
а .
*1Л
|^т
удовлетворяется
а . м.
*Л гі
а і
1гк
а ■ и.
І / I
Ы г
Непосредственная проверка условий этой теоремы не представляется возможной. В нашем случае мы имеем приближенную информацию о параметрах задачи. При этом дело не только в неточном задании компонент вектор-функционала, но и в отсутствии точной фундаментальной матрицы х(г). Приводимые ниже соображения позволяют утверждать, что при выполнении условий (6) гарантируется выполнение условий теоремы С.Н.Черникова для любых возможных значений компонент а и
ь.
Построим сначала двухстороннюю оценку приближения производной Ха(/) фундаментальной матрицы по ее невязке а(г). Из (ьха) (г) = А(г) следует, что
X(() а -) = —, і)А(я^я .
Учитывая свойства
матрицы Коши С((,я) [7, 11], имеем
п
Ъ
а
а
У
Ъ
V
г
0
Щ-Xa(t) = A(t)+fc't(t,s)A(s)ds. Таким образом,
Ix(t)-xa (o| „ IA(0|+J0'|c;(f, s)\|a(5)|*. Как показано в [7, 11, 1, c. 242], матрица Коши C(t, s) связана с резольвентным ядром R(t, s) ядра к(t, s) соотношениями
R(t, s) = C',(t, s), 0 „ s „ t „ T. (В)
C(t, s) = E + J R(r, s)dr, 0 „ s „ t „ T. (9)
J s
Спектральный радиус интегрального оператора Вольтерра с ядром к(t, s) равен нулю [6, c. 153], поэтому
R(t, s) = Xj (/, s) + K2(t,s) н---------------1- (/, s) ч-,
где к (t, s) = к(t, s),
Kq (t,s) = f‘K(t, £,)Kq_, (£,,s)d£ q = 2,3,.... Покажем, что из неравенства
[ x(z)dz \Kq(t,s)\„ K(t) s
q-1
(q -1)!
следует неравенство
K+1(M)|„ *(0
Действительно,
f x(z)dz
j s
\^ic(z)dz
Js
q-1
(q -1)!
f (ri \\\y{z)dz
-d£ =
9-1
f tc(z)dz
J s
= к(0
9!
Отсюда
|Д(М)|„ £(/)ехр |£ АГ(Г)Й?Г |,
и для х(/), учитывая введенные обозначения, получаем
Ш)=Х\Г)-МГ)„ Х{Г)„ ха{Т)+мг) = пг).
Теперь, по построению, а9 „ а „ а»,
Ь:,, Ь:,, й, / = и все й и Ъ.
—; у у 7 ,/ у у у г
допускают эффективное вычисление.
Для завершения доказательства остается заметить, что из определения интервальных операций [4, с. 16] следует, что условия теоремы С.Н.Черникова выполнены для всех возможных значений параметров системы неравенств.
Приведем иллюстрирующий пример. Рассмотрим краевую задачу
x(t)+t[ (s-2)2x(s)ds = ——y(f\- — z(t\,
200 10
1 г> 1
y(t) = --x(t)-josy(s)ds + —z(/);
1 ft 1 1 t
z(t) H----f z(s)ds =-------x{t) H-----y(t)-------z(t);
10J° 10 100 9П w
t є [0,2].
2
J x(s)ds + y(1) + z(2) = 1;
x(1) + J y(s)ds + z(0) = 2;
2
x(0) + y(2) + J sz(s)ds = 3;
(•2 ,.2
J x(s)ds +J y(s)ds +
2
+J z(s)ds = 0
Первые три уравнения определяют динамику развития трех отраслей с состояниями х(г), у(г), z(t) соответственно. Отрасли оказывают взаимное влияние друг на друга, сила этого влияния определяется соответствующими коэффициентами, и обладают полной памятью об изменении своих состояний. Интегральные и точечные характеристики функционирования системы определены левыми частями четырех краевых условий. Требуется дать ответ на вопрос, достижимы ли значения этих показателей, заданные в правых частях, при заданном допустимом уровне погрешности.
В этом примере е -приближенная разрешимость установлена с помощью вычислительного эксперимента, реализованного в системе Мар1е, для
9 9 9 9
є = col(—, — , — ,—). 250 250 250 250
Матрица
используемая в процессе вычислительного эксперимента, имеет элементы
K(t ) =
t є [0,
10і
2t, t є [—,2].
10
Для представления о порядке коэффициентов а приведем их приближенные значения:
Ч
q
ч
au ' • au ' 1.659910494 0.8539527760 0.2829884432'
0.6987867571 1.427084461 0.8931947935
fl41 ' • й43 0.5638229580 0.005162564559 1.840289503
1.526003338 1.447018074 1.380024455
Список литературы
1. Азбелев Н.В., Максимов В.П.,
Рахматуллина Л.Ф. Элементы современной теории функционально-дифференциальных уравнений. М.: Ин-т компьютерных
исследований, 2002. 384 с.
2. Азбелев Н.В., Максимов В.П.,
Рахматуллина Л.Ф. Введение в теорию функционально -дифференциальных уравнений. М.: Наука, 1991. 280 с.
3. Азбелев Н.В., Максимов В.П.,
Симонов П.М.Функционально-дифференциальные уравнения и их приложения // Вестн. Удмуртского ун-та. Математика. Механика. Компьютерные науки. 2009. № 1. С. 3-23.
4. Алефельд Г., Херцбергер Ю. Введение в интервальные вычисления. М.: Мир, 1987. 356 с.
5. Анохин А.В. О линейных импульсных
системах для функционально-
дифференциальных уравнений // Доклады АН СССР. 1986. Т. 286. № 5. С. 1037-1040.
6. Забрейко П.П., Кошелев А.И., Красносельский М.А., Михлин С.Г. и др. Интегральные уравнения. М.: Наука, 1968. 448 с.
7. Максимов В.П. О формуле Коши для функционально-дифференциального уравнения // Дифференциальные уравнения. 1977. Т.13. №4. С. 601-606.
8. Максимов В.П. Арифметика
рациональных чисел и компьютерное
исследование интегральных уравнений // Соросовский образовательный журнал. 1999. № 3. С. 121-126.
9. Максимов В.П. Об одном подходе к задаче наведения системы в окрестность нормативной траектории // Вестн. Перм. ун-та. Серия: Экономика. 2008. № 8. С. 108-112.
10. Максимов В.П. Импульсная
коррекция управления для динамических
моделей с последействием // Вестн. Перм. унта. Серия: Экономика. 2009. № 1. С. 91- 95.
11. Максимов В.П., Рахматуллина Л.Ф. О представлении решений линейного функционально-дифференциального уравнения // Дифференциальные уравнения. 1973. Т. 9. № 6. С. 1026-1036.
12. Максимов В.П., Румянцев А.Н. Краевые задачи и задачи импульсного управления в экономической динамике. Конструктивное исследование // Известия
высших учебных заведений. Математика. 1993. № 5. С. 56-71.
13. Максимов В.П., Поносов Д.А., Чадов А.Л.Некоторые задачи экономикоматематического моделирования // Вестн. Перм. ун-та. Серия: Экономика. 2010. № 2. С. 45-50.
14. Максимов В.П., Чадов А.Л. О конструктивном исследовании краевых задач с приближенным выполнением краевых условий // Известия высших учебных заведений. Математика. 2010. № 10. С. 82-86.
15. Максимов В.П., Чадов А.Л. Гибридные модели в задачах экономической динамики // Вестн. Перм. ун-та. Серия: Экономика. 2011. № 2. С. 13-23.
16. Румянцев А.Н. Доказательный вычислительный эксперимент в исследовании краевых задач. Пермь: Изд-во Перм. ун-та, 1999. 172 с.
17. Треногин В.А. Функциональный анализ. М.: Наука, 1980. 496 с.
18. Черников С.Н. Линейные
неравенства. М.: Наука, 1968. 488 с.
19. Azbelev N.V., Maksimov V.P.,
Rakhmatullina L.F. Introduction to the theory of functional differential equations: Methods and applications. N.Y.: Hindawi Publishing
Corporation, 2007. 314 p.
20. Azbelev N.V., Maksimov V.P., Simonov P.M. Theory of functional differential equations and applications // International Journal of Pure and Applied Mathematics. 2011. V. 69. № 2. P. 203-235.
21. Kaucher E.W., Miranker W.L. Selfvalidating numerics for functional space problems. N. Y.: Academic Press, 1988. 256 p.
22. Kolev L.V. Outer interval solution of the eigenvalue problem under general form parametric dependences // Reliable Computing. 2006. V. 12. P.121-140.
23. Lin Y., Stadther M.A. Validated solutions of initial value problems for parametric ODE's // Applied Numerical Mathematics. 2007. V. 57. P. 1145-1162.
24. Nakao M.T., Hashimoto K., Watanabe Y. A numerical method to verify the invertibility of linear elliptic operators with applications to nonlinear problems // Computing. 2005. V. 75. P. 1-14.