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

Серия «Математика»

Том 1 (2007), № 1, С. 291-302

Онлайн-доступ к журналу: http://isu.ru/izvestia

ИЗВЕСТИЯ

Иркутского государственного университета

УДК 517.956

Обобщенное решение смешанной задачи для полулинейной гиперболической системы с суммируемыми по Лебегу входными данными *

В. А. Терлецкий (terletsky@math.isu.ru)

Иркутский государственный университет, г. Иркутск

Аннотация. Рассматривается гиперболическая система полулинейных (линейный дифференциальный оператор, нелинейная правая часть) дифференциальных уравнений со смешанными условиями достаточно общего вида. Предлагается сравнительно простой способ построения обобщенного решения.

Ключевые слова: обобщенное решение, гиперболическая система, интегральный эквивалент, характеристики

Введение

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

* Работа выполнена при финансовой поддержке РФФИ, проекты 05-01-00187, 0601-81016.

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

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

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

1. Постановка смешанной задачи

В прямоугольнике П = 5 х Т, Б = («0,51), Т = (¿о,^) рассмотрим систему дифференциальных уравнений с частными производными

Хг + А(в, ¿)Хз = / (х,в,г), (1.1)

в которой х = х(в,Ь) - искомое п-мерное вещественное решение, А = А(в, ¿) - заданная матричная функция размера п х п; / = /(х, в, ¿) -заданная п-мерная векторная функция.

Систему (1.1) будем считать гиперболической. Для этого, как известно, ([1], с. 23), следует предположить, что все собственные значения Х^ = Х^в,Ь), г = 1, 2,...,п матрицы А являются вещественными в П функциями с одинаковой алгебраической и геометрической кратностью. Последнее равносильно ([2], с.58) существованию в каждой точке (в,Ь) € П базиса в пространстве И"" из левых и правых

V = (р(1),р(2),...,р(п)) собственных векторов матрицы А со свойствами СА = ЛС, АР = VЛ, СР = Е. Здесь т - знак транспонирования, Л = diag{Хl, Х2,..., Хп}, Е = diag{1,1,..., 1}.

Дополнительно предположим, что все элементы матриц Л, С и V являются непрерывными и непрерывно дифференцируемыми функциями в замыкании П области П.

Перейдем к постановке смешанных условий для системы (1.1). С этой целью определим на границе 9П = П \ П прямоугольника П функции Хг = Хг(в,Ь), положив

Хг(в,Ь) = Уо(в,Ь) + Хг(в,Ь)г(в,Ь), (в,Ь) € ЙП, г = 1, 2,...,п,

где г>о = го(в,Ь), г = г(в,Ь) - компоненты вектора единичной внешней нормали к границе ЙП. За исключением угловых точек прямоугольника П они определяются однозначно: го = (—1)]+1, г = 0, в € Б,Ь = Ь]; го = 0; г = ( — 1)]+1, в = в], Ь € Т, ] = 0,1. По функциям Хг построим двузначные функции х±, положив

Х±(в, Ь) = sign хг(в, хг(в,Ь) ± 1]/2, г = 1, 2,...,п, (в,Ь) € ЙП.

Сконструируем диагональные матрицы 2- = diag{x-, Х-,..., Х-}, 2+ = diag {х+,х+, ... ,х+}, 2± = 2- + 2+. Тогда [3,4] для системы (1.1) будут корректными смешанные условия

г

2~(в,г)[Ст(в,г)х(в,г)-д(2+Стх, I ^(2±Стх,г)ёг,вл)] =0, (в,Ь) € ЙП,

го

(1.2)

в которых вектор-функция ц = д(х,у,в,Ь) задана в области И" х И1 х <9П.

Поясним смысл смешанных условий (1.2). Во-первых, отметим, что матрица 2- всюду на верхней границе прямоугольника П, т.е. при Ь = Ь1, равна нулю и, напротив, на нижней границе Ь = ¿о совпадает с единичной. Это означает, что решение х свободно от каких-либо условий на верхней границе П, но полностью закреплено связью х = Vц на нижней границе. На боковых границах в = во и в = 51 прямоугольника П картина более разнообразна. Здесь в точке (в,Ь) скалярное произведение < ¡(г\х > должно совпадать с значением функции цг, если Хг(в,Ь) > 0, и ничем не регламентируется в случае Хг(в,Ь) < 0, г = 1, 2,... ,п.

Во-вторых, подчеркнем, что первые два аргумента вектор-функции ц не существенны на нижней границе области П. Действительно, здесь влияние первого аргумента нейтрализуется тождеством 2+(в,Ьо) = 0, а второго - равенством Ь = ¿о. На боковых же границах они могут моделировать разнообразные функциональные и динамические связи для решения х. Первый аргумент 2+Стх формируется линейными комбинациями < ¡(г\х >, соответствующими отрицательным собственным значениям Хг, а второй аргумент позволяет подчинять линейные комбинации решения х обыкновенным дифференциальным уравнениям,

правая часть которых заведомо не зависит от координат образа Стх, отвечающих нулевым собственным значениям.

Укажем предположения на входные данные /, ц при которых будет строится обобщенное решение х задачи (1.1),(1.2).

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

\\/(х + Ах,в,г) - /(х,в,г)\\< ь\\Ах\\, (в,г) е п, \\д(х + Ах, у + Ау, г) - я(х,у,г)\\< ЩАх\\ + \\Ду\\), (1.4) \\ф(х + Ах,г) - ф(х,г)\\ < ь\\Ах\\, г е т,

для произвольных х, Ах е И™, у, Ау е И™ с константой Ь > 0.

При фиксированных значениях х е И™, у е И™ вектор-функции /, ц, и ф суммируемы по Лебегу с некоторой степенью р > 1 :

/(х, •) е ь™(П), д(х, у, •, •) е ь™(дп), ф(х, •) е ьт(т).

2. Интегральный эквивалент

Введем в рассмотрение инварианты Римана г = г(в,г), г(в,г) е И™, воспользовавшись линейным невырожденным преобразованием

г(в,г) = ст(в,г)х(в,г), х(в,г) = р(в,г)г(в,г). (2.1)

Если х - классическое (непрерывное и непрерывно-дифференцируемое в П) решение системы (1.1), то г удовлетворяет [3,4] инвариантной системе

гг + Лг3 = д(г,в,г), (2.2)

в которой д(г, в, г) = Ст (/(Рг, в, г) — (Рг + ЛР3)г).

Преимуществом инвариантной системы (2.2) по сравнению с исходной системой является простота ее дифференциального оператора. В каждое уравнение системы (2.2) ввиду диагональности матрицы Л входят производные только одного инварианта гг, г = 1, 2,..., п. Это обстоятельство позволяет рассматривать выражение г^ + \ггга как полную производную функции гг(в(г\£,т; г), г) по переменной г при условии, что функция в(г\£,т; г) является интегральной кривой обыкновенного дифференциального уравнения

ис

— = Хг(в,г), г = 1,2,..,п (2.3)

аг

со свойством в(г^(£,т; г) = £ при т = г.

Уравнения (2.3) принято называть характеристическими уравнениями системы (1.1), а интегральные кривые в(г\^,т;Ь) - характеристиками.

Пусть (§(г'),§(г')) = (§(г\в,г),§(г\в,г)) € 9П - начальная точка г-ой характеристики, выпущенной из точки (в, Ь) € П в обратном времени так, что Ь(г\в,Ь) < Ь. Тогда результатом интегрирования каждого дифференциального уравнения инвариантной системы (2.2) вдоль соответствующей характеристики будет система интегральных уравнений типа Вольтерра

г

Гг(в,Ь) = Гг(§ (г),Ь(г)) + I дг(г,в(г)(в,Ь,т ),т )(1т, г = 1, 2,...,п, (в, г) € П.

(2.4)

Здесь первое слагаемое определяется по правилу

(ф, 0,§(г),Ьо), §г)(в,Ь) = Ьо, г.(§(г) §(г)) = ^ №

Л ' ; ^±Стх(§(г),т),т)йт, §(гг)Л(г)), §(г)(в,г) >Ьо.

К го

(2.5)

В свою очередь, гладкие решения интегральной системы (2.4),(2.5) удовлетворяют [3,4] инвариантной системе (2.2).

Решение интегральной системы (2.4), (2.5) можно объявить обобщенным решением инвариантной системы (2.2) с соответствующими смешанными условиями, а обобщенным решением исходной задачи (1.1), (1.2) считать связанную с ним преобразованием (2.1) вектор-функцию х. Именно такая схема построения обобщенного решения использовалась в [3]. Однако, ничто не мешает рассматривать инвариантную систему (2.2) и вытекающую из нее интегральную систему (2.4), (2.5) как промежуточные результаты и пройти через них транзитом от задачи (1.1),(1.2) к ее интегральному варианту, имеющему вид

п „

х(в,1) = ^Р(г)(в,Ш(§(г)(в,Ь),§(г)(в,Ь))+ у дг(Стх, в(г)(вт),т)йт]. г=1 №(з,г)

(2.6)

Эквивалентность интегральной системы (2.5),(2.6) смешанной задаче (1.1)-(1.3) на гладких решениях и отсутствие в ней производных от решения х позволяет дать

Определение 1. Обобщенным решением задачи (1.1), (1.2) назовем вектор-функцию х, удовлетворяющую почти всюду в П интегральной системе (2.6) и почти всюду на 9П равенствам (2.5).

Идея доказательства существования и единственности обобщенного решения заключается в установлении сжимаемости отображения,

стоящего в правой части равенства (2.6). Обозначим его оператором X = X(д,ф,/; х). Тогда обобщенное решение есть неподвижная точка отображения X.

Существование и единственность такой точки х обоснуем методом последовательных приближений. Для этого построим последовательность функций {х(к^}, к = 0,1,..., положив

х(к+1) = X (д,ф,/; х(к)), (2.7)

начиная с х(0 = 0. Понятно, что для сходимости последовательных

приближений достаточно доказать сходимость ряда ^ 5(к), если поло-

к=0

жить

¿к+1)(<1 г)=\\\х(к+1)(в,г) -х(к)(в,г)\\, (в,г) е п, 0 \\\г±(в,г)(х(к+1\в,г) -х(к)(в,г))\\, (в,г) е дп.

Отправной точкой здесь служит неравенство

5(к+1)(в,г) < к]Т || 5(к)( в(г)(в, г; т),т)Ит+

г=1 Ь« (з,г) ™

+ I [¿)z+(s(■г),i(■г))6(k)(в(j)(s(■г)(в,г),г(г)(в,г);т),т)+ го 1=1

+ 5(к)(^(г)(в, г); т)]ат}, (2.8)

являющееся следствием условий Липшица (1.4) с константой к > 0, зависящей только от константы Ь и матрицы А. Отметим, что неравенство (2.8) справедливо лишь для "малых" интервалов Т, т.к. предполагает наличие не более, чем одного "отражения" характеристик от боковых границ. Для этого достаточно считать, что

г1 - г0 < ( тах тах \Аг(в, г)\)_1(в1 - во). г=1,2,...,™ (з,г)&и

Однако, такое ограничение не обременительно. Его всегда можно выполнить, разбив интервал т на "малые" интервалы. Понятно, что построение обобщенного решения на каждом последующем интервале проводится в точности так же, как и на самом первом.

3. Построение обобщенного решения строго гипреболической

системы

Упростим рассматриваемую задачу, предположив, что система (1.1) является строго гиперболической ([5], с.13), или гиперболической в узком смысле ([1], с.26, [6], с.77). Тогда существует число й > 0, при котором справедливы неравенства

\Хг(в,г) - хз (в,г)\ > й> о, г = э, г,э = 1,2,...,п, (в,г) е П. (3.1)

В принципе, сделанного предположения достаточно для получения всех результатов данного параграфа. Но, стремясь максимально просто и выпукло изложить сжимающее свойство отображения X и тем самым обосновать сходимость метода последовательных приближений (2.7), дополнительно упростим и смешанные условия (1.2). Будем считать вектор-функцию ц зависящей только от переменных (в, Ь) е 9П. Фактически это означает, что смешанное условие (1.2) превращается в обычное условие Коши для соответствующего инварианта Римана, которое задано в начальной точке характеристики, находящейся на боковой границе области П. Тогда неравенство (2.8) примет более простой вид

п £

5(к+1)(в,Ь) < К]Т / 5(к)( в(г)(в,Ь; т),т)йт, к = 0,1,..., (в,Ь) е П.

г=1т{з,г)

(3.2)

Применив здесь однократное рекуррентное усиление, будем иметь при к = 1, 2,..., (в,Ь) е П неравенства

п п £ т

6(к+1) (в,Ь) < К2 £ ]Т I ( ! 6(к-1)(в(з)((,т; а),а)йа )\^Чз,;т) йт. г=1 з=1№(з,г)

(3.3)

Последующие подстановки оценки (3.2) в неравенство (3.3), по-видимому, лишены смысла, т.к. приводят лишь к наращиванию кратности интегралов в нем. Однако, преодолеть эту трудность позволяет свойство

(3.1).

Заменим переменную т на новую переменную (, положив при г = Э, г,э = 1, 2,... ,п

( = в(з)(в(г)(в, Ь; т),т; а). (3.4)

Тогда

й( = [вТз)(С, т; а) + в3 (С, т; а)вУ (в, Ь; т)]\5=^(з,р,а)йт. Воспользуемся здесь тождеством [3, с.10] в(тз)(С,т; а) + Хз (С,т )в[з)(£,т; а) = 0, э = 1, 2,...,к, (С,т) е П, а е Т

и уравнением (2.3) для характеристики £ = в(г\в,Ь; т). Получим равенство

Лт = [(Ш, т) - А,(£, т))в(](£, т; а)]-1Л(. (3.5)

Оно показывает, что для непрерывности и гладкости замены перемен-ных(3.4) свойство (3.1) является необходимым условием. Второе такое необходимое условие, а именно, отделимость от нуля частной производной вЯ)(£, т; а) при заданных предположениях на собственные значения Аз матрицы А системы (1.1) выполняется всегда. В [3] показано, что существует константа Ь > 0, обеспечивающая при всех ] = 1, 2,... ,п и всех (£,т) € П, (в(з\£,т; а), а) € П неравенства

¿¿\£„т; а) > Ь. (3.6)

Нетрудно видеть, что новая переменная ( при фиксированном значении а изменяется на интервале

со{в(г)(в,Ь; т),в(з) (в, г; т)}П 5.

Далее для определенности будем считать, что

Аг > Хз и (в(г)(в,Ь; т),в(з)(в,Ь; т)) С 5

при всех т € (Ьо, Ь). Тогда в силу (3.4)-(3.6) при всех г,] = 1, 2,... ,п, г = ] справедливы соотношения

} ( } 6(к-1)(в(з)(£,т; а),а)Ла ^т =

№)(з,г) №)(£,Т)

г з^^Ца) -

= Г Г 1>а,а)йСйа < (3 _)

г з(Я(з, г;а)

< 1 I I $(к-1)(СаЖЛа.

№)(з,г) 8(*)(з,г;а)

Любой треугольник Огя(в,Ь), служащий областью интегрирования в правой части неравенств (3.7), всегда входит в область С(в,Ь) определенности решения х в точке (в,Ь), которая строится, как известно [1], по правилу

С(в, Ь) = {(£, т) € П : втф, Ь; т) < £ < втах(в, Ь; т), т € Т}.

Здесь кривые £ = вт1П(в, Ь; т) и £ = втах(в,Ь; т) являются интегральными кривыми дифференциальных уравнений

Л£ Г. Аг(£,т), £ € 5

Лт 1 0, £ = во, £ = в1,

dZ f. min Ai(£>T), Z e S,

dT l 0, Z — so, Z — si

соответственно.

Вследствие указанного свойства области G и неотрицательности подынтегральной функции справедливы неравенства

t sL)(s,t;a)

j j 5(k-i)(Z,a)dZda < jj 5(k~l)(Z,T)dZdT (3.8)

i(i)(s,t) s(i)(s,t;a) G(s,t)

при всех i — j, i,j — 1, 2,... ,n. Неравенства (3.8) устанавливают единообразную оценку для всех двойных интегралов из неравенства (3.3) при i — j, i, j — 1, 2, . . . , n.

Обратимся теперь к случаю i — j. Тогда в силу очевидного свойства s(i (s(i) (s, t; t),t; a) — s(i\s,t; a) двойной интеграл превращается в повторный по интервалу (t(i')(s, t),t). Нетрудно показать, что для любого номера j — 1, 2,... выполняется тождество

I f (t - a)j-iö(k-j) (s(i) (s, t; a),a)dadT <

i(i) (s,t) m(s,t) (,339)

< } U-OLö(k-j)(s(i)(s,t; a),a)da.

t(i">(s,t)

На основании формул (3.7)-(3.9) неравенство (3.3) при к = 1 может быть преобразовано к виду

п £

5(2)(в,Ь) < К2[£ / 5(0)(в(г)(в,Ь; т),т)йт+

г=1 №(з,£) (3.10)

+I 5(0)(С,тШт].

о(з,г)

Описанная схема позволяет "подниматься" по номерам к = 2,3,...

в неравенстве (3.2), опираясь на тождество (3.9) и оценку (3.10). На произвольном шаге к будем иметь

П £ /< \к — 1

5(к)(в,Ь) < Кк [ Е I ^—-ГТ 5(0)(в(г)(в,Ь; т),т)йт+

г=Ч0)(з,£)2 (3.11)

+I *(0) (С,тШт].

В свою очередь, приращение 5° удовлетворяет оценке

5(0)(в,Ь) = \\^(1)(в,Ь)\\ < К1Е[\\д(0, 0,ё (г)(в,Ь),№(в,Щ\ +

£ г=1 (3.12)

+ / \\/(0,в(г)(в,Ь; т),т)\\йт].

Здесь константа К1 > 0 зависит только от матрицы А системы (1.1). Пусть вектор-функция х есть предел последовательности {хк}, или,

что то же х(в,Ь) = (х(к+1^(в, Ь) — х(к)(в,Ь)). Тогда, очевидно,

k=0

11®(м)П<£ t(k)(s,t).

k=0

Оценив сумму ряда с помощью неравенств (3.11),(3.12), получим итоговую оценку для решения х в точках П в виде

г

НхМК К2^2[\\д(0,0,§ (г)(в,Ь),§(г)(в,т\ + I\\/(0,в(г)(в,Ь; т),т)\\йт] +

г=1 №(з, г)

+ I \\д(0,0,£,т)\\йш + Ц \\1(0,£,т)\\й£йт). (3.13)

9ппда(8,г) а(в,ь)

Здесь константа К2 > 0 зависит только от констант К, К1,п, а, Ь, Ь1 — Ьо, а интеграл по участку границы дП П дС является интегралом первого рода, т.е. йш = \]йт2 + й£2.

Нетрудно показать, что функция, стоящая в правой части неравенства (3.13) суммируема по Лебегу в области П со степенью р и имеет суммируемые в 5 со степенью р следы при любом фиксированном Ь € Т. Поэтому обобщенное решение х существует и единственно в ^(П) и удовлетворяет оценке роста (3.13).

4. Построение обобщенного решения при существенно ограниченных входных данных

Как мы видели, результаты предыдущего параграфа существенно используют предположение (3.1) о строгой гиперболичности системы (1.1). Здесь мы откажемся от ограничения (3.1), заменив его на условие существенной ограниченности входных данных задачи (1.1),(1.2). А именно, будем считать, что для любых x € Rn, y € Rm

f (x, ■, ■) € L^U), q(x, y, ■, ■) € 1&(дП), ф, ■) € L™(T). (4.1)

В этом случае все последовательные приближения xk, k = 1,2,..., построенные методом (2.7) будут также существенно ограниченными в П, что позволяет ввести в рассмотрение функции

w(k)(t) = ess sup 5(k) (s,t), определенные на интервале T при всех k = 0,1,....

Из неравенства (2.8) следует, что

г

из(к+1) (Ь) < К^ и)(к) (т)йт, (4.2)

го

где К3 = К(п2/2 + 2).

Условия (4.1) и (4.2) обеспечивают существование и единственность обобщенного решения х, удовлетворяющего оценке

\\х\\ЬПоо{п) < К4(М0, •)\\Ьгп{т) + Ш0,0, •, (От) + \\1 (0, •, •)\\ь^0(П))

с некоторой константой К4, зависящей только от К\, К3 и (Ь\ — Ьо).

Таким образом, для существенно ограниченных входных данных задачи (1.1), (1.2) построение обобщенного решения выглядит намного проще, чем в случае входных данных лишь суммируемых по Лебегу с некоторой степенью р > 1. Более того, здесь не требуется отделимости всюду в замкнутом прямоугольнике П любых двух собственных значений Хг и Хз, г = ], г,] = 1, 2,...,п, т.е. не используется условие

(3.1).

5. Обоснование существования обобщенного решения в

общем случае

Займемся построение обобщенного решения задачи (1.1),(1.2), считая выполненными только предположения первого параграфа.

Введем в рассмотрение последовательности вектор-функций {/к}, {д(к^} и {у(к} к = 1, 2,..., определив их с помощью покоординатных срезок величины к соответствующих вектор-функций:

(к)(„ в +) = \ fi(x,в,Ь), \Мх,в,Ь)\ < k,

Цк)(х,8,Ь) =

к 8щи^(х,в,Ь), \^(х,в,Ь)\ > к, х € И"", (в,Ь) € П;

( Яг(х,у,в,Ь), \дг(х,у,в,Ь)\ < к, д\ (х, у, в, Ь) = < к 8щщг(х,у,в,Ь), \ дг(х,у,в,Ь) \ > к, х € И",у € И™,

I (в,Ь) € Ш;

г = 1,2,...,п;

щ(к)(х Ь) = 1 Уз(x,Ь), \(х,Ь) \ < k,

ух,Ь) \ к вщи^з (х,Ь), \ щ (х,Ь) \ > к, х € И", Ь € Т; ] = 1, 2,...,ш.

Понятно, что при любом к = 1, 2,... так построенные вектор-функции f(к), д(к) и у(к) являются существенно ограниченными в области своего определения, т.е. удовлетворяют условию (4.1). Следовательно, для любого к = 1, 2,... существует единственное обобщенное решение х(к € Ь"(П) задачи (1.1),(1.2) с входными данными f(к, д(к, у(к).

Покажем, что последовательность {хк} является фундаментальной в ^(П). Пусть символ А(к'1 означает разность элементов соответствующей последовательности с номерами к и I. В частности, А(к'1^х = х(к) — х(1\ А(к,1) / = /— /и т.п. По определению обобщенного решения имеем равенства х(к) = X(ц(к),ф(к), /(к); х(к)), х(1) = X(ц(1), ф(1), /(1); х(1)), в которых оператор X строится по правилу (2.6). Отсюда после несложных, но громоздких процедур можно получить оценку

\\А(к>1)х(;Щь^) < К5[(}(\\Ак1)х(,т)\\ь^) +

£о

+ Е \\г±(вз,т)А(к'1)х(вз,т)\\)йт + \\Ак)/(х(к), •, он™) + 3=0

+\\Ак1)ф(к), ОН^сг) + \\Ак1)д(х(к),у(к), •, ОН^ш)]

с константой К5 той же природы, что и предыдущие. Из данного неравенства и фундаментальности последовательностей {/(к^}, {ф^} и {д(к^}, следует фундаментальность последовательности {х(к^}, а значит и наличие предельного элемента х этой последовательности, являющегося обобщенным решением задачи (1.1),(1.2).

Список литературы

1. Рождественский Б.Л., Яненко Н.А. Системы квазилинейных уравнений и их приложения к газовой динамике. — М.: Наука, 1978. — 687 с.

2. Ланкастер П. Теория матриц. — М.: Наука, 1982. —270 с.

3. Терлецкий В.А. Обобщенное решение гиперболических систем одномерных полулинейных дифференциальных уравнений // Серия: Оптимизация и управление . Вып. 11. — Иркутск: Изд-во Иркут. ун-та, 2004. — 48 с.

4. Терлецкий В.А. Обобщенное решение одномерных полулинейных гиперболических систем со смешанными условиями // Изв. Вузов. Математика. — 2004. — № 12. — С. 75-83.

5. Куликовский А.Г., Погорелов Н.В., Семенов А.Ю. Математические вопросы численного решения гиперболических систем уравнений. — М.: ФИЗМАТЛИТ, 2001. — 608 с.

6. Петровский И.Г. Лекции об уравнениях с частными производными. — М.: ФИЗМАТГИЗ, 1961. —303 с.

V. A. Terletsky

Generalized solution of a mixed problem for semi-linear hyperbolic system with Lebesgue summable input data

Abstract. A hyperbolic system of semilinear differential equations with boundary conditions of a very general type is considered and a rather simple way of constructing its generalized solution is proposed. Key-words: generalized solution, hyperbolic system, integral equivalent, characteristics.