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

0387986995Basic TheoryC

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

118

V. SINGULARITIES OF THE FIRST KIND

Step 6. Setting

 

 

#(x) = xN+1 ("liM xJ`'+1) = ku> nik(x),

it is easy to show that n'(x) satisfies integral equation (V.2.12). It is also evident that i7(x) is analytic for IxI < 6. Thus, the proof of Theorem V-2-5 is completed. 0 Now, finally, by using the argument given in Observations V-2-2 and V-2-3, we

obtain the following theorem.

Theorem V-2-7. Every formal solution E xC[[zjJ" of system (V.2.1) is conver- gent, i.e., E xC{x}".

Remark V-2-8. In general, (V.2.1) may not have any formal solutions. However,

Theorem V-2-7 states that if (V.2.1) has formal solutions, then every formal solution is convergent.

Remark V-2-9. If yi(x) = xA f (x) E XAC[[xfj" is a formal solution of a linear system X!LY = A(x)y, the formal power series f is a formal solution of xd _

[A(x) - A!, }i. Therefore, f E C{x}" by virtue of Theorem V-2-7. This proves convergence of the formal solution constructed in Theorem V-1-2.

V-3. The S-N decomposition of a matrix of infinite order

In §V-4, we shall define the S-N decomposition of a linear differential operator.

As a linear differential operator will be represented by a lower block-triangular matrix of infinite order, we derive, in this section, the S-N decomposition of such

a matrix

0

0 ... ... ...

All

A21

A22

0 ... ... ...

Amt Am2

...

Am,n

0 . .

where for each (j, k), the quantity AJk is an n, x nk constant matrix. Set

Al = All,

0

 

 

 

All

O

..

0

A21

A22

O

..

0

Am =

 

 

 

(m > 2).

Am l

Am2

...

...

Am,

Form > 2, we write Am -

Am-!

0

Set Nm = m nt. Then, A. is an

B.

 

[

Amm J

l-1

 

 

 

 

N", x Nm matrix, while Bm is an nm x N,"_ 1 matrix. Also, let Amm = Smm +Nmm

3. THE S-N DECOMPOSITION OF A MATRIX OF INFINITE ORDER 119

and Am = Sm + Aim be the S-N decompositions of Amm and Am, respectively, where Smm is an nm x nm diagonalizable matrix, Sm is an Nm x Nm diagonalizable matrix, Nmm is an nm x nm nilpotent matrix, Aim is an N. x Nn nilpotent matrix,

SmmNmm = NmmSmm, and SmAm = NmSm. The following lemma shows how the two matrices Sm and Aim look.

Lemma V-3-1. The matrices Sm and Aim have the following forms:

S1 = s11,

S"' =

rI Sm_1

0 l

(m > 2),

J

 

L Cm

Smm

 

 

[Aim_i

 

 

N, = Ni1,

Aim =

 

(m > 2),

fm

O

 

 

Nmm 1

 

where Cm and fm are am x N,_1 matrices.

Proof

Consider the case m > 2. Since the matrices Sm and Aim are polynomials in Am with constant coefficients, it follows that

Sm = [B m1

and Aim = I fm1 01,

 

ILM0

where 13m-1 and Dm-1 are Nm-1 x Nm_1 matrices, Cm and fm are nm x Nm_1 matrices, and µm and vm are nm x nm matrices. Furthermore, Dm_1 and Vm are nilpotent. Also, Bm-1Dm-1 = Dm-1Bm-i and tlmvm = 1.m µm. Hence, it suffices to show that Bm-1 and µm are diagonalizable.

Note that since Sm is diagonalizable, Sm has Nm linearly independent eigen

vectors. An eigenvector of Sm has one of two forms k and [?.] , where p is an

P J

eigenvector of Bm_1i whereas q is an eigenvector of µm. Therefore, if we count those independent eigenvectors, it can be shown that Bm-1 has N,, -I linearly independent eigenvectors, while Pm has nm linearly independent eigenvectors. Note that Nm = Nm_ 1 + nm. This completes the proof of Lemma V-3-1. 0

Lemma V-3-1 implies that the matrices Sm and Nn have the following forms:

S1 = SI1,

S11

0

0 ...

0

C21

S22

0 ...

0

Sm =

 

 

(m>2)

Cmi

Cm2 ... ...

Smm

and

 

 

 

N1 = N11,

0

D ...

0

Nil

Y21

JN22

0

0

Aim =

 

 

(m>2),

Fm 1 Fm2 ... ... Nmm

120

V. SINGULARITIES OF THE FIRST KIND

where C,k and Fjk are n3 x nk matrices.

Set

 

 

S11

0

0

 

 

 

C21

S22

0

 

 

S =

 

 

 

 

 

Cml

Cm2

"'

Smm 0

 

f

N11

0

0

"'

 

f21

N22

0

 

N=

 

 

...

.Wmm 0

 

 

1m 1

Fm2

Then,

 

 

 

 

 

(V.3.1)

A=S+N

and

SN=NS.

We call (V.3.1) the S-N decomposition of the matrix A.

V-4. The S-N decomposition of a differential operator

Consider a differential operator

(V.4.1)

C(yl = x

+ f2(x)Y,

w

 

 

 

where f2(x) _ > xt1

and 1

E

Let us identify a formal power series

t=o

 

ao

 

 

 

 

xt dt with the vector p =

a1

E C. Then, the operator C is represented

a2

t>o

by the matrix

flo

0

0

0 ...

f2l

In + f2o

0

0 ...

f12

l1

21, + f2o

0 ...

A =

11m f4.-l

fZn-2

... f21 mI, + fZa 0 .

5. A NORMAL FORM OF A DIFFERENTIAL OPERATOR

121

where I is the n x n identity matrix. Let A = S + N be the S-N decomposition

of A. Since A =

j) A1,11

where (i = R2 and I,> is the oo x oo identity

L

matrix, the matrices S and N have the forms

 

So

0

 

0

O

 

al

I. + So

0

0

 

a2

al

21n + So

0

(V.4.2)

S =

 

 

 

 

 

am

am-l

 

am-2

 

and

 

 

 

 

 

 

 

No

O

O

 

 

 

V1

No

0

 

 

N =

V2

V1

No

 

(V.4.3)

 

 

 

 

 

 

Vm

Vm-1

Vm-2 ...

... ...

... ...

... ...

at min + SO 0

VI No 0 ...

respectively, where the a. and v, are n x n matrices and S2o = So +No is the S-N decomposition of the matrix S2o.

00

00

Set a(x) = So + E xtat and v(x) = No +

trot. Then, S represents a

t_1

tvl

differential operator Lo[y1 = x + a(x)y, while N represents multiplication by v(x) (i.e., the operator: y" -+ v(x)y-). Since A = S + N and SN = NS, it follows

that

(V.4.4)

L[gf = Lo[Yj + v(x)y"

and

Lo[v(x)y1 = v(x)Lo[yl-

We call (V.4.4) the S-N decomposition of operator (V.4.1). We shall show in th next section that Lo[yj is diagonalizable and v(x)n = O.

V-5. A normal form of a differential operator

Let us again consider the differential operator

(V.5.1)

L[yj = Xf + Q(x)U,

00

where f2(x) = Ex1fle and Sgt E WC). In §V-4, we derived the S-N decomposi-

t=o

tion

(V.5.2)

L[y-] = Lo[yl + v(x)y"

and

Co[v(x)yl = v(x)Lo[yj.

122

V. SINGULARITIES OF THE FIRST KIND

 

where

 

 

 

 

(V.5.3)

Co[YI = x

+ a(x)y

 

and

 

 

 

Co

 

 

 

 

 

a(x) = So + 00 xtat

and

v(x) = No +

xtvt

 

t=1

 

 

t_1

(cf. (V.4.4)). Notice that N = So + No is the S-N decomposition of S2o and that the operator Go[ul and multiplication by v(x) are represented respectively by matrices (V.4.2) and (V.4.3).

Set

 

0

0

0

 

 

 

 

I. +,%

0

0

 

Sm =

al

21 + So

0

 

 

 

 

 

 

am_ 1

am-2

 

Since So is diagonalizable, there exist an invertible n x n matrix P0 and a diagonal matrix Ao such that

SoPo = PoAo, NO = diag[al, A2, ... , . \n)-

Note that A, (j = 1, 2,... , n) are eigenvalues of So. Hence, A. (j = 1, 2,... , n) are also eigenvalues of flo. For every positive integer m, Sm is diagonalizable. Hence,

 

Po

there exists an (m + 1)n x n matrix P.. =

Pi

such that SmPm = PmA0. If

 

Pm

m is sufficiently large, we can further determine matrices Pt for t > m + 1 by the equations

 

t

(V.5.4)

(tl + S0)Pt + E ahP1_h = PtAo.

h=1

Equation (V.5.4) can be solved with respect to Pt, since the linear operator Pt -+

(tl + So)P1 - PtAo is invertible if m is sufficiently large. In this way, we can find

A

P1

an oo x n matrix P =

I such that SP = PAa. Set P(x) = >x'P,. Then,

Pm

t=D

entries of P(x)-1 are also formal power series in x and

(V.5.5)

Co[P(x)j = P(x)Ao.

5. A NORMAL FORM OF A DIFFERENTIAL OPERATOR

123

Define two differential operators and an n x n matrix by

X [U1 = P(x)-iC(P(x)uZ,

/Co(u1 = P(x)-'Co(P(x)Uj,

{ vo(x) = P(x)-'v(x)P(x),

respectively. Then 1C(u'j = X 0(u'] + vo(x)u'. Observe that

.- - _.

_.. _, f _. , l du' \ _ . _. 1

(V.5.6)

=x du +Aou".

This shows that the operator Co(y1 is diagonalizable. Furthermore,

 

 

dvo

P(x)-'Lo(P(x)vo(x)1 = P(x)`Co[v(x)P(x)1

(V.5.7)

x

dxx) + Aovo(x) =

 

 

= P(x)-'v(x)Lo(P(x)1 =vo(x)P(x)-'Co(P(x)1= vo(x)Ao

This shows that the entry i,k(x) on the j-th row and k-th column of the matrix vo(x) must have the form v,k(x) = 7jkxAk-a', where ryjk is a constant. Since Lt(x) is a formal power series in x, it follows that

(V.5.8)

vjk(x) = 0

if

Ak - )j is not a non-negative integer.

Observe further that the matrix &. (x) can be written in the form

'711

'Yin

vo(x) = x-^°rx^°,

where

r =

'7ni Inn

Hence, for any non-negative integer p, we have vo(x)P = x-A0rPxAo, where

z^O = exp ((log x)Ao] = diag(za' , Z112'... , za^ I.

On the other hand, since No is nilpotent, vo(x)P can be written in a form i (x)P =

X'pQP(x), where mp is a non-negative integer such that lim mP = +oo and

P+oo

the entries of the matrix Qp(x) are power series in x with constant coefficients. Therefore, L0(x)P = 0 if p is sufficiently large. This implies that the matrix r is nilpotent. Since v(x) = P(x)vo(x)P(x)-1, we obtain v(x)^ = O.

Thus, we arrive at the following conclusion.

Theorem V-5-1. For a given differential operator (V.5.1), let 0o = So +No be the S-N decomposition of the matrix N. Then, there exists an n x n matrix P(x) such that

(1)the entries of P(x) are formal power series in x with constant coefficients,

(2)P(O) is invertible and SoP(0) = P(0)Ao, where A0 is a diagonal matrix whose diagonal entries are eigenvalues of Q0r

124 V. SINGULARITIES OF THE FIRST KIND

(5) the transformation

(V.5.9) P(x)u

changes the differential operator (V.5.1) to another differential operator

(V.5.10)

P(x)-1G(P(x)ui = KOK + vo(x)u,

where

Ko(61 = x dil + Aou,

 

 

vo(x) = x-^Orx^0,

and

r=

Ynl Inn

with constants "ljk such that

(V.5.11)

Vj k = 0

if

Ak - Aj is not a non-negative integer.

Furthermore, the matrix r as nilpotent.

Remark V-5-2. It is easily verified that the matrix P(x) is a formal solution of the system

xd) = P(x)(Ao + vo(x)) - Sl(x)P(x).

Since the entries of vo(x) are polynomials in x, the power series P(x) is convergent if f2(x) is convergent (cf. Theorem V-2-7). Therefore, in such a case, v(x) is convergent and, hence, o(x) is convergent.

Observation V-5-3. Choose integers l1, ... , to so that A.,+f j = Ak +fk if Aj -Ak is an integer. Then, vo(x) = x4rx-L, where L = diag(e1, f2, ... , en]. If we set

7t{V1 = x'! 1C(xf and ho{v") = x-LICo(x' i , it follows that

dv'

?{o(vl = x'L fXLTf + zLLv + Aoz' 7} = xdx + (Ao + L)v

Hence,

(V.5.12)

R{v = x + (Ao + L + r)v.

Note that

Ao + L = diag(A1 + P1, A2 + e2, ... , An + enj

and

(V.5.13)

(Ao + L)r = r(A0 + L).

It was already shower in §V-4 that r is nilpotent. Thus, the following theorem is

proved.

5. A NORMAL FORM OF A DIFFERENTIAL OPERATOR

125

Theorem V-5-4. The transformation y = P(x)xLiT changes the system

 

(V.5.14)

G[yi =

xdy + fl(x)yl = 0

 

 

to

 

 

 

 

(V.5.15)

VV! = x d1U

+ (Ao + L + r)u = 0.

 

Observation V-5-5. The matrix

 

 

 

 

 

n-1

h

 

o{x) =

x-Ao-L-r = x-Ao-L I +

(-1)h[logx[ rh

 

h!

 

 

 

h-1

 

 

 

 

 

is a fundamental matrix solution of system (V.5.15). Note that if (A3+tj)-(Ak+tk) is an integer, then A, + ej = Ak + tk in Ao + L.

Remark V-5-6. A fundamental matrix solution 4D(x) of system (V.5.14) is given by 4>(x) = P(x)x4o(x), which can be written in the form

45(x) = P(x)xLx-M-L In

n-1

(lo x h

rh

 

+ F{-1)h h!

 

 

 

 

IL

h=1

 

 

 

Since L - A0 - L = -A0, the matrix 44(x) can be also written in the form

 

 

 

 

 

n-1

 

 

 

4i(x) =

P(x)x-ne In +

(-1)h (lohx)h rh

 

 

 

 

 

 

h=1

 

 

 

 

 

I n-i

h

 

 

 

 

However, x-Ao and

I

nr >(-1)h [lohx[

rh do not commute, whereas

x-A°-L

 

h=1

 

 

 

 

 

n-I

 

 

 

 

 

 

 

[log x[ h

 

 

 

 

 

and In + E(-1)h

rh commute.

 

 

 

h=1

 

h!

I

 

 

 

 

 

 

 

 

 

 

Remark V-5-7. The methods used in §§V-3, V-4, and V-5 are based on the

original idea given in [GerL[.

Example V-5-8.

In order to illustrate the results of this section, consider the differential equation of the Bessel functions

(V.5.16)

za (zT) + (z2

- a2)y = 0,

 

 

where a is a non-negative integer. If we change the independent variable z by x = z2, (V.5.16) becomes

x d

2y=0.

lxdx

126

V. SINGULARITIES OF THE FIRST KIND

This equation is equivalent to the system

(V.5.17)

x

+ Q(x)y = 0,

y =

 

0

-1

 

where 11(x) =

x - a2

=10+x11i c o =

 

4

0

 

 

 

 

0 -1

 

0

0

a2

and 0i =

1

 

 

 

0

 

4

0

4

 

To begin with, let us remark that the S-N decomposition of the matrix 11 is given by S2o = So + No, where

O

if

80 = {Q

if

a = 0,

and

N

a > 0,

N0 - to

 

if

if

a = 0,

a>0.

 

a

0

 

Also, set Ao =

2

a

Note that two eigenvalues of Q are

0

2

Now, we calculate in the following steps:

Step 1. Fix a non-negative integer m so that m > 2 - (_a) = a.

Step 2. Find three 2(m + 1) x 2(m + 1) matrices A,,,, Sm, and Nm (cf. §§V-4 and

V-5).

Po

Step 3. Find a 2(m + 1) x 2 matrix Pm =

P1

 

, where the Pr are 2 x 2 matrices

such that SmPm = Pn&.

Pm

 

 

 

 

m

in

Step 4. Find two 2 x 2 matrices Mm(x) = No + >2xtvt and Qm(x) _ >xtPt

t=1

t-o

(cf. §§V-4 and V-5).

Step 5. We must obtain Qm(x)-1Mm(x)Qm(x) = vo(x) +O(xm+i), where

 

Po InoP0

if

a = 0,

vo(x)

r 0

a(a)xn

 

a>0,

 

if

 

0

0

 

J

 

where a(a) is a real constant depending on a (cf. Theorem V-5-1).

After these calculations have been completed, we come to the following conclu-

sion.

5. A NORMAL FORM OF A DIFFERENTIAL OPERATOR

127

+00

Conclusion V-5-9. There exists a unique 2 x 2 matrix P(x) = >xtPe such that

1=o

(1)the matrices Pl for t = 0,1, ... , m are given in Step 8,

(2)the power series P(x) converges for every x,

(8)the transformation y = P(x)u" changes (V.5.17) to

(V.5.18)

dii

xaj + (Ao + vo(x)) i = 0'.

We illustrate the scheme given above in the case when a = 0. First we fix m = 2. Then,

0 -1

0

0 0 0

 

 

 

0

0

0

0 0 0

0

0

0

0

0

0

 

 

 

0

0

0

0

0 0

0

0

1

-1

0 0

 

 

 

1

1

 

0

0 0

 

 

 

4

2

1

A2 =

 

 

 

 

 

 

 

 

 

 

 

 

0

0

1

0

0

 

S2 =

 

1

1

0

 

0

0

1

 

 

 

 

 

1

4

 

 

 

 

 

 

 

 

4

4

 

 

 

 

0

0

0

0 2 -1

 

 

 

3

3

1

1

2

0

 

 

 

 

 

 

 

 

-

32

32

4

2

 

 

 

 

 

 

 

 

 

 

0

0

4

0

0

2

 

 

1

3

1

1

0

2

 

 

 

 

 

 

 

 

16

32

4

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

and

 

 

 

 

 

-1

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0

0

0

0

 

 

 

 

 

 

 

 

0

 

0

0

0

0

0

 

 

 

 

 

 

 

 

1

 

1

0

-1

0

0

 

 

 

 

 

 

 

 

4

 

2

 

 

 

 

 

 

N2=

 

 

 

 

 

 

 

 

 

 

 

0

 

1

0

0

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

3

1

1

0

-1

 

 

 

 

 

 

 

 

32

 

32

4

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

16

 

T2-

0

9

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1192

-320

 

 

 

 

 

 

 

 

 

 

 

 

64

-128

 

Calculating eigenvectors of S2, we find a 6 x 2 matrix P2 =

-16

 

16

such

-32

 

48

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

1

 

 

 

 

 

 

 

 

 

 

 

 

1

 

0

 

that S2P2 = 0. Note that, in this case, Ao = 0.