Tiga Metode, Tiga Daerah
Untuk persamaan uji dy/dx = λy, tiga metode ODE eksplisit memiliki daerah stabilitas berikut dalam bidang hλ kompleks:
Metode Euler (orde pertama): daerah stabilitas adalah disk |1 + hλ| ≤ 1, lingkaran berjari-jari 1 berpusat di (-1, 0). Nilai hλ negatif real harus berada di [-2, 0].
Runge-Kutta 2 (metode titik tengah) (orde kedua): daerah stabilitas adalah |1 + hλ + (hλ)²/2| ≤ 1. Lebih besar dari disk Euler, tetapi masih terbatas.
Runge-Kutta 4 (orde keempat): daerah stabilitas memenuhi |1 + hλ + (hλ)²/2 + (hλ)³/6 + (hλ)⁴/24| ≤ 1. Nilai hλ negatif real mencapai sekitar -2.785. Daerah ini jauh lebih besar daripada Euler.
Backward Euler (implisit): daerah stabilitas adalah seluruh bidang kompleks kecuali disk |1 - hλ|⁻¹ > 1, setara dengan |1/(1-hλ)| ≤ 1. Untuk λ dalam setengah bidang kiri (Re(λ) < 0), ini secara bersyarat stabil — tidak ada batasan pada h dari stabilitas.
Fungsi Amplifikasi
Untuk setiap metode Runge-Kutta, faktor amplifikasi per langkah R(hλ) adalah aproksimasi polinomial untuk 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 sejati: |e^z| = e^(Re(z)). Untuk Re(z) < 0 (ODE stabil), solusi sejati meluruh. Metode numerik stabil jika |R(z)| ≤ 1 — cocok dengan perilaku meluruh.
Nilai Eigen Imajiner Murni: Sistem Osilatif
Banyak sistem fisik memiliki nilai eigen imajiner murni: λ = iω (osilasi tanpa redaman). Sistem pegas-massa, mekanika orbital, dinamika pendulum.
Untuk λ = iω: hλ = ihω terletak pada sumbu imajiner.
Stabilitas Euler pada sumbu imajiner: |1 + ihω|² = 1 + (hω)² > 1 untuk setiap h > 0. Euler tidak stabil untuk ukuran langkah apa pun pada nilai eigen imajiner murni. 'Osilasi' yang dihitung tumbuh tanpa batas.
Stabilitas RK4 pada sumbu imajiner: daerah stabilitas meluas ke sekitar |hω| ≤ 2.83 pada sumbu imajiner. Untuk h yang cukup kecil, RK4 menangani osilasi tanpa redaman. Euler tidak bisa.
Geometri ini menjelaskan mengapa Euler gagal pada sistem konservatif (pegas-massa, orbit, persamaan gelombang) bahkan dengan h kecil, sementara RK4 menanganinya dengan baik.
Geometri Masalah Kaku
Sistem ODE kaku memiliki nilai eigen dengan besarnya sangat berbeda. Rasio kekakuan: κ = max|Re(λᵢ)| / min|Re(λᵢ)| >> 1.
Mengapa kekakuan mahal untuk pemecah eksplisit:
Stabilitas memerlukan h·max|λᵢ| ≤ C (di mana C bergantung pada metode). Nilai eigen paling negatif menetapkan batasnya.
Keakuratan untuk dinamika lambat memerlukan h·min|λᵢ| ≥ ε (menyelesaikan mode paling lambat dengan cukup).
Jika κ besar, kedua persyaratan ini memaksa h yang sangat kecil: cukup kecil untuk stabilitas mode cepat, cukup besar untuk mencontohkan mode lambat. Jumlah langkah skala dengan κ.
Gambaran geometris dalam spektrum nilai eigen: nilai eigen Jacobian ∂f/∂y membentuk sekumpulan titik dalam bidang kompleks. Daerah stabilitas pemecah eksplisit harus berisi semua titik h·λᵢ. Jika nilai eigen berkisar dari -1 hingga -1000, daerah stabilitas harus mencakup rentang 1000 sepanjang sumbu real — memerlukan h yang sangat kecil.
Pemecah implisit: daerah stabilitas backward Euler mencakup seluruh setengah bidang kiri. Semua nilai eigen dengan Re(λ) < 0 secara otomatis berada di dalam daerah stabilitas terlepas dari h. Batasan pada h hanya berasal dari keakuratan, bukan stabilitas.
Rasio Kekakuan & Biaya
Pertimbangkan jaringan reaksi kimia dengan reaksi cepat (skala waktu 10⁻⁶ s) dan reaksi lambat (skala waktu 1 s).
Rasio kekakuan: κ = 10⁶ / 1 = 10⁶.
Dengan RK4 (batas stabilitas h·|λ_max| ≤ 2.785): h_max = 2.785 / 10⁶ ≈ 2,8 × 10⁻⁶ s.
Untuk mengintegrasikan selama 10 s waktu reaksi: langkah = 10 / (2,8 × 10⁻⁶) ≈ 3,6 × 10⁶.
Dengan backward Euler (tidak bersyarat stabil): h dapat dipilih untuk keakuratan reaksi lambat. h = 10⁻² s (100 sampel di atas skala 1 s). Langkah = 10 / 10⁻² = 1000.
Rasio biaya: eksplisit 3,6 juta langkah vs implisit 1000 langkah — faktor 3600. Setiap langkah implisit memerlukan penyelesaian sistem linier (biaya per langkah lebih tinggi), tetapi biaya total jauh lebih rendah untuk masalah yang sangat kaku.
Mengapa Tabung n-Dimensi Bukan Apa yang Anda Pikirkan
Dalam 2D, 'tabung' dengan jari-jari ε di sekitar kurva C adalah himpunan titik dalam jarak ε dari C. Penampang silang adalah lingkaran dengan jari-jari ε. Volume tabung tumbuh sebanding dengan panjangnya.
Dalam n dimensi, geometri tabung berubah secara fundamental, karena fenomena dari Bab 9:
Paradoks sudut berdimensi-n: dalam ruang berdimensi-n, hampir semua volume hiperkubus berdimensi-n terletak di sudut — bukan di wilayah pusat. Ketika n meningkat, fraksi volume dalam jarak ε dari pusat mendekati nol untuk setiap ε yang tetap.
Diterapkan pada tabung solusi ODE:
Dalam 2D: jika solusi sejati melewati pusat tabung, sebagian besar titik di dekatnya dekat dengan kurva. Gangguan kecil membuat Anda tetap dekat dengan solusi sejati.
Dalam dimensi tinggi: sebagian besar titik dalam kotak batas tabung sebenarnya jauh dari kurva solusi sejati. 'Volume' tabung didominasi oleh sudut — wilayah yang jauh dari pusat dalam beberapa dimensi secara bersamaan.
Konsekuensi untuk simulasi: dengan 28 ODE berpasangan (masalah intersepsi Angkatan Laut Hamming), gangguan ukuran ε dalam setiap dimensi dapat menghasilkan perpindahan total ε√28 ≈ 5,3ε dari solusi sejati. Tabung harus dipahami dalam hal norma L2 di semua dimensi, bukan hanya perpindahan maksimum dalam dimensi apa pun.
Stabilitas dalam dimensi tinggi: sistem di mana setiap komponen meluruh secara independen (setiap nilai eigen memiliki bagian real negatif) mungkin masih menunjukkan perpindahan gabungan yang besar karena kesalahan komponen menambah dalam norma L2. Tabung 28-dimensi bukan hanya 28 tabung 1-dimensi independen — geometri menggandengnya.
Dari Geometri ke Desain
Wawasan geometris dari Bab 18-20 bersatu sebagai serangkaian prinsip desain untuk simulasi numerik:
Pemilihan ukuran langkah: h harus menempatkan h·λ di dalam daerah stabilitas untuk setiap nilai eigen. Untuk sistem kaku, metode implisit menghilangkan batasan stabilitas, meninggalkan hanya persyaratan keakuratan.
Akumulasi kesalahan dalam dimensi tinggi: kesalahan global adalah vektor dalam ruang berdimensi-n. Normanya tumbuh sebagai √n kali kesalahan per-komponen. Simulasi berdimensi tinggi membutuhkan persyaratan keakuratan per-langkah yang lebih ketat.
Umpan balik sebagai penstabil: jika simulasi menggabungkan umpan balik (output yang dihitung mempengaruhi input berikutnya, seperti dalam sistem panduan), umpan balik konvergen meredam kesalahan. Simulasi dapat menoleransi input yang tidak tepat untuk jumlah di dalam loop umpan balik.
Instabilitas sebagai sinyal: untuk masalah dengan bidang arah divergen, instabilitas dapat dimanfaatkan: arah divergensi membawa informasi tentang kesalahan kondisi awal, memungkinkan penyesuaian korektif.