Kuliah 6 Analisis Numerik Lanjut

Sistem Persamaan Linear

Persamaan Linear

Persamaan linear adalah sebuah persamaan yang tidak ada perkalian antar variabel bebas, seperti 3x + 2y = 1.

Sistem Persamaan Linear

Sebuah sistem yang berisi persamaan-persamaan linear, bentuk umumnya adalah

a₁₁ x₁ + a₁₂ x₂ +...+a_1n x_n = b₁

a₂₁ x₁ + a₂₂ x₂ +...+a_2n x_n = b₂

_...

a_m1 x₁ + a_m2 x₂ +...+a_mn x_n = b_m

jika ditulis dalam bentuk matriks akan menghasilkan

Ax = b

Atau

Pangkat(A) = pangkat(A|b) = n  à  SPL mempunyai solusi tunggal.

Pangkat(A) ≠ pangkat(A|b) à SPL tidak mempunyai solusi.

Bila suatu SPL mempunyai solusi tunggal, maka terdapat banyak cara untuk mencari penyelesaian SPL tersebut, di antaranya adalah : x = A^-1b. 

Metode Substitusi Maju

Metode ini dimulai dengan mengubah terlebih dahulu atriks A menjadi matriks segitiga bawah

kemudian menyelesaikannya mulai dari x₁

Metode Substitusi Mundur

Metode ini dimulai dengan mengubah terlebih dahulu atriks A menjadi matriks segitiga atas

kemudian menyelesaikannya mulai dari x_n

Eliminasi Gauss Naif

Metode ini adalah kembangan dari metode eliminasi.

Langkah penyelesaian:

1.      Tulis SPL dalam bentuk matriks diperbesar

2.      Ubah matriks tersebut menjadi matriks segitiga atas atau segitiga bawah dengan operasi baris dasar (OBD)

Contoh : selesaikan SPL berikut dengan metode eliminasi Gauss

6x₁ – 2x₂ + 2x₃ + 4x₄ = 16

12x₁ – 8x₂ + 6x₃ + 10x₄ = 26

3x₁ – 13x₂ + 9x₃ + 3x₄ = -19

-6x₁ + 4x₂ + x₃ - 18x₄ = -34

Jawaban

kemudian solusi diperoleh dengan menerapkan substitusi mundur

Vektor Eror

Untuk suatu SPL, didefinisikan vektor eror sebagai berikut

yaitu hampiran dikurangi nilai eksak

Vektor Residu

Untuk suatu SPL, didefinisikan vektor residu sebagai berikut

yaitu sisa yang dihasilkan oleh hampiran

Contoh : misalkan diberikan matriks A, vektor b, hampiran, dan nilai eksak sebagai berikut

tentukan eror dan residu

Jawaban

Eliminasi Gauss dengan Pivoting

Metode ini digunakan untuk mengantisipasi gagalnya proses komputasi Eliminasi Gauss tanpa pivoting. Gagalnya eliminasi Gauss disebabkan oleh adanya elemen diagonal utama matriks A yang bernilai nol.

Pivot Parsial

Langkah :

1.      Tentukan r sehingga

2.      Tukar baris i dengan r, jika i=r tidak ditukar

3.      Buat nol elemen di bawah a_ii

4.      Ulangi langkah 1-3 sampai membentuk matriks segitiga atas

5.      Lakukan subtsitusi mundur

Pivot Parsial Terskala

1.      Definisikan vektor indeks l = [l₁,l₂,...,l_n] = [1,2,...,n]

2.      Definisikan vektor skala s = [s₁,s₂,...,s_n] dengan

3.      Tentukan rasio masing-masing baris menggunakan

4.      Pilih j, yaitu indeks dengan rasio maksimum. Baris j adalah pivot untuk iterasi k (k = 1, 2, …, n-1). Jika banyaknya rasio maksimum lebih dari satu, maka pilih indeks terkecil.

5.      Tukar l_k dengan l_j pada vektor indeks

6.      Tukarkan baris pada matriks sesuai dengan vektor indeks.

7.      Buat nol elemen di bawah a_kk.

8.      Kembali ke langkah 3. Vektor indeks yang digunakan adalah yang terbentuk pada langkah 5.

Contoh : Tentukan solusi dari SPL berikut menggunakan eliminasi Gauss dengan pivot parsial terskala

Jawaban (hanya 1 iterasi)

Pertama definisikan vektor indeks dan vektor skala

l = [1,2,3] = [l₁,l₂,l₃]

s = [2,1,4] = [s₁,s₂,s₃]

Untuk iterasi pertama (k = 1), tentukan j, yaitu indeks dengan rasio maksimum. Jika banyaknya rasio maksimum lebih dari satu, maka dipilih indeks terkecil.

Tukarkan l_k  dengan l_j  pada vektor indeks. Artinya, untuk k=1 dan j=2, tukarkan l₁ dan l₂ , sehingga diperoleh vektor indeks baru l = [2,1,3]

Tukarkan baris pada matriks sesuai dengan vektor indeks sehingga diperoleh

Buat nol elemen-elemen di bawah a₁₁

Iterasi pertama selesai

Struktur Banded

Sebuah matriks berukuran n x n disebut mempunyai struktur banded jika terdapat bilangan bulat k (k ≤ n) sehingga a_ij = 0  ketika | i – j | ≥ k.

Tridiagonal System

Suatu sistem yang direpresentasikan dengan matriks yang memenuhi:

1. Terdapat 3 diagonal, yaitu diagonal utama, superdiagonal, dan subdiagonal.
2. Elemen-elemen a_ij ≠ 0 jika |i–j|≤ 1  dan a_ij=0 jika |i–j|≥ 2.

Contoh :

Strictly Diagonal Domiance

Matriks A = (a_ij)_nxn adalah strictly diagonally dominant jika memenuhi syarat berikut

Pentadiagonal System

Suatu sistem yang direpresentasikan dengan matriks yang memenuhi:

1. Terdapat 5 diagonal, yaitu diagonal utama,  2 super-diagonal, dan 2 subdiagonal.
2. Elemen-elemen a_ij ≠ 0 jika |i–j|≥ 2  dan a_ij = 0 jika|i–j|≥ 3, untuk setiap i, j.

Contoh :