dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 1 dari 19
PENYELESAIAN SISTEM PERSAMAAN ALJABAR LINIER
Materi Kuliah:
Sistem Persamaan Aljabar Linier; Overview Aljabar Matriks;
Metode Eliminasi Gauss; Metode Gauss-Jordan;
Metode Iteratif (Gauss-Seidel & Jacobi); Metode Thomas
# SISTEM PERSAMAAN ALJABAR LINIER #
Secara umum, sistem persamaan aljabar linier dapat dinyatakan dalam bentuk:
a11 x1 + a12 x2 + ....... + a1n xn = b1
a21 x1 + a22 x2 + ....... + a2n xn = b2
.....................................................
an1 x1 + an2 x2 + ....... + ann xn = bn
atau:
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
=
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
n n nn n n
n
n
b
b
b
x
x
x
a a a
a a a
a a a
2
1
2
1
1 2
21 22 2
11 12 1
atau: A x = b
(A dan b diketahui, x harus ditentukan)
Penyelesaian sistem persamaan linier merupakan salah satu masalah sentral dalam ilmu rekayasa
modern, termasuk ilmu rekayasa kimia. Sistem persamaan linier Ax = b memiliki hanya satu
penyelesaian x, jika A berupa matriks wajar (det (A) ≠ 0)
# OVERVIEW ALJABAR MATRIKS #
Vektor adalah himpunan beraturan dari (lambang) bilangan-bilangan yang disebut elemen-elemen
(dari vektor tersebut); banyaknya elemen di dalam himpunan dinamakan order atau dimensi dari
vektor dimaksud. Dari analogi geometrik, sebuah vektor berorder n mendefinisikan (posisi) suatu
titik di dalam ruang berdimensi n.
Vektor kolom:
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
=
n x
x
x
x
2
1
atau: [ ] i x = x
Vektor baris: [ ] n
xT x1 x2 x = , dengan: T ≡ transposisi
Matriks: ≡ himpunan beraturan dari vektor-vektor yang berdimensi sama
≡ jajaran segi-empat dari elemen-elemen (lambang/bilangan)
Matriks A:
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
=
m m mn
n
n
a a a
a a a
a a a
A
1 2
21 22 2
11 12 1
berorder/berdimensi m x n
[ ] ij A = a
aij ≡ elemen matriks A yang posisinya di baris ke-i dan kolom ke-j
aii (i = 1 s.d. n) ≡ elemen diagonal utama
Jika m = n, matriks A disebut matriks bujur sangkar berorder n
Jika m = 1 (hanya satu baris) vektor baris
Jika n = 1 (hanya satu kolom) vektor kolom
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 2 dari 19
Matriks transposisi AT dari matriks A:
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
=
n n mn
m
m
T
a a a
a a a
a a a
A
1 2
12 22 2
11 21 1
berorder/berdimensi n x m
Matriks simetrik: matriks bujur sangkar yang bersifat: aij = aji
Jika A matriks simetrik, maka: AT = A
Contoh:
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
= = −
1 2 1
1 3 2
2 1 1
A AT
Matriks antisimetrik: matriks bujur sangkar yang bersifat: aij = -aji
Jika A matriks antisimetrik, maka: AT = -A
Contoh:
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
−
−
−
=
1 2 0
1 0 2
0 1 1
A
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
−
−
−
= − =
1 2 0
1 0 2
0 1 0
AT A
Matriks kerancang (sparse matrix):
matriks (berdimensi besar) yang sebagian besar
elemennya adalah nol.
Matriks segitiga atas (upper triangular matrix,
U): matriks bujur sangkar yang semua elemen
subdiagonalnya nol.
Contoh:
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
=
0 0 0 5
0 0 6 9
0 3 5 8
1 2 4 7
U
Jika A adalah matriks segitiga atas, maka:
aij = 0 jika i > j
Matriks segitiga bawah (lower triangular matrix,
L): matriks bujur sangkar yang semua elemen
superdiagonalnya nol.
Contoh:
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
=
8 9 2 5
6 1 7 0
3 4 0 0
1 0 0 0
L
Jika B adalah matriks segitiga atas, maka:
bij = 0 jika i < j
Matriks pita (band matrix):
matriks kerancang dengan elemen-elemen tak
nol yang membentuk pita sepanjang diagonal.
Contoh:
⎥ ⎥ ⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢ ⎢ ⎢
⎣
⎡
0 0 0 2 3
0 0 1 2 3
0 4 3 1 2
3 2 1 1 0
4 2 1 0 0
Lebar pita: atas = 2, bawah = 1, total = 4
Matriks tridiagonal: matriks pita dengan
lebar pita atas = 1, bawah = 1, total = 1
Jika A matriks tridiagonal, maka:
1
0 1
− ≤
− >
⎩ ⎨ ⎧ = i j
i j
a
a
ij
ij
Matriks diagonal (D):
matriks bujur sangkar yang bersifat:
i j
i j
d
i
ij =
≠
⎩ ⎨ ⎧
=
λ
0
Contoh:
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
=
3
2
1
0 0
0 0
0 0
λ
λ
λ
D
Matriks identitas (I):
matriks diagonal dengan λi = 1 untuk
semua nilai i
Sifat-Sifat Aljabar Matriks dan Vektor
1. Penjumlahan dan pengurangan matriks
Dua matriks A dan B bisa saling ditambahkan atau dikurangkan jika mempunyai dimensi yang
sama (misal, matriks A dan B berdimensi m x n)
2. Perkalian matriks
Dua matriks A dan B hanya bisa dikalikan menghasilkan matriks produk C = AB, jika:
banyaknya kolom matriks A sama dengan banyaknya baris matriks B.
Jika A berdimensi m x n dan B berdimensi n x p, maka:
C = AB berdimensi m x p dan
j s d p
i s d m
c a b
n
k
ij ik kj 1 . .
1 . .
1 =
=
=Σ=
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 3 dari 19
Pada umumnya:
• Keterdefinisian AB tidak berarti keterdefinisian BA (banyaknya kolom B harus sama
dengan banyaknya baris A)
• AB ≠ BA (jika sama: kebetulan)
Contoh: Jika
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
= −
3 1
2 1
1 0
A dan ⎥⎦
⎤
⎢⎣
⎡
− −
=
1 2 1
2 3 2
B , maka:
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
= ⎥⎦
⎤
⎢⎣
⎡
− −
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
= −
5 11 5
5 4 5
2 3 2
1 2 1
2 3 2
3 1
2 1
1 0
AB
3. Hukum asosiatif dan distributif terhadap perkalian dengan skalar:
αA = Aα αAB = AαB = ABα α(A+B) = αA + αB
4. Hukum asosiatif untuk penjumlahan/pengurangan dan perkalian:
A ± (B ± C) = (A ± B) ± C A(BC) = (AB)C
5. Hukum komutatif untuk penjumlahan/pengurangan: A + B = B + A
6. Hukum distributif pekalian terhadap penjumlahan:
A(B + C) = AB + AC (A + B)C = AC + BC
7. (A + B)T = AT + BT (AB)T = BTAT (ABC)T = CT(AB)T = CTBTAT
Matriks Kebalikan (Invers)
Jika A dan B adalah dua matriks bujur sangkar yang sama-sama berdimensi n dan BA = AB = I,
maka B disebut matriks kebalikan (invers) dari matriks A dan diberi simbol A-1, sehingga:
A-1A = AA-1 = I
Contoh:
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
=
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− −
−
−
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− 0 0 1
0 1 0
1 0 0
5/ 4 3/ 4 1/ 4
1/ 2 1/ 2 1/ 2
1/8 1/8 1/8
2 2 1
4 1 0
2 1 1
A A-1 I
Matriks bujur sangkar yang tidak memiliki kebalikan disebut matriks tak wajar (singular
matrix). Matriks ini mempunyai determinan = 0.
Aturan-aturan untuk matriks kebalikan:
• (αA)-1 = α-1A-1 dengan: α ≡ skalar
• (A-1)T = (AT)-1
• (AB)-1 = B-1A-1
• Jika D adalah matriks diagonal dengan elemen-elemen diagonal dii (i = 1 s.d. n), maka D-1
adalah juga matriks diagonal dengan elemen-elemen diagonalnya adalah
ii d
1
Determinan Matriks
Determinan adalah fungsi bernilai skalar dari semua elemen sebuah matriks bujur sangkar.
Contoh: ⎥⎦
⎤
⎢⎣
⎡
=
c d
a b
A det (A) = ad – bc
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
=
p q r
d e f
a b c
A det (A) = aer +bfp + cdq – pec – qfa – rdb
Determinan dari suatu skalar (alias matriks berorder 1 x 1) adalah nilai skalar itu sendiri.
Determinan matriks bujur sangkar berorder > 1 dapat dievaluasi dengan menggunakan formula
pengembangan dari Laplace, yakni:
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 4 dari 19
Jika A adalah matriks bujur sangkar berorder n, maka:
Σ=
= − +
n
k
ik
i k
ik A a M
1
det ( ) ( 1) dengan: i = 1 atau 2 ... atau n
atau:
Σ=
= − +
n
k
kj
k j
kj A a M
1
det ( ) ( 1) dengan: j = 1 atau 2 ... atau n
di mana: rs
(−1)r + s M adalah kofaktor dari ars dan Mrs adalah determinan minor ars
Minor
Jika Q adalah matriks bujur sangkar berorder (n-1) yang diperoleh dengan cara menghapus
baris ke-i dan kolom ke-j dari matriks bujur sangkar A (yang berorder n), maka matriks Q
dinamakan minor aij dari matriks A.
Contoh:
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
=
31 32 33
21 22 23
11 12 13
a a a
a a a
a a a
A
Maka: Minor a11 = ⎥⎦
⎤
⎢⎣
⎡
32 33
22 23
a a
a a
Minor a12 = ⎥⎦
⎤
⎢⎣
⎡
31 33
21 23
a a
a a
Minor a13 = ⎥⎦
⎤
⎢⎣
⎡
31 32
21 22
a a
a a
Minor a21 = ⎥⎦
⎤
⎢⎣
⎡
32 33
12 13
a a
a a
Minor a22 = ⎥⎦
⎤
⎢⎣
⎡
31 33
11 13
a a
a a
Minor a23 = ⎥⎦
⎤
⎢⎣
⎡
31 32
11 12
a a
a a
Minor a31 = ⎥⎦
⎤
⎢⎣
⎡
22 23
12 13
a a
a a
Minor a32 = ⎥⎦
⎤
⎢⎣
⎡
21 23
11 13
a a
a a
Minor a33 = ⎥⎦
⎤
⎢⎣
⎡
21 22
11 12
a a
a a
Tiap matriks bujur sangkar berorder n mempunyai n2 buah minor.
Contoh:
ad bc cb da ad cb bc da
c d
a b
= − = − + = − = − +
p q
d e
c
p r
d f
b
q r
e f
a
p q r
d e f
a b c
= (1) + (−1) + (1) = aer – aqf – bdr + bpf + cdq – cpe
= aer + bpf + cdq – pec – qfa – rdb
7 3
4 5
3
7 7
4 6
2
3 7
5 6
1
7 3 7
4 5 6
1 2 3
= − + = 35 – 18 – 56 + 84 + 36 – 105 = -24
Sifat-sifat determinan:
1. det (I) = 1
2. det (AT) = det (A)
3. Jika A dan B sama-sama berorder n, maka: det (AB) = det (A).det (B)
4. det (AB) = det (BA)
5. Jika A merupakan matriks tak wajar, maka: det (A) = 0
6. Jika elemen-elemen dari suatu baris (atau kolom) dari A semuanya bernilai nol, maka:
det A) = 0
7. Jika A merupakan matriks segitiga atas/bawah atau matriks diagonal, maka:
det (A) = a11.a22. .... . ann
8. Jika dua baris (atau kolom) dari A sama, maka: det (A) = 0
9. Jika matriks B diperoleh dengan mempertukarkan dua baris (atau kolom) pada matriks A,
maka: det (B) = - det (A)
10. Pengurangan kelipatan suatu baris atau kolom dari baris (kolom) tidak mengubah nilai
determinan.
11. det (αA) = αn . det (A)
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 5 dari 19
Aturan Cramer untuk Penyelesaian Sistem Persamaan Linier Ax = b
Jika diketahui matriks bujur sangkar A yang berorder n dan vektor b yang juga berorder n, maka
penyelesaian sistem persamaan linier Ax = b adalah:
det (A)
x i
i
Δ
= dengan: i = 1 s.d. n
di mana Δi adalah determinan dari matriks yang diperoleh dengan mengganti kolom ke-i dari
matriks A dengan vektor b.
Contoh: ⎥⎦
⎤
⎢⎣
⎡
= ⎥⎦
⎤
⎢⎣
⎡
⎥⎦
⎤
⎢⎣
⎡
0
1
3 4
2 3
2
1
x
x
det(A) = 8 – 9 = -1
A x b
4 0 4
0 4
1 3
1 Δ = = − = dan 0 3 3
3 0
2 1
2 Δ = = − = −
Maka: 4
1
4
det ( )
1
1 = −
−
=
Δ
=
A
x dan 3
1
3
det ( )
2
2 =
−
−
=
Δ
=
A
x
sehingga: x = [ -4 3 ]T
# METODE ELIMINASI GAUSS #
Eliminasi Gauss merupakan salah satu metode penyelesaian sistem persamaan aljabar linier. Strategi
penyelesaian sistem n buah persamaan linier Ax = b, di mana A adalah matriks bujur sangkar biasa
dilakukan dengan:
1. Mereduksi Ax = b menjadi Ux = c, di mana U adalah matriks segitiga atas berorder n, dan
2. Menyelesaikan sistem Ux = c tersebut di atas dengan metode substitusi balik (backward
substitution), sehingga diperoleh x.
Metode pereduksian tersebut di atas disebut proses eliminasi Gauss, dan akan diilustrasikan dalam
contoh berikut ini:
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
−
=
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
1
1
1
1
1 2 3 4
2 3 4 3
3 4 3 2
4 3 2 1
4
3
2
1
x
x
x
x
A x b
Bentuk matriks perbesarannya:
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
−
1
1
1
1
1 2 3 4
2 3 4 3
3 4 3 2
4 3 2 1
A b
Langkah-langkahnya:
1. Buatlah elemen-elemen di bawah a11 menjadi bernilai nol, dengan cara:
a) Hitung:
11
21
21 a
l = a ,
11
31
31 a
l = a , dan
11
41
41 a
l = a
4
3
21 l =
2
1
4
2
31 l = =
4
1
41 l =
b) Baris ke-2 baru = baris ke-2 lama – l21 (baris ke-1)
Baris ke-3 baru = baris ke-3 lama – l31 (baris ke-1)
Baris ke-4 baru = baris ke-4 lama – l41 (baris ke-1)
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 6 dari 19
sehingga diperoleh:
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
−
5/ 4
3/ 2
1/ 4
1
0 5/ 4 5/ 2 15/ 4
0 3/ 2 3 5/ 2
0 7 / 4 3/ 2 5/ 4
4 3 2 1
2. Buatlah elemen-elemen di bawah a22 (yang baru) menjadi bernilai nol, dengan cara:
a) Hitung:
22
32
32 a
l = a dan
22
42
42 a
l = a
7
6
7 / 4
3/ 2
32 l = =
7
5
7 / 4
5/ 4
42 l = =
b) Baris ke-3 baru = baris ke-3 lama – l32 (baris ke-2)
Baris ke-4 baru = baris ke-4 lama – l42 (baris ke-2)
sehingga diperoleh:
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
−
10/ 7
12/ 7
1/ 4
1
0 0 10/ 7 20/ 7
0 0 12/ 7 10/ 7
0 7 / 4 3/ 2 5/ 4
4 3 2 1
3. Buatlah elemen-elemen di bawah a33 (yang baru) menjadi bernilai nol, dengan cara:
a) Hitung:
33
43
43 a
l = a
6
5
12/ 7
10/ 7
43 l = =
b) Baris ke-4 baru = baris ke-4 lama – l43 (baris ke-3)
sehingga diperoleh:
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
0
12/ 7
1/ 4
1
0 0 0 5/ 3
0 0 12/ 7 10/ 7
0 7 / 4 3/ 2 5/ 4
4 3 2 1
U c
Proses eliminasi selesai setelah (n-1) tahap/langkah, di mana n adalah order dari matriks koefisien.
Penyelesaian persamaan:
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
0 0 0 5/ 3
0 0 12/ 7 10/ 7
0 7 / 4 3/ 2 5/ 4
4 3 2 1
⎥ ⎥ ⎥ ⎥ ⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
4
3
2
1
x
x
x
x
=
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
0
12 / 7
1/ 4
1
U x = c
selanjutnya dapat dilakukan dengan cara substitusi balik:
0
5/3
0
4 x = =
1
12/ 7
(( 12/ 7) (10/ 7)(0))
3 = −
− −
x =
1
7 / 4
((1/ 4) (3/ 2)( 1) (5/ 4)(0))
2 =
− − −
x =
0
4
(1 (3)(1) (2)( 1) (1)(0))
1 =
− − − −
x =
sehingga diperoleh: x = [0 1 -1 0]T
Verifikasi hasil: = b
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
−
=
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
1
1
1
1
0
1
1
0
1 2 3 4
2 3 4 3
3 4 3 2
4 3 2 1
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 7 dari 19
Persoalan yang Muncul dan Pemecahannya
Persoalan-persoalan yang dapat dijumpai dalam penyelesaian sistem persamaan linier dengan metode
eliminasi Gauss meliputi:
1. Division by Zero
Pembagian dengan bilangan nol sangat mungkin terjadi dalam penyelesaian sistem persamaan
linier dengan eliminasi Gauss.
Contoh: Gunakan eliminasi Gauss untuk menyelesaikan kumpulan persamaan:
2 x2 + 3 x3 = 8
4 x1 + 6 x2 + 7 x3 = -3
2 x1 + x2 + 6 x3 = 5
Persoalan yang sama dapat dijumpai jika elemen diagonal matriks koefisiennya mendekati nol
(sangat kecil).
2. Round-Off Errors
Error yang disebabkan oleh pembulatan bilangan yang terjadi selama proses eliminasi
3. Ill-Conditioned Systems
Penyelesaian sebuah sistem persamaan ditentukan oleh kondisi sistem tersebut. Pada wellconditioned
systems, perubahan kecil dalam satu atau lebih koefisiennya mengakibatkan
perubahan yang kecil dalam penyelesaiannya. Dan sebaliknya, pada ill-conditioned systems,
perubahan kecil dalam satu atau lebih koefisiennya mengakibatkan perubahan yang besar
dalam penyelesaiannya. Ill-conditioned systems juga memungkinkan terjadinya jawaban
dengan rentang yang lebar, yang merupakan penyelesaian dari sistem persamaan tersebut.
4. Singular Systems
Kasus ini terjadi jika determinan matriks koefisien sistem persamaan bernilai nol (atau
merupakan matriks tak wajar), sehingga tidak mempunyai penyelesaian.
Persoalan tersebut di atas (1-3) dapat dipecahkan melalui beberapa cara, seperti (1) Penggunaan
significant figures yang lebih banyak, (2) Pivoting, maupun (3) Scaling.
Contoh-contoh:
# Ill-Conditioned Systems
Selesaikan sistem persamaan: x1 + 2 x2 = 10
1,1 x1 + 2 x2 = 10,4
Ulangi lagi, tetapi jika koefisien x1 dalam persamaan kedua diganti menjadi 1,05.
Penyelesaian:
Dengan menggunakan aturan Cramer, diperoleh:
4
1(2) 2(1,1)
2(10) 2(10,4)
1 =
−
−
x = dan 3
1(2) 2(1,1)
1(10,4) 1,1(10)
2 =
−
−
x =
Dengan cara yang sama, untuk koefisien x1 dalam persamaan kedua diganti menjadi 1,05:
8
1(2) 2(1,05)
2(10) 2(10,4)
1 =
−
−
x = dan 1
1(2) 2(1,05)
1(10,4) 1,05(10)
2 =
−
−
x =
[Perhatikanlah bahwa perubahan yang kecil pada salah satu koefisien mengakibatkan
penyelesaian (x1 dan x2) yang jauh berbeda.]
# Pengaruh Scaling terhadap Determinan
Lakukan evaluasi nilai determinan (det) pada sistem persamaan berikut ini:
(a) 3 x1 + 2 x2 = 18 (b) x1 + 2 x2 = 10
- x1 + 2 x2 = 2 1,1 x1 + 2 x2 = 10,4
(c) Ulangilah untuk sistem (b), tetapi jika kedua persamaan tersebut dikalikan faktor 10.
Penyelesaian:
(a) det = 3(2) – 2(-1) = 8
Sistem ini tergolong well-conditioned system.
(b) det = 1(2) – 2(1,1) = -0,2
Sistem ini tergolong ill-conditioned system (karena mempunyai determinan yang
nilainya mendekati nol)
(Perhatikan bahwa tahap reduksi atau eliminasi
yang pertama akan gagal, karena a11 = 0)
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 8 dari 19
(c) Sistem (b) dikalikan 10 menjadi: 10 x1 + 20 x2 = 100
11 x1 + 20 x2 = 104
Perkalian sebuah persamaan dengan sebuah konstanta (scaling) tidak mempengaruhi
penyelesaiannya. Dengan kata lain, sistem persamaan tetap bersifat ill-conditioned.
Namun demikian, nilai determinannya jauh berubah, yakni:
det = 10(20) – 20(11) = -20
# Partial Pivoting
Gunakan eliminasi Gauss untuk menyelesaikan: 0,0003 x1 + 3,0000 x2 = 2,0001 (i)
1,0000 x1 + 1,0000 x2 = 1,0000 (ii)
(Catatan: Penyelesaian eksak sistem persamaan ini adalah: x1 = 1/3 dan x2 = 2/3)
Penyelesaian:
Perhatikan bahwa elemen pivot yang pertama, a11 = 0,0003, mendekati nol...!!
Cara pertama:
Jika persamaan (i) dikalikan 1/(0,0003), maka diperoleh: x1 + 10000 x2 = 6667
yang dapat digunakan untuk mengeliminasi x1 dari persamaan (ii):
-9999 x2 = -6666 x2 = 2/3
x2 di atas selanjutnya dapat disubstitusikan balik ke persamaan (i) untuk mengevaluasi x1:
0,0003
2,0001 3(2/ 3)
1
−
x =
Namun demikian, hasil yang diperoleh sangat dipengaruhi oleh banyaknya significant figures
yang digunakan, seperti ditabelkan berikut ini:
Significant
figures x2 x1 Nilai mutlak percent
relative error untuk x1
3 0,667 -3,33 1099
4 0,6667 0,0000 100
5 0,66667 0,30000 10
6 0,666667 0,330000 1
7 0,6666667 0,3330000 0,1
Cara lain: Pivoting
Jika kedua persamaan saling dipertukarkan posisinya (supaya elemen pivot pertama
merupakan nilai yang terbesar dalam satu kolomnya), maka:
1,0000 x1 + 1,0000 x2 = 1,0000 (i.a)
0,0003 x1 + 3,0000 x2 = 2,0001 (ii.a)
Melalui eliminasi dan substitusi, diperoleh: x2 = 2/3.
Selanjutnya, x2 ini disubstitusikan balik sehingga diperoleh x1 sebesar:
1
1 (2/3)
1
−
x =
Pada beberapa nilai significant figures yang digunakan, diperoleh hasil perhitungan sbb.:
Significant
figures x2 x1 Nilai mutlak percent
relative error untuk x1
3 0,667 0,333 0,1
4 0,6667 0,3333 0,01
5 0,66667 0,33333 0,001
6 0,666667 0,333333 0,0001
7 0,6666667 0,3333333 0,00001
Berdasarkan hasil di atas, terlihat bahwa significant-figures-yang-digunakan mempunyai
pengaruh yang jauh kurang sensitif terhadap perhitungan (bandingkan dengan cara
sebelumnya...!!).
Kesimpulan: Pada kasus ini, strategi pivoting memberikan hasil yang jauh lebih baik.
# METODE GAUSS-JORDAN #
Metode eliminasi Gauss merupakan “separuh jalan” dari metode Gauss-Jordan, yakni metode yang
dilakukan untuk mengubah matriks bujur sangkar A menjadi sebuah matriks diagonal D. Proses
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 9 dari 19
eliminasi yang dilakukan terhadap matriks segitiga atas U agar menjadi matriks diagonal D dilakukan
melalui langkah-langkah yang serupa. Ilustrasi dalam contoh sebelumnya akan diproses lebih lanjut
menggunakan metode Gauss-Jordan. (lihat kembali contoh sebelumnya)
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
−
=
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
1
1
1
1
1 2 3 4
2 3 4 3
3 4 3 2
4 3 2 1
4
3
2
1
x
x
x
x
A x b
Matriks perbesaran:
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
−
1
1
1
1
1 2 3 4
2 3 4 3
3 4 3 2
4 3 2 1
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
0
12/ 7
1/ 4
1
0 0 0 5/ 3
0 0 12/ 7 10/ 7
0 7 / 4 3/ 2 5/ 4
4 3 2 1
A b U c
7
6
5/3
10/ 7
34 g = =
4
3
5/3
5/ 4
24 g = =
5
3
5/3
1
14 g = =
8
7
12/ 7
3/ 2
23 g = =
6
7
12/ 7
2
13 g = =
7
12
7 / 4
3
12 g = =
D s
Dengan demikian: 0
4
0
1 x = =
1
7 / 4
7 / 4
2 x = =
1
12/ 7
12/ 7
3 = −
−
x =
0
5/3
0
4 x = =
Karena penyelesaian Ax = b sudah dapat dicapai secara efisien dengan metode eliminasi Gauss,
penggunaan metode Gauss-Jordan bukanlah dalam penyelesaian sistem persamaan linier, melainkan
dalam perhitungan manual penentuan matriks kebalikan (atau invers) A-1 dari matriks bujur
sangkar A.
A . A-1 = I
(dapat ditentukan/dicari dengan metode Gauss Jordan)
Eliminasi Gauss
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
0
12 / 7
1/ 4
1
0 0 0 5/ 3
0 0 12 / 7 0
0 7 / 4 3/ 2 0
4 3 2 0
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
0
12/ 7
7 / 4
3
0 0 0 5/ 3
0 0 12/ 7 0
0 7 / 4 0 0
4 3 0 0
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
0
12/ 7
7 / 4
0
0 0 0 5/ 3
0 0 12/ 7 0
0 7 / 4 0 0
4 0 0 0
x = [0 1 -1 0]T
(sama dengan hasil yang diperoleh sebelumnya,
dengan metode eliminasi Gauss)
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 10 dari 19
Contoh: Tentukan invers dari matriks A =
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− 2 2 1
4 1 0
2 1 1
, dengan metode Gauss-Jordan!
Penyelesaian:
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− 2 2 1 0 0 1
4 1 0 0 1 0
2 1 1 1 0 0
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− − −
0 3 2 1 0 1
0 1 2 2 1 0
2 1 1 1 0 0
A I
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− −
− − −
0 0 4 5 3 1
0 1 2 2 1 0
2 1 1 1 0 0
(keadaan pada tahap akhir eliminasi Gauss)
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− −
− − −
0 0 4 5 3 1
0 1 2 2 1 0
2 1 1 1 0 0
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− −
− − −
−
0 0 4 5 3 1
0 1 0 1/ 2 1/ 2 1/ 2
2 1 0 1/ 4 3/ 4 1/ 4
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− −
− − −
−
0 0 4 5 3 1
0 1 0 1/ 2 1/ 2 1/ 2
2 0 0 1/ 4 1/ 4 1/ 4
Selanjutnya, baris pertama dibagi dengan 2, baris kedua dibagi dengan (-1), dan baris ketiga
dibagi dengan (-4), sehingga diperoleh:
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− −
−
−
0 0 1 5/ 4 3/ 4 1/ 4
0 1 0 1/ 2 1/ 2 1/ 2
1 0 0 1/8 1/8 1/8
I A-1
Verifikasi hasil:
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− 2 2 1
4 1 0
2 1 1
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− −
−
−
5/ 4 3/ 4 1/ 4
1/ 2 1/ 2 1/ 2
1/8 1/8 1/8
=
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− −
−
−
5/ 4 3/ 4 1/ 4
1/ 2 1/ 2 1/ 2
1/8 1/8 1/8
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− 2 2 1
4 1 0
2 1 1
=
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
0 0 1
0 1 0
1 0 0
A A-1 = A-1 A = I
Jadi, matriks invers A adalah: A-1 =
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎣
⎡
− −
−
−
5/ 4 3/ 4 1/ 4
1/ 2 1/ 2 1/ 2
1/8 1/8 1/8
# METODE ITERATIF GAUSS-SEIDEL #
Metode Gauss-Seidel merupakan metode penyelesaian sistem persamaan aljabar linier yang dilakukan
secara iteratif. Berawal dari nilai-nilai tebakan awal untuk x, nilai awal ini selanjutnya diproses secara
berulang-ulang (atau iteratif), hingga diperoleh nilai x yang cukup dekat dengan nilai sesungguhnya
(atau memenuhi batas toleransi yang telah ditetapkan). Untuk sistem persamaan linier yang melibatkan
matriks A dan vektor b yang berorder n, diperlukan nilai tebakan awal x sejumlah n-1 buah.
Tinjaulah sistem persamaan linier (berorder 3) yang dinyatakan dalam bentuk matriks sbb.:
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
31 32 33
21 22 23
11 12 13
a a a
a a a
a a a
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
3
2
1
x
x
x
=
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
3
2
1
b
b
b
A x = b
Tahap pertama
Tahap pertama
Tahap kedua
Tahap kedua
eliminasi Gauss
eliminasi Gauss
Gauss-Jordan
Gauss-Jordan
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 11 dari 19
Dari baris 1, diperoleh: 3
11
13
2
11
12
11
1
1 x
a
x a
a
a
a
x = b − −
Dengan cara yang sama, x2 dan x3 dapat ditentukan dari baris 2 dan baris 3:
3
22
23
1
22
21
22
2
2 x
a
x a
a
a
a
x = b − − dan 2
33
32
1
33
31
33
3
3 x
a
x a
a
a
a
x = b − −
Pendekatan tersebut di atas digunakan dalam metode iteratif Gauss-Seidel. Jika dipilih 2 nilai tebakan
awal untuk x2 dan x3, maka x1 dapat dihitung. Berdasarkan x1 (baru) yang dihitung x3 (lama), maka x2
(baru) dapat dihitung. Berdasarkan x1 (baru) dan x2 (baru), maka x3 (baru) dapat dihitung. Demikian
seterusnya, sampai diperoleh harga-harga x1, x2, dan x3 yang konvergen. Prosedur iteratif ini
digambarkan dalam diagram alir berikut ini. Prosedur yang analog dapat diterapkan terhadap sistem
persamaan linier berorder n.
Persentase error perhitungan dapat dinyatakan sebagai:
100%
,
, , x
x
x x
Error
i iterasi sekarang
i iterasi sekarang i iterasi sebelumnya
i
−
=
Contoh:
Gunakan metode Gauss-Seidel untuk menyelesaikan sistem persamaan:
3 x1 – 0,1 x2 – 0,2 x3 = 7,85
0,1 x1 + 7 x2 – 0,3 x3 = -19,3
0,3 x1 – 0,2 x2 + 10 x3 = 71,4
Catatan: Nilai sebenarnya dari sistem persamaan ini adalah: x1 = 3; x2 = -2,5; dan x3 = 7
Penyelesaian:
Dari setiap persamaan, dapat diperoleh x1, x2, dan x3 sebesar:
3
7,85 0,1 0,2 2 3
1
x x x + +
=
7
19,3 0,1 0,3 1 3
2
x x x − − +
=
10
71,4 0,3 0,2 1 2
3
x x x − +
=
MULAI
Tebakan awal: x2,lama dan x3,lama
Toleransi: tol
baru lama lama x
a
x a
a
a
a
x b 3,
11
13
2,
11
12
11
1
1, = − −
baru baru lama x
a
x a
a
a
a
x b 3,
22
23
1,
22
21
22
2
2, = − −
baru baru baru x
a
x a
a
a
a
x b 2,
33
32
1,
33
31
33
3
3, = − −
Error ≤ tol ?
x1 = x1,baru
x2 = x2,baru
x3 = x3,baru
SELESAI
x2,lama = x2,baru
x3,lama = x3,baru
Ya
Tidak
ke iterasi
berikutnya
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 12 dari 19
Iterasi pertama:
Dengan menggunakan nilai tebakan awal untuk x2 dan x3 sebesar: x2 = x3 = 0, maka
diperoleh: 2,616667
3
7,85 0,1 (0) 0,2 (0)
1 =
+ +
x =
Gunakan x1 = 2,616667 ini dan x3 = 0 untuk menghitung x2:
2,794524
7
19,3 0,1 (2,616667) 0,3 (0)
2 = −
− − +
x =
Gunakan x1 = 2,616667 dan x2 = -2,794524 untuk menghitung x3:
7,005610
10
71,4 0,3 (2,616667) 0,2 ( 2,794524)
3 =
− + −
x =
Iterasi kedua:
Dengan cara yang sama, maka diperoleh:
2,990557
3
7,85 0,1 ( 2,794524) 0,2 (7,005610)
1 =
+ − +
x = = 0,31%; = 12,5% t a ε ε
2,499625
7
19,3 0,1 (2,990557) 0,3 (7,005610)
2 = −
− − +
x = = 0,015%; = 11,8% t a ε ε
7,000291
10
71,4 0,3 (2,990557) 0,2 ( 2,499625)
3 =
− + −
x = = 0,0042%; = 0,076% t a ε ε
Komentar: Perhatikan bahwa metode ini bersifat konvergen (menuju nilai-nilai x sebenarnya).
Langkah-langkah iterasi berikutnya dapat dilakukan untuk meningkatkan ketelitian dan
kedekatan dengan nilai sebenarnya.
Kriteria Konvergensi untuk Metode Gauss-Seidel
Metode Gauss-Seidel mirip dengan metode iterasi satu titik (dalam penentuan akar persamaan tak
linier), yang dapat memunculkan 2 kemungkinan persoalan, yakni: (1) kadang-kadang tidak
konvergen, atau (2) jika konvergen, maka laju konvergensinya lambat.
Untuk sistem persamaan linier Ax = b yang melibatkan n buah persamaan, konvergensi metode ini
akan selalu tercapai, jika: Σ
≠ =
>
n
i j
j
ii ij a a
1
(diagonally dominant)
Metode relaksasi merupakan modifikasi dari metode Gauss-Seidel yang dirancang untuk
meningkatkan konvergensi. Setelah masing-masing xi dihitung dengan persamaan-persamaan di
atas, nilai xi ini selanjutnya dimodifikasi dengan memberikan faktor pembobot terhadap hasil
iterasi sekarang dan iterasi sebelumnya; atau:
i baru i baru i lama x x x , , , = λ + (1− λ )
di mana: λ ≡ faktor pembobot (bernilai antara 0 – 2; biasanya ditentukan secara empirik)
Jika λ = 1, maka metode ini tidak dimodifikasi/direlaksasi
Jika 0 ≤ λ ≤ 1, maka modifikasinya tergolong dalam underrelaxation (biasanya digunakan untuk
membuat sistem yang tidak konvergen menjadi konvergen)
Jika 1 ≤ λ ≤ 2, maka modifikasinya tergolong dalam overrelaxation (biasanya digunakan untuk
mempercepat konvergensi sistem yang sudah konvergen). Pendekatan ini biasa disebut juga
successive or simultaneous overrelaxation (SOR).
# METODE ITERATIF JACOBI #
Metode Jacobi merupakan metode penyelesaian sistem persamaan aljabar linier secara iteratif yang
mirip dengan metode Gauss-Seidel. Berawal dari nilai-nilai tebakan awal untuk x, nilai awal ini
selanjutnya diproses secara berulang-ulang (atau iteratif), hingga diperoleh nilai x yang cukup dekat
dengan nilai sesungguhnya (atau memenuhi batas toleransi yang telah ditetapkan). Untuk sistem
persamaan linier yang melibatkan matriks A dan vektor b yang berorder n, diperlukan nilai tebakan
awal x sejumlah n buah.
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 13 dari 19
Tinjaulah sistem persamaan linier (berorder 3) yang dinyatakan dalam bentuk matriks sbb.:
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
31 32 33
21 22 23
11 12 13
a a a
a a a
a a a
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
3
2
1
x
x
x
=
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
3
2
1
b
b
b
A x = b
Seperti halnya metode Gauss-Seidel, metode Jacobi juga akan selalu menghasilkan penyelesaian yang
konvergen jika matriks koefisien A bersifat diagonally dominant. Jika A merupakan matriks bujur
sangkar yang bersifat diagonally dominant, maka dari baris 1, diperoleh: 3
11
13
2
11
12
11
1
1 x
a
x a
a
a
a
x = b − −
Dengan cara yang sama, x2 dan x3 dapat ditentukan dari baris 2 dan baris 3:
3
22
23
1
22
21
22
2
2 x
a
x a
a
a
a
x = b − − dan 2
33
32
1
33
31
33
3
3 x
a
x a
a
a
a
x = b − −
Jika dipilih 3 nilai tebakan awal untuk x1, x2, dan x3, maka x1 (baru), x2 (baru), dan x3 (baru) dapat
dihitung secara bersamaan. Demikian seterusnya, sampai diperoleh harga-harga x1, x2, dan x3 yang
konvergen. Prosedur iteratif ini digambarkan dalam diagram alir berikut ini. Prosedur yang analog
dapat diterapkan terhadap sistem persamaan linier berorder n.
Persentase error perhitungan dapat dinyatakan sebagai:
100%
,
, , x
x
x x
Error
i iterasi sekarang
i iterasi sekarang i iterasi sebelumnya
i
−
=
Contoh:
Ulangi soal sebelumnya dengan metode iteratif Jacobi!
Penyelesaian:
Jika digunakan nilai-nilai tebakan awal x, yaitu: x1 = x2 = x3 = 0
dan dengan cara yang mirip dengan sebelumnya, maka diperoleh hasil-hasil di bawah ini.
MULAI
Tebakan awal: x1,lama, x2,lama,
dan x3,lama ; Toleransi: tol
baru lama lama x
a
x a
a
a
a
x b 3,
11
13
2,
11
12
11
1
1, = − −
baru lama lama x
a
x a
a
a
a
x b 3,
22
23
1,
22
21
22
2
2, = − −
baru lama lama x
a
x a
a
a
a
x b 2,
33
32
1,
33
31
33
3
3, = − −
Error ≤ tol ?
x1 = x1,baru
x2 = x2,baru
x3 = x3,baru
SELESAI
x1,lama = x1,baru
x2,lama = x2,baru
x3,lama = x3,baru
Ya
Tidak
ke iterasi
berikutnya
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 14 dari 19
Tabulasi hasil perhitungan (dengan metode Jacobi):
Bandingkan hasil ini dengan hasil yang diperoleh dengan metode Gauss-Seidel sbb:
# METODE THOMAS UNTUK MENYELESAIKAN SISTEM PERSAMAAN
LINIER DENGAN MATRIKS KOEFISIEN TRIDIAGONAL #
Metode eliminasi Gauss dapat diterapkan untuk sistem persamaan linier yang matriks koefisiennya
berupa matriks tridiagonal. (Ingatlah kembali, apa yang dimaksud dengan matriks tridiagonal...!).
Metode terapan ini biasa disebut metode Thomas, yang alur penyelesaiannya akan digambarkan
dalam uraian berikut ini.
Tinjaulah sistem persamaan linier:
a11 x1 + a12 x2 = b1
a21 x1 + a22 x2 + a23 x3 = b2
a32 x2 + a33 x3 + a34 x4 = b3
.........................................................................
an,n-1 xn-1 + ann xn = bn
Atau, dapat dituliskan menjadi:
d1 x1 + e1 x2 = b1
c2 x1 + d2 x2 + e2 x3 = b2
c3 x2 + d3 x3 + e3 x4 = b3
.........................................................................
cn xn-1 + dn xn = bn
Dalam bentuk perkalian matriks dan vektor, persamaan di atas dapat dituliskan sebagai:
⎥ ⎥ ⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢ ⎢ ⎢
⎣
⎡
=
⎥ ⎥ ⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢ ⎢ ⎢
⎣
⎡
⎥ ⎥ ⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢ ⎢ ⎢
⎣
⎡
n n n n b
b
b
b
x
x
x
x
c d
c d e
c d e
d e
... ...
0 0 0
... ... ... ... ...
0 0
0 0
0 0 0
3
2
1
3
2
1
3 3 3
2 2 2
1 1
... (*)
Matriks tridiagonal
Dengan menerapkan langkah-langkah yang telah dipelajari dalam metode eliminasi Gauss (pada
materi sebelumnya), penyelesaian dari sistem persamaan linier (*) adalah sbb.:
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 15 dari 19
n n x = γ
dan:
i
i i
i i
x e x
β
= γ − +1 (i = n-1, n-2, ..., 1)
dengan γi dan βi yang ditentukan melalui serangkaian perhitungan berurutan sbb.:
1 1 β = d
1
1
1 β
γ = b
1
1
−
= − −
i
i i
i i
d c e
β
β (i = 2, 3, ..., n)
i
i i i
i
b c
β
γ
γ −1 −
= (i = 2, 3, ..., n)
Contoh:
Selesaikan sistem persamaan berikut:
2 x1 – 3 x2 = -4
x1 + 2 x2 – x3 = 2
4 x2 – x3 + x4 = 9
2 x3 – x4 = 2
Penyelesaian:
Sistem persamaan linier di atas memuat matriks koefisien tridiagonal. Mula-mula besaranbesaran
ci, di, ei, dan bi akan didefinisikan lebih dahulu, untuk menyesuaikan dengan notasi
yang diuraikan dalam metode Thomas tersebut di atas.
i ci di ei bi
1 0 2 -3 -4
2 1 2 -1 2
3 4 -1 1 9
4 2 -1 0 2
Dengan metode Thomas, maka diperoleh hasil-hasil sbb.:
β1
= d1 = 2
γ1
=
1
1β
b =
2
− 4 = -2
β2
=
1
2 1
2 β
d − c e =
2
2 1. ( 3)
−
− = 3,5
γ2
=
2
2 2 1
β
b − c γ = 3,5
2 −1. (−2) = 1,14
β3
=
2
3 2
3 β
d − c e = 3,5
1 4 . ( 1)
−
− − = 0,143
γ3
=
3
3 3 2
β
b − c γ = 0,143
9 − 4 .1,14 = 31
β4
=
3
4 3
4 β
d − c e = 0,143
−1− 2 .1 = -15
γ4
=
4
4 4 3
β
b − c γ =
15
2 2 . 31
−
− = 4
Dengan demikian:
x4 = γ
4 = 4
x3 =
3
3 4
3 β
γ − e x = 0,143
31− 1. 4 = 3
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 16 dari 19
x2 =
2
2 3
2 β
γ − e x =
3,5
1,14 ( 1) (3)
−
− = 2
x1 =
1
1 2
1 β
γ − e x =
2
2 ( 3) (2)
−
− − = 1
Jadi, penyelesaian sistem persamaan linier tersebut di atas adalah:
x1 = 1; x2 = 2; x3 = 3; dan x4 = 4
# CONTOH APLIKASI #
Neraca Massa Rangkaian Proses
Umpan berupa zat A murni dengan laju 100 kmol/jam.
Kendala: 1. 80% dari A dan 40% dari B di dalam alur 2 didaur ulang (recycle).
2. Perbandingan mol A terhadap mol B di dalam alur 1 adalah 5:1
Neraca massa (dalam kmol/jam):
Di sekitar pencampur:
NA1 = NA3 + 100 atau: NA1 − NA3 = 100 ... (1)
NB1 = NB3 atau: NB1 − NB3 = 0 ... (2)
(NA1 menyatakan laju alir molar A di dalam alur 1, dst.)
Di sekitar reaktor:
NA2 = NA1 − r atau: − NA1 + NA2 + r = 0 ... (3)
NB2 = NB1 + r atau: − NB1 + NB2 − r = 0 ... (4)
(r menyatakan laju reaksi)
Di sekitar pemisah:
NA3 + NA4 = NA2 atau: − NA2 + NA3 + NA4 = 0 ... (5)
NB3 + NB4 = NB2 atau: − NB2 + NB3 + NB4 = 0 ... (6)
Berdasarkan kendala 1:
NA3 = 0,8 NA2 atau: − 0,8 NA2 + NA3 = 0 ... (7)
NB3 = 0,4 NB2 atau: − 0,4 NB2 + NB3 = 0 ... (8)
Berdasarkan kendala 2:
NA1 = 5 NB1 atau: NA1 − 5 NB1 = 0 ... (9)
Berdasarkan penjabaran neraca massa di atas, dihasilkan 9 buah persamaan linier dengan 9 variabel
yang tak diketahui (yakni NA1, NB1, NA2, NB2, NA3, NB3, NA4, NB4, dan r)
Dengan demikian, terbentuk sistem persamaan linier yang dapat diselesaikan secara simultan!
Pemisah
Produk
Umpan
A → B
Pencampur Reaktor
4
3
1 2
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 17 dari 19
(Sebagai latihan, susunlah ke-9 persamaan di atas dalam bentuk perkalian matriks dan
vektor, dan selanjutnya silakan Anda selesaikan sendiri sistem persamaan linier tersebut.
Untuk mengecek jawaban Anda, hasil perhitungan dari persoalan neraca massa ini adalah:
NA1 = 227,3 kmol/jam; NB1 = 45,45 kmol/jam; NA2 = 159,1 kmol/jam; NB2 = 113,6 kmol/jam;
NA3 = 127,3 kmol/jam; NB3 = 45,45 kmol/jam; NA4 = 31,82 kmol/jam; NB4 = 68,18 kmol/jam;
dan r = 68,18 kmol/jam)
# LATIHAN SOAL #
1. Tuliskan sekumpulan persamaan berikut ini dalam bentuk matriks:
2 3 30 = 2 x + 6 x
2 1 20 = 3 x + 8 x
1 3 10 = x + x
Tuliskan tranpose matriks koefisiennya!
2. Gunakan metode grafis untuk menyelesaikan sistem persamaan berikut:
2 6 18 1 2 x − x = −
8 40 1 2 − x + x =
Cek jawaban Anda dengan mensubstitusikannya kembali ke dalam persamaan-persamaan di atas.
3. Untuk sekumpulan persamaan berikut:
2 5 1 2 3 x + x =
2 2 1 1 2 3 x + x + x =
3 2 1 2 x + x =
(a) Hitunglah nilai determinannya
(b) Gunakan cara Cramer untuk menentukan harga x1, x2, dan x3.
(c) Substitusikan jawaban Anda (pada butir (b)) ke dalam persamaan-persamaan semula, untuk
mengecek kebenarannya.
4. Gunakan eliminasi Gauss untuk menyelesaikan sistem persamaan berikut:
4 2 1 2 3 x + x − x = −
5 2 4 1 2 3 x + x + x =
6 6 1 2 3 x + x + x =
Lakukan partial column pivoting dan cek kebenaran jawaban Anda dengan mensubstitusikan
kembali ke dalam persamaan-persamaan semula.
5. Selesaikan sistem persamaan:
3 1 2 3 x + x − x = −
6 2 2 2 1 2 3 x + x + x =
3 4 1 1 2 3 − x + x + x =
dengan metode:
(a) eliminasi Gauss naive
(b) eliminasi Gauss dengan partial pivoting
(c) Gauss-Jordan naive
(d) Gauss-Jordan dengan partial pivoting
6. Sistem persamaan berikut ini dirancang untuk menentukan konsentrasi (c, dalam g/m3) dalam
serangkaian reaktor sebagai fungsi jumlah massa yang diumpankan ke dalam masing-masing
reaktor (pada ruas kanan persamaan, dalam g/hari):
17 2 3 500 1 2 3 c − c − c =
5 21 2 200 1 2 3 − c + c − c =
5 5 22 30 1 2 3 − c − c + c =
(a) Gunakan matriks invers untuk menentukan harga konsentrasi c1, c2, dan c3.
(b) Tentukan c1, c2, dan c3 dengan metode Gauss-Seidel, hingga tercapai εs = 5%.
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 18 dari 19
7. Gunakan metode Gauss-Seidel dengan relaksasi (dengan λ = 0,90 dan εs = 5%) untuk menyelesaikan
sistem persamaan:
5 12 80 1 3 − x + x =
4 2 1 2 3 x − x − x = −
6 8 45 1 2x + x =
(Jika perlu, lakukan penyusunan ulang terhadap persamaan untuk mencapai konvergensi).
8. Selesaikan 3 buah persamaan linier berikut ini dengan metode eliminasi Gauss, dan tunjukkan
semua langkah perhitungannya.
3 2 3 1 2 3 x − x + x =
2 4 2 2 1 2 3 x + x − x =
4 2 3 12 1 2 3 x − x − x = −
Ulangilah jika menggunakan metode iteratif Gauss-Seidel (atau Jacobi).
9. Selesaikan sistem persamaan linier berikut ini dengan metode eliminasi Gauss, dan tunjukkan
semua langkah perhitungannya.
2 3 2 3 x − x =
2 3 1 1 2 3 x − x + x =
3 2 5 10 1 2 3 x − x + x =
10. Selesaikan 3 buah persamaan linier berikut ini dengan metode eliminasi Gauss, dan tunjukkan
semua langkah perhitungannya.
3 3 2 1 2 3 x − x + x =
5 3 6 1 2 3 x + x + x =
2 2 1 2 3 x + x − x =
11. Selesaikan sistem persamaan linier berikut ini dengan metode Thomas, dan tunjukkan semua
langkah perhitungannya.
1 1 2 x − x =
2 3 1 1 2 3 x + x − x =
3 4 2 3 2 3 4 x − x + x =
5 3 3 4 x − x =
Ulangilah jika menggunakan metode iteratif Gauss-Seidel (atau Jacobi).
12. Sebuah kolom ekstraksi yang memiliki 4 tahap ideal/teoretik akan digunakan untuk mengekstraksi
asam format yang terlarut dalam metil isobutil keton (MIBK) dengan menggunakan air sebagai
pengekstrak. Operasi akan dilakukan dengan mode berlawanan arah (counter current); fasa organik
diumpankan dari dasar kolom dengan laju 5 m3 MIBK/jam dan konsentrasi asam format 1,0
mol/liter MIBK, sedangkan fasa akuatik dialirkan dari puncak kolom dengan laju alir 2,45 m3
air/jam serta sama sekali tidak mengandung asam format. Operasi berlangsung isotermal dan harga
rasio kesetimbangan konsentrasi asam format di dalam MIBK terhadap konsentrasi asam format di
dalam air adalah K = 0,445. MIBK dan air dapat dianggap sama sekali tidak larut satu sama lain. 1
m3 = 1000 liter. Tentukan konsentrasi asam format pada aliran-aliran yang keluar dan masuk tiap
pelat/tahap ideal dalam kolom ekstraksi tersebut.
13. Sebuah kilang minyak menghasilkan 3 jenis elpiji: A, B, dan C, dengan komposisi seperti disajikan
pada tabel berikut ini (angka-angka yang ditunjukkan menyatakan %-volume komponen):
Komponen A ElBpi ji C
C2 5,0 - -
C3 90,0 10,0 -
i-C4 5,0 85,0 8,0
n-C4 - 5,0 80,0
i-C5 - - 12,0
Jumlah 100,0 100,0 100,0
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 19 dari 19
Dari sebuah perusahaan calon pelanggan, kilang minyak tersebut menerima tawaran kontrak jual
beli 106 m3 elpiji yang disyaratkan harus mengandung: 31,2%-v propana, 53,4%-v isobutana, dan
12,6%-v n-butana. Kadar etana dan i-C5 tidak disyaratkan. Kilang minyak tersebut bermaksud
memenuhi kontrak ini dengan cara mencampurkan elpiji A, elpiji B, dan elpiji C.
Tentukan kebutuhan masing-masing elpiji A, B, dan C (dalam m3) agar persyaratan kontrak dapat
dipenuhi dengan baik!
14. Gambar berikut ini memperlihatkan sistem tiga buah reaktor yang terhubungkan oleh pipa-pipa.
Laju perpindahan setiap komponen/zat melalui setiap pipa merupakan perkalian dari laju alir
volumetriknya (Q, dalam satuan m3/detik) dengan konsentrasinya yang keluar dari masing-masing
reaktor (C, dalam satuan mg/m3). Sistem proses berada dalam kondisi steady, yang berarti bahwa
laju massa masuk ke setiap reaktor sama dengan laju massa keluar.
Susunlah neraca massa pada masing-masing reaktor dan selanjutnya selesaikanlah sistem
persamaan aljabar linier yang terbentuk (untuk menentukan harga-harga C1, C2, dan C3).
15. Tinjaulah proses perpindahan panas secara konduksi melalui dinding sebuah pipa. Suhu dinding
bagian dalam: 200oC, suhu dinding bagian luar: 80oC, dan ketebalan dinding pipa: 0,05 meter. Jarijari
bagian dalam pipa: 0,05 meter (ro) dan jari-jari bagian luar pipa: 0,1 meter. Persamaan
diferensial yang dijabarkan berdasarkan neraca panas dalam dinding pipa dan menggambarkan
distribusi suhu sepanjang dinding pipa dinyatakan sbb.:
0 2
2
= ⎟
⎟⎠
⎞
⎜ ⎜⎝
⎛
+
dr
r d T
dr
dT
Pendekatan beda hingga (finite difference) yang diterapkan untuk menyelesaikan persamaan di atas
menghasilkan serangkaian persamaan aljabar linier berikut:
0
2
2 1
2
1
1 1 = ⎟ ⎟⎠
⎞
⎜ ⎜⎝
⎛
Δ
+ − +
Δ
− ⎟
⎟⎠
⎞
⎜ ⎜⎝
⎛
Δ
+ i + i i − T
r
T r
r
T r
r
r
untuk: i = 2, 3, ..., 10
dengan:
r = r + Δ r (i −1) = 0,05 + 0,005 (i −1) o
dan
T1 = 200oC dan T11 = 80oC
Tentukan distribusi suhu di sepanjang dinding pipa menggunakan metode Thomas! (Jangan lupa
bahwa r tidaklah konstan)
☺☺☺ Selamat Belajar!!! ☺☺☺
400
mg/detik
200
mg/detik
3
2
1
Q13C1
Q12C1 Q23C2
Q21C2
Q33C3
Q33 = 120
Q13 = 40
Q12 = 80
Q23 = 60
Q21 = 20
Keterangan:
PENYELESAIAN SISTEM PERSAMAAN ALJABAR LINIER
Materi Kuliah:
Sistem Persamaan Aljabar Linier; Overview Aljabar Matriks;
Metode Eliminasi Gauss; Metode Gauss-Jordan;
Metode Iteratif (Gauss-Seidel & Jacobi); Metode Thomas
# SISTEM PERSAMAAN ALJABAR LINIER #
Secara umum, sistem persamaan aljabar linier dapat dinyatakan dalam bentuk:
a11 x1 + a12 x2 + ....... + a1n xn = b1
a21 x1 + a22 x2 + ....... + a2n xn = b2
.....................................................
an1 x1 + an2 x2 + ....... + ann xn = bn
atau:
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
=
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
n n nn n n
n
n
b
b
b
x
x
x
a a a
a a a
a a a
2
1
2
1
1 2
21 22 2
11 12 1
atau: A x = b
(A dan b diketahui, x harus ditentukan)
Penyelesaian sistem persamaan linier merupakan salah satu masalah sentral dalam ilmu rekayasa
modern, termasuk ilmu rekayasa kimia. Sistem persamaan linier Ax = b memiliki hanya satu
penyelesaian x, jika A berupa matriks wajar (det (A) ≠ 0)
# OVERVIEW ALJABAR MATRIKS #
Vektor adalah himpunan beraturan dari (lambang) bilangan-bilangan yang disebut elemen-elemen
(dari vektor tersebut); banyaknya elemen di dalam himpunan dinamakan order atau dimensi dari
vektor dimaksud. Dari analogi geometrik, sebuah vektor berorder n mendefinisikan (posisi) suatu
titik di dalam ruang berdimensi n.
Vektor kolom:
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
=
n x
x
x
x
2
1
atau: [ ] i x = x
Vektor baris: [ ] n
xT x1 x2 x = , dengan: T ≡ transposisi
Matriks: ≡ himpunan beraturan dari vektor-vektor yang berdimensi sama
≡ jajaran segi-empat dari elemen-elemen (lambang/bilangan)
Matriks A:
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
=
m m mn
n
n
a a a
a a a
a a a
A
1 2
21 22 2
11 12 1
berorder/berdimensi m x n
[ ] ij A = a
aij ≡ elemen matriks A yang posisinya di baris ke-i dan kolom ke-j
aii (i = 1 s.d. n) ≡ elemen diagonal utama
Jika m = n, matriks A disebut matriks bujur sangkar berorder n
Jika m = 1 (hanya satu baris) vektor baris
Jika n = 1 (hanya satu kolom) vektor kolom
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 2 dari 19
Matriks transposisi AT dari matriks A:
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
=
n n mn
m
m
T
a a a
a a a
a a a
A
1 2
12 22 2
11 21 1
berorder/berdimensi n x m
Matriks simetrik: matriks bujur sangkar yang bersifat: aij = aji
Jika A matriks simetrik, maka: AT = A
Contoh:
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
= = −
1 2 1
1 3 2
2 1 1
A AT
Matriks antisimetrik: matriks bujur sangkar yang bersifat: aij = -aji
Jika A matriks antisimetrik, maka: AT = -A
Contoh:
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
−
−
−
=
1 2 0
1 0 2
0 1 1
A
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
−
−
−
= − =
1 2 0
1 0 2
0 1 0
AT A
Matriks kerancang (sparse matrix):
matriks (berdimensi besar) yang sebagian besar
elemennya adalah nol.
Matriks segitiga atas (upper triangular matrix,
U): matriks bujur sangkar yang semua elemen
subdiagonalnya nol.
Contoh:
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
=
0 0 0 5
0 0 6 9
0 3 5 8
1 2 4 7
U
Jika A adalah matriks segitiga atas, maka:
aij = 0 jika i > j
Matriks segitiga bawah (lower triangular matrix,
L): matriks bujur sangkar yang semua elemen
superdiagonalnya nol.
Contoh:
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
=
8 9 2 5
6 1 7 0
3 4 0 0
1 0 0 0
L
Jika B adalah matriks segitiga atas, maka:
bij = 0 jika i < j
Matriks pita (band matrix):
matriks kerancang dengan elemen-elemen tak
nol yang membentuk pita sepanjang diagonal.
Contoh:
⎥ ⎥ ⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢ ⎢ ⎢
⎣
⎡
0 0 0 2 3
0 0 1 2 3
0 4 3 1 2
3 2 1 1 0
4 2 1 0 0
Lebar pita: atas = 2, bawah = 1, total = 4
Matriks tridiagonal: matriks pita dengan
lebar pita atas = 1, bawah = 1, total = 1
Jika A matriks tridiagonal, maka:
1
0 1
− ≤
− >
⎩ ⎨ ⎧ = i j
i j
a
a
ij
ij
Matriks diagonal (D):
matriks bujur sangkar yang bersifat:
i j
i j
d
i
ij =
≠
⎩ ⎨ ⎧
=
λ
0
Contoh:
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
=
3
2
1
0 0
0 0
0 0
λ
λ
λ
D
Matriks identitas (I):
matriks diagonal dengan λi = 1 untuk
semua nilai i
Sifat-Sifat Aljabar Matriks dan Vektor
1. Penjumlahan dan pengurangan matriks
Dua matriks A dan B bisa saling ditambahkan atau dikurangkan jika mempunyai dimensi yang
sama (misal, matriks A dan B berdimensi m x n)
2. Perkalian matriks
Dua matriks A dan B hanya bisa dikalikan menghasilkan matriks produk C = AB, jika:
banyaknya kolom matriks A sama dengan banyaknya baris matriks B.
Jika A berdimensi m x n dan B berdimensi n x p, maka:
C = AB berdimensi m x p dan
j s d p
i s d m
c a b
n
k
ij ik kj 1 . .
1 . .
1 =
=
=Σ=
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 3 dari 19
Pada umumnya:
• Keterdefinisian AB tidak berarti keterdefinisian BA (banyaknya kolom B harus sama
dengan banyaknya baris A)
• AB ≠ BA (jika sama: kebetulan)
Contoh: Jika
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
= −
3 1
2 1
1 0
A dan ⎥⎦
⎤
⎢⎣
⎡
− −
=
1 2 1
2 3 2
B , maka:
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
= ⎥⎦
⎤
⎢⎣
⎡
− −
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
= −
5 11 5
5 4 5
2 3 2
1 2 1
2 3 2
3 1
2 1
1 0
AB
3. Hukum asosiatif dan distributif terhadap perkalian dengan skalar:
αA = Aα αAB = AαB = ABα α(A+B) = αA + αB
4. Hukum asosiatif untuk penjumlahan/pengurangan dan perkalian:
A ± (B ± C) = (A ± B) ± C A(BC) = (AB)C
5. Hukum komutatif untuk penjumlahan/pengurangan: A + B = B + A
6. Hukum distributif pekalian terhadap penjumlahan:
A(B + C) = AB + AC (A + B)C = AC + BC
7. (A + B)T = AT + BT (AB)T = BTAT (ABC)T = CT(AB)T = CTBTAT
Matriks Kebalikan (Invers)
Jika A dan B adalah dua matriks bujur sangkar yang sama-sama berdimensi n dan BA = AB = I,
maka B disebut matriks kebalikan (invers) dari matriks A dan diberi simbol A-1, sehingga:
A-1A = AA-1 = I
Contoh:
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
=
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− −
−
−
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− 0 0 1
0 1 0
1 0 0
5/ 4 3/ 4 1/ 4
1/ 2 1/ 2 1/ 2
1/8 1/8 1/8
2 2 1
4 1 0
2 1 1
A A-1 I
Matriks bujur sangkar yang tidak memiliki kebalikan disebut matriks tak wajar (singular
matrix). Matriks ini mempunyai determinan = 0.
Aturan-aturan untuk matriks kebalikan:
• (αA)-1 = α-1A-1 dengan: α ≡ skalar
• (A-1)T = (AT)-1
• (AB)-1 = B-1A-1
• Jika D adalah matriks diagonal dengan elemen-elemen diagonal dii (i = 1 s.d. n), maka D-1
adalah juga matriks diagonal dengan elemen-elemen diagonalnya adalah
ii d
1
Determinan Matriks
Determinan adalah fungsi bernilai skalar dari semua elemen sebuah matriks bujur sangkar.
Contoh: ⎥⎦
⎤
⎢⎣
⎡
=
c d
a b
A det (A) = ad – bc
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
=
p q r
d e f
a b c
A det (A) = aer +bfp + cdq – pec – qfa – rdb
Determinan dari suatu skalar (alias matriks berorder 1 x 1) adalah nilai skalar itu sendiri.
Determinan matriks bujur sangkar berorder > 1 dapat dievaluasi dengan menggunakan formula
pengembangan dari Laplace, yakni:
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 4 dari 19
Jika A adalah matriks bujur sangkar berorder n, maka:
Σ=
= − +
n
k
ik
i k
ik A a M
1
det ( ) ( 1) dengan: i = 1 atau 2 ... atau n
atau:
Σ=
= − +
n
k
kj
k j
kj A a M
1
det ( ) ( 1) dengan: j = 1 atau 2 ... atau n
di mana: rs
(−1)r + s M adalah kofaktor dari ars dan Mrs adalah determinan minor ars
Minor
Jika Q adalah matriks bujur sangkar berorder (n-1) yang diperoleh dengan cara menghapus
baris ke-i dan kolom ke-j dari matriks bujur sangkar A (yang berorder n), maka matriks Q
dinamakan minor aij dari matriks A.
Contoh:
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
=
31 32 33
21 22 23
11 12 13
a a a
a a a
a a a
A
Maka: Minor a11 = ⎥⎦
⎤
⎢⎣
⎡
32 33
22 23
a a
a a
Minor a12 = ⎥⎦
⎤
⎢⎣
⎡
31 33
21 23
a a
a a
Minor a13 = ⎥⎦
⎤
⎢⎣
⎡
31 32
21 22
a a
a a
Minor a21 = ⎥⎦
⎤
⎢⎣
⎡
32 33
12 13
a a
a a
Minor a22 = ⎥⎦
⎤
⎢⎣
⎡
31 33
11 13
a a
a a
Minor a23 = ⎥⎦
⎤
⎢⎣
⎡
31 32
11 12
a a
a a
Minor a31 = ⎥⎦
⎤
⎢⎣
⎡
22 23
12 13
a a
a a
Minor a32 = ⎥⎦
⎤
⎢⎣
⎡
21 23
11 13
a a
a a
Minor a33 = ⎥⎦
⎤
⎢⎣
⎡
21 22
11 12
a a
a a
Tiap matriks bujur sangkar berorder n mempunyai n2 buah minor.
Contoh:
ad bc cb da ad cb bc da
c d
a b
= − = − + = − = − +
p q
d e
c
p r
d f
b
q r
e f
a
p q r
d e f
a b c
= (1) + (−1) + (1) = aer – aqf – bdr + bpf + cdq – cpe
= aer + bpf + cdq – pec – qfa – rdb
7 3
4 5
3
7 7
4 6
2
3 7
5 6
1
7 3 7
4 5 6
1 2 3
= − + = 35 – 18 – 56 + 84 + 36 – 105 = -24
Sifat-sifat determinan:
1. det (I) = 1
2. det (AT) = det (A)
3. Jika A dan B sama-sama berorder n, maka: det (AB) = det (A).det (B)
4. det (AB) = det (BA)
5. Jika A merupakan matriks tak wajar, maka: det (A) = 0
6. Jika elemen-elemen dari suatu baris (atau kolom) dari A semuanya bernilai nol, maka:
det A) = 0
7. Jika A merupakan matriks segitiga atas/bawah atau matriks diagonal, maka:
det (A) = a11.a22. .... . ann
8. Jika dua baris (atau kolom) dari A sama, maka: det (A) = 0
9. Jika matriks B diperoleh dengan mempertukarkan dua baris (atau kolom) pada matriks A,
maka: det (B) = - det (A)
10. Pengurangan kelipatan suatu baris atau kolom dari baris (kolom) tidak mengubah nilai
determinan.
11. det (αA) = αn . det (A)
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 5 dari 19
Aturan Cramer untuk Penyelesaian Sistem Persamaan Linier Ax = b
Jika diketahui matriks bujur sangkar A yang berorder n dan vektor b yang juga berorder n, maka
penyelesaian sistem persamaan linier Ax = b adalah:
det (A)
x i
i
Δ
= dengan: i = 1 s.d. n
di mana Δi adalah determinan dari matriks yang diperoleh dengan mengganti kolom ke-i dari
matriks A dengan vektor b.
Contoh: ⎥⎦
⎤
⎢⎣
⎡
= ⎥⎦
⎤
⎢⎣
⎡
⎥⎦
⎤
⎢⎣
⎡
0
1
3 4
2 3
2
1
x
x
det(A) = 8 – 9 = -1
A x b
4 0 4
0 4
1 3
1 Δ = = − = dan 0 3 3
3 0
2 1
2 Δ = = − = −
Maka: 4
1
4
det ( )
1
1 = −
−
=
Δ
=
A
x dan 3
1
3
det ( )
2
2 =
−
−
=
Δ
=
A
x
sehingga: x = [ -4 3 ]T
# METODE ELIMINASI GAUSS #
Eliminasi Gauss merupakan salah satu metode penyelesaian sistem persamaan aljabar linier. Strategi
penyelesaian sistem n buah persamaan linier Ax = b, di mana A adalah matriks bujur sangkar biasa
dilakukan dengan:
1. Mereduksi Ax = b menjadi Ux = c, di mana U adalah matriks segitiga atas berorder n, dan
2. Menyelesaikan sistem Ux = c tersebut di atas dengan metode substitusi balik (backward
substitution), sehingga diperoleh x.
Metode pereduksian tersebut di atas disebut proses eliminasi Gauss, dan akan diilustrasikan dalam
contoh berikut ini:
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
−
=
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
1
1
1
1
1 2 3 4
2 3 4 3
3 4 3 2
4 3 2 1
4
3
2
1
x
x
x
x
A x b
Bentuk matriks perbesarannya:
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
−
1
1
1
1
1 2 3 4
2 3 4 3
3 4 3 2
4 3 2 1
A b
Langkah-langkahnya:
1. Buatlah elemen-elemen di bawah a11 menjadi bernilai nol, dengan cara:
a) Hitung:
11
21
21 a
l = a ,
11
31
31 a
l = a , dan
11
41
41 a
l = a
4
3
21 l =
2
1
4
2
31 l = =
4
1
41 l =
b) Baris ke-2 baru = baris ke-2 lama – l21 (baris ke-1)
Baris ke-3 baru = baris ke-3 lama – l31 (baris ke-1)
Baris ke-4 baru = baris ke-4 lama – l41 (baris ke-1)
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 6 dari 19
sehingga diperoleh:
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
−
5/ 4
3/ 2
1/ 4
1
0 5/ 4 5/ 2 15/ 4
0 3/ 2 3 5/ 2
0 7 / 4 3/ 2 5/ 4
4 3 2 1
2. Buatlah elemen-elemen di bawah a22 (yang baru) menjadi bernilai nol, dengan cara:
a) Hitung:
22
32
32 a
l = a dan
22
42
42 a
l = a
7
6
7 / 4
3/ 2
32 l = =
7
5
7 / 4
5/ 4
42 l = =
b) Baris ke-3 baru = baris ke-3 lama – l32 (baris ke-2)
Baris ke-4 baru = baris ke-4 lama – l42 (baris ke-2)
sehingga diperoleh:
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
−
10/ 7
12/ 7
1/ 4
1
0 0 10/ 7 20/ 7
0 0 12/ 7 10/ 7
0 7 / 4 3/ 2 5/ 4
4 3 2 1
3. Buatlah elemen-elemen di bawah a33 (yang baru) menjadi bernilai nol, dengan cara:
a) Hitung:
33
43
43 a
l = a
6
5
12/ 7
10/ 7
43 l = =
b) Baris ke-4 baru = baris ke-4 lama – l43 (baris ke-3)
sehingga diperoleh:
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
0
12/ 7
1/ 4
1
0 0 0 5/ 3
0 0 12/ 7 10/ 7
0 7 / 4 3/ 2 5/ 4
4 3 2 1
U c
Proses eliminasi selesai setelah (n-1) tahap/langkah, di mana n adalah order dari matriks koefisien.
Penyelesaian persamaan:
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
0 0 0 5/ 3
0 0 12/ 7 10/ 7
0 7 / 4 3/ 2 5/ 4
4 3 2 1
⎥ ⎥ ⎥ ⎥ ⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
4
3
2
1
x
x
x
x
=
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
0
12 / 7
1/ 4
1
U x = c
selanjutnya dapat dilakukan dengan cara substitusi balik:
0
5/3
0
4 x = =
1
12/ 7
(( 12/ 7) (10/ 7)(0))
3 = −
− −
x =
1
7 / 4
((1/ 4) (3/ 2)( 1) (5/ 4)(0))
2 =
− − −
x =
0
4
(1 (3)(1) (2)( 1) (1)(0))
1 =
− − − −
x =
sehingga diperoleh: x = [0 1 -1 0]T
Verifikasi hasil: = b
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
−
=
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
1
1
1
1
0
1
1
0
1 2 3 4
2 3 4 3
3 4 3 2
4 3 2 1
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 7 dari 19
Persoalan yang Muncul dan Pemecahannya
Persoalan-persoalan yang dapat dijumpai dalam penyelesaian sistem persamaan linier dengan metode
eliminasi Gauss meliputi:
1. Division by Zero
Pembagian dengan bilangan nol sangat mungkin terjadi dalam penyelesaian sistem persamaan
linier dengan eliminasi Gauss.
Contoh: Gunakan eliminasi Gauss untuk menyelesaikan kumpulan persamaan:
2 x2 + 3 x3 = 8
4 x1 + 6 x2 + 7 x3 = -3
2 x1 + x2 + 6 x3 = 5
Persoalan yang sama dapat dijumpai jika elemen diagonal matriks koefisiennya mendekati nol
(sangat kecil).
2. Round-Off Errors
Error yang disebabkan oleh pembulatan bilangan yang terjadi selama proses eliminasi
3. Ill-Conditioned Systems
Penyelesaian sebuah sistem persamaan ditentukan oleh kondisi sistem tersebut. Pada wellconditioned
systems, perubahan kecil dalam satu atau lebih koefisiennya mengakibatkan
perubahan yang kecil dalam penyelesaiannya. Dan sebaliknya, pada ill-conditioned systems,
perubahan kecil dalam satu atau lebih koefisiennya mengakibatkan perubahan yang besar
dalam penyelesaiannya. Ill-conditioned systems juga memungkinkan terjadinya jawaban
dengan rentang yang lebar, yang merupakan penyelesaian dari sistem persamaan tersebut.
4. Singular Systems
Kasus ini terjadi jika determinan matriks koefisien sistem persamaan bernilai nol (atau
merupakan matriks tak wajar), sehingga tidak mempunyai penyelesaian.
Persoalan tersebut di atas (1-3) dapat dipecahkan melalui beberapa cara, seperti (1) Penggunaan
significant figures yang lebih banyak, (2) Pivoting, maupun (3) Scaling.
Contoh-contoh:
# Ill-Conditioned Systems
Selesaikan sistem persamaan: x1 + 2 x2 = 10
1,1 x1 + 2 x2 = 10,4
Ulangi lagi, tetapi jika koefisien x1 dalam persamaan kedua diganti menjadi 1,05.
Penyelesaian:
Dengan menggunakan aturan Cramer, diperoleh:
4
1(2) 2(1,1)
2(10) 2(10,4)
1 =
−
−
x = dan 3
1(2) 2(1,1)
1(10,4) 1,1(10)
2 =
−
−
x =
Dengan cara yang sama, untuk koefisien x1 dalam persamaan kedua diganti menjadi 1,05:
8
1(2) 2(1,05)
2(10) 2(10,4)
1 =
−
−
x = dan 1
1(2) 2(1,05)
1(10,4) 1,05(10)
2 =
−
−
x =
[Perhatikanlah bahwa perubahan yang kecil pada salah satu koefisien mengakibatkan
penyelesaian (x1 dan x2) yang jauh berbeda.]
# Pengaruh Scaling terhadap Determinan
Lakukan evaluasi nilai determinan (det) pada sistem persamaan berikut ini:
(a) 3 x1 + 2 x2 = 18 (b) x1 + 2 x2 = 10
- x1 + 2 x2 = 2 1,1 x1 + 2 x2 = 10,4
(c) Ulangilah untuk sistem (b), tetapi jika kedua persamaan tersebut dikalikan faktor 10.
Penyelesaian:
(a) det = 3(2) – 2(-1) = 8
Sistem ini tergolong well-conditioned system.
(b) det = 1(2) – 2(1,1) = -0,2
Sistem ini tergolong ill-conditioned system (karena mempunyai determinan yang
nilainya mendekati nol)
(Perhatikan bahwa tahap reduksi atau eliminasi
yang pertama akan gagal, karena a11 = 0)
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 8 dari 19
(c) Sistem (b) dikalikan 10 menjadi: 10 x1 + 20 x2 = 100
11 x1 + 20 x2 = 104
Perkalian sebuah persamaan dengan sebuah konstanta (scaling) tidak mempengaruhi
penyelesaiannya. Dengan kata lain, sistem persamaan tetap bersifat ill-conditioned.
Namun demikian, nilai determinannya jauh berubah, yakni:
det = 10(20) – 20(11) = -20
# Partial Pivoting
Gunakan eliminasi Gauss untuk menyelesaikan: 0,0003 x1 + 3,0000 x2 = 2,0001 (i)
1,0000 x1 + 1,0000 x2 = 1,0000 (ii)
(Catatan: Penyelesaian eksak sistem persamaan ini adalah: x1 = 1/3 dan x2 = 2/3)
Penyelesaian:
Perhatikan bahwa elemen pivot yang pertama, a11 = 0,0003, mendekati nol...!!
Cara pertama:
Jika persamaan (i) dikalikan 1/(0,0003), maka diperoleh: x1 + 10000 x2 = 6667
yang dapat digunakan untuk mengeliminasi x1 dari persamaan (ii):
-9999 x2 = -6666 x2 = 2/3
x2 di atas selanjutnya dapat disubstitusikan balik ke persamaan (i) untuk mengevaluasi x1:
0,0003
2,0001 3(2/ 3)
1
−
x =
Namun demikian, hasil yang diperoleh sangat dipengaruhi oleh banyaknya significant figures
yang digunakan, seperti ditabelkan berikut ini:
Significant
figures x2 x1 Nilai mutlak percent
relative error untuk x1
3 0,667 -3,33 1099
4 0,6667 0,0000 100
5 0,66667 0,30000 10
6 0,666667 0,330000 1
7 0,6666667 0,3330000 0,1
Cara lain: Pivoting
Jika kedua persamaan saling dipertukarkan posisinya (supaya elemen pivot pertama
merupakan nilai yang terbesar dalam satu kolomnya), maka:
1,0000 x1 + 1,0000 x2 = 1,0000 (i.a)
0,0003 x1 + 3,0000 x2 = 2,0001 (ii.a)
Melalui eliminasi dan substitusi, diperoleh: x2 = 2/3.
Selanjutnya, x2 ini disubstitusikan balik sehingga diperoleh x1 sebesar:
1
1 (2/3)
1
−
x =
Pada beberapa nilai significant figures yang digunakan, diperoleh hasil perhitungan sbb.:
Significant
figures x2 x1 Nilai mutlak percent
relative error untuk x1
3 0,667 0,333 0,1
4 0,6667 0,3333 0,01
5 0,66667 0,33333 0,001
6 0,666667 0,333333 0,0001
7 0,6666667 0,3333333 0,00001
Berdasarkan hasil di atas, terlihat bahwa significant-figures-yang-digunakan mempunyai
pengaruh yang jauh kurang sensitif terhadap perhitungan (bandingkan dengan cara
sebelumnya...!!).
Kesimpulan: Pada kasus ini, strategi pivoting memberikan hasil yang jauh lebih baik.
# METODE GAUSS-JORDAN #
Metode eliminasi Gauss merupakan “separuh jalan” dari metode Gauss-Jordan, yakni metode yang
dilakukan untuk mengubah matriks bujur sangkar A menjadi sebuah matriks diagonal D. Proses
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 9 dari 19
eliminasi yang dilakukan terhadap matriks segitiga atas U agar menjadi matriks diagonal D dilakukan
melalui langkah-langkah yang serupa. Ilustrasi dalam contoh sebelumnya akan diproses lebih lanjut
menggunakan metode Gauss-Jordan. (lihat kembali contoh sebelumnya)
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
−
=
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
1
1
1
1
1 2 3 4
2 3 4 3
3 4 3 2
4 3 2 1
4
3
2
1
x
x
x
x
A x b
Matriks perbesaran:
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
−
1
1
1
1
1 2 3 4
2 3 4 3
3 4 3 2
4 3 2 1
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
0
12/ 7
1/ 4
1
0 0 0 5/ 3
0 0 12/ 7 10/ 7
0 7 / 4 3/ 2 5/ 4
4 3 2 1
A b U c
7
6
5/3
10/ 7
34 g = =
4
3
5/3
5/ 4
24 g = =
5
3
5/3
1
14 g = =
8
7
12/ 7
3/ 2
23 g = =
6
7
12/ 7
2
13 g = =
7
12
7 / 4
3
12 g = =
D s
Dengan demikian: 0
4
0
1 x = =
1
7 / 4
7 / 4
2 x = =
1
12/ 7
12/ 7
3 = −
−
x =
0
5/3
0
4 x = =
Karena penyelesaian Ax = b sudah dapat dicapai secara efisien dengan metode eliminasi Gauss,
penggunaan metode Gauss-Jordan bukanlah dalam penyelesaian sistem persamaan linier, melainkan
dalam perhitungan manual penentuan matriks kebalikan (atau invers) A-1 dari matriks bujur
sangkar A.
A . A-1 = I
(dapat ditentukan/dicari dengan metode Gauss Jordan)
Eliminasi Gauss
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
0
12 / 7
1/ 4
1
0 0 0 5/ 3
0 0 12 / 7 0
0 7 / 4 3/ 2 0
4 3 2 0
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
0
12/ 7
7 / 4
3
0 0 0 5/ 3
0 0 12/ 7 0
0 7 / 4 0 0
4 3 0 0
⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢
⎣
⎡
−
0
12/ 7
7 / 4
0
0 0 0 5/ 3
0 0 12/ 7 0
0 7 / 4 0 0
4 0 0 0
x = [0 1 -1 0]T
(sama dengan hasil yang diperoleh sebelumnya,
dengan metode eliminasi Gauss)
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 10 dari 19
Contoh: Tentukan invers dari matriks A =
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− 2 2 1
4 1 0
2 1 1
, dengan metode Gauss-Jordan!
Penyelesaian:
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− 2 2 1 0 0 1
4 1 0 0 1 0
2 1 1 1 0 0
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− − −
0 3 2 1 0 1
0 1 2 2 1 0
2 1 1 1 0 0
A I
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− −
− − −
0 0 4 5 3 1
0 1 2 2 1 0
2 1 1 1 0 0
(keadaan pada tahap akhir eliminasi Gauss)
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− −
− − −
0 0 4 5 3 1
0 1 2 2 1 0
2 1 1 1 0 0
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− −
− − −
−
0 0 4 5 3 1
0 1 0 1/ 2 1/ 2 1/ 2
2 1 0 1/ 4 3/ 4 1/ 4
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− −
− − −
−
0 0 4 5 3 1
0 1 0 1/ 2 1/ 2 1/ 2
2 0 0 1/ 4 1/ 4 1/ 4
Selanjutnya, baris pertama dibagi dengan 2, baris kedua dibagi dengan (-1), dan baris ketiga
dibagi dengan (-4), sehingga diperoleh:
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− −
−
−
0 0 1 5/ 4 3/ 4 1/ 4
0 1 0 1/ 2 1/ 2 1/ 2
1 0 0 1/8 1/8 1/8
I A-1
Verifikasi hasil:
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− 2 2 1
4 1 0
2 1 1
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− −
−
−
5/ 4 3/ 4 1/ 4
1/ 2 1/ 2 1/ 2
1/8 1/8 1/8
=
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− −
−
−
5/ 4 3/ 4 1/ 4
1/ 2 1/ 2 1/ 2
1/8 1/8 1/8
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
− 2 2 1
4 1 0
2 1 1
=
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
0 0 1
0 1 0
1 0 0
A A-1 = A-1 A = I
Jadi, matriks invers A adalah: A-1 =
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎣
⎡
− −
−
−
5/ 4 3/ 4 1/ 4
1/ 2 1/ 2 1/ 2
1/8 1/8 1/8
# METODE ITERATIF GAUSS-SEIDEL #
Metode Gauss-Seidel merupakan metode penyelesaian sistem persamaan aljabar linier yang dilakukan
secara iteratif. Berawal dari nilai-nilai tebakan awal untuk x, nilai awal ini selanjutnya diproses secara
berulang-ulang (atau iteratif), hingga diperoleh nilai x yang cukup dekat dengan nilai sesungguhnya
(atau memenuhi batas toleransi yang telah ditetapkan). Untuk sistem persamaan linier yang melibatkan
matriks A dan vektor b yang berorder n, diperlukan nilai tebakan awal x sejumlah n-1 buah.
Tinjaulah sistem persamaan linier (berorder 3) yang dinyatakan dalam bentuk matriks sbb.:
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
31 32 33
21 22 23
11 12 13
a a a
a a a
a a a
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
3
2
1
x
x
x
=
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
3
2
1
b
b
b
A x = b
Tahap pertama
Tahap pertama
Tahap kedua
Tahap kedua
eliminasi Gauss
eliminasi Gauss
Gauss-Jordan
Gauss-Jordan
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 11 dari 19
Dari baris 1, diperoleh: 3
11
13
2
11
12
11
1
1 x
a
x a
a
a
a
x = b − −
Dengan cara yang sama, x2 dan x3 dapat ditentukan dari baris 2 dan baris 3:
3
22
23
1
22
21
22
2
2 x
a
x a
a
a
a
x = b − − dan 2
33
32
1
33
31
33
3
3 x
a
x a
a
a
a
x = b − −
Pendekatan tersebut di atas digunakan dalam metode iteratif Gauss-Seidel. Jika dipilih 2 nilai tebakan
awal untuk x2 dan x3, maka x1 dapat dihitung. Berdasarkan x1 (baru) yang dihitung x3 (lama), maka x2
(baru) dapat dihitung. Berdasarkan x1 (baru) dan x2 (baru), maka x3 (baru) dapat dihitung. Demikian
seterusnya, sampai diperoleh harga-harga x1, x2, dan x3 yang konvergen. Prosedur iteratif ini
digambarkan dalam diagram alir berikut ini. Prosedur yang analog dapat diterapkan terhadap sistem
persamaan linier berorder n.
Persentase error perhitungan dapat dinyatakan sebagai:
100%
,
, , x
x
x x
Error
i iterasi sekarang
i iterasi sekarang i iterasi sebelumnya
i
−
=
Contoh:
Gunakan metode Gauss-Seidel untuk menyelesaikan sistem persamaan:
3 x1 – 0,1 x2 – 0,2 x3 = 7,85
0,1 x1 + 7 x2 – 0,3 x3 = -19,3
0,3 x1 – 0,2 x2 + 10 x3 = 71,4
Catatan: Nilai sebenarnya dari sistem persamaan ini adalah: x1 = 3; x2 = -2,5; dan x3 = 7
Penyelesaian:
Dari setiap persamaan, dapat diperoleh x1, x2, dan x3 sebesar:
3
7,85 0,1 0,2 2 3
1
x x x + +
=
7
19,3 0,1 0,3 1 3
2
x x x − − +
=
10
71,4 0,3 0,2 1 2
3
x x x − +
=
MULAI
Tebakan awal: x2,lama dan x3,lama
Toleransi: tol
baru lama lama x
a
x a
a
a
a
x b 3,
11
13
2,
11
12
11
1
1, = − −
baru baru lama x
a
x a
a
a
a
x b 3,
22
23
1,
22
21
22
2
2, = − −
baru baru baru x
a
x a
a
a
a
x b 2,
33
32
1,
33
31
33
3
3, = − −
Error ≤ tol ?
x1 = x1,baru
x2 = x2,baru
x3 = x3,baru
SELESAI
x2,lama = x2,baru
x3,lama = x3,baru
Ya
Tidak
ke iterasi
berikutnya
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 12 dari 19
Iterasi pertama:
Dengan menggunakan nilai tebakan awal untuk x2 dan x3 sebesar: x2 = x3 = 0, maka
diperoleh: 2,616667
3
7,85 0,1 (0) 0,2 (0)
1 =
+ +
x =
Gunakan x1 = 2,616667 ini dan x3 = 0 untuk menghitung x2:
2,794524
7
19,3 0,1 (2,616667) 0,3 (0)
2 = −
− − +
x =
Gunakan x1 = 2,616667 dan x2 = -2,794524 untuk menghitung x3:
7,005610
10
71,4 0,3 (2,616667) 0,2 ( 2,794524)
3 =
− + −
x =
Iterasi kedua:
Dengan cara yang sama, maka diperoleh:
2,990557
3
7,85 0,1 ( 2,794524) 0,2 (7,005610)
1 =
+ − +
x = = 0,31%; = 12,5% t a ε ε
2,499625
7
19,3 0,1 (2,990557) 0,3 (7,005610)
2 = −
− − +
x = = 0,015%; = 11,8% t a ε ε
7,000291
10
71,4 0,3 (2,990557) 0,2 ( 2,499625)
3 =
− + −
x = = 0,0042%; = 0,076% t a ε ε
Komentar: Perhatikan bahwa metode ini bersifat konvergen (menuju nilai-nilai x sebenarnya).
Langkah-langkah iterasi berikutnya dapat dilakukan untuk meningkatkan ketelitian dan
kedekatan dengan nilai sebenarnya.
Kriteria Konvergensi untuk Metode Gauss-Seidel
Metode Gauss-Seidel mirip dengan metode iterasi satu titik (dalam penentuan akar persamaan tak
linier), yang dapat memunculkan 2 kemungkinan persoalan, yakni: (1) kadang-kadang tidak
konvergen, atau (2) jika konvergen, maka laju konvergensinya lambat.
Untuk sistem persamaan linier Ax = b yang melibatkan n buah persamaan, konvergensi metode ini
akan selalu tercapai, jika: Σ
≠ =
>
n
i j
j
ii ij a a
1
(diagonally dominant)
Metode relaksasi merupakan modifikasi dari metode Gauss-Seidel yang dirancang untuk
meningkatkan konvergensi. Setelah masing-masing xi dihitung dengan persamaan-persamaan di
atas, nilai xi ini selanjutnya dimodifikasi dengan memberikan faktor pembobot terhadap hasil
iterasi sekarang dan iterasi sebelumnya; atau:
i baru i baru i lama x x x , , , = λ + (1− λ )
di mana: λ ≡ faktor pembobot (bernilai antara 0 – 2; biasanya ditentukan secara empirik)
Jika λ = 1, maka metode ini tidak dimodifikasi/direlaksasi
Jika 0 ≤ λ ≤ 1, maka modifikasinya tergolong dalam underrelaxation (biasanya digunakan untuk
membuat sistem yang tidak konvergen menjadi konvergen)
Jika 1 ≤ λ ≤ 2, maka modifikasinya tergolong dalam overrelaxation (biasanya digunakan untuk
mempercepat konvergensi sistem yang sudah konvergen). Pendekatan ini biasa disebut juga
successive or simultaneous overrelaxation (SOR).
# METODE ITERATIF JACOBI #
Metode Jacobi merupakan metode penyelesaian sistem persamaan aljabar linier secara iteratif yang
mirip dengan metode Gauss-Seidel. Berawal dari nilai-nilai tebakan awal untuk x, nilai awal ini
selanjutnya diproses secara berulang-ulang (atau iteratif), hingga diperoleh nilai x yang cukup dekat
dengan nilai sesungguhnya (atau memenuhi batas toleransi yang telah ditetapkan). Untuk sistem
persamaan linier yang melibatkan matriks A dan vektor b yang berorder n, diperlukan nilai tebakan
awal x sejumlah n buah.
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 13 dari 19
Tinjaulah sistem persamaan linier (berorder 3) yang dinyatakan dalam bentuk matriks sbb.:
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
31 32 33
21 22 23
11 12 13
a a a
a a a
a a a
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
3
2
1
x
x
x
=
⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢
⎣
⎡
3
2
1
b
b
b
A x = b
Seperti halnya metode Gauss-Seidel, metode Jacobi juga akan selalu menghasilkan penyelesaian yang
konvergen jika matriks koefisien A bersifat diagonally dominant. Jika A merupakan matriks bujur
sangkar yang bersifat diagonally dominant, maka dari baris 1, diperoleh: 3
11
13
2
11
12
11
1
1 x
a
x a
a
a
a
x = b − −
Dengan cara yang sama, x2 dan x3 dapat ditentukan dari baris 2 dan baris 3:
3
22
23
1
22
21
22
2
2 x
a
x a
a
a
a
x = b − − dan 2
33
32
1
33
31
33
3
3 x
a
x a
a
a
a
x = b − −
Jika dipilih 3 nilai tebakan awal untuk x1, x2, dan x3, maka x1 (baru), x2 (baru), dan x3 (baru) dapat
dihitung secara bersamaan. Demikian seterusnya, sampai diperoleh harga-harga x1, x2, dan x3 yang
konvergen. Prosedur iteratif ini digambarkan dalam diagram alir berikut ini. Prosedur yang analog
dapat diterapkan terhadap sistem persamaan linier berorder n.
Persentase error perhitungan dapat dinyatakan sebagai:
100%
,
, , x
x
x x
Error
i iterasi sekarang
i iterasi sekarang i iterasi sebelumnya
i
−
=
Contoh:
Ulangi soal sebelumnya dengan metode iteratif Jacobi!
Penyelesaian:
Jika digunakan nilai-nilai tebakan awal x, yaitu: x1 = x2 = x3 = 0
dan dengan cara yang mirip dengan sebelumnya, maka diperoleh hasil-hasil di bawah ini.
MULAI
Tebakan awal: x1,lama, x2,lama,
dan x3,lama ; Toleransi: tol
baru lama lama x
a
x a
a
a
a
x b 3,
11
13
2,
11
12
11
1
1, = − −
baru lama lama x
a
x a
a
a
a
x b 3,
22
23
1,
22
21
22
2
2, = − −
baru lama lama x
a
x a
a
a
a
x b 2,
33
32
1,
33
31
33
3
3, = − −
Error ≤ tol ?
x1 = x1,baru
x2 = x2,baru
x3 = x3,baru
SELESAI
x1,lama = x1,baru
x2,lama = x2,baru
x3,lama = x3,baru
Ya
Tidak
ke iterasi
berikutnya
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 14 dari 19
Tabulasi hasil perhitungan (dengan metode Jacobi):
Bandingkan hasil ini dengan hasil yang diperoleh dengan metode Gauss-Seidel sbb:
# METODE THOMAS UNTUK MENYELESAIKAN SISTEM PERSAMAAN
LINIER DENGAN MATRIKS KOEFISIEN TRIDIAGONAL #
Metode eliminasi Gauss dapat diterapkan untuk sistem persamaan linier yang matriks koefisiennya
berupa matriks tridiagonal. (Ingatlah kembali, apa yang dimaksud dengan matriks tridiagonal...!).
Metode terapan ini biasa disebut metode Thomas, yang alur penyelesaiannya akan digambarkan
dalam uraian berikut ini.
Tinjaulah sistem persamaan linier:
a11 x1 + a12 x2 = b1
a21 x1 + a22 x2 + a23 x3 = b2
a32 x2 + a33 x3 + a34 x4 = b3
.........................................................................
an,n-1 xn-1 + ann xn = bn
Atau, dapat dituliskan menjadi:
d1 x1 + e1 x2 = b1
c2 x1 + d2 x2 + e2 x3 = b2
c3 x2 + d3 x3 + e3 x4 = b3
.........................................................................
cn xn-1 + dn xn = bn
Dalam bentuk perkalian matriks dan vektor, persamaan di atas dapat dituliskan sebagai:
⎥ ⎥ ⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢ ⎢ ⎢
⎣
⎡
=
⎥ ⎥ ⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢ ⎢ ⎢
⎣
⎡
⎥ ⎥ ⎥ ⎥ ⎥ ⎥
⎦
⎤
⎢ ⎢ ⎢ ⎢ ⎢ ⎢
⎣
⎡
n n n n b
b
b
b
x
x
x
x
c d
c d e
c d e
d e
... ...
0 0 0
... ... ... ... ...
0 0
0 0
0 0 0
3
2
1
3
2
1
3 3 3
2 2 2
1 1
... (*)
Matriks tridiagonal
Dengan menerapkan langkah-langkah yang telah dipelajari dalam metode eliminasi Gauss (pada
materi sebelumnya), penyelesaian dari sistem persamaan linier (*) adalah sbb.:
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 15 dari 19
n n x = γ
dan:
i
i i
i i
x e x
β
= γ − +1 (i = n-1, n-2, ..., 1)
dengan γi dan βi yang ditentukan melalui serangkaian perhitungan berurutan sbb.:
1 1 β = d
1
1
1 β
γ = b
1
1
−
= − −
i
i i
i i
d c e
β
β (i = 2, 3, ..., n)
i
i i i
i
b c
β
γ
γ −1 −
= (i = 2, 3, ..., n)
Contoh:
Selesaikan sistem persamaan berikut:
2 x1 – 3 x2 = -4
x1 + 2 x2 – x3 = 2
4 x2 – x3 + x4 = 9
2 x3 – x4 = 2
Penyelesaian:
Sistem persamaan linier di atas memuat matriks koefisien tridiagonal. Mula-mula besaranbesaran
ci, di, ei, dan bi akan didefinisikan lebih dahulu, untuk menyesuaikan dengan notasi
yang diuraikan dalam metode Thomas tersebut di atas.
i ci di ei bi
1 0 2 -3 -4
2 1 2 -1 2
3 4 -1 1 9
4 2 -1 0 2
Dengan metode Thomas, maka diperoleh hasil-hasil sbb.:
β1
= d1 = 2
γ1
=
1
1β
b =
2
− 4 = -2
β2
=
1
2 1
2 β
d − c e =
2
2 1. ( 3)
−
− = 3,5
γ2
=
2
2 2 1
β
b − c γ = 3,5
2 −1. (−2) = 1,14
β3
=
2
3 2
3 β
d − c e = 3,5
1 4 . ( 1)
−
− − = 0,143
γ3
=
3
3 3 2
β
b − c γ = 0,143
9 − 4 .1,14 = 31
β4
=
3
4 3
4 β
d − c e = 0,143
−1− 2 .1 = -15
γ4
=
4
4 4 3
β
b − c γ =
15
2 2 . 31
−
− = 4
Dengan demikian:
x4 = γ
4 = 4
x3 =
3
3 4
3 β
γ − e x = 0,143
31− 1. 4 = 3
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 16 dari 19
x2 =
2
2 3
2 β
γ − e x =
3,5
1,14 ( 1) (3)
−
− = 2
x1 =
1
1 2
1 β
γ − e x =
2
2 ( 3) (2)
−
− − = 1
Jadi, penyelesaian sistem persamaan linier tersebut di atas adalah:
x1 = 1; x2 = 2; x3 = 3; dan x4 = 4
# CONTOH APLIKASI #
Neraca Massa Rangkaian Proses
Umpan berupa zat A murni dengan laju 100 kmol/jam.
Kendala: 1. 80% dari A dan 40% dari B di dalam alur 2 didaur ulang (recycle).
2. Perbandingan mol A terhadap mol B di dalam alur 1 adalah 5:1
Neraca massa (dalam kmol/jam):
Di sekitar pencampur:
NA1 = NA3 + 100 atau: NA1 − NA3 = 100 ... (1)
NB1 = NB3 atau: NB1 − NB3 = 0 ... (2)
(NA1 menyatakan laju alir molar A di dalam alur 1, dst.)
Di sekitar reaktor:
NA2 = NA1 − r atau: − NA1 + NA2 + r = 0 ... (3)
NB2 = NB1 + r atau: − NB1 + NB2 − r = 0 ... (4)
(r menyatakan laju reaksi)
Di sekitar pemisah:
NA3 + NA4 = NA2 atau: − NA2 + NA3 + NA4 = 0 ... (5)
NB3 + NB4 = NB2 atau: − NB2 + NB3 + NB4 = 0 ... (6)
Berdasarkan kendala 1:
NA3 = 0,8 NA2 atau: − 0,8 NA2 + NA3 = 0 ... (7)
NB3 = 0,4 NB2 atau: − 0,4 NB2 + NB3 = 0 ... (8)
Berdasarkan kendala 2:
NA1 = 5 NB1 atau: NA1 − 5 NB1 = 0 ... (9)
Berdasarkan penjabaran neraca massa di atas, dihasilkan 9 buah persamaan linier dengan 9 variabel
yang tak diketahui (yakni NA1, NB1, NA2, NB2, NA3, NB3, NA4, NB4, dan r)
Dengan demikian, terbentuk sistem persamaan linier yang dapat diselesaikan secara simultan!
Pemisah
Produk
Umpan
A → B
Pencampur Reaktor
4
3
1 2
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 17 dari 19
(Sebagai latihan, susunlah ke-9 persamaan di atas dalam bentuk perkalian matriks dan
vektor, dan selanjutnya silakan Anda selesaikan sendiri sistem persamaan linier tersebut.
Untuk mengecek jawaban Anda, hasil perhitungan dari persoalan neraca massa ini adalah:
NA1 = 227,3 kmol/jam; NB1 = 45,45 kmol/jam; NA2 = 159,1 kmol/jam; NB2 = 113,6 kmol/jam;
NA3 = 127,3 kmol/jam; NB3 = 45,45 kmol/jam; NA4 = 31,82 kmol/jam; NB4 = 68,18 kmol/jam;
dan r = 68,18 kmol/jam)
# LATIHAN SOAL #
1. Tuliskan sekumpulan persamaan berikut ini dalam bentuk matriks:
2 3 30 = 2 x + 6 x
2 1 20 = 3 x + 8 x
1 3 10 = x + x
Tuliskan tranpose matriks koefisiennya!
2. Gunakan metode grafis untuk menyelesaikan sistem persamaan berikut:
2 6 18 1 2 x − x = −
8 40 1 2 − x + x =
Cek jawaban Anda dengan mensubstitusikannya kembali ke dalam persamaan-persamaan di atas.
3. Untuk sekumpulan persamaan berikut:
2 5 1 2 3 x + x =
2 2 1 1 2 3 x + x + x =
3 2 1 2 x + x =
(a) Hitunglah nilai determinannya
(b) Gunakan cara Cramer untuk menentukan harga x1, x2, dan x3.
(c) Substitusikan jawaban Anda (pada butir (b)) ke dalam persamaan-persamaan semula, untuk
mengecek kebenarannya.
4. Gunakan eliminasi Gauss untuk menyelesaikan sistem persamaan berikut:
4 2 1 2 3 x + x − x = −
5 2 4 1 2 3 x + x + x =
6 6 1 2 3 x + x + x =
Lakukan partial column pivoting dan cek kebenaran jawaban Anda dengan mensubstitusikan
kembali ke dalam persamaan-persamaan semula.
5. Selesaikan sistem persamaan:
3 1 2 3 x + x − x = −
6 2 2 2 1 2 3 x + x + x =
3 4 1 1 2 3 − x + x + x =
dengan metode:
(a) eliminasi Gauss naive
(b) eliminasi Gauss dengan partial pivoting
(c) Gauss-Jordan naive
(d) Gauss-Jordan dengan partial pivoting
6. Sistem persamaan berikut ini dirancang untuk menentukan konsentrasi (c, dalam g/m3) dalam
serangkaian reaktor sebagai fungsi jumlah massa yang diumpankan ke dalam masing-masing
reaktor (pada ruas kanan persamaan, dalam g/hari):
17 2 3 500 1 2 3 c − c − c =
5 21 2 200 1 2 3 − c + c − c =
5 5 22 30 1 2 3 − c − c + c =
(a) Gunakan matriks invers untuk menentukan harga konsentrasi c1, c2, dan c3.
(b) Tentukan c1, c2, dan c3 dengan metode Gauss-Seidel, hingga tercapai εs = 5%.
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 18 dari 19
7. Gunakan metode Gauss-Seidel dengan relaksasi (dengan λ = 0,90 dan εs = 5%) untuk menyelesaikan
sistem persamaan:
5 12 80 1 3 − x + x =
4 2 1 2 3 x − x − x = −
6 8 45 1 2x + x =
(Jika perlu, lakukan penyusunan ulang terhadap persamaan untuk mencapai konvergensi).
8. Selesaikan 3 buah persamaan linier berikut ini dengan metode eliminasi Gauss, dan tunjukkan
semua langkah perhitungannya.
3 2 3 1 2 3 x − x + x =
2 4 2 2 1 2 3 x + x − x =
4 2 3 12 1 2 3 x − x − x = −
Ulangilah jika menggunakan metode iteratif Gauss-Seidel (atau Jacobi).
9. Selesaikan sistem persamaan linier berikut ini dengan metode eliminasi Gauss, dan tunjukkan
semua langkah perhitungannya.
2 3 2 3 x − x =
2 3 1 1 2 3 x − x + x =
3 2 5 10 1 2 3 x − x + x =
10. Selesaikan 3 buah persamaan linier berikut ini dengan metode eliminasi Gauss, dan tunjukkan
semua langkah perhitungannya.
3 3 2 1 2 3 x − x + x =
5 3 6 1 2 3 x + x + x =
2 2 1 2 3 x + x − x =
11. Selesaikan sistem persamaan linier berikut ini dengan metode Thomas, dan tunjukkan semua
langkah perhitungannya.
1 1 2 x − x =
2 3 1 1 2 3 x + x − x =
3 4 2 3 2 3 4 x − x + x =
5 3 3 4 x − x =
Ulangilah jika menggunakan metode iteratif Gauss-Seidel (atau Jacobi).
12. Sebuah kolom ekstraksi yang memiliki 4 tahap ideal/teoretik akan digunakan untuk mengekstraksi
asam format yang terlarut dalam metil isobutil keton (MIBK) dengan menggunakan air sebagai
pengekstrak. Operasi akan dilakukan dengan mode berlawanan arah (counter current); fasa organik
diumpankan dari dasar kolom dengan laju 5 m3 MIBK/jam dan konsentrasi asam format 1,0
mol/liter MIBK, sedangkan fasa akuatik dialirkan dari puncak kolom dengan laju alir 2,45 m3
air/jam serta sama sekali tidak mengandung asam format. Operasi berlangsung isotermal dan harga
rasio kesetimbangan konsentrasi asam format di dalam MIBK terhadap konsentrasi asam format di
dalam air adalah K = 0,445. MIBK dan air dapat dianggap sama sekali tidak larut satu sama lain. 1
m3 = 1000 liter. Tentukan konsentrasi asam format pada aliran-aliran yang keluar dan masuk tiap
pelat/tahap ideal dalam kolom ekstraksi tersebut.
13. Sebuah kilang minyak menghasilkan 3 jenis elpiji: A, B, dan C, dengan komposisi seperti disajikan
pada tabel berikut ini (angka-angka yang ditunjukkan menyatakan %-volume komponen):
Komponen A ElBpi ji C
C2 5,0 - -
C3 90,0 10,0 -
i-C4 5,0 85,0 8,0
n-C4 - 5,0 80,0
i-C5 - - 12,0
Jumlah 100,0 100,0 100,0
dy/analisis numerik/penyelesaian sistem persamaan aljabar linier/februari 2008/halaman 19 dari 19
Dari sebuah perusahaan calon pelanggan, kilang minyak tersebut menerima tawaran kontrak jual
beli 106 m3 elpiji yang disyaratkan harus mengandung: 31,2%-v propana, 53,4%-v isobutana, dan
12,6%-v n-butana. Kadar etana dan i-C5 tidak disyaratkan. Kilang minyak tersebut bermaksud
memenuhi kontrak ini dengan cara mencampurkan elpiji A, elpiji B, dan elpiji C.
Tentukan kebutuhan masing-masing elpiji A, B, dan C (dalam m3) agar persyaratan kontrak dapat
dipenuhi dengan baik!
14. Gambar berikut ini memperlihatkan sistem tiga buah reaktor yang terhubungkan oleh pipa-pipa.
Laju perpindahan setiap komponen/zat melalui setiap pipa merupakan perkalian dari laju alir
volumetriknya (Q, dalam satuan m3/detik) dengan konsentrasinya yang keluar dari masing-masing
reaktor (C, dalam satuan mg/m3). Sistem proses berada dalam kondisi steady, yang berarti bahwa
laju massa masuk ke setiap reaktor sama dengan laju massa keluar.
Susunlah neraca massa pada masing-masing reaktor dan selanjutnya selesaikanlah sistem
persamaan aljabar linier yang terbentuk (untuk menentukan harga-harga C1, C2, dan C3).
15. Tinjaulah proses perpindahan panas secara konduksi melalui dinding sebuah pipa. Suhu dinding
bagian dalam: 200oC, suhu dinding bagian luar: 80oC, dan ketebalan dinding pipa: 0,05 meter. Jarijari
bagian dalam pipa: 0,05 meter (ro) dan jari-jari bagian luar pipa: 0,1 meter. Persamaan
diferensial yang dijabarkan berdasarkan neraca panas dalam dinding pipa dan menggambarkan
distribusi suhu sepanjang dinding pipa dinyatakan sbb.:
0 2
2
= ⎟
⎟⎠
⎞
⎜ ⎜⎝
⎛
+
dr
r d T
dr
dT
Pendekatan beda hingga (finite difference) yang diterapkan untuk menyelesaikan persamaan di atas
menghasilkan serangkaian persamaan aljabar linier berikut:
0
2
2 1
2
1
1 1 = ⎟ ⎟⎠
⎞
⎜ ⎜⎝
⎛
Δ
+ − +
Δ
− ⎟
⎟⎠
⎞
⎜ ⎜⎝
⎛
Δ
+ i + i i − T
r
T r
r
T r
r
r
untuk: i = 2, 3, ..., 10
dengan:
r = r + Δ r (i −1) = 0,05 + 0,005 (i −1) o
dan
T1 = 200oC dan T11 = 80oC
Tentukan distribusi suhu di sepanjang dinding pipa menggunakan metode Thomas! (Jangan lupa
bahwa r tidaklah konstan)
☺☺☺ Selamat Belajar!!! ☺☺☺
400
mg/detik
200
mg/detik
3
2
1
Q13C1
Q12C1 Q23C2
Q21C2
Q33C3
Q33 = 120
Q13 = 40
Q12 = 80
Q23 = 60
Q21 = 20
Keterangan:
Tidak ada komentar:
Posting Komentar