Petunjuk Praktikum - Metode Numerik 19
2. Metode Gauss-Jordan
Tugas 2: Buatlah Program untuk menyelesaikan SPL diatas dengan
algoritma Gauss-Jordan.
Praktikum 4. Metode Hitungan langsung - Dekomposisi LU
Praktikum 4. Metode Hitungan langsung - Dekomposisi LU
B. Metode Dekomposisi
Salah satu metode penyelesaian SPL AX=B yang lain adalah yang
biasa dikenal dengan sebutan Dekomposisi LU. Prinsip Metode ini adalah
menfaktorkan Matriks A menjadi Suatu perkalian 2 matrik yaitu matriks L
(Matriks segitiga bawah) dan U (matriks segitiga atas). (A=LU)
a11 a12 a1n x1 b1
a21 a22 a2 n x2 b2
a31 a32 a3n x3 = b3
a an 2 ann xn bn
n1
Matriks A dapat difaktorkan sebagai berikut:
a11 a12 a13 a14 l11 0 0 0 u11 u21 u31 u41
a a24 l21 l22 0 0
21 a22 a23
= 0 u22 u32 u42
a31 a32 a33 a34 l31 l32 l33 0 0 0 u33 u43
a41 a42 a43 a44 l41 l42 l43 l44 0 0 0 u44
M. Ziaul Arif, [Link] , Rusli Hidayat, [Link] – Universitas Jember
Petunjuk Praktikum - Metode Numerik 20
Langkah Penyelesaian SPL AX=B sebagai berikut:
AX=B
LUX=B
Misal: UX=Y
Maka,
LY=B, mencari matriks Y dengan penyulihan maju
UX=Y, mencari matriks X atau penyelesaian SPL dengan penyulihan mundur.
1. Metode Doolittle
Prinsip Metode Dekomposisi LU ini adalah bahwa diagonal utama
matriks L-nya bernilai 1 dan diagonal matriks U tak nol. sehingga
bentuknya sebagai berikut:
a11 a12 a13 a14 1 0 0 0 u11 u12 u13 u14
a a24 l21 1 0 0
21 a22 a23
= 0 u22 u23 u24
a31 a32 a33 a34 l31 l32 1 0 0 0 u33 u34
a41 a42 a43 a44 l41 l42 l43 1 0 0 0 u44
u11 u12 u13 u14
l u
21 11 l21u12 + u22 l21u13 + u23 l21u14 + u24
l31u11 l31u12 + l32u22 l31u13 + l32u23 + u33 l31u14 + l32u24 + u34
l41u11 l41u12 + l42u22 l41u13 + l42u23 + l43u33 l41u14 + l42u24 + l43u34 + u44
M. Ziaul Arif, [Link] , Rusli Hidayat, [Link] – Universitas Jember
Petunjuk Praktikum - Metode Numerik 21
Contoh 4: Buatlah program untuk menyelesaikan SPL berikut dengan metode
Doolittle.
x1 − x2 + 2 x3 − x4 = − 8
x1 + x2 + x3 =− 2
2 x1 − 2 x2 + 3x3 − 3x4 = − 20
x1 − x2 + 4 x3 + 3x4 = 4
M. Ziaul Arif, [Link] , Rusli Hidayat, [Link] – Universitas Jember
Petunjuk Praktikum - Metode Numerik 22
Penyelesaian:
Berikut Hasil running program pada MATLAB.
Command Window
>> Doolittle
masukkan matriks A yang akan diproses (Beri []) :[1 -1 2 -1;1 1 1
0;2 -2 3 -3;1 -1 4 3]
Masukkan matriks B (beri []) :[-8 -2 -20 4]
Matriks A
A =
1 -1 2 -1
1 1 1 0
2 -2 3 -3
1 -1 4 3
Matriks L
L =
1 0 0 0
1 1 0 0
2 0 1 0
1 0 -2 1
Matriks U
U =
1 -1 2 -1
0 2 -1 1
0 0 -1 -1
0 0 0 2
Matriks Y
Y =
-8
6
-4
4
Matriks X
ans =
-7
3
2
2
M. Ziaul Arif, [Link] , Rusli Hidayat, [Link] – Universitas Jember
Petunjuk Praktikum - Metode Numerik 23
Listing program
Berikut Listing program metode Doolittle. Simpan pada MATLAB dengan
nama ”doolittle”.
%%
clear
clc
disp(' METODE DOOLITTLE ')
disp('Press Enter to continue')
pause
clc
%%
A=input('masukkan matriks A yang akan diproses (Beri []) :');
B1=input('Masukkan matriks B (beri []) :');
B=B1';
[m,n]=size(A);
ordo=length(A);
if m~=n
disp('Matriks harus persegi')
break
end
U=zeros(m);
L=zeros(m);
for j=1:m
L(j,j)=1;
end
for j=1:m
U(1,j)=A(1,j);
end
for i=2:m
for j=1:m
for k=1:i-1
s1=0;
if k==1
s1=0;
else
Algoritma pencarian L
for p=1:k-1
s1=s1+L(i,p)*U(p,k);
end
end
L(i,k)=(A(i,k)-s1)/U(k,k);
end
M. Ziaul Arif, [Link] , Rusli Hidayat, [Link] – Universitas Jember
Petunjuk Praktikum - Metode Numerik 24
end
for k=i:m
s2=0;
for p=1:i-1
s2=s2+L(i,p)*U(p,k);
end Algoritma pencarian U
U(i,k)=A(i,k)-s2;
end
end
end
Y1(1) = B(1);
for k=2:ordo
sigma1=0;
for j=1:k-1
sigma1=sigma1+L(k,j)*Y1(j); Penyulihan Maju
Y1(k)= B(k)-sigma1;
end
end
Y=Y1';
X(ordo)=Y(ordo)/U(ordo,ordo);
for k=ordo-1:-1:1
sigma2=0;
for j=k+1:ordo; Penyulihan mundur
sigma2=sigma2+U(k,j)*X(j);
X(k)=(Y(k)-sigma2)/U(k,k);
end
end
disp('Matriks A')
A
disp('Matriks L')
L
disp('Matriks U')
U
disp('Matriks Y')
Y
disp('Matriks X')
X'
M. Ziaul Arif, [Link] , Rusli Hidayat, [Link] – Universitas Jember
Petunjuk Praktikum - Metode Numerik 25
2. Metode Crout
Perbedaan Metode Crout dengan Metode Doolitte adalah ada pada
matriks U. Matriks U pada metode Crout diagonal utamanya bernilai 1 dan
diagonal matriks L tak nol. sehingga bentuknya sebagai berikut:
Tugas 3: Buatlah Program untuk menyelesaikan SPL berikut dengan Metode
Crout.
M. Ziaul Arif, [Link] , Rusli Hidayat, [Link] – Universitas Jember