• Sonuç bulunamadı

Toeplitz ve Hankel matrislerin en küçük singüler değerleri için alt sınırlar

N/A
N/A
Protected

Academic year: 2021

Share "Toeplitz ve Hankel matrislerin en küçük singüler değerleri için alt sınırlar"

Copied!
48
0
0

Yükleniyor.... (view fulltext now)

Tam metin

(1)

FEN BİLİMLERİ ENSTİTÜSÜ

TOEPLITZ VE HANKEL MATRİSLERİN

EN KÜÇÜK SİNGÜLER DEĞERLERİ İÇİN ALT SINIRLAR Mesut KIRICI

YÜKSEK LİSANS TEZİ

MATEMATİK ANABİLİM DALI

(2)
(3)

İÇİNDEKİLER

1. GİRİŞ……….………..…...1

2. TOEPLITZ VE HANKEL MATRİSLERİ……… 6

3. MATRİS ÇARPIMLARININ SİNGÜLER DEĞERİ……….... 8

4. MATRİSLERİN SİNGÜLER DEĞERLERİ İÇİN TAHMİNLER…………. 15

5. MATRİSLERİN EN KÜÇÜK SİNGÜLER DEĞERLERİ İÇİN ALT SINIRLAR

5.1. Gersgorin Tipi Alt Sınırlar……….25 5.2. Ostrowski Tipi Alt Sınırlar………....26

5.3. Brauer Tipi Alt Sınırlar………....27

5.4. Gudkov Tipi Alt Sınırlar………...31 6. UYGULAMALAR………....33 7. SONUÇ VE ÖNERİLER………..41 8. KAYNAKLAR………..42

(4)

ÖZET

YÜKSEK LİSANS TEZİ

TOEPLITZ VE HANKEL MATRİSLERİN EN KÜÇÜK SİNGÜLER DEĞERLERİ İÇİN ALT SINIRLAR

Mesut KIRICI

Selçuk Üniversitesi Fen Bilimleri Enstitüsü Matematik Anabilim Dalı

Danışman: Doç. Dr. Ramazan TÜRKMEN 2010, 43 Sayfa

Jüri: Prof. Dr. Durmuş BOZKURT : Doç. Dr. Ramazan TÜRKMEN

: Doç. Dr. Süleyman SOLAK

Bu çalışmamızda, matris ve matris çarpımlarının Singüler değerleri, matrislerin Singüler değerleri için tahminler ve matrislerin Singüler değerleri için alt sınır belirlemeyi amaç edinen Gersgorin, Ostrowski, Brauer ve Gudkov tipi özel yöntemlerin Cauchy-Toeplitz ve Cauchy-Hankel tipi matrislerde E. E. Tryteskinov

2 1 

g ve =1 özel durumu için en küçük singüler değerlerlerinde alt sınır bulma

yöntemleri uygulanmıştır. Cauchy-Toeplitz ve Cauchy-Hankel tipi matrislerde

h

2 1 

g , =1 ve =3 özel durumu için en yaklaşık sonucu veren yöntem tespit

edilmiştir.

h n

Anahtar Kelimeler: Cauchy-Hankel Matris, Cauchy-Toeplitz Matris, Singüler Değer, En Küçük Singüler Değer İçin Alt Sınır

(5)

ABSTRACT

THE POST GRADUATE THESIS

THE LOWER BOUND FOR SMALLEST SINGULAR VALUES OF TOEPLITZ AND HANKEL MATRICES

Mesut KIRICI

Selçuk Üniversity Graduate School of Natural and Aplied Science Department of Mathematics

Adviser: Assistant Prof. Dr. Ramazan TÜRKMEN 2010, 43 Page

Jury: Prof. Dr. Durmuş BOZKURT

: Assoc. Prof. Dr. Ramazan TÜRKMEN

: Assoc. Prof. Dr. Süleyman SOLAK

In this study, matrix and matrices multiplication’s Singular values, the guesses for the Singular values of matrices and the special methods of Gersgorin, Ostrowski,Brauer and Gudkow which aimed to determine the lower bounds for singular values of matrices are applied on Couchy-Toeplitz and Couchy-Hankel matrices in the E. E. Tryteskinov’s g 1 2 and h=1 special condition. the method

which gave the closest result for g 12, =1 and h n3 special condition on

Cauchy-Toeplitz and Cauchy-Hankel matrices was determined.

Key words: Cauchy-Hankel Matrix, Cauchy-Toeplitz Matrix, Singular Value, Lower Bound For the Smallest Singular Value.

(6)

ÖNSÖZ

Bu çalışma, Selçuk Üniversitesi Fen-Edebiyat Fakültesi Matematik Bölümü öğretim üyesi Doç. Dr. Ramazan TÜRKMEN yönetiminde yapılarak, Selçuk Üniversitesi Fen Bilimleri Enstitüsü’ ne Yüksek Lisans tezi olarak sunulmuştur. Bu çalışmamızda matrislerin Singüler değerleri için alt sınır belirlemeyi amaç edinen Gersgorin, Ostrowski, Brauer ve Gudkov tipi özel yöntemlerin Cauchy-Toeplitz ve Cauchy-Hankel matrisleri üzerinde uygulamaların çalışması yapılmıştır.

Bu amaçla öncelikli olarak, Singüler değerin tanımı, Toeplitz, Hankel, Cauchy-Toeplitz ve Cauchy-Hankel matrislerin tanımı ve özellikleri, matrislerin en küçük singüler değerleri için alt sınır bulma yöntemleri olan Gersgorin, Ostrowski, Brauer ve Gudkov özel metodların çıkarımları verilmiştir.

Son olarak ise Cauchy-Toeplitz ve Cauchy-Hankel tipi özel E. E. Trytyshnikov 2

1 

g ve özel durumu için n=3 şartında Gersgorin, Ostrowski, Brauer ve

Gudkov en küçük singüler alt değer bulma özel metodlarının uygulaması yapılmıştır. 1

h

Bu çalışmanın hazırlanmasında hiçbir zaman yardımlarını esirgemeyen tez danışmanım Sayın Doç. Dr. Ramazan TÜRKMEN’ e teşekkürü borç bilirim.

Mesut KIRICI KONYA, 2010

(7)

1. GİRİŞ

Günümüzde elemanter matris cebiri, matematik biliminin yanında elektrik ve bilgisayar mühendisliği gibi çeşitli teknik alanlar içinde kullanılan matematiksel bilginin ayrılmaz parçası haline gelmiştir. İngiliz matematikçi Sylvester, 1850 yılında matris kavramını kullanmıştır. Yine, İngiliz matematikçi Cayley, 1858 yılında Matris Cebirinin modern esaslarını kurmuştur. Daha sonraları ise Alman matematikçi Frobenius, matrislerle ilgili yepyeni kavramlar ve teoremler üzerinde durmuşlardır.

Tanım 1. 1. m,n Zve 1im , 1 jn olmak üzere

 

i,j ikililerin cümlesi M olsun. herhangi bir cisim ve F

 

i j f

 

i j aij F M f    , , : F aij  olmak üzere

mn m m n n

a

a

a

a

a

a

a

a

a

n m

.

.

.

.

.

.

.

2 1 2 22 21 1 12 11

şeklinde düzenlenen tabloya matris adı verilir.

Bir matrisi oluşturan değerlere matrisin elemanları denir. Burada matrisin satır (yatay) sayısını de matrisin sütun (dikey) sayısını belirtmektedir. Genel olarak bir matrisin gösterimi m n ij m n A a   

  şekliyle ifade edilir. i ve j indisleri matris elemanının konumunu gösterir. Örneğin; aija23 elemanı matrisin 2. satır ve 3. sütunundaki elemanını göstermektedir.

Eğer matrisin satır sayısı sütun sayısına eşitse, yani mn ise bu matrise kare matris adı verilir.

(8)

A matrisinin kare matris olması durumunda elemanları esas köşegen üzerinde bulunurlar.

nn a a a a11, 22, 33,..., Tanım 1. 2              n A

. . 0 0 . 0 0 . 0 0 . 0 3 2 1

tipindeki nnmatrisine köşegen matris denir ve Aköş(1,2,...,n) ile gösterilir. Tanım 1. 3. Esas köşegeni üzerinde elemanları 1 ve diğer bütün elemanları 0 olan köşegen matrise birim matris denir ve veya In I şekliyle gösterilir.

Tanım 1.4. A   aij nxn, nntipinde bir matris olsun. A.B .BAI olacak şekilde tipinde

n

nB matrisi var ise, 'B ye A matrisinin tersi denir ve B A1 şeklinde

gösterilir.

Tanım 1. 5. A   aij mxn

nxm

tipinde bir matris olsun. A matrisinin satır ve sütunlarının yer değiştirilmesiyle oluşan ve ile gösterilen matrise matrisinin transpozesi denir. T A AT A   aji şekliyle gösterilir.

Tanım 1.6. A, n-kare reel matris olmak üzere AT  ise A matrisine simetrik A matris (her i, j için aijaji ), eğer AT A ise A matrisine ters simetrik matris denir. ( her i, j için aij aji)

Tanım 1.7. . A, n-kare kompleks matris olsun. , A matrisinin elemanlarının kompleks eşleniğini göstermek üzere ise A matrisine hermityen matris denir. _ A AA AH (_)T

Tanım 1. 8. T : V → V lineer dönüşümü verilsin. x ∈ V , olan sıfırdan farklı bir x vektörü için T(x) = λ.x eşitliğini sağlayan bir λ sayısı varsa, λ sayısına T dönüşümünün özdeğeri, x vektörüne de λ özdeğerine karşılık gelen özvektörü denir.

(9)

Tanım 1.9. tipindeki herhangi bir matrisinin kendisi ile eşlenik transpozunun çarpımının öz değerinin kare köküne matrisinin singüler değeri denir. n mA A Örneğin;         4 1 2 3 5 1 A

matrisinin singüler değerlerini bulmaya çalışalım.

            4 3 1 5 2 1 T A olur ki buradan;       1 4 2 3 5 1 . =            4 3 1 5 2 1       21 9 9 35  T A A.

o halde son elde edilen matrisin öz değerini bulmaya çalışalım,            21 9 9 35 ) . det(   I AAT = (35).(21)81 = 2 56654 buradan devam edecek olursak,

654 56 2  =0 denkleminin köklerini, 130 28 1    ve 2 28 130 olarak buluruz.

O halde A matrisinin singüler değerlerini; 130

28

1  

s ve s2  28 130 olarak elde etmiş oluruz. Tanım 1. 10. V , cismi üzerinde vektör uzayı olsun. Buna göre; fonksiyonu; F I:VVF ) i xV için (x,x)0,(x,x)0 x0, V y x F a ii)  , ,  için (ax,y)a(x,y),

(10)

V y x F a iii)  , ,  (x,y)(y,x),(x,ay)a(x,y), ( ,a ' nın eşleniği) a V z y x iiii) , ,  için (xy,z)(x,z)(y,z),

şartlarını sağlıyorsa bu fonksiyona iç çarpım fonksiyonu denir.

Tanım 1. 11. Bir vektör uzayı üzerinde iç çarpım tanımlanmışsa bu uzaya iç çarpım uzayı denir. Ayrıca herhangi bir vektörü için uzunluk (veya norm) v v olarak gösterilir.

Eğer x y,  n ise bu halde (1in) olmak üzere ve olarak alabiliriz. Bu durumda

) ,..., , (x1 x2 xn x ) (

yy1,y2,...,yn x ve y vektörlerinin iç çarpımı

1 1 2 2 1 ( , ) ... n n n i i i x y x y x y x y x y      

şeklinde tanımlanır.

VEKTÖR VE MATRİS NORMLARI 1. 2.1 Normlar

üzerinde tanımlanan mutlak değer fonksiyonu ile sayıların büyüklükleri dizilerin yakınsaklığı, fonksiyonların sürekliliği, limitleri ve verilen bir reel sayı için bu sayıya en yakın asal ve tamsayıyı bulma gibi yaklaşım problemleri çözülebilir. Aynı şeyler bir vektör uzayı üzerinde tanımlanan norm için de geçerlidir. V reel vektör uzayı olmak üzere bu uzayda tanımlanan bir norm ile vektör normlarını karşılaştırabiliriz. Vektör dizilerinin yakınsaklığı irdelenebilir, dönüşümlerin sürekliliği ve limitleri çalışılabilir. Ayrıca, normlar singüler değer ayrışımında probleminin analizinde önemli rol oynar. Bu bölümde öncelikle vektör normları üzerinde daha sonra da matris normları üzerinde duracağız.

b

Ax

Genel olarak norm, . sembolü ile gösterilmektedir.

Tanım 1. 12. F reel ya da kompleks sayılar cismi ve V, F cismi üzerinde tanımlanmış vektör uzayı olmak üzere;

(11)

şeklinde ifade edilen ve

i) Her vV için, a) v0ise v 0 dır

b) =0 olması için gerek ve şart v v 0 olmasıdır, ii) aFve vV için, ava v dir,

iii)u,vV için uvuv dir

aksiyomlarını sağlayan . dönüşümüne, vektör normu denir.

Herhangi bir x vektörünün pozitif bir sayıya dönüştürülmesi işlemine norm denir.

1. 2. 2 Matris Normları

Tanım 1.13. F reel yada kompleks sayılar cismi ve ; bileşenleri F cisminin elemanları olan -kare matrisler kümesi olmak üzere;

) (F Mn n . :M Fn( )   {0} A A

şeklinde ifade edilen ve i) AMn(F) için, a)A0 ise A >0 dır.

b) A =0 olması için gerek ve yeter şart =0 olmasıdır, A ii) aFve AMn(F)için, aAa A dır,

iii) için, A,BMn(F) ABAB dir, iv) A,BMn(F) için, ABA B dir,

(12)

2. TOEPLITZ VE HANKEL MATRİSLERİ

Tanım 2. 1. n1 olmak üzere

     1 0 , 1( , ) n j i j i j i n x x h x x H

kuadratik formuna Hankel formu denir. Bu forma uyan matrise de Hankel matrisi denir ve 1 0 , 1 ( )      i j nij n h H olarak gösterilir.

Bir Hankel matrisinin açık gösterimi,

                                  2 2 3 2 1 1 3 2 4 2 1 1 1 2 1 1 1 0 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . n n n n n n n n n n n n n h h h h h h h h h h h h h h h h H

şeklindedir. Buradan görüldüğü gibi Hankel matrisleri simetriktir. Ayrıca sonsuz mertebeden bir Hankel matrisi de

    (hi j)i,j 0 H olarak tanımlanır.

Tanım 2. 2. n1 ve ti,j ler kompleks sayılar olmak üzere

1 1 , 0 ( , ) n n i i j T x xt x x  

j i j 0

kuadratik formuna Toeplitz formu denir. Bu forma tekabül eden

1

1 ( ),

n n i j i j

Tt

(13)

Bir Toeplitz matrisi açık olarak,                                      0 1 2 1 1 0 1 2 3 0 1 1 2 1 0 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . t t t t t t t t t t t t t t t t T n n n n n n n n n

şeklinde yazılır. Buradan görüldüğü üzere bir Toeplitz matrisinin elemanları esas köşegene paralel köşegenler boyunca aynıdır. Dolayısıyla bir Toeplitz matrisini, matrisin ilk satır vektörü ile ilk sütun vektörü temsil eder diyebiliriz.

Tanım 2. 3. g, h herhangi sayılar ve g/h Z  olmak üzere; n j i n jh ih g T 1 , 1          

matrisine Cauchy-Toeplitz matrisi denir.

Bu matris hem Cauchy hem Toeplitz matrisidir. Tanım 2. 4. g,h herhangi sayılar ve g/h Z

n j i n g ih jh H 1 , 1          

matrisine Cauchy-Hankel matrisi denir. Burada matrisi, hem Cauchy hem de Hankel matrisidir.

n H

(14)

3. MATRİS ÇARPIMLARININ SİNGÜLER DEĞERLERİ

n m , kompleks matrislerin uzayı olsun. H hermityen matrisinin n n

özdeğerleri 1( ) ...H  n( )H ve Cn m nin singüler değerleri

1( )C 2( )C ... n( )C     olmak üzere 12(CC) i  ) (Ci  , ( ) dir. ve n i1,2,..., p m

C qmax

 

p,mn olmak üzere, q1( )C  ... n( ) 0'C  dır. olacak şekilde tamsayılarını alalım ve aşağıdaki Lemmaları verelim. ni i1 2 1  ... i ,...,1 in Lemma 3.1: W ,...,1 Wk, kümeleri n k V V ,...,1 boyVt   n it 1, , t i t

boyWi 1i1 ...ik ve V1  ...Vk olmak üzere nin alt uzayı olsun. Bu takdirde bir ortonormal baza sahip bir alt uzay ve

' n

x1,...,xk

nWxt ve Wt

, olmak üzere diğer ortonormal bazı t

t V

zt1,...,k

z 1, ,z dir. [Wang ve Xi, k

1997]

Lemma 3.2: G H,  , pozitif yarı tanımlı hermityen matrisler ve n n

olsun. Buradan n i ik   ... 1 1

   k t k t i i i GH G H 1 1 ) ( ) ( ) (    ve ) ( ) ( ) ( 1 1 1 H G GH n t k t k i i it t   

 

 

olup kn için eşitlik vardır. p n A ,  Bn m olsun. Bu takdirde 1 2 1 1 ( ) ( t t k k i i t t AB ABB A    )   

(15)

1 2 1 ( ) t k i t A ABB          

1 2 1 1 ( ( t k i n t t AA BB             

) )

   k t i n t B A t 1 1 ) ( ) (   dir. Buradan, 1 1 1 ( ) ( ). ( t t k k i i n t t t AB A B        

(3.1)   k t t n it A B 1 1( ) ) (  

elde edilir. [Wang ve Xi, 1997]

Lemma 3.3: Ap n, Bn m ve i i n k     ... 1 1 olsun. Bu takdirde,

(3.2)      k t i n i k t t AB t A t B 1 1 1 ) ( ) ( ) (    ve

(3.3)      k t k t i n i t AB t A t B 1 1 1( ) ) ( ) (    dir. Sonuç 3.1: A , p nBn m, i i n k     ... 1 1 ve 0 r   olsun. Bu takdirde

     k t t n r r i k t r it AB t A B 1 1 1 ) ( ) ( ) (   

     k t i n r r i k t r it AB t A t B 1 1 1 ) ( ) ( ) (   

(3.4)      k t i n r i k t r it AB t A t B 1 1 1 ) ( ) ( ) (   

(16)

İspat:. olmak üzere A matrisinin özdeğerleri, özdeğerlerine karşılık gelen ortonormal özvektörleri olsun. Yani,

n n A( ),..., ( ) 1 AA n AA  n u u u1, 2,..., i i i i i AA u Au u AA  ( )  2( ) , i1,2,...,n dir. n ile t

R  R t

u1,u2,...uit

vektörleri tarafından gerilen bir alt uzay ve de St

uit,uit1,...,un

, t 1,2,...,k

vektörleri tarafından gerilen alt uzayı verilsin. 1, ,...2 n için k x x x 

n k k k x x x C X1, 2,...,   , k k k diag C I  (1,1,...,1)  olur. Lemma 3.4: An m olsun ve , , 1, 2,..., t t t t yR zS tk n k k olmak üzere; , ) k y ,..., , ( 1 2 k y y YZk ( , ,..., )z z1 2 z  olsun. ytRt,ztSt,t 1,2,...,k olmak üzere;

i) Eğer YkYkIk ise buradan (Y A) (A),

t

i k

t

   t 1,2,...,k' dır, ii) Eğer ZkZkIk ise buradan (Z A) (A),

t i k t     t 1,2,...,k' dır. [Wang ve Xi, 1997]

İspat: Sabit bir t (1tk) için Rt yi geren

y1,y2,...yt,vt1,...,vit

A ve nın aynı singüler değere sahip olduğunu görmek kolaydır.

1 2

( , ,..., ) '

t

i

u u uA

y y1, ,... ,2 y vt t1,...,v it

Rt yi geren ortonormal bir baz olacak şekilde vt1,...,vitRt

seçelim. Bu takdirde BBq q hermityen matrisinin özdeğerleri ve alt

matrisinden p CC p ) ( ) (B i C i    ve qi1(B)pi1(C) (i 1,2,...,p) olup t(YkA)t

(y1,y2,...,yt)A

it

(y1,y2,...yt,vt1,...,vit)A

i

u u ut A

t   ( 1, 2,... ) it(A)

olur. Benzer şekilde

t(ZkA)

(zt,...,yk)A

1

(17)

ui un A

t  1 ( ,... ) = (A) t i  dir. Teorem 3.1: ACnm, ve i i n k     ... 1 1 ), ,..., k

ve her bir değişken aralığında artan olacak şekilde

0,

(1

k değişkenin herhangi bir fonksiyonu olsun. Bu taktirde

k I k X k X t W t k t t t k x k k k i W i i A X A X         max min ( ),...,,..., ( ) ),... ( 1 dim ,..., 1 1      (3.5)

dir. [Wang ve Xi, 1997]

İspat: t 1,2,...,k boyWtit ile Cn in herhangi bir alt uzayı olsun. t W boyStnit 1,t 1,2,...,k ve S1S2 ...Sk' dır. ' de Cn ve t t W xzt  , St olmak üzere

ve k ,..., 2 , 1 x ,1

tx2,...,xk

z1,z2,...,zk

ortonormal tabanlara sahip bir Cn alt uzayı vardır. Lemma 3.4 (ii) ve  artan bir fonksiyon olmadığından

1(XkA),...,,...,k(Xk)

 

1(ZkA),...,,...,k(Zk)

 ( ( ),..., ( )) 1 A ik A i    

olur. Dolayısıyla (3.5)' in sol tarafının sağ taraftan daha küçük olmadığı görülür. Diğer taraftan WtRt (t 1,2,...,k) ve

y1,y2,...,yk

, olacak şekilde herhangi bir ortonormal bir küme olsun. Lemma 3.4 (i)’ den

t t R y

1(Y Ak ),..., k(Y Ak )

i1( ),...,A ik( )A

       

olur. Buradan (3.5) in sol tarafının sağ taraftan daha büyük olmadığı görülür. '

Teorem 3. 2. A,BCnxn, i i n k     , .... 1 ve 0rR olmak üzere 1 1 (3. 6) 1 1 ( ) ( ) ( t k k r r r i n t n t t t AB A B        

)

(18)

İspat:

k ve t r t k 1 1,...., ) (  

WtRt' yi alalım. Teorem 3. 1, Özellik 3. 1 ve Lemma 3. 4 (i) yi kullanırsak '

     k t k r t I Y Yy R k t r i AB Y AB k k k t t 1 * 1 ) * ( min ) (  

     k t r t n k r t I Y Yy R B A Y k k k t 1 1 * ) ( ) * ( min  

   k t r t n r it A B 1 1( ) ) (   elde edilir .

Sonuç 3.2. , pozitif yarı tanımlı hermityen matrisler ve

ve ise nxn C H G,  n0 r  i i   k  ... 1 1

(3. 7)     k t r t r t n r i r it GH t G H 1 1 1 ) ( ) ( ) (   

     k t r t r t n r i r it GH t G t H 1 1 1( ) ) ( ) (    dır. [Wang ve Xi, 1997]

İspat: A,B' nin yer değiştirmesi ve r üzerinde aşağıdaki sonuçlar G12, H 12 ve 2r ile

(3. 4) ve (3. 6)' da verildi. (3. 7)’ den, A,BCnxn için Teorem 3. 2 sağlanır.

Teorem 3. 3. m3, 1i1 ....ikn ve 1,..., nxn olsun. Bu takdirde, m A AC i) 1 ( ... ) ( ) ( )... ( ), 1 2 1 1

  k t m r t k t r t r i m r it A At AAA  ii )

tk r m t r t r i k t m r it A A 1 t A1 A2 A 1 ( 1... )  ( ) ( ).... ( ), iii) ( ... ) ( ) ( )... ( ), 1 1 1 2 1 1 1

    k t m r t n k t r t n r i m r it A At AAA  iv) 1 ( ... ) ( ) ( )... 1( ), 1 2 1 1 1

     k t m r t n k t r t n r i m r it A At At AA

(19)

İspat: (i) ve (ii) aşikârdır. (iii) ve (iv), Sonuç 3. 1 den m üzerinden tümevarımla, '

    k t k t m r t n m r i m r it A A t A A A 1 1 1 1 1 1... ) ( ... ) ( ) (    =      

k t m r it A A 1 1... ) ( 

  k t m r t n A 1 1( )  ( ) ( )... ( ) . ( ) 1 1 1 1 1 2 1 1

             k t m r t n m r t n k t r t n r it AAAA  = ( ) ( )... 1( ), 1 2 1 1 m r t n k t r t n r it A A  A  

   1 ( ... ) ( ) ( )... 1( ) 1 3 1 2 1 1

     k t m r t n k t r t n r i m r it A At AAt AA1 ( ... ) ( ) ( )... 1( ) 1 2 1 1 1

     k t m r t n k t r t n r i m r it A At At AA  dir.

Not: (3. 6) ve (3. 7) den ' m3 için (iii) ve (iv) doğal konjektürü sağlanır.

Örnek 3. 1. , 1 0 0 0 1 0 0 0 3 1            A ,            1 0 0 0 3 0 0 0 1 2 A            3 0 0 0 1 0 0 0 1 3 A ve k=n, r=1 alalım. Buradan ( ) ( ) ( ) 1( 3), 1 1 2 1 1 3 2 1A A A A A A n t k t k t t n t t    

 

   (913) dır. Sonuç 3. 3. A A A Cnxn m it ik n, m , 3,1 ....  .... 2 1 0 r   olsun. Buradan; i ) 1 ( ... ) ( ) ( ) ... 1( ), 1 2 1 1 1

     k t m r t n r n k t r t n r i m r it A At At A   A  ii) 1 ( ... ) ( ) ( ) ... 1( ) 1 2 1 1 1

     k t m r t n r n k t r ti n r m r it A Aikt At A   A

(20)

dir. [Wang ve Xi, 1997]

İspat: (3. 6) ve (3. 4) ' ü kullanarak herhangi A B C, nxn için ) ( ) ( ) (AB A r B n r n r i     dir.

(21)

4. MATRİSLERİN SİNGÜLER DEĞERLERİ İÇİN TAHMİNLER

Teorem 4. 1. , . mertebeden bir kompleks matris olsun.

matrisi, dan büyük olsun. Yani, ) (aij AA ' n ( ) nxn ij Bb   i,j için bij

 

aij olsun. Buradan A' nın bütün öz değerleri; { : ( ) } 1 ii n i ii B b a z C z     

(4. 1) nin içindedir.

Singüler değerlerin tahmini, pratikte önemlidir. A ' nın singüler değerleri için tahminler elde etmek maksadıyla matrisi yukardaki teoremde kullanılır. Bu yaklaşım bizi aşağıdaki teoreme götürür. [Li, 1999]

* A A Teorem 4. 2. nxn olsun. ij C a A( )

    i j j i ji i ij i a c a r , u max(0,u),si max(ri,ci) ve aiaii,i 1,2...n (4. 2)

olmak üzere ' nın bütün singüler değerleri, A

B Bi, i [(aisi) , aisi]  (4. 3) içindedir. [Li, 1999]

Teorem 4. 3. ve olacak şekilde negatif olmayan

matris olsun. nxn ij C a A( ) ( ) nxn ij Bb   0 ) ( }, , max{a a i j bijij ji  

(22)

ise A ' nın bütün singüler değerleri; 1 { : ( ) } n i ii i z z aB b       

(4. 4) birleşimi içindedir. Burada; B ' nin spektral yarıçapı ve (i1,2....,n) için aiaii

ve , pozitif reel sayılardır. [Li, 1999]

İspat: ,A' nın bir Singüler değeri olmak üzere

x A y * ,yAx (4. 5)

olacak şekilde T ve vektörleri vardır.

n x x x x( 1, 2,..., ) T n y y y y( 1, 2,..., )

Farz edelim ki; ve Teorem 4. 3 ü sağlasın. İlk olarak ı alalım. Bu durumda aşağıdaki teorem verilebilir.

0 )

( 

bij

B ' B0'

Lemma 4. 1. c,, 1 ve   olsun. Buradan;

 c max{ c, c} (4. 6) dır. [Li, 1999] İspat: i)   c ise; c c c          ii)  c ise;   c   c  c olup ispat tamamdır.

Teorem 4. 4. nxn olsun. ij C

a

(23)

     a z a s s R z z j i j n j i j i i. , } : { 1 ,

(4. 7) içerisindedir. [Li, 1999]

İspat: u (u1,u2,...,u )T,u 0 olmak üzere n

u1 0,(i) ve Bu(B)u şartını sağlayan bir u vektörü vardır.

i i i i u y y u x xˆ ,ˆ  olsun. (4. 5) den; '

     n j n j i j j ij i j j ji i i n u u x a y u u y a x 1 1 ... 2 , 1 , ˆ ˆ , ˆ ˆ   olup i j i j ji j ji i u u y a y a xˆ ˆ ˆ 1

    (4. 8) i j j j ij i ii i u u x a x a yˆ ˆ ˆ 1

    (4. 9)

dır. zi max{xˆi, yˆi} dır. }zp max{zj :1 jn ise olduğu açıktır. Genelliği bozmadan 0  p z p p p y

z  ˆ olsun. (zpyˆpxˆp olduğu zaman benzer

argümanlar kullanılır.) p p y x ˆ ˆ   alalım. (4.8) ve (4.9)` da eşitsizlikler; p j j p j jp p pp u u y a z a

ˆ     p j j p j pj p pp u u x a z a

ˆ     

(24)

olur. Buradan; p j p j jp p pp u u a z a

    (4. 10) p j p j pj p pp u u a z a

     (4. 11)

olur. ajpbpj olsun. Bu(B)u ' nın p. bileşeni;

     p j p pp p j pj j p j jpu b u B u b u a ( ) (4. 12)

dir. (4. 10) ve (4. 12) eşitsizlikleri gereği;

aˆpp (B)bpp (4. 13)

elde edilir. Benzer şekilde; apjbpj ve (4. 11)’ den

 app (B)bpp (4. 14)

elde edilir.  1 idi. (4. 6), (4. 13) ve (4. 14) ile

 ap (B)bpp (4. 15)

(25)

Son olarak B0 ve benzer durumları sağlasın.  0 keyfi sayısı için

) ( )

(  bij 

B pozitif ve diğer kabulleri sağlasın. )B( ' a yukardaki durumları

uygularsak; ) 0 ( ) ( ) ( )) ( (         B bii B bii

olur. Bu ise ispatı tamamlar.

Teorem 4. 5. Teorem 4. 3' teki kabulleri alalım. max{0, },dd d olmak üzere R

ii i B B b l ( )( ) ve Vi [(aili(B)),aili(B)] (4. 16) olsun. buradan, a) Bütün Singüler değerler

n i i dedir. i V V 1 :

b) Eğer, V aralığının n bileşeni varsa ve n-k intervallerinden bağımsız ise,V ' de A' i nın Singüler değeri vardır. [Li, 1999]

İspat: (a)' nın ispatı, Teorem 4. 3' ten görülür. (b)’ nin ispatı aşağıdaki gözleme dayanır.

(b. 1) [0,1], Ddiag(a11,a22,...,ann) , A()DC ve C=A-D olmak üzere , ) ,...,bnn B G F , (b11 b22 diag

G ()  olsun, F=B-G olmak üzere;

B B A A G B D A(0) , (0) , (1) , (1)

(26)

dir. B(),A()' a göre daha geniş bir spektral yarıçapa sahiptir. (a) dan A()' nun bütün Singüler değerleri,

n i i i i i i i a l B a l B V V 1 ) ( : ))] ( ( ))) ( ( [( ) (          aralığındadır. s t, [0,1], s t için; )) ( ( )) ( (B sB t  

dır. ye göre daha geniş bir spektral yarıçapa sahip olduğundan B(t),B(s)' )) ( ( )) ( (B s l B t lii ve Vi(s)Vi(t)Vi(1)Vi dır.

(b. 2) A' nın Singüler değerleri sürekli fonksiyonlardır. A(0) ' ın Singüler değerleri, örneğin ) ,..., 2 , 1 ( ), ,..., 2 , 1 (i n a V i n aiii  sağlanır.

(b. 3) , 0' dan 1 e değiştiği zaman, sürekli olarak, ' A() ' un i() Singüler değeri, a den ji  (1) ' e değişir.

Teorem 4.4. ün farklı bir ispatı:  , ' nın bir Singüler değeri ve (4. 7)' yi sağlasın. A , ) ,..., , ( , ) ,..., , ( 1 2 1 2 T n T n y y y y x x x x  i max{xi, yi} ve i n i p      1

max ise p 0 olduğu açıktır i p i

q

 max olsun. (4.5) ' de p . denklem;

p pp p jp j p j x a y a y    

(4.17)

   p j j jp p pp p a x a x y  (4.18)

olur. Genelliği bozmadan, pypxp olsun . Buradan; yp 0 dır.

p p y x   yazalım. (4.17) ve (4.18) eşitsizliklerinden; pp p q jp q p j p a a s       

(27)

p q p j pj q p pp a s a      

 

olur. Lemma 4.1 i uygularsak; '

p q p p s a     (4.19)

elde edilir. Eğer q 0 ise  ap ve  , (4.6) ' de verilen küme içindedir. q 0 olsun. Benzer argümanlarla (4.19) ' da

q p p q s a      (4.20) olur. Sonuçta: q p q p a s s a     

elde edilir. Böylece Teorem 4.4 ispat edilmiş olur.

Örnek 4.1.         2 1 1 6

A ise 1(A)6.1231, 2 2.1231 olur. negatif olmayan matris ve        2 1 1 6 B 2361 . 6 ) B (

 , matrisinden daha geniş bir spektral yarıçapa sahiptir. Teorem 4.2 uygulanırsa; A ' nın bütün singüler değerleri

A

0,6.2361

aralığında olur. Diğer taraftan 2361 . 6 ) ( max A   dir. Bu ise 7 ) ( 1 max AA A  

(28)

Teorem 4.2 ' yi uygulayarak ise, bütün singüler değerlerin

 

      6 1 1 6 1 B 1,3 

 

5,7 '

de olduğunu biliyoruz. Buradan; min 1' dir. Bu ise;

         2 min ) ( min i i i c r a A

ile aynıdır. Buradan

  i j i r aij ,

  i j ji a i

c dir. Şart sayısı için bir üst sınır

2 / ) ( min ) ( ) ( 12 1 min max i i i r c a A A A A K        (4.21) olarak verilir.

Sonuç olarak; A ' nın şart sayısı için üst sınır ( ) 6.2361K A  olarak bulunur. Bu ise; (4.21) ile verilen sınırdan daha iyidir.

Örnek 4.2.            3 1 0 1 10 1 0 0 6

A ise, A nın singüler değerleri 10.2137, 5.9588 ve 2.8590 dır. ' '

3649 . 10 ) (            3 1 0 1 10 1 0 1 6

B negatif olmayan matris ve p B olup A ' ya göre daha

geniş bir spektral yarıçapa sahiptir. Teorem 4.3' ten biliyoruz ki; A ' nın bütün singüler değerleri

0, 10.3649

aralığındadır. Buradan Teorem 4.3 ' ü uygulayarak

           12 1 0 1 10 1 0 1 9 1 B 4605 . 12 ) (B1p

matrisi daha geniş bir spektral yarıçapa sahip bir matris ise;

(Perron köktür) tir.

Bu durumda A ' nın bütün singüler değerleri;

2.5395, 12.4605

aralığındadır. Son olarak 4.5 i '

(29)

           10 1 0 1 10 1 0 1 10 2

B negatif olmayan matrisine uygularsak (Perron

kök) olur. Bu durumda bütün singüler değerler

4142 , 11 ) (B2p

1.5858, 4.4142

 

 4.5858,7.4142

 

 8.5858,11.4142

aralığındadır. Bu üç aralığın her biri ' nın bir singüler değerini içerir. (4.3) ' ü uygulayarak A 4891 . 11 ) ( max A   ve 2 ) ( min A  

elde edilir. Teorem 4.3 ve Teorem 4.5 ' den max(A)10.3649 ve min(A)2.5395 olur. A ' nın şart sayısı için bir üst sınır K(A)4.0815 olup, (4.21) ' den

sağlanır. Teorem 4.4 ' ü uygulayalım. (4.5) ' ten; 7446 . ) (A K 5

2.6972,3.3820

 

 5.5505,6.5858

 

 9.4142,10.4495

ve 8742 . 3 ) (AK

dir. A' nın şart sayısının daha dar bir alt sınırı elde edilir. Bizim nümerik deneylerimizden; Teorem 4.4' ün Teorem 4.3 ' ten daha iyi bir sınır verdiği görülür. Fakat bu durumda Kuadratik eşitsizlikleri çözmek için çaba sarf etmek gerekir. Örnek 4.3.         1 0 1 1 ) (aij

A olsun. A' nın en küçük singüler değeri;

2( )A 5 1 / 2 0.6180

   

dir. (4.3) veya (4.4) ten '

( )

/2 0.5 min ) ( 2 Aairici  

elde edilir. Ancak Teorem 4.4 ten ' 2(A) için 0' ın alt sınır olduğu açıktır. Teorem 4.3 ve daha geniş spektral yarıçaplı negatif olmayan matristen;

(30)

(4.22)        y x b B ij 1 1 ) ( alarak. p(B)=

x y (xy)2 4

/2 elde edilir. )' ( 2 A

nın keyfi bir alt sınırını elde etmek için;

a11p(B)x , a22p(B)y (4.23) elde edilmelidir. Fakat bu (4.23) ile çelişir; x ve y y olur. x

Teorem 4.6. n n, indirgenemez ve C(A), açık olmayan ij C

a

A( )(A)

çemberinin kümesi  olsun. Buradan A ' nın bütün özdeğerleri, ) ( : A C pi ii pi i r a z C z           

(4.24)

bölgesi içindedir. [Li, 1999] Örnek 4.4.            0 0 1 . 0 1 . 0 0 0 0 1 0

A indirgenemez bir matris olsun. Teorem 4.2, Teorem 4.3 ve

Teorem 4.5 ' i alalım. A ' nın bütün singüler değerleri

z0: z03 s1s2s3

0,0.4642

(4.25) kümesi içindedir. Ancak A' nın singüler değerleri, 1 1, 12 3 0. ' dir ve (4.25),  ' i içermez. 1

(31)

5. MATRİSLERİN EN KÜÇÜK SİNGÜLER DEĞERLERİ İÇİN ALT SINIRLAR

A(aij) biçiminde nn tipinde kompleks bir matris olsun, k 1,...,n için

  k j kj k A a P ( ) ve

  k j jk k A a Q ( ) şeklinde tanımlayalım.

A ' nın en küçük singüler değerini )n(A ile göstereceğiz. Gersgorin Teoremi kullanılarak ;

[ ( ) ( )]} 2 1 { min ) ( 1 a P A Q A A kk k k n k n    (5.1) olduğu gösterilmiştir. [4]

Burada, Gersgorin Teoreminin yerine Ostrowski ([2] , Teorem 6.4.1), Brauer ([2] , Teorem 6.4.7) ve Gudkov ([1]) metotları kullanılarak )n(A için yeni sınırlar elde edilmiştir.

Yeni sınırlar daha karmaşık ifadeler içerebilir, ancak genellikle daha büyük, daha iyi sınırlar vermektedir ve karşılıklı olarak karşılaştırılamazlar. Gersgorin Teoremi' nin Singüler olmama göstergesinin optimumluğuna rağmen bunun gibi benzer durumlarda değişimlerin olmasının hiçbir önemi yoktur.

Bu sonuçlar kare matrislerin elemanları için gösterilmiştir, fakat [3]` ün 3. 7 bölümündeki metot ise dikdörtgen matrisler içinde genelleştirilebilir.

Burada Ostrowski ve Brauer Teoremi kadar bilindik olmayan Gudkov` un Singülerlik testini inceleyeceğiz.

Teorem 5. 1. R1(A)P1(A) olsun ve i2,...,n için

      n i k ik kk k i k ik i a a A R a A R 1 1 1 ) ( ) (

(32)

tanımlansın. Eğer; i1,...,n için aiiRi(A) ise Singüler değildir. [Horn ve A Johnson, 1985]

5. 2. n(A) İçin Ostrowski Tipi Alt Sınırlar

Teorem 5. 2 :

A

' nın en küçük Singüler değeri, )n(A ise;

({4 [ ( ) ( )] } [ ( ) ( )])} 2 1 { min ) ( 2 1 2 2 1 a P A Q A P A Q A A kk k k k k n k n      (5.2)

şeklinde bir alt sınır elde edilmiştir. [Horn ve Johnson, 1985]

İspat : Öncelikli olarak, [2] deki Teorem 7. 3. 7 kullanılarak         I A A A I A n n ) ( ) ( *

singüler olduğu görülür ve buradan;

(5.3)        I A A A I A n n ) ( ) ( *                 * ) ( ) ( 0 0 A I A I A A I I n n   matrisi de singülerdir.

Ostrowski Teoremi ile matrisin singülerliliğini birleştirerek en az bir için ( ); i n i  1 aii 2 [Pi(A)n(A)][Qi(A)n(A)] (5.4) eşitsizliği elde edilir. (5.2) yardımıyla (5.4) düzenlenirse ;

2 2 ) ( ). ( ) ( )] ( ) ( [ ) ( 0n APi AQi An APi A Qi Aaii

olduğu görülür. (5.2) ile (5.1) sınırları arasında ilişki aşağıdaki gibi doğrudan görülür.

Sonuç 5. 1 : (5.2) sınırı her zaman (5.1) sınırından büyüktür. Üstelik Pr(A)Qr(A) olmak üzere için ; A

)] ( ) ( [ 2 1 )]} ( ) ( [ 2 1 { min 0 1 k n akkPk AQk AarrPr AQr A   

(33)

5. 3. n(A) İçin Brauer Tipi Alt Sınırlar

Önceki bölümde, [2] deki Teorem 6. 4. 7 (Brauer Teoremi) ile Ostrowski Teoremini değiştirirsek ;              )]) ~ ( ) ~ ( [ } )] ~ ( ) ~ ( [ ~ . ~ 4 ({ 2 1 min ) ( 2 1 2 2 , 1 a a P A P A P A P A A ii jj i j i j n j i n j i  (5.5)

olmak üzere (5.5) sınırı elde edilebilir.

* 0 0 ~ A A A

(5.5) ve (5.2) karşılaştırıldığında, önceki sınırın sağ tarafı en fazla yeni sınırın sağ tarafı kadar büyük olabilir.

Bu nedenle biz, (5.1) ' in uygun sınırları arasına (5.5) ' yi ilave etmeyeceğiz. Bununla birlikte [3] ' deki 3. 7 bölümü ve Ostrowski ve Brauer öz değer bölgesi arasındaki ilişki, en azından (5.1) ' den daha büyük ve (5.2) ile kıyaslanamayan “Brauer-Tipi” türü bir sınır belirtir.

Bu fikri doğrulamak için;

         2[Re Re (Re Re ) ( ) ( )] 1 min ) ( 2 * * , 1 a a a a P A A P A A A kk jj kk jj k j n j k n j k

ile tanımlı n(A) değeri ile A' yı ilişkilendirmek uygun olacaktır. Bir lemma ile başlayalım.

Lemma 5. 1 C(cij), nn tipinde bir hermityen matris olsun. O halde I

C

(34)

İspat :Cn(C)I hermityen olduğunu biliyoruz. Bu yüzden bütün öz değerler reeldir. Öz değerlerin negatif olmadığını göstermek için Brauer teoremini uygulayalım.( [2] , 6. 4. 11 )

Bunu takiben Cn(C)I için yukarda bahsedilen spektral özelik olan bir yeter koşul cqq n(C)0 (q1,...,n) (5.6) ve [css n(C)][crr n(C)]Ps(C)Pr(C) (s,r 1,...,n;sr) (5.7) elde edilecektir. n q 

1 olmak üzere Cn(C)I ' nın köşegen elemanını düşünerek q.

             2 ( ) 4 ( ) ( ) 1 min ] ) ( ) ( 4 ) ( [ 2 1 ) ( 2 , 1 2 P CP C c c c c P CP C c c c c C c kk ij kk jj k j n j k r q rr qq rr qq n qq j k

0 2 1 ) ( 2 1 ) ( ) ( 4 ) ( 2 1 2 cqq crr cqq crr Pq C Pr C cqq crr cqq crr

elde edilir. (5.7) yi elde ederken, ' a2  b4 0 için f() reel fonksiyonunun )

2 ,

( a aralığında kesin azalan ve f()2 ab olduğuna dikkat edelim.

Böylece herhangi bir ( 4 ) 2

1 a a2 b

 değeri için f()0 olur. Bu ifadeye dayanarak (5.7) elde edilir.

Artık, lemma yardımıyla önerilen sınırları elde edebiliriz.

Teorem 5. 3. A ` nın en küçük singüler değeri için ) ( ) (A n A n   

sağlanır. [Horn ve Johnson, 1985]

İspat: [4]' de Teorem (1) ' in ispatındaki yol izlenerek sadece BAn(A)I için

*

B

B öz değerlerinin negatif olmadığını göstereceğiz. ) ( ) ( ) (A A* A A* n n n     dikkate alınarak;

(35)

A A

I A A A A I A

A B

B * * n( )n( *) * n( *) elde edilir. Kabul edilen spektral özelik Lemma 2. 1 den görülür. '

Teorem 2. 2 ' nin gelecek sonucu n(A) için (1) ve (2) dekilerin karışımı olan bir sınırı verir.

Sonuç 5. 2. A ' nın en küçük Singüler değeri



     

   2 ( ) [ ( ) ( )][ ( ) ( ) 1 min ) ( 2 , 1 a a a a P A Q A P A Q A A kk jj kk jj k k j j j k n j k n  (5.8) dır. [Horn ve Johnson, 1985]

İspat : [3], 3. 7. 17 ' nın Sonuç 2. 1 İspatı, bu ispatla benzerdir ve bu sebepten dahil edilmemiştir.

(5.8) sınırı ve (5.1) sınırı arasındaki ilişkiyi açıklayarak bu bölümü kapatalım.

Sonuç 5. 3. (5.8) sınırı her zaman (en azından) , (5.1) sınırı kadar büyüktür. Üstelik A için;

( ) ( )

} 2 1 { min 0 1 k n akkPk AQk A    [ ( ) ( )] 2 1 A Q A P arrrr  [ ( ) ( )] 2 1 A Q A P ajjjj  (5.9) eşitsizliği herhangi bir \ için sağlanır. Önceki sınır açık değildir ve bir sonrakinden daha iyidir. [Horn ve Johnson, 1985]

} ,..., 1 { n j{r} İspat :

     2 ( ) [ ( ) ( )][ ( ) ( )] 1 min 2 , 1 akk ajj akk ajj Pk A Qk A Pj A Qj A j k n j k

( ) [ ( ) ( )][ ( ) ( )]

2 1 a a a a 2 P A Q A P A Q A t t s s tt ss tt ss        olduğunu varsayalım.

Buradan Sonuç 2. 2 yardımı ile

n(A)

( ) [ ( ) ( )][ ( ) ( )]

2 1 a a a a 2 P A Q A P A Q A t t s s tt ss tt ss       (5.10)

Referanslar

Benzer Belgeler

Variorum bu y~l, Osmanl~~ kent tarihi alan~ndaki ara~t~rmalanyla dikkati çeken ve makaleleri aranan ODTÜ Tarih bölümü ö~retim üyelerinden Suraiya Faroqhi'nin; Bekta~i

- Adli t ıp uzmanlığı, doğal olarak tüm dünyada patoloji- otop- si a ğırlıldı bir eğitim ve uzmanlık yap ısına sahiptir. Ama bizde ne tuhaft ır ki, akademisyenler de

Sera koşullarında bitki büyümesini artırıcı rizobakter uygulamalarının marulda kök yaş ağırlığına etkisi.. Çizelge 4.1 ve Şekil 4.1.’den de görülebileceği

Paksoy (2005), malzeme ihtiyaç kısıtı altında çok aşamalı bir tedarik zinciri ağı tasarımı için karma tamsayılı bir doğrusal programlama modeli

Üçüncü monomer olarak toplam mol miktarına göre %15 mol NHMAAm kullanılan hidrojelde şişme büzülme değerlerinden LCST değeri 24°C ve toplam mol miktarına göre %15

Bunlar: Devlet Başkanının halk tarafından seçilmesi, Devlet Başkanının önemli (anayasal) yetkilere sahip olması ve yürütme görevini yerine getiren bir

Bu anlamda bu çalışmada öncelikle bir ilişki türü olarak uluslararası ilişkilerin tanımlanmasına yönelik sorunlara kısaca değinilecek daha sonra disiplinin tarihi

Şirin Devrimin, Joseph Conrad'a benzeyen amcası ile rakı içtiğini, Devrimin Yale Tiyatro Okulunda okuduğunu İngiltere ve Amerika’da sah­ neye çıkmış