Bir Sınıf Jacobi Matrisi İçin Özdeğer Problemi
1Ozan ÖZKAN∗
Selçuk Üniversitesi, Fen-Edebiyat Fakültesi, Matematik Bölümü 42079 Kampüs, Konya
Özet : Bu çalışma;
n
×
n
simetrik Jacobi matrislerinin özdeğerleri üzerine bazı yeni sonuçları içermektedir. Ele alınan problem; simetrik Jacobi matrisler ailesinin özel bir halidir. Burada ele alınan simetrik Jacobi matrisi, bir sınıf hiperbolik tip diferensiyel denklemin fark denklemi hale getirilmesi sonucu oluşan katsayılar matrisi ile aynıdır [4]. Elde edilen sonuçlar; bazı diferensiyel denklem sistemlerinin çözümünün davranışını irdelemeye imkan verir.Anahtar sözcükler: Özdeğer, simetrik Jacobi matrisi, Diferensiyel Denklemler
Eigenvalue Problem For A Class Of Jacobi Matrices
Abstract: This study contains some new results about the eigenvalues of a
n
×
n
symmetric Jacobi matrix. The problem is a special kind for the of family of the symmetric Jacobi matrices. The symmetric Jacobi matrix in this paper, is the same as the coefficient matrix obtained by converting one class of hyperbolic type differential equation into difference equation [4]. The obtained results enable to analyses the behavior of the solution of the system of some differential equations.
Key words: Eigenvalue, symmetric Jacobi matrix, Differential equations
1.Giriş
Bu çalışmada bir sınıf simetrik Jacobi matrisinin özdeğerleri ve bu özdeğerlerin davranışı incelenecektir. Jacobi matrislerinin çözümlerinin varlık ve tekliklerinin temel teorisi litaratürde mevcuttur [5,6]. Simetrik Jacobi matrisleri genel olarak,
J = nxn n n n c b b c b c b − − a . . . 0 . . . . . . . . . . . 0 ... 0 a 0 0 ... 0 a 0 ... 0 a 1 1 3 3 2 2 2 1 1 1 (1.1)
biçiminde verilen ve
c
i=
b
i,
i
=
1
,
n
−
1
şartını sağlayan matrisler olarak tanımlanmaktadır [7]. İkinci bölümde; (1.1) ile gösterilen simetrik Jacobi matrisilerinden özel bir matris ele alınacaktır. Bu matris; ƒ(x) ve g(x) fonksiyonları sürekli fonksiyonlar, a, b katsayıları reel sayı vel
> 0 için 2 2t
U
∂
∂
- a2 2 2x
U
∂
∂
+ bt
U
∂
∂
= 0 , (t, x)∈Q≡R+ × [0,l
] (1.2) U(t, 0) = U(t,l
) = 0 , t ≥ 0 (1.3) U(0, x) = ƒ(x) Ut (0, x) = g(x) , x∈ [0,l
] (1.4)biçiminde verilen lineer hiperbolik diferensiyel denkleminin aşağıda gösterileceği gibi fark denklemi halinde yazılması sonucu elde edilen katsayılar matristir [4]. Yeterli kadar büyük n doğal sayısı için [0,
l
] aralığında aşağıdaki gibi noktalar belirleyip uzaklığa bağlı sonlu farklar yöntemi kullanılarak; (1.2)-(1.4) problemi sonlu farklar problemine dönüştürülsün [1,8].xj = j
1
+
n
l
, j =0
,
n
+
1
olmak üzere (1.2)-(1.4) problemine uygun sonlu farklar yöntemiyleoluşturulmuş problem; 2 2 2
)
1
(
)
,
(
)
,
(
l
a
n
x
t
U
b
x
t
U
′′
i+
′
i−
+
(U
(
t
,
x
i+1)
−
2
U
(
t
,
x
i)
+
U
(
t
,
x
i−1)
)= 0, i =1
,
n
(1.5) U(t,x0)=U(t, xn+1) =0 (1.6)U(0, xi ) = f(xi) ,
U′
(0, xi) = g(xi ) i =1
,
n
(1.7)olur. Yukarıda her bir i değeri için n tane denklemden elde edilen sistemde;
= ) ( . . . ) ( 1 n x U x U U , = ) ( . . . ) ( 1 n x f x f f ve = ) ( . . . ) ( 1 n x g x g g
eşitlikleri ile 1×n boyutlu vektörleri,
A = − − − − 2 1 . . . 0 0 1 . . . . . . . . . . . . . . . . . . . . 0 . 0 1 2 1 0 0 . . 0 1 2 1 0 . . . 0 1 2 (1.8)
ile de n×n boyutlu matrisi gösterilirse (1.5)-(1.7) probleminden aşağıdaki ifadeyi elde edilir [9,10,11].
(
)
0
1
2 2 2=
+
−
′
+
″
U
A
a
n
U
b
U
l
(1.9)( )
f
U
0
=
,U
′
( )
0
=
g
(1.10)Bu çalışmada (1.9)-(1.10) probleminin çözümünden ziyade (1.9) denklemindeki (1.8) ile gösterilen A matrisi incelenecektir. Ele alınan matrisin karakteristik polinomu vasıtasıyla sırasıyla; özdeğerlerin tanımlanmış oldukları aralık, bu aralıkta nasıl sıralandıkları ve davranışları incelenmiştir. Elde edilen sonuçlar matrisin elemanlarına bağlı olduğundan, bu çalışma matrisin elemanlarına bağlı olarak özdeğerlerin sınıflandırılması için yeni bir karakterasyon vermiş olacaktır.
2. Simetrik üçlü bant matrisin öz değerleri
Bu bölümde n-mertebeden (1.8) matrisinin spektral yapısı incelenecektir.
A = nxn − − − − 2 1 . . . 0 1 . . . . . . . . . . . 0 ... 0 1 2 1 0 0 ... 0 1 2 1 0 ... 0 1 2
Biliniyor ki keyfi A matrisinin özdeğerlerini ve özvektörlerini bulmak için,
A ϕ = λ ϕ
denkleminden faydalanılır [1,2,3]. Bu sistem homojen sistem olduğundan; bilindiği gibi sıfırdan faklı çözümünün olabilmesi için A-λE matrisinin determinantının sıfıra eşit olması lazımdır. Bunun için aşağıdaki formda tanımlanmış Pn(λ) polinomunu incelemek yeterlidir.
Pn (λ) = λ λ λ λ − − − − − − − − 2 1 0 . . . 0 1 . . . . . . 0 . 0 1 2 1 0 0 . . 0 1 2 1 0 . . . 0 1 2 (2.1) Lemma 2.1. ∀ n ≥ 3 için; Pn (λ) = ( -2 - λ) Pn-1 (λ) - Pn-2 (λ) (2.2) eşitliği doğrudur.
İspat : Yukarıdaki (2.1) determinantının değeri birinci satırı kullanarak hesap edilirse;
Pn (λ ) = (-2 - λ ) ( - 1)1+1 ) 1 ( ) 1 ( 2 1 0 0 1 0 . . . 0 ... 1 2 1 0 ... 0 1 2 − − − − − − − − n x n λ λ λ + 1( -1)1+2 ( 1) ( 1) 2 1 0 ... 0 1 . . . . . . . . . . . . . . 0 ... 0 1 2 0 0 ... 0 1 1 − − − − − − n x n λ λ
olur. Elde edilen yukarıdaki eşitlikteki toplamın birinci determinantı Pn–1(λ)’ ya eşittir. Toplamın ikinci determinantının değeri ise; birinci sütun kullanılarak hesap edilirse değerinin Pn–2(λ)’ ya eşit olduğunu görülür. Böylelikle (2.2) bağıntısının doğruluğu ispatlanır.
Sonuç 2.2.
∀ n ≥ 1 için Pn(- 4 - λ) = (-1)n Pn (λ) (2.3) dır.
İspat :
P
n(
λ
)
’ nın tanımından ,P1 (λ) = -2-λ, P2 (λ) = (2 + λ)2 –1
değerleri kolaylıkla bulunabilir. Bu değerler indüksiyon metodunun ilk aşaması olarak (2.3)’ ün doğruluğunu kontrol için kullanırsa;
eşitlikleri bulunur. Bu ise n = 1 ve n = 2 için (2.3) ün doğru olduğunu gösterir.İkinci aşama olarak varsayalım ki; n = k-1 ve n = k için (2.3) eşitliği doğru olsun. n = k+1 için doğru olduğu gösterilmelidir. (2.2) formülüne göre,
Pk+1 (-4-λ) = (2+λ) Pk (-4-λ)-Pk - 1 (-4-λ) = (2+λ) (-1) Pk (λ)-(-1)k - 1Pk - 1 (λ)
= (-1)k + 1[ (-2-λ) Pk (λ)-Pk - 1 (λ)] = (-1)k + 1Pk + 1 (λ)
olur. Bu ∀ n ≥ 1 için (2.3) eşitliğinin doğru olduğunu gösterir. Sonuç 2.3.
∀ n ≥ 1 için aşağıdaki eşitlikler doğrudur
Pn (0) = (-1)n (n + 1) (2.4)
Pn (- 4) = n + 1 (2.5)
İspat: (2.5) eşitliğinin doğruluğu (2.3) ve (2.4) den kolayca görülür. Bu nedenle sadece (2.4) eşitliğinin doğru olduğunun ispatlanması yeterli olacaktır. Yine indüksiyon metodunu kullanılırsa,
n = 1 ve n = 2 için; P1 (0) = -2 , P2 (0) = 3
olur. (2.4) eşitliğinin n = k-1 ve n = k için doğru kabul edip, n = k+1 için doğru olduğu gösterilmelidir. (2.2) den dolayı ;
Pk+1 (0)= -2 Pk (0) - Pk-1 (0) = (-2)(-1)k (k+1)-(-1)k-1 k = (-1)k+1 (2k+2-k)
= (-1)k+1 ((k+1)+1)
olur. Bu ise (2.4) ün n = k+1 için doğru olması demektir. Böylelikle Sonuç 2.3’ ün ispatı yapılmış olur.
Teorem 2.4.
Pn(λ) ve Pn-1 (λ) polinomlarının kökleri aşağıdaki formda sıralanmışlardır.
0
...
4
<
λ
i(n)<
λ
i(n 1)<
<
λ
n i(n 1)<
λ
n(n)<
−
− − − , i =1
,
n
n=2k+1, ∀k ≥ 0 için , Pn (-2) = 0ve kökler
λ
=
−
2
noktasına göre simetriktir.İspat: Köklerin
λ
=
−
2
noktasına göre simetrikliği Sonuç 2.2.’den görülebilir. Lemma 2.1’ deki (2.2) eşitliğinden dePn (-2) = - Pn-2 (-2)
dir. Buradan her k ≥ 1 için
P2k-1(-2) = - P2k-3(-2) = ... = (-1)k-1P1(-2) = 0
olur.Teoremin esas hükmünü ispatlamak için burada da indüksiyon metodu kullanılacaktır. n = 1 ve n = 2 için;
P1(λ) = 0 ⇒ -2-λ = 0 ⇒ λ = -2
P2(λ) = 0 ⇒ (-2-λ)2 = 1 ⇒ λ1 = -3, λ2 = -1
dir. Pk-1 (λ) ve Pk (λ) için teoremin hükmü doğru olsun. Yani, bu polinomların sırasıyla (k-1) ve (k) sayıda negatif, aynı zamanda
λ
=
−
2
noktasına göre simetrik olan kökleri mevcut olsun vePk-1(λ) nın kökleri Pk(λ) nın kökleri arasına girsin. Köklerin negatif ve
λ
=
−
2
’ ye göre simetrikliğinden anlaşılıyor ki, kökler (-4, 0) aralığında tanımlıdır. Bu kökler sırasıyla,λi(k-1) , i=
1
,
k
−
1
λi(k) , i =1
,
k
şeklinde gösterilirse,0
...
4
<
( )<
( 1)<
<
( 1)<
( )<
−
− − − k k k i k k i k iλ
λ
λ
λ
olur. Bu ise ;)
1
-k
1,
(i
4
)
k
1,
(i
4
) 1 ( ) 1 ( ) ( 1 ) (=
−
−
=
=
−
−
=
− − − − + k i k k i k i k k iλ
λ
λ
λ
demektir.Buradan
P
k+1(
λ
)
için teoremin hükmünün doğruluğu görülür. Kökler simetrik olduğundan (-4,-2) veya (-2,0) aralıklarından sadece birinde inceleme yapmak yeterli olacaktır.P
k(
λ
)
veP
k−1(
λ
)
polinomlarını çarpanlarına ayrıldığında,
−
−
−
−
=
−
−
−
−
=
− − − − − −)
)...(
)(
(
)
1
(
)
(
)
)...(
)(
(
)
1
(
)
(
2 1 1 1 1 2 1 1 1 1 k k k k k k k k k k k kP
P
λ
λ
λ
λ
λ
λ
λ
λ
λ
λ
λ
λ
λ
λ
(2.6)olur.
[
−
4
,
−
2
]
aralığındaP
k+1(
λ
)
=
0
denkleminin köklerini bulunsun.[
−
4
,
−
λ
1(k)]
aralığının uçlarındaP
k+1(
λ
)
’nın işareti incelenirse, (2.5)’ den0
2
)
4
(
1−
=
+
>
+k
P
k olur. (2.2) formülünden ;)
(
)
(
)
2
(
)
(
1( ) 1( ) 1( ) 1 1( ) 1 k k k k k k kP
P
P
+λ
=
−
−
λ
λ
−
−λ
olup, burada (2.6) kullanırsa;
0
)
(
1(k)=
kP
λ
veP
k−1(
λ
1(k))
>
0
olur ki, bu0
)
(
1( ) 1<
+ k kP
λ
olmasını gerektirir. Bu ise Cauchy teoreminden;
∃
λ
1(k+1)∈
(
−
4
,
λ
1(k))
∋
P
k+1(
λ
1(k+1))
=
0
olması demektir. Şimdi[
λ
i(k−1),
λ
i+1(k)]
=
−
2
1
,...,
2
,
1
k
i
parçalarının uçlarındaP
k+1(
λ
1)
nın işaretleri bakılmalıdır. (2.2) formülünden;0
)
(
)
(
)
2
(
)
(
( 1) ( 1) ( 1) 1 ( 1) 1 −=
−
−
− −−
− −=
+ i k i k k i k k i k kP
P
P
λ
λ
λ
λ
dır. (2.6) den dolayı ;P
k+1(
λ
i(k−1))
=
(
−
2
−
λ
i(k−1))
P
k(
λ
i(k−1))
(2.7) ) 6 . 2( daki ikinci eşitlikten dolayı
P
k(
λ
i(k−1))
nin işareti ,(
−
1
)
k.(
−
1
)
k−i=
(
−
1
)
i olur. Bu sonuç veλ
i(k−1)∈
(
−
4
,
−
2
)
olduğu (2.7) de göz önüne alınırsa;P
k+1(
λ
i(k−1))
nin işaretii
)
1
(−
olur. (2.6) ve(
2
.
2
)
ye göre
P
k+1(
λ
i+1(k))
=
−
P
k−1(
λ
i+1(k))
(2.8) dır. (2.6) ’e göre deP
k−1(
λ
i+1(k))
nın işareti(
−
1
)
k−1(
−
1
)
l−1−i=
(
−
1
)
i olur. Bu (2.8)’de göz önüne alınırsa;P
k+1(
λ
i+1(k))
nın işareti(
−
1
)
i+1 olur.Yine Cauchy teoreminden;
(
( ))
1 ) 1 ( ) 1 ( 1k i k,
i k i+ +∈
− +∃
λ
λ
λ
=
−
2
1
,...,
2
,
1
k
i
∋
P
k+1(
λ
i+1(k+1))
=
0
olur. Böylelikle (-4,-2) aralığında
P
k+1(
λ
)
nın toplam olarak1
2
1 +
+
k
tane kökü olduğunu bulunmuş olur. k ’nın tek veya çift sayı oluşuna göre köklerin sayısına bakılırsa;
I) Eğer k tek ise; (-4,-2) aralığında
P
k+1(
λ
)
nın köklerinin sayısı2
1
+
k
olur. Simetri özelliğinden (-4,0) aralığında köklerin toplam sayısı (k+1) dir.
II) Eğer k çift olursa, (-4,-2) aralığında
P
k+1(
λ
)
nın köklerinin sayısı2
k
olur.
Simetri özelliğinden (-2,0) aralığındaki köklerin sayısı
2
k
olur.2
−
=
λ
;P
k+1(
λ
)
nın kökü olduğundan (-4,0) aralığında köklerin toplam sayısı (k+1) bulunur.Teoremin ispatından anlaşılır ki
P
n(
λ
)
nın kökleriP
k+1(
λ
)
nın kökleri arasına yerleşir.Teorem 2.5.
)
(
λ
nP
polinomunun kökleri için,0
4
) ( ) ( 1
→
−
→
∞ → ∞ → n n n n nλ
λ
, limitler doğrudur.Bu teoremin ispatı aşağıda ispatı yapılacak olan lemmaların bir sonucudur. Ayrıca bu lemmalarda Teorem 2.5’ de söylenen limitlerin yakınsama hızları da incelenecektir.
İlk olarak;
λ
n(n) nin sıfıra yakınsama hızına bakalım. Öncelikle genel terimleri ;1
2
1
2
1
+
−
=
n
n
nµ
ve1
2
2
−
=
n
n
nχ
şeklinde olan dizileri göz önüne alalım∃m
0>
1
seçelim ki[
]
0
4
,
4
µ
mλ
∈
−
−
+
∀
için 0)
2
(
)
1
)(
3
(
mχ
λ
λ
λ
>
−
−
+
+
(2.9)olsun. Bu
m
0’ların varlığı aşağıdaki limitlerden açıktır.2
3
)
2
(
)
3
)(
1
(
lim
=
−
+
−
+
−
∞ → n n n nµ
µ
µ
1
lim
n n→∞χ
=
Lemma 2.6.[
−
4
,
−
4
+
µ
n+m0−1]
∈
λ
∀
için, 1(
)
1(
)
0λ
χ
λ
n m n nP
P
+>
+ − eşitsizliği vardır. İspat: İspat indüksiyon metoduyla yapılırsa;n =1 için,
P
1(
λ
)
=
−
2
−
λ
P
2(
λ
)
=
(
2
+
λ
)
2−
1
olur. (2.9) eşitsizliğinden[
4
,
4
]
0 mµ
λ
∈
−
−
+
∀
için , 2(
)
1(
)
0λ
λ
x
P
P
>
m olur. Lemmanınhükmünün n=k için doğru kabul edip n=k+1 için bakılırsa;
]
4
,
4
[
1 0− ++
−
−
∈
∀
λ
µ
k m için,P
k+1(
λ
)
>
χ
k+m0−1P
k(
λ
)
(2.10) olup,∀
λ
∈
[
−
4
,
−
4
+
µ
k+m0]
için (2.2) formülünden)
(
)
2
(
)
(
1 2λ
λ
k kλ
kP
P
P
+=
−
−
+−
(2.11) 1 +>
n nµ
µ
olduğundan (2.10)-(2.11) den ( ) 2 1 1( ) 1 2 0λ
λ
λ
+ − + + − − − > k m k k P x P 2 1 1( ) 1 1 0 0λ
µ
+ − + − + − − ≥ k m k m k P x
)
(
)
1
(
2
1
)
1
(
2
1
)
1
(
2
1
)
1
(
2
1
2
1 0 0 0 0λ
+
−
+
−
−
+
−
+
−
+
+
−
+
−
=
P
km
k
m
k
m
k
m
k
)
(
)
(
1
)
(
2
)
(
2
)
(
1
)
1
(
2
1
1
1 1 0 0 1 0 0λ
λ
λ
+ + + ++
−
=
+
=
+
−
+
+
=
kP
kx
k mP
km
k
m
k
P
m
k
dır. Lemma 2.7.(
)
4
1
lim
( ) 2≤
−
∞ →n
n n nλ
limiti doğrudur. İspat:[
−
4
,
−
4
+
+ 0−2]
∈
∀
λ
µ
n m için,0
1
)
2
(
2
1
)
2
(
2
1
2
2
2
)
(
0 0 2 1=
−
−
≥
−
+ 0−=
−
n
+
m
−
+
n
+
m
−
+
>
P
λ
λ
µ
n mve
χ
n>
1
olduğundan Lemma 2.6.’ ya göre0
)
(
...
)
(
)
(
λ
>
P
−1λ
>
>
P
1λ
>
P
n nolur. Buradan anlaşılıyor ki
[
4
,
4
2]
0−
+
+
−
−
µ
n m aralığındaP
n(
λ
)
nın kökü yoktur.Sonuç 2.2’ den anlaşılır ki,
[
−
µ
n+m0−2,
0
]
parçasında daP
n(
λ
)
nın kökü yoktur. Buradan da 2 ) ( 0− +−
<
n m n nµ
λ
elde edilir. Böylelikle
(
)
(
)
4
1
lim
lim
( ) 2 2 2 0=
−
−
≤
+ − ∞ → ∞ →n
n n mn
n n nλ
µ
olur ∴Teorem 2.5 İçin İspat Teorem 2.4.’e göre
4
...
) 2 ( 1 ) 1 ( 1 ) ( 1n>
λ
n+>
λ
n+>
>
−
λ
yani;
{ }
λ
1(n) dizisi monoton azalan ve alttan sınırlı olduğundan bir noktaya yakınsar. Yani;4
) ( 1
→
→∞
α
≥
−
λ
n n , Benzer olarak0
...
) 1 ( 11 ) ( ) 1 ( 1 −<
<
+ +<
<
− n n n n n nλ
λ
λ
yani;
{ }
λ
n(n) dizisi monoton artan ve üstten sınırlı olduğundan bir noktaya yakınsar. Başka bir değişleλ
n(n)
→
n→
∞β
≤
0
olur. Ancak Lemma 2.7’ den β= 0 bulunur. Aynı zamanda) ( 1 ) (n
4
n nλ
λ
=
−
−
olduğundan α= -4 dür. Böylelikle sadece teoremin ispatı yapılmakla kalınmayıp, köklerin yakınsama hızları da bulunmuş olur. Yani,{ }
λ
1(n) ve{ }
λ
n(n) dizilerinin kendi limitlerine hangi hızla yaklaştığı görülür.Jacobi matrisinin elemanlarına bağlı olarak özdeğerleri hakkında bilgi edinilmesine imkan tanımıştır.
3. Kaynaklar
1. Alan Jeffrey, “Linear Algebra And Ordinary Differential Equations” CRC press, Inc., Boca Raton Ann Arbor. London, Tokyo,( 1993).
2. Kurosh, “Higher Algebra”, Mır Publıshers, Moscow, (1975).
3. John T. Moore, “Elements Of Linear Algebra And Matrix Theory”, New York, (1968). 4. O. Özkan, “İkinci Mertebeden Lineer Hiperbolik Denklemler Üzerine Bazı Karışık
Problemler”, Yüksek Lisans Tezi, S.Ü. Fen Bilimleri Enst., Konya, (1999).
5. Hochstandt H., “On Construction Of A Jacobi Matrices”, Lin. Alg. Appl., 8, 435-446, (1974).
6. Hald O., “Inverse Eigenvalue Problems For Jacobi Matrices”, Lin. Alg. Appl.,14, 63-85, (1976).
7. M. Marcus, H. Minc, “A Survey Of Matrix Theory And Matrix İnequalities”, Dover Publications, NewYork, 166-167, (1964) .
8. Courant-Hilbert, “Methods of Mathematical Physics”, Interscıence Publıshers, Inc., New York, (1953).
9. B. Aliev and A. Kh. Khanmamedov, “Energy Estimates for Solutions of the Mixed
Problem for Lineer second-order Hperbolic Equations”, Mathematical Notes, vol. 59,
No.4, (1996).
10. S. G. Kreın, “Linear Differential equations in Banach spaces (Russian )”. Ed. Nauka, Moscow, (1969).
11. J.-L. Lions, E. Magenes, “Problemes aux limites nonhomogenes et applications”, vol.1, Dunod, Paris, (1968).