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

МАТЕМАТИКА

УДК 517.957; 532.526.2

В. В. Феоктистов, О. О. Мякинник

СТРУКТУРА РЯДА ДЛЯ РЕШЕНИЯ СИСТЕМЫ УРАВНЕНИЙ В ЧАСТНЫХ ПРОИЗВОДНЫХ 1-ГО ПОРЯДКА

Предложена структура ряда с неопределенными числовыми коэффициентами для представления решения линейной и нелинейной систем уравнений в частных производных 1-го порядка и связанных с этими системами задач. Для линейной системы доказана теорема о разложении аналитического решения задачи Коши по бегущим волнам с матричными коэффициентами. E-mail: olga.mknk@ipmtel.ru

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

В 1968г. Ф.И.Федоров (Институт физики АН БССР) записал в матричной форме и предложил к изучению систему нелинейных дифференциальных уравнений в частных производных 1-го порядка, нелинейные члены которой представляют собой квадратичную форму [1]:

Ак — = Со + С В + БТС2Б. (1)

к= 1

Число уравнений системы (1) равно числу компонент искомой п-мерной вектор-функции В независимых переменных хк, к = 1, в:

Вт = (В1(Х1, Х2, . . . , х3), В2(Х1, Х2, . . . , х3), ..., Вп(х 1, Х2, . . . , Хs)) .

Коэффициенты Ак, к = 1, в, представляют собой квадратные числовые матрицы размера п х п с элементами а^. Каждый элемент п-мерного столбца в правой части системы представляет собой полином 2-й степени относительно компонент неизвестной вектор-функции В: числовой вектор С0 составлен из свободных членов уравнений системы с0г}, г = 1, п; числовая квадратная матрица С1 размера п х п содержит коэффициенты при неизвестных функциях В^, г,] = 1,п; числовая кубическая матрица С2 размера п х п х п состоит из коэффициентов 0г2т при квадратичных членах В^В], г,],1 = 1,п. Таким образом, г-е уравнение системы (1) имеет вид

£ Е дХ = 4"+Е +Е Е c2{j!> в в-

k=i j=i k j=i j=i 1=1

Поясним структуру матрицы С2. Кубическую матрицу размера п х х п х п можно представить как книгу из п страниц, каждой из которых соответствует квадратная матрица размера п х п. Умножив вектор-строку Вт справа на кубическую матрицу С2, получим квадратную матрицу. Каждая ¿-я строка этой матрицы является результатом умножения Вт на г-ю "страницу", г = 1 ,п. Умножив затем квадратную матрицу ВтС2 справа на вектор-столбец В, получим столбец ВтС2В, г-й элемент которого имеет вид j=1¿^ 1=1 с2 Вj В;.

Свойства системы (1) определяются алгебраическими свойствами задающих систему числовых матриц-коэффициентов [1]. Матрицы могут быть вырожденными, состоять только из нулевых элементов или содержать нулевые строки, что позволяет рассматривать в рамках одной системы как нелинейные, так и линейные дифференциальные и алгебраические уравнения.

Выделим некоторые важные частные случаи системы (1). Это, прежде всего, система

5 О п

V Лк — = Со + СД (2)

к= дхк

которая определяется в работах самого Ф.И. Федорова и В.Я. Скоро-богатько (Институт прикладных проблем механики и математики АН УССР) как универсальная форма линейных дифференциальных уравнений [1, 2]. Действительно, при условии равенства числа неизвестных функций и числа уравнений система линейных дифференциальных уравнений с постоянными коэффициентами любого порядка может быть приведена к виду (2) при помощи введения в рассмотрение дополнительных функций.

При п =1 и 8 = 1 система (1) представляет собой дифференциальное уравнение Риккати

Л— = Со + С1В + С2В2, п = 1, 8 = 1. ах

Поставив задачу отыскания решений системы (1), которые удовлетворяют условию

дВ = 0, з = М, (3)

дхк

получим систему квадратичных алгебраических уравнений

п п п

о = с?} + V с{гЛBj + V V с?0 BjВ;, г = М, (4)

¿=1 j=l ;=1

относительно неизвестного числового вектора В. В работе [3] решения системы (1), удовлетворяющие условию (3), выделены в класс эле-

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

Система (1) допускает случай, когда все матрицы Лк, к = 1,8, задающие левую (дифференциальную) часть, являются нулевыми (содержат только нулевые элементы). Тогда система (1) является алгебраической системой 2-го порядка (4) и имеет только элементарные решения.

К виду (1) можно привести любую систему дифференциальных уравнений с постоянными коэффициентами, нелинейные члены которой представляют собой произведение неизвестных функций и их производных. Этому требованию удовлетворяют многие уравнения математической физики. Система (1) получила название системы Ф.И. Федорова (доказавшего фундаментальное значение этой системы в теоретической физике) или системы универсальных нелинейных уравнений.

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

Задача Коши для системы (1). Предположим, что хотя бы одна из матриц-коэффициентов Лк, к = 1,8, задающих левую часть системы (1), является обратимой. Пусть

Введем обозначение х3 = £ для независимого переменного х3, связанного с выделенной матрицей Л3. Умножим систему (1) слева на -Л-1 и, вернувшись после выполнения этой операции к прежним обозначениям, получим систему уравнений

где I — единичная матрица размера n х n, A — линейный дифферен-циалный оператор. Будем рассматривать t как временное переменное, а xk, k = 1, (s — 1), как пространственные переменные. Система (5) представляет собой систему нелинейных нестационарных уравнений в матричной форме и является важным частным случаем системы (1) [2].

Вместе с тем запись (5) представляет собой нормальную форму системы дифференциальных уравнений в частных производных, или систему Ковалевской [4].

det Ak = 0 при k = s.

(5)

Для системы (5) рассмотрим начальное условие по переменному

В(£, х) ^ = 0(х), X = (хЬ . . . , £,5-1), 0т(х) = (^1(х), . . . , ^п(х)) ,

где вектор-функция 0(х) — аналитическая функция всех своих аргументов в некоторой области С0 гиперплоскости £ = 0, содержащей точку х0, х0 = (х1,..., х<0-1). Таким образом, задача Коши (5)-(6) для системы Ф.И. Федорова (1) в случае ее приведения к нормальной форме (5) удовлетворяет условиям теоремы Коши-Ковалевской, следовательно, задача (5)-(6) имеет единственное аналитическое решение в некоторой достаточно малой области С, С Э С0, пространства переменных (£, х1,..., х5-1) [5].

При исследовании системы Ф.И. Федорова для сохранения ее матричной формы В.Я. Скоробогатько предложил метод, названный им матричным [2]. В частности, решение задачи (5)-(6) в классе функций Сж представлено им как решение нелинейного интегрального уравнения

г

В = ехр (Л£) 0(х) - У ехр Л(£ - в) (<Со + С1В + 0тС20) (7) 0

При этом первое слагаемое интегрального уравнения (7) удовлетворяет как линейному дифференциальному уравнению

дВ=

которое является частным случаем нелинейного уравнения (5), так и начальному условию (6) [2].

Рассмотрим задачу Коши (5)-(6) при нулевых матрицах С1 , С2 и нулевом векторе С0. Докажем теорему, из которой следует представление выражения ехр (Л£) 0(х), являющегося решением задачи, в виде ряда специального вида. Теорема 1. Ряд

ж

%х) = ^ {Х^Ха ••• Х,--1} • 7«, х = (хь...,х,_1), (8)

а, ||а|| =0

в_1

Хк = Рхк + £Лк, а = («1,..., а,_1), ||а|| = ^ ак,

к=1

где В — п-мерная вектор-функция, 7а — числовой п-мерный вектор-коэффициент, Лк, к = 1, (в — 1), — квадратные числовые матрицы размера п х п, I — единичная матрица, является частным решением

системы уравнений в частных производных 1-го порядка

дБ ^ „ дВ

= Е Лк эх? <9>

к=1

Прежде чем перейти к доказательству теоремы, дадим пояснения, относящиеся к построению ряда (8). Суммирование в (8) проводится по мультииндексу а размерности (з — 1). Напомним, что мультиин-декс представляет собой числовой вектор заданной размерности, компоненты которого принимают значения из множества N0. Сумма этих значений называется длиной мультииндекса а и обозначается ||а||. Выпишем, например, все двумерные мультииндексы а, а = (а1,а2), длина которых равна трем:

(3, 0), (2, 1), (1, 2), (0, 3).

Рассмотрим отдельные сомножители Хк, Хк = /хк + £Лк, заключенные в фигурные скобки выражения (8). В скалярном случае (п =1, Лк — числовые коэффициенты) верно утверждение: гладкая функция В(£, х) удовлетворяет уравнению в частных производных 1-го порядка вида (9) тогда и только тогда, когда функция В(£, х) является первым интегралом системы обыкновенных дифференциальных уравнений

^ = —Лк, к = М—Г. (10)

Выражения

Хк = /хк + £Лк, к = 1, (з — 1), (11)

представляют собой (в — 1) независимых первых интегралов системы (10) и являются решениями линейного однородного уравнения (9), имеющими форму плоской бегущей волны [4].

Запишем систему обыкновенных дифференциальных уравнений (10) в симметричной форме, предположив, что (9) является системой (п > 2), задаваемой матрицами Лк :

/¿х1 №2 1

„ „ , = —■(12)

Л1 Л2 Л5_1

По аналогии со скалярным случаем рассмотрим выражения (11) для линейной однородной системы уравнений в частных производных 1-го порядка (9) (п > 2) и связанной с ней системы обыкновенных дифференциальных уравнений с матричными коэффициентами (12). Очевидно, что при формальной замене скалярных коэффициентов на матричные использование бегущих волн (11) для изучения решения системы (9) связано с проблемой некоммутативности умножения матриц.

В качестве возможного подхода к данной проблеме В.Я. Скоробо-гатько в рамках матричного метода ввел операцию симметризации [2]. В записи (8) при формулировке доказываемой теоремы эта операция обозначена фигурными скобками и означает, что следует взять сумму всех возможных существенных произведений, построенных из заключенных в скобки сомножителей. Пусть, например, а является двумерным мультииндексом длины 3 с компонентами а1 = 2 и а2 = 1. Тогда из двух сомножителей Х^ и Х2, стоящих в фигурных скобках соответствующего члена ряда (8), можно составить три различных произведения:

{х2Х2} = (Х1)2Х2 + Х1Х2Х1 + Х2(Х1)2.

Таким образом, ряд (8) представляет собой разложение искомой вектор-функции В(£, х) по соответствующим мультииндексу а, а = (а1,..., ав-1), степеням Х^, к = 1,8 — 1, выражений вида (11), к результату некоммутативного умножения которых при каждом а применена операция симметризации.

Описательное определение операции симметризации, данное в работе [2], неудобно при проведении преобразований. Построим алгоритм применения этой операции к произведению р сомножителей, которые имеют вид Х^. Представим результат симметризации как сумму р слагаемых, в каждое из которых сгруппированы члены, имеющие первым левым сомножителем элемент Хк, к = 1,р. При этом порядок следования остальных сомножителей Х^, ] = к, является произвольным, что соответствует определению операции симметризации.

Тогда, вынося Хк, к = 1,р, за скобки, получаем

р

{Х?1 Х2*2 • • • Хра^ = ^ Хк • {Х?1 • • • Х--1 • • • Хра^ , ак > 1. (13)

к=1

Применим алгоритм (13) для понижения порядка заключенных в фигурные скобки сомножителей. На некотором шаге под знаком операции симметризации получим произведение т сомножителей первой степени, где т < р. Следуя алгоритму далее, будем последовательно уменьшать число сомножителей в фигурных скобках до полного их раскрытия:

{Х1Х2 • • • Хт} =

т т т т т

= ^2 П Хк2\ = ^2 Хк1 ^2 Хк^ П Хк^ =

к1 = 1 к2 = к1 = 1 к2 = кз=1

к2=к1 к2=к1 к3=к2=к1

т т т

= • • • = ^ Хк1 Хк2 • • • Хкт-1 • Хкт . (14)

к1 = 1 к2=1 кт-1 = 1

к2 = к1 кт-1=кт-2 = -=к1

Построенный алгоритм (13)-(14) может быть рассмотрен как определение оператора (операторной скобки) на элементах счетного множества, для которых введены операции умножения и сложения.

Поставим в соответствие системе линейных уравнений в частных производных 1-го порядка (9) множество К_ь элементами которого являются матрицы Лк, задающие эту систему, и построенные на их основе бегущие волны Хк:

Ak е Vs-1, Xk е Vs-1, Xk = Jxfe + tAfe, k = 1, s - 1. (15)

Введем на элементах множества Vs-1 действие оператора W по мультииндексу а, а = (а1, ... , ар), размерности (s — 1):

W(Xk)=Xk, k = p = 1, s — 1;

W (Xf1, X2a2,. ..,Xa )= У Xk W (Xf1,

k=1

X afc-1

(16)

,Xpap); (17)

p

£

k=1

W (XbX2,...,Xp) = > Xk W (X1,...,Xk-1,Xk+1,...,Xp). (18)

Сотношения (16)-(18) определяют некоторый закон взаимодействия бегущих волн. Назовем введенный оператор Ш оператором (операторной скобкой) волнового взаимодействия или волновой Ш-скобкой по мультииндексу а размерности р. Длину мультииндекса а назовем порядком волнового взаимодействия (или порядком Ш-скобки). Соотношение (17) позволяет понижать порядок оператора, соотношение (18) связывает операторы разной размерности. При этом для множеств р =1,8 — 1, на которых рассматривается волновая Ш-скобка по мультииндексу а размерности р, имеет место соотношение

V! С V с •••с V, С •••С К_1. (19)

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

Для аргументов скобки волнового взаимодейстия Ш имеют место следующие перестановочные свойства:

Лк (ХкГ = (Хк Г Лк; (20)

Лк Хг + ЛгХк = ХгЛк + Хк л, к = г, к,г = 1,8 — 1. (21)

Свойство (20) является очевидным, поскольку матрица Лк — перестановочная с (Лк )р, единичной матрицей I и скалярными сомножителями хк. Для проверки свойства (21) нужно раскрыть скобки в

выражениях Хк и вынести матрицы Ак и А правым сомножителем. Свойство (21) означает, что некоммутативные матрицы, задающие систему уравнений (9), коммутируют по сумме с бегущими волнами, построенными на основе этих матриц. Отметим, что

дХк дХк дХк т

- — Ак, - — - — I,

д£ дхк дХк

и рассмотрим операцию дифференцирования оператора волнового взаимодействия Ш. Из (18) непосредственно следует, что

д

-Ш (Х1,...,Хк,..., Хр ) =

р

= ^ Ш(Х1,... ,Хк-1, Ак, Хк+1,..., Хр) =

к=1

р

= р ^ Ак Ш (Х1,...,Хк-1,Хк+1,...,Хр), р =178—1. (22)

к=1

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

Рассмотрим дифференцирование оператора волнового взаимодействия Ш по его аргументу. Выполнение этой операции связано с понижением порядка Ш:

д

Ш (Ха ,...,ХГ ,...,Ха) =

= ||а||Ш (ХГ1 ,...,ХкГк-1 ,...,Ха), р =1,8 — 1. (23)

Доказательство правила дифференцирования (23) приведено в ходе доказательства теоремы.

Полученные правила дифференцирования оператора волнового взаимодействия Ш и перестановочные свойства его аргументов (20)-(21) используются в ходе доказательства теоремы.

Доказательство теоремы. Введем обозначение Ра для общего члена ряда (8):

-Г = Ш(ХГ1,..., Ха-1) • 7а, а = (а1,..., а3_1). (24)

Волновая Ш-скобка в (24) рассматривается на связанном с системой (9) множестве (см. (15)), обладающем свойством (19).

Очевидно, что при ||а|| =0 вектор Ра представляет собой числовой вектор Р0, -ро = I•То, и является свободным членом ряда (8). При этом

вектор Р0 является элементарным решением системы (9), содержащей произвольное число (в — 1) пространственных переменных хк.

Пусть (в —1)=1. Тогда система (9) содержит одно пространственное переменное х1 и, соответственно, одну матрицу Л1:

дБ дВ

1дВ = Л ах? (25)

В этом случае отсутствует проблема некоммутативности умножения матриц. Действие Ш-скобки по одномерному мультииндексу а, а = а1, на элементах множества V является, согласно определению, действием единичного оператора. При этом ряд (8) принимает вид

те те те

В(*, х) = ^ Рв1 = ^ Ш(Х?1) • 7«! = ^ (1х1 + ¿Л^1 • 7«!. (26)

а1>0 а1>0 а1>0

Непосредственная подстановка вектор-функции РГ1 в систему (25) показывает, что каждый отдельный член ряда (26) является решением этой системы. В силу линейности системы и аналитичности рассматриваемых функций ряд (26) и любая его частичная сумма также являются решением системы (25).

Рассмотрим случай (в — 1) = 2. Тогда система (9) принимает вид

г аВ , дВ , дВ

1дв = Л1 ахг+Л2 ах2. (27)

Для построении решения системы (27) в качестве аргументов Ш-скобки по двумерному мультииндексу используются элементы множества У^, V э У1.

Выделим сначала те члены ряда (8), которые содержат только одно пространственное переменное:

Р(а1, 0) = Ш(Ха1 )^7(„1, 0), Р(0,а2) = Ш(Х2*2 )^7(0,а2), а1, а2 > 0.

Проверка того, что векторы Р(Г1, 0) и Р(0,а2), зависящие от одного пространственного переменного, являются решением (27), совпадает с соответствующей проверкой при (в — 1) = 1 для элементов множества У1.

Рассмотрим теперь те члены ряда (8), которые содержат оба пространственных переменных. Тогда (24) принимает вид

Ра = Р(а1,а2) = Ш(Х?1 ,Х2*2) • 7(^2), аг = 0, г = 1, 2, (28)

где действие оператора Ш-скобка по двумерному мультииндексу не является тривиальным.

Методом математической индукции докажем, что вектор Ра с любыми ненулевыми значениями компонент двумерного мультииндек-

са а является решением системы (27). В качестве параметра индукции выберем порядок волнового взаимодействия | а| , соответствующего вектору -р(Г1 ,12). Выполним первый шаг математической индукции для ||а|| = 2, т.е. для случая, когда в (28) порядок оператора Ш равен его размерности:

||а|| =2, -(1,1) = Ш(Х1,Х2) • 7(1,1). Вычислив производные вектор-функции -р(1,1) по правилам (22) и (23):

д-р(1,1)

dt

= 2AiX2Y(i,i) + 2A2XiY(i,i),

дР(1,1) = 2Х;7 дР(1,1) =2Х;7 —- = 2Х27(1,1)^ —- = ^^Д^

дх1 дх2

убеждаемся, что -р(1,1) является решением системы (27).

Пусть утверждение математической индукции верно, когда параметр индукции ||а|| равен т, где т > 2. Это означает, что вектор-функция, имеющая вид (28), где а1 + а2 = т, также является решением системы (27). Шаг индукции от т к (т + 1) соответствует переходу к оператору волнового взаимодействия большего порядка и осуществляется увеличением на единицу степени ак, к = 1, 2, одного из двух его аргументов. Пусть это аргумент Х1:

-р(а1+1,а2) = Ш(ХГ + \ХГ ) • 7(п+1,а2), (29)

||а|| = (а1 + 1) + а2 = т +1.

Покажем, что из утверждения "вектор-функция -р(а1,Г2) является решением системы (27)" следует утверждение "вектор-функция -р(а1+1,Г2) является решением системы (27)".

В соответствии с введенным определением рассмотрим результат действия Ш-скобки порядка (т +1) по двумерному мультииндексу как сумму двух слагаемых, содержащих Ш-скобки порядка т по мульти-индексу такой же размерности:

-р(а1+1,а2) = Х1Ш(ХГ1 ,Х2"2 ) • 7(п+1,Г2) +

+ Х2Ш(Х1а1+1,Х2Г2_1) • 7(а1 + 1,«2). (30)

Введем обозначение

Ц/т = Ш(ХГ1 ,Х212) • 7(11+1,12), Кт = Ш(ХГ1+1,Х2а2_1) • 7(11 + 1,12).

(31)

Следовательно, вектор (29) имеет вид

-р(11 + 1,12) = Х1ит + X2рm, (32)

где для вектор-функций ит и утверждение математической индукции выполнено. Вычислим для построенной вектор-функции Р(а1+1,а2) составляющие системы (27):

А дР(д1+1,а2) = + АЛ ^дХт + аЛ дд^т;

A2

d-P(oi + 1,02)

d—(oi + 1,02)

= A2X1 ^ + ^Km + A2X

dx

dVm .

dx '

ÖUn

- = А1С/т + ^ + А2Кт + Х2 .

Воспользуемся в полученных выражениях свойством (20) для перестановки матрицы А и одноименной волны Х^ и свойством коммутативности по сумме (21) для того, чтобы выразить произведения А^Х? через X? А (г = ^, г,^ = 1, 2). Подставим полученный результат в систему (27):

d-

(01+1,02)

= X1

^ A1 d-P(oi+1,o2)

dx1

dUm dUm 1 dx1

dt

dvm A dUm 1 dx1

dt

"Ao

д-U

+ X2

+ t (A1A2 - A2A1)

dUm 1 д^2

dvm

dvm dUm

(01+1,02) dX2

+

+

dx1 dx2

(33)

Применив правило дифференцирования (23) к вектор-функциям ?7„ и Р™, заданным выражениями (31), получим

dvm öu„

= («1+«2)w (ха1 ,x2a2-1).

dx1 dx2

(34)

Тогда из выражения (33) следует утверждение: если вектор-функции Ц7т и являются решением системы (27), то вектор-функция -Р(о1+1,о2) также является решением системы (27). Утверждение математической индукции для случая (в — 1) = 2 доказано.

По построению рассмотренные вектор-функции -Р(о,о), —(а1,0), .Р(о,а2) и Р(а1 ,а2) исчерпывают все возможные члены ряда (8) при (в — 1) = 2. Таким образом, для случая (в — 1) = 2 доказано, что каждый член ряда (8), представляющий собой результат взаимодействия волн, является решением системы (27). В силу линейности системы и

2

аналитичности рассматриваемых функций ряд

те те

р = Е -(11,12) = Е Ш(Х11 ,Х212) • 7« (35)

||1||>о 1Н1>о

и любая его частичная сумма также являются решением системы (27). Утверждение теоремы для (в — 1) = 2 доказано.

Доказательство утверждения теоремы для системы (9), содержащей произвольное фиксированное число (в — 1), (в — 1) > 3, пространственных переменных, проводится аналогично. Убедимся, что каждый член ряда (8), имеющий вид (24), является решением системы (9). Пусть в (24) мультииндекс а размерности (в — 1) таков, что содержит р ненулевых компонент. Перенумеруем ненулевые компоненты мультииндекса в порядке их следования при произвольном положении (в — 1 — р) нулевых компонент:

а = (0,..., а1, 0,..., 0, а2, 0,..., 0, ар, 0,..., 0), р < (в — 1).

Тогда вектор (24) принимает вид

-1 = Ш(Х11 ,...,Хр1р) • 71, р = 1,8 — 1. (36)

Пусть р < (в — 1), т.е. число ненулевых компонент мультииндекса меньше размерности соответствующего этому мультииндексу оператора волнового взаимодействия Ш. Тогда проверка того, что вектор (36) является решением системы (9), повторяет рассуждение, проведенное для вектора, построенного при помощи Ш-скобки по мультииндексу меньшей размерности.

Рассмотрим, например, для системы (9), содержащей три пространственные переменные, член ряда (8), соответствующий мультииндексу (а1,а2, 0). Это означает, что р = 2, а волна (!ж3 + ¿А3) не входит в число аргументов Ш-скобки по трехмерному мультииндексу. Тогда производные вектора (36) по х3 обращаются в нуль и случай сводится к рассмотренным ранее случаям действия Ш-скобки на элементах множеств V и V и соответствующих им систем (25) и (27).

Пусть р = (в — 1), т.е. при построении вектора (36) используются волны Хк, содержащие все пространственные переменные системы (9). Доказательство, как и в случае (в — 1) = 2, проведем методом математической индукции с параметром индукции, равным порядку волнового взаимодействия. Первый шаг индукции выполним при ||а|| = р (все компоненты мультииндекса а равны единице), применив к вектор-функции (36) правила дифференцирования (22) и (23). Пусть является решением системы (9) при ||а|| = т, т > р. Пусть далее (как и в (29)) шаг индукции от т к (т + 1) осуществляется

увеличением на единицу значения аь Мультииндекс с измененной первой компонентой обозначим через а':

= а1 + 1, а'к = ак, к = 2, р.

Рекуррентное соотношение, аналогичное соотношению (32), получим в соответствии с (17), последовательно группируя в одно слагаемое все члены ряда, имеющие первым левым сомножителем волну Хк по каждому хк, к = 1,р:

-Ра' —

£ XkРk, Р1 — w(ха1,..., xpp)7а

k=1

Рk — W (X°

ai + 1

Xак -1

..хар )Ya', k — 2,p.

(37)

При этом для функций ик, к = 1,р, связанных по построению с мультииндексом длины т, утверждение индукции выполнено.

Рассмотрим систему (9) при В = Ра/. Воспользовавшись после вычисления Ак • дРа/ /джк, к = 1,р, свойствами (20) и (21), получим

д-Ра dt

S- 1

д-

2> ^ = £ X

k=1

k=1

д Р

k

dt

£ Ak

дРk

p

+ t > v (AkAj - Aj Ak)

k,j=1

k=1 др j др k

дж

+

дж

k

дж,

p < (s - 1). (38)

Применив правило дифференцирования (23) к функциям ик, заданным формулой (37), получим

дик м = гр, к=

дРj

джк

дж,

Из (38) следует утверждение: если функции ик являются решением системы (9), то и функция Ра/ является решением системы (9). Утверждение индукции доказано для произвольного (в — 1), 8 > 3.

Таким образом, для произвольного 8, 8 > 3, рассмотрены все возможные мультииндексы а, по которым выполняется суммирование в ряде (8), и соответствующие им члены ряда, каждый из которых является решением системы (27). В силу линейности системы и аналитичности рассматриваемых функций ряд (8) и любая его частичная сумма, которая представляет собой нелинейную комбинацию бегущих волн, также являются решением системы (27). ■

Замечание. Выведем для двумерного мультииндекса правило дифференцирования Ш-скобки волнового взаимодейстия (23), которое было использовано в ходе доказательства теоремы. Применим метод математической индукции с параметром индукции, равным порядку вол-

нового взаимодействия. Для случая а1 = а2 = 1 утверждение индукции принимает вид

(Х1,Х2) = 2Х2.

дХ1

Выполним шаг индукции от (а1 + а2) к (а1+а2 + 1), воспользовавшись свойством (17) для вынесения (внесения) аргументов за (под) знак W-скобки:

д

—w (х1в1+1,хаа ) =

д

= [Х^(Х?1 ,Ха) + X2W(Х1а1+1,Х2аг-1Я = дХ1

дд

= W (Х?1 ,Х2"2) + Х^-—W (Х^1 ,Х2"2) + Х^-—W (Х1а1+1,Х2а2-1) = 1 2 дХ1 1 2 дХ1 1 2

= W(Х^1 ,Х2а2)+(«1+«2) Х^(Х?1-1^2)+X2W(Х^1 ,Х2"2-1)] = = W (Х^1 ,Х2а2) + («!+а2^ (Х^1 ,Х2а2) = («1+«2+1^ (Х^1 ,Х2а2).

Полученное соотношение доказывает правило дифференцирования (23) для двумерного мультииндекса а.

Следствие 1. Поскольку общий член ряда (8) представляет собой результат действия W-скобки, заданной произвольным мультииндек-сом а произвольной фиксированной размерности (в — 1), и удовлетворяет матричному уравнению (9), это уравнение можно рассматривать как правило дифференцирования по £ оператора W-скобка, который задан мультииндексом а размерности (в — 1):

д -Г-а,. "-аа_1\

(Х- ,...,Х^ ,...,Х5а!!1 ) =

в-1 в-1

= ||аП £ АW(ХГ1,..., Х^-1,..., Х5а-_1), ||а|| = £ ак. (39)

к=1 к=1

Правая часть формулы (39) получена в результате применения в правой части уравнения (9) формулы (23), которая задает правило дифференцирования W-скобки по пространственному переменному. Использованная при доказательстве теоремы формула (22) является частным случаем (39) для ак = 1, к = 1, в — 1.

Следствие 2. Пусть 1 < р < д < (в — 1). Рассмотрим системы дифференциальных уравнений в частных производных 1-го порядка

4(40

и

di f-f dxk

Тогда каждый член ряда (8), построенного для системы (40), и любая частичная сумма этого ряда представляют собой и соответствующие элементы ряда (8), построенного для системы (41), являясь при этом решениями обеих систем.

Следствие 3. Из доказанной теоремы непосредственно следует, что ряд (8) с коэффициентами Ya есть решение поставленной задачи Коши, если

1 dai+a2+-----ha3_i

Ya = й! dxa...ас! ^ о• (42)

II II 1 2 S—1 ж=ж0

Следовательно, ряд (8) c коэффициентами (42) является сходящимся в некоторой достаточно малой области G, содержащей точку x0, и представляет собой иную форму записи ряда, фигурирующего в теореме Коши-Ковалевской.

Общий случай. Пусть система дифференциальных уравнений в частных производных 1-го порядка (1) с произвольными коэффициентами получена в результате рассмотрения некоторой линейной или нелинейной модели математической физики. К таким моделям можно отнести, например, систему уравнений Максвелла, системы уравнений газовой динамики и гидродинамики.

Некоторыми исследователями ряды различной структуры предлагаются в качестве асимптотического метода аналитического исследования не только линейных, но и нелинейных задач. В качестве примера приведем серию работ А.Ф. Сидорова (Институт математики и механики УО РАН) [14, 15]. Вместе с тем общие методы построения рядов сразу для широкого класса нелинейных уравнений и связанных с ними краевых задач отсутствуют.

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

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

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

В работе [6] предложена структура ряда с неопределенными числовыми коэффициентами, задающего решение системы уравнений (1):

В(х) = Е Та• ха, ха = х?1 х^2 ...х?', 7Г = (т!, ...,-£) . (43)

а, ||а||=0

Суммирование в ряде (43), записанном относительно произведения в независимых переменных системы, проводится по мультииндексу а размерности в (ср. с рядом (8)). Коэффициент 7а представляет собой неопределенный числовой вектор, верхний индекс при компонентах которого указывает на соответствующую компоненту вектора В:

Р(х)т =(B1(x),...,Bn(x))= ( £ т1 • xa,..., £ тП • ха

. а, liai =0 а, ||а|| =0

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

Рассмотрим уравнение

£ A

д , ж ,

■ ж - ха1 ха2...xas i = с + c^o+ттc^, (44)

^дх I / ^ 'а ^

fe=1

.0.

получающееся в результате сравнения коэффициентов при х для ||а|| = 0. Левая часть уравнения (44) представляет собой полином 1-й степени относительно Та, где ||а|| = 1, правая часть этого уравнения — полином 2-й степени относительно 70. Запишем уравнение (44)

в виде соотношения

1Л(7«, а : ||а|| = 1) = 2Д(Та, а : ||а|| = 0), (45)

где полиномы 1-й и 2-й степени обозначены соответственно символами 1Л и 2Д. В круглых скобках указана длина мультииндекса векторных коэффициентов 7а, относительно которых эти полиномы записаны.

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

1 Д(7а, а : ||а|| = 2) = 2Д(7а, а : ||а|| = 0,1); (46)

1 Л(7а, а : ||а|| = 3) = 2Д(7а, а : ||а|| = 0, 2);

1 Л(7а, а : ||а|| = N — 1) = 2Д(7а, а : ||а|| = 0,Х — 2); (47) 1 Л(7а, а : ||а|| = N) = 2Д(7а, а : ||а|| = 0^ — 1); (48) 1 Л(7а, а : ||а|| = N + 1) = 2Д(7а, а : ||а|| = 0^); (49) ..........................................................................................(50)

В линейной части каждого уравнения системы все коэффициенты отвечают мультииндексам одной длины /, где I > 1. В нелинейной правой части содержится коэффициент с мультииндексом длины (/ — 1) и произведения коэффициентов, сумма индексов которых не превышает (/ — 1).

Отметим, что если в виде ряда (43) представить решение системы (2), то алгебраическая система для определения коэффициентов разложения имеет вид

1Л(7а, а : ||а|| = I + 1) = ^(т«, а : ||а|| = /), I = 0,1,..., то. (51)

Таким образом, если для линейного случая первый коэффициент разложения 7а, ||а|| = 0, содержится в качестве неизвестного только в первом уравнении бесконечной линейной алгебраической системы (51), то в нелинейном случае — во всех уравнениях бесконечной алгебраической системы 2-го порядка (45)-(50).

В результате представления решения системы дифференциальных уравнений в частных производных 1 -го порядка с квадратичной нелинейностью (1) в виде ряда (43) с неопределенными числовыми коэффициентами получена бесконечная алгебраическая система 2-го порядка (45)-(50) для нахождения векторных коэффициентов 7а. В отличие от линейного случая при увеличении числа рассматриваемых уравнений от N до ^ + 1) появляются не только коэффициенты 7а,

||а|| = N + 1, но и меняются значения коэффициентов 7а, ||а|| = 0,Х. Это означает, что последовательное рекуррентно точное [14] нахождение коэффициентов ряда становится невозможным, поскольку уравнения для их определения являются взаимозависимыми.

Допустим, что, начиная с мультииндекса а определенной длины, все векторы 7а нулевые. Тогда верно следующее утверждение.

Утверждение 1. Пусть для коэффициентов ряда (43), задающего решение системы дифференциальных уравнений в частных производных 1-го порядка с квадратичной нелинейностью (1), выполнено условие

= 0, Уа ||а|| > N. (52)

Тогда в системе алгебраических уравнений 2-го порядка относительно коэффициентов этого ряда, получаемой при подстановке (43) в (1), число уравнений равно числу неизвестных.

Действительно, для каждого искомого коэффициента 7а с фиксированным мультииндексом а, 0 < ||а|| < N — 1, существует уравнение, которое содержит этот коэффициент в правой части в следующей форме: во-первых, линейно — как результат умножения на матрицу С1; во-вторых, в виде произведений 70 т • 7а и 7а т • т0 — как результат построения квадратичной формы ВтС2В. Все другие искомые коэффициенты в правой части этого уравнения отвечают мультииндексам меньшей длины, а все искомые коэффициенты в левой части — муль-тииндексам большей длины. В предположении 7а = 0 при ||а|| > N в уравнениях вида (49) и левая, и правая часть обращаются в нуль. В уравнениях вида (48) правая часть обладает описанным свойством относительно 7а, ||а|| = N — 1. При этом, поскольку левая часть таких уравнений обращается в нуль, число искомых коэффициентов не увеличивается. Следовательно, в системе уравнений (45)-(48), выделенной из бесконечной системы (45)-(50) при помощи условия (52), число уравнений равно числу неизвестных 7а, ||а|| = N — 1.

Допустим, что при ограничении (52) на коэффициенты разложения ряда (43) алгебраическая система 2-го порядка решена тем или иным методом. Тогда известны коэфффициенты разложения 7а, ||а|| < N. Построенная частичная сумма ряда (43) может быть рассмотрена как представление частного решения системы (1).

В качестве одного из способов разрешения нелинейности в системе уравнений (1) и связанной с ней алгебраической системе 2-го порядка (45)-(48) предлагаем использовать постановку нелинейных краевых условий. Такой подход реализован в работах авторов [8-11] на модели уравнений пограничного слоя [12, 13]. Для полученной системы уравнений в частных производных 1-го порядка с квадратичной

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

Таким образом, для нелинейной системы (1) и ее частного случая — линейной системы (9) построены ряды (8) и (43) с числовыми векторными коэффициентами. Ряды предложенной структуры можно использовать для аналитического исследования моделей математической физики.

СПИСОК ЛИТЕРАТУРЫ

1. Федоров Ф. И. Уравнения первого порядка для гравитационного поля // ДАН СССР. - 1968. - T. 179, № 4. - C. 802-805.

2. Скоробогатько В. Я. Решение системы дифференциальных уравнений матричным методом // ДАН УССР. - 1988. T. 3. - С. 28-31.

3. М я к и н н и к О. О. Нахождение элементарных решений системы уравнений типа Ф.И. Федорова // Докл. Акад. наук Украины. - 1994. T. 11. - С. 45-50.

4. К у р а н т Р. Уравнения с частными производными. - М.: Мир, 1964. - 829 с.

5. Петровский И. Г. Лекции об уравнениях с частными производными. -М.: Физматлит, 1961. -400 с.

6. Skorobogat'koV., MyakinnikO. Ona Power Series Representation of the General Solution of Fedorov's Set of Equations // Gravitation and Cosmology. - 1995. - Vol. 1. -No. 4. -P. 315-318.

7. Скоробогатько В. Я., Мякинник О. О. О решении задачи Коши для системы уравнений Ф.И. Федорова // Тез. докл. конф. Современные методы нелинейного анализа. 25-29 апреля 1995, Воронеж.

8. Ф е о к т и с т о в В. В., М я к и н н и к О. О. Решение краевой задачи для нелинейной системы, сводящейся к системе с квадратичной нелинейностью // Вестник МГТУ Сер. "Естественные науки". - 2000. - № 2 (5). - С. 3-9.

9. Феоктистов В. В., Мякинник О. О. Реализация квадратичной нелинейности в уравнениях с частными производными // Материалы конф. Современные методы теории функций и смежные проблемы. 26 января-2 февраля 2003. - Воронеж, 2003. - С. 267-268.

10. Ф е о к т и с т о в В. В., М я к и н н и к О. О. О решении задач пограничного слоя, преобразованных к системе уравнений в частных производных первого порядка с квадратичной нелинейностью // Вестник МГТУ. Сер. "Естественные науки". - 2004. - № 1 (12). - С. 54-71.

11. Феоктистов В. В., Мякинник О. О. Специальные ряды для решения систем уравнений в частных производных 1-го порядка // Современные проблемы прикладной математики и математического моделирования: Материалы 3-й Междун. научн. конф., Воронеж, 2-7 февраля 2009 г. - Воронеж: ВГУ, 2009.

12. ШлихтингГ. Теория пограничного слоя. - М.: Наука, 1969. - 744 с.

13. О л е й н и к О. А., С а м о х и н В. Н. Математические методы в теории пограничного слоя. - М.: Наука. Физматлит, 1997. - 512 с.

14. С и д о р о в А. Ф. О некоторых аналитических представлениях решений нелинейного уравнения нестационарной фильтрации // В сб. Сидоров А.Ф. Избранные труды: Математика. Механика. - М.: Физматлит, 2001. - С. 281-288.

15. Сидоров А. Ф. Специальные конструкции рядов для решения нелинейных уравнений с частными производными // В сб. Сидоров А.Ф. Избранные труды: Математика. Механика. - М.: Физматлит, 2001. - С. 217-225.

Статья поступила в редакцию 20.03.2009

Владимир Васильевич Феоктистов родился в 1944 г., окончил в 1967 г. МГУ им. М.В. Ломоносова. Д-р техн. наук, профессор кафедры "Математическое моделирование" МГТУ им. Н.Э. Баумана. Академик Российской академии космонавтики им. К.Э. Циолковского, лауреат премии им. Н.Е. Жуковского РАН. Автор более 70 научных работ в области прикладной математики.

V.V. Feoktistov (b. 1944) graduated from the Lomonosov Moscow State University in 1967. D. Sc. (Eng), professor of "Mathematical Modeling" Department of the Bauman Moscow State Technical University. Academician of the Russian Academy of Cosmonautics named after K.E. Tsiolkovsky, winner of N.Ye. Zhukovsky prize of Russian Academy of Sciences. Author of more then 70 publications in the field of applied mathematics.

Ольга Олеговна Мякинник окончила в 1984 г. Львовский государственный университет. В 1997-2000 гг. проходила стажировку на кафедре "Математическое моделирование" МГТУ им. Н.Э. Баумана. Соискатель кафедры "Математическое моделирование" МГТУ им. Н.Э. Баумана. Автор 15 научных работ в области нелинейных уравнений математической физики.

O.O. Miakinnik graduated in 1984 from Lvov State University. Trainee of "Mathematical Simulation" Department of the Bauman Moscow State Technical University (1997-2000). Applicant for academic degree of "Mathematical Simulation" Department of the Bauman Moscow State Technical University. Author of 15 publications in the field of non-linear equations of mathematical physics.