14
Tekil Nokta Etrafında Çözüm: Frobenius Yöntemi
Önceki kısımda verilen diferansiyel denklemlerin çözümleri, bu diferansiyel denklemin adi noktaları civarında seri çözümleri bulunmuştur. Bu kısımda diferansiyel denklemin tekil noktası civarında çözümler ele alınacaktır. Buna göre
1
0( ) 1( ) 1 ... ( ) 0
n n
n n n
d y d y
a x a x a x y
dx dx
şeklinde verilen diferansiyel denklemin
x x
0 noktası tekil nokta olsun. Genelliklex
0 civarında kuvvet serisi biçiminde çözümü yoktur. Belirli koşullar çerçevesinde çözüm0 0
0
( ) ( )r n( )n
n
y x x x C x x
şeklinde aranacaktır. Burada
Cn ve r sabitlerdir (reel veya kompleks sayı olabilir).
Tekil Noktaların Sınıflandırılması 0 1
( x x P x ) ( )
ve( x x
0)
2P x
2( )
fonksiyonları analitik isex
0 noktasına regüler tekil nokta denir. Aksi halde irregüler tekil nokta denir.Örnek 1.
2 2
2 d y2 dy 5 0
x x x y
dx dx diferansiyel denklemi için x=0 noktasının özelliğini belirleyiniz.
2
2 2
1 5
2 2 0
d y dy x x dx y
dx x
için 1
( ) 1 P x 2
x
ve
2 2
( ) 5 2 P x x
x
olur.
1
( ) 1 xP x 2
ve 2 2( ) 52
x P x x fonksiyonları x=0 noktasında analitiktir. Dolayısıyla x=0 noktası diferansiyel denklemin regüler tekil noktasıdır.
Örnek 2.
2
2 2 2 2
2 1
0
( 2) ( 2)
d y dy x
dx y
dx x x x x
diferansiyel denklemi için x=0; 2 noktalarının özelliklerini belirleyiniz
15
1 2
( ) 2
2 P x
x x
ve
2 2 2
( ) 1
2 P x x
x x
olur. Buradan,
1 2
(x 2)P x( ) 2 x
ve
2 2 22 ( ) x 1
x P x
x
için bu fonksiyonlar x=2 noktasında süreklidir,
dolayısıyla x=2 noktası diferansiyel denklemin regüler tekil noktasıdır.
1
( ) 2
( 2) xP x x x
ve
2
2 2
( ) 1
2 x P x x
x
için x=0 noktasında
xP x
1( )
tanımsızdır, dolayısıyla x=0 noktası diferansiyel denklemin irregüler tekil noktasıdır.Teorem.
x x
0 diferansiyel denklemin regüler tekil noktası olsun. Bu durumda,0 0
0
( )r n( )n
n
x x C x x
,0 x x
0 R
serisi şeklinde
x
0 civarında sıfır olmayan yakınsak en az bir çözümü vardır.Örnek 3.
2 2
2 d y2 dy 5 0
x x x y
dx dx
,
0 x R
için denkleminix x
00
noktası civarında seri çözümünü bulunuz.
2
2 2
1 5
2 2 0
d y dy x x dx y
dx x
; 1 1
( ) 2
P x x
ve 2 25
( ) 2 P x x
x
için x=0 regülar tekil noktadır.
Dolayısıyla çözüm için
0 0
( ) r n n n n r
n n
y x x C x C x
şeklinde aranır.Frobenius Yöntemi:
Farzedelim ki,
x
0 regüler tekil noktası olsun.
0
0
0
0 0
( ) r n n n n r
n n
y x x x C x x C x x
,C
00
(5)16
0
10
( )
n n rn
dy x n r C x x
dx ,
2 2
2 0
( ) 1
n n rn
d y x
n r n r C x dx
Diferansiyel denklemde yerine yazılır ve
xx0
nin derecesine göre gruplandırılır.
1
20 0 r k 1 0 r k 2 0 r k
...
K x x
K x x
K x x
=0 bulunur. Burada k belirli bir tam sayıdır.Sağ taraf sıfır olduğu için bu durumda
K
0 K
1 K
2 ... 0
olmalıdır.Bu eşitliklerden
K
00
den r için ikinci dereceden denklem elde edilir, buna indeks denklemi denir. Bu denklemin köklerir
1 ver
2(r
1>r
2 olsun) bulunur ve bu köklere diferansiyel denklemin üstelleri (ekponentleri) denir.1)
K K
1,
2,...
ler sıfıra eşitlenerek, bunların ifadelerinde yer alanC
n ler için denklemler elde edilir.2)
r r
1, (5) de elde edilen denklemlerde yerine yazılarakC
n ler bulunur.3) Eğer
r
2 r
1ise, (6) deki işlemlerr r
2 için tekrarlanarak, diğer seri çözüm bulunur.Ancak bu çözüm
r r
1için belirlenen çözümden lineer bağımsız olmayabilir. Bu durum daha sonra araştırılacaktır.Örnek 4.
2 2
2 d y2 dy 5 0
x x x y
dx dx
0
( )
n n rn
y x C x
,C
00
;
10
( ) n n r
n
y x n r C x
;
20
( ) 1 n n r
n
y x n r n r C x
için diferansiyel denklemde yerlerine yazılırsa,
0
0 1
1
2 ( 1) 5 2 1 ( ) 5 0
n
r n r
n n
K n K
r r r C x n r n r n r C C x
olur.17
0 2 ( 1) 5 0
K r r , r 1 5
r
2 ,r
21
n 0
K dan,
1 52
r
için 1(2 7)
n n
C C
n n
,n 1
,y x
1( ) C x
0 5 2 1
19x
1981x
2
77221x
3 ...
2
1
r
için 1(2 7)
n n
C C
n n
,n 1 y
2( ) x C x
0 1 1
15x
301x
2
901x
3 ...
Her iki çözümdeki C birbirinden farklı sabitlerdir. Genel çözüm, 0
1 2
5 2 1 1 2 1 3 1 1 1 2 1 3
01 9 198 7722 02 5 30 90
( ) ( )
( ) 1 ... 1 ...
y x y x
y x C x x x x C x
x x x
,yazılabilir. Burada
C
01 veC
02keyfi sabitlerdir.Burada “her zaman birbirinden bağımsız iki seri çözüm bulunabilir mi?” sorusu ortaya çıkar veya 1. Hangi koşullar çerçevesinde verilen
1
0( ) 1( ) 1 ... ( ) 0
n n
n n n
d y d y
a x a x a x y
dx dx
diferansiyel denklemin
0
0
0
r n
n n
x x C x x
, x civarında lineer bağımsız iki çözümü vardır? x02. Eğer bu diferansiyel denklemin x civarında yukarıda verilen formda lineer bağımsız iki x0 çözümü yoksa, lineer bağımsız iki çözüm nasıl bulunur?
Teorem.
Farzedelim ki,
r
1 r
2 bilinen denklemin kökleridir.i.
r
1 r
2N
(N=0,1,2,… pozitif tam sayı) bu durumda,a. 1
0
1
0
0
( ) r n n
n
y x x x C x x
b. 2
0
2 *
0
0
( ) r n n
n
y x x x C x x
18
biçiminde iki lineer bağımsız çözümü vardır. Burada,
C
00
,C
0*0
dır.ii. Eğer,
r
1 r
2N
(N=0,1,2,… pozitif tam sayı) ise, bu durumdaa. 1
0
1
0
0
( ) r n n
n
y x x x C x x
,b. 2
0
2
0
1 00
( ) r n n ( ) ln
n
y x x x C x x C y x x x
biçiminde lineer bağımsız ikiçözümü vardır. Burada,
C
00
,C
0*0
ve C sabit/sıfır olabilir.iii. Eğer,
r
1 r
20
ise bu durumdaa. 1
0
1
0
0
( ) r n n
n
y x x x C x x
,C
00
veb. 2
0
2 1
0
1 00
( ) r n n ( ) ln
n
y x x x C x x Cy x x x
biçiminde lineer bağımsız ikiçözümü vardır. Burada,
C
00
,C
0*0
ve C sabit/sıfır olabilir.Örnek 5.
2 2 2
2 d y2 dy 3 0
x x x y
dx dx
diferansiyel denkleminin x=0 regüler tekil noktası civarında çözümünü bulunuz.
0
( ) n n r
n
y x C x
;
10
( ) n n r
n
y x n r C x
;
20
( ) 1 n n r
n
y x n r n r C x
için
20 0 0 0
2 1 3 0
n n r
n n r
n n r
n n rn n n n
n r n r C x n r C x C x C x
1
0 1
0 0
2
2 0
2 1 3 2 1 1 3
2 1 3 0
r r
n r
n n
n
r r r C x r r r C x
n r n r n r C C x
19
2r r için 1 r 3 0
2 r
2 r 3 0
buradan 1 3 22
, 1
r r
1 2 520
r r
ve 5/2 tam sayı ( N 0,1, 2,...
) değil.
12r r 1 r 1 3 C 0
22 n r n r 1 n r 3 CnCn 0
,
n 2
1 32
r r
için7 C
10
ise,C
10
2 5
n n 2 0n n C C ,
n 2
2
n 25
n
C C
n n
,n 2
için 2 018
C C
; 3 10 33
C C
; 4 0936 C C
,…1 32
r r
için çözüm;3 2 2 4
1 0
1 1
( ) 1 ...
18 936
y x C x x x bulunur.
2
1
r r
için*
3 C
10
ise,C
1*0
2 5
n n 2 0n n C C ,
n 2
2
n 25
n
C C
n n
,n 2
için* 0*
2
2
C C
;* 1*
3
0
3
C C
;*
* *2 0
4
12 24
C C C
,…2
1
r r
için çözüm;* 1 2 4
2 0
1 1
( ) 1 ...
2 24
y x C x x x bulunur. Dolayısıyla genel çözüm,
1 1 2 2
( ) ( ) ( )
y x K y x K y x
olur. K1C0 ve K2C*0 alınmıştır.20 Örnek 6.
2 2 2 5
4
2 0
d y dy
x x x y
dx dx diferansiyel denkleminin x=0 regüler tekil noktası civarında çözümünü bulunuz.
0
( ) n n r
n
y x C x
;
10
( ) n n r
n
y x n r C x
;
20
( ) 1 n n r
n
y x n r n r C x
için
20 0 0 0
1 5 0
4
n n r
n n r
n n r
n n rn n n n
n r n r C x n r C x C x C x
20 2
1 5 0
4
n r n r
n n
n n
n r n r n r C x C x
5 5 1
0 1
4 4
0 0
5 2
4 2
0
1 1 1
1 0
r r
n r
n n
n
r r r C x r r r C x
n r n r n r C C x
2 5
2 4 0
r r için 1 5 2 1
2
,
2r r
1 2
3
r r
tam pozitif sayıdır. Dolayısıyla yukarıda verildiği üzere:ii. Eğer,
r
1 r
2N
(N=0,1,2,… pozitif tam sayı) ise, bu durumdaa. 1
0
1
0
0
( ) r n n
n
y x x x C x x
,b. 2
0
2
0
1 00
( ) r n n ( ) ln
n
y x x x C x x C y x x x
biçiminde lineer bağımsız ikiçözümü vardır. (Burada C sıfır/sıfırdan farklı olabilir) 1 52
r r
içinC
10
3
n n 2 0n n C C ,
n 2
n 23
n
C C
n n
,n 2
için 2 02 5 C C
; 3 10 3 6 C C
; 44 7
22 4 5 7
0C C C
,21
0
2n
2 4 6 ... 2 5 7 9 ... 2 3 C C
n n
,n 1
1 52
r r
için çözüm;
5 2 2
1 0
1
( ) 1
2 4 6 ... 2 5 7 9 ... 2 3
n
n
y x C x x
n n
bulunur.2 21
r r
için (Bu kısımdaki C katsayıları üst indis * lı olup, yazım kolaylığı açısından yazımda ihmal edilmiştir.) y2(x) çözümünü ilk kökte olduğu gibi arayalım. Sonuçta lineer bağımsız çözüm bulunursa, çözüm tamamlanmış olur, aksi halde (ii) de belirtildiği gibi bir ek terim gelmelidir.1
0 C
3
n n 2 0n n C C ,
n 2
n 23
n
C C
n n
,n 2
ven 3
(C
3 başka yerden bulunacak)n=2 için 2 0
2
C C
; 4 2 04 2 4
C C C
; 52 5
3C C
,…
0
2n
2 4 6 ... 2 3 5 7 ... 2 3 C C
n n
,n 3
3
2 1
2 4 6 ... 2 2 5 7 9 ... 2 1
n
C C
n n
,n 2
2
1 2
r r
için çözüm;
2 4 2
1 2
2 0
3
2 1 1 2 3
3
2
( ) 1
2 2 4 2 4 6 ... 2 3 5 7 ... 2 3
2 4 6 ... 2 2 5 7 9 ... 2 1
n
n
n
n
x x x
y x C x
n n
C x x x
n n
bulunur. Burada
C
0( C
0*)
veC
3( C
3*)
keyfi sabitlerdir.y x
1( )
vey
2( ) x
aralarında lineer bağımlıdır fakat iki tane sabit içermektedir. Bu nedenle çözümlerden biri (örneğiny
2( ) x
tek22
başına) genel çözüm olabilir. Dolayısıyla, yukarıda (ii) maddesinde yer alan
y
2( ) x
çözümünde doğal logaritmanın katsayısı C=0 olur.Örnek 7.
2 2 2
2 3 3 0
d y dy
x x x y
dx dx
diferansiyel denkleminin x=0 civarında Frobenius yöntemi ile çözünüz .
0
( ) n n r
n
y x C x
;
10
( ) n n r
n
y x n r C x
;
20
( ) 1 n n r
n
y x n r n r C x
için
10 1
1 3 3 n n r 1 n n r 0
n n
n r n r n r C x n r C x
0
1
0 1 0
1 3 3 r 1 3 3 n 1 n n r 0
n
r r r C x n r n r n r C n r C x
2 4 3 0
r r için
r
1 3, r
2 1
1 2
2
r r
tam pozitif sayı (2 N
) 13
r r
için
2
n
2
n 1 0n n C n C ,
n 1
1n n
C C
n
,n 1
içinC
1 C
0 ; 2 1 02 2!
C
C C
; 3 2 03 3!
C C C
,...2 0
( 1)
!
n n
C C
n
,n 1
1
3
r r
için çözüm;2 3
3 3
1 0 0
( ) 1 ... ( 1) ...
2! 3! !
n n x
x x x
y x C x x C x e
n
bulunur.
2
1
r r
için (not edelim ki sabitler üst indis * içermektedir, yazım kolaylığı açısından yazımda ihmal edilmiştir)
2
n n 1 0n n C nC ,
n 1
23
n2
1
n
C C
n
,n 1
ven 2
n=1 için
C
1 C
0 ;n=2 için
0 C
2 2 C
1 0
içinC
10
dolayısıyla,C
00
bulunur. Ancak çözüm gereği 00
C
olmalıdır. Buradan,r r
21
için bir özel çözüm yoktur. Böylece C2 keyfi kabul edilerek,3 2
C C
; 4 33
C C
; 5 4 23 3!
C C
C
,… 2( 1)
2!
n n
C C
n
,n 1
2
1
r r
için çözüm;4 5 2 2 3
2 3 3
2 2 2
( 1) ( 1)
( ) ... ... 1 ... ...
2! 3! ! 2! 3! !
n n n n
x x x x x x
y x C x x x C x x
n n
2
( )
2 3 x 1( )
y x C x e
y x
bulunur. Teoremden görüldüğü gibi ikinci çözümy
2( ) x
* 1
1 0
( ) ln
n n n
C x Cy x x
,C
0*0
veC 0
olmalıdır. Bu çözümü belirlemek içi pek çok yöntem vardır, bunlardan biri de türev mertebesini düşürmedir.MERTEBE DÜŞÜRME METODU
Bu metoda göre diferansiyel denklemin bir çözümü y1(x) yardımıyla çözüm,
( )
3 x( )
y x x e
u x
, (y x
1( ) x e
3 x) şeklinde aranır. Burada u(x)=?
3 x
3
2 x 3 x( )
dy du
x e x e x e u x
dx dx
2 2
3 2 3 3 2
2 x 2
2 3
x x x6
x6
x( )
d y d u du
x e x e x e x e x e xe u x
dx dx dx
Bu ifadeler 2d y22
2 3
dy 3 0x x x y
dx dx
denkleminde yerine yazılır düzenlenirse,
2
2 3 0
d u du
x x
dx dx
;
( ) du x dx
değişken dönüştürmesi kullanılarak,24
3
( ) 0xd x x
dx
için
( ) x x e
3 x veu x ( ) ( ) x dx x e dx
3 x3 3
( )
2( )
x xy x y x x e
x e dx
olur. Bu fonksiyony x
1( )
ile lineer bağımsızdır.2 3 4
3 3
2
( ) 1 ...
2 6 24
x
x x x
y x x e
x
x dx
İntegral alınırsa,
3 2
2 2
1 1 1 1 1
( ) ln ...
2 6 48
2
y x x e
xx x x
x x
3 2 3
2 2
1 1 1 1 1
( ) ... ln
6 48 2
2
x x
y x x e x x x e x
x x
bulunur. Dolayısıyla genel çözüm,
1 1 2 2
( ) ( ) ( )
y x K y x K y x
olur.
Yukarıdaki örneklerde teoremde verilen
r
1 r
20
veyar
1 r
2 hariç tüm haller ele alınmıştır.Hatırlatalım ki, bu durumda çözümler
1 0 0
0
( )
r n nn
y x x x C x x
ve
1 *
2 0 0 1 0
0
( )
r n n( ) ln
n
y x x x C x x y x x x
,(
C
0 0 ve C
0* 0
) şeklinde olur.ÖRNEKLER
Aşağıdaki denklemlerin tekil noktalarını sınıflandırınız.
i.
2 3
d y22
2
dy 0x x x y
dx dx
ii.
4 2 3 2
d y22 2
1
dy 2 0x x x x x y
dx dx
25 iii.
5 4 6 3
d y22 2dy
2
0x x x x x y
dx dx
Aşağıdaki denklemleri, regüler tekil noktaları civarında Frobenius yöntemi ile çözünüz.
iv. 2 2d y22 dy
2 1 0x x x y
dx dx
v. 2d y22 dy
2 19
0x x x y
dx dx
vi. 2d y22 dy
2 14
0x x x y
dx dx vii.
2 2 2
2 2 0
d y dy
x x y
dx dx
viii. 2d y22
3
dy 3 0x x x y
dx dx
ix.
2
2 2 0
d y dy
x y
dx dx
BESSEL DENKLEMİ VE BESSEL FONKSİYONLARI A. Sıfırıncı Dereceden Bessel Denklemi
p bir parametre olmak üzere,
2 2 2 2
2 0
d y dy
x x x p y
dx dx
(1)
bu denkleme p. dereceden Bessel denklemi denir. Denklemin çözümüne ise p. dereceden Bessel fonksiyonu denir.
Eğer, p=0 olursa (1) den,
2
2 0
d y dy
x xy
dx dx
(2)
denklemi elde edilir. (2) denklemi sıfırıncı dereceden Bessel denklemi adlanır. Bu denklem için x=0 regüler tekil noktadır ve çözüm,
0
( )
n n rn
y x C x
,C
00
(3)26 şeklinde aranır. (2) ve (3) den
1
1 10 0 0
1 n n r n n r n n r 0
n n n
n r n r C x n r C x C x
2 1 2 10 2
n r n r 0
n n
n n
n r C x C x
2
22 1 1
0 1 2
0 0 2
0
1 0
r r n r
n n
n
r C x r C x n r C C x
2
0
r
içinr
1 r
20
bulunur. 1 r
2C
1 0
2 n n 2
0
n C C
,n 2
22 n n
C C
n
,n 2
; 2 02
2C C
, 3 12
0
3
C C
, 4 2 02 2 2
4 2 4
C C C
, …2 2 20
( 1) ( !) 2
n
n n
C C
n
,n 1
r=0
0
2
1 0 2 0 0
0
( )
( ) ( 1) ( )
( !) 2
n n
n
J x
y x C x C J x
n
;2
0 2
0
( ) ( 1)
( !) 2
n n
n
J x x
n
toplamına 0. dereceden1. tür Bessel fonksiyonu denir.
r
1 r
20
olduğu için yukarıda verilen teoreme göre diğer çözüm,*
2 0
0
( )
n n( ) ln
n
y x x C x J x x
dır.Şimdi, mertebe düşürerek bu çözümü bulmaya çalışalım:
y x
2( ) y x u x
1( ) ( )
2 1 1
( ) ( )
( ) dy x ( ) ( ) du x
y x u x y x
dx dx
2 2
1 1
2 2 1 2
( ) ( ) ( ) ( )
( ) d y x ( ) 2 dy x du x ( ) d u x
y x u x y x
dx dx
dx dx
27 Diferansiyel denklemde yerine yazlırsa: 2
2 0
d y dy
x xy
dx dx
2
1
d u
2 12 dy
1du 0
xy y x
dx dx dx
içndu
dx
ve 22
d u d dx dx
dönüşümü sonunda,
0
( ) d
0( ) 2 d
0( ) 0
xJ x J x x J x
dx dx
0 0
0
0 0
( ) 2 ( )
( ( ))
( ) 2 ( )
J x x d J x
d J x
d dx dx dx
xJ x x J x
0
2ln ln x ln J ( ) x
0
21 ( )
du x J x dx
için
0
2
0
2( ) 1
( ) ( )
dx
e
xu x dx dx
x J x J x
için,
2 0 2 0 2
0 0
( ) ( ) ( )
( ) ( )
dx
e
xdx
y x J x J x
J x x J x
0
2 2 4 63 5
( ) 1 ...
2 32 576
x x x
J x
2 4 6
2 0
1 5 23
1 ...
2 32 576 ( )
x x x
J x
3 5
2 0
1 5 23
( ) ( ) ...
2 32 576
x x x
y x J x dx
x
2 4 6 2 4 6
0
5 23
( )ln 1 ... ...
4 64 2304 4 128 3456
x x x x x x
J x x
2 4 6
0