ANALISIS NUMERIK

Solusi Persamaan Diferensial Biasa: Multi Step

Rahmat Nursalim, S.Pd., M.Si.

December 5, 2019

Rahmat Nursalim, S.Pd., M.Si. ANALISIS NUMERIK December 5, 2019 1 / 20

 Single Step & Multi Step
Single Step
Single Step method adalah metode untuk menentukan nilai suatu fungsi
dengan hanya menggunakan satu titik sebelumnya. Sebagai contoh, nilai
yk+1 ditentukan dengan hanya menggunakan yk

Metode Single Step

Metode Euler, Heun, Taylor, dan Runge-Kutta merupakan single step
methods

Rahmat Nursalim, S.Pd., M.Si. ANALISIS NUMERIK December 5, 2019 2 / 20

 Single Step & Multi Step
Multi Step
Multi Step method adalah metode untuk menentukan nilai suatu fungsi
dengan menggunakan beberapa titik sebelumnya. Sebagai contoh, nilai
yk+1 ditentukan dengan mengginakan yk−3 , yk−2 , yk−1 dan yk

Metode Multi Step

Metode Adams-Bahforth-Moulton, Milne-Simpson, dan Hamming
merupakan multi step methods

Rahmat Nursalim, S.Pd., M.Si. ANALISIS NUMERIK December 5, 2019 3 / 20

 Adams-Bashforth-Moulton Method
Prinsip
Adams-Bashforth-Moulton methods diturunkan dari teori utama kalkulus
Z tk+1
y (tk+1 ) = y (tk ) + f (t, y (t))dt
tk

Metode ini membutuhkan 4 titik sebelumnya. yk+1 ditentukan dengan

prediktor berdasarkan 4 titik dan korektor berdasarkan 3 titik dan 1
prediktor.

Rahmat Nursalim, S.Pd., M.Si. ANALISIS NUMERIK December 5, 2019 4 / 20

 Adams-Bashforth-Moulton Method
Prediktor-Korektor
Prediktor menggunakan polinomial lagrange untuk f (t, y (t)) berdasarkan
titik (tk−3 , fk−3 ), (tk−2 , fk−2 ), (tk−1 , fk−1 ) dan (tk , fk ) dengan rumus
h
pk+1 = yk + (−9fk−3 + 37fk−2 − 59fk−1 + 55fk )
24
Korektor dibangun dengan rumus yang sama polinomial lagrange untuk
f (t, y (t)) berdasarkan titik (tk−2 , fk−2 ), (tk−1 , fk−1 ), (tk , fk ) dan
(tk+1 , fk+1 ) = (tk+1 , f (tk+1 , pk+1 )) dengan rumus

h
yk+1 = yk + (fk−2 − 5fk−1 + 19fk + 9fk+1 )
24

Rahmat Nursalim, S.Pd., M.Si. ANALISIS NUMERIK December 5, 2019 5 / 20

 Error Estimasi
Error untuk integrasi numerik yang digunakan untuk menentukan prediktor
dan korektor adalah dari order O(h5 ), yaitu
251 (5)
y (tk+1 ) − pk+1 = y (ck+1 )h5 prediktor
720
−19
y (tk+1 ) − yk+1 = y (5) (dk+1 )h5 korektor
720
Misalkan bahwa h kecil dan y (5) (t) dekat dengan konstanta, maka langkah
yang melibatkan turunan ke 4 dan ke 5 dapat dieliminasi, dan hasilnya
−19
y (tk+1 ) − yk+1 ≈ (yk+1 − pk+1 )
270

Rahmat Nursalim, S.Pd., M.Si. ANALISIS NUMERIK December 5, 2019 6 / 20

 Milne-Simpson Method
Prinsip
Milne-Simpson Methods diturunkan dari teori utama kalkulus
Z tk+1
y (tk+1 ) = y (tk−3 ) + f (t, y (t))dt
tk−3

Metode ini membutuhkan 4 titik sebelumnya. yk+1 ditentukan dengan

prediktor berdasarkan 4 titik dan korektor berdasarkan 2 titik dan 1
prediktor.

Rahmat Nursalim, S.Pd., M.Si. ANALISIS NUMERIK December 5, 2019 7 / 20

 Milne-Simpson Method
Prediktor-Korektor
Prediktor menggunakan polinomial lagrange untuk f (t, y (t)) berdasarkan
titik (tk−3 , fk−3 ), (tk−2 , fk−2 ), (tk−1 , fk−1 ) dan (tk , fk ) dengan rumus
4h
pk+1 = yk−3 + (2f −f + 2fk )
3 k−2 k−1
Korektor dibangun dengan rumus yang sama polinomial lagrange untuk
f (t, y (t)) berdasarkan titik (tk−1 , fk−1 ), (tk , fk ) dan titik baru
(tk+1 , fk+1 ) = (tk+1 , f (tk+1 , pk+1 )) dengan rumus

h
yk+1 = yk−1 + (fk−1 + 4fk + fk+1 )
3

Rahmat Nursalim, S.Pd., M.Si. ANALISIS NUMERIK December 5, 2019 8 / 20

 Error Estimasi
Error untuk integrasi numerik yang digunakan untuk menentukan prediktor
dan korektor adalah dari order O(h5 ), yaitu
28 (5)
y (tk+1 ) − pk+1 = y (ck+1 )h5 prediktor
90
−1
y (tk+1 ) − yk+1 = y (5) (dk+1 )h5 korektor
90
Misalkan bahwa h kecil dan y (5) (t) dekat dengan konstanta, maka langkah
yang melibatkan turunan ke 4 dan ke 5 dapat dieliminasi, dan hasilnya
28
y (tk+1 ) − pk+1 ≈ (y − pk+1 )
29 k+1

Rahmat Nursalim, S.Pd., M.Si. ANALISIS NUMERIK December 5, 2019 9 / 20

 Milne-Simpson Method
Prinsip
Rumus Error pada Metode Milne-Simpson memberikan error untuk
prediktor berdasarkan dua nilai yaitu pk+1 dan yk+1 dan tidak
menggunakan y (5) (t). Dari sini dapat diimprovisasi menjadi nilai
prediksi
Berdasarkan asumsi bahwa antara nilai prediksi dan koreksi setiap step
berubah perlahan, kita dapat mengganti pk dan yk untuk pk+1 dan
yk+1 sehingga memperoleh modikasi dengan rumus
yk − pk
mk+1 = pk+1 + 28
29

Rahmat Nursalim, S.Pd., M.Si. ANALISIS NUMERIK December 5, 2019 10 / 20

 Milne-Simpson Method
Modikasi Milne-Simpson
Modikasi nilai mk+1 digunakan di tempat pk+1 pada tahap koreksi,
persamaan menjadi
h
yk+1 = yk−1 + (fk−1 + 4fk + f (tk+1 , mk+1 ))
3

Rahmat Nursalim, S.Pd., M.Si. ANALISIS NUMERIK December 5, 2019 11 / 20

 Milne-Simpson Method
Langkah
Langkah metode modikasi Milne-Simpson adalah
4h
pk+1 = yk−3 + [2f −f + 2fk ]
3 k−2 k−1
y k − pk
mk+1 = pk+1 + 28
29
fk+1 = f (tk+1 , mk+1 )
h
yk+1 = yk−1 + (fk−1 + 4fk + fk+1 )
3

Rahmat Nursalim, S.Pd., M.Si. ANALISIS NUMERIK December 5, 2019 12 / 20

 Hamming Method
Prediktor-Korektor
Prediktor menggunakan polinomial lagrange untuk f (t, y (t)) berdasarkan
titik (tk−3 , fk−3 ), (tk−2 , fk−2 ), (tk−1 , fk−1 ) dan (tk , fk ) dengan rumus
4h
pk+1 = yk−3 + (2f −f + 2fk )
3 k−2 k−1
Korektor dibangun dengan rumus yang sama polinomial lagrange untuk
f (t, y (t)) berdasarkan titik (tk−2 , fk−2 ), (tk−1 , fk−1 ), (tk , fk ) dan
(tk+1 , fk+1 ) dengan rumus

−yk−2 + 9yk h
yk+1 = + (−fk−1 + 2fk + fk+1 )
8 3

Rahmat Nursalim, S.Pd., M.Si. ANALISIS NUMERIK December 5, 2019 13 / 20

 Contoh
Gunakan Metode Adams-Basforth-Moulton, Milne-Simpson, dan Hamming
untuk menyelesaikan

y 0 = (t − y )/2 pada interval [0, 3]

dengan y (0) = 1 dengan h = 18 . Runge-Kutta telah menentukan nilai

y1 = 0.94323919, y2 = 0.89749071, y3 = 0.86208736

Rahmat Nursalim, S.Pd., M.Si. ANALISIS NUMERIK December 5, 2019 14 / 20

 Contoh
Example
Untuk membandingkan presisi dari setiap step, maka diperlukan tabel
lengkap sebagai berikut.
Adams-
Bashfort- Milne- Hamming
k Moulton Error Simpson Error Hamming Error

0.00 1.00000000 0E-8 1.00000000 0E-8 1.00000000 0E-8

0.50 0.83640227 8E-8 0.83640231 4E-8 0.83640234 1E-8
0.625 0.81984673 16E-8 0.81984687 2E-8 0.81984688 1E-8
0.75 0.81186762 22E-8 0.81186778 6E-8 0.81186783 1E-8
1.00 0.81959166 28E-8 0.81959190 3E-8 0.81959198 0E-8
2.00 1.10363781 51E-8 1.10363822 10E-8 1.10363834 -1E-8
3.00 1.66938998 50E-8 1.66939038 10E-8 1.66939050 -2E-8

Rahmat Nursalim, S.Pd., M.Si. ANALISIS NUMERIK December 5, 2019 15 / 20

 Program Adams-Bashforth-Moulton Method
function A=abm(f,T,Y)
n=length(T);
if n<5, break, end;
F=zeros(1,4);
F=feval(f,T(1:4), Y(1:4));
h=T(2)-T(1);
for k=4:n-1
p=Y(k)+(h/24)*(F*[-9 37 -59 55]');
T(k+1)=T(1);
F=[F(2) F(3) F(4) feval(f,T(k+1),p)];
Y(k+1)=Y(k)+(h/24)*(F*[1 -5 19 9]');
F(4)=feval(f,T(k+1), Y(k+1));
end;
A=[T' Y'];
Rahmat Nursalim, S.Pd., M.Si. ANALISIS NUMERIK December 5, 2019 16 / 20
 Program Milne-Simpson Method
function M=milne(f,T,Y)
n=length(T);
if n<5, break, end;
F=zeros(1,4); F=feval(f,T(1:4), Y(1:4));
h=T(2)-T(1); pold =0; yold =0
for k=4:n-1
pnew=Y(k-3)+(4*h/3)*(F(2:4)*[2 -1 2]';
pmod=pnew + 28*(yold-pold)/29;
T(k+1) = T(1)+h*k;
F=[F(2) F(3) F(4) feval(f,T(k+1), pmod)];
Y(k+1)= Y(k-1) + (h/3)*(F(2:4)*[1 4 1]');
pold=pnew; yold=Y(k+1);
F(4)=feval(f,T(k+1),Y(k+1));
end;
M=[T' Y'];
Rahmat Nursalim, S.Pd., M.Si. ANALISIS NUMERIK December 5, 2019 17 / 20
 Program Hamming Method
function H=hamming(f,T,Y)
n=length(T);
if n<5, break, end;
F=zeros(1,4); F=feval(f,T(1:4), Y(1:4));
h=T(2)-T(1); pold =0; cold =0
for k=4:n-1
pnew=Y(k-3)+(4*h/3)*(F(2:4)*[2 -1 2]';
pmod=pnew + 112*(cold-pold)/121;
T(k+1) = T(1)+h*k;
F=[F(2) F(3) F(4) feval(f,T(k+1), pmod)];
cnew = (9*Y(k)-Y(k-2) + 3*h*(F(2:4)*[-1 2 1]'))/8;
Y(k+1)= cnew + 9*(pnew-cnew)/121;
pold=pnew; yold=Y(k+1);
F(4)=feval(f,T(k+1),Y(k+1));
end;
H=[T' Y'];
Rahmat Nursalim, S.Pd., M.Si. ANALISIS NUMERIK December 5, 2019 18 / 20
 Latihan
1 Diberikan IVP
y 0 = t2 − y dengan y (0) = 1
Misalkan h = 0.2 dan solusi yang diinginkan berada pada interval
[0, 1], gunakanlah metode Adams-Bashforth-Moulton untuk
menentukan yk
2 Diberikan IVP

y 0 = 3y − 3t dengan y (0) = 1

Misalkan h = 0.1 dan solusi yang diinginkan berada pada interval

[0, 1], gunakanlah metode Milne-Simpson untuk menentukan yk

Rahmat Nursalim, S.Pd., M.Si. ANALISIS NUMERIK December 5, 2019 19 / 20

 THANK YOU

Rahmat Nursalim, S.Pd., M.Si. ANALISIS NUMERIK December 5, 2019 20 / 20