Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

0387986995Basic TheoryC

.pdf
Скачиваний:
23
Добавлен:
10.06.2015
Размер:
11.75 Mб
Скачать

428

XIII. SINGULARITIES OF THE SECOND KIND

so that A' is a diagonal matrix with n distinct diagonal entries. Then, a transformation similar to (T') can be found so that (E") is changed to (E') with suitable constants ao, al, ... , an_ I .

XIII-6. The Hukuhara-Turrittin theorem

In this section, we explain a theorem due to M. Hukuhara and H. L. TLrrittin that clearly shows the structure of solutions of a system of linear differential equations of the form

(XIII.6.1)

where the entries of the n x n matrix A are in K (cf. §XIII-5). In order to state this theorem, we must introduce a field extension f- of K. To define L, we first set

+a

E a.nxm1" : a.., E C and

M E Z

I

,

M=M

 

 

+a

where 7L is the set of all integers. For any element a = F a,nxn'/° of K,,, we

M=M

+00

define xji by x da = F \ v) amx'nl '. Then, K, is a differential field. The field

M=M

+oo

L is given by L = U K which is also a differential field containing K as a subfield.

L=1

Furthermore, L is algebraically closed. The Hukuhara-T rrittin theorem is given as follows.

Theorem XIII-6-1 ({Huk4j and [Tu1j). There exists a transformation

(XIII.6.2)

y' = UI

such that

(i)the entries of the matrix U are in L and det U ;j-1 0,

(ii)transformation (XIII.6.2) changes system (X111.6.1) to

(XIII.6.3)

xd = Bz',

where B is an n x n matrix in the Jordan canonical form

B=diag[Bl,B2,...,Bpj, Bj=diag[BJi,Bj2i...,B,,n,],

(XIII.6.4)

Bjk =AI In,,. + Jn,,,.

6. THE HUKUHARA-TURRITTIN THEOREM

429

Here, I,,,,

is the njk x njk identity matrix, J.,,, is an n 1k x n3k nilpotent

matrix of the form

 

 

 

0

1

0

 

0

0

1

(XIII.6.5)

Jn",

 

 

 

0

0

0

 

0

0

0

and the Aj are polynomials in xfor some positive integer s, i.e.,

 

d,

 

 

 

(XIII.6.6)

Aj = > A

x-'18

where all E C

(j = 1, 2, ... , p)

 

r=O

 

 

 

and

 

 

 

 

(XII1.6.7)

0 if dJ > 0

 

A., - A, are not integers if i 0 j.

A,d,

and

Proof.

Without loss of generality, assume that the matrix A of system (XIII.6.1) has the form

0

1

0

0

0

 

0

0

1

0 ...

0

 

A

 

 

 

 

 

0

0

0

0 ...

1

 

00

a1

a2

a3 "'

an-1

 

 

 

 

 

an- 1

In.

where ak E 1C (k = 0, 1, ... , n - 1) (cf. Theorem XIII-5-4). Set E= A - n

Then,

 

 

 

 

 

(XIIi.6.8)

trace [E] = 0.

 

 

Consider the system

 

 

 

 

 

(XIII.6.9)

X E = E.

 

 

Case 1. If there exists an n x n matrix S with the entries in X such that det S 0 0 and the transformation w = Sii changes (XIII.6.9) to x dii = A(x)u", where the

entries of A are in C[[x]], then there exists another n x n matrix S with the entries in X such that det S 96 0 and the transformation

(XIII.6.10)

w = Suu

changes (XIII.6.9) to xdig = Aoii, where the entries of the matrix A0 are in C.

F irthermore, any two distinct eigenvalues of A0 do not differ by an integer (cf.

Theorem V-5-4). Hence, in this case, system (XIII.6.1) is changed by transformation (XIII.6.10) to

dil

xaj = [an 1 In+Ao]u.

This proved Theorem XIII-6-1 in this case.

430 XIII. SINGULARITIES OF THE SECOND KIND

Case 2. Assume that there is no n x n matrix S with entries in X such that det.S 0 and the transformation w = Su' changes (XM.6.9) to xdu = A(x)u, where the entries of A are in C([x]]. Since

 

1

0

0

 

0

 

- lan_1

1

0

 

0

 

n

 

 

 

 

E =

 

 

 

-lan_1

 

0

0

0

 

1

 

 

 

 

n

 

ap

al

a2

a3

an-2

an-1 -- -an-1,

a cyclic vector can be found for system (XIII.6.9) by using the matrix W defined

by

 

w11

win

 

 

 

 

with

 

writ

Wnn

 

 

 

 

 

[wll ... Win] = (10 ... 01,

 

 

[w,,1 ... wJn] = V1-1[1 0

01,

where V([c1

cn]) = x[c1

[c1

- - cn]E. The matrix W is lower-

triangular and the diagonal entries are {1, ... ,1}, i.e.,

(XIII.6.11)

W =

1

0

0

0

 

1

0

0

 

 

 

 

1

 

 

If (XIII.6.9) is changed by the transformation v" = Wtv", then

 

(XIU.6.12)

X!LV _ x

+ WEJ W-Y.

 

 

It follows from (XIII.6.8) and (XIII.6.11) that

 

 

 

 

( XIII . 6. 13)

tracel x

+ WE]W- l

= 0.

 

 

Also,

 

 

 

 

 

 

 

 

0

1

0

0

0

 

 

0

0

1

0

0

x

+ WE] W-1 = V[ W]W-1

0

0

0

l

11

 

0

 

 

W i31

/32

2 2

A

 

 

 

An-1

6. THE HUKUHARA-TURRITTIN THEOREM

431

where 1k E )C (k = 0, 1,... , n - 1). In particular, from (XIII.6.13), it follows that On-1 = 0. Under our assumption, not all Qt are in C[[x]]. Set 9 = (t: at V C[[x]]}

+00

 

and set Also, /3t = x-u" E

Qtmxm (t E 3), where, for each t, the quantity µt is a

m=0

 

+oo

 

positive integer, 1: /3tmxm E C[[xJ] and #to 54 0. Set m=0

k=max(n µt t :IE91.

J

Then, ut < k(n - t) for every t E J and µt = k(n - t) for some t E J. This implies

that

(I)

and

(11)

Of = L. (t=0,1,... n- 1)

Qt,mxm

m>-k(n-t)

Qt = x-k(n-t) (Ct + xqt)

for some t such that k(n - t) is a positive integer, ct is a nonzero number in C, qt E C[[x]], and Qt,m E C. We may assume without any loss of generality that

k = h for some positive integers h and q.

9

Let us change system (XIII.6.12) by the transformation

v" = diag [1, x-k, ... , x-(n-1)kI u".

Then,

(XIII.6.14)

where

 

0

1

0

0

 

0

0

1

0

(XIII.6.15) F =

 

 

+ kxkdiag [0,1, ... , n - 1]

 

0

0

0

0

 

70

71

72

73

and

= :

 

 

 

7t = xk(n-t)o,

 

$,,mxm+k(n-t) E C[ [x'1911 (0 < t < n - 1).

 

 

m>-k(n-t)

In particular, 7n-1 = 0.

432

XIII. SINGULARITIES OF THE SECOND KIND

+oo

Setting F = E xm/9Fm, where the entries of F,,, are in C, we obtain m=o

0 1 0 0 0 01

Fo=

C1 C2 C3 ... Cn_2 0

CO

where the constants co, cl ... , cn_2 are not all zero. This implies that the matrix

F0 must have at least two distinct eigenvalues. Hence, there exists an n x n matrix

T such that

(1) T = XM19Tm, where the entries of the matrices Tm are in C and To is

m=o

invertible,

(2)the transformation

(XIII.6.16)

y = Ti

changes system (XIII.6.14) to a system xd = x-kGxi with a matrix G in

a block-diagonal form G = [ 0 0 , where Gl and G2 are respectively

G2 J

n1 x nl and n2 x n2 matrices with entries in C[[x1/91and that nl + n2 = n and n., > 0 (j = 1, 2) (cf. §XIII-5). Therefore, the proof of Theorem XIII-6-1 can be completed recursively on n. 0

Observation XIII-6-2. In order to find a fundamental matrix solution of (XIII.6.1), let us construct a fundamental matrix solution of (XIII.6.3) in the following way:

Step 1. For each (j, k), set

4'1k = xA,0 eXp[AJ (x)] exp[(log x)Jf, ],

where

if dj = 0,

-t/s if d1>0.

Step 2. For each j, set

Dj = diag ['j1,

xAJ0 exp[A,(x)] exp[(logx)Jj],

where J. = diag [J,,,,, J,,72 , ... J,mJ ] .

6. THE HUKUHARA-TURRITTIN THEOREM

433

Step 3. Set

A = ding [Al(x)In A2(x)I,,,, ... ,

C = diag [,\101., + J1, 1\20I., + J2, ... ,

Then,

(XIII.6.17) 4) = diag [451, 4i2, ... , 4ip] = xC exp[A]

is a fundamental matrix solution of (XIII.6.3), where

xC = diag [xA1OxJ, xl\2OxJ2' ... 'X aPOxJP] ,

xJ, = exp[(logx)JjJ.

The matrix

 

 

(XIII.6.18)

U4i = Uxc exp[A]

 

is a formal fundamental matrix solution of system (XIII.6.1), where U is the matrix of transformation (XIII.6.2) of Theorem XIII-6-1. The two matrices exp[A] and xC commute.

Observation XIII-6-3. Theorem XIII-6-1 is given totally in terms of formal power series. However, even if the matrix A(x) of system (XIII.6.1) is given analytically, the entries of U of transformation (XIII.6.2) are, in general, formal power series in xl1", since the entries of the matrix T(x) of transformation (XIII.6.16) are formal power series in general. Transformation (XIII.6.16) changes system (XIII.6.14) to a block-diagonal form. Therefore, in a situation to which Theorem XIII-4-1 applies, transformation (XIII.6.2) can be justified analytically. The following theorem gives such a result.

Theorem XIII-6-4. Assume that the entries of an n x it matrix A(x) are holomorphic in a sector So = {x E C 0 < lxJ < ro, I argxi < ao} and admit asymptotic expansions in powers of x as x 0 in So, where ro and ao are positive numbers.

Assume also that d is a positive integer and y" E C'. Let S be a subsector of So whose opening is sufficiently small. Then, Theorem XIII-6-1 applies to the system

(XIII.6.19)

xd+1dy = A(x)f

 

dx

with transformation (XIII.6.2) such that the entries of the matrix U of (XIII.6.2) are holomorphic in S and each of them is in a form x,00(x) where p is a rational number and O(x) admits an asymptotic expansion in powers of xl"° as x 0 in S, where s is a positive integer.

Observation XIII-6-5. In the case when the entries of the matrix A(x) on the right-hand side of (XIII.6.19) are in C[[x]l and A(0) has n distinct eigenvalues, the matrix A(x) also has n distinct eigenvalues A1(x), A2(x), ... , which are in

C((xJJ. Furthermore, the corresponding eigenvectors p"1(x), p"2(x), ... ,15n(x) can be constructed in such a way that their entries are in Chill and that p""1(0),7"2(0), ... ,

434

XIII. SINGULARITIES OF THE SECOND KIND

p",,(0) are n eigenvectors of A(O). Denote by P(x) the n x n matrix whose column vectors are p1(x), p2(x), ... , p,+(x). Then, detP(O) 36 0 and P(x)-lA(x)P(x) = diag[A1(x),A2(x),... ,An(x)]. This implies that the transformation y = P(x)ii changes system (XIII.6.19) to

(XIII.6.20) xd+1 dx {diagiAi(x)A2(x).... , An(x)] - xd+1P(x)-1 d ( )

It is easy to construct another n x n matrix Q(x) so that (a) the entries of Q(x) are in C[(x]], (b) Q(0) = 4,, and (c) the transformation i = Q(x)v changes system

(XIII.6.20) to

(XIII.6.21)

xd+1 dv

= diag [ft1(x), µ2(x), ... , µn(x)] v,

dx

 

 

where µ1(x), µ2(x), ... , i (x) are polynomials in x of degree at most d such that A, (x) = it) (x) + O(xd+1) ( j = 1, 2, ... , n). Therefore, in this case, the entries of the matrix U of transformation (XIII.6.2) are in 1C.

Observation XIII-6-6. Assume that the entries of A(x) of (XIII.6.19) are in

C([x]]. Assume also that A(O) is invertible. Then, upon applying Theorem XIII-6-1

to system (XIII.6.19), we obtain

following theorem.

d1 = d for all j. Using this fact, we can prove the

s

Theorem XIII-6-7. Let Q;i1(x) and Qi2(x) be two solutions of a system

(XIII.6.22)

xd+1 dy

= A(x)yf + x f (x),

 

dx

 

where d is a positive integer, the entries of the n x n matrix A(x) and the C"-valued function 1 *(x) are holomorphic in a neighborhood of x = 0, and A(O) is invertible.

Assume that for each j = 1, 2, the solution ¢,(x) admits an asymptotic expansion in powers of x as x - 0 in a sector S3 = {x E C : lxl < ro, aj < arg x < b3 }, where ro is a positive number, while aJ and bi are neat numbers. Suppose also that S1 n S2 0. Then, there exist positive numbers K and A and a closed subsector S =

{x : lxl < R,a < argx < b} of Sl nS2 such that

K exp[-Alxl -d]

in S.

Proof.

Since the matrix A(O) is invertible, the asymptotic expansions of1(x) and ¢2(x) are identical. Set 1 (x) = ¢1(x) -$2(x). Then, the C"-valued function >G(x) satisfies system (XIII.6.19) in S, n$2 and ii(x) ^_- 6 as x 0 in S1nS2. By virtue of Theorem

XIII-6-4, a constant vector 66 E C" can be found so that r%i(x) = U4D(x)6, where lb(x) is given by (XIII.6.17). Now, using Observation XIII-6-6, we can complete the proof of Theorem XIII-6-7.

6. THE HUKUHARA-TURRITTIN THEOREM

435

Observation XIII-6-8. The matrix A = diag [AI In,, A2In3, ... , API,y] on the right-hand side of (XIII.6.18) is unique in the following sense. Assume that an-

other formal fundamental matrix solution Uxcexp[A] of system (XIII.6.1) is con- structed with three matrices U, C, and A similar to U, C, and A. Since the

matrices Uxc exp[A] and UxO exp[A] are two formal fundamental matrices of sys- tem (XHI.6.1), there exists a constant n x n matrix r E GL(n, C) such that

Uxc exp[A] = U5C exp[A]I' (cf. Remark IV-2-7(1)). Hence, exp[A]T exp[-A] = x-CU-IUxC. Using the fact that r is invertible, it can be easily shown that A = A if the diagonal entries of A are arranged suitably. For more information concerning the uniqueness of the Jordan form (XIII.6.3) and transformation (XIII.6.2), see, for example, [BJL], [Ju], and [Leve].

Observation XIII-6-9. The quantities A.,(x) are polynomials in x1l'. Set w =

2a[!]

and x

1/a =

wx

1/a

Then, i = x. Therefore, if z

1/e

C

exp

 

 

 

in Ur exp[A] is

replaced by zI/', then another formal fundamental matrix of (XIII.6.1) is obtained.

This implies that the two sets {Aj (i) : j = 1, 2,... , p} and (A.,(x) : j = 1, 2,... , p) are identical by virtue of Observation XIII-6-8.

Observation XIII-6-10. A power series p(x) in x1/' can be written in a form a-1

Ax) = Exh1'gh(x), where ql(x) E C[[x]] (j = 0, 1, ... , s - 1). Using this fact and

h=0

Observation XIII-6-9, we can derive the following result from Theorem XIII-6-1.

Theorem XIII-6-11. There exist an integer q and an n x n matnx T(x) whose entries are in C[[x]] such that

(a)det T (x) 96 0 as a formal power series in x,

(b)the transformation y" = T(x)t changes system (XII1.6.1) to x dil = E(x)iZ with an it x it matrix E(x) such that entries of x9E(x) are polynomials in x.

The main issue here is to construct, starting from Theorem XIII-6-1, a formal transformation whose matrix does not involve any fractional powers of x in such a way that the given system is reduced to another system with a matrix as simple as possible. A proof of Theorem XIII-6-11 is found in [BJL].

Changing the independent variable x by x-1, we can apply Theorem XIII-6-1 to singularities at x = oo. The following example illustrates such a case.

Example XIII-6-12. A system of the form

 

P(x) 0]b,

y = lyd,

where P(x) = xm+

ahxin is a positive odd integer, and the ah are complex

h=1

numbers, has a formal fundamental matrix solution of the form

x

F(x)

1

0

f 1

1

e

0

 

 

 

 

 

0

x-1/2 ]

l I

-111 0

eE(t,a) I

set th =
q

436 XIII. SINGULARITIES OF THE SECOND KIND

where

m

1/2

+00 bk(a)

 

ak

 

 

+ E xk

 

+ E

xk

 

k=1

 

 

k=1

1<h<(m+2)/2i

 

E(x, a) =

2

2) x(m+2)/2 +

(a)x(m+2-2h)/2,

(m +

(.m + 2 - 2h) bh

+00

and F(x) = xq> x-h Fr, with an integer q and 2x2 constant matrices Fh such that

h--O

+00

det h=O x-h Fh ,-6 0 as a formal power series in x-l. For details of construction, see [HsSj and [Sil3j.

XIII-7. An n-th-order linear differential equation at a singular point of the second kind

Let us look at the formal fundamental matrix solution (XIII.6.18) of system

(XIII.6.1). First, notice that if we set k = max { a :j = 1,2,... pthen k is the

 

s

JJJ

 

order of singularity of system (XIII.6.1) at x = 0 (cf. Definition V-7-8). Let

 

(XIII.7.1)

0<k1 <k2 <... <k9_1 <kq

 

 

 

 

d1 d2

 

be the distinct values among p non-negative rational numbers

s s , ... , .s

Also,

n, (h = 1, 2,... , q). It is easy to see that lh > 0 (h = 1, 2,... , q)

d /s=kh

P

and Elh = En, = n.

h=1 7=1

Observation XIII-7-1. System (XIII.6.1) has th linearly independent formal so- lutions of the form

(XIII.7.2)

'ih,v(x) = x'1h - exp[Qh,v(x)Wh,,,(x) (v = 1, ... ,;h = 1, ... , q),

where -yh,,, E C, Qh,,,(x) is either equal to 0 or a polynomial in x-1/' of the form

Qh.v(x) =

I2h,vx-kh(I +O(x1/'))

(11h,v E C and I2h,, # 0),

and the entries of lh,,,(x) are polynomials in logx with coefficients in

Define q + 1 points (Xh, Yh) (h = 0,1, ... , q) recursively by

(Xo,Y0)

(010),

 

(Xh, Yh)== (Xh-1 + eh, Yh-t + kheh)

(h = 1, 2,... , q).

Let us denote by N the polygon whose vertices are q + 1 points (Xh, Yh) (h =

0,1, ... , q). The polygon N has q distinct slopes kh given in (XIII.7.1).

7. AN N-TH-ORDER LINEAR DIFFERENTIAL EQUATION

437

Definition XIII-7-2. The polygon N is called the Newton polygon of system

(XIII.6.1) at x = 0.

Observation XIII-7-3. In §XIII-5, it was shown that system (XIII.6.1) is equiv- alent to an n-th-order linear differential equation

(XIII.7.3)

where b = x operator

(XIII.7.4)

n-1

anb" r) + E atbtrl = 0,

t=o

, at E C((xj), and an # 0 (cf. Theorem XIII 5 4). For the differential

n-1

C(>,J = anb"r) + j:atdtrt. t=o

the Newton polygon N(C) is defined in the following way.

If a power series a = E c,,x'n E C((xj] is not 0, we set v(a) = min{m : cn

m=0

0}. If a = 0, set v(0) = +oo. For operator (XIII.7.4), consider n + 1 points

(Q, v(at)) (t = 0, 1, ... , n) on an (X, Y)-plane. Set

Pt = ((X, Y) : 0 < X < e, Y > v(at)}, and P = UP,.

t=o

Definition XIII-7-4. The boundary curve C of the smallest convex set containing P is called the Newton polygon of the operator C at x = 0.

Denote by N(C) the Newton polygon of C at x = 0.

Definition XIII-7-5. Two Newton polygons are said to be identical if the two polygons become the same by moving one or the other upward in the direction of the Y-axis.

Now, we prove the following theorem.

Theorem XIII-7-6. If system (XIII.6.1) and differential equation (X111.7.3) are equivalent in the sense of Theorem X111-5-4. then the two Newton polygons N and

N(C) are identical.

Proof.

The proof of this theorem will be completed if the following three statements are verified:

(a) If./V(C) has only one nonvertical side with slope k, then differential equation

(XIII.7.3) is equivalent to a system xk+1

= A(x)g with a matrix A(x)

whose entries are power series in x11", and A(0) is invertible if k > 0, where s is a positive integer such that sk is an integer.