Cheat Sheet Materi 7

TI5001 / Optimisasi & Intelejensia Komputasional · Metode Numerik Multi

1 · Peta metode multi variabel

Model tak berpembatas: x=(x1,…,xn)T, opt f(x).

BasisMetode
evaluasi fungsi (nilai f)cyclic coordinate, Hooke-Jeeves, Rosenbrock
turunanNewton, gradient search

Slide merinci cyclic coordinate dan gradient search; yang lain hanya disebut.

2 · Cyclic coordinate: prinsip

  • n arah pencarian = vektor satuan sumbu: dj = 1 di baris j, 0 lainnya. Untuk n=3: d2=(0,1,0)T
  • y(1)=x(k) (awal iterasi k)
y(j+1) = y(j) + λjdjλj = opt f(y(j)+λjdj), fungsi satu variabel f(λj) → metode single (Materi 6)
  • Satu iterasi k: sumbu j=1..n berurutan; x(k+1)=y(n+1) (n=3 → y(4))
  • Pencarian satu variabel setelah K iterasi = n·K (n=2, 3 iterasi → 6)
  • λj boleh negatif

3 · Cyclic coordinate: langkah

  1. Langkah 0: ε>0, d1…dn, x(1), y(1)=x(1), k=1, j=1
  2. Langkah 1: cari λj; y(j+1)=y(j)+λjdj. j<n → j=j+1, ulang; j=n → langkah 2
  3. Langkah 2: x(k+1)=y(n+1). Jika ‖x(k+1)−x(k)‖ < ε berhenti; else y(1)=x(k+1), j=1, k=k+1, ke langkah 1

‖v‖ = √(v1²+v2²). Tidak memakai gradien.

4 · Contoh cyclic coordinate (slide)

(1,2) optimumx(1)=(0,0)kontur f=(x₁−1)²+(2x₁−x₂)²jalur: langkah sejajar sumbu(0,0)→(0,2;0)→(0,2;0,4)→(0,36;0,4)→(0,36;0,72)→…

min f=(x1−1)²+(2x1−x2)², ε=0,001, x(1)=(0,0)

Sumbu 1: f(λ1)=(λ1−1)²+4λ1² → f′=10λ1−2=0 → λ1=0,2, y(2)=(0,2;0)

kx(k)λ1λ2x(k+1)‖Δ‖f
1(0;0)0,20,4(0,2;0,4)0,44720,64
2(0,2;0,4)0,160,32(0,36;0,72)0,35780,4096
3†(0,36;0,72)0,1280,256(0,488;0,976)0,28620,2621

Lanjut sampai ‖Δ‖<ε; optimum (1,2), f=0. λ negatif: dari (2,0), f(2+λ;0)=(1+λ)²+(4+2λ)² → λ1=−1,8, y=(0,2;0). Dari (0,1;0,5): λ1=0,3. Sumbu 2 di (0,5;1): f=0,25+λ2² → λ2=0. Catatan: slide menulis 0,707 untuk ‖Δ‖ k=1, tabel ringkasan 0,4472 (benar).

5 · Gradient search: prinsip & langkah

  • Gradien ∇f(x) = (∂f/∂xj) ; arah penurunan tercepat = −∇f (min) → steepest descent; maks: +∇f → steepest ascent †
d(k) = −∇f(x(k)) / ‖∇f(x(k))‖(atau −∇f tanpa normalisasi). Satuan: panjang d = 1
x(k+1) = x(k) + λ(k)d(k) , λ>0λ(k) = opt f(x(k)+λd(k)) (satu variabel)
  1. 0: ε>0, x(1), k=1
  2. 1: jika ‖∇f(x(k))‖ < ε berhenti; else hitung d(k)
  3. 2: cari λ
  4. 3: x(k+1), k=k+1, ke 1

6 · Contoh gradient search (slide)

f=(x1−1)²+(2x1−x2)²; ∇f=(10x1−4x2−2, −4x1+2x2)T; ε=0,001; x(1)=(0,0)

kx(k)∇f‖∇f‖dλx(k+1)
1(0;0)(−2;0)2(1;0)0,2(0,2;0)
2(0,2;0)(0;−0,8)0,8(0;1)0,4(0,2;0,4)
3(0,2;0,4)(−1,6;0)1,6(1;0)0,16(0,36;0,4)
4†(0,36;0,4)(0;−0,64)0,64(0;1)0,32(0,36;0,72)

k=3: f(λ)=(λ−0,8)²+4λ² → 10λ−1,6=0 → λ=0,16

Di (1,1): ∇f=(4,−2), ‖∇f‖=√20≈4,472, d=(−0,894; 0,447). Di (1,2): ∇f=(0,0) → berhenti

Jalur sama dengan cyclic coordinate karena tiap ∇f hanya punya satu komponen ≠0 (−∇f searah sumbu); fungsi lain umumnya berbeda. Slide menulis √((−2)²+0²)=1,6 pada x(3): seharusnya √((−1,6)²+0)=1,6.

7 · Cyclic coordinate vs gradient search

Cyclic coordinateGradient search
basisnilai fturunan (∇f)
arahtetap: d1…dn (sumbu)−∇f/‖∇f‖, berubah tiap iterasi
λboleh negatif> 0
pencarian 1 var / iterasin1
titik barux(k+1)=y(n+1)x(k)+λd(k)
henti‖x(k+1)−x(k)‖<ε‖∇f(x(k))‖<ε

Dirangkum dari slide kuliah TI5001. Cocokkan lagi dengan slide untuk detail lengkap.