Buders Sayısal Analiz

1 SAYISAL ANALİZ 1

2 SAYISAL ANALİZ 3. Hafta SAYISAL ANALİZDE HATA KAVRAMI VE HATA TÜRLERİ 2

3 İÇİNDEKİLER 1. de Problem Çözümünde İzlenilecek Adımlar 2. Matematiksel Modelleme 3. de Hata Kavramı 4. de Hataların Sebepleri 5. Sayısal Hata ve Hata Türleri Mutlak Hata Bağıl Hata Yaklaşım Hataları Yüzde Hata Hata Sınırlaması 6. Anlamlı Basamak 7. Sayıların Bilgisayarda Gösterilmesi 8. Kesme ve Yuvarlama Hataları 3

4 Giriş Sayısal analiz, matematik problemlerinin bilgisayar yardımı ile çözümlenme tekniğidir. Genellikle analitik olarak çözümleri çok zor veya imkansız olan matematiksel problemleri belirli bir hata aralığında çözmek için kullanılır. Giriş Bilgisi Hesap Yöntemi Çıkış Bilgisi DATA ALGORİTMA SONUÇ Sayısal analiz işleyiş şeması Sayısal analiz yöntemleri ile bulunan/hesaplanan sonuçlar kesin değerler değil yaklaşık değerlerdir. 4

5 de Problem Çözümünde İzlenilecek Adımlar (Detaylı) Problem Mühendislik problemi tanımlanır. Matematiksel Modelleme Problemin çözümü için matematiksel terimlerle ifade edilen model oluşturulur. Uygun matematiksel model tespit edilir. Tasarım/Algoritma Geliştirme Programlama için çözüme en uygun algoritma/akış diyagramı tasarlanır. Kodlama Algoritma/akış diyagramı istenilen bir programlama dili ile kodlanır. Test Program çalıştırılır, sonuçlar, grafikler elde edilir. Yazılım/Uygulama En uygun (optimum) çözümün elde edildiği program (İterasyonlardan sonra ) Not: Hatalara göre süreç, tekrar gözden geçirilir, işlemler tekrar edilir. 5

6 Matematiksel Modelleme Bağımlı değişken = f ( bağımsız değişkenler, parametreler, zorlayıcı fonksiyonlar) Bağımlı değişken, sistemin davranışını ya da konumunu gösteren karakteristik Bağımsız değişkenler, sistemin davranışının incelendiği boyutlar (zaman, konum vb.) Parametreler, sistemin özelliklerini, yapısını ya da bileşenlerini yansıtan parametreler Zorlayıcı fonksiyonlar, sisteme etkiyen dış etkiler Örnek: Bir paraşütçünün düşme hızının hesabı (Newton 2. yasası) v Analitik Çözüm Modeli gm = c ( ( c/ mt ) 1 e ) v: hız (diğer zorlayıcı kuvvetlere, parametrelere ve bağımsız değişkenlere bağlı olarak değişen bir bağımlı değişkendir. t: zaman (bağımsız değişken) g: yer çekimi sabiti (zorlayıcı kuvvet) c: havanın direnç katsayısı (sistemin fiziksel özelliği) m: kütle (sistemin fiziksel özelliği) Sayısal Çözüm Modeli Yeni değer = Eski değer + eğim * adım uzunluğu (Euler Yöntemi) 6

7 de Hata Kavramı, matematiksel hesaplamaları ve problemleri tekrarlı işlemler ile bilgisayarlar (programlar) aracılığı ile çözmektir. Özetle,, bir problemi hesap makinesi, bilgisayar gibi araçlarla yaklaşık olarak çözmek için kullanılan yöntemlerdir. Bilgisayarlı sayısal analizde çeşitli sebeplerle hatalı sonuçlar elde edilebilir. Kullanıcıdan kaynaklanan hatalar Modelleme hatası Kodlama hatası Bilgisayar ya da kullanılan yazılımdan kaynaklanan hatalar Belirli bir uzunlukta sayı depolayabilme, Yuvarlatmadan kaynaklanan hata, Kesmeden kaynaklanan hata 7

8 de Hataların Sebepleri Fiziksel veya sosyal olayların matematiksel olarak çözülmelerinde yapılan hatalar genellikle üç ana başlıkta toplanır. Modelleme Hatası, bir olayın formüle edilmesi esnasında varsayımlardan kaynaklanan hatalardır. Örnek: Serbest düşme problemlerinin modellenmesinde, hava ile cisim arasındaki sürtünme kuvvetinin ihmal edilmesinden dolayı meydana gelen hatalar bu tür hatalar grubuna girer. Ölçme Hatası, deney ve gözlemede ölçmelerden dolayı meydana gelen hatalardır. Örnekte, eğer serbest düşme yapan cismin, düştüğü mesafe veya havada düşerken gecen süre eğer yanlış ölçülürse bu tür hatalar ölçme hatası olarak tanımlanabilir. Sayısal hatalar veya diğer bir deyimle modelin çözümlemesinde yapılan hatalardır. Dersimizde sayısal hatalar ile ilgileneceğiz. 8

9 Sayısal Hata Sayısal hatalar, matematiksel işlemler ve değerlerin yaklaşık kullanımlarından ortaya çıkan farklardır. Doğruluk: Ölçme veya hesaplama sonuçlarının gerçek değere yakınlığıdır. Çoğu durumda kesinliğin değeri bilinmediğinden en iyi doğru cevap genellikle en iyi tahmin edilen değere göre ölçülür. Kesinlik (hassasiyet): Bir büyüklüğün defalarca ölçülmesi sonucu ölçümlerin birbiriyle ne kadar yakınlıkta olduğudur. Bilgisayarlı hesaplama doğruluğu, bilgisayarın kelime uzunluğu (bir kerede işlenen bit sayısı), ile doğrudan ilişkilidir. Bilgisayarlar sonlu sayıda rakamı saklayabilirler. Bu sebeple, hesaplamalar, tam değil, yaklaştırmalar ile yapılabilir. Hesaplamalarda ne kadar hata vardır. Ne kadarlık bir hata kabul edilebilir 9

10 Sayısal Hata Türleri Mutlak Hata Bir büyüklüğün, analitik olarak (doğru/gerçek) bilinen değeri ile sayısal hesaplamalarda elde edilen yaklaşık (hesaplanan) değeri arasındaki mutlak farktır Hesaplamalarda, genellikle hatanınişaretiyle değil, büyüklüğüyle ilgilenilir. Bu sebeple hata kavramında mutlak değer kullanılır. 10

11 Sayısal Hata Türleri BBaağıııı = (ggggggçeeee ddddğeeee) (yyyyyyyyyyşıııı ddddğeeee) ggggggçeeee ddddğeeee MMMMMMMMMMMM HHHHHHHH = GGGGGGçeeee DDDDğeeee 11

12 Sayısal Hata Türleri Yaklaşım Hataları Uygulamalarda genellikle, gerçek sonuç baştan bilinemeyeceğinden, bu durumda hata, gerçek değerin bilinen en iyi tahminine göre normalize edilir. Sayısal yöntemler bir sonucu bulmak için genellikle iteratif yöntemler kullanır. İterasyon, bir işlev ya da yöntemin tekrar tekrar uygulanmasıdır. Yöntem doğru ise her adımda gerçek değere biraz daha yaklaşılır. Dizi ya da tekrarlı işlemlerde kullanılır. Yaklaşık Mutlak Hata aa = (eeee iiiiii tttthmmmmmm) (yyyyyyyyyyşıkk ddddğeeee) Yaklaşık Bağıl Hata aa = (eeee iiiiii tttthmmmmmm) (yyyyyyyyyyşıkk ddddğeeee) eeee iiiiii tttthmmmmmm 12

13 Sayısal Hata Türleri Yaklaşım Hataları Bir niceliğin tam değeri verilmediği zaman gerçek bir hatanın hesaplanması mümkün değildir. Yaklaşık hata, hatanın sınırlarını belirlemede kullanılır. Bu durumda Scarborough Kriteri kullanılabilinir. Eğer yaklaşık bağıl hata e a < m ise o zaman sonuç m nin en küçük basamağı için doğrudur. ya da Aşağıdaki Scarborough kriteri gerçekleştiğinde, sonucun en az n anlamlı basamak için kesinlikle doğrudur denilebilir. ss = %( nn ) 13

14 Sayısal Hata Türleri Yüzde Hata Bağıl, mutlak, yaklaşım hatası gibi tüm hata ifadelerinin % ile çarpılması ile % değerleri elde edilir. Buna göre hatalar, bağıl yüzde hata, yüzde yaklaşım hatası şeklinde adlandırılır. 14

15 Sayısal Hata Türleri Örnek 1: Gerçek değeri cm olan bir elektrik kablosunun uzunluğu ölçüldüğünde cm bulunuyor. Mutlak hatayı, bağıl hatayı ve yüzde bağıl hatayı hesaplayınız? Mutlak Hata: εε mmmmmmmmmmmm = (ggggggçeeee ddddğeeee) (yyyyyyyyyyşıkk ddddğeeee) = = 1 cccc Bağıl Hata: BBaağıııı = (ggggggçeeee ddddğeeee) (yyyyyyyyyyşıııı ddddğeeee) ggggggçeeee ddddğeeee = MMMMMMMMMMMM HHHHHHHH GGGGGGçeeee DDDDğeeee Yüzde Bağıl Hata: BBBBğıll= = 0, % BBBBğıll== BBBBğıll = 0, = % 0,5 15

16 Anlamlı Basamak Aracın hızı tam olarak kaç km/s Aracın kilometresi kaçtır Bir sayının anlamlı basamakları, güvenle kullanılabilecek basamaklarıdır. Chapra ve Canale nin Kitabından Alınmıştır. 16

17 Anlamlı Basamak Mutlak ve Bağıl hata hesaplamasında gerçek değere ihtiyaç vardır. Ancak gerçek değer bilinmediğinde, özellikle dizi (seri) şeklindeki hesaplamalarda yaklaşık hata hesaplaması kullanılır. Bu işlemde son iki yineleme arasındaki rakamların kaç tanesinin tekrar ettiğine bakılarak fonksiyonun gerçek değerine ne kadar yaklaştığı (yakınsadığı) kontrol edilebilir. Her yeni tekrardaki değer, bir önceki değer ile yakınsa sonuç ya da ölçüm hassastır denilebilir. Ama doğruluk tartışılır. Örnek: Tekrar Sayısı Sonuç Anlamlı Basamak Sayısı 1 f(x)=abcdef 0 2 f(x)=abghik 2 17

18 Anlamlı Basamak Örnek 2: Sin(x) in seri açılımı aşağıda görülmektedir. sin(π/6) değerinin gerçek değerini bulup, yukarıdaki seri açılıma göre sin(π/6) için terim sayısına göre, mutlak ve bağıl hataları ile anlamlı basamak sayılarını gösteriniz? 18

19 Sayısal Hata Türleri Örnek 3: Scarborough kriteri ile hataları anlamlı sayılar ile ilişkilendirmek? ~ π = π = π = n ~ π = π π = s ( 2 ) % *10 n Serhat Yılmaz ın Sunusundan Alınmıştır. 19

20 Sayısal Hata Türleri Matematikte fonksiyonlar/işlevler çoğu zaman sonsuz seriler ile gösterilir. Matematiksel bir ifadenin seri açılımında terim sayısı arttıkça, hesaplanan değer gerçeğe daha yaklaşır (Mac Laurin Serisi Açılımı). Örnek 4: e x ifadesinin seri açılımı aşağıda görülmektedir. e 2 3 n x x x x = 1 + x ! 3! n! x= değeri için e in gerçek değeri Serinin birinci teriminden başlayarak, tek tek terim ekleyerek, e sayısını tahmin edin. Her adımda (terim eklemede) gerçek ve yüzde bağıl hataları bulun. Yaklaşık mutlak hata değeri, 3 anlamlı basamak veren belirli bir ε s kriterinden daha küçük oluncaya kadar terim eklemeye devam edin. hata Serhat Yılmaz ın Sunusundan Alınmıştır. 20

21 Sayısal Hata Türleri Çözüm 4: İlk olarak sonucun en az 3 üç anlamlı basamak için doğru olmasını garanti eden hata kriteri belirlenebilir. ss = %( nn ) ( 2 3) = % *10 = % s Bu seviyenin altına inene kadar terim eklenmeye devam edilir. Serhat Yılmaz ın Sunusundan Alınmıştır. 21

22 Sayısal Hata Türleri Çözüm 4: İlk terimden başlayarak terim ekleme, İlk terim için; e x = e = 1, Bağıl yüzde hata, % r = 1, *% 1, = % İlk iki terim için; e x = 1 + x e = = Bağıl yüzde hata, % r 1, = * % 1, = % Yaklaşık bağıl yüzde hatanın hesaplanması: 1. adımda elimizdeki tek tahmin değeri bulunduğundan yaklaşık hata hesaplanamaz. 2. adıma geçtiğimizde, ilk adımda bir önceki tahmin değerine sahip olduğumuzdan yaklaşık hata hesaplanabilir = a *% = % Bu değer, % hata sınırından küçük olmadığından yeni terim eklenmeye devam edilir a > s Serhat Yılmaz ın Sunusundan Alınmıştır. e = =

23 Sayısal Hata Türleri Örnek 4 ün Çözüm Algoritması ve MATLAB Programı Başla İlk Değerleri Ata e in en son yaklaşık değerini göster H % > a? s E Döngü İçindeki Hesaplamaları (Terim ekle, Toplamı bul, hatayı hesapla ) Serhat Yılmaz ın Sunusundan Alınmıştır. 23

24 Sayıların Bilgisayarda Gösterilmesi Matematiksel hesaplamalarda farklı sayı sistemleri kullanılabilir. Hexadecimal(16), Decimal(10), Octal(8) ve binary(2) Decimal olarak sayısının gösterimi; = 1* *10 + 3*1 = 1* * *10 0 Bilgisayarlar ise ikilik (binary) sayı sistemine göre çalışır. sayısının bilgisayarda binary olarak gösterimi; = 1* * * * * * *2 0 = ( ) 2 Elektronik hesaplayıcılarda sayılar iki tabanında ancak belirli uzunlukta ifade edilebilir. Örneğin, reel sayılar için normal hassasiyette 32 bitlik bir yer ayrılan hesaplayıcıda 7 ondalık basamağa, çift hassasiyette ise 64 bitlik yer ayrılır ve buda yaklaşık 15 ondalık basamağa karşılık gelir. Bu nedenle değerler için hesaplayıcılardaki ayrılan yerler veri tipine göre değişmektedir. Buda farklı bir türde hataya neden olabilmektedir. 24

25 Sayıların Bilgisayarda Gösterilmesi Sayıların temsili; x = m E < m < 1 x 0 için m mantis, E üs tür. 0 sayısı toplam olarak normal notasyona sahiptir Bir sayının normal ya da kayan noktalı formda temsili x = 0.d 1 d 2 d d k 10 n d 1 0 ve d k 0 la bir sayının önemli basamakları olduğunu ve sayının gerçek değerine göre güvenle kullanılabilen basamaklarıdır. Ondalık noktaları kaydırmak için kullanılan 0 noktaları önemli basamaklar olarak sayılmaz. Kalan sıfırlar önemli olabilirde olmayabilirde. Örnek; x=0, Önemli figüre sahiptir. x=1, 4,3 veya 2 önemli figüre sahip olabilir. 25

26 Sayıların Bilgisayarda Gösterilmesi Sayıların Kayan Noktalı Gösterilimi (floating Point) x= (işaret) m*b (işaret)e Burada m mantis, b taban (binary sistem için b=2) ve E üs dür. 8 bitlik işlemciler için kayan noktalı gösterim, Mantis işareti Mantis basamakları Üs işareti Üs basamakları d 1 d 2 d 3 d 4 d 5 d 6 d 7 d 8 Örnek: Mantis işareti Mantis basamakları Üs işareti Üs basamakları m = +( ) m = +( ) = 1 E = -[(1 2 1 ) + (1 2 0 )] = -(2+1) = - 3 sayı = (desimal sistemde 0, sayısına eşittir. monash.pw Asaf Varol un Sunusundan Alınmıştır. 26

27 Kesme ve Yuvarlama Hataları Bilgisayarlar, işleyebildikleri kelime uzunluğunun sınırlı sebebiyle (sınırlı sayıda bilgiyi saklayabilmek) sayı değerlerinin sonlu olması gerekir. Hesaplama sırasında bir sayının izin verilemeyen (saklanamayan) kısmının atılması, hesaplanan değerin hatalı olarak bulunmasına sebep olur. Bu durumda karşılaşılan hatalara kesme ve yuvarlama hataları denir. Kesme hatası, olarak algılanır. Bu durumda hata= olarak algılanır. Bu durumda hata= Yuvarlama hatası, olarak algılanır. Bu durumda hata= olarak algılanır. Bu durumda hata=

28 Kesme ve Yuvarlama Hataları Yuvarlama Hatası: Örnek olarak 1/3 kesrini bilgisayar gibi belli adet hane kullanarak yazar. Sayıların tanımlanması için kaç hane kullanılacağı rakamların nasıl tanımlandığı ve bilgisayarın mimarisi ile ilgilidir. Bu tür hatalara yuvarlama hatası (round-off error) denir. 28

29 Kesme ve Yuvarlama Hataları Kesme hatası, verilerin sayısal işlemlere girmesinden kaynaklanan hatadır. Sonsuz terimli bir seriyi uygun şekilde keserek sayısal sonuçlar elde edilir. Belirli terimden sonra gelen terimlerin ihmal edilmesi kesme hatası olarak bilinir. Burada yapılan hata atılan terimlerin toplamı kadar olur. Hesaplamada ihmal edilen terimlerin toplamı yapılan kesme hatasına eşit olur. Örnek: Kesme hatasının oluşması: 29

30 Kesme ve Yuvarlama Hataları Örnek: Körfez savaşı esnasında ( ), kesme hatası sebebiyle patriot füzesi, skud füzesini algılayamadı. Bilgisayar saatinden alınan değer 1/10 ile çarpılarak geçen süre hesaplanmaktaydı. Bilgisayar bit sabit kayıtçı kullanmaktaydı. Serhat Yılmaz ın Sunusundan Alınmıştır. 30

31 Kesme ve Yuvarlama Hataları 1/10 sayısı 2 li tabanda 1/2 +1/2 +1/2 +1/2 +1/2 +1/ şeklinde sonsuz terimlidir şeklinde Sayı: bit: Kesme hatası: saatlik batarya süresince bu hata: = sn Saniyede 1, metre yol alan skud füzesi saniyede metreden fazla gider. Bu kadarlık hata ise Skud un menzil dışında görülmesine yol açar. Serhat Yılmaz ın Sunusundan Alınmıştır. 31

32 Örnek Bir matematikçi dairenin alanını analitik olarak A D = πr 2 bulmuştur. Sadece ikiz kenar üçgenlerin alanını hesaplayabilen bir bilgisayarla bu dairenin alanını % s hata ile bulabilecek bir algoritma oluşturun ve programını MATLAB programlama dilinde yazınız. Program, her adımda daireyi üçgenlerden oluşan n eşit parçaya bölsün, istenen doğruluk sağlanmadıkça n i 1 arttırsın (r=1 cm olsun, n, 5 ten başlasın). Daireyi yaklaşık olarak oluşturacak n tane üçgenin alanı ; n*a ü çgen ( 2 r = n* *sin α ) 2 A D α n α = 2π n 32

33 Başla Çözüm İlk Değerleri Ata: n, r, hata, tol, A0 Son Alan değerini göster H hata tol mı? E 2 r = A 2 A1 n* *sin ( α ) D % s A = 1 A 0 * % A 1 A0 A1 n=n+1 33

34 Çalışma Sorusu f(x) = e -1,5 fonksiyonu için; İlk terimden başlayarak, 5 terim için adım adım f(x) = e -1,5 fonksiyonunu hesaplayınız. Gerçek değer olarak kabul edip, her adım da yüzde mutlak hata, yüzde bağıl hatayı hesaplayınız. Gerçek değeri bilmediğinizi kabul edip her adım da yüzde yaklaşım hata ve yüzde yaklaşım bağıl hatayı hesaplayınız. Her iki çözümde de varsa anlamlı basamak sayısını gösteriniz. e x = 1 x + x x 3! + e x = 1 x e = 1+ x x x 3!

35 KAYNAKLAR Serhat YILMAZ, Bilgisayar İle Sayısal Çözümleme, Kocaeli Üniv. Yayınları, No, Kocaeli, Steven C. Chapra, Raymond P. Canale (Çev. H. Heperkan ve U. Kesgin), Yazılım ve Programlama Uygulamalarıyla Mühendisler İçin Sayısal Yöntemler, Literatür Yayıncılık. İlyas ÇANKAYA, Devrim AKGÜN, Sezgin KAÇAR Mühendislik Uygulamaları İçin MATLAB, Seçkin Yayıncılık Yüksel YURTAY, Ders Notları, Sakarya Üniversitesi monash.pw Asaf Varol, Ders Notları, Fırat Üniversitesi 35