04.05.2013 Views

Pertemuan 11 : Pengintegralan Numerik - Matematika IPB

Pertemuan 11 : Pengintegralan Numerik - Matematika IPB

Pertemuan 11 : Pengintegralan Numerik - Matematika IPB

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 1<br />

(P<strong>11</strong>) <strong>Pengintegralan</strong><br />

numerik<br />

Rujukan:<br />

1.Buku 1 bab 5<br />

2.Buku 2 bab 7<br />

• Pengantar Kuadratur<br />

• Metode <strong>Pengintegralan</strong><br />

<strong>Numerik</strong> :<br />

• Aturan Trapesium<br />

(Trapezoid)<br />

• Aturan Trapesium<br />

Majemuk<br />

• Aturan Simpson<br />

• Aturan Simpson<br />

Majemuk<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 2<br />

Pengantar Kuadratur<br />

• Pendekatan analitik untuk menghitung integral<br />

tentu menggunakan TDK II<br />

dengan<br />

b<br />

<br />

a<br />

f ( x)<br />

dx F(<br />

b)<br />

F(<br />

a)<br />

F( x)<br />

f ( x)<br />

• Bila anti turunan dari suatu fungsi f(x), F(x),<br />

tidak dpt ditentukan, maka TDK II tidak dapat<br />

digunakan.<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 3<br />

• <strong>Pengintegralan</strong> numerik adalah alat utama yang<br />

digunakan oleh para ilmuwan untuk memperoleh<br />

hampiran jawaban bagi nilai integral tentu yang<br />

tidak dapat diselesaikan secara analitik.<br />

• Untuk beberapa kasus, meskipun cara analitik dapat<br />

digunakan, namun bentuknya sangat kompleks,<br />

sehingga pendekatan numerik dirasa lebih mudah<br />

dengan hasil yang cukup akurat.<br />

• Pendekatan numerik dilakukan dengan mencari<br />

fungsi hampiran bagi f(x), yaitu dengan polinomial<br />

Taylor atau interpolasi polinomial.<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 4<br />

• Contoh, untuk mengevaluasi<br />

e<br />

e<br />

dapat digunakan hampiran fungsi berikut :<br />

t<br />

x<br />

2<br />

1 2<br />

1<br />

t t <br />

<br />

2!<br />

2 1 4<br />

1<br />

x x <br />

<br />

2!<br />

dengan:<br />

0 <br />

dx x<br />

2<br />

1<br />

t<br />

n!<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong><br />

n<br />

1<br />

x<br />

n!<br />

<br />

2n<br />

1<br />

<br />

0<br />

1<br />

( n <br />

e x<br />

1)!<br />

t<br />

2<br />

dx<br />

n1<br />

1<br />

x<br />

( n 1)!<br />

e<br />

ct<br />

2n2<br />

e<br />

dx


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 5<br />

1<br />

<br />

0<br />

e<br />

x<br />

2<br />

dx<br />

<br />

1<br />

<br />

0<br />

( 1<br />

<br />

x<br />

2<br />

<br />

1 3 1 5 1 1<br />

x x x <br />

x<br />

3 10 n!<br />

2n<br />

1<br />

1 1 1 1<br />

1<br />

<br />

E<br />

3 10 n!<br />

2n<br />

1<br />

Dengan mengambil n=3, maka :<br />

1<br />

2!<br />

x<br />

4<br />

<br />

<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong><br />

1<br />

n!<br />

x<br />

2n<br />

1 1 1<br />

1<br />

E 1.<br />

4571<br />

E<br />

3 10 42<br />

1<br />

e 8 e<br />

0 E<br />

x dx 0.<br />

0126<br />

24 216<br />

0<br />

) dx <br />

E<br />

2n1<br />

<br />

1<br />

0<br />

<br />

E


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 6<br />

• Meskipun terlihat mudah, tapi pada umumnya<br />

pengintegralan menggunakan pendekatan deret<br />

Tayor tersebut tidak selalu mudah, sehingga cara<br />

yang kedua, yaitu menggunakan fungsi interpolasi<br />

polinomial lebih disenangi.<br />

• Integral tentu pada interval [a, b] didefinisikan sbg:<br />

b<br />

n<br />

f ( x)<br />

dx <br />

a<br />

n<br />

i1<br />

I(<br />

f ) lim<br />

<br />

• Aturan kuadratur n titik:<br />

f ( x ) x<br />

• Titik x i disebut titik kuadratur, w i disebut pembobot<br />

Q<br />

i<br />

n<br />

(<br />

f<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong><br />

i<br />

)<br />

<br />

n<br />

<br />

i1<br />

w<br />

i<br />

f<br />

( x<br />

i<br />

)


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 7<br />

Definisi<br />

• Andaikan bahwa a = x 0 < x 1 < ... < x M = b. Suatu<br />

formula dalam bentuk<br />

[<br />

M<br />

<br />

k0<br />

k k 0 0 1 1<br />

Q f ] w f ( x ) w f ( x ) w f ( x ) <br />

w<br />

b<br />

<br />

a<br />

f<br />

( x)<br />

dx Q[<br />

f ] E[<br />

f<br />

• disebut suatu pengintegralan numerik atau rumus<br />

kuadratur.<br />

• suku E[f] disebut kesalahan akibat pemangkasan<br />

selama proses interpolasi.<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong><br />

]<br />

M<br />

f<br />

( x<br />

M<br />

)


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 8<br />

M<br />

x k0<br />

• k disebut titik kuadratur<br />

• disebut faktor pembobot<br />

M<br />

wk k0<br />

• Tergantung pada aplikasinya, titik-titik {x k}<br />

dapat dipilih dalam berbagai cara.<br />

• Pada beberapa metode, seperti metode<br />

Trapesium dan Simpson, titik-titik yang dipilih<br />

berjarak sama.<br />

• Dalam hal ini juga penting untuk diketahui<br />

akurasi dari solusi numerik yang diperoleh.<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 9<br />

Teorema (Aturan Trapesium)<br />

• Pandang y = f(x) pada [x 0, x 1], dengan<br />

x 1= x 0 + h. Aturan Trapesium<br />

merupakan suatu hampiran numerik untuk f(x)<br />

pada [x 0, x 1]. Jadi :<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 10<br />

Sisaan dari aturan trapesium ini berbentuk<br />

dengan c terletak antara x 0 dan x 1, serta<br />

mempunyai persamaan<br />

Contoh:<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> <strong>11</strong><br />

Andaikan<br />

maka bila sekatan dapat dibuat lebih kecil,<br />

akan diperoleh:<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 12<br />

Aturan Trapesium Majemuk<br />

• Bila beberapa trapesium digunakan maka disebut<br />

aturan trapesium majemuk (composite<br />

trapezoidal). Ini dilakukan untuk mendapatkan<br />

tingkat akurasi yang lebih tinggi.<br />

• Semakin banyak jumlah sekatan yang dibuat,<br />

akan semakin akurat hasil yang diperoleh.<br />

• Contoh, bila dibuat 3 sekatan.<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 13<br />

Bandingkan kesalahan yang terjadi<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 14<br />

Teorema Trapesium Majemuk<br />

• Pandang y=f(x) pada [a,b]. Andaikan interval [a,<br />

b] dipartisi menjadi m subinterval<br />

berjarak sama yaitu dengan titik-titik<br />

x k= x 0 + kh untuk k = 1,2,…m.<br />

• Aturan Trapesium Majemuk untuk m subinterval<br />

merupakan suatu hampiran numerik untuk f(x)<br />

pada [a, b], sehingga dapat dituliskan sebagai<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 15<br />

Sisaan dari Aturan Trapesium Majemuk<br />

• Andaikan [a, b] dibagi menjadi m<br />

subinterval dengan jarak<br />

• Aturan Trapesium Majemuk<br />

adalah suatu hampiran numerik bagi integral<br />

dan<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 16<br />

• Lebih jauh, jika maka terdapat<br />

suatu nilai c dengan a < c < b sedemikian<br />

sehingga berbentuk :<br />

• Atau bila dituliskan dalam notasi "big O"<br />

menjadi<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 17<br />

• Pada contoh di atas:<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 18<br />

≈ T n(f)<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 19<br />

Contoh:<br />

• Tentukan hampiran numerik dari<br />

menggunakan aturan Trapesium m = 1, 4, 8, dan 16<br />

Jawab:<br />

m Tm(f)<br />

1 4.04864<br />

4 3.49710<br />

8 3.46928<br />

16 3.46232<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 20<br />

Contoh:<br />

• Tentukan hampiran numerik dari<br />

dgn aturan Trapesium m = 100, 200, 400, dan 800<br />

Jawab:<br />

m<br />

Tm(f)<br />

100 3.46005707746<br />

200 3.46001252350<br />

400 3.46000138500<br />

800 3.45999860038<br />

Solusi analitik menghasilkan nilai integral = 3.45999767217<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 21<br />

kesalahan dari hampiran numerik:<br />

m Kesalahan hampiran<br />

100 -0.0000594053<br />

200 -0.0000148513<br />

400 -0.00000371283<br />

800 -0.000000928209<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 22<br />

Aturan Simpson<br />

• Bila aturan trapesium menggunakan interpolasi<br />

polinomial derajat linear, maka aturan Simpson<br />

menggunakan interpolasi polinomial derajat 2 :<br />

Sehingga diperoleh :<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 23<br />

Gambarannya seperti berikut :<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 24<br />

Teorema Aturan Simpson<br />

• Pandang y=f(x) pada [x 0, x 2], dengan x 1= x 0 + h, dan<br />

x 2= x 0 + 2h. Aturan Simpson adalah<br />

• Sisaan dari aturan Simpson ini adalah<br />

dengan c terletak di antara x 0 dan x 2 serta<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 25<br />

Contoh:<br />

Sebagai pembanding, dgn aturan Trapesium:<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 26<br />

Aturan Simpson Majemuk<br />

• Bila beberapa parabolik digunakan, metode ini<br />

disebut aturan Simpson Majemuk (composite<br />

Simpson). Pendekatan ini digunakan untuk<br />

memperoleh tingkat akurasi yang lebih tinggi.<br />

• Untuk menghitung :<br />

digunakan cara berikut :<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 27<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 28<br />

yang dapat disederhanakan menjadi<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 29<br />

Teorema Aturan Simpson Majemuk<br />

• Pandang y = f(x) pada [a, b]. Andaikan [a, b] disekat<br />

menjadi 2m subinterval berjarak sama<br />

yaitu dengan titik-titik x k= x 0 + kh untuk<br />

k = 0,1,2,…2m, maka aturan Simpson Majemuk<br />

merupakan suatu hampiran numerik untuk f(x) pada<br />

[a, b], sehingga dapat dituliskan sebagai<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 30<br />

Sisaan dari Aturan Simpson Majemuk<br />

• Andaikan [a, b] disekat menjadi 2m<br />

subinterval berjarak sama yaitu<br />

. Aturan Simpson Majemuk<br />

merupakan suatu hampiran numerik untuk<br />

integral, dan<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 31<br />

• Lebih jauh, jika , maka terdapat<br />

suatu nilai c di antara a dan b sedemikian<br />

sehingga kesalahan berbentuk.<br />

Bila diekspresikan menggunakan notasi "big O"<br />

menjadi<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 32<br />

Contoh: ≈ Sn(f)<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 33<br />

Contoh:<br />

• Tentukan hampiran numerik dari<br />

menggunakan aturan Simpson m = 1, 4, 8, dan 16<br />

Jawab:<br />

m Sm(f)<br />

1 3.4613498419<br />

4 3.46008250981<br />

8 3.46000297964<br />

16 3.45999800397<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 34<br />

Contoh:<br />

• Tentukan hampiran numerik dari<br />

dgn aturan Trapesium m = 20, 40, 80, dan 160<br />

Jawab:<br />

m<br />

20<br />

40<br />

80<br />

160<br />

Solusi analitik menghasilkan nilai integral =<br />

Sm(f)<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 35<br />

kesalahan dari hampiran numerik:<br />

m Kesalahan hampiran<br />

20<br />

40<br />

80<br />

160<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>


<strong>Pertemuan</strong> <strong>11</strong> : <strong>Pengintegralan</strong> <strong>Numerik</strong> 36<br />

Semester Pendek 2012 Departemen <strong>Matematika</strong> FMIPA <strong>IPB</strong>

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!