Lineer Olmayan Denklem Sistemlerinde Newton Yöntemi

Konunun Amacı

Birden çok lineer olmayan denklemi Jacobian matrisiyle yerel olarak lineerleştirip düzeltme vektörü üzerinden eşzamanlı çözme. 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?

Çok bilinmeyenli denge, akım ağı ve malzeme durum denklemleri Newton lineerleştirmesiyle çözülü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. Fonksiyon vektörü

Tanım: Bütün doğrusal olmayan denklemlerin artıklarını bir arada tutan F vektörü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. Jacobian matrisi

Tanım: Fonksiyonların bilinmeyenlere göre birinci kısmi türevlerinden oluşan matristir.

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.

3. Düzeltme vektörü

Tanım: Mevcut tahmine eklenerek yeni tahmini oluşturan Δx vektörü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.

4. Artık normu

Tanım: Denklem takımının ne ölçüde sağlandığını tek sayıyla ölçen ||F|| değeridir.

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.

Temel Bağıntı

J(xᵢ)·Δxᵢ=−F(xᵢ); xᵢ₊₁=xᵢ+Δxᵢ

Her iterasyonda fonksiyon vektörü ve kısmi türevlerden Jacobian kurulur. Düzeltme, Jacobian'ı ters almak yerine lineer sistem çözülerek bulunur. 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ı

1. F vektörünü ve bilinmeyen sırasını tanımla.
2. Başlangıç vektöründe F ve Jacobian J'yi hesapla.
3. J·Δx=−F sistemini güvenilir bir lineer çözücüyle çöz.
4. Tahmini Δx kadar güncelle.
5. Vektör normu ve artık normu toleranslarını kontrol et.

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

J=[[2,1],[1,3]] ve F=[1,2] ise JΔ=−F çözülür. Denklem çözümünden Δx=−0,2 ve Δy=−0,6 bulunarak mevcut tahmin güncellenir.

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

• Jacobian'ı yanlış bilinmeyen sırasıyla kurmak.
• Her iterasyonda Jacobian ve artığı güncellememek.
• Tekil veya kötü koşullu Jacobian ile kontrolsüz büyük adım atmak.

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

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

Özet

Lineer Olmayan Denklem Sistemlerinde Newton Yöntemi, 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.