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

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

ЛИТЕРАТУРА

1. Тарасевич Л.С., Гребенников П.И., Леусский А.И. Макроэкономика: Учебник. 6-е издание, испр. и доп. М.: Высшее образование, 2006. 654 с.

БЛАГОДАРНОСТИ: Работа выполнена при финансовой поддержке ЗАО «ПРОГНОЗ».

Поступила в редакцию 10 июня 2015 г.

Simonov P.M., Shults M.N. MODELING HOUSEHOLDS IN GENERAL EQUILIBRIUM THEORY

The article discusses the use of the general equilibrium theory in constructing and analyzing a model of the household functioning.

Key words: general equilibrium theory; computable general equilibrium model; households.

Симонов Петр Михайлович, Пермский государственный национальный исследовательский университет, г. Пермь, Российская Федерация, доктор физико-математических наук, профессор, профессор кафедры информационных систем и математических методов в экономике, e-mail: simpm@mail.ru

Simonov Pyotr Mihkailovich, Perm State National Research University, Perm, the Russian Federation, Doctor of Physical and Mathematical Sciences, Professor of the Information Systems and Mathematical Methods in Economics Department, e-mail: simpm@mail.ru

Шульц Михаил Николаевич, Пермский государственный национальный исследовательский университет, г. Пермь, Российская Федерация, аспирант кафедры информационных систем и математических методов в экономике, e-mail: mshults@mail.ru

Shults Mikhail Nikolaevich, Perm State National Research University, Perm, the Russian Federation, Post-graduate Student of the Information Systems and Mathematical Methods in Economics Department, e-mail: mshults@mail.ru

УДК 517.911.5

АППРОКСИМАЦИЯ ГИПЕРБОЛИЧЕСКИХ ДИФФЕРЕНЦИАЛЬНЫХ ВКЛЮЧЕНИЙ С ВНЕШНИМИ ВОЗМУЩЕНИЯМИ

© В.В. Скоморохов

Ключевые слова: гиперболические дифференциальные включения с импульсными воздействиями; аппроксимирующее отображение; радиус внешних возмущений; модуль непрерывности отображения; 6 -решение.

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

Пусть Мга — п -мерное векторное пространство с нормой | ■ | ; сотр[Мга] — множество всех непустых, компактов Мга; Л,[-, ■] — расстояние по Хаусдорфу между множествами, содержащимися в пространстве Мга . Обозначим Сп([а, Ь] х [с, й]) пространство непрерывных функций и: [а, Ь] х [с, й] — Мга с нормой ||и|| = тах{|и(£,ж)|: (¿, ж) € [а, Ь] х [с, й]} .

х [с, й]) — пространство суммируемых по Лебегу функций и: [а, Ь] х [с, й] —► Мга с ь й

нормой ||и|| = // |и(^, ж) | .

а с

Будем говорить, что ^: [а, Ь] х [с, й] х Мга — сотр[Мга] удовлетворяет условиям Кара-теодори, если выполняются следующие условия:

1) при каждом и € Мга отображение ^(■,-,и) измеримо;

2) при почти всех (¿, ж) € [а, Ь] х [с, й] отображение ^(¿, ж, ■) непрерывно;

3) для каждого ограниченного множества V С Кга найдется такая функция шу (■, ■) € € ^([а, Ь] х [с, й]), что при почти всех (¿, ж) € [а, Ь] х [с, й] и всех и € V выполняется неравенство (¿,ж,и)|| ^ шу(¿,ж).

Пусть € [а, Ь] (а < ¿1 < ... < ¿т < Ь) — конечный набор точек. Обозначим через Сп([а, Ь] х [с, й]) множество всех непрерывных на каждом из промежутков [а, ¿1] х [с, й] , (¿1, ¿2] х [с, й] , ..., (¿т,Ь] х [с, й] ограниченных функций и: [а,Ь] х [с, й] — Мга, имеющих пределы справа в точках , к = 1, 2,...,ш с нормой НжЦ^, П[аЬ] = sup{|u(t, ж)|: £ € [а, Ь] х х [с, й]} .

Рассмотрим задачу

д2и

€ ^ (¿, ж, и(£, ж)), (¿,ж) € [а, Ь] х [с,й], (1)

д£дж

А(и(^,ж)) = 4(и(4, ж)), к = 1,2,..., ш, (2)

и(£, 0) = а(г), и(0,ж) = в(ж), (3)

где а(-) , в(') , непрерывны и а(0) = в(0) , отображение ^: [а, Ь] х [с, й] х Мга — сотр[Мга] удовлетворяет условиям Каратеодори. Отображения : Мга — Мга , к = 1, 2,...,ш , непрерывны, А(и(^, ж)) = и(^ + 0, ж) — и(^, ж), к = 1,2,..., ш.

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

£ Х т

и(£,ж) = а(£) + в (ж) — а(0) + / д(«,т) йзйт + ^ ,ь]/& (и(4 ,ж)). (4)

ас

Обозначим через К ([а, Ь] х [с, й] х Мга х [0, то)) множество всех функций п: [а, Ь] х [с, й] х х Мга х [0, то) — [0, то), обладающих следующими свойствами:

1) при каждых (и, 5) € Мга х [0, то) функция п(', ',и, 5) измерима;

2) при почти всех (¿, ж) € [а, Ь] х [с, й] и всех 5 € [0, то) функция п(£, ж, ■, 5) непрерывна;

3) для каждых и € сотр[Мга] и 5 € [0, то) существует такая суммируемая функция : [а, Ь] х [с, й] — [0, то), что при почти всех (¿, ж) € [а, Ь] х [с, й] и всех и € и и т € [0, 5]

выполняется неравенство п(^,ж,и, т) ^ (¿, ж);

4) при почти всех (¿, ж) € [а, Ь] х [с, й] и каждого и € Кга выполняются равенства Ит п(^,ж,^,5)= п(^,ж,и, 0) = 0.

5^0+0

Обозначим через P([a, b] х [c, d] х R" х [0, то)) множество всех функций n: [a, b] х [c, d] х х R" х [0, то) ^ [0, то), обладающих свойствами из класса функций K([a, b] х [c, d] х R" х х [0, то)), а также удовлетворяющих следующим условиям: для каждых U € comp[R"] и 5 € (0, то) найдутся такие числа r(U, 5) > 0 и y(U, 5) ^ 0, что при почти всех (t, ж) € € [a, b] х [c, d] и всех u € U число r(U, 5) удовлетворяет неравенству r(U, 5) ^ n(t,x,u, 5), а для числа y(U, 5) при почти всех (t, ж) € [a, b] х [c, d] , всех u € U и т € [0, 5] имеет место оценка n(t, ж, u, т) ^ 7(U, 5).

Пусть ■, ■, ■) € P([a, b] х [c, d] х R" х [0, то)) . Определим функцию : [a, b] х [c, d] х х R" х [0, то) ^ [0, то) равенством

^(■0)(t, ж, u, 5) = sup h[F(t, ж, u), F(t, ж, v)]. (5)

Значения функции ■, ■, ■) в точке (t, ж,u, 5) будем называть модулем непрерыв-

ности отображения F: [a, b] х [c, d] х R" ^ comp[Rn] в точке (t, ж,u) по переменной u в шаре B[u, ^(t, ж, u, 5)], функцию ■, ■, ■) — функцией радиуса модуля непрерывности или просто радиусом непрерывности, а саму функцию ■, ■, ■) — функцией модуля непре-

рывности или просто модулем непрерывности отображения F: [a, b] х [c, d] хRra ^ comp[Rn] относительно радиуса непрерывности ■, ■, ■).

Будем говорить, что многозначное отображение F: [a, b] х [c, d] х R" х [0, то) ^ comp[Rn] аппроксимирует отображение F: [a, b] х [c, d] х R" ^ comp[Rn], если найдется такая функция £(•, ■, ■, ■) € K([a, b] х [c, d] х R" х [0, то)), что при почти всех (t, ж) € [a, b] х [c, d] и всех (u, 5) € R" х [0, то) выполняется оценка

h[F(t,ж,u),F(t,ж,u, 5)] ^ £(t,ж, u, 5). (6)

Отображение F(-, ■, ■, ■) будем называть аппроксимирующим отображение F(■, ■, ■) или просто аппроксимирующим. Функция £(•, ■, ■, ■) € K([a, b] х [c, d] х R" х [0, то)) в неравенстве (6) определяет степень близости значения F(t, ж, u, 5) в точке (t, ж, u) € [a, b] х [c, d] хR" к значению F(t^,u) для каждого фиксированного 5 € [0, то). Эту функцию {(•, ■, ■, ■) будем называть степенью аппроксимации отображения F: [a, b] х [c, d] х R" ^ comp[Rn] отображением F: [a, b] х [c, d] х R" х [0, то) ^ comp[Rn] или просто степенью аппроксимации.

Пару (F(-, ■, ■, •),{(•, ■, ■, ■)) будем называть аппроксимацией отображения F(■, ■, ■) или просто аппроксимацией, а если при почти всех (t, ж) € [a, b] х [c, d] и всех (u, 5) € R" х [0, то) выполняется включение F(t, ж, u) С FF(t, ж, u, 5) , то аппроксимацией вложением.

Значения аппроксимирующего отображения iF(-, ■, ■, ■) могут вычисляться с некоторой степенью точности, которую можно задать некоторой функцией n(', ■, ■, ■) € K([a, b] х [c, d] х х R" х [0, то)) .

В связи с этим рассмотрим отображение Qn: [a, b] х [c, d] х R" х [0, то) ^ comp[Rn], определенное равенством

Qn (t, ж, u, 5) = F(t, ж, u, 5)n(t'x'u's), (7)

где функция n(', ■, ■, ■) € K([a, b] х [c, d] х R" х [0, то)) в каждой точке (t, ж^) € [a, b] х [c, d] х хК" при каждом фиксированном 5 € [0, то) определяет погрешность вычисления значений аппроксимирующего отображения F(-, ■, ■, ■) . Далее, функцию n(',',',') будем называть радиусом внешних возмущений аппроксимирующего отображения F(-, ■, ■, ■) или просто радиусом внешних возмущений.

Пусть n(', ■, ■, ■) € K([a, b] х [c, d] х R" х [0, то)) . Рассмотрим при каждом фиксированном 5 € [0, то) дифференциальное включение

d2u

—- € Qn(t,ж,u(t,ж),5), (8)

где отображение : [a, b] х [c, d] х Rn х [0, то) ^ comp[Rn] задано равенством (7). Дифференциальное включение (8) будем называть аппроксимирующим дифференциальное включение (1) с внешними возмущениями.

Каждое решение включения (8) с импульсными воздействиями (2) и условиями (3) при фиксированном 5 > 0 будем называть 5 -решением (приближенным решением с точностью до 5 или просто приближенным решением) включения (1).

Пусть отображение F : [a, b] х [c, d] х Rn ^ comp[Rn] удовлетворяет условиям Каратео-дори. Рассмотрим задачу

д 2u

7—— € coF(i,x,u(i,x)), (9)

diax

A(u(ifc ,ж)) = 4 (u(ifc ,ж)), k = 1,2,..., m, (10)

u(t, 0) = a(t), u(0,x)= в(х), (11)

где co F(■, -,u(-, ■)) — выпуклая оболочка множества F(■, -,u(-, ■)) , отображения : Rn ^ ^ Rn , k = 1,2,..., m , непрерывны, A(u(tk, ж)) = u(tk + 0, ж) — u(tk, ж), k = 1,2,..., m , а(-),в(-) , непрерывны и a(0) = в(0) .

Обозначим через H(V), Hco(V) множества решений задач (1)-(3) и (9)-(11), соответственно, принадлежащих множеству V С Cn([a, b] х [c, d]) .

Пусть V С Cn([a, b] х [c, d]) , •, •, •) € K([a,b] х [c, d] х Rn х [0, то)) . Обозначим через Hrç(5)(V) множество всех решений задачи (8), (2), (3) с заданным радиусом внешних возмущений, принадлежащих множеству V .

Теорема. Пусть V — ограниченное замкнутое множество пространства С ([a, b] х [c, d]) и пусть ■, ■, ■) € P ([a, b] х [c, d] х Rn х [0, то)) . Далее, пусть пара (F(-, ■, ■, •),£(•, ■, ■, ■)) аппроксимирует отображения F(■, ■, ■) вложением. Тогда для любой функции п(', ■, ■, ■) € K ([a, b] х [c, d] х Rn х [0, то)), для которой существует такое число е > 0, что при почти всех (t, ж) € [a, b] х [c, d] , всех u € (U(V))£ и 5 € [0, то) имеет место неравенство

^(■0)(i, ж, u, 5) ^ n(t, ж, u, 5),

где ж, u, 5) — модуль непрерывности отображения F(■, ■, ■), выполняется соотно-

шение

Hco(V) = f| Hn(ô)(V*), ¿>0

где Hn(^)(Vô) — замыкание в пространстве Cn([a, b] х [c, d]) множества Hn(^)(Vô), Vô — замкнутая в пространстве Cn([a, b] х [c, d]) 5 -окрестность множества V.

ЛИТЕРАТУРА

1. Филиппов А.Ф. Дифференциальные уравнения с разрывной правой частью. М.: Наука, 1985.

2. Самойленко А.М., Перестюк Н.А. Дифференциальные уравнения с импульсным воздействием. К.: Вища шк., 1987.

3. Завалищин С.Т., Сесекин А.Н. Импульсные процессы. Модели и приложения. М.: Наука, 1991.

4. Булгаков А.И., Скоморохов В.В., Филиппова О.В. Асимптотические свойства множества S -решений функционально-дифференциального включения с импульсными воздействиями // Вестник Тамбовского университета. Серия Естественные и технические науки. Тамбов, 2011. Т. 16. Вып. 4. С. 1039-1043.

5. Булгаков А.И., Скоморохов В.В. Аппроксимация дифференциальных включений // Матем. сборник. 2002. Т. 193. № 2. С. 35-52.

Поступила в редакцию 9 июня 2015 г.

Skomorokhov V.V. APPROXIMATION OF HYPERBOLIC DIFFERENTIAL INCLUSIONS WITH IMPULSES AND WITH EXTERNAL PERTURBATIONS

In this paper, hyperbolic differential inclusions with external perturbations and with impulses are considered. Here we represent the concept of approximate solution ( S -solution) for a hyperbolic differential inclusion with impulses. The asymptotic properties of solutions sets to approximating differential inclusions with external disturbance are derived.

Key words: hyperbolic differential inclusions with impulses; approximating map; radius of external perturbations; modulus of continuity; S -solution.

Скоморохов Виктор Викторович, Тамбовский государственный технический университет, г. Тамбов, Российская Федерация, кандидат физико-математических наук, доцент кафедры высшей математики, e-mail: uaa@nnn.tstu.ru

Skomorokhov Viktor Viktorovich, Tambov State Technical University, Tambov, the Russian Federation, Candidate of Physics and Mathematics, Associate Professor of the Higher Mathematics Department, e-mail: uaa@nnn.tstu.ru

УДК 517.968

О ВЫЧИСЛЕНИИ ЯДРА ИНТЕГРАЛЬНОГО УРАВНЕНИЯ, ВОЗНИКАЮЩЕГО В ОБРАТНОЙ ГРАНИЧНОЙ ЗАДАЧЕ ТЕПЛОПРОВОДНОСТИ

© С.В. Солодуша, И.В. Мокрый

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

В статье рассматривается интегральное уравнение Вольтерра I рода типа свертки:

t

J KN(t - s)^(s)ds = y(t), 0 ^ s ^ t ^ T, (1)

0

N 1

KN(t - s) = £ (-1)q+1 q2e-n2q2(t-s), y(t) = go (t) , (2)

q=i

введенное в [1] в связи с поиском решения u(1,t) обратной граничной задачи

ut = uxx, ж € (0,1), t ^ 0,

u (ж, 0) = 0, u (0, t) = 0, ux (0, t) = go (t).

При численном решении (1), (2) возникают погрешности, связанные не только с погрешностью метода, но и с ошибками при выполнении операций машинной арифметики над вещественными числами с плавающей точкой.

Цель данной работы — разработать алгоритм вычисления Kn , учитывающий особенности машинной арифметики и обеспечивающий желаемое (заданное) число достоверных