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

УДК 517.9

ОБ АСИМПТОТИЧЕСКОМ ПОВЕДЕНИИ РЕШЕНИЙ НЕЛИНЕЙНЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ С КОМПЛЕКСНЫМИ КОЭФФИЦИЕНТАМИ

И.В. Асташова 4)

Московский государственный университет им. М.В.Ломоносова, Ленинские горы, ГСП-1, Москва, 119991, Россия, e-mail: ast@diffiety.ac.ru

Аннотация. Для нелинейных дифференциальных уравнений типа Эмдена-Фаулера с комплексными коэффициентами получены асимптотические формулы для модуля и аргумента решений и равномерные оценки решений.

Ключевые слова: асимптотические формулы, равномерные оценки решений, комплексные коэффициенты.

1. Введение. Рассматривается дифференциальное уравнение

y"(x) = p(x)\y(x)\my(x) , (1)

где m > 0, x Е R, а p(x) - непрерывная комплекснозначная функция. Если p(x)

- действительнозначная функция, то это уравнение превращается в хорошо известное уравнение типа Эмдена-Фаулера, асимптотические свойства которого детально исследовались в работах Ф.Аткинсона, Р.Беллмана, И.Кигурадзе, А.Кнезера, В.Кондратьева,

A.Мышкиса, Дж.Сансоне и других авторов. Подробную библиографию см. в [1]. С другой стороны, (1) - это одномерное уравнение Шредингера. Качественные свойства решений различных задач, связанных с этим уравнением в n-мерном случае (п ^ 2), были описаны М.Ф.Бидо-Верон, Х.Брезисом, Л.Вероном, Б.Гершем, С.Дои, Т.Като,

B.Кондратьевым, П.Константином, Н.Хаяси, М.Шубиным и другими авторами, см. [2, 3, 4, 5, 6, 7, 8, 9, 10].

В настоящей работе получены асимптотические формулы для модуля и аргумента решений и равномерные оценки решений.

2. Основные результаты. При p(x) = p0 Е C \ R существует решение Y(x), определенное на (0, +то), которое имеет вид

\Y(x)\ = Cix-2/m, arg Y(x) = C2 lnx

4Работа выполнена при поддержке РФФИ (грант № 11-01-00989).

с постоянными

C — Q 1 + 4/m

C2 — — Q ----------------

Im p0

-Re p0 + t (Re p0)2 + 8(m,+, 2 (Im Po)2 \/ (m + 4)2

Q — V

Теорема 1. Пусть m > 0 и p(x) = p0 — const G C \ R. Тогда все нетривиальные решения уравнения (1) исчерпывающе описываются следующим образом:

1. Все непродолжаемые решения, определенные на полуоси (—то, x0) или (x0, +то), которые имеют точный вид:

\y(x)\ = I Y(|x — x0|) I ■ arg y(x) — arg Y (\x — x0\) + Ф0 с произвольными вещественными x0 и ф0.

2. Для любого непродолжаемого решения, определенного на ограниченном интервале (xi, x2), справедливо представление

\y(x)\ — |Y(|x — xk1)1 (1 + о(1)) ■

argy(x) — arg Y(\x — xk\) (1 + o(1)),

где x — xk, k =1, 2.

Теорема 2. Пусть p(x) - непрерывная комплекснозначная функция, m > 0 иp(x0) = p0 Е C\R. Пусть y(x) -непродолжаемоерешение уравнения (1), определенное на (xi,x0) или (x0,x2) при — ж ^ xi < x0 < x2 ^ +ж. Тогда

\y(x)\ = \Y(\x — xo\)\ (1 + o(1)),

argy(x) = argY(\x — x0\) (1 + o(1)),

при x —— xo.

Теорема 3. Пусть p(x) - непрерывная комплекснозначная функция, е = ±1, m> 0, p(x) — p0 Е C\R при x — еж. Пусть, далее, y(x) - решение уравнения (1), определенное в окрестности еж. Тогда

\y(x)\ = \Y(\x\)\ (1 + о(1)),

argy(x) — argY(\x\) (1 + o(1))

при х ^ ж.

Теорема 4. Пусть И,е р(х) > р* > 0. Тогда для любого решения у(х) уравнения (1), определенного на (х0 — е, х0 + е) и такого, что у(х0) = 0, справедлива оценка

е2 < С |у(Хо)|-т

Р*

с постоянной С > 0, зависящей только от т.

Следствие. Пусть для функции р(х) выполняются условия теоремы 4. Тогда для любого решения у(х) уравнения (1), определенного на [а, Ь], выполнено

ш\ <’"' C

е2р*

для всех x Е [а + е,Ь — е].

Следствие. Пусть для функции p(x) выполняются условия теоремы 4. Тогда для любого решения y(x) уравнения (1), определенного на (—ж, x0) или (x0, +ж), на всей области определения выполняется неравенство

\y(x)\ < \x — xo\-2/m VC/p* •

Следствие. Если И,е р(х) > д*х г, д* > 0, г > 0, то для любого решения у(х) уравнения (1), определенного на (0, +ж), для всех х > 0 выполнено

|у(х)| <х(г-2)/т ч/СЛ* .

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

Следствие. Если функция р(х) удовлетворяет условиям теоремы 4, то единственным решением уравнения (1), определенным на (—ж, +ж), является тривиальное решение у(х) = 0.

3. Доказательство основных результатов. Всюду ниже для уменьшения количества дробей в формулах вместо т/4 используется ц.

3.1. Фазовое пространство. Заметим, что если функция р(х) является постоянной и у(х) является решением (1), то и функция г(х) = А у (1А12^(х — х0)) при произвольных константах А Е С и х0 Е К также является решением этого уравнения. Это позволяет понизить размерность задачи, отождествляя решения, связанные приведенным соотношением.

Пара функций (у(х),у;(х)) порождает кривую в С2. Кривые, порожденные нетривиальными решениями, лежат в С2 \ {0}. Решения у(х) и у(х — х0) порождают одну и ту же кривую (с точностью до параметризации).

Рассмотрим отношение эквивалентности в С2 \ {0}, при котором решения у^х) и у2(х) = Ау1(|А|2^х) порождают одну и ту же кривую в факторпространстве. Это отношение может быть задано формулой

(■20, %1) ~ (А^0, А| А|2^^1>

для произвольного комплексного А = 0.

Обозначим через Ф факторпространство С2 \{0} по этому отношению эквивалентности. Его можно снабдить структурой действительного двумерного многообразия класса С1 с помощью атласа, состоящего из двух карт. Обе карты являются биекциями подмножеств Ф на С.

Первая карта определена на классах эквивалентности комплексных пар (г0, 21), для которых 20 = 0, то есть на всем Ф, кроме точки-классе эквивалентности пары (0,1). Биекция определяется комплекснозначной функцией

и ;[(.°.21)]^ ^ •

Вторая карта определена для классов пар (20,21), г1 = 0, следующим образом:

z0|z1|‘2^/(2^+1')

и : [(^0,^1)]

z\

Непосредственно проверяется, что эти функции корректно определены и являются биекциями. Замены координат задаются соотношениями

|и|2^/(2^+1)

и =

и

и и ^ ’ и

и принадлежат классу С1 как отображения К2 \ {0} —> К2 \ {0}.

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

Рис. 1.

3.2. Динамическая система на фазовом пространстве для постоянной р(х). Опишем в координатах кривые, порождаемые на Ф решениями (1) с постоянной р(х) = р0.

На первой карте имеем

= у'

и у|у|2м ’

откуда непосредственными вычислениями получаем

и' = р - (ш + 1)и2 - ш|и|2) .

Следовательно, выбирая в качестве параметра переменную і, для которой М = |y|2^dx, получим внутреннее описание кривой:

du , ч 2 і 12

и = — = р0 — ш + 1)и — шщ . dt

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

= 1 + }и|, ( №о\и|2 — (Ш + 1)Рои2

(^шро^и |2 — (ш + 1)р0и 2^

¿т 2л + 1

с другим параметром т, для которого ¿т = ^'|2м/(2м+1)^х.

Правые части обоих уравнений принадлежат классу С1 в действительном смысле. Из двух параметров с помощью разбиения единицы можно сделать один так, чтобы все кривые, порожденные на Ф решениями (1), были траекториями автономной динамической системы с новым параметром в качестве независимой переменной. Ввиду компактности Ф, любая траектория системы продолжена на всю ось (—ж, +ж), причем именно такие полные траектории, а не их части порождаются непродолжаемыми решениями (1).

У системы есть ровно две неподвижные точки (при условии, что р0 = 0). Они обе находятся в первой карте и отличаются только знаком. Уравнение г/ = 0, записанное в терминах V = Ие и0 и = 1т и0:

{

(2ш + 1)^0 — ^0 = Ие Р0 , 2(ш + 1)^0^0 = 1т р0 ,

можно легко решить, получив два решения: и0 = v0 + т0{, где

У0 =

^0

\

(

2ш + 1

Ие р0 + ^ / (Ие Р0)2 + (ш + 1)2 (Іт Р0)2 4ш + 2

Іт р0

2(ш + 1) У0

и второе решение —и0 = —ь0 — т0ї.

Иногда удобнее записывать систему в терминах неподвижной точки и0, а не р0:

и = (л + 1)(и0 — и2) + л(Iи012 — |и|2).

3.3. Случай и0 = ± г. Замкнутые траектории. Хотя случай и0 = ± г соответствует действительному значению р0, его исследование помогает понять поведение траекторий и для комплексных р0.

В этом случае система записывается следующим образом:

и = —1 — (л + 1)и — ли . (4)

Среди ее решений легко находится одно действительнозначное, меняющееся от +ж до — ж. На самом деле, это только часть замкнутой траектории на Ф, проходящей через единственную не покрытую первой картой точку. Так как другие траектории не могут проходить через эту же точку, они все полностью лежат в первой карте. Точнее, в полуплоскости 1т и > 0 или 1т и < 0. Ввиду инвариантности системы относительно комплексного сопряжения достаточно рассмотреть только первый случай. Для любой такой траектории, не являющейся неподвижной точкой, исследуем поведение а^(и — г), используя обозначения V = Ие и, = 1т и :

^ и

— а^(и — г) = 1т

dt и — і

(—1 — (ш + 1)и2 — ш|и|2) (и + і)

Іт

щ — і!2

—т + |и|2т + 1 + (ш + 1) Ие (и2) +

щ — і|2 —т(т2 — 1) — v2(w + 1) — 2^2 + т2 — 1

V2 + (т — 1)2

— (т2 — 1)(т — 1) — V2 (т + 1) — 2щу2

V2 + (т — 1)2

( V \

—т — 1 — 2и - г I —1 .

\|и — А)

Отсюда следует, что траектория обходит точку і по часовой стрелке, регулярно меняя знак Ие и.

Из (4) также вытекает, что все траектории симметричны относительно мнимой оси. Значит, все они, кроме действительной оси и двух неподвижных точек, представляют собой овалы, окружающие одну из неподвижных точек, причем обход г происходит по, а —г — против часовой стрелки.

Глобально Ф, покрытое траекториями, выглядит, как глобус с двумя полюсами и множеством параллелей (см. рис. 2).

3.4. Случай комплексных р0. Перейдем к случаю р0 с ненулевой мнимой частью, при этом и0 будет иметь и ненулевую действительную часть. Для использования предыдущего результата повернем и сожмем/растянем картину траекторий так, чтобы неподвижная точка, для которой Ие и0 > 0, попала в г. Это преобразование записывается в первой карте в виде и м ги/и0 и легко продолжается до глобального диффеоморфизма пространства Ф.

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

Рис. 2.

То, как ее траектории проходят через описанные выше овалы, можно выяснить, оценив знак мнимой части произведения и из (2) на и из (5).

положительно для всех и £ С, кроме ±г. Это значит, что вне неподвижных точек все траектории системы (5) последовательно покидают все овалы, лежащие в полуплоскости 1т и < 0, пересекают действительную ось, после чего последовательно проникают в овалы из полуплоскости 1т и > 0 (см. рис. 3). Эти траектории не могут иметь предельную точку, отличную от —г при Ь м —ж и отличную от г при Ь м +ж. Поэтому, ввиду компактности многообразия Ф, для всех нетривиальных траекторий эти точки являются пределами.

Таким образом, траектории системы (3) - это две неподвижные точки и траектории, стремящиеся от одной из них к другой.

Зная неподвижные точки системы, можно явно выписать семейство решений уравнения (1). Используя полярную форму для у = регф, можно записать уравнение и = и0

Рис. 3.

Это произведение равно

+ (ш2 + ш) (1 — |и|2) ((1 + и2) ио — (1 + и2) и^ .

Его мнимая часть выглядит менее громоздко

^(ш + 1)2 11 + и2|2 — ¡л2 (1 — |и2|)2^ Ие и0 .

Согласно неравенству треугольника для векторов — 1 и и2, последнее выражение строго

в виде

Решая его отдельно для действительной и мнимой частей, получим

р ^ — —2^Re Uo (x — Xo),

а затем

Im Uo і I I ,

Ф — —^------------ In Iх — Xo| + Ф0•

2^Re u0

Из неотрицательности p следует, что это решение определено на (—ж, хо). Аналогичные формулы для —и0 описывают решение, заданное на (х0, +ж).

Для остальных траекторий имеют место соотношения и ~ и0 при t ^ +ж и и ~ —u0 при t ^ — ж, приводящие к асимптотическим формулам для соответствующих решений уравнения (1), определенным на конечных интервалах Х,х2).

Учитывая (2) и возвращаясь в обозначениях к m = 4^, получим следующее описание решений.

В рассматриваемом случае p(x) = p0 = const G C \ R точное решение Y(x), определенное на (0, +ж), имеет вид

IY(х)| — C\x 2/m , arg Y(x) — C2 ln x

с постоянными

f>( !й?)2 ■

C — Q 1 + 4/m

C2 — — Q~,----------

Im p0

—Re po + 4 (Re po)2 + 8(m,+22 (Im po)2 у (m + 4)2

Q

Оказывается, что все решения уравнения (1) за исключением тривиального у = 0 имеют такую же асимптотику, как Y(х).

Таким образом, доказана Теорема 1 (см. также рис. 4).

Рис. 4.

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

й = р(£) — + 1)й2 — ц\и\2

с различными функциями р(£) для разных у(х). Однако, некоторые свойства р(х) наследуются р(£), что помогает исследовать асимптотическое поведение траекторий на Ф и решений (1).

Пусть у(х) - непродолжаемое решение уравнения (1), определенное на интервале (х1 ,х2) (возможно, неограниченном). Поскольку уравнение не изменяется при преобразовании х М —х, достаточно исследовать поведение решения у(х) только вблизи х2. Пусть р(х) М р0 при х М х2 (этот предел автоматически существует при х2 < +ж). Рассмотрим кривую, порождаемую у(х) на Ф и преобразуем Ф (как и в случае р(х) = р0) таким образом, чтобы точка й0 = У0 + /ш01, определяемая (2), перешла в г.

Преобразованная кривая не является траекторией (5). Но вне сколь угодно малых окрестностей точек ±г при р(£) близких к р0 с учетом непрерывности эта кривая пересекает овалы описанные для й0 = г в том же порядке, что и траектории (5). Таким образом, единственной возможной причиной для кривой не стремиться к ±г (или к ±йо

перед преобразованием Ф) является ограниченность t и т при х ^ х2. Докажем, что этого не может быть.

Пусть х2 < +ж. Тогда, в силу непродолжаемости, хотя бы одна из функций y(x) или у1 (х) должна быть неограниченной. Но это невозможно если t ограничено:

d

— Ь(|у' |2 + уГ+2)

< dx _ 2|у 1 \у''\ + (4^ + 2) Ы4^11у1 \

< dt \у'\2 + |у|4м+2

2|у\ (sup \p(x)| + 2^ + 1)

Ы2^ (\у'\2 + Ы4"+2)

2|у\ М2^1 (sup \p(x)| + 2^ + 1)

<

< \у\2 + \у\^+2

< sup |p(x)| + 2^ + 1.

Докажем, что случай х2 = +ж также невозможен. Пусть s G Ф - предел траектории при х ^ +ж. Рассмотрим произвольное решение уравнения (1) с p(x) = p0 и начальными условиями (у3,у3) при х3 G (х1,х2) порождающим s. Поскольку s = ±u0, существует х4 > х3 такое, что (у(х4),у'(х4)) порождает другую точку в Ф. Так как окрестность V точки s может быть выбрана так, что для любого решения (1) с p(x) достаточно близкой к p0 и начальными данными достаточно близкими к (у3,у3) соответствующая кривая в Ф будет покидать V перед х = х4. Из этого множества решений, используя подстановку z(x) = A у (|A|2^(x — х')), можно получить любое решение (1) с p( х) достаточно близкой к p0 и начальными данными (для всех х) порождающими точку в Ф, достаточно близкую к s. Следовательно, любая кривая в Ф, порождаемая решением (1) и определенная в окрестности +ж не может иметь предела, отличного от ±u0. Это доказывает Теорему 2 и Теорему 3.

3.6. Оценки. Докажем Теорему 4.

Рассмотрим вещественнозначную непрерывную функцию, связанную с решением соотношением

v (х) = 2|у(х)|'\у(х)|-2^-1 и определенную на максимальном интервале (х0 — 8*, х0 + 8*), где определено и отлично от нуля у( х) .

Воспользовавшись очевидным соотношением

2|у|' = (\у\2)'\у|-1 = (у'у + уу ')\у\-1

можно представить функцию V(х) в виде

v(х) = (у'у + уу') \у\-211-2

Дифференцируя это соотношение, получаем

V (х) =

= 2\у\2^И,е р(х) + 2\у'\2\у\_2м_2 —

— (^ + 1)У(х)2\у\2м >

> \у\2м (2р* — + 1) У(х)2 ) •

Вначале рассмотрим случай У(х0) > 0. В этом случае, поскольку У'|у 0 > 0, функция У(х) остается положительной для всех х Е (х0, х0 + 8*). Из определения У следует, что \у(х)\ возрастает на этом интервале. Следовательно, \у(х)\ > \у(х0)\ = 0 для х Е (х0, х0 + 8*) и функция У (х) определена для всех х Е (х0, х0 + е).

Обозначим у/р*/(ц + 1) через У*. Предположим, что У(х0) Е [0,У*]. Пока У(х) остается на этом отрезке, для х > х0 выполняется:

У'(х) > \у\2» (2р* — (» + 1)У*2) =

= \у\2^р* > \у(х0)\2мр* •

Следовательно, У(х1) становится равной У* для некоторых х1 Е (х0, х0 + е) таких что

< У* = \у(х0)\_2м (6)

х1 — х0 < 1—тр;-----------------------------------------------— -тгтт" • (6)

\у(х0)\2^р* (^ + 1)У*

и поэтому \у(х 1) \ > \у(х0)\.

Теперь предположим, что У(х0) > У*. Поскольку У'|у_у > \у(х0)\2мр* > 0, неравенство У (х) > У* остается справедливым х Е (х0, х0 + е). Поэтому для таких х получаем

(\у\ 2^) = —^У < —^У* и следовательно \у(х)\_2м < \у(х0)\_2м — цУ*(х — х0). Таким образом,

< \у(х0)\_2м (7)

х — х0 < ---—----- • (7)

цУ*

Оценки (6) и (7) показывают что ни одна из функций У(х) и у(х) не может быть определена при х > х0 + е для

г2 > \у(х0)\ 1 + Л 2 = 16(т + 2)2

^ У2 + 1 ц) т2(т + 4) р*

Случай У(х0) < 0 исследуется аналогичным образом, но левее х0. А именно, для отрицательной У(х0) ни одна из функций У(х) и у(х) не может быть определена при х > х0 + е для е, удовлетворяющих (8). Это завершает доказательство.

Из доказанной теоремы немедленно вытекают все четыре следствия.

Замечание. Результаты об асимптотическом поведении решений уравнения с действительными коэффициентами содержатся в [11] (см. также библиографию), об асимптотическом поведении решений уравнения (1) с постоянным комплексным коэффициентом содержатся в [12], [13], о равномерных оценках положительных решений уравнений с действительными коэффициентами - в [14], некоторые результаты об оценках решений для уравнения (1) были опубликованы в [15], см. также [16].

Литература

1. Кигурадзе И.Т., Чантурия Т.А. Асимптотические свойства решений неавтономных обыкновенных дифференциальных уравнений //М.: Наука, 1990. - 432 с.

2. Bidaut-Véron M.F. Local and global behaviour of solutions of quasilinear elliptic equations of Emden-Fowler type // Arch. Rat. Mech. Anal. - 1989. - 107. - P.293-324.

3. Brezis H., Kato T. Remarks on the Shrödinger operator with singular complex potential // J.Math. pures et appl. - 1979. - 58. - P.137-151.

4. Constantin P. Decay estimates of Schrödinger equations // Commun. Math. Phys. - 1990. -127. - P.101-108.

5. Doi S. On the Cauchy problem for Schrodinger type equations and the regularity of solutions // J. Math. Kyoto Univ. - 1994. - 34. - P.319-328.

6. Guerch B., Veron L. Local properties of stationnary solutions of some nonlinear singular Schrödinger equation // Rev. Mat. Iberoamericana. - 1991. - 7. - P.65-114.

7. Hayashi N. Global existence of small solutions to quadratic nonlinear Schrödinger equations // Comm. P.D.E.. - 1993. - 18. - P.1109-1124.

8. Kato T. Shrodinger operators with singular potentials // Israel J. Math. - 1972. - 13. -P.135-148.

9. Kato T. On some Shrodinger operators with a singular complex potential // Ann. Sc. Norm. Sup. Pisa, Ser. IV. - 1978. - 5. - P.105-114.

10. Kondrat’ev V., Shubin M. Discreteness of spectrum for the Schrodinger operators on manifolds of bounded geometry / Operator Theory: Advances and Applications, Vol. 110 / Birkhauser Verlag: Basel/Switzerland, 1999.

11. Асташова И. В. Применение динамических систем к исследованию асимптотических свойств решений нелинейных дифференциальных уравнений высоких порядков // Современная математика и ее приложения. - 2003. - 8. - C.3-33.

12. Асташова И. В. Об асимптотическом поведении решений уравнения типа Эмдена-Фаулера с комплексным коэффициентом // Современная математика и ее приложения. - 2005. -29. - C.14-18.

13. Astashova I.V. On asymptotic properties of the one-dimensional Schrodinger equation / Operator Theory: Advances and Applications, Vol. 114 / Birkhauser Verlag: Basel/Switzerland, 2000. - P.15-19.

14. Асташова И.В. Равномерные оценки положительных решений квазилинейных дифференциальных уравнений // Известия РАН. - 2008. - 72;6. - C.103-124.

15. Astashova I.V. Estimates of Solutions to One-dimensional Schrodinger Equation / World Scientific: Progress in Analysis, v. II / Proceedings of the 3rd International ISAAC Congress. -Singapore, 2003. - P.955-960.

16. Асташова И.В. Качественные свойства решений квазилинейных обыкновенных дифференциальных уравнений / М: МЭСИ, 2010. - 242 с. (ISBN 978-5-7764-0647-8)

ON ASYMPTOTIC BEHAVIOR OF SOLUTIONS TO NONLINEAR EQUATIONS

WITH COMPLEX COEFFICIENTS

I.V. Astashova

Moscow Lomonosov State University,

Leninskie gory, GSP-1, Moscow, 119991, Russia, e-mail: ast@diffiety.ac.ru

Abstract. Asymptotic formulas for modulus and argument of solutions and uniform estimates of solutions are obtained to nonlinear differential equations of Emden-Fowler’s type with complex coefficients.

Key words: asymptotic behavior, uniform estimates of solutions, complex coefficients.