Iterācijas metode sakņu atrašanā
Lietišķajā matemātikā, fizikā, inženierzinātnēs un datorzinātnēs ļoti bieži rodas "saknes atrašanas" problēma. Sakne ir \(x\) vērtība, kas padara funkciju par nulli, tas ir, vienādojuma risinājums:
\[
f(x)=0
\]
Ne visiem vienādojumiem ir risinājumi, ko var izteikt slēgtas formas formulās, piemēram, kvadrātvienādojumos. Daudziem reālās pasaules gadījumiem, piemēram, sarežģītiem nelineāriem vienādojumiem, mums ir nepieciešamas skaitliskas pieejas. Viena no svarīgākajām pieejām ir iteratīvā metode — procedūra, kas rada virkni aptuvenu risinājumu, kas, izmantojot iterāciju, tuvojas saknei.
Šajā rakstā ir aplūkoti iterācijas metožu pamatjēdzieni, to konverģences nosacījumi un dažas bieži izmantotas iteratīvās metodes sakņu atrašanai.
Sākot no
1. Iterācijas metodes pamatideja
Iterācijas metode darbojas, veicot sākotnējo minējumu \(x_0\) un pēc tam pakāpeniski to uzlabojot, lai iegūtu secību:
\[
x_0, x_1, x_2, \punkti, x_n
\]
ar cerībām:
\[
x_n \to \alpha
\]
kur \(\alpha\) ir vienādojuma \(f(x)=0\) patiesā sakne.
Vispārīgi runājot, iterācijas metode pārveido problēmu \(f(x)=0\) līdzvērtīgā formā:
\[
x = g(x)
\]
Pēc tam tiek veikta iterācija:
\[
x_{n+1} = g(x_n)
\]
Ja šis process konverģē, tad fiksētais punkts \(g(x)\) ir sākotnējā vienādojuma saknes risinājums.
Sākot no
2. Konverģence: Kad iterācija ir veiksmīga?
Ne visas funkcijas \(g(x)\) ģenerē stabilas iterācijas. Lai iterācija \(x_{n+1}=g(x_n)\) konverģētu uz sakni \(\alpha\), bieži tiek izmantoti šādi vispārīgie nosacījumi:
1. \(g(\alpha)=\alpha\) (sakne ir fiksēts punkts)
2. \(|g'(\alpha)| < 1\) (kontraksi lokal)
Intuisi dari \(|g'(\alpha)| < 1\) adalah: di sekitar solusi, fungsi \(g\) “tidak terlalu curam”, sehingga setiap iterasi membawa nilai \(x_n\) semakin dekat, bukan menjauh.
Konvergensi juga dipengaruhi oleh tebakan awal. Dua metode yang sama dapat berhasil atau gagal tergantung pada \(x_0\).
---
3. Metode Bagi Dua (Bisection) sebagai Iterasi Sederhana
Walaupun sering diklasifikasikan terpisah, metode bagi dua dapat dilihat sebagai metode iteratif yang sangat andal. Syaratnya: fungsi \(f(x)\) kontinu pada interval \([a,b]\) dan terjadi perubahan tanda:
\[
f(a)\cdot f(b) < 0
\]
Artinya, ada akar di antara \(a\) dan \(b\). Algoritmanya:
1. Hitung titik tengah \(c=\frac{a+b}{2}\)
2. Tentukan subinterval yang masih mengurung akar (berdasarkan perubahan tanda)
3. Ulangi sampai toleransi tercapai
Kelebihan metode ini: pasti konvergen jika syarat perubahan tanda terpenuhi. Kekurangannya: konvergensinya relatif lambat karena galat berkurang kira-kira setengah tiap iterasi (konvergensi linear).
---
4. Metode Iterasi Titik Tetap (Fixed-Point Iteration)
Inilah bentuk iterasi paling langsung:
\[
x_{n+1} = g(x_n)
\]
Langkahnya:
1. Ubah \(f(x)=0\) menjadi \(x=g(x)\)
2. Pilih tebakan awal \(x_0\)
3. Lakukan iterasi sampai \(|x_{n+1}-x_n|\) atau \(|f(x_n)|\) lebih kecil dari toleransi
Keunggulannya adalah kesederhanaan. Namun, metode ini sangat sensitif terhadap pilihan \(g(x)\). Untuk persamaan yang sama, terdapat banyak cara menulis \(x=g(x)\), tetapi hanya sebagian yang konvergen.
Sebagai contoh, jika ingin mencari akar dari \(f(x)=x^3-2x-5\), kita bisa menulis:
- \(x = \sqrt[3]{2x+5}\) sehingga \(g(x)=\sqrt[3]{2x+5}\)
Lalu lakukan iterasi \(x_{n+1}=\sqrt[3]{2x_n+5}\). Keberhasilan iterasi bergantung pada apakah \(|g'(x)|<1\) di sekitar akar.
---
5. Metode Newton-Raphson: Iterasi Cepat Berbasis Turunan
Metode Newton-Raphson adalah salah satu metode paling populer karena konvergensinya biasanya sangat cepat. Rumus iterasinya:
\[
x_{n+1} = x_n - \frac{f(x_n)}{f'(x_n)}
\]
Interpretasinya: pada \(x_n\), kita buat garis singgung fungsi \(f(x)\). Perpotongan garis singgung dengan sumbu-\(x\) dipakai sebagai perkiraan berikutnya.
Keunggulan:
- Konvergensi kuadratik (sangat cepat) jika sudah cukup dekat dengan akar dan \(f'(\alpha)\neq 0\).
Kekurangan:
- Membutuhkan turunan \(f'(x)\).
- Bisa gagal jika tebakan awal buruk, atau jika \(f'(x_n)\) mendekati nol, sehingga langkah iterasi menjadi tidak stabil.
Metode ini banyak digunakan dalam optimasi, pemodelan fisika, hingga komputasi teknik karena efisiensinya ketika kondisi mendukung.
---
6. Metode Secant: Alternatif Newton Tanpa Turunan
Jika turunan sulit dihitung, metode secant menawarkan kompromi. Ide utamanya adalah mendekati turunan dengan selisih hingga:
\[
f'(x_n)\approx \frac{f(x_n)-f(x_{n-1})}{x_n-x_{n-1}}
\]
Sehingga rumus iterasinya:
\[
x_{n+1}=x_n - f(x_n)\,\frac{x_n-x_{n-1}}{f(x_n)-f(x_{n-1})}
\]
Metode ini memerlukan dua tebakan awal: \(x_0\) dan \(x_1\). Kecepatan konvergensinya umumnya lebih baik daripada bagi dua dan fixed-point sederhana, meskipun biasanya sedikit lebih lambat daripada Newton. Namun, karena tidak perlu turunan, secant sering lebih praktis.
---
7. Kriteria Berhenti (Stopping Criteria)
Dalam komputasi numerik, iterasi harus dihentikan ketika sudah cukup akurat atau jika dicurigai tidak konvergen. Kriteria umum:
1. Galat antar iterasi kecil :
\[
|x_{n+1}-x_n|<\varepsilon
\]
2. Nilai fungsi mendekati nol :
\[
|f(x_n)|<\varepsilon
\]
3. Batas iterasi maksimum untuk mencegah loop tak berujung:
\[
n \le n_{\max}
\]
Pemilihan toleransi \(\varepsilon\) bergantung pada kebutuhan: simulasi teknik mungkin memerlukan toleransi ketat, sementara perhitungan kasar cukup longgar.
---
8. Perbandingan Singkat Metode Iterasi
Secara ringkas:
- Bisection : paling stabil, pasti konvergen (dengan syarat perubahan tanda), tapi lambat.
- Fixed-point : sangat sederhana, tetapi konvergensi tidak selalu terjamin.
- Newton-Raphson : sangat cepat, tetapi butuh turunan dan sensitif terhadap tebakan awal.
- Secant : tidak perlu turunan, cukup cepat, tetapi bisa kurang stabil dibanding bisection.
Dalam praktik, pemilihan metode bergantung pada sifat fungsi, ketersediaan turunan, kebutuhan kecepatan, dan kestabilan.
---
Kesimpulan
Metode iterasi adalah tulang punggung pencarian akar secara numerik untuk persamaan nonlinier. Dengan membangun urutan perkiraan yang diperbarui secara berulang, kita dapat mendekati solusi ketika metode analitik tidak tersedia. Pemahaman tentang konvergensi, pemilihan tebakan awal, dan kriteria berhenti sangat penting agar iterasi menghasilkan akar yang benar dan efisien.
Pada aplikasi nyata, sering kali digunakan strategi gabungan: mulai dengan metode stabil seperti bisection untuk “mengunci” interval akar, lalu beralih ke Newton atau secant untuk mempercepat konvergensi. Dengan demikian, kita memperoleh keseimbangan antara keandalan dan kecepatan—dua aspek yang sangat berharga dalam komputasi numerik.
---
Jika Anda ingin, saya bisa menambahkan contoh perhitungan langkah demi langkah (numerik) untuk salah satu metode di atas agar artikel lebih konkret.
