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

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

М. В. Коровина, В.Ю. Смирнов

Аннотация — Целью статьи является исследование асимптотик решений вырожденных эллиптических уравнений. Технический аппарат, который применяется для интерпретации полученных асимптотик, называется ресургентным анализом. Его основанием является интегральное преобразование Лапласа-Бореля. Обратное преобразование Лапласа-Бореля дает метод регулярного суммирования рядов. В этой статье мы решаем проблему построения асимптотики обратного преобразования Лапласа-Бореля одной функции экспоненциального роста. Это преобразование необходимо для построения асимптотик решений дифференциальных уравнений с вырождением в коэффициентах методом повторного квантования. Этот метод используется для изучения асимптотик решений дифференциальных уравнений с голоморфными коэффициентами. В статье приводится пример использования этого метода для одного уравнения 4-го порядка. С использованием метода повторного квантования и обратного преобразования функции экспоненциального типа, полученного в этой статье, строятся асимптотики решения этого уравнения.

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

I. Введение

Работа является продолжением исследований, касающихся построений асимптотик решений линейных дифференциальных уравнений с голоморфными коэффициентами для случая, когда коэффициенты при старших производных вырождаются. Для обыкновенных дифференциальных уравнений эта классическая задача хорошо изучена в случае, когда точка вырождения является регулярной особой точкой, например, в работах Кондратьева [1], [2]. В этих работах показано, что в случае регулярной особой точки асимптотики решений являются конормальными. Однако, в том случае, когда особая точка является

Статья получена 5 июля 2018.

М. B. Коровина, Московский государственный университет им. М. В. Ломоносова, факультет Вычислительной математики и кибернетики, Москва, РФ (e-mail:[email protected]).

В. Ю. Смирнов, Московский авиационный институт (национальный исследовательский университет), Российский университет транспорта (МИИТ), Москва, РФ (e-mail: [email protected]).

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

Обыкновенные дифференциальные уравнения с голоморфными коэффициентами могут быть представлены в виде

где

HI г,-1 rk+l d)u = f, 1 k dr )

h (r, p) = £ a (r ) p,

(i)

(2)

здесь через а1 (г) обозначены голоморфные функции,

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

Если в представлении уравнений (1) минимальное к равно нулю, то такое вырождение называется коническим, это случай регулярной особой точки. Асимптотика решения уравнения (1) в некоторой окрестности нуля является конормальной, а именно имеет вид

т ю

^га' \П г^ а^г' .

j=1 '=0

ю

Здесь а' - некоторые числа, ^ а'г1 - сходящийся

1=0

ряд. В случае, когда минимальное к > 0, вырождение называется вырождением типа клюва порядка к+1. Это случай иррегулярной особенности. Рассмотрению частных случаях этой задачи посвящены главы в таких классических книгах, как [4], [5].

К примеру, в работе [4] строятся асимптотики уравнения вида

— | и(х) + а1 (х}-—\и(х) + а0(х)и(х) = 0 (3)

—х ) ^ —х )

в окрестности бесконечности, считая, что коэффициенты а1 (х) голоморфны в некоторой ее

i=0

окрестности. Эта задача, путем замены г =— сводится к

х

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

В книге [4] так же, как и в работах [5], [6] построены асимптотические разложения решений некоторых обыкновенных дифференциальных уравнений, они были получены в виде произведений соответствующих экспонент на расходящиеся ряды, однако вопрос об интерпретации полученных расходящихся рядов был оставлен открытым, иными словами метод суммирования этих расходящихся рядов отсутствует. Такие асимптотики в дальнейшем получили названия ВКБ асимптотик. Это название появилось при решении некоторых задач квантовой механики, где были получены асимптотические разложения такого типа.

В конце 80-х годов прошлого века был получен аппарат пригодный для суммирования подобных рядов, основанный на преобразовании Лапласа-Бореля и понятии ресургентной функции, впервые введенном французским математиком Ж. Экалем [7].

В работах [8], [9] была доказана теорема о ресургентности решений обыкновенных

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

Благодаря этому результату в работах [8], [10] были построены регулярные асимптотики решений для случая, когда корни старшего символа H о (p )= H (0, p ) имеют первый порядок. В этом случае асимптотики решений, имеют вид

akrk,

(4)

ak -

где а, i = 1,...,п - корни полинома Н0(р) и а] и некоторые комплексные числа.

Однако методы, которые использовались для построения асимптотик решений в случае, когда основной символ имел простые корни, оказались неприменимыми в случае кратных корней, кроме случая уравнений второго порядка. В работах [11], [12] вопрос о построении асимптотик решений для уравнения второго порядка с произвольными голоморфными коэффициентами решен. Также в этой работе решена задача о построении асимптотик решений для уравнения Лапласа на многообразии с особенностью типа клюва.

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

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

Иными словами, суть метода повторного квантования состоит в том, что преобразование Лапласа-Бореля делается два раза. То есть, сдвинув корень р. с помощью замены

PJ_

u(г) = e г UJ(г)

в точку ноль, делается преобразование Лапласа-Бореля и получается интегро-дифференциальное уравнение относительно функции и. (p), которая будет иметь

особенность в нуле. Для того чтобы найти асимптотику этой функции в нуле делается преобразование Лапласа-Бореля второй раз. При этом мы получаем интегральное уравнение относительно функции и. ^), которая

является преобразованием Лапласа-Бореля функции и. (р). Для полученного уравнения, с помощью метода

последовательных приближений строится асимптотика его решения в окрестности корней главного символа. Найдя асимптотику функции и. ^), и затем, сделав

обратное преобразование Лапласа-Бореля, строится асимптотика функции и. (р) в окрестности нуля. Эта

асимптотика может быть ВКБ-асимптотикой или иметь более общий вид и содержать экспоненты с дробными степенями, поэтому возникает необходимость вычисления обратного преобразования Лапласа-Бореля

а

от функций вида ер g(р), где п,к е N, а е Я. Проблеме построения асимптотики этого преобразования и посвящена данная статья.

II. Определение преобразования Лапласа-Бореля и его СВОЙСТВА

Дадим определения некоторых понятий ресургентного анализа, которые понадобятся нам ниже.

Обозначим

через

S„

сектор

^£={"1 -s< argr <£H <R}. Будем говорить, что

функция f аналитическая на SR £ имеет не более чем к-

экспоненциальный рост, если существуют такие неотрицательные константы С и а , что в секторе SR £

выполнено неравенство

< Ce

Через Ek(Srs) обозначим пространство функций

голоморфных в области

S

которые k-

экспоненциально растут в нуле.

Основным техническим средством для построения теории ресургентных функций является преобразование Лапласа-Бореля. В монографии [14] описано определение и свойства этого преобразования в классе гиперфункций к-экспоненциального роста при к=1.

u

i=l

k=0

a

Г

Обозначим через E(рЛе) пространство голоморфных функций экспоненциального роста в области

= \ p

-П-е< arg p <П + е\p| > R J-, а через e(c )

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

k -преобразованием Лапласа-Бореля функции f (г) е Ek называется отображение

Bk: Ek Ц Е( О^ )/Е(С), такое что

'0

Bkf=J

= i p' rf (r) 7+Г •

о '

dr

Обратное k-преобразование Лапласа-Бореля определяется формулой

k

Bk-lf = J ep 'rk f ( p) dp •

2m J

Контур у изображен на рисунке 1 в работе [8]. Заметим, что для ^-преобразования Лапласа-Бореля верны формулы

Bk ° (-1 rk+1 |)f (r) = pBf

dp ° Bkf=-Bk ((r)

III. Основной результат В этом параграфе мы найдем обратное преобразование

Лапласа- Бореля от функции e р" g (p ).

а

к

ТЕОРЕМА. Асимптотика функции В~1е р g(р) в окрестности нуля имеет вид

Ё exP

ZO-

n+k \ r-j f n+k

Cr"

(5)

v r„™

dr

2 2 + С r I - r

u + cn-irI - r

k+n-2

2

dr

u +

dr

n-U 2 +... + cr I - r

dr

u +

. n . n 2^1 n

+ cr u + cnr I - r — I u = 0, dr

здесь через c¡, 1 = 0,..., п -1 и е обозначены некоторые константы . Сделаем преобразование Бореля

рп+к~ + С0 (- 1)" р42 р^ )—Р!...—Ря +

1 1

+ С1 (- 1)И-1 |...| Р1к+'и(Р1 )—Р1...—Ри-1 + 11

2 Рр2

с2 V- Ч 11 Р1 "\Р1 —Р—2

+ С2 (- i)n-2 J J pik+2f (pi )dpidp2 + ... + 1 1

+ Cn-i (- i)p pi+n-if (pi )dpi + i

p p2

+ c(- l)n J... J f (pi )dpi...dpn = f (p).

i i

Через f (р) обозначена произвольная голоморфная функция. Продифференцируем п раз

' — "

dp

n+k . / л \n k . {л V-1

p u + c0 (- i pu + сД- i

+ C2 (- i)n

dp

fd-T pk+2u+...+ dp )

k+i . | p u +

(7)

( d Y-i

+ Cn-i(- \dp} pk+n-iu +(- i)ncu = f(p) ■

Покажем, что можно выбрать константы

Здесь

OU j = k + «

корни полинома

е,е1,1 = 0,...,п-1 так, чтобы функция и(р) = рав р являлась решением уравнения (7), для этого подставим эту функцию в уравнение.

Тогда получим уравнение относительно констант е, е.,' = 0,...,п -1

n+k . p +

n+k

n + k !

k

(- i)n+i I ak I , а коэффициенты aj при

i < n - i, С, - некоторые константы.

ДОКАЗАТЕЛЬСТВО.

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

k+n-i

d

d

2

r

+

d

k

а

r

n

( , k , 2k , 3k ^

, a+k— a+k----a+k--

n a+k , n „ , n n . — n —

a0 p + a1 p n +a2 p n +a3 p n +... + anp

\

+ c_

/

k 2k 3k

, , , a+k— , a+k----, a+k--

n-1 a+k . n-1 n . n-1 n . n-1 n .

a0 p + a1 p n +a2 p n +a3 p n +

v

, (n-l)k Л

, a+k---'—

+...+a::; p n

+ c„

/

f , k , (n-2)k Л

1 a+k — a+k-!—

a>1 pa+k +ain-2 p n +... + ann-:22 p n

(

+... + c1

v

, kA

1 ir 1 a+k—

al p +a1 p n

j

+ C0 pa + k + C(: i)npa= 0,

здесь

a/, i = 1,...,-, j = 1,...,n -

соответствующие

ak

константы. Заметим, что a- = I — | . Из последнего

I -

уравнения следует

c(- l)n =-an"

a"-1c , = -a"

an-2c +an-2c =-an

"n-1 °n-2 T Lin-2L'n-1 n-2

n—1 n—2 1 n

a0 Cn-1 +a0 Cn-2 + ... + a0C1 + C0 = -a0 .

1 k+1 d --X -

k dx

n + k k f 1 k +1 d

+ n + k | x2k f 1 x^+1 d " 1 k J V k dx

n-1

u + cn-1-x" I--x"

k V k dx

k+n-2

u+

u+

+ ... + C,

n+k k

n+k

(n-1)k ( - I k+1

1 k+1 d

k+1

xv" *'"I--x""-| U +

k dx

n + k

+1-| CU + C,

k J 01 k

-+1+i

n + k л" nk f 1 k+1 d N*

"f n + k

+sl —

x)kI-1xk+1 d | u = 0. k dx J

k dx

k-1-i

u+

Так как основной символ этого уравнения равен

н 0 (р)=рп+к+(- 1)п+1 Г Т+к с=рп+к +

(-1)-

7 \ n+k / 7

-+k Л f ak

k

и все корни полинома Н0 (р) имеют первый порядок, то асимптотика решения уравнения (6) в нуле имеет вид

n+k f k

* i^w la

u(x)=S exp| Kj s

j=1 V i=1

x

x j S CJ,x'

Очевидно, что эта система имеет единственное решение, отсюда следует, что можно выбрать константы

а

к

с,с,i = 0,...,п -1 так, чтобы функция м(р) = рае р" была решением уравнения (7). Решение и( х) уравнения (6) является обратным преобразованием Лапласа-Бореля этой функции.

Вернемся к уравнению (6) и найдем его решение.

Разделим это уравнение на гп и запишем его в виде

где aj, j = 1,..., - + k являются корнями полинома H 0 (p ). Так как

n+k s a

V xj xa S e x s Ci j=1 ,

(

S Cix = S exp

j=1

kj

л

s

^S c

V r J

то окончательно получили, что асимптотика функции

a

k

-1цa„ p"

7~>—1 a d"

B1 p e v в окрестности нуля имеет вид

^ 2 - d _ Г n+k _

dr

k f 2—— d _ Г n+k _

dr

U +

k n+k

B-1 paep «Sexp

j=1

S-

V r"™ J

r -+k S cj r -+k l

k

- r

7

2--+k d

dr

+... +

(n-1)

- r

V

nk

'' 7

2--+k d

dr

u +

- r

1 (n+i)k f 2 -

'' 7

2-n+k d

k

+ Sair

i=1

-r

dr

_d_ dr

u+

/

k-1-i

u = 0,

Теорема доказана.

IV. Пример

В качестве примера применения формулы (5) построим асимптотику решения уравнения 4-го порядка с кратными корнями, а именно уравнения с индексом сингулярности равным трем.

Пусть индекс сингулярности равен трем, тогда рассмотрим уравнение

2 d Л4 , f 2 d Л2 2f 2 d Л

- r — | u + c1r| - r — | u + c2r | - r — |u + dr J V dr J V dr J

здесь через а1, с', i = 1,..., к -1,. = 0,..., п -1 обозначены соответствующие константы. Сделаем в (8) подстановку

x = rn+k, тогда получим уравнение

3 21 2 d

+ c3r u + c4r | - r — | u +

3 4 dr J

+ c5r| - r2 u + r3 H f r J - r2 d | |u = 0.

(8)

+

k+--1

+

+

-

-

и

i=1

k+n—1

u + Cn-1r

i=1

k+n-2

2

+ Cn- 2r

k

+ c.r

+ cu + c0 r

Здесь через е1,' = 1,...,5 обозначены некоторые константы и выполнено условие е1 Ф 0. Через

2 d

+ a \- r2 + r Z b (rí-r2 d dr J i=o i dr

(9)

+ e2 (- 1)2 j j Pi~(Pl )dPl dP2 + 11

+ e4 (- 1)2 j j Pl2y(Pl )dPldP2 +

(10)

y

С5 (- !)j P3«(P)dP-

метод последовательных приближений аналогично тому, как это сделано, к примеру, в работе [8], легко

H \ r, r — II обозначен дифференциальный оператор с гладкими коэффициентами.

r ' d)]" ai i- r ' d ) + a, (-r" dr I +

показать, что слагаемое

Bj BH \ r, \ r2 d|I u(r )dp

dr

дает

младшие члены асимптотики функции и(д) в окрестностях точек д = 0 и д = е1. Это следует из равенства

B

где а1,1 = 1,2,3 некоторые константы, а через Ь1 (г) -соответствующие голоморфные функции. Сделаем преобразование Лапласа-Бореля уравнения (8)

р 4и(р) + е1 (-1)| р 2и(р )—р +

+ (-1)3 р2 ВН(г, (г2 — )) и (г)—р,—р2—рз = f (р).

Через f (р) обозначена произвольная голоморфная функция. Продифференцируем это уравнения два раза:

(- р2 —) ~(р) + С' (- р2 —) ~(р) +

+ е2 рм(р) + е1 р 2м (р ) +

+е4 р(- р2 — ) р)-

(11)

—р

~\ВН ( г,( г' — )) и(г *= (— )2 /(р>

Здесь через е1,' - 2,3,4 обозначены соответствующие константы. Применим метод повторного квантования. Основной символ имеет два корня первого порядка д = 0 и д = е1. Легко показать, что слагаемое

|ВН (г, (г2 —)\ и(г )—р дает младшие члены

асимптотики функции и(р) в окрестности нуля. Это

следует из того, что, сделав преобразование Лапласа-

Бореля уравнения (11), получим уравнение:

д

д(д - е1) и (д) + |(е2 + е4д)и(д)—д +

1

+ е\ Ц и(д2) —д2 —д1 + (12)

11

+ В |ВН ('', (' ■2 — \\и(г — = В( —\2 f(ру.

Выразив из уравнения (12) функцию и(д) и проведя

р 1 д

|м(р)—р =—Цг5(д1)—д1—д , (13)

) д <?0 которое было доказано в работе [13].

Таким образом, применения метод последовательных приближений к уравнению (12) с учетом того, что слагаемое (13) дает младшие члены асимптотики, а корни главного символа являются простыми, получим асимптотику решения уравнения (11) в окрестности точки р - 0. Она будет иметь вид

ю ю

и(р)* е рр^-Xа,1 р' + р^2Xа'р .

1=0 1=0 Чтобы получить асимптотику решения уравнения (8) сделаем обратное преобразование Лапласа-Бореля. Из формулы (5) получим

B-1

p"1 za) p'+p"2 z aiPi

= e^r r^ Zejr2 + "2 ze/r2 + r'3 Ze/r.

l=0 l=0 l=0

Здесь a;, i = 1,2 - корни многочлена p2 + 4e,.

БИБЛИОГРАФИЯ

[1] Кондратьев В. А. Краевые задачи для эллиптических уравнений в областях с коническими или угловыми точками // Труды Московского математического общества. - 1967. - Т. 16. - С. 209-292.

[2] Кондратьев В. А. Краевые задачи для эллиптических уравнений в конических областях. Докл. АН СССР, 153:1 (1963), 27-29

[3] Кац Д. С. Вычисление асимптотик решений уравнений с полиномиальными вырождениями коэффициентов // Дифференциальные уравнения. — 2015. — Т. 51, № 12. — С. 1612-1617.

[4] Олвер Ф. Асимптотика и специальные функции. Пер. с англ. под ред. А. П. Прудникова. - М.: Наука. Гл. ред. физ. -мат. лит,1990. - 528 с.

[5] Коддингтон Э.А., Левипсоп Н. Теория обыкновенных дифференциальных уравнений. М.: Иностранная литература, 1958. — 475 с.

[6] Л. Чезари. Асимптотическое поведение и устойчивость решений обыкновенных дифференциальных уравнений. М.: Мир, 1964. — 477 с.

[7] J. Ecalle. Cinq applications des fonctions résurgentes. // Prepub. Math. d'Orsay, 1984, 84T62, # 110 pp.

[8] Коровина М. В., Шаталов В. Е. Дифференциальные уравнения с вырождением и ресургентный анализ // Дифференциальные уравнения. — 2010. — Т. 46, № 9. — С. 1259-1277.

[9] Коровина М. В. Существование ресургентного решения для уравнений с вырождением высших порядков // Дифференциальные уравнения. — 2011. — Т. 47, № 3. — С. 349357.

[10] Коровина М. В. Асимптотики решений уравнений с высшими вырождениями // Доклады Академии паук. — 2011. — Т. 437, № 3. — С. 302-304.

[11] Коровина М. В. Асимптотики решений уравнений второго порядка со старшими вырождениями и уравнение Лапласа на многообразии с особенностью типа клюва // Доклады Академии паук. — 2014. — Т. 456, № 4. — С. 396-399.

e

[12] Korovina M. V. Asymptotic solutions of second order equations with holomorphic coefficients with degeneracies and laplace's equations on a manifold with a cuspidal singularity // Global Journal of Science Frontier Research (GJSFR): F Mathematics & Decision Sciences. — 2017. — Vol. 17, no. 6. — P. 57-71.

[13] Коровина М. В. Метод повторного квантования и его применения к построению асимптотик решений уравнений с вырождениями // Дифференциальные уравнения. — 2016. — Т. 52, № 1. — С. 60-77.

About one problem that occurs when applying the repeated quantization method to linear differential equation with holomorphic

coefficients

M. V. Korovina, V.Y. Smirnov

Abstract — The present paper deals with asymptotic expansions for solution degenerate elliptic differential equations. The technique for the interpretation and construction of asymptotic expansions on the basis of the Laplace-Borel transform is referred to as resurgent analysis. The inverse Laplace-Borel transform provides a regular method for the summation of the series. In this paper, we solve the problem of constructing asymptotic expansions the inverse-transform Laplace-Borel of one exponential type functions. This transform is necessary for construction of asymptotic of solutions of differential equations with degeneration in the coefficients with repeated quantization method. This method is used to study the asymptotic of solutions of equations with holomorphic coefficients. By using the repeated quantization method, we will consider an example fourth order differential equation and used inverse-transform Laplace-Borel of one exponential type functions, which receive in the paper, we construct asymptotic expansions of solution of this equation.

Keywords — differential equations with singularities, singular points, degenerate, Laplace-Borel transformation, resurgent function, principle operator symbol, asymptotic expansion.