Ana Sayfa » Ankara Üniversitesi ANKADEM » Mühendislik Fakültesi » Computer Simulation and Modelling

Computer Simulation and Modelling

Ünite 7: Numericanalysis

Ankara Üniversitesi Mühendislik Fakültesi

Bu ders materyali, Ankara Üniversitesi Bilgisayar Mühendisliği Bölümü'nde 2014-2015 Bahar Dönemi'nde Semra Gündüç tarafından verilen 'Computer Simulations' dersi kapsamında hazırlanmıştır.

İleri Yönlü Fark Formülü Hatası

İki noktalı ileri yönlü fark formülü kullanılarak hesaplanan sayısal birinci türevin hata mertebesi O(Delta x) düzeyindedir. Bu durum, adım boyutu yarıya indirildiğinde hatanın da yaklaşık yarıya ineceği anlamına gelir.

Üç Noktalı Merkezi Fark Formülü

x_i+1 ve x_i-1 noktalarındaki Taylor açılımlarının birbirinden çıkarılmasıyla elde edilen birinci türev formülü, O(Delta x^2) mertebesinde yüksek doğruluğa sahiptir.

Beş Noktalı Birinci Türev Formülü

x_i-2, x_i-1, x_i+1 ve x_i+2 noktalarındaki değerleri kullanan formül, türevi 1/(8*Delta x) * [f(x_i-2) - 8f(x_i-1) + 8f(x_i+1) + f(x_i+2)] + O(Delta x^4) şeklinde hesaplayarak çok yüksek doğruluk sunar.

Anahtar Kavramlar

Taylor Serisi AçılımıSürekli bir fonksiyonun, belirli bir nokta etrafındaki sonsuz türevleri kullanılarak bir polinom şeklinde ifade edilmesidir. Sayısal analizde fonksiyonları basitleştirmek ve yaklaşıklık hesabı yapmak için kullanılır.

Ayrık Uzay Ağı (Discrete Space-Mesh)Sürekli bir uzay veya aralığın, bilgisayarda işlenebilmesi için belirli sayıda sonlu ve ayrık noktalara bölünmesiyle oluşan yapıdır. İki nokta arasındaki mesafe adım boyutu olarak adlandırılır.

İleri Yönlü Fark (Forward Difference)Bir fonksiyonun türevini, mevcut nokta ile bir sonraki nokta arasındaki farkı adım boyutuna bölerek tahmin eden birinci derece sayısal türev yöntemidir.

Geri Yönlü Fark (Backward Difference)Bir fonksiyonun türevini, mevcut nokta ile bir önceki nokta arasındaki farkı adım boyutuna bölerek hesaplayan sayısal türev yaklaşımıdır.

Merkezi Fark (Central Difference)Türev hesabında mevcut noktanın hem bir adım ilerisindeki hem de bir adım gerisindeki noktaları kullanarak daha yüksek doğruluk (ikinci derece) sağlayan yöntemdir.

Kesme Hatası (Truncation Error)Sonsuz bir serinin (örneğin Taylor serisinin) belirli bir terimden sonra kesilerek sonlu sayıda terimle hesaplanması sonucu oluşan matematiksel hatadır.

Makine Hassasiyeti (Machine Accuracy)Bir bilgisayarın sayısal veri tipleriyle temsil edebileceği en küçük anlamlı fark veya limit değerdir. Bu limitin altındaki adımlar yuvarlama hatasına yol açar.

İnterpolasyon (Interpolation)Değerleri sadece belirli ayrık noktalarda bilinen bir fonksiyonun, bu noktaların arasında kalan ara değerlerini tahmin etmek için kullanılan matematiksel yöntemdir.

Kuadratik Polinomİkinci dereceden (ax^2 + bx + c formunda) olan polinomdur. Üç veri noktasından geçen eğrileri modellemek ve daha hassas interpolasyon yapmak için kullanılır.

Fonksiyonun KöküBir f(x) fonksiyonunu sıfıra eşitleyen, yani f(x) = 0 denklemini sağlayan x değeridir. Sayısal analizde bu değer tolerans sınırları dahilinde aranır.

Tolerans (Epsilon)Sayısal hesaplamalarda, bir sonucun kabul edilebilir derecede doğru olduğunu belirten, sıfıra çok yakın seçilen eşik değerdir (örneğin 10^-6).

Newton-Raphson YöntemiBir fonksiyonun kökünü bulmak için teğet doğrularını kullanan, başlangıç noktasına bağlı olarak çok hızlı (kuadratik) yakınsayan yinelemeli bir algoritmadır.

Secant YöntemiNewton-Raphson yöntemindeki türev hesaplama zorunluluğunu ortadan kaldırmak için türev yerine iki noktalı fark yaklaşımı kullanan kök bulma yöntemidir.

Kuadratik YakınsamaBir yinelemeli yöntemin her adımında hata miktarının karesiyle orantılı olarak azalması, yani doğru basamak sayısının her iterasyonda yaklaşık iki katına çıkması durumudur.

Kök Katlılığı (Multiplicity)Bir kökün denklemde birden fazla kez tekrarlanması durumudur (örneğin (x-a)^P = 0 için P katlılık derecesidir). Katlılık arttıkça standart yöntemlerin yakınsama hızı düşer.

Bisection (İkiye Bölme) AlgoritmasıFonksiyonun işaret değiştirdiği bir aralığı sürekli ikiye bölerek kökü arayan, yavaş ama her zaman yakınsayan garantili bir kök bulma yöntemidir.

Simpson İntegral KuralıEntegrasyon aralığını çift sayıda eşit parçaya bölüp, her bölgede fonksiyonu parabolik eğrilerle yaklaştırarak alan hesaplayan yüksek doğruluklu sayısal integral yöntemidir.

Sıradan Diferansiyel Denklem (ODE)Tek bir bağımsız değişkene ve bu değişkene bağlı fonksiyonun türevlerine sahip olan matematiksel denklemdir.

Euler YöntemiBirinci mertebeden diferansiyel denklemleri çözmek için kullanılan, sonraki adımı sadece mevcut noktadaki eğimle tahmin eden en basit ve ilk derece sayısal entegrasyon yöntemidir.

Runge-Kutta YöntemiDiferansiyel denklemlerin çözümünde, aralık içindeki farklı noktalarda eğim hesaplamaları yaparak yüksek doğruluk sağlayan tek adımlı yinelemeli yöntemler ailesidir.

Kısmi TürevÇok değişkenli bir fonksiyonun, diğer tüm değişkenleri sabit tutulurken sadece tek bir değişkene göre alınan türevidir.

Runge FenomeniEşit aralıklı noktalarda çok yüksek dereceli polinomlarla interpolasyon yapıldığında, aralığın uç noktalarında meydana gelen aşırı büyük salınım ve hata durumudur.

Adım Boyutu (Delta x)Sayısal analizde ayrıklaştırılmış noktalar arasındaki doğrusal mesafedir. Adım boyutu küçüldükçe doğruluk artar ancak hesaplama süresi uzar.

Aşkın (Transcendental) DenklemlerTrigonometrik, logaritmik veya üstel fonksiyonlar içeren, cebirsel yöntemlerle doğrudan çözülemeyen ve sayısal yöntemlerle kökü aranan denklemlerdir.

Garantili YakınsamaBaşlangıç noktasından bağımsız olarak, belirli koşullar altında algoritmanın kesinlikle çözüme (köke veya değere) ulaşacağını ifade eden özelliktir.

Diğer Önemli Bilgiler

İkinci Derece Türevin Beş Noktalı Formülü

İkinci türevin doğruluğunu artırmak için kullanılan beş noktalı formül, O(Delta x^4) hata payı ile çalışır ve katsayıları sırasıyla -1, -16, -30, 16, -1 şeklindedir.

Sinüs Fonksiyonu İnterpolasyon Örneği

Metinde sunulan örnekte, sin(0.80)=0.71736 ve sin(0.82)=0.73115 değerleri kullanılarak yapılan interpolasyonla sin(0.81) değeri 0.72426 olarak hesaplanmış ve gerçek değere çok yakın bulunmuştur.

Newton-Raphson Maksimum Yineleme Sınırı

Sunulan C kodlarında, sonsuz döngüleri ve ıraksamayı engellemek amacıyla maksimum yineleme sayısı (MAXCOUNT) makro tanımıyla 100 olarak sınırlandırılmıştır.

Katlı Köklerde Yakınsama Hızı

Katlılığı P > 1 olan köklerde standart Newton-Raphson yönteminin hata yakınsama hızı e_n+1 = [(P-1)/P] * e_n formülüyle doğrusal (lineer) hıza düşer.

x^4 - x = 10 Denkleminin Çözümü

Newton-Raphson yöntemiyle x_0=2 başlangıç değeri seçildiğinde, kök değerleri sırasıyla x_1=1.871, x_2=1.85578 ve x_3=1.855585 olarak hızlıca yakınsar.

Karekök Bulma Algoritması

Bir 'a' sayısının karekökünü bulmak için f(x) = x^2 - a = 0 fonksiyonu tanımlanır. Newton-Raphson formülü uygulandığında x_n+1 = x_n - (x_n^2 - a)/(2*x_n) yinelemesi elde edilir.

Simpson Yöntemi Hata Mertebesi

Simpson entegrasyon kuralında, simetrik aralıklar nedeniyle tek dereceli terimler birbirini yok eder ve hata mertebesi beklenmedik şekilde yüksek olan O(Delta x^5) düzeyine ulaşır.

Euler Yöntemi Hata Sınırı

Diferansiyel denklemlerin çözümünde kullanılan en basit yöntem olan Euler algoritmasının yerel kesme hatası O(Delta x^2), birikimli küresel hatası ise O(Delta x) düzeyindedir.

Dördüncü Derece Runge-Kutta (RK4) Doğruluğu

RK4 yöntemi, her adımda dört farklı eğim tahmini (k1, k2, k3, k4) hesaplayarak diferansiyel denklemleri O(h^5) mertebesinde yerel hata ile çözer.

Çift Sayıda Aralık Zorunluluğu

Simpson 1/3 entegrasyon kuralının uygulanabilmesi için toplam alt aralık sayısının (N) mutlaka çift bir sayı olması gerekmektedir.

Karmaşık Köklerin Başlangıç Koşulu

Gerçek katsayılı cebirsel denklemlerde, başlangıç değeri (x_0) gerçek sayı seçilirse Newton-Raphson yöntemiyle karmaşık köklere ulaşmak matematiksel olarak imkansızdır.

Sınavda Dikkat Et

Türev formüllerinin hata mertebelerine dikkat edin; merkezi fark O(Delta x^2) iken ileri/geri fark O(Delta x) mertebesindedir. Sınavda formül karşılaştırması sorulabilir.
Simpson yönteminin uygulanabilmesi için aralık sayısının (N) çift olması gerektiğini unutmayın; bu kural doğrudan doğru-yanlış veya boşluk doldurma sorusu olarak gelebilir.
Newton-Raphson yönteminin yakınsama hızının kuadratik olduğunu, ancak katlı köklerde (P > 1) bu hızın doğrusala düştüğünü ve P katsayısı ile modifiye edilmesi gerektiğini aklınızda bulundurun.
Gerçek katsayılı denklemlerde karmaşık köklere ulaşmak için başlangıç değerinin de karmaşık seçilmesi gerektiği bilgisini sınavda karıştırmayın.
Diferansiyel denklemlerin çözümünde RK4 yönteminin hata mertebesinin O(h^5) olduğunu ve Euler yöntemine göre çok daha kararlı sonuçlar verdiğini unutmayın.
Sayısal analizde adım boyutu Delta x'in aşırı küçültülmesinin makine hassasiyeti sınırları nedeniyle her zaman doğruluğu artırmayacağını, aksine kararsızlığa yol açabileceğini bilin.