Teğet Doğru Yaklaşıklığı
Geometrik İdea
Sıradan bir diferansiyel denklem dy/dx = f(x,y) (x,y) düzlemindeki her noktaya bir eğim atar — bir yön alanı. Gerçek çözüm y(x), bu atanmış eğimleri her yerde izleyen bir eğridir.
Euler'in yöntemi, sürekli yön alanını ayrık bir yürüyüşe dönüştürür:
> (xₙ, yₙ) → (xₙ + h, yₙ + h·f(xₙ, yₙ))
(xₙ, yₙ) noktasından, teğet doğru boyunca h mesafesi hareket edin. Yaklaşık bir sonraki noktaya ulaşın. Tekrarlayın.
Geometrik hata: (xₙ, yₙ) noktasındaki teğet doğrunun eğimi f(xₙ, yₙ)'dir, ancak gerçek eğri [xₙ, xₙ + h] aralığı boyunca her noktada farklı bir eğime sahiptir. Euler adımı, sol uç noktadaki eğimi baştan sona kullanır — 'olan eğim.' Adım başına hata h² olarak büyür.
Birikmiş Hata
Sabit bir uç nokta x = a'ya ulaşmak için N adımda, h = a/N ile:
- Adım başına yerel kesme hatası: O(h²)
- Adım sayısı: N = a/h
- Küresel hata: O(h²) × (a/h) = O(h) — birinci dereceden doğruluk
Euler'in yöntemi birinci mertebedendir: h'yi yarıya indirmek küresel hatayı yarıya indirir.
Euler'in Yöntemi Çalıştırma
dy/dx = y'yi, başlangıç koşulu y(0) = 1 ile düşünün. Gerçek çözüm: y(x) = eˣ, yani y(1) = e ≈ 2.71828.
Euler'in yöntemi h = 0.5 ile, x = 0'dan x = 1'e kadar (2 adım) uygulayın:
Adım 1: y₁ = y₀ + h·f(x₀, y₀) = 1 + 0.5·(1) = 1.5. Yeni nokta: (0.5, 1.5).
Adım 2: y₂ = y₁ + h·f(x₁, y₁) = 1.5 + 0.5·(1.5) = 2.25. Yeni nokta: (1.0, 2.25).
Euler 2.25 verirken gerçek değer 2.71828'dir. Hata: 0.468. Bağıl hata: ~17%.
Euler'in Stabilite Bölgesini Türetme
Test denklemi dy/dx = λy için (λ bir karmaşık sayı olmak üzere), Euler'in yöntemi şunları verir:
> yₙ₊₁ = yₙ + h·λ·yₙ = yₙ·(1 + hλ)
Adım başına amplifikasyon faktörü: z = 1 + hλ.
Stabilite koşulu: hesaplanan çözüm, ancak ve ancak |z| ≤ 1, yani |1 + hλ| ≤ 1 ise sınırlı kalır.
Bu, karmaşık hλ düzleminde geometrik bir koşuldur: hλ noktası (-1, 0) merkezli yarıçapı 1 olan daire içinde yer almalıdır.
Euler'in stabilite bölgesi: { hλ ∈ ℂ : |1 + hλ| ≤ 1 }
Gerçek, negatif λ için (dy/dx = -2y gibi azalan bir ODE): hλ, reel eksen üzerinde (-2, 0] aralığında yer almalıdır. λ = -2 ve h = 0.5 ile: hλ = -1. Bu, stabilite sınırının tam üzerindedir — yöntem marjinal olarak stabildir, bu da önceki örnekteki niteliksel başarısızlığı açıklar.
h = 1 ve λ = -2 ile: hλ = -2, bizi stabilite bölgesinin dışına koyar. Çözüm büyüyen amplitüttle salınır.
Stabilite Sınırını Bulma
Runge-Kutta 4 (RK4), Euler'den daha geniş bir stabilite bölgesine sahiptir, bu da çoğu problem için tercih edilme nedenlerinden biridir.
Gerçek negatif λ için, RK4 hλ'yı reel eksen üzerinde yaklaşık -2.785'e kadar (Euler'in -2 sınırına karşı) izin verir.
Özdeğerleri çok farklı büyüklüklerde olan rijit denklemler için — say λ₁ = -1 ve λ₂ = -1000 — stabilite hλ₂'nin bölge içinde kalmasını gerektirir. RK4 için reel eksen üzerinde: h·(-1000) ≥ -2.785, yani h ≤ 0.002785.
Bu küçük adım boyutu, rijit özdeğer λ₂ tarafından belirlenmiş, yavaş bileşen λ₁ h = 2 kullanabilmesine rağmen simülasyonu pahalı hale getirir.
Sabit Noktalar & Çekim Havzaları
Euler'in yöntemi dy/dx = f(y)'ye uygulandığında, ayrık bir harita tanımlar: yₙ₊₁ = g(yₙ) = yₙ + h·f(yₙ).
Bu haritanın bir sabit noktası: g(y) = y olacak şekilde y. Euler dy/dx = f(y) üzerinde için, sabit noktalar f(y) = 0 sağlar — ODE'nin dengeli noktaları.
Sabit noktanın stabilitesi: eğer |g'(y)| < 1 ise, yakın yinelemeler y'e yakınsar. Eğer |g'(y*)| > 1 ise, ıraksarlar.
g'(y) = 1 + h·f'(y). Sabit bir nokta y'de: stabilite için |1 + h·f'(y)| < 1.
Bu, tam olarak λ = f'(y*) ile Euler stabilite koşuludur — denge noktasında ODE'nin doğrusallaştırılması.
Çekim havzası: Euler haritası altında y*'e yakınsamayan başlangıç koşullarının kümesi. Doğrusal olmayan sistemler için, havza sınırı simülasyonun ODE dengesini güvenilir bir şekilde izleyeceği yeri vs başka bir çekiciye sapacağını tanımlar.
Simülasyon döngüsü ayrık dinamik bir sistemdir. Niteliksel davranışı — yakınsama, salınım, sapma — adım boyutu h'ye ODE'nin yön alanının geometrisine göre bağlıdır.
Geometriyi Simülasyon Tasarımına Bağlama
Sayısal simülasyonun geometrisi üç soruya iner:
1. Stabilite bölgesi nerede? Euler için: disk |1 + hλ| ≤ 1. RK4 için daha geniş. İmplisit yöntemler için sınırsız (tüm sol yarı-düzlem).
2. ODE'nin özdeğerleri nerede? Her noktada f'nin Jacobian'ının özdeğerleri λ, hangi stabilite bölgesinin hλ içermesi gerektiğini belirler.
3. Hangi h tüm hλ'yı bölge içinde tutmuş? Maksimum izin verilen h = (stabilite bölgesi yarıçapı) / max|λ|.
Rijit sistemler için: max|λ| muasaldır, açık yöntemler için minuscul h'yi zorlar. İmplisit yöntemler adım başına pahalı ancak büyük h'ye izin verir — adım başına maliyeti stabilite için değişirler.
Hamming'in içgörüsü çevirisi: sayısal yöntem seçimi, ODE'nin özdeğer spektrumunun geometrisi hakkında bir bahis kodlar. Bu bahsi açık hale getirmek herhangi bir simülasyondaki ilk tasarım kararıdır.