JURNAL PRAKTIKUM

KOMPUTASI NUMERIK TEKNIK KIMIA

MODUL 3
(METODE BAIRSTOW)

NAMA : ADITYA MARDIANSYAH

NRP : 02211740000058
ASISTEN : HANA PUTRI TUADAYANI
HARI/TANGGAL PRAKTIKUM : SENIN/ 07 OKTOBER 2019

LABORATORIUM SIMULASI DAN KOMPUTASI

DEPARTEMEN TEKNIK KIMIA
INSTITUT TEKNOLOGI SEPULUH NOPEMBER
SURABAYA
2019
METODE BAIRSTOW 2

I. Metode Bisection
I.1. Tujuan Praktikum
Praktikum ini bertujuan untuk menyelesaikan persamaan non linear dengan metode
bairstow

I.2. Dasar Teori

Dalam menyelesaikan fungsi kuadrat,sudah terdapat rumus untuk mencari akar akar
nya.Masalah yang dihadapi adalah bagaimana mencari akar akar polinomial derajat
n.Dengan metode Bairstow yabg didasari atas pembagian polinomial dengan suatu faktor
kuadrat,akar persamaan kuadrat hasil operasi fungsi bairstow dapay diselesaikan dengan
metde solusi persamaan kuadrat.Metode ini digunakan ntu mencari akar akar persamaan
polinomial dengan menentukan faktor faktor kuadratnya.Dalam hal ini polinomial ditulis
sebagai berikut :
Pn(x) = a1xn+a2xn-1+.......+anx+an+1
Bila x2-rx-s adalah faktor kuadratis yang dicari,maka :
Pn(x) = a1xn+a2xn-1+.......+anx+an+1
Pn(x) = (x2-rx-s)( b1xn-2 + b2xn-3+.......+ bn-1)( bn(x-r) + bn+1)
Apabila akar kuadratis sudah benar,maka bn = 0 dan bn+1 = 0
∂bn ∂bn
0 = b n+ ∂ r ∆ r + ∂ s ∆ s..........................................................................................
(1)

+ ∂ b n+1 ∂ b n+ 1
0 = b n+1 ∆ r+ ∆ s..............................................................................
∂r ∂s
(2)
∂bn ∂bn
Dimana, =¿ cn-1 =¿ cn-2
∂r ∂s
∂ b n+1 ∂ b n+1
=¿ cn =¿ cn-1
∂r ∂s
Maka persamaan (1) dan (2) bisa ditulis :

cn-1 ∆ r + cn-2 ∆ s = - bn.............................................................................................

(3)

cn ∆ r + cn-1 ∆ s = - bn+1.............................................................................................
(4)
LABORATORIUM SIMULASI DAN KOMPUTASI
DEPARTEMEN TEKNIK KIMIA
INSTITUT TEKNOLOGI SEPULUH NOPEMBER
2019
METODE BAIRSTOW 3

Dari Persamaan (3) dan 4 ∆ r dan ∆ s dapat diperoleh dengan :

−bn c n−1+ bn+1 c n−2
∆r=
( cn −1 )2−c n c n−2
−bn+ 1 c n−1 +bn c n
∆ s=
(c n−1)2−c n c n−2
Dengan didaptkan nya ∆ r dan ∆ s ,maka r* dan s* dapat diisi dengan
r* = r + ∆ r
s* = s + ∆ s
I.3. Algoritma
Prosedur metode bairstow secara sistematis adalah :
1. Pilih harga pendekatan awal r1,s1 dan toleransi
2. Tentukan b(i) dan c(i)
b = a1 c1 = b1
b2 = a2 + rb1 c2 = b2 + rc1
b3 = a3 + rb2 + sb1 c3 = b3 + rc2 + sc1
bn = an + rbn-1 + sbn-2 cn = bn + rcn-1 + scn-2
bn+1 = an+1 + rbn-1 + sbn-2 cn+1 = bn+1 + rcn-1 + scn-2
3. Tentukan DENOM
DENOM = (Cn-1)2 – CnCn-2
4. Bila DENOM = 0 maka R = R+1, S = S+1 dan kembali ke tahap 3
Bila DENOM ≠ 0 mak lanjut ke tahap 5
5. Tentukan DEL R dan DEL S,yaitu :
DEL R = [-bncn-1 + bn+1cn-2] / DENOM
DEL S = [-bn+1cn-1 + bncn] / DENOM
6. Tentukan R baru atau S baru
R baru = R lama + DEL R
S baru = S lama + DEL S
7. Bila |DEL S| + |DEL R|≤ Toleransi, maka perhitungan dihentikan
Serta,harga S dan R terakhir adalah R = S yang dicari.
Bila |DEL S| + |DEL R|>¿ Toleransi, kembali ke tahap 2

LABORATORIUM SIMULASI DAN KOMPUTASI

DEPARTEMEN TEKNIK KIMIA
INSTITUT TEKNOLOGI SEPULUH NOPEMBER
2019
METODE BAIRSTOW 4

I.4. Flowchart

START

n ; tol

P = zeros(7, n+1)

i = n : -1 :0

P(1, n-i+1) = input ([‘konstanta untuk derajat ke –‘, num2str(i)’= ’])

r;s

disp
Disp p(1,:)
p(1,:) ;; nn

IteIte==00;;dr=1 ds=1; ;k=1

dr=1 ; ds=1 k=1

N ¿2

|dr|+|ds|>tol

A
B C

LABORATORIUM SIMULASI DAN KOMPUTASI

DEPARTEMEN TEKNIK KIMIA
INSTITUT TEKNOLOGI SEPULUH NOPEMBER
2019
METODE BAIRSTOW 5

B C A

P(4,1) = P(1,1)

i = 2 ; n+1

P(2,i) = r*p(4,i-1)

tidak i>2

ya

P(3,i) = s*p(4,i-2)

P(4,i) = p(1,i) + p(2,i) +

p(3,i)

p(7,1) = p(4,1)

i = 2 ; n+1

P(5,i) = r*p(7,i-1)

i>2 tidak

ya

P(6,i) = s*(p7,i-2)

B C D
E

LABORATORIUM SIMULASI DAN KOMPUTASI

DEPARTEMEN TEKNIK KIMIA
INSTITUT TEKNOLOGI SEPULUH NOPEMBER
2019
METODE BAIRSTOW 6

B C D E

p(7,i) =p(4,i)+p(5,i)+p(6,i)

DENOM = p(7,n-1)2-p(7,n)*p(7,n-2)

DENOM ≈ 0 tidak

ya

dr = (-p(4,n)*p(7,n-1))+(p(4,n+1)*p(7,n-2))/DENOM r = r + dr
ds = (-p(4,n+1)*p(7,n-1))+(p(4,n)*p(7,n))/DENOM
s = s + ds
r = r + dr
s = s + ds
ite = ite +1

root=ones(1,3)
root(1,2)= -r
root(1,3)=-s
D=root(1,2)^2-(4*root(1,1)*root(1,3))
x1=(-root(1,2)-D^0.5)/(2*root(1,1))
x2=(-root(1,2)+D^0.5)/(2*root(1,1))

disp k = x1 ; k+1 = x2, ite

A=zeros(1,n+1);
A(1,:)=p(1,:);
pnew = deconv(A,root)
p=zeros(7,n-1); n=n-2;
ite=0; k=k+2
dr=1; ds=1

LABORATORIUM SIMULASI DAN KOMPUTASI

DEPARTEMEN TEKNIK KIMIA
INSTITUT TEKNOLOGI SEPULUH NOPEMBER
2019
METODE BAIRSTOW 7

yeF

n == 2
tidak

ya

root=ones(1,3)
root(1,:) = p(1,:)
D=root(1,2)^2-(4*root(1,1)*root(1,3))
x1=(-root(1,2)-D^0.5)/(2*root(1,1))
x2=(-root(1,2)-D^0.5)/(2*root(1,1))

disp k=x1, k+1=x2

n == 1

tidak
ya

disp k = -p(1,2)/p(1,1)

end

LABORATORIUM SIMULASI DAN KOMPUTASI

DEPARTEMEN TEKNIK KIMIA
INSTITUT TEKNOLOGI SEPULUH NOPEMBER
2019
METODE BAIRSTOW 8

TUGAS KHUSUS

X3 – 20LX2 + A2SX + FA5.103 = 0

Dimana:
L = 167
A=3
S = 266939
F = -257

X3 – 20(167)X2 + (32)(266939)X + (-257)(35)(103) = 0

X3 – 3340X2 + 2402451X + 62451000 = 0
Setelah proses iterasi, didapat R= 1027 dan S= -27000
Faktor kuadratis: X2 – RX – S = X2 -1027X +27000
Polinomial hasil bagi: X – 2313
Maka:
(X3 – 3340X2 + 2402451X + 62451000) = (X2 – 1027X +27000) (X – 2313) = 0
(X – 1000) (X – 27) (X – 2313) = 0
X1 = 27 ; X2 = 1000 ; X3 = 2313

Metode Tol R S Error Ite X1 X2 X3

114558,4
-1019.7 -31391.7 4 31.775 987.924 2320.3

Bairstow 0.001 1026.9 -26997 4401.148 5 27.187 998.975 2314

1027 -27000 0.0009 6 27 1000 2313

LABORATORIUM SIMULASI DAN KOMPUTASI

DEPARTEMEN TEKNIK KIMIA
INSTITUT TEKNOLOGI SEPULUH NOPEMBER
2019