Tugas Metode Eliminasi Gauss Dengan Matlab 2q392r

Tugas 2

Teknik Komputasi

Metode Eliminasi Gauss dengan Matlab

Dosen

Nama NIM Kelas

: Dr. Ir. Nazori Az, M.T

: Fransiscus Xaverius Eko Budi Kristanto : 1111600126 : XA

MAGISTER ILMU KOMPUTER UNIVERSITAS BUDI LUHUR JAKARTA 2012

Metode Eliminasi Gauss dengan Matlab 2012 Soal: Lakukan solusi untuk 2 contoh persamaan linear berikut dengan menggunakan eliminasi Gauss-Jordan dan buatlah script dalam MATLAB, dari pemrograman sampai keluarannya dengan metode eliminasi Gauss. 1. Persamaan linear 1: P1 :

x1

+

x2

P2 :

2 x1

+

x2

P3 :

3x1

− x2

P4 :

− x1

+ 2 x2

+

3x 4 = 4

− x3

+

x4

− x3

+ 2 x4 = − 3

+ 3 x3

= 1

− x4

= 4

2. Persamaan linear 2: P1 :

x1

P2 :

2 x1

P3 : P4 :

−

x2

+ 2 x3 −

− 2 x2

+ 3 x3

x1

+

x2

+

x1

−

x2

+ 4 x3

x4 = − 8

− 3 x4 = − 20 = −2

x3 + 3 x4

= 4

Jawab: 1. Penyelesaian Persamaan Linear 1 dengan Metode Eliminasi Gauss: Persamaan linear: P1 :

x1

+

x2

P2 :

2 x1

+

x2

P3 :

3x1

− x2

P4 :

− x1

+ 2 x2

+

3x 4 = 4

− x3

+

x4

− x3

+ 2 x4 = − 3

+ 3 x3

= 1

− x4

= 4

Bentuk Matriks lengkapnya:

Script Program Matlab:

1 2 � 3 −1

1 1 −1 2

0 −1 −1 3

3 1 2 −1

4 1 � −3 4

% Penyelesaian Persamaan Linear 1 dengan Metode Eliminasi Gauss clear all clc A=[1 1 0 3 4;2 1 -1 1 1;3 -1 -1 2 -3;-1 2 3 -1 4]; % Data matriks disp('Matriks A:') A disp('Jumlah Persamaan:') 1

Metode Eliminasi Gauss dengan Matlab 2012 n=4 % jumlah persamaan pause %===Proses Triangularisasi=== for j=1:(n-1) %---mulai proses pivot--if (A(j,j)==0) for p=1:n+1 u=A(j,p); v=A(j+1,p); A(j+1,p)=u; A(j,p)=v; end end %---akhir proses pivot--jj=j+1; for i=jj:n m=A(i,j)/A(j,j); for k=1:(n+1) A(i,k)=A(i,k)-(m*A(j,k)); end end end disp('Matriks A hasil Proses Triangularisasi:') A pause %===Akhir Proses Triangularisasi=== %---Proses Substitusi mundur--x(n,1)=A(n,n+1)/A(n,n);

x

for i=n-1:-1:1 S=0; for j=n:-1:i+1 S=S+A(i,j)*x(j,1); end x(i,1)=(A(i,n+1)-S)/A(i,i); end

2

Metode Eliminasi Gauss dengan Matlab 2012 Hasil: Matriks A: A= 1 1 0 3 4 2 1 -1 1 1 3 -1 -1 2 -3 -1 2 3 -1 4 Jumlah Persamaan: n= 4 Matriks A hasil Proses Triangularisasi: A= 1 1 0 3 4 0 -1 -1 -5 -7 0 0 3 13 13 0 0 0 -13 -13 x= -1 2 0 1

Dengan demikian, maka penyelesaiannya adalah: x1 = -1, x2 = 2, x3 = 0 dan x4 = 1

3

Metode Eliminasi Gauss dengan Matlab 2012 2. Penyelesaian Persamaan Linear 2 dengan Metode Eliminasi Gauss: Persamaan linear: P1 :

x1

P2 :

2 x1

P3 : P4 :

−

x2

+ 2 x3 −

− 2 x2

+ 3 x3

x1

+

x2

+

x1

−

x2

+ 4 x3

x4 = − 8

− 3 x4 = − 20 = −2

x3 + 3 x4

= 4

Bentuk Matriks lengkapnya: 1 −1 2 −2 � 1 1 1 −1

2 −3 1 4

−1 −3 0 3

−8 −20 � −2 4

Script Program Matlab: % Penyelesaian Persamaan Linear 2 dengan Metode Eliminasi Gauss clear all clc A=[1 -1 2 -1 -8;2 -2 3 -3 -20;1 1 1 0 -2;1 -1 4 3 4]; % Data matriks disp('Matriks A:') A disp('Jumlah Persamaan:') n=4 % jumlah persamaan pause %===Proses Triangularisasi=== for j=1:(n-1) %---mulai proses pivot--if (A(j,j)==0) for p=1:n+1 u=A(j,p); v=A(j+1,p); A(j+1,p)=u; A(j,p)=v; end end %---akhir proses pivot--jj=j+1; for i=jj:n m=A(i,j)/A(j,j); for k=1:(n+1) 4

Metode Eliminasi Gauss dengan Matlab 2012 A(i,k)=A(i,k)-(m*A(j,k)); end

end end disp('Matriks A hasil Proses Triangularisasi:') A pause %===Akhir Proses Triangularisasi=== %---Proses Substitusi mundur--x(n,1)=A(n,n+1)/A(n,n);

x

for i=n-1:-1:1 S=0; for j=n:-1:i+1 S=S+A(i,j)*x(j,1); end x(i,1)=(A(i,n+1)-S)/A(i,i); end

5

Metode Eliminasi Gauss dengan Matlab 2012 Hasil: Matriks A: A= 1 -1 2 -2 1 1 1 -1

2 -1 -8 3 -3 -20 1 0 -2 4 3 4

Jumlah Persamaan: n= 4 Matriks A hasil Proses Triangularisasi: A= 1 -1 2 -1 -8 0 2 -1 1 6 0 0 -1 -1 -4 0 0 0 2 4 x= -7 3 2 2 Dengan demikian, maka penyelesaiannya adalah: x1 = -7, x2 = 3, x3 = 2 dan x4 = 2

6