Tiga Metode, Tiga Daerah
Untuk persamaan uji dy/dx = λy, tiga metode ODE eksplisit memiliki daerah stabilitas berikut dalam ruang kompleks hλ:
Metode Euler (pertama): daerah stabilitas adalah lingkaran |1 + hλ| ≤ 1, sebuah lingkaran dengan jari-jari 1 yang berpusat di (-1, 0). hλ negatif sebenarnya harus berada di [-2, 0].
Runge-Kutta 2 (metode tengah) (kedua): daerah stabilitas adalah |1 + hλ + (hλ)²/2| ≤ 1. Lebih besar dari lingkaran Euler, tetapi masih terbatas.
Runge-Kutta 4 (keempat): daerah stabilitas memenuhi |1 + hλ + (hλ)²/2 + (hλ)³/6 + (hλ)⁴/24| ≤ 1. hλ negatif sebenarnya sekitar -2.785. Daerah tersebut jauh lebih besar daripada Euler.
Backward Euler (implisit): daerah stabilitas adalah seluruh ruang kompleks kecuali lingkaran |1 - hλ|⁻¹ > 1, atau setara dengan |1/(1-hλ)| ≤ 1. Untuk λ di setengah ruang kompleks (Re(λ) < 0), ini stabil secara unconditional — tidak ada batasan pada h dari segi stabilitas.
Fungsi Amplifikasi
Untuk setiap metode Runge-Kutta, faktor amplifikasi per langkah R(hλ) adalah pendekatan polinomial terhadap e^(hλ):
- Euler: R(z) = 1 + z (dipotong pada derajat 1)
- RK2: R(z) = 1 + z + z²/2 (dipotong pada derajat 2)
- RK4: R(z) = 1 + z + z²/2 + z³/6 + z⁴/24 (dipotong pada derajat 4)
Daerah stabilitas adalah {z ∈ ℂ : |R(z)| ≤ 1}. Amplifikasi solusi sebenarnya: |e^z| = e^(Re(z)). Untuk Re(z) < 0 (ODE stabil), solusi sebenarnya menurun. Metode numerik stabil jika |R(z)| ≤ 1 — mencocokkan perilaku menurun.
Nilai Eigen Hanya Imaginer: Sistem Osilasi
Banyak sistem fisika yang memiliki nilai-nilai eigen hanya imaginernya: λ = iω (osilasi tanpa pengurangan). Sistem pegas-massa, mekanika orbital, dinamika lonceng.
Untuk λ = iω: hλ = ihω berada di sumbu imaginernya.
Stabilitas Euler pada sumbu imajiner: |1 + ihω|² = 1 + (hω)² > 1 untuk setiap h > 0. Euler tidak stabil untuk setiap ukuran langkah pada eigenvalue yang sepenuhnya imajiner. 'Gesekan' yang dihitung tumbuh tanpa batas.
Stabilitas RK4 pada sumbu imajiner: daerah kestabilan mencakup sekitar |hω| ≤ 2.83 pada sumbu imajiner. Untuk h cukup kecil, RK4 mengatasi gelombang yang tidak diredam. Euler tidak bisa.
Geometri ini adalah mengapa Euler gagal pada sistem konservatif (spring-mass, orbit, persamaan gelombang) meskipun dengan h kecil, sementara RK4 mengatasi mereka dengan baik.
Geometri Masalah Keras
Sistem ODE keras memiliki eigenvalue dengan skala besar yang sangat berbeda. Rasio kekakuan: κ = max|Re(λᵢ)| / min|Re(λᵢ)| >> 1.
Mengapa kekakuan mahal untuk pemecah eksplisit:
Kebutuhan kestabilan adalah h·max|λᵢ| ≤ C (di mana C tergantung pada metode). Eigenvalue yang paling negatif yang menentukan batas.
Akurasi untuk dinamika lambat membutuhkan h·min|λᵢ| ≥ ε (menyelesaikan mode yang paling lambat dengan baik).
Jika κ besar, kedua persyaratan ini memaksa h kecil: cukup kecil untuk kestabilan mode cepat, cukup besar untuk mengambil sampel mode lambat. Jumlah langkah berkorelasi dengan κ.
Gambar geometri di spektrum eigenvalue: eigenvalue dari Jacobian ∂f/∂y membentuk sebuah himpunan titik di ruang kompleks. Daerah kestabilan solver eksplisit harus mencakup semua titik h·λᵢ. Jika eigenvalue mencakup rentang dari -1 hingga -1000, daerah kestabilan harus menutupi rentang 1000 di sumbu real — membutuhkan h sangat kecil.
Solver implisit: daerah kestabilan Euler belakang menutupi seluruh setengah ruang kiri. Semua eigenvalue dengan Re(λ) < 0 secara otomatis berada di dalam daerah kestabilan, terlepas dari h. Ketentuan h hanya datang dari akurasitas, bukan kestabilan.
Rasio Ketegangan & Biaya
Sertakan sebuah jaringan reaksi kimia dengan reaksi cepat (skala waktu 10⁻⁶ s) dan reaksi lambat (skala waktu 1 s).
Rasio ketegangan: κ = 10⁶ / 1 = 10⁶.
Dengan RK4 (batas kestabilan h·|λ_max| ≤ 2.785): h_max = 2.785 / 10⁶ ≈ 2,8 × 10⁻⁶ s.
Untuk mengintegrasikan selama 10 s waktu reaksi: langkah-langkah = 10 / (2,8 × 10⁻⁶) ≈ 3,6 × 10⁶.
Dengan Euler belakang (kestabilan unconditional): h dapat dipilih untuk akurasitas reaksi lambat. h = 10⁻² s (100 sampel selama 1 s skala). Langkah-langkah = 10 / 10⁻² = 1000.
Rasio biaya: 3,6 juta langkah eksplisit vs 1000 langkah implisit — faktor 3600. Setiap langkah implisit membutuhkan penyelesaian sistem linear (biaya per langkah lebih tinggi), tetapi total biaya jauh lebih rendah untuk masalah sangat kaku.
Mengapa Tabung Berdimensi n Bukan Sesuatu yang Anda Pikirkan
Dalam 2D, 'tabung' dengan jari-jari ε di sekitar kurva C adalah himpunan titik dalam jarak ε dari C. Silinder potongannya adalah lingkaran dengan jari-jari ε. Volume tabung berkembang proporsional dengan panjangnya.
Dalam n dimensi, geometri tabung berubah secara fundamental, karena fenomena dari Bab 9:
Paradoks sudut n-dimensional: di ruang n-dimensional, hampir seluruh volume kubus n-dimensional terletak di sudut-sudut - bukan di region tengah. Seiring bertambahnya n, fraksi volume dalam jarak ε dari pusat turun ke nol untuk ε tetap.
Diterapkan pada tabung solusi ODE:
Dalam 2D: jika solusi sebenarnya melewati pusat tabung, sebagian besar titik tetangga dekat dengan kurva. Gangguan kecil menjaga Anda tetap dekat dengan solusi sebenarnya.
Dalam dimensi tinggi: sebagian besar titik dalam kotak batas tabung sebenarnya jauh dari kurva solusi sebenarnya. Volume tabung didominasi oleh sudut-sudut - region yang jauh dari pusat dalam dimensi beberapa dimensi secara bersamaan.
Konsekuensi untuk simulasi: dengan 28 persamaan ODE yang berkaitan (masalah intersepsi Angkatan Laut Hamming), gangguan ukuran ε dalam setiap dimensi dapat menghasilkan perpindahan total ε√28 ≈ 5,3ε dari solusi sebenarnya. Tabung harus dipahami dalam hal norma L2 di seluruh dimensi, bukan hanya perpindahan maksimum dalam satu dimensi.
Stabilitas dalam dimensi tinggi: suatu sistem di mana setiap komponen mengalami penurunan independen (setiap nilai eigen memiliki bagian real negatif) mungkin menunjukkan penurunan kombinasi besar karena kesalahan komponen yang menambahkan norma L2. Tabung 28-dimensional bukan hanya 28 tabung 1-dimensional yang independen - geometri mengkaitkan mereka.
Dari Geometri ke Desain
Insight geometri Bab 18-20 bergabung menjadi sebuah set prinsip desain untuk simulasi numerik:
Pemilihan langkah: h harus menempatkan h·λ di dalam daerah stabil untuk setiap nilai eigen. Untuk sistem yang kaku, metode implisit menghapus ketentuan kestabilan, meninggalkan hanya persyaratan akurasi.
Akumulasi kesalahan dalam dimensi tinggi: kesalahan global adalah vektor dalam ruang n-dimensional. Normnya tumbuh sebesar √n kali kesalahan per-komponen. Simulasi berdimensi tinggi membutuhkan persyaratan ketelitian per-langkah yang lebih ketat.
Balikkan sebagai stabilizer: jika simulasi menggabungkan umpan balik (output yang dihitung mempengaruhi input berikutnya, seperti dalam sistem panduan), umpan balik konvergen menurunkan kesalahan. Simulasi dapat tolerir input yang tidak akurat untuk jumlah dalam lingkaran umpan balik.
Instabilitas sebagai sinyal: untuk masalah dengan bidang arah divergen, instabilitas dapat diambil alih: arah divergensi mengangkut informasi tentang kesalahan kondisi awal, memungkinkan penyesuaian korektif.