Metode iterasi invers dengan pergeseran dirancang untuk menghitung nilai eigen spesifik yang mendekati tebakan awal. Algoritma ini menggunakan parameter pergeseran μ sebagai panduan pencarian menuju nilai eigen tertentu.
Bayangkan kamu mencari stasiun radio di tengah banyak frekuensi. Tanpa pergeseran, kamu mendengar semua stasiun sekaligus (tidak jelas). Dengan pergeseran, kamu mengarahkan dial ke frekuensi tertentu untuk mendapatkan satu stasiun yang jernih.
Algoritma dimulai dengan matriks A∈Rn×n dan parameter pergeseran . Vektor awal dinormalisasi menjadi dan iterasi dimulai dari .
Langkah pertama adalah menghitung dekomposisi LU dari matriks A−μ⋅I. Dekomposisi ini dilakukan sekali di awal dan digunakan berulang dalam setiap iterasi untuk efisiensi komputasi.
wk+1:=(A−μ⋅I)−1⋅w^k
Praktiknya, alih-alih menghitung invers matriks langsung, kita menyelesaikan sistem persamaan linear menggunakan dekomposisi LU
Setelah memperoleh wk+1, lakukan normalisasi untuk mencegah pertumbuhan tak terkendali
w^k+1:=wk+1/∥wk+1∥
Estimasi nilai eigen menggunakan rasio komponen vektor. Untuk indeks j dengan (w^k)j=0
λ:=(wk+1)j/(w^k)j
Iterasi berlanjut hingga memenuhi kriteria konvergensi ∣1/λ−1/λlama∣≤toleransi. Parameter toleransi adalah nilai ambang batas yang menentukan seberapa akurat hasil yang diinginkan, misalnya 10−6 untuk akurasi enam digit desimal.
Bayangkan seperti mengukur tinggi badan dengan penggaris. Toleransi menentukan seberapa presisi pengukuran yang kamu terima (apakah cukup akurat sampai centimeter, atau perlu sampai milimeter). Semakin kecil nilai toleransi, semakin akurat hasilnya, tetapi memerlukan lebih banyak iterasi.
Hasil akhir memberikan nilai eigen 1/λ+μ dan vektor eigen w^k.
Metode iterasi vektor von Mises menemukan nilai eigen dengan magnitudo terbesar (nilai eigen dominan). Algoritma ini menggunakan proses iterasi sederhana dengan perkalian matriks berulang.
Algoritma dimulai dengan vektor awal w0∈Rn yang dinormalisasi w^0:=w0/∥w0∥ dan iterasi dimulai dari k:=0.
Estimasi nilai eigen menggunakan rasio komponen untuk indeks j dengan (w^k)j=0
λ:=(wk+1)j/(w^k)j
Iterasi berlanjut hingga ∣λ−λlama∣≤toleransi. Nilai toleransi menentukan tingkat ketelitian yang dibutuhkan, biasanya berkisar antara 10−8 hingga 10−12 untuk perhitungan presisi tinggi. Hasil akhir adalah nilai eigen dominan λ dan vektor eigen w^k.
Metode ini berhasil jika ∣λ1∣>∣λ2∣ dan vektor awal w0 memiliki komponen tidak nol dalam arah vektor eigen dominan. Dengan asumsi tersebut, iterasi konvergen ke nilai eigen dominan λ1 dan vektor eigen terkait.
Untuk matriks yang dapat diinvertibel, terdapat hubungan penting antara nilai eigen matriks dan inversnya. Jika A⋅vn=λn⋅vn, maka berlaku
A−1⋅vn=λn1⋅vn
Artinya, nilai eigen terkecil dari A menjadi nilai eigen terbesar dari A−1. Dengan menerapkan iterasi vektor pada A−1, kita memperoleh nilai eigen terkecil dari matriks asli.
Praktiknya, setiap iterasi menyelesaikan sistem linear A⋅wk+1=w^k menggunakan dekomposisi LU, menghindari perhitungan invers eksplisit yang tidak efisien.
