Algoritma Thomas

Unduh sebagai docx, pdf, atau txt
Unduh sebagai docx, pdf, atau txt
Anda di halaman 1dari 2

Algoritma Thomas

Diberikan sistem persamaan linier berikut

Tu=q
Dimana T adalah matriks tridiagonal, u adalah vektor variabel, dan q adalah vektor
konstanta.
SPL tersebut dapat diselesaikan dengan menggunakan Algoritma Thomas yang
merupakan penyederhanaan dari Eliminasi Gauss.
Algoritma Thomas adalah sebagai berikut.
1. Memfaktorkan matriks T menjadi matriks segitiga bawah B dan matriks
segitiga atas A, dimana B dan A adalah matriks bidiagonal. Langkah ini
dikenal dengan LU decomposition.

t 11 t 12
t 21 t 22

0 0
t 23 0


0 0 t nn1

][

1 0 0 0
b
1 0 0
= 2

0 0 bn 1
t nn

][

d 1 a1
0 d2


0 0 0 dn

T
Dimana

0 0
a2 0

d 1=t 11 , a1=t 12

Untuk

i=2, 3, , n1

bi=t i ,i 1 /d i1
d i=t i ,ibi ai1
ai=t i ,i +1

bn =t n , n1 /d n1
d n=t n ,n bn an 1
2. Substitusi maju dengan melibatkan matriks B;

1 0 0 0
b2 1 0 0

0 0 bn 1
Dimana

][ ] [ ]
r1
q1
r2
q
= 2

rn
qn

r 1=q 1
r i=qib i r i1 ,

untuk

i=2, 3, , n

3. Substitusi mundur dengan melibatkan matriks A;

d 1 a1
0 d2

0 0
a2 0

Br =q

Au=r

][ ] [ ]

u1
r1
u2
r
= 2


u
r
0 0 0 dn n
n
Dimana

un=r n /d n
ui=(r iai ui +1)/d i ,

untuk

i=n1, ,3, 2,1

FALSE

Reff: http://www.cfm.brown.edu/people/gk/chap6/node13.html

Diagram Alir Program


TRUE

Start

Deklarasi variabel, dsb


nx=10
nt=20
k=nx/2

Input matriks koefisien

Input nilai awal

Input syarat batas

LU decomposition

j=0
j=j+1

j <= nt-1

Kalkulasi nilai vektor konstanta


sebanyak k

Anda mungkin juga menyukai