CC BY
12
МАТЕМАТИКА MATHEMATICS
УДК 517.55
DOI 10.18413/2075-4639-2018-50-3-245-253
О РАЗРЕШИМОСТИ И УСТОЙЧИВОСТИ ЗАДАЧИ КОШИ ДЛЯ ПОЛИНОМИАЛЬНОГО РАЗНОСТНОГО ОПЕРАТОРА С ПЕРЕМЕННЫМИ КОЭФФИЦИЕНТАМИ
ON THE SOLVABILITY AND STABILITY OF THE CAUCHY PROBLEM FOR A POLYNOMIAL DIFFERENCE OPERATOR WITH VARIABLE COEFFICIENTS
М.С. Апанович M.S. Apanovich
Красноярский государственный медицинский университет
имени профессора В.Ф. Войно-Ясенецкого, Россия, 660022, г. Красноярск, ул. Партизана Железняка, 1
Krasnoyarsk State Medical University named after Prof. V. F. Voino-Yasenetsky, 1 Partizana Zheleznyaka St., Krasnoyarsk, 660022, Russian,
Аннотация
Рассматривается проблема разрешимости и устойчивости полиномиальных разностных операторов, основным источником появления которых является теория разностных схем. Основной результат работы сформулирован в теореме, в которой приведены условия на переменные коэффициенты полиномиального разностного оператора, обеспечивающие его разрешимость и устойчивость. Общая идея доказательства состоит в том, что мы сводим задачу Коши для полиномиального разностного уравнения к вопросу о разрешимости бесконечной системы линейных уравнений с бесконечным числом неизвестных.
Abstract
We study the solvability and stability of polynomial difference operators coming mainly from the theory of difference schemes. The main result of the paper is formulated in Theorem, which contains conditions on the variable coefficients of a polynomial difference operator that ensure its solvability and stability. The general idea of the proof is that we reduce the Cauchy problem for a polynomial difference equation to the question of the solvability of an infinite system of linear equations with an infinite number of unknowns.
Ключевые слова: полиномиальный разностный оператор, задача Коши, разрешимость, устойчивость.
Keywords: polynomial difference operator, Cauchy problem, solvability, stability.
1. Постановка задачи и основной результат
Для комплекснозначных функций /(х) целочисленных переменных х = (х1,..., Хп ) определим операторы 8 ^ сдвига по переменным х^:
83/(х) = / (х1,...,х ,_1,х, +1 х м,..., хп)
и рассмотрим полиномиальный разностный оператор порядка т
Р(б) =2са(*У ,
|а|<т
где С = (а1,...,ап) - мультииндекс, |а|=с +... + ап, да =8"1...8"п, са(х) - коэффициенты разностного оператора, представляющие собой функции целочисленных переменных
х = (Х1'...' Хп )•
Будем рассматривать разностные уравнения вида
Р(ё)/(х) = g (х), х ег:, (1)
где /(х) - неизвестная, а g(х) - заданная на г+ = г+ х... функция и г+ - множество целых неотрицательных чисел.
Для п > 1 пространство решений уравнения (1) бесконечномерно, для выделения из него единственного решения требуются дополнительные условия. Эти условия можно сформулировать различными способами (см., например, [1]). В данной работе определим их следующим образом.
Для двух точек х, у целочисленной решетки гп неравенство х > у означает, что ^ у{ для / = 1,..., п, а запись х £ у означает, что найдется ¿0 е {1,...,п\ такое, что х^ < у^ . Фиксируем мультииндекс \ такой, что
¡\ = т и ср(х) Ф 0, для всех х е Ъп+. (*)
Обозначим Х0\ = {х е г" : х £ ¡} и сформулируем задачу: найти решение /(х) уравнения (1), которое для х е Х0\ совпадает с заданной функцией р(х), т. е. удовлетворяет условию
/(х) = р(х\ х е Х05 . (2)
Если \ = (т,0,...,0) или \ = (0,...,0, т), то с точки зрения теории разностных схем мы имеем явную разностную схему (см., например, [2]). В этом случае разрешимость и единственность задачи (1) - (2) очевидны.
Для других \ таких, что \\\ = т может оказаться, что решение не единственное или
что задача не имеет решений (см. [3]).
Задачу (1) - (2) будем называть задачей Коши для полиномиального разностного оператора Р(б), а функцию р(х) - начальными данными этой задачи.
Нас интересуют условия на разностный оператор Р(р), обеспечивающие разрешимость задачи (1) - (2), т. е. существование и единственность решения для любых начальных данных ррх) и правых частей g(х), и устойчивость задачи (1) - (2), т. е. существование постоянной м такой, что при любых р(х) и g(x) для нормы соответствующего решения
/(х) задачи (1) - (2) справедлива оценка 1/(ху|ш <м(|g(хуЦ + Цр^у|ш).
Теорема. Если для коэффициентов полиномиального разностного оператора Р(б) выполнено условие (*) и
1) выполнено условие
|с\(х)| > 21 Сс(х),х ег+, (3)
\а\ = т,аФ \
то задача (1)-(2) имеет единственное решение;
2) выполнено условие
\с\(х)>2\са(х), х ег+, (4)
аф\
то задача (1)-(2) устойчива.
Замечание 1. Очевидно, что если выполнено условие (4), то выполнено и условие (3).
2. Анализ проблемы и методы исследования
Отметим, что условие, аналогичное условию (3), использовалось в [4] для доказательства разрешимости в классе аналитических функций варианта обобщенной задачи Коши для полиномиального дифференциального оператора р(б) с начально-краевыми условиями типа Рикье. Коэффициенты разложения в степенной ряд аналитических решений этой задачи удовлетворяют соотношениям вида (1)-(2). Отметим также работы [5], [6], в которых отыскание аналитических решений дифференциальных уравнений приводит к соотношениям вида (1)-(2) для коэффициентов разложения в ряд этих решений.
Случай постоянных коэффициентов рассматривается в работе [7], в которой исследована устойчивость однородной двухслойной линейной разностной схемы, а также в работах [8]—[12], где к исследованию устойчивости многослойных однородных разностных схем применяется теория амеб алгебраических гиперповерхностей и получена формула для решения задачи Коши через ее фундаментальное решение. Случай переменных коэффициентов рассматривается в [2], где исследуются однородные разностные схемы для стационарного уравнения и для уравнения теплопроводности с переменными коэффициентами.
Уравнения (1)-(2) мы рассматриваем как бесконечную систему линейных уравнений относительно бесконечного числа неизвестных /(у), у £ . После упорядочения она примет специфический вид, а именно: в каждое уравнение системы войдет только конечное число неизвестных. Такая система совместна, если всякое конечное число уравнений из этой системы совместно (см. [13], лемма 6.3.7). Построим последовательность подсистем системы (1)-(2), которые состоят из конечного числа уравнений и в каждую следующую входят все уравнения предыдущей. Совместность каждой подсистемы из этой последовательности означает, что совместно и всякое конечное число уравнений из (1)-(2).
Возьмем произвольное р е^ . Неизвестные /(у) будем "нумеровать" элементами множества Jр = {у еТ" : |у| < р| и упорядочим это множество однородно-лексикографическим способом. Уравнения "занумеруем" элементами двух множеств 1р = {х еТ" : |х| < р -т} и Iд, = | е Х0,д : || < р}. Если обозначить #м - число элементов конечного множества м , то нетрудно видеть, что # 1р +# 1р р =# Jp . Так как Iр,р " {Р" 1р }= Jр, то элементам множества Iр,р присвоим те же "номера", с которыми они входят в множество Jр, а элементам х множества Iр - те "номера", с которыми Д + х входят в Jр .
Обозначим Ьр систему уравнений относительно конечного числа неизвестных /(у), у е вида
Ес« (х)/(х + а) = 8(х) х £р , (5)
|а|< т
/ = р. (6)
Число уравнений Ы IД р) этой системы равно числу неизвестных # Jp (и равно числу Мр целых неотрицательных решений у неравенства |у| < р ). Отметим, что все уравнения из системы Ьр входят в систему уравнений Ьр+1.
Обозначим через Ад,р матрицу системы уравнений (5)-(6), а через Ад,р - ее определитель. Его порядок равен Nр и по построению в строках, соответствующих уравнениям (6), все элементы, кроме стоящей на главной диагонали единицы, равны нулю. В строках, соответствующих уравнениям (5), на главной диагонали стоит коэффициент Сд(х).
Пример 1. Рассмотрим для п = 2 разностный оператор
Р(81,82 ) = С2,0 (хр у8 " Си(х у)8182 " С0,2 (хр У)822 " С1,0 Ь у)81 " ^^ У)82 " С0,0 Ь У) ,
где
m = 2, ^ = (1,1),
с2 0(x, y) = x4 + 3 У + 3, c X(x, y) = 10x5 + 5 y + 20, c0 2(x, У) = x3 + X + 5 y3 + 4, •10(x, y) = 2x2 + x + y2 +1, c0 X(x, y) = x3 + 2x2 + x + y3 + 2, c0,0 (x, y) = 2x + 3y +1.
Для p = 3 система уравнений (5)-(6) примет вид
C2,0(x, y)/(x + 2, y)+cu(x y)f (x +1, y +0+ C0,2(x, y)f (x, y + 2)+c1,0(x, y)f (x +1, y)+
+ C»(x, y)/(x, y +1)+C00 (x, y )/(x, y ) = g(x, y), (x, y )e I3 (5')
/ (x, y) = ф, y), (x, y) e I(UX3 (6')
В ней десять неизвестных /(y, y2 ), (yx, y2 ) e J3,
J = {(0,0), (1,0), (0,1), (2,0), (1,1), (0,2), (3,0), (2,1), (1,2), (0,3)}.
Уравнения (5') нумеруются элементами множества 13= {0,0), (1,0) (0,1)}, а уравнения (6') - элементами множества I(1,1),3 = {(0,0), (1,0), (0,1), (2,0), (0,2), (3,0), (0,3)}. Так как объединение I3 U I(j j) 3 дизъюнктное, то точки с координатами (x, y) и (x, y) считаются различными. Определитель Л(1,1),3 этой системы равен
00 000 0000
Л,
(1,1 ),3
1 0 0 1
00 1 2 00
00 3 0 0 4 00
000 000
1 0 0
3 20 4 0 0 1
000
4 6 0 023 000
0 0 0 0 0 0 0 0
00 00 00
00 0 0 00
00 60
10 4 30 0 6 25 9 0 0 0 1
3. Вспомогательные утверждения
Для доказательства теоремы 1 нам потребуется несколько вспомогательных утверждений.
Лемма 1. Задача (1)-(2) для всех р(х) и g(х) имеет единственное решение тогда и только тогда, когда для всех р = 0,1,2,... определители Д\,р Ф 0.
В определителях Д\ р , р = 0,1,2,..., отвечающих за разрешимость задачи Коши, присутствуют все коэффициенты Са(х) характеристического многочлена Р(г). Покажем, что разрешимость задачи (5)-(6) зависит в действительности только от коэффициентов однородной составляющей старшей степени этого многочлена, т. е. определяется главным символом разностного оператора.
Обозначим I\ = {х е г+ : |х| = ч - т} , I= \ е : \\ = Ч} и 3\ = {у е г+ : |у| = ч} . Нетрудно видеть, что #I' ч +#I\ ч =#Т ' ч = Nч, где Nч - число целых неотрицательного решений уравнения у1 +... + уп = Ч. Заметим, что N0 +... + Nр = Ыр - числу целых неотрицательных решений неравенства у] +... + уп < р.
Для ч = 0,1,..., р обозначим Dpq миноры определителя Д\,р, составленные из его строк, соответствующих уравнениям
Р(5)/(х) = g(х), х е Г ,
/(и) = р(\), \е Г
и столбцов, соответствующих неизвестным /(у), где у е Тч .
Отметим, что в определителях присутствуют только коэффициенты однородной составляющей старшей степени многочлена Р(х).
Пример 2. Для разностного оператора Р(рх ,£2), рассмотренного в примере 1, имеем для ч = 0,1,2,3
ц(1,1),0 = 1, D(],]),]'-
1 0 0 1
Ц
1 0 0 3 20 4 0 0 1
Ц
(1Д).3
1 0 0 0
4 30 6 0
0 6 25 9
0 0 0 1
Замечание 2. Если в определителе ч вычеркнуть строки и столбцы, соответствующие уравнениям (6), то получим определитель ленточной матрицы (см, например, [14]), которая использовалась в [3] для доказательства разрешимости задачи (1)-(2) в случае п = 2.
Связь между определителями Д\р и ч дается следующей леммой.
Лемма 2. Для всякого р е г+ справедливо равенство Д\,р = •
Пример 3. Для оператора из примеров 1, 2 в соответствии с леммой 2
имеем
Д
(1,1), 0
1 0 0 1
00 00
0 0
00 00
00 00
0 0 1 0 0 0 0 0 0
1 2 3 20 4 0 0 0 0
1 , А1(1,1),3 = 0 0 0 0 1 0 0 0 0
0 0 0 0 0 1 0 0 0
3 0 4 6 0 4 30 6 0
0 4 0 2 3 0 6 25 9
0 0 0 0 0 0 0 0 0
Д(1,1),1=
1 0 0 1
А
(1,1),3
1 0 0 0 0 0 0
3 20 4 0 0 0 0
0 0 1 0 0 0 0
0 0 0 1 0 0 0
4 6 0 4 30 6 0
0 2 3 0 6 25 9
0 0 0 0 0 0 1
Д(1,1),2 =
1 0 0 3 20 4 0 0 1
1 0 0 0
А3(1,1),3 = 4 30 6 0
0 6 25 9
0 0 0 1
=Д
(1,1),3 .
А(1,1),3 = Д(1,1),0 ■ Д(1,1),1 ■ Д(1,1),2 ■ Д(1,1),3 .
Из лемм 1 и 2 сразу следует
Лемма 3. Задача (5)-(6) имеет единственное решение тогда и только тогда, когда Дрк Ф 0 для всех к = 0,1,2,....
4. Доказательства
Доказательство леммы 1.
Необходимость. Предположим, что для некоторого р0 е определитель Ар,^ = 0. Это означает, что найдется нетривиальная линейная комбинация левых частей уравнений системы Ьр равная нулю, т. е.
£ Ах £ са(х)/(х + а)+ £ В/М = 0.
хе1 \а\<т ^еI
Р, рО
Но тогда и для любых правых частей g (х) и начальных данных р(х) этой системы получим
2 Ых)+ 2 врМ=0,
хе:!р0 ^¡^0
а в силу произвольности выбора g(х)и р(х)все коэффициенты Ах и В{ равны нулю.
Достаточность. Так как для любого р е г+ определители Д\,р Ф 0, то любая подсистема системы уравнений (1)-(2) из конечного числа уравнений совместна. Используя упомянутую выше лемму из [13] (лемма 6.3.7) получим, что решение задачи (1)-(2) существует. Если предположить, что оно не единственное, то для некоторого р окажется, что
система уравнений Ьр с правой частью равной нулю имеет нетривиальное решение, т. е. Д\ = 0. Противоречие.
Доказательство леммы 2.
Будем последовательно преобразовывать определитель Д\,р , пользуясь теоремой Лапласа.
На первом шаге (ч = 0) разложим определитель Д\,р по строкам с номерами х е I'0, \е I\ 0. Единственный отличный от нуля минор порядка N'0 = 1 обозначим 0, очевидно = 1.
Определитель, полученный из Д\ р вычеркиванием этой строки и столбца, соответствующего неизвестной /(0), обозначим Д\,р . Очевидно, что Д\, = Ц\0 • Д\.
На втором шаге (ч = 1) разложим определитель Д\, р по строкам с номерами из множеств I \ и I'¡,1. Среди миноров порядка N\ не содержит нулевых строк только определитель, составленный из столбцов с номерами у е Т\. Обозначим этот минор 1, тогда Д\р = Цц • Д\,р , где определитель Д\,р получен из определителя Д\, р вычеркиванием строк с номерами х е I\ и \ е I'¡д и столбцов с номерами у е Т\.
Продолжая процедуру, на (ч +1)-ом шаге среди миноров порядка Nу определителя Д\ р , соответствующих строкам с номерами из множеств I и I\ ч не содержит нулевой строки только определитель, составленный из столбцов с номерами у е Т ' ч , т. е. цр,ч . Таким образом, Д\,р = • Д\+р,} .
Так как на последнем р -ом шаге Д\ = Ц\,р, то Д\,р = • Цц •...• Ц\,р.
Доказательство теоремы 1.
1. У определителя Д\ р на главной диагонали стоят единицы и выделенный коэффициент с\(х). Согласно лемме 2 имеем Д\,р = • Цц •... • , где - главные миноры определителя Д\,р порядка N' ч . Определители Ц\,ч зависят только от коэффициентов са(х) однородной составляющей старшей степени. Если выполнено условие (3), то являются определителями матриц с диагональным преобладанием, поэтому ч Ф 0 . По лемме 3 задача (1)-(2) имеет единственное решение.
2. Обозначим ЦАЦ^ = max £ а максимум-норму матрицы a , а через R(A) = |ай| -
' j=i
|aj , i = 1,2,..., n , обозначим величину диагонального преобладания в i -ой строке, а
j *i
также положим R (a) = min R (a). Если R* (A)> 0 , то A - матрица с диагональным преоб-
1<i<n
ладанием (см., например, [15],[16]).
Для доказательства устойчивости потребуется оценка нормы матрицы Aßß, обратной к матрице Aß подсистемы (5)-(6).
Для любого p е Z+ величины диагонального преобладания в строках с "номерами" /е Ißp равны Rü(Aß,p)= 1, а в строках с "номерами" x е Ip равны R(a^ р) = С(x)-- £ С (x) и эти величины не зависят от р. В силу условия (4)
а у
a^ß
R Wp H
1, |ca(xH = R* ф 0. Для матрицы APpp с диагональным преоблада-
1 { а*р ) ,
нием справедлива оценка (см., например, [15],[16])
KIIи- R. <7>
Рассмотрим вектор fp с координатами f (y), y £ J и вектор hp с координатами g (x ), ç(m) , где x £ Ip, U £ Ip,p , KL - (l g(x)L + Mx)| J для всех p . Систему (5)-(6) запишем в виде Appfp = hp, где A^,p Ф 0, отсюда найдем fp = A-php. Оценим норму fp с учетом (7):
iwl =IaIaL -KLKL - RpX, - «(|g(x)L+lKx)iJ.
Поскольку неравенство справедливо для любого p, то
llfll U- M (l g(x)L+lk(x)L).
Таким образом, задача (1)-(2) устойчива.
Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта №18-31-00232.
Список литературы References
1. Bousquet-Mélou M., Petkovsek M. 2000. Linear recurrences with constant coefficients: the multivariate case. Discrete Mathematics, 225: 51-75.
2. Самарский А.А. 1977. Теория разностных схем. М., Наука, 656.
Samarskiy A.A. 1977. Teoriya raznostnykh skhem [Theory of difference schemes]. Moscow, Nauka, 656. (in Russian)
3. Рогозина М.С. 2014. О разрешимости задачи Коши для полиномиального разностного оператора. Вестн. НГУ. Сер. матем., мех., информ., 14(3): 83-94.
Rogozina M.S. 2016. Solvability of the Cauchy Problem with a Polynomial Difference Operator. Journal of Mathematical Sciences, 213 (6): 887-896.
4. Хермандер Л. 1965. Линейные дифференциальные операторы с частными производными. Перев. с англ. М., Мир, 379. (Hormander L. 1963. Linear partial differential operators. Springer, 292.)
Khermander L. 1965. Lineynye differentsial'nye operatory s chastnymi proizvodnymi [Linear partial differential operators]. Moscow, Mir, 379. (Hormander L. 1963. Linear partial differential operators. Springer, 292.)
5. Гюнтер Н.М. 1925. О распространении теоремы Коши на любую систему уравнений в частных производных. Математический сборник, 32 (2): 367-447.
Gyunter N.M. 1925. O rasprostranenii teoremy Koshi na lyubuyu sistemu uravneniy v chastnykh proizvodnykh [On the extension of Cauchy's theorem to any system of partial differential equations]. Sbornik: Mathematics, 32 (2): 367-447. (in Russian)
6. Казаков А.Л. 2007. Обобщенная задача Коши с данными на двух поверхностях для квазилинейной аналитической системы. Сибирский математический журнал, 48 (5): 1041-1055.
Kazakov A.L. 2007. The generalized Cauchy problem with data on two surfaces for a quasilinear analytic system. Siberian Mathematical Journal, 48(5): 837-848.
7. Федорюк М.В. 1987. Асимптотика: интегралы и ряды. М., Наука, 546.
Fedoryuk M.V. 1987. Asimptotika: integraly i ryady [Asymptotics. Integrals and series]. Moscow, Nauka, 546. (in Russian)
8. Лейнартас Е.К. 2011. Устойчивость задачи Коши для многомерного разностного оператора и амеба характеристического множества. Сибирский математический журнал, 52 (5): 10871095.
Leinartas E.K. 2011. Stability of the Cauchy problem for a multidimensional difference operator and the amoeba of the characteristic set. Siberian Mathematical Journal, 52 (5): 864-870.
9. Рогозина М.С. 2012. Устойчивость многослойных разностных схем и амебы алгебраических гиперповерхностей. Журнал СФУ. Математика и физика, 5 (2): 256-263.
Rogozina M.S. 2012. Stability of multilayer finite difference schemes and amoebas of algebraic hypersurfaces. Journal of Siberian Federal University, 5 (2): 256-263. (in Russian, with English summary)
10. Некрасова Т.И. 2014. Об иерархии производящих функций решений многомерных разностных уравнений. Известия Иркутского государственного университета, 9: 91-102.
Nekrasova T.I. 2014. On the Hierarchy of Generating Functions for Solutions of Multidimensional Difference Equations. The Bulletin of Irkutsk State University, 9: 91-102. (in Russian, with English summary)
11. Лейнартас Е.К., Рогозина М.С. 2015. Разрешимость задачи Коши для полиномиального разностного оператора и мономиальные базисы факторов в кольце полиномов. Сибирский математический журнал, 56 (1): 111-121.
Leinartas E.K., Rogozina M.S. 2015. Solvability of the Cauchy problem for a polynomial difference operator and monomial bases for the quotients of a polynomial ring. Siberian Mathematical Journal, 56 (1): 92-100.
12. Лейнартас Е.К., Ляпин А.П. 2009. О рациональности многомерных возвратных степенных рядов. Журнал СФУ. Математика и физика, 2 (4) 449-455.
Leinartas E.K., Lyapin A.P. 2009. On rationality multidimentional recursive power series. Journal of Siberian Federal University, 2 (4): 449-455. (in Russian, with English summary)
13. Хермандер Л. 1968. Введение в теорию функций нескольких комплексных переменных. Перев. с англ. М., Мир, 280. (Hormander L. 1966. An introduction to complex analysis in several variables. Van Nostrand, 208.)
Khermander L. 1968. Vvedenie v teoriyu funktsiy neskol'kikh kompleksnykh peremennykh [An introduction to complex analysis in several variables]. Moscow, Mir, 280. (Hormander L. 1966. An introduction to complex analysis in several variables. Van Nostrand, 208.)
14. Ильин В.П., Лиснянский И.М. 1978. О решении алгебраических уравнений с ленточными теплицевыми матрицами. Сибирский математический журнал, 19 (1): 44-48.
Il'in V.P., Lisnyanskiy I.M. 1978. O reshenii algebraicheskikh uravneniy s lentochnymi tep-litsevymi matritsami [On the solution of algebraic equations with ribbon Toeplitz matrices]. Siberian Mathematical Journal, 19 (1): 44-48. (in Russian)
15. Taussky O. 1949. A recurring theorem on determinants. The American Mathematical Monthly, 56 (1): 672-676.
16. Ahlberg J.H., Nilson E.N. 1963. Convergence properties of the spline fit. J.SIAM, 11 (1): 95104.
