Praktikum 12 Integrasi Numerik Ruang lingkup bahasan Kuadratur
Praktikum 12 Integrasi Numerik
Ruang lingkup bahasan Ø Kuadratur Ø Aturan Trapesium Ø Aturan Simpson Ø Integrasi Romberg
Aturan Simpson • Aturan Simpson Majemuk
Contoh • Dengan menggunakan aturan Simpson, hampiri integral fungsi f(x)=1/x pada [1, 9] dengan delapan subinterval seragam.
Numerik function hasil=simpmaj(f, x 0, x 1, n) h=(x 1 -x 0)/n; hasil=[]; c=abs((x 1 -x 0)/n); x=x 0: c: x 1; a=f(x 0); b=f(x 1); c=0; d=0; for i=1: n/2, d=d+4*(f(2*i)); end for i=2: n/2, c=c+2*(f(2*i-1)); end hasil=(a+b+c+d)*h/3
Integrasi Romberg • Aturan rekursif trapesium • Misalkan f suatu fungsi yang terdefinisi pada [a, b] dan h=(b-a). Untuk n=1, 2, 4, 8, … atau n=20, 21, 22, 23, …, kita definisikan barisan aturan trapesium dengan Barisan aturan trapesium tersebut memenuhi hubungan
Integrasi Romberg • Dalam menghitung hampiran suatu integral dengan menggunakan aturan trapesium rekursif, kita lakukan langkah-langkah sebagai berikut.
Secara numerik • • • function Tn=trreku(f, n, a, b) h=b-a; if n==0, Tn=h*(f(a)+f(b))/2; else if n>0 index=[1: 2: 2^n-1]; x=a+h*index/(2^n); F=f(x); Jf=sum(F); Tn=trreku(f, n-1, a, b)/2+Jf*h/(2^n); end endfunction
Contoh • Hitunglah hampiran integral dengan menggunakan aturan trapesium rekursif dengan cacah interval 1, 2, 3, 4, dan 5. hitunglah pula galat masing-masing hampiran. •
Integrasi Romberg • Aturan rekursif simpson • Misalkan {Tn; n=0, 1, 2, …} adalah barisan aturan trapesium majemuk dan Sn adalah aturan simpson majemuk untuk fungsi f dengan 2 n subinterval pada interval [a, b]. Hubungan antara aturan simpson dangan aturan trapesium majemuk adalah • Aturan rekursif boole • Misalkan {Sn; n=0, 1, 2, …} adalah barisan aturan simpson majemuk dan Bn adalah aturan boole majemuk untuk fungsi f dengan 2 n subinterval pada interval [a, b]. Hubungan antara aturan boole dangan aturan simpson majemuk adalah
Integrasi Romberg • Metode Romberg Integrasi Romberg dengan ekstrapolasi Richardson
Integrasi Romberg
Contoh • Dengan menggunakan metode Romberg carilah hampiran integral tentu
Contoh • Secara analitik
Dalam Scilab deff('y=f(x)', 'y=(x. ^2+x+1). *cos(x)') T=[]; for n=0: 5, Tn=traperekursif(f, n, 0, %pi/2); T=[T; Tn]; end i=[2: 6]'; S(i)=(4*T(i, 1)-T(i-1, 1))/3; i=[3: 6]'; B(i)=(16*S(i)-S(i-1))/15; i=[4: 6]'; R 4(i)=(64*B(i)-B(i-1))/63; i=[5: 6]'; R 5(i)=(256*R 4(i)-R 4(i-1))/255; R 6(6)=(4^5*R 5(6)-R 5(5))/(4^5 -1); [T S B R 4 R 5 R 6]
- Slides: 15