Gauss-Legendre Kuadratürü

Konunun Amacı

İntegrasyon noktalarını eşit aralık yerine en yüksek polinom doğruluğunu verecek apsis ve ağırlıklarda seçerek belirli integrali hesaplama. Bu ders Sayısal Analiz ve Numerik Analiz adlarıyla da anılır. Amaç yalnız bir formülü uygulamak değil; yöntemin hangi varsayımla çalıştığını, ne zaman yakınsadığını ve yaklaşık sonucun nasıl doğrulanacağını öğrenmektir.

Mühendislikte Neden Gereklidir?

Sonlu elemanlarda eleman rijitliği ve yük integralleri az sayıda fonksiyon değerlendirmesiyle Gauss noktalarında hesaplanır. Gerçek bir mühendislik hesabında veri hassasiyeti, model kabulleri ve bilgisayarın sonlu basamakla çalışması sonucu etkiler. Bu nedenle sonuç; birimi, hata ölçüsü ve kullanılan toleransla birlikte raporlanır.

Bir algoritmanın hızlı olması tek başına yeterli değildir. Kararlı olmayan bir yöntem küçük veri veya yuvarlama değişimlerini büyütebilir. Aynı problem farklı adım büyüklüğü, başlangıç tahmini veya ağ yoğunluğuyla yeniden çözüldüğünde sonuçların birbirine yaklaşması beklenir.

Temel Kavramlar

1. Gauss apsisi

Tanım: Fonksiyonun değerlendirildiği özel tᵢ integrasyon noktasıdır.

Bu kavram hesap tablosunda yalnız bir etiket olarak bırakılmaz. Kullanılan değer, birim, işaret kabulü ve geçerli olduğu nokta açıkça belirtilir. Sonuç yorumlanırken yöntemin varsayımıyla birlikte okunur.

2. Gauss ağırlığı

Tanım: Her apsisteki fonksiyon değerini çarpan Wᵢ katsayısıdır.

3. Standart aralık

Tanım: Gauss-Legendre tablolarının tanımlandığı [-1,1] aralığıdır.

4. Koordinat dönüşümü

Tanım: [a,b] aralığını standart t aralığına taşıyan doğrusal eşlemedir.

Temel Bağıntı

∫[-1,1]F(t)dt≈Σ WᵢF(tᵢ); x=(b−a)t/2+(a+b)/2

Standart aralık dışındaki integral önce [-1,1]'e dönüştürülür ve Jacobian (b−a)/2 çarpanı eklenir. N noktalı Gauss kuralı 2N−1 dereceye kadar polinomlarda kesindir. Bağıntıya değer yerleştirmeden önce indislerin, vektör sırasının ve birimlerin tutarlı olduğu kontrol edilir. Ara değerler mümkün olduğunca tam hassasiyetle taşınır; yalnız nihai sonuç uygun basamakta yuvarlanır.

Çözüm Algoritması

[a,b] aralığını t∈[-1,1] aralığına dönüştür.
Dönüşüm Jacobian'ını belirle.
Seçilen N için Gauss apsis ve ağırlıklarını al.
Dönüştürülmüş fonksiyonu apsislerde hesapla.
Ağırlıklı toplamı Jacobian ile çarp.

Her adımın çıktısı bir sonraki adımın girdisidir. Bu nedenle ara tabloda iterasyon numarası, tahmin, artık, hata ölçüsü ve varsa ölçekleme katsayısı ayrı sütunlarda tutulur. Algoritma durduğunda hangi ölçütün sağlandığı açıkça yazılır.

Çözümlü Örnek

İki noktalı kuralda t=±1/√3 ve ağırlıklar 1'dir. F(t)=t² için toplam 1/3+1/3=2/3 olur; bu değer integralin kesin sonucudur.

Bu örnekte hesaplanan sayının yanında işlemin yapısı da önemlidir. Sonucun makul mertebede olduğu yaklaşık bir hesapla kontrol edilir. Bir iterasyon yöntemi kullanılıyorsa son tahmin denkleme geri yazılarak artık bulunur; bir integrasyon veya türev yaklaşımında ise daha küçük adımla ikinci hesap yapılır.

Hata, Yakınsama ve Kararlılık Kontrolü

Artık kontrolü: Yaklaşık çözümü özgün denklemde yerine koy ve sağlanmayan kısmı ölç.
Adım veya ağ kontrolü: h küçüldüğünde ya da nokta sayısı arttığında sonucun kararlı bir değere yaklaşıp yaklaşmadığını incele.
Duyarlılık kontrolü: Başlangıç tahmini veya giriş verisindeki küçük değişimin sonucu orantısız büyütmediğini doğrula.
Bağımsız kontrol: Mümkünse düşük mertebeli yaklaşık hesap, denge eşitliği veya farklı yöntemle sonucu karşılaştır.

Yakınsama ile doğruluk aynı değildir. Yanlış kurulmuş modele veya yanlış denkleme uygulanan algoritma çok düzenli yakınsasa bile fiziksel olarak yanlış bir sonuca ulaşabilir. Önce model, sonra algoritma ve son olarak uygulama doğrulanır.

Sık Yapılan Hatalar

[a,b] integralinde standart apsisleri dönüşümsüz kullanmak.
Jacobian çarpanını unutmak.
Apsisleri eşit aralıklı noktalar sanmak.

Teknik Görsel Nasıl Okunur?

Şemadaki mavi öğeler bilinen verileri, kırmızı öğeler hesaplanan veya güncellenen büyüklükleri, yeşil oklar ise algoritmanın ilerleme yönünü gösterir. Grafik üzerindeki noktalar tam ölçekli ölçüm değil, yöntemin geometrik anlamını açıklayan eğitim gösterimidir.

Uygulama Kontrol Listesi

Problem, değişkenler ve birimler açık mı?
Kullanılan formülün varsayımları sağlanıyor mu?
Başlangıç değeri, adım büyüklüğü ve tolerans raporlandı mı?
Ara hesaplar erken yuvarlanmadan taşındı mı?
Sonuçta artık, hata veya ağ yakınsaması kontrol edildi mi?
Bulunan değer fiziksel sınırlar içinde mi?

Özet

Gauss-Legendre Kuadratürü, yaklaşık çözümü denetlenebilir bir algoritmaya dönüştürür. Güvenilir sayısal analiz; doğru formül, açık durdurma ölçütü, hata değerlendirmesi ve mühendislik yorumu birlikte verildiğinde tamamlanır.