Sayısal Metotlar
Kaynak Kitaplar :
Numerical Analysis
Richard L.Burden, J.Douglas Faires
International Thomson Publishing Company
Numerical Methods Using Matlab John H. Mathews, Kurtis D. Fink Prentice Hall
Mühendisler için Sayısal Yöntemler
Steven C. Chapra, Raymond P.Canale
Çeviri: Hasan Heperkan, Uğur Kesgin
Literatür Yayınevi
İÇERİK
-Tek değişkenli fonksiyonların köklerinin bulunması -Interpolasyon(ara nokta tayini)
-Sayısal integral -Sayısal türev
-Adi Diferansiyel denklemlerin köklerinin bulunması -LU ayrıştırması
-Öz değer, Öz vektör
-Lineer olamayan denklemlerin köklerinin bulunması
Tek değişkenli fonksiyonlar
Yarılama Metodu :
[a,b] aralığında köke bakılır p1 bu aralığın orta noktasıdır.
f(p1)=0 ise p=p1 dir.
f(p1) ve f(a1) aynı işaretliyse kök [p1,b1], a2=p1, b2=b1 dir.
f(p1) ve f(a1) farklı işaretliyse kök [a1,p1], a2=a1, b2=p1 dir.
Algoritması:
Adım 1: i=1 FA=f(a)
Adım 2: while i<=N0 do steps 3-6 Adım 3: p=(a+b)/2
FP=f(p)
Adım 4: if FP=0 or (b-a)/2<TOL then OUTPUT(p)
STOP Adım 5: i=i+1
Adım 6: if FA.FP>0 then a=p; FA=FP else b=p
Adım 7: OUTPUT(‘Method failed after N0 iterations, N0=‘,N0) STOP
10 4
)
( x x
3 x
2 f
Örnek:
fonksiyonunun [1,2] aralığında kökünün olup olmadığına bakılacak:
f(1)= - 5 , f(2)= 14
Sabit Nokta Metodu :
10 4
)
( x x
3 x
2 f
Örnek:
fonksiyonunun [1,2] aralığında kökünün olup olmadığına bakılacak:
f(1)= - 5 , f(2)= 14
10 4
)
( x x
3 x
2
f
Newton-Raphson Metodu :
Secant Metodu :
Regula-Falsi Metodu :
Interpolasyon
Interpolasyon
Lagrange Interpolasyonu :
Newton’un Bölünmüş Farklar Metodu :
Newton’un Bölünmüş Farklar Metodu :
Örnek : Aşağıdaki nokta değerleri verilmiş olsun.
Newton’un İleri Bölünmüş Farklar Metodu :
Newton’un geri Bölünmüş Farklar Metodu :
Newton’un orta Bölünmüş Farklar Metodu :
Geri Fark Formülü : İleri Fark Formülü :
Orta Fark Formülü :
İleri=
Geri=
Orta=
Hermit Metodu :
Sayısal Türev
Sayısal İntegral
Adi Diferansiyel Denklemler
Euler Metodu :
y’=f(t,y) , [t0,tM] , y(t0)=y0
tk=a+kh , k=0,1,2,…,M , h=(b-a)/M
y(t) fonksiyonu t=t0 noktasında Taylor serisine açıldığında :
y’(t0)=f(t0,y(t0)), h=t1-t0 , değerlerini yerine yazarak
elde edilir.
2
) )(
) ( )(
( )
( )
(
2 0 1
0 0
0
t t c t y
t t y t
y t
y
) 2 ( ))
( , ( )
( )
(
2 1 0
0 0
1
c h y t
y t hf t
y t
y
h yeterince küçük seçilirse ikinci türeve sahip terim ihmal edilebilir.
Bu durumda :
Heun Metodu :
Taylor Metodu :
Runge-Kutta Metodu :
Örnek :