b belajar.myudak

Mencari akar persamaan nonlinier

Metode tertutup (biseksi) dan terbuka (Newton–Raphson, secant), galat, kriteria berhenti, dan laju konvergensi.

dasar #metode numerik#akar#iterasi#galat

Banyak persamaan tidak bisa diselesaikan secara analitik, misalnya cos⁡x=x\cos x = x atau ex=3xe^x = 3x. Metode numerik mencari hampiran akarnya, yaitu nilai xx yang membuat f(x)=0f(x) = 0, lewat iterasi sampai cukup teliti.

Metode tertutup: biseksi

Jika ff kontinu dan f(a)⋅f(b)<0f(a) \cdot f(b) < 0, maka (Teorema Nilai Antara) ada akar di [a,b][a, b]. Biseksi terus membagi dua selang:

  1. c=(a+b)/2c = (a + b)/2
  2. Jika f(a)⋅f(c)<0f(a) \cdot f(c) < 0, akar di [a,c][a, c], sehingga b←cb \leftarrow c. Jika tidak, a←ca \leftarrow c.
  3. Ulangi sampai (b−a)/2<ε(b - a)/2 < \varepsilon.

Biseksi selalu konvergen, dan jumlah iterasinya bisa dihitung di muka: n≥log⁡2 ⁣(b−aε)n \ge \log_2\!\left(\frac{b - a}{\varepsilon}\right). Kekurangannya, ia lambat (konvergensi linier).

Metode terbuka: Newton–Raphson & secant

Newton–Raphson memakai garis singgung: tarik garis singgung di xix_i, potongkan dengan sumbu-x, lalu jadikan titik potong itu tebakan berikutnya.

xi+1=xi−f(xi)f′(xi)x_{i+1} = x_i - \frac{f(x_i)}{f'(x_i)}

Konvergensinya kuadratik: jumlah digit benar kira-kira berlipat dua setiap iterasi. Tapi butuh turunan f′f', dan bisa divergen jika tebakan awal jauh atau f′(xi)≈0f'(x_i) \approx 0.

Secant mengganti turunan dengan kemiringan dari dua titik terakhir, sehingga tidak perlu menghitung f′f':

xi+1=xi−f(xi) xi−xi−1f(xi)−f(xi−1)x_{i+1} = x_i - f(x_i)\,\frac{x_i - x_{i-1}}{f(x_i) - f(x_{i-1})}
Mencari akar f(x) = 0
iterasi 1/20 · x ≈ 1.50000000
iabxᵢf(xᵢ)galat
11.0000002.0000001.50000000-1.25e-15.00e-1

Toleransi 10⁻⁶. Biseksi selalu konvergen tapi linier (galat setengahnya tiap iterasi); Newton konvergen kuadratik (digit benar kira-kira berlipat dua) asal tebakan awal dekat dan f′(x) ≠ 0.

Bandingkan jumlah iterasi ketiga metode untuk fungsi yang sama dengan toleransi 10−610^{-6}. Newton biasanya selesai dalam 4–6 iterasi, sementara biseksi butuh sekitar 20.

Galat dan kriteria berhenti

  • Galat absolut: ∣xi+1−xi∣|x_{i+1} - x_i|
  • Galat relatif: ∣xi+1−xixi+1∣\left|\frac{x_{i+1} - x_i}{x_{i+1}}\right|
  • Berhenti jika galat <ε< \varepsilon, atau ∣f(xi)∣<ε|f(x_i)| < \varepsilon, atau iterasi melebihi batas maksimum (wajib ada, untuk berjaga-jaga jika divergen).
Newton–Raphson untuk √a
Cek pemahaman
  1. 1. Syarat awal metode biseksi pada selang [a, b]?

  2. 2. Berapa iterasi biseksi minimal untuk selang panjang 1 dan toleransi 10⁻³?

  3. 3. Kelemahan utama Newton–Raphson?