F


Степенные ряды и нестепенные асимптотики решений второго уравнения Пенлеве

( Power series and nonpower asymptotics of solutions to the second Painlevé equation
Preprint, Inst. Appl. Math., the Russian Academy of Science)

Брюно А.Д., Завгородняя Ю.В.
(A.D.Bruno, Yu.V.Zavgorodnyaya)

ИПМ им. М.В.Келдыша РАН

Москва, 2003
Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (проект № 02-01-01067)

Аннотация

С помощью степенной геометрии, кратко изложенной в § 1, в случае общего положения вычислены все степенные разложения решений второго уравнения Пенлеве вблизи точек z=0, z = Ґ (§ 2) и z = z0 0 (§ 3). Аналогично при a=0 вычислены все степенные разложения решений и логарифма решений (§ 4). Обнаружены новые тонкие свойства некоторых из этих разложений.

Abstract

By means of Power Geometry, shortly presented in § 1, in the generic case we compute all power expansions of solutions to the second Painlevé equation at points z=0, z = _ (§ 2) and z = z0 0 (§ 3). Analogously for a=0 we compute all power expansions of solutions and of logarithm of solutions (§ 4). We have found new fine properties of some of these expansions.

E-mail:     bruno@keldysh.ru, zavgorodnyaya@rambler.ru

О степенной геометрии
§ 0. О степенной геометрии
Здесь приводятся те сведения из степенной геометрии [1-7], которые используются в этой работе.
1.1. Постановка задачи. Пусть x - независимая и y - зависимая переменные, x,y О \Bbb C. Положим X=(x,y). Дифференциальным мономом a(x,y) называется произведение обычного монома
cxr1yr2    def
=
 
  cXR,
(1)
где c=const О \Bbb C, R=(r1,r2) О \Bbb R2, и конечного числа производных вида
dly/dxl,    l О \Bbb N.
(2)
Сумма дифференциальных мономов
f(X)= е
ai(X)
(3)
называется дифференциальной суммой.
Основная задача. Пусть задано обыкновенное дифференциальное уравнение
f(X)=0,
(4)
где f(X) - дифференциальная сумма, в которую y входит в целых неотрицательных степенях. При x® 0 и при x®Ґ для решений y=j(x) уравнения (4) найти все степенные разложения
y=crxr+ е
cs xs,
(5)
где коэффициенты cr,cs=const О \Bbb C, cr 0, показатели степени r,s О \Bbb R и не имеют в \Bbb R точек накопления; r < s, если x® 0, и r > s, если x®Ґ.
Положим
w = м
п
н
п
о
-1,
если
x® 0,
1,
если
x®Ґ.
(6)
Степенная геометрия дает теорию и алгоритмы для решения этой задачи. Процедура вычисления разложений (5) состоит из двух этапов: вычисление первых приближений
y=crxr,    cr 0
(7)
и вычисление дальнейших членов разложения (5), которые удобно разбить на несколько ветвящихся шагов, описанных ниже.
1.2. Выделение укороченных уравнений. Каждому дифференциальному моному a(X) ставится в соответствие его (векторный) показатель степени Q(a)=(q1,q2) О \Bbb R2 по следующим правилам. Для монома (1) Q(cxr1yr2)=(r1,  r2); для производной (2) Q(dly/dxl)=(-l,1); при умножении дифференциальных мономов их показатели степени складываются как векторы Q(a1a2)=Q(a1)+Q(a2). Множество S(f) показателей степеней Q(ai) всех дифференциальных мономов ai(X), входящих в дифференциальную сумму (3), называется носителем суммы f(X). Очевидно, S(f) О \Bbb R2. Через fQ(X) обозначим сумму тех мономов ai(X) из (3), у которых Q(ai)=Q. Тогда дифференциальную сумму (3) можно записать в виде
f(X)= е
fQ(X) по Q О S(f).
(8)
Выпуклая оболочка G(f) носителя S(f) называется многоугольником суммы f(X). Граница G(f) многоугольника G(f) состоит из вершин Gj(0) и ребер Gj(1). Их называют (обобщенными) гранями Gj(d), где верхний индекс указывает размерность грани, а нижний - ее номер. Каждой грани Gj(d) соответствуют граничное подмножество Sj(d)=S(f)ЗGj(d) множества S и укороченная сумма
^
f
 
(d)
j 
(X)= е
ai(X)    по    Q(ai) О Sj(d).
(9)
Пусть плоскость \Bbb R*2 сопряжена плоскости \Bbb R2 так, что для P=(p1,p2) О \Bbb R*2 и Q=(q1,q2) О \Bbb R2 определено скалярное произведение
бP,Qс    def
=
 
  p1q1+p2q2.
(10)
Каждой грани Gj(d) в \Bbb R*2 соответствуют: свое нормальное подпространство
N(d)j={PбP,Qс = бP,Qўс    по    Q,Qў О Sj(d)}
и свой нормальный конус
Uj(d) = м
п
н
п
о
P:
бP,Qс = бP,Qўс,
Q,Qў О Sj(d),
бP,Qс > бP,Q"с,
Q" О S(f)\Sj(d)
ь
п
э
п
ю
.
Для ребра Gj(1) нормальное подпространство N(1)j это прямая, ортогональная ребру Gj(1) и проходящая через начало координат P=0, а нормальный конус это тот луч прямой Nj(1), который направлен от ребра Gj(1) наружу многоугольника G(f). Для вершины Gj(0) нормальное подпространство N(0)j=\Bbb R*2, а нормальный конус Uj(0) это открытый сектор (угол) на плоскости \Bbb R*2 с вершиной в нуле P=0 и ограниченный лучами, являющимися нормальными конусами ребер, примыкающих к вершине Gj(0).
Предложение 1.1. Если носитель S(f) состоит из одной вершины G1(0), то нормальный конус U1(0)=\Bbb R*2 т.е. в U1(0) лежат оба вектора ±R. В противном случае раствор угла Uj(0) меньше p, т.е. из двух векторов ±R в нормальный конус Uj(0) попадает только один. Здесь вектор R 0, R О \Bbb R*2.
Предложение 1.2. Если многоугольник G(f) состоит из одного ребра G1(1), то нормальный конус U1(1) это прямая, т. е. в U1(1) лежат оба вектора ±R. В противном случае в нормальном конусе Uj(1) лежит только один из векторов ±R. Здесь вектор R 0 является нормалью к ребру Gj(1).
Итак, каждой грани Gj(d) соответствуют: нормальное подпространство Nj(d) и нормальный конус Uj(d) в плоскости \Bbb R*2 и укороченное уравнение
^
f
 
(d)
j 
(X)=0.
(11)
Для разложения решения (5), которое запишем как y-j(x)=0, имеются аналогичные объекты, которые обозначаются малыми буквами: носитель s(y-j), многоугольник g(y-j). Из граней многоугольника g(y-j) нас интересует только ребро g1(1), которое соответствует укорочению (7) разложения (5), т.е. натянуто на точки Q1=(0,1) и Q2=(r,0). Через n и u обозначаются нормальное подпространство и нормальный конус ребра g1(1). Очевидно, n - это прямая, натянутая на вектор R=(1,  r), а u - это ее луч {P=lw(1,r), l > 0}, направленный наружу многоугольника g(y-j) и ортогональный ребру g1(1)
Теорема 1.1 [2, гл. VI, теорема 1.1]. Если разложение (5) удовлетворяет уравнению (4) и u М Uj(d), то укорочение (7) решения (5) является решением укороченного уравнения (9), (11).
Поэтому для нахождения всех укороченных решений (7) уравнения (4) надо вычислить: носитель S(f), многоугольник G(f), все его грани Gj(d) и их нормальные конусы Uj(d). Затем для каждого укороченного уравнения (9), (11) надо найти все его решения (7), у которых вектор R=(1,r) лежит в Nj(d), и из них отобрать те, у которых один из векторов ±(1,  r) лежит в нормальном конусе Uj(d). Если d=0, то это означает, что один из векторов ±(1,r)=wR лежит в Uj(d). Если d=1, то это свойство всегда выполнено. При этом согласно предложениям 1.1 и 1.2 определяется также значение w.
1.3. Решение укороченного уравнения. Здесь рассмотрим по отдельности два случая: вершины Gj(0) и ребра Gj(1).
Вершине Gj(0)={Q} соответствует укороченное уравнение (11) с точечным носителем Q. Положим g(X)=X-Q[^(f)]j(0)(X), тогда решение (7) уравнения (11) удовлетворяет уравнению g(X)=0. Подставляя y=cxr в g(X), получаем, что g(x,cxr) не зависит от x и c, и является многочленом от r, т.е. g(x,cxr) є c(r), где c(r) - характеристический многочлен дифференциальной суммы [^(f)]j(0)(X). Следовательно, для решения (7) уравнения (11) показатель r является корнем характеристического уравнения
c(r)    def
=
 
  g(xxr)=0,
(12)
а коэффициент cr - произвольный. Из корней ri уравнения (12) надо отобрать только те, для которых один из векторов w(1,r), где w = ±1, лежит в нормальном конусе Uj(0) вершины Gj(0). При этом согласно предложению 1.1 значение w определяется однозначно. Соответствующие выражения (7) с произвольной константой cr являются кандидатами на роль укороченных решений уравнения (4). При этом если w = -1, то x® 0 и в разложении (5) s > r, а если w = 1, то x®Ґ и в (5) s < r.
Укороченное уравнение (11) называется алгебраическим, если оно не содержит производных.
Замечание 1.1 [2]. Если укороченное уравнение (11) с d=0 является алгебраическим, то оно не имеет решений вида (7). Поэтому укорочения, состоящие из одного алгебраического монома, можно не рассматривать.
Ребру Gj(1) соответствует укороченное уравнение (11), нормальный конус Uj(1) которого является лучом { P = lwў(1,rў), l > 0}, где вектор wў(1,rў) является внешней нормалью ребра Gj(1). Нормальный конус u={P=lw(1,r), l > 0} решения (5) пересекается с Uj(1) только если эти лучи совпадают, т.е. w = wў и r=rў. Согласно предложению 1.2 этим однозначно определяется показатель степени r укороченного решения (7) и значение w. Для определения коэффициента cr надо выражение (7) подставить в укороченное уравнение (11). После сокращения на некоторую степень x получаем алгебраическое уравнение для коэффициента cr
~
~
f  
 
(c)    def
=
 
  x-s
^
f
 
(1)
j 
(x,crxr)=0.
(13)
Каждому его корню cr=cr(i) 0 соответствует свое выражение (7), которое является кандидатом на роль укороченного решения уравнения (4). При этом если получилось, что w = -1, то x® 0 и в (5) s > r, а если w = 1, то x®Ґ и в (5) s < r.
Итак, каждое укороченное уравнение (11) имеет несколько подходящих решений (7) с u М Uj(d). Объединим их в непрерывные по w, r, cr и параметрам уравнения (4) семейства, которые обозначим Fj(d)k, где k=1, 2,.... Если нас интересуют не все решения (5) уравнения (4), а только те решения (5), у которых нормальный конус u лежит в некотором заданном конусе K, то K называется конусом задачи [2, гл. I, § 6]. Например, для укороченного уравнения (11) нормальный конус Uj(d) является конусом задачи, если нет других ограничений.
1.4. Критические числа укороченного решения. Если найдено укороченное решение (7), то замена
y=crxr+z,
(14)
приводит уравнение (4) к виду
f(x,cxr+z)    def
=
 
  
~
f
 
(x,z)=0,
(15)
где [(f)\tilde](x,z) - дифференциальная сумма, все точки Q=(q1,q2) ее носителя S([(f)\tilde]) имеют целую неотрицательную координату q2. К уравнению (15) можно применить описанные выше вычисления (т.е. носителя, многоугольника, укорочений и т.д.) и получить следующий член cs0xs0 разложения (5), у которого s0 > r, если x® 0, и s0 < r, если x® Ґ. Следовательно, получилась основная задача для уравнения (15), но теперь с конусом задачи
K= м
п
н
п
о
P=(p1p2):
p2/p1 > r,  p1 < 0,
если
x® 0,
p2/p1 < r,  p1 > 0,
если
x®Ґ
ь
п
э
п
ю
.
Согласно (6) для s=p2/p1 этот конус задачи K можно записать в виде
K={s=p2/p1:  sw < rw}.
(16)
Однако во многих случаях дифференциальная сумма [(f)\tilde](x,z) имеет специальный вид, что позволяет существенно сократить вычисления разложений (5). Предположим, что уравнение (15) имеет вид
~
f
 
(x,z)    def
=
 
  L(x)z+h(x,z)=0,
(17)
где L(x) - линейный дифференциальный оператор и носитель S(Lz) состоит из одной точки (v,1), являющейся вершиной [(G)\tilde]1(0) многоугольника G([(f)\tilde]), у всех точек Q=(q1,q2) носителя S(h) координата q2 і 0 и нет точки Q=(v,1), нормальный конус вершины [(G)\tilde]1(0) содержит вектор P=(p1,p2) с p1w > 0. Напомним, что дифференциальная сумма f(x,y) имеет первую вариацию [8] df(x,y)/dy, которая обладает следующими свойствами.
d(c xq1yq2)/dy=cq2xq1yq2-1,    d(dl y/d xl)/dy=dl/dxl,

d(f+g)/dy=df/dy+dg/dy,    d(fg)/dy=(df/dy)g+f(dg/dy).
Теорема 1.2 [4, 5, 6]. Пусть (7) - решение укороченного уравнения (11) с u О U(d)j. Тогда в уравнении (17) оператор
L(x) = d
^
f
 
(d)
j 
(x,y)/dy     на    y=crxr,
(18)
т.е. равен первой вариации, вычисленной на кривой (7).
Следовательно, после подстановки (15) уравнение (4) принимает вид (17), если L(x)\not є 0. Пусть n(k) - характеристический многочлен дифференциальной суммы L(x)z, т.е.
n(k)=x-v-kL(x)xk.
(19)
Корни k1,...,ks многочлена n(k) называются собственными значениями укороченного решения (7). Те из вещественных собственных чисел ki, которые лежат в конусе задачи, т.е. удовлетворяют неравенствам (16), называются критическими числами. Они играют важную роль при нахождении разложения (5), что будет показано в следующем пункте.
Замечание 1.2 [5]. Степенное решение (7) алгебраического укороченного уравнения (11) с d=1 не имеет собственных значений и критических чисел, ибо для него
n(k) є n0=const=
^
f
 
(1)
j 

y
(1,cr).
(20)
Если cr - простой корень уравнения (13), то n0 0. Если cr - кратный корень уравнения (13, то n0=0.
Если L(x)\not є 0, то n(k)\not є 0. Если же L(x) є 0, то для уравнения (15) надо решать основную задачу как описано в пп. 1.2 и 1.3 с учетом конуса задачи (16).
1.5. Носитель разложения решения (5). Рассмотрим уравнение (17) в невырожденном случае, когда n(k)\not є 0. Параллельно сдвинем носитель S([(f)\tilde]) на вектор -(v,1). Тогда вершина (v,1), соответствующая члену L(x)z, сдвинется в начало координат. Пусть Sў+ - множество конечных сумм векторов Q О Sў[(   def) || ( = )]  S([(f)\tilde])-(v,1) и пусть
K=Sў+ З{q2=-1}.
(21)
Пусть æ - число критических значений ki укороченного решения (7).
Рассмотрим по отдельности случаи с æ = 0 и æ = 1.
Случай æ = 0.
Теорема 1.3 [7]. В случае æ = 0 имеется разложение (5) решения уравнения (4), в котором все показатели s лежат в множестве K и все коэффициенты cs определяются однозначно.
Лемма 1.1. Пусть G(0)j - такая вершина многоугольника G(f) уравнения (4), что укороченное уравнение [^(f)](0)j(X)=0 имеет решение (7), вектор w(1,  r) О U(0)j и все точки сдвинутого носителя S(f)-G(0)j представляются в виде lR*+mR*, где целые l,m і 0, R* и R* - некоторые вектора. Тогда для множества (21) уравнения (15) справедливо включение
K М {k=r+lr*+mr*, целые l,m і 0, l+m > 0 },
где r*=б(1,  r),  R*с, r*=б(1,  r),  R*с.
Точечное множество Z плоскости \Bbb R2 называется решеткой, если для любых двух векторов Q,  Qў О Z их сумма Q+Qў и разность Q-Qў принадлежат Z. Векторы B1 и B2 образуют базис решетки Z, если любая точка Q решетки Z представляется в виде Q=n1B1+n2B2, где n1 и n2 О \Bbb Z, т.е. целые числа. Пусть S - некоторое множество в \Bbb R2 и точка Q О S. Множество точек S-Q будем называть сдвинутым множеством S. Для укороченного решения (1.6) носитель s состоит из двух точек Q1=(0, 1) и Q2=(r, 0). Поэтому сдвинутый носитель s состоит из точек 0 и Q2-Q1=(r-1).
Замечание 1.3. Если сдвинутые носители S(f) уравнения (4) и s укороченного решения (7) лежат в некоторой решетке Z, то множество
K М ZЗ{q2 = -1}ЗK.
(22)
Случай æ = 1. Пусть k1 - единственное критическое значение. Обозначим через Sў+(k1) множество конечных сумм векторов Q О Sў[(   def) || ( = )]  S([(f)\tilde])-(v,1) и вектора (k1,-1), и положим
K(k1)=Sў+(k1)З{q2=-1}.
(23)
Теорема 1.4 [7]. Если в случае æ = 1 единственное критическое значение k1 не лежит в множестве K, то имеется разложение (5) решения уравнения (4), где показатели s лежат в множестве K(k1), при этом коэффициент ck1 произволен, а остальные коэффициенты cs определяются однозначно.
Лемма 1.2. Если множество (21) имеет вид
K= м
н
о
k=r+l
~
r
 
, целые l > 0 ь
э
ю
,
то множество
(23) имеет вид
K(k1)= м
н
о
k=r+l
~
r
 
+m(k1-r), целые l,m і 0,  l+m > 0 ь
э
ю
.
   Замечание 1.4. Если множество Sў=S([(f)\tilde])-(v,1) и точка (k1,-1) лежат в некоторой решетке Z(k1), то
K(k1) М Z(k1)З{q2=-1}ЗK.
(24)
Теорема 1.5 [7]. Множества (21) и (23) не имеют точек накопления на прямой \Bbb R, если множество S([(f)\tilde]) не имеет точек накопления на плоскости \Bbb R2.
Если в ситуации теоремы 1.4 условие на критическое значение k1 не выполнено (т. е. k1 О K), то надо вычислить начальный отрезок
y=crxr+ m
е
 
ckxk+u
(25)
разложения (5) до такого значения m, что |m-r| і |k1-r|. В уравнении для u условие теоремы 1.4 будет выполнено.
1.6. Вычисление коэффициентов разложения. Упорядочим точки s множества K или K(k1) по росту |s-r|. Двигаясь по точкам s из соответствующего множества в соответствии с этой упорядоченностью, для каждого коэффициента cs из (5) получаем линейное уравнение L cs xs+bsxs+v=0, т. е.
n(s)cs+bs=0,
(26)
где bs зависит от предыдущих коэффициентов cl с индексами l, предшествующими индексу s. Если n(s) 0, то cs=-bs/n(s). Если n(s)=0 и bs=0, то cs - произвольно. Если n(s)=0 и bs 0, то уравнение (26) не имеет решения.
   Замечание 1.5. Если n(s)=0, то получается условие разрешимости (совместности) уравнения (26)
bs=0.
(27)
Если исходное уравнение (4) зависит от параметров A, то, вообще говоря, коэффициент bs также зависит от этих параметров, т.е. bs=bs(A). При этом возможны три случая: (a) bs(A) є 0 при любых значениях параметров A; (b) условие (27) является некоторым условием на параметры A; (c) bs(A) 0 при любых допустимых значениях параметров A.
Для каждого семейства Fj(d)k получается некоторое семейство разложений (5) с определенным носителем, в котором для критических коэффициентов cki могут существовать условия разрешимости. Если они выполнены, то эти коэффициенты произвольны, а остальные определяются однозначно. Такие семейства разложений будем обозначать Gj(d)k.
1.7. Вычисление второго приближения. Пусть известны укороченное уравнение (11), его подходящее решение (7) и его критические числа ki. Покажем, как находится ближайшая к r точка s0 множества K. Пусть вектор P=w(1,  r), где w = ±1 (см. (1.5ў)), направлен наружу многоугольника G(f). Пусть m = бP,  Qс и mў=maxбP,  Qўс, где Q О S(d)j,  Qў О S(f)\S(d)j.
Теорема 1.6 [4, 5]. Ближайшая к r точка s0 такова, что |s0-r|=m-mў.
Если критические значения ki отсутствуют или таковы, что |ki-r| > m-mў, то можно вычислить коэффициент cs0 следующим образом. Выделим в сумме f(x) второе приближение [^([^(f  )])] =еfQ(X) по Q:  бP,  Qс = mў, делаем подстановку y=crxr+cs0xs0 в сумму
^
f
 
(d)
j 
+
^
^
f  
 
(28)
и приведем подобные члены. Приравнивая нулю коэффициент при xmўw в сумме (28), получаем линейное уравнение (26) для коэффициента cs0, где bs это коэффициент суммы [^([^(f  )])](x,crxr) при xmўw и s=s0.
Если имеется критическое значения k1 такое, что |k1-r| < m-mў, то s0=k1 и коэффициент cs0 произвольный.
1.8. Вопросы сходимости. Максимум порядков l всех производных (2), входящих в дифференциальный моном a(X), называется его порядком дифференцирования и обозначается p(a). Порядком дифференцирования p(f) суммы (3) называется наибольший из порядков дифференцирования p(ai) ее мономов ai(X).
Условие 1.1. В сумме (17) порядок дифференцирования p([(f)\tilde]) Ј p(Lz).
Теорема 1.7 [4, 5]. Пусть в ситуациях теорем 1.3, 1.4 выполнено условие 1.1, тогда разложение (5) абсолютно сходится для достаточно малых |x| 0, если x® 0,   и |x|-1 0, если x® Ґ (т.е. для |x|-w 0).
1.9. Нестепенные асимптотики решений. В п.п. 1.1 - 1.4 изложен метод вычисления всех степенных асимптотик, т.е. вида (7). Поэтому здесь сосредоточимся на вычислении всех нестепенных асимптотик решений, т.е. не имеющих вида (7). Порядком функции j(x) называется число r=lim(ln|j(x)|/ln|x|) О [-Ґ,+Ґ], если предел существут. Две функции j(x) и y(x) сильно (асимптотически) эквивалентны, если j(x)=y(x)[1+o(|x|we|y(x)|w2e)] для некоторого e < 0, где
w2 = м
п
н
п
о
-1,
если    y(x)® 0,
0,
если    y(x)®const 0,
1,
если    y(x)®Ґ.
При этом функция j(x) является сильной асимптотикой функции y(x) и наоборот. Пусть y=j(x) - решение уравнения (4), и функция j(x) имеет порядок r. Тогда нормальным конусом u этого решения является луч
u = м
п
н
п
о
lw(0,-1),
если    r=-Ґ;
lw(1,r),
если    r О \Bbb R;
lw(0,1),
если    r=+Ґ,    l > 0.
Понятие конус задачи применимо и к решениям с так определенным нормальным конусом u.
Теорема 1.8 [2, гл. VI, теорема 1.1]. Если решение y=j(x) уравнения (4) имеет нормальный конус u, то укороченное уравнение (11), для которого u М Uj(d), имеет сильно асимптотически эквивалентное решение y = y(x).
1.10. Случай горизонтального ребра G(1)j. В этом случае у всех точек Q=(q1,q2) ребра Gj(1) координата q2 одинакова. Положим
g(x,y)    def
=
 
   y-q2
^
f
 
(1)
j 
(x,y)
(29)
и сделаем логарифмическое преобразование
h = dlny/dx.
(30)
Согласно теореме 2.4 из [2, гл. VI] при этом преобразовании дифференциальная сумма g(x,y) перейдет в дифференциальную сумму h(x,h)[(   def) || ( = )]  g(x,y) и уравнение (11) примет вид
h(x,h)=0.
(31)
Теперь для уравнения (31) получаем основную задачу с конусом задачи p1+p2 і 0. Сумму порядков всех производных, входящих в дифференциальный моном a(x,y), назовем суммарным порядком дифференцирования монома a(x,y) и обозначим D(a). Для дифференциальной суммы (3) суммарный порядок дифференцирования D(f)=maxiD(ai). Пусть ребро Gj(1) соединяет вершины Gj-1(0)=(q1ў,q2ў) и Gj(0)=(q1ўў,q2ў) с q1ў < q1ўў.
Пусть Wj - степень характеристического многочлена cj(r) укорочения [^(f)](0)j, соответствующего вершине G(0)j.
Теорема 1.9 [7, п. 5.4]. Если Wj-1 = D([^(f)]j-1(0))=D([^(f)]j(1)) (или D([^(f)]j(0))=D([^(f)]j(1)) = Wj), то уравнение (11) не имеет подходящих решений при x® 0 (или x®Ґ).
Следовательно, если Wj-1 = D([^(f)]j-1(0))=D([^(f)]j(1)) = D([^(f)]j(0)) = Wj, то преобразование (30) и дальнейшее исследование делать не надо.
§ 0. Степенные разложения решений в окрестности нуля и бесконечности при a 0
Для второго уравнения Пенлеве
f(z,w) def
=
 
-w" + 2w3 + zw + a = 0,
(32)
где a - комплексный параметр, при a 0 найдем все степенные разложения
w = crzr + е
cszs
(33)
его решений, где коэффициенты cr,cs = const О \Bbb C, cr 0, показатели степени r,s О \Bbb R: r < s, если z® 0, и r > s, если z®Ґ.
2.1. Основные объекты степенной геометрии. Пусть параметр a 0, тогда носитель S(f) уравнения (32) состоит из четырех точек:
Q1 = (-2,1),     Q2 = (0,3),     Q3 = (1,1),     Q4 = (0,0).
(34)
Их выпуклая оболочка G(f) это четырехугольник с вершинами G(0)j = Qj и ребрами G(1)j = {Qj,Qj+1}, где j=1,2,3,4, Q5=Q1 (см. рис. 1).
Внешние нормали Nj к ребрам G(1)j суть N1 = (-1,1),    N2 = (1,1/2),     N3 = (1,-1),     N4 = (-1,-2). Нормальные конусы U(1)j ребер G(1)j суть (см. рис. 2):

U(1)1 = {P:p2 = 2p1 < 0} = {lN1},
U(1)2 = {P: p2 = -p1 > 0} = {lN2},
U(1)3 = {P: p2 = p1/2 > 0} = {lN3},
U(1)4 = {P: p2 = -p1 < 0} = {lN4}, l > 0.
(35)
Нормальные конусы Uj(0) вершин Gj(0) это углы между нормальными конусами прилегающих к вершине ребер (см. рис. 2):

U(0)1 = {P: p2 < -p1,    p2 > 2p1 },
U(0)2 = {P: p2 > p1/2,    p2 > -p1 },
U(0)3 = {P: p2 > -p1,    p2 < p1/2 },
U(0)4 = {P: p2 < -p1,    p2 < 2p1 }.
(36)
Теорема 2.1. Носитель S(f) уравнения (32) лежит в решетке Z с базисом B1 = (1,1), B2 = (2,-1).
Доказательство. Согласно (34) имеем Q1 = -B1,    Q2 = 2B1 - B2,     Q3 = B1,    Q4 = 0, т.е. все точки носителя S(f) имеют вид mB1 + nB2 с целыми m и n. При этом точки -Q1 и Q3 совпадают с базисными векторами решетки. [¯]
2.2. Решения, соответствующие вершинам. Рассмотрим укороченное уравнение, соответствующее вершине G(0)1:
^
f
 
(0)
1 
def
=
 
- w" = 0
(37)
и найдем все его подходящие укороченные решения
w = crzr,    cr 0
(38)
с w(1,r) О U(0)1. Согласно (36) p1 < 0 в конусе U(0)1. Следовательно, w = -1, z® 0 и разложения будут по возрастающим степеням z.
Размерность грани d=0, поэтому согласно п. 1.3 положим
g(z,w) = X-Q4
^
f
 
(0)
4 
= z2w
^
f
 
(0)
4 
= z2w-1(-w").
Для w = zr получим характеристический многочлен c(r)[( def) || =]g(z,zr) = r(1 -r), т.е. r1 = 0 и r2 = 1 суть корни уравнения c = 0. Рассмотрим их по отдельности.
Корню r=r1=0 соответствует вектор R=(1,0) и согласно (36), вектор wR О U(0)1. Итак, получили семейство F(0)11 укороченных решений w=c0 с произвольной постоянной c0 0 и w = -1.
Найдем первую вариацию:
d
^
f
 
(0)
4 
/dw = -d2/dz2.
(39)
Согласно теореме 1.2 получаем оператор
L(z) = -d2/dz2\not є 0.
(40)
Его характеристический многочлен n(k) имеет вид
n(k) = z2-kL(z)zk = k(1-k).
(41)
Уравнение n(k) = 0 имеет 2 корня: k1 = 0 и k2 = 1. В нашем случае w = -1 и r=0, следовательно, конус задачи (16) есть K = {s > 0}. Поэтому имеется одно критическое значение k2 = 1 > 0 и æ = 1. Согласно п. 1.5. найдем множество K. Сдвинутый носитель решения w=c0 это вектор (0,1). Вместе с векторами B1 = (1,1) и B2 = (2,-1) он порождает целочисленную решетку \Bbb Z2. Но множество Sў лежит в секторе, ограниченном лучами, натянутыми на векторы B1 и B2. Его целочисленные точки, лежащие на прямой q2 = -1, образуют множество K. Следовательно,
K = {s = l, целые l > 1}.
(42)
Так как число критических чисел æ = 1 и k2 П K, то по теореме 1.4 показатели s лежат в множестве K(k2) = K(1), при этом коэффициент c1 произволен, а остальные коэффициенты cs определяются однозначно. Множество K(1) находим, используя лемму 1.2:
K(1) = {s = l, целые l і 1}.
(43)
Таким образом, получаем, что разложение (33) в случае r=0 есть
w = c0 + c1z + Ґ
е
k=2 
ckzk = Ґ
е
k=0 
ckzk,
(44)
где c0 = const 0 и c1 - произвольны, а ck, k > 1, однозначно определены. Это семейство разложений G(0)11.
Рассмотрим теперь случай r=r2=1. Здесь вектор P=w(1,1) = -(1,1) О U(0)1. Конус задачи K = {s > 1}. В этом случае первая вариация снова есть (39) и не зависит от укороченного решения (38), тогда по теореме 1.2, учитывая наше решение w=cz, получаем, что оператор L(z) снова есть (40). Характеристический многочлен n(k) имеет вид (41). Опять корни уравнения n(k) = 0 равны k=0 и k=1. Но ни одно из этих собственных значений не лежит в конусе задачи K = {s > 1}. Следовательно, число критических чисел æ = 0. По теореме 1.3 все показатели s разложения (33) лежат в множестве K и все коэффициенты cs определяются однозначно. Используя замечание 1.3, найдем множество K. Теперь сдвинутый носитель укороченного решения w=cz это вектор (1,-1). Он не лежит в решетке теоремы 2.1. Вместе с векторами B1, B2 он порождает целочисленную решетку \Bbb Z2. Согласно замечанию 1.3 множество K есть
K = {s = 1+ m, целые m > 0}.
(45)
Следовательно, решение будет выглядеть следующим образом:
w = cz + Ґ
е
k=2 
ckzk = Ґ
е
k=1 
ckzk,
(46)
где произвольное c = const = c1 0, а остальные ck определяются однозначно. Вычислим второе приближение. Так как вектор P = w(1,r) = -(1,1), то m = 1, а mў = 0. Поэтому ближайшей точкой s0 к r=1 будет s0 = 2. В данном случае критических чисел нет, поэтому можем вычислить коэффициент c2 по формуле (26). Выделим в сумме f(X) второе приближение [^([^(f  )])] = еfQ(X)    по     Q:бP,Qс = mў, т.е. [^([^(f  )])] = a = b2. Согласно (26)
c2 = -
^
^
f  
 
/n(2) = -a/2,
(47)
т.к. по (41) n(2)=2. Таким образом, получили семейство G(0)12.
Вершинам G(0)2, G(0)3 и G(0)4 соответствуют алгебраические укороченные уравнения. Согласно замечанию 1.1 они не имеют решений вида (38).
2.3. Решения, соответствующие ребру G(1)1. Ребру G(1)1 соответствует укороченное уравнение
^
f
 
(1)
1 
def
=
 
-w" + 2w3
(48)
с нормальным конусом U(1)1 = {P = l(-1,1) = -l(1,-1),     l > 0}. Поэтому здесь r=-1 и w = -1. Для уравнения (48) будем искать решение w=cz-1, c 0. Подставляя это выражение в уравнение (48), для постоянной c получаем алгебраическое уравнение 2c3 -2c = 0. Оно имеет два ненулевых решения c = ±1. Итак, получили два укороченных решения w = ±z-1. Вычислим их критические числа.
Первая вариация по w
d
^
f
 
(1)
3 
/dw = -d2/dz2 + 6w2,
(49)
тогда, учитывая наше решение w = cz-1 и то, что c2 = 1, получим:
L(z) = -d2/dz2 + 6z-2.
(50)
Характеристический многочлен n(k) имеет вид
n(k) = z 2-kL(z)zk = -k(k-1)+6 = -(k+2)(k-3).
(51)
Он имеет два корня k1 = -2 и k2 = 3. Конус задачи здесь
K = {s > -1}.
(52)
Следовательно, только k2 = 3 О K, т.е. имеется только одно критическое число и æ = 1. Согласно п. 1.5 вычислим множество K показателей s в разложении (33). Для укороченного уравнения w = c-1z-1 сдвинутый носитель содержит точку (-1,-1) = -B1, которая лежит в решетке Z, описанной в теореме 2.1. Эта решетка пересекается с прямой q2=-1 по точкам q1 = s вида
s = -1+3l,
(53)
где l - целое число. Согласно (52) и замечанию 1.3 получаем, что множество K состоит из точек (53) c l > 0. Критическое число k2 = 3 не лежит в множестве K, тогда по теореме 1.4 все показатели решения s лежат в множестве K(k2) = K(3) , причем коэффициент ck3 = c3 - произволен, а остальные cs однозначно определяются. По лемме 1.2 множество
K(3) = {s = 3l + 4m -1, целые     l,m і 0,    l+m > 0}.
(54)
Найдем второе приближение решения. Для этого сначала определим значения m и mў. В нашем случае имеем вектор P = w(1,r) = (-1,1), следовательно, получаем, что m = 3, а mў = 0. Тогда получаем, что s0 = 2. Так как у нас единственное критическое число k2 = 3 и условие |k2 -r | > m-mў выполнено, то коэффициент c2 можем найти по формуле (26). Здесь [^([^(f  )])] = zw+a и b2 = [^([^(f  )])](z,cz-1) = a+c. Согласно (51) n(2)=4, следовательно, согласно (26) c2 = -(c+a)/4, где c = ±1.
Множество (54) состоит из точек s=2,3,5,6,..., т.е. 2,3 и целые s і 5, поэтому ребру G(1)1 соответствует степенное разложение решений
w = cz-1 + c2z2 + c3z3 + Ґ
е
s=5 
cszs,
(55)
где коэффициенты c = ±1, c2 = -(c+a)/4, c3 - произвольное, а остальные cs - однозначно определяются. Они образуют два однопараметрических семейства G(1)11 и G(1)12.
Однако, при c3 = 0 получаем разложения с показателями (53), т.е.
w = cz-1 + z-1 Ґ
е
k=1 
ckz3k.
(56)
2.4. Решения, соответствующие ребру G(1)2. Ребру G(1)2 соответствует укороченное уравнение
^
f
 
(1)
2 
def
=
 
zw + 2w3 = 0
(57)
с нормальным конусом U(1)2 = {P=l(1,1/2)}. Поэтому здесь r=1/2 и w = 1 и ищем решения уравнения (57) в виде w = cz1/2. Для c получаем алгебраическое уравнение c + 2c3 = 0, которое имеет два ненулевых корня c = ±i/Ц2. Т.к. w = 1, то z®Ґ, поэтому степенные разложения вида (33) являются разложениями по убывающим степеням zs с s < 1/2. Укороченное уравнение (57) является алгебраическим. По замечанию 1.2 его решения w=cz1/2 не имеют критических чисел (т.е. æ = 0) и для них
n(k) = n0 = z-1
^
f
 
(1)
2 

w
= z-1(z+6w2) = -2.
(58)
Найдем множество K значений показателей s в разложении (33). Сдвинутый носитель укороченного решения w = cz1/2 содержит точку B3 = (1/2,-1). Вместе с векторами B1 и B2 теоремы 2.1 она порождает решетку с базисом B1, B3. Она пересекается с прямой q2=-1 по точкам
q1 = s =  1

2
+  3

2
m,
(59)
где m - целые. Поскольку K = {s < 1/2}, то получаем
K = {  1

2
-  3

2
m, целые m > 0}.
(60)
По теореме 1.3 все показатели s лежат в множестве (60) и все коэффициенты cs определяются однозначно. Найдем коэффициент при втором приближении решения. Здесь [^([^(f  )])] = a, согласно (58) n(-1) = -2 и по (26) имеем c-1 = a/2.
Окончательно получаем, что ребру G(1)2 соответствуют два разложения
w = cz1/2 + c-1z-1 +z1/2 Ґ
е
s=2 
csz-3s/2,
(61)
где c=±i/Ц2, c-1 = -a/2, а все остальные cs однозначно определяются. Они образуют семейства G(1)21 и G(1)22.
2.5. Решения, соответствующие ребру G3(1). Ему соответствует укороченное уравнение
^
f
 
(1)
3 
= zw + a = 0
(62)
с нормальным конусом U(1)3 = {P=l(1,-1),   l > 0}. Следовательно, r=-1 и w = 1. Ищем решение w = c-1z-1. Подставляя его в (62), получаем c-1 = -a и c-1 0, т.к. по преположению a 0.
Т.к. w = 1, то z ® Ґ и степенные разложения вида (33) являются разложениями по убывающим степеням zs, s < -1, и конус задачи K выглядит следующим образом
K = {s < -1}.
(63)
Укороченное уравнение (62) является алгебраическим. По замечанию 1.2 его решения w=c-1z-1 не имеют критических чисел и для них
n(k) = n0 = z-1(
^
f
 
(1)
3 

w
) = 1.
(64)
По теореме 1.3 все показатели степени s разложения вида (33) лежат в множестве K и все коэффициенты cs определяются однозначно. Найдем множество K. Для укороченного уравнения w = c-1z-1 сдвинутый носитель содержит точку (-1,-1) = -B1, которая лежит в решетке Z, описанной в теореме 2.1. Эта решетка пересекается с прямой q2=-1 по точкам q1 = s вида (53). Согласно (63) и замечанию 1.3 множество
K = {-1-3l, целые     l > 0}.
(65)
Таким образом, получили, что разложение
w = c-1z-1 + Ґ
е
k=1 
ck z-3k-1
(66)
является одним из решений уравнения (32), где c-1 = -a, а остальные cs однозначно определяются. Найдем коэффициент во втором приближении. Для этого найдем m и mў. Т.к. P = w(1,r) = (1,-1), то m = 0, а mў = -3. Следовательно, ближайшей точкой к r=-1 является точка s = -4. Т.к. критических чисел нет, то мы можем найти коэффициент c-4 тем же способом, который использовался при нахождении коэффициента второго приближения в случае ребра G(1)2. Здесь [^([^(f  )])] = -w" + 2w3. Поэтому b-4 = -2a-2a3 и c-4 = 2a+2a3, согласно (64).
Окончательно получаем, что уравнение (32) имеет одно решение
w = c-1z-1 + Ґ
е
k=1 
c-3k-1z-3k-1,
(67)
где c-1 = -a, c-4 = 2a(a2 +1) и все остальные cs однозначно определяются, образующее семейство G(1)3.
2.6. Решения, соответствующие ребру G(1)4. Ребру G(1)4 соответствует укороченное уравнение
^
f
 
(1)
4 
def
=
 
-w" + a = 0
(68)
с нормальным конусом U(1)4 = {P: p2 = 2p1 < 0}, т.е. r=2 и w = -1. Будем искать решение w=cz2. Для постоянной c получим следующее значение: c = a/2 и c 0, т.к. по предположению, a 0. Т.к. w = -1, то z ® 0, и разложения вида (33) будут по возрастающим степеням zs с s > 2. Множество
K = {s > 2}
(69)
это конус задачи для ребра G(1)4. Первая вариация d[^(f)] / dw, оператор L(z) и функция n(k) здесь имеют вид (39), (40) и (41) соответственно. Поскольку корни k1=0 и k2 = 1 уравнения n(k)=0 не лежат в конусе задачи (69), то критических чисел нет. Найдем множество K показателей s разложения (33). Сдвинутый носитель укороченного решения w=cz2 это точка (2,-1) = B2, которая лежит в решетке Z теоремы 2.1. Эта решетка пересекается с прямой q2=-1 по точкам q1=s вида (53). Согласно (69) и замечанию 1.3 имеем
K = {s = 2+3l, целые l > 0}.
(70)
Причем, т.к. нет критических чисел, то по теoреме 1.4 все показатели s лежат в множестве (70). Найдем второе приближение. У нас m = 0, mў = -3, тогда ближайшей точкой к r = 2 будет s0 = 5. Т.к. æ = 0, то вычислим коэффициент c5 по формуле (26). Здесь [^([^(f  )])] = zw,     b5 = c, по (41) n(5)=-20, откуда по (26) имеем c5 = c/20 = a/40. Таким образом, ребру G(1)4 соответствует степенное разложение
w = cz2 + c5z5 + Ґ
е
l=2 
c2+3lz2+3l,
(71)
где c = a/2, c5 = a/40 и все остальные cs однозначно определяются. Это семейство обозначим G(1)4.
2.7. Сводка результатов и их обсуждение. Были найдены следующие разложения решений уравнения (32).
  1. Двупараметрическое (по c 0 и c1) семейство G(0)11 степенных разложений (44).

  2. Однопараметрическое (по c) семейство G(0)12 степенных разложений (46), которое включается в предыдущее семейство при c0 = 0.

  3. Два однопараметрических (по c3) семейства G(1)11 и G(1)12 степенных разложений (55), которые при c3 = 0 имеют вид (56).

  4. Два степенных разложения G(1)21 и G(1)22 (61).

  5. Cемейство G(1)3, состоящее из одного степенного разложения (66).

  6. Cемейство G(1)4 - одно степенное разложение (71).

Существование решений 1, 2 семейств G(0)11 и G(0)12 следует из теоремы Коши [11, гл. I, § 5, с. 18]. Решения 3 семейств G(1)11 и G(1)12 были найдены в [12]. Эти решения также были получены в [13, гл. I, § 4, пример 4.6] с помощью перехода к системе двух уравнений и с использованием специальной сложной программы, проверяющей выполнение условия совместности. Решения 4 и 5 семейств G(1)21, G(1)22 найдены в [14]. Решение в семействе G(1)4 и разложения (56) найдены, по-видимому, впервые. В [11, гл. III, § 10, с. 183, 184] утверждается, что второе уравнение Пенлеве не имеет решений, которые вблизи точки z=Ґ разлагаются в ряд вида w = akzk+ ak+1zk-1 + ... Но это не верно, ибо разложение 5 (т.е. (66)) является именно таким с k=-1. Отметим еще, что согласно теореме 1.7 разложения 1, 2, 3, 6 сходятся для достаточно малых |z|, а 4, 5 могут расходиться при любых z.
§ 0. Степенные разложения решений в окрестности точки z0 0
Рассмотрим второе уравнение Пенлеве (32) в окрестности точки z0 0, т.е. в исходном уравнении (32) сделаем замену z = [(z)\tilde] + z0:
f(
~
z
 
,x) def
=
 
-w" + 2w3 +
~
z
 
w + z0w +a,
(72)
где a - комплексный параметр и a 0. Найдем все степенные разложения вида (33) решений уравнения (72), где z заменено на [(z)\tilde].
3.1. Основные объекты степенной геометрии. Носитель S(f) уравнения (72) состоит из пяти точек:
Q1 = (-2,1), Q2 = (0,3), Q3 = (1,1), Q4=(0,0), Q5 = (0,1)
(73)
и показан на рис. 3. В данном случае ребра, вершины (см. рис. 3) и нормальные конусы (см. рис. 2) для них те же, что были в случае § 2. Носитель S(f) уравнения (72) лежит в целочисленной решетке \Bbb Z2 и не лежит в более редкой решетке. Поскольку нас интересуют только решения с [(z)\tilde] ® 0, то конус задачи p1 < 0. С ним пересекаются только нормальные конусы U(0)1, U(0)2, U(0)4 и U(1)1, U(1)4. Поэтому далее рассматриваем только решения, соответствующие граням G(0)1, G(0)2, G(0)4, и G(1)1, G(1)4.
3.2. Решения, соответствующие вершинам. Случай вершины G(0)1 аналогичен уже рассмотренному в пункте 2.2. Т.о. получаем, что вершине G(0)1 соответствуют семейство G(0)11 (формула (44)) и семейство G(0)12 (формулы (46) и (47)) c заменой z на [(z)\tilde].
Вершинам G(0)2 и G(0)4 соответствуют алгебраические укороченные уравнения [^(f)](0)j=0. Согласно замечанию 1.1 они не имеют решений вида (38).
3.3. Решения, соответствующие ребру G(1)1. Ему соответствует укороченное уравнение
^
f
 
(1)
1 
def
=
 
-w" + 2w3.
(74)
Этот случай аналогичен уже рассмотренному в п. 2.3. Т.о. ребру G(1)1 соответствует нормальный конус U(1)1 = {P=l(-1,1),    l > 0}, поэтому w = -1 и r=-1. Будем рассматривать решение w=c[(z)\tilde]-1 уравнения (74), где c=±1. Первая вариация d[^(f)](1)1 /dw, оператор L([(z)\tilde]) и функция n(k) здесь имеют вид (49), (50) и (51) соответственно. Т.к. из корней k1 = -2 и k2 = 3 уравнения n(k)=0 в конусе задачи K = {s > -1} лежит только k2 = 3, то он и является критическим числом и æ = 1. Множество K в данном случае есть
K = {2l-1, целые     l > 0},
(75)
где l - целое. Т.к. критическое число k2 = 3 лежит в множестве (75), то надо вычислить начальный отрезок разложения
w = c
~
z
 
-1
 
+ c1
~
z
 
+ c2
~
z
 
2
 
+ c3
~
z
 
3
 
.
(76)
Найдем значение коэффициента c1. В этом случае, m = 3, mў=1, тогда ближайшей точкой к r=-1 будет точка s0 = 1. Тогда по формуле (26), получаем c1 = -cz0/6 = z0/6, т.к. c=±1. Чтобы найти коэффициент c2, в уравнение (72) подставим w = c[(z)\tilde]-1 + c1[(z)\tilde] + c2[(z)\tilde]2 и приравняем к нулю коэффициент при [(z)\tilde]0. Получаем сумму
-2c2 + 2(c
~
z
 
-1
 
+ c1
~
z
 
+ c2
~
z
 
2
 
)3 + c +c1
~
z
 
2
 
+ c2
~
z
 
3
 
+ z0c
~
z
 
-1
 
+ z0c1
~
z
 
+ z0c2
~
z
 
2
 
+a.
(77)
Ее свободный член есть -2c2 + 6c2c2 + c. Следовательно, c2 = -c/4 = 1/4.
Вычислим в сумме (77) коэффициент при [(z)\tilde]1. Получаем
b3 = 6cc12+z0c1 = c1(6cc1+z0) = c1(-c2z0+z0)=0.
Следовательно, коэффициент c3 произволен. Это случай (a) замечания 1.5. Итак, ребру G(1)1 соответствуют степенные разложения решений
w = c
~
z
 
-1
 
+ c1
~
z
 
+ c2
~
z
 
2
 
+ c3
~
z
 
3
 
+ Ґ
е
k=4 
ck
~
z
 
k
 
,
(78)
где c = ±1, c1 = ±z0/6, c2 = 1/4, c3 - произвольно, а все остальные ck однозначно определяются. Они образуют два однопараметрических семейства G(1)11 и G(1)12.
3.4. Решения, соответствующие ребру G(1)4. Ребру G(1)4 соответствует укороченное уравнение
^
f
 
(1)
4 
def
=
 
-w" + a = 0
(79)
с нормальным конусом U(1)4 = {P: p2 = 2p1 < 0}, т.е. r=2 и w = -1. Будем искать решение w=c[(z)\tilde]2. Для постоянной c получим следующее значение: c = a/2 и c 0, т.к. по предположению, a 0. Т.к. w = -1, то z ®0, и разложения вида (33) у нас будут по возрастающим степеням [(z)\tilde]s с s > 2. Множество
K = {s > 2}
(80)
это конус задачи для ребра G(1)4. Первая вариация d[^(f)] / dw, оператор L(z) и функция n(k) здесь имеют вид (39), (40) и (41) соответственно. Поскольку корни k1=0 и k2 = 1 уравнения n(k)=0 не лежат в конусе задачи (80), то критических чисел нет. Т.к. нет критических чисел, то по теореме 1.4 все показатели s лежат в множестве
K = {4,5,6,7...}.
(81)
Найдем второе приближение. У нас m = 0, mў = -2, тогда ближайшей точкой к r = 2 будет s0 = 4. Т.к. æ = 0, то вычислим коэффициент c4 по формуле (26). Здесь [^([^(f  )])] = z0w,     b4 = z0c, по (41) n(4)=-12, откуда по (26) имеем c4 = z0c/12 = az0/24. Таким образом, ребру G(1)4 соответствует степенное разложение
w = c
~
z
 
2
 
+ c4
~
z
 
4
 
+ Ґ
е
l=5 
cl
~
z
 
l
 
,
(82)
где c = a/2, c4 = (az0)/24 и все остальные cs однозначно определяются. Это семейство обозначим G(1)4.
3.5. Сводка результатов и их обсуждение. Были найдены следующие разложения решений уравнения (72):
  1. Двупараметрическое (по c0 и c1) семейство G(0)11 степенных разложений вида (44).

  2. Однопараметрическое (по c1) семейство G(0)12 степенных разложений вида (46), которое включается в предыдущее семейство при c0 = 0.

  3. Два однопараметрических (по c3) семейства G(1)11 и G(1)12 решений (78).

  4. Семейство G(1)4, состоящее из одного решения (82).

По теореме 1.7 все они сходятся при достаточно малых |[(z)\tilde]|. Существование разложений 1 и 2 семейств G(0)11 и G(0)12 следует из теоремы Коши. Разложения 3 семейств G(1)11 и G(1)12 найдены в [12]. Разложение 4, т.е. семейство G(1)4, найдено, по-видимому, впервые.


§ 0. Разложения решений при a=0
При a=0 второе уравнение Пенлеве есть
f(z,w) def
=
 
-w" + 2w3 + zw = 0.
(83)
4.1. Объекты степенной геометрии. Носитель S(f) уравнения (83) состоит из трех точек:
Q1 = (-2,1),     Q2 = (0,3),     Q3 = (1,1).
(84)
Их выпуклая оболочка G это треугольник с вершинами G(0)j = Qj и ребрами G(1)j = {Qj,Qj+1}, где j=1,2,3,Q4=Q1 (см. рис. 4).
Внешние нормали Nj(1) к ребрам G(1)j cуть N1(1) = (-1,1),    N2(1)=(1,1/2),    N3(1) = (0,1). Нормальные конусы ребер суть
U(1)1 = {P: p1 = 0, p1 > -p2} = {lN1},
U(1)2 = {P: p1 = -p2, p1 < 2p2 } = {lN2},
U(1)3 = {P: p2=p1/2 > 0} = {lN3}, l > 0.
(85)
Нормальные конусы Uj(0) вершин Gj(0) - это углы между нормальными конусами прилегающих к вершине ребер (см. рис. 5):

U(0)1 = {P: 2p2 < p1,    p1 > 0 },
U(0)2 = {P: p1 > -p2,    p1 < 2p2 },
U(0)3 = {P: p1 < 0,    p1 < -p2 }.
(86)
Теорема 4.1. Сдвинутый носитель Sў(f) лежит в решетке Z3 с базисом B4 = (3,0), B5 = (-1,2).


4.2. Решения, соответствующие вершинам. Рассмотрим укороченное уравнение, соответствующее вершине G(0)1:
^
f
 
(0)
1 
def
=
 
- w" = 0.
(87)
Этот случай аналогичен тому, который описан в пункте 2.2.
Следовательно, вершине G(0)1 соответствуют два семейства решений: семейство G(0)11 (формула (44)) и семейство G(0)12 (формула (46)).
Укороченные уравнения, соответствующие вершинам G(0)2 и G(0)3:
^
f
 
(0)
2 
def
=
 
2w3 = 0,    
^
f
 
(0)
3 
def
=
 
zw = 0
(88)
являются алгебраическими. По замечанию 1.1 они не имеют решений вида w = cr zr,     cr 0. Поэтому укорочения, состоящие из одного алгебраического монома, можно не рассматривать.
4.3. Решения, соответствующие ребру G(1)1. Ребру G(1)1 соответствует укороченное уравнение
^
f
 
(1)
1 
def
=
 
-w" + 2w3 = 0.
(89)
Этот случай аналогичен тому, который был описан в п. 2.3. Таким образом, ребру G(1)1 соответствуют степенные разложения (55):
w = cz-1 + c2z2 + c3z3 + Ґ
е
s=5 
cszs,
(90)
где коэффициенты c = ±1, c2 = -c/6, c3 - произвольное, а остальные cs однозначно определяются. Они образуют два однопараметрических семействa G(1)11 и G(1)12. При c3 = 0 получается разложение (56).
4.4. Решения, соответствующие ребру G(1)2. Ему соответствует укороченное уравнение
^
f
 
(1)
2 
def
=
 
zw + w3 = 0.
(91)
Уравнение (91) имеет два решения w = cz1/2, где c2 = 1/2,     c = ±iЦ{1/2}, которые нам подходят. Этот случай аналогичен тому, который описан в п. 2.4. Исследовав данное ребро по той же схеме, получим, что ребру G(1)2 соответствует степенное разложение
w = cz1/2+
е
s 
cszs,
(92)
где все показатели степени s О K:
K = {-3l-5/2, целые     l і 0},
(93)
где l - целое и c = ±i/Ц2, все cs однозначно определены. Найдем коэффициент второго приближения. Сначала найдем значения m и mў. В нашем случае, т.к. P = w(1,r) = (1, 1/2), m = 3/2, а mў = -3/2. Тогда s0 = -5/2 и по формуле (26) находим, что c-5/2 = -c/16. Итак, ребру G(1)2 соответствуют степенные разложения
w = cz1/2 + Ґ
е
l 
clz-3l-5/2,
(94)
причем c = ±Ц2/2,     c-5/2 = -c/16. Они образуют семейства G(1)21 и G(1)22.
4.5. Решения, соответствующие ребру G(1)3. Ему соответствует укороченное уравнение
^
f
 
(1)
3 
def
=
 
-w" + zw = 0.
(95)
Уравнение (95) имеет решение w є 0. Но оно не имеет решений вида w = crzr,     cr 0. Ребро G(1)3 горизонтально, поэтому для нахождения решений укороченного уравнения (95) воспользуемся методами п.п. 1.9 и 1.10. Согласно п.1.10 сделаем логарифмическое преобразование вида (30), т.е. y = d lnw/dz. Тогда wў = yw, w" = yўw+y2w, и уравнение (95) переходит в уравнение
w(-yў- y2 + z) = 0.
(96)
После сокращения на w получаем уравнение
h(z,y) def
=
 
-yў- y2 + z = 0.
(97)
Для уравнения (97) строим его носитель [(S)\tilde] = S(h), многоугольник [(G)\tilde] = G(h), его грани [(G)\tilde](d)j, их нормальные конусы [(U)\tilde](d)j. Носитель [(S)\tilde] (см. рис. 6) состоит из трех точек [(Q)\tilde]1 = (-1,1),     [(Q)\tilde]2 = (0,2),    [(Q)\tilde]3 = (1,0). Многоугольник [( G)\tilde] = G(h) это треугольник с вершинами [(G)\tilde](0)j = [(Q)\tilde]j (см. рис. 6). Ребра (см. рис. 6) [(G)\tilde](1)1 = {[(Q)\tilde]1, [(Q)\tilde]2},    [(G)\tilde](1)2 = {[(Q)\tilde]2,[(Q)\tilde]3},     [(G)\tilde](1)3 = {[(Q)\tilde]3, [(Q)\tilde]1}.
Внешние нормали [(N)\tilde]j к ребрам [(G)\tilde](1)j суть [(N)\tilde]1 = (-1,1),     [(N)\tilde]2 = (2,1),     [(N)\tilde]3 = (-1,-2).
Нормальные конусы [(U)\tilde](d)j (см. рис. 7) суть
~
U
 
(1)
1 
= {p1 = -p2 < 0} = {l
~
N
 

1 
},
~
U
 
(1)
2 
= {p2 = 2p1 > 0} = {l
~
N
 

2 
},
~
U
 
(1)
3 
= {p1 = 2p2 < 0} = {l
~
N
 

3 
},l > 0.
(98)

~
U
 
(0)
1 
= {p1 < -p2, p1 < 2p2},
~
U
 
(0)
2 
= {p2 > -p1,p2 > 2p1},
~
U
 
(0)
3 
= {p1 > 2p2, p2 < 2p1}.
(99)
Теорема 4.2. Носитель S(h) лежит в решетке с Z6 базисом B6 = (3,-1), B7 = (0,1).
Согласно п. 1.10 для уравнения (97) получаем конус задачи [(p)\tilde]1 + [(p)\tilde]2 > 0. На рис.7 показана его граница. С этим конусом задачи пересекаются только нормальные конуса [(U)\tilde](0)2, [(U)\tilde](0)3 и [(U)\tilde](1)2. Поэтому в дальнейшем достаточно рассмотреть только решения, соответствующие граням [(G)\tilde](0)2, [(G)\tilde](0)3 и [(G)\tilde](1)2. Вершинам [(G)\tilde](0)2 и [(G)\tilde](0)3 соответствуют алгебраические уравнения -y2 = 0 и z=0 соответственно. По замечанию 1.1 они не дают подходящих решений. Ребру [(G)\tilde](1)2 соответствует алгебраическое укороченное уравнение
-y2 + z=0.
(100)
Оно имеет два подходящих решения y = ±z1/2, т.е. r=1/2 и w = 1. Согласно замечанию 1.2 эти укороченные решения не имеют критических чисел (æ = 0) и для них
n(k) = 2.
(101)
Здесь вектор сдвинутого носителя укороченного уравнения есть (1/2, -1). Вместе с векторами теоремы 4.2 он порождает решетку Z6 с базисом (3/2, 0), (1/2, -1). Она пересекается с прямой q2=-1 по точкам q1=1/2 + 3l/2, где l - целое. Поскольку конус задачи для искомого разложения есть s < 1/2, то по замечанию 1.3 множество
K = {s = 1/2 - 3l/2, целые l > 0}.
(102)
По теореме 1.3 уравнение (97) имеет решение в виде ряда
y = ±z1/2 + z1/2 Ґ
е
l=1 
clz-3l/2,
(103)
где все коэффициенты cl определяются однозначно при фиксированном знаке y корня из z. Вычислим коэффициент c1 в разложении (103). Здесь [^([^(h  )])] = -yў. Поэтому b1 = 1/2. Согласно (26) и (101) имеем c1 = -1/4. Пусть h = lnw, тогда hў = y и согласно (103)
h = у
х
y dz = ±  2

3
z3/2+c0 -  1

4
lnz+ Ґ
е
l=2 
 2cl

3(1-l)
z3(1-l)/2,
(104)
где c0 - произвольная постоянная. Для решений w(z) укороченного уравнения (97) это дает разложения
w =
~
c

z1/4
exp[±  2

3
z3/2+ Ґ
е
l=2 
 2cl

3(1-l)
z3(1-l)/2]
(105)
при z® Ґ. Согласно нормальному конусу [(U)\tilde](1)3 имеем w ® 0 при z® Ґ. Следовательно, в квадратной скобке в (105) знак берется так, чтобы
±Rez3/2 < 0.
(106)
По теореме 1.8 решение укороченного уравнения (95) является асимптотикой решения исходного уравнения, если w® 0, т.е. при выполнении условия (106) для z. Таким образом, получили два однопараметрических семейства G(1)31 и G(1)32 сильных асимптотик.
4.6. Сводка результатов и их обсуждение. Были найдены следующие разложения решений уравнения (83).
  1. Двупараметрическое (по c0 и c1) семейство G(0)11 степенных разложений (44).

  2. Однопараметрическое (по c) семейство G(0)12 степенных разложений (46).

  3. Два однопараметрических (по c3) семейства G(1)21 и G(1)22 степенных разложений (90). При c3 = 0 получаются разложения (56).

  4. Два семейства G(1)31 и G(1)32 степенных разложений (92).

  5. Два однопараметрических (по [(c)\tilde]) семейства G(1)31 и G(1)32 асимптотик (105).

По теореме 1.7 разложения 1 - 3 сходятся при достаточно малых |z|. Разложения 4 и 5 могут расходиться при любом z. Интересно, что носители разложений 2 - 4 отличаются от носителей соответствующих разложений в случае a 0. Это, по-видимому, новый результат. Укороченное уравнение (95) это известное уравнение Эйри [15] и разложения (105) для его решений известны [16].


Литература
  1. Брюно А.Д. Локальный метод нелинейного анализа дифферециальных уравнений. М.: Наука, 1979. 256 c.

  2. Брюно А.Д. Степенная геометрия в алгебраических и дифференциальных уравнениях. М.: Физматлит, 1998. 288 c.

  3. Брюно А.Д. Автомодельные решения и степенная геометрия // Успехи математических наук, 2000, т. 55, N 1, с. 3-44.

  4. Брюно А.Д. Степенные разложения решений одного алгебраического или дифференциального уравнения // Доклады АН, 2001, т. 380, N 2, с. 155-159.

  5. Брюно А.Д. Степенные разложения решений одного алгебраического или дифференциального уравнения. Препринт N 63. М.: ИПМ им. М.В. Келдыша, 2000. 22 с.

  6. Bruno A.D. Power geometry as a new calculus // Proceeding of ISAAC 2001. Amsterdam: Kluwer, 2003.

  7. Брюно А.Д. Асимптотики и разложения решений обыкновенного дифференциального уравнения. Препринт N 9. М.: ИПМ им. М.В. Келдыша, 2003. 39 с.

  8. Тихомиров В.М. Фреше производная // Математическая Энциклопедия. М.: Советская Энциклопедия, 1985, т. 5, с. 666.

  9. Розов Н.Х. Пенлеве уравнение // Математическая Энциклопедия. М.: Советская энциклопедия, 1984, т. 4, с. 233-234.

  10. Gromak V.I., Laine I., Shimomura S., Painleve Differential Equations in the Complex Plane, Walter de Gruyter. Berlin, New York, 2002. 303 p.

  11. Голубев В.В. Лекции по аналитической теории дифференциальных уравнений. М.-Л.: ГИТТЛ, 1941, 400 с. 2-е изд., М.-Л.: ГИТТЛ, 1950.

  12. Голубев В.В. К теории уравнений Пенлеве. // Матем. сборник, 1912, т. 28, N 2, с. 323-349.

  13. Козлов В.В., Фурта С.Д. Асимптотики решений сильно нелинейных систем дифференциальных уравнений. М.: МГУ, 1996. 244 с.

  14. Яблонский А.И. Асимтотическое разложение правильных решений некоторых классов дифференциальных уравнений // ДАН БССР, 1964, т. 8, N 2, с. 77-80.

  15. Федорюк М.В. Эйри уравнение // Математическая энциклопедия, т. 5, М.: Советская энциклопедия, 1985, с. 939.

  16. Федорюк М.В. Эйри функция // Математическая энциклопедия, т. 5, М.: Советская энциклопедия, 1985, с. 939 - 941.




File translated from TEX by TTH, version 3.40.
On 02 Mar 2005, 19:41.