Penyelesai ODE Orde Kedua

Tentang Penyelesai ODE Orde Kedua

Penyelesai ODE Orde Kedua menerima persamaan diferensial biasa linier dalam bentuk a·y″ + b·y′ + c·y = g(x) dengan koefisien riil konstan, secara otomatis menurunkan persamaan karakteristiknya, mengklasifikasikan rezim redaman (overdamped, critically damped, underdamped, undamped, atau unstable), dan menghasilkan solusi simbolik bentuk tertutup serta solusi numerik akurasi tinggi. Output interaktif menyandingkan plot waktu kurva ganda y(x) dan y′(x) dengan lintasan bidang fase (y, y′) — tampilan yang menunjukkan rezim secara sekilas: spiral-masuk untuk underdamped, simpul-masuk untuk overdamped, loop tertutup untuk undamped, spiral-keluar untuk unstable.

Apa Itu ODE Linier Orde Kedua Dengan Koefisien Konstan?

Sebuah persamaan diferensial biasa linier orde kedua dengan koefisien konstanta riil adalah persamaan dalam bentuk

di mana a ≠ 0, b, c adalah konstanta riil dan g(x) adalah suku pemaksa. Dua kondisi awal y(x₀) = y₀ dan y′(x₀) = y′₀ mengubah ini menjadi masalah nilai awal dengan solusi unik pada lingkungan x₀ — hal ini mengikuti teorema Picard-Lindelöf yang diterapkan pada sistem orde pertama yang setara.

Jika g(x) = 0 persamaannya adalah homogen. Jika tidak, itu adalah non-homogen, dan solusi lengkapnya terurai sebagai

di mana y_h adalah solusi umum dari persamaan homogen terkait (berisi dua konstanta bebas) dan y_p adalah solusi partikular apa pun dari persamaan lengkap. Menerapkan dua kondisi awal akan menetapkan nilai kedua konstanta bebas tersebut.

Persamaan Karakteristik

Menebak y = e^(r·x) dalam persamaan homogen menghasilkan persamaan karakteristik (atau tambahan)

sebuah kuadrat yang diskriminannya Δ = b² − 4ac mengontrol seluruh perilaku kualitatif:

Tiga Kasus Akar & Rezim Redaman

Metode Koefisien Tak Tentu (Kasus Non-homogen)

Ketika g(x) mengambil salah satu dari bentuk sederhana berikut, metode koefisien tak tentu memberikan solusi partikular dengan mengasumsikan percobaan dari bentuk yang sama dengan koefisien yang tidak diketahui dan menyelesaikannya:

• Konstanta g(x) = k. Percobaan: y_p = K. Jika c = 0 kalikan dengan x; jika b = 0 juga, kalikan dengan x lagi.
• Polinomial derajat n. Percobaan: polinomial umum derajat n. Kalikan dengan x atau x² jika suku konstan atau linier beresonansi.
• Eksponensial g(x) = A·e^(k·x). Percobaan: y_p = K·e^(k·x). Jika k bertepatan dengan akar karakteristik, kalikan dengan x (akar sederhana) atau x² (akar ganda) — ini adalah resonansi.
• Sinusoidal g(x) = A·cos(ω·x) + B·sin(ω·x). Percobaan: y_p = K₁·cos(ω·x) + K₂·sin(ω·x). Kalikan dengan x jika iω adalah akar (resonansi frekuensi murni).
• Produk dan jumlah mengikuti prinsip linieritas dan aturan perkalian.

Membaca Bidang Fase

Sistem orde pertama yang setara adalah u = y, v = y′ dengan u′ = v dan v′ = (g(x) − b·v − c·u)/a. Memplot v terhadap u secara parametrik dalam x memberikan lintasan bidang fase. Untuk sistem otonom homogen (tidak ada x dalam g), orbit ditentukan secara unik oleh titik awalnya (y₀, y′₀) dan menunjukkan rezim secara sekilas:

• Underdamped: lintasan melingkar ke dalam menuju titik asal.
• Overdamped: lintasan mendekati titik asal sepanjang garis invarian (eigenvektor lambat).
• Critically damped: simpul degenerasi, lintasan bersinggungan dengan eigenvektor tunggal.
• Undamped: elips tertutup mengelilingi titik asal — osilasi abadi.
• Unstable: lintasan melingkar atau lari keluar menuju tak terhingga.

Contoh Pengerjaan: Osilator Harmonik Teredam Paksa

Pertimbangkan persamaan y″ + 2·y′ + 5·y = 10 dengan y(0) = 0, y′(0) = 0 — sebuah sistem underdamped yang dipaksa.

1. Persamaan karakteristik: r² + 2r + 5 = 0 → Δ = 4 − 20 = −16 → r = −1 ± 2i.
2. Solusi homogen: y_h = e^(−x)·(C₁·cos 2x + C₂·sin 2x).
3. Solusi partikular untuk pemaksa konstanta g = 10: coba y_p = K, sehingga 5K = 10, memberikan y_p = 2.
4. Terapkan IC: y(0) = 0 → C₁ + 2 = 0 → C₁ = −2. y′(0) = 0 → −C₁ + 2C₂ = 0 → C₂ = −1.
5. Jawaban akhir: y(x) = 2 − e^(−x)·(2·cos 2x + sin 2x) — berosilasi dengan amplop yang meluruh dan batas y → 2.

Cara Menggunakan Kalkulator Ini

1. Masukkan koefisien a, b, c di baris atas. a harus non-zero (jika tidak, persamaan tersebut adalah orde pertama).
2. Ketik suku pemaksa g(x), atau biarkan 0 untuk masalah homogen. Solusi partikular bentuk tertutup diturunkan untuk konstanta, polinomial hingga derajat 2, dan eksponensial tunggal A·e^(k·x) termasuk kasus resonansi.
3. Berikan kondisi awal (x₀, y₀, y′₀). Baik y maupun y′ pada x₀ harus ditentukan karena persamaannya adalah orde kedua.
4. Pilih rentang x untuk plot. Solver berintegrasi keluar dari x₀ ke kedua arah x menggunakan RK4.
5. Klik Selesaikan & Visualisasikan. Anda mendapatkan persamaan karakteristik dengan akar-akarnya pada bidang kompleks, klasifikasi rezim redaman, solusi bentuk tertutup homogen dan partikular, plot waktu kurva ganda y dan y′, serta lintasan bidang fase.

Aplikasi Umum

• Sistem pegas-massa-redaman mekanis: m·x″ + c·x′ + k·x = F(t). Overdamped, critically damped, dan underdamped sesuai dengan rasio redaman yang berbeda ζ = c/(2·√(m·k)).
• Sirkuit listrik RLC: sirkuit RLC seri mematuhi L·Q″ + R·Q′ + Q/C = V(t) — struktur identik, simbol berbeda.
• Pendulum (sudut kecil): θ″ + (g/L)·θ = 0 memberikan gerak harmonik sederhana; menambahkan hambatan udara memberikan osilasi teredam.
• Respon bangunan terhadap gempa bumi: struktur derajat-kebebasan-tunggal dengan percepatan dasar sebagai suku pemaksa.
• Sistem servo yang dikendalikan PID: dinamika kesalahan loop tertutup mereduksi menjadi ODE orde kedua yang rasio redamannya mengatur overshoot.
• Model populasi dengan inersia: pertumbuhan ekonomi dengan jeda akumulasi modal, atau model ekologis dengan respon yang tertunda.

Metode Numerik — Runge-Kutta Klasik (RK4) pada Sistem 2D

Alat ini mereduksi a·y″ + b·y′ + c·y = g(x) menjadi sistem orde pertama

dengan u(x₀) = y₀, v(x₀) = y′₀. Runge-Kutta empat tahap kemudian diterapkan pada status vektor (u, v). RK4 memiliki kesalahan pemotongan lokal O(h⁵) dan kesalahan global O(h⁴); default 400 sub-langkah di setiap arah memberikan akurasi sekitar enam digit untuk masalah non-kaku.

Pertanyaan yang Sering Diajukan

Apa itu ODE linier orde kedua dengan koefisien konstan?

ODE linier orde kedua dengan koefisien konstan memiliki bentuk a·y″ + b·y′ + c·y = g(x), di mana a, b, c adalah konstanta riil dan g(x) adalah suku pemaksa (non-homogen). Dengan dua kondisi awal y(x₀) = y₀ dan y′(x₀) = y′₀ solusinya adalah unik. Kasus homogen g(x) = 0 selalu memungkinkan solusi bentuk tertutup melalui persamaan karakteristik a·r² + b·r + c = 0; kasus non-homogen diselesaikan sebagai y(x) = y_h(x) + y_p(x).

Apa itu persamaan karakteristik?

Untuk a·y″ + b·y′ + c·y = 0, substitusi ansatz y = e^(r·x) menghasilkan a·r² + b·r + c = 0 — persamaan karakteristik atau tambahan. Akar-akarnya menentukan bentuk solusi homogen: dua akar riil berbeda memberikan y_h = C₁·e^(r₁·x) + C₂·e^(r₂·x); akar berulang r memberikan y_h = (C₁ + C₂·x)·e^(r·x); akar konjugasi kompleks α ± β·i memberikan y_h = e^(α·x)·(C₁·cos(β·x) + C₂·sin(β·x)).

Apa arti under-, critically, dan overdamped?

Terminologi ini berasal dari model pegas-massa-redaman m·x″ + c·x′ + k·x = 0. Overdamped (diskriminan > 0, dua akar riil) berarti sistem kembali ke keseimbangan perlahan tanpa osilasi. Critically damped (diskriminan = 0, akar berulang) adalah pengembalian tercepat tanpa overshoot. Underdamped (diskriminan < 0, akar kompleks) memberikan osilasi yang meluruh. Undamped (b = 0, c/a > 0) memberikan osilasi sinusoidal murni selamanya.

Apa itu metode koefisien tak tentu?

Untuk pemaksa sederhana g(x) — konstanta, polinomial, eksponensial, sinus, kosinus, dan produknya — solusi partikular y_p diasumsikan memiliki bentuk yang sama dengan g dengan koefisien yang tidak diketahui, yang ditentukan dengan mensubstitusi ke dalam ODE dan mencocokkan suku-sukunya. Percobaan harus dikalikan dengan x (atau x² untuk akar ganda) ketika g(x) beresonansi dengan akar karakteristik.

Apa itu bidang fase?

Untuk persamaan orde kedua yang direduksi menjadi sistem 2D (y, y') bidang fase memplot y' terhadap y seiring bertambahnya x. Kurva solusi dalam bidang fase menunjukkan rezim secara sekilas: spiral meluruh untuk underdamped, simpul masuk untuk overdamped, elips tertutup untuk gerak harmonik tak teredam, dan spiral keluar untuk osilasi tidak stabil. Ini adalah rekan geometris dari diagram akar persamaan karakteristik.

Metode numerik apa yang digunakan alat ini?

Metode Runge-Kutta orde keempat (RK4) klasik diterapkan pada sistem orde pertama yang setara u = y, v = y′, dengan u′ = v dan v′ = (g(x) − b·v − c·u)/a. RK4 memiliki kesalahan pemotongan lokal O(h⁵) dan default 400 sub-langkah per arah memberikan akurasi sekitar enam digit untuk persamaan non-kaku pada jendela yang dipilih.

Bacaan Lebih Lanjut

• Persamaan diferensial linier — Wikipedia
• Persamaan karakteristik — Wikipedia
• Metode koefisien tak tentu — Wikipedia
• Osilator harmonik — Wikipedia
• Bidang fase — Wikipedia
• Metode Runge-Kutta — Wikipedia