0387986995Basic TheoryC
.pdf118 |
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.
