Dikdörtgen Hesaplamaları
Geometrik Fikir
Bir sıradan diferansiyel denklem dy/dx = f(x,y) her nokta (x,y) düzleminde bir eğim atar - bir yön alanı. Gerçek çözüm y(x) her zaman o atanan eğimleri takip eden bir eğri olmalıdır.
Euler'ın yöntemi sürekli yön alanını diskrete bir yürüyüşe dönüştürür:
> (xₙ, yₙ) → (xₙ + h, yₙ + h·f(xₙ, yₙ))
(xₙ, yₙ) noktasından, h mesafesine tangent düzlemin üzerinde hareket et. Yaklaşık bir sonraki noktaya var. Tekrarla.
Geometrik hata: tangent düzlemin (xₙ, yₙ) noktasındaki eğim f(xₙ, yₙ) olsa da, gerçek eğrinin her [xₙ, xₙ + h] aralığında noktasında farklı bir eğim vardır. Euler adımının sol uçtaki eğimleri kullanmasıdır - 'geçmişteki eğim.' Her adım hata büyür ve h² ile orantılıdır.
Biriktirilmiş Hata
N adımlık sabit bir uç nokta x = a'ya ulaşmak için, h = a/N:
- Yerel kesme hata: O(h²)
- Adım sayısı: N = a/h
- Toplamlı hata: O(h²) × (a/h) = O(h) - birinci dereceden doğruluk
Euler'ın yöntemi birinci dereceden doğruluktur: h'yi ikiye katlarsanız, toplam hata da ikiye katlanır.
Euler'ın Yöntemini Çalıştırma
dy/dx = y, başlangıç koşulu y(0) = 1. Gerçek çözüm: y(x) = eˣ, böylece y(1) = e ≈ 2.71828.
Euler'ın yöntemini h = 0.5 ile x = 0'dan x = 1'e (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. Hata: 0.468. Oransal hata: ~17%.
Euler'un Stabilite Bölgesi Derlemesi
Test denklemi dy/dx = λy (burada λ karmaşık bir sayıdır) için Euler'ın yöntemi verecek:
> yₙ₊₁ = yₙ + h·λ·yₙ = yₙ·(1 + hλ)
Adım başına büyüme faktörü: z = 1 + hλ.
Stabilite koşulu: hesaplanan çözümün sınırsız kalması için ve sadece |z| ≤ 1, yani |1 + hλ| ≤ 1 olması gerekmektedir.
Bu, karmaşık hλ düzleminde bir geometrik koşuldur: hλ noktası, (-1, 0) merkezli ve yarıçap 1 olan daire içinde olmalıdır.
Euler'un stabilite bölgesi: { hλ ∈ ℂ : |1 + hλ| ≤ 1 }
Gerçek, negatif λ için (azalan ODE gibi) dy/dx = -2y: hλ, gerçek eksen üzerinde (-2, 0] aralığında olmalıdır. λ = -2 ve h = 0.5 ile: hλ = -1. Bu, metodun niceliksel başarısızlık örneğinde niceliksel kararlılık açıklar - tam olarak stabilite sınırında. h = 1 ve λ = -2 ile: hλ = -2, stabilite bölgesinden dışarıda. Çözüm büyüyen amplitüde ile sallanır.
Stabilite Sınırını Bulma
Runge-Kutta 4 (RK4), Euler'den daha büyük bir stabilite bölgesi sağlar, bu da genellikle problemler için tercih edilmesinin bir nedenidır.
Gerçek negatif λ için, RK4 gerçek eksen üzerinde -2.785'a kadar hλ'yi sağlar (Euler'un -2 sınırlayıcı limiti).
Eigen değerler λ arasında çok farklı büyüklüklerde olan kasıtsız denklemeler için (ör. λ₁ = -1 ve λ₂ = -1000) stabilite, hλ₂'nin stabilite bölgesinde kalması gerektiği anlamına gelir. Gerçek eksen üzerinde RK4 için: h·(-1000) ≥ -2.785, böylece h ≤ 0.002785.
Bu küçük adım boyutu, yavaş bileşen λ₁ için h = 2 kullanabileceği rağmen, kasıtsız denklemeler için hλ₂'nin stabil olması gerektiği için pahalı bir simülasyon sağlar.
Sabit Noktalar & Çekirdek Alanları
Euler'ın yöntemi, dy/dx = f(y) diferansiyel denklemini çözerek diskret bir harita oluşturur: yₙ₊₁ = g(yₙ) = yₙ + h·f(yₙ).
Bu haritanın sabit noktası: y such that g(y) = y. Euler'ın dy/dx = f(y) üzerindeki uygulamasında, diferansiyel denklemin sabit noktaları f(y) = 0'u karşılar - diferansiyel denklemin dengi.
Sabit nokta istikrarı: eğer |g'(y)| < 1, y etrafındaki iteratlar y'ye doğru yaklaşırlar. Eğer |g'(y)| > 1, uzaklaşırlar.
g'(y) = 1 + h·f'(y). Sabit nokta y üzerinde: |1 + h·f'(y)| < 1 için istikrar.
Bu, ODE'nin lineerleştirilmesinin denklemi olan λ = f'(y*) ile Euler istikrar koşulu ile tam olarak aynıdır.
Çekirdek alanı: Euler haritası üzerinde y*'ye doğru yaklaşılan başlangıç koşullarının kümesidir. Nonlinear sistemler için, simülasyonun ODE denklemi dengeyi izlemeye devam etmesi veya başka bir cazip nokta etrafında dalgalı olması sınırı belirler.
Simülasyon döngüsü, bir diskret dinamik sistemdir. Nitelikleri - yaklaşımlar, dalgalı hareket, uzaklaşma - h adımının ODE'nin yönlendirme alanının geometrisine göre değişir.
Cebirsel Şekillere Uygulanan Geometri
Numarik simülasyonun geometrisi üç soruya indirgenebilir:
1. Stabilite bölgesinin nerede? Euler için: disk |1 + hλ| ≤ 1. Büyük için RK4. İmplisit yöntemler için sınırsız (her iki yarıplanekteki sol kısmı kapsar).
2. ODE'nin eksen değerleri nerede? ODE'nin Jacobian'ın f'nin her noktasındaki eksen değerleri λ, hλ'nin hangi stabilite bölgesinde olması gerektiğini belirler.
3. h'nin tüm hλ'yi bölgede tutmak için ne olmalı? En izin verilebilir h = (stabilite bölgesinin yarıçapları) / max|λ|.
Sert sistemler için: max|λ| devasa, açık yöntemler için küçük h'yi zorlar. İmplisit yöntemler her adımda maliyetli olsa da büyük h'ye izin verir — onlar, stabilite için hızı adım maliyeti için değiştirir.
Hamming'in keşfi çevirir: sayısal metodun seçimi, ODE'nin eksen değerleri spektrumunun geometrisi hakkında bir bahse dönüşür. Bu bahsi açık hale getirmek, herhangi bir simülasyonun tasarım kararının ilk aşamasıdır.