Pertemuan 11 : Pengintegralan Numerik - Matematika IPB
Pertemuan 11 : Pengintegralan Numerik - Matematika IPB
Pertemuan 11 : Pengintegralan Numerik - Matematika IPB
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>