Simpsonin sääntö ja trapetsoidisäännöllä Numerical Integration laskee viljelyala toiminto rikkomalla toiminto moniksi pienempiä alueita, jotka on helpompi laskea.
>> puolisuunnikkaan ('sin (x) /x', 1,3,8)
ans =
0,902337806742469
I = ( ((3-1) /8) /2) [f
(1) +2 {f (1,25) + f (1,5) + f (1,75) + f
(2) + f (2,25) + f (2,5) + f (2.
75)} + f
(3)]
= 0,125 [0,841471 + 2 (3,165097) 0,047040)
= 0,902338
trapetsoidisäännöllä varten numeerinen integrointi käyttäen trapazoids laskea ala-toiminto. Jokainen trapazoid malleja lineaarinen viiva kahden pisteen funktio ja laskee ala linja.
>> simprule ('sin (x) /x', 1,3 , 8)
ans =
0,902568788567005
Sen sijaan käyttää lineaarista linjat mallin toiminto, Simpsonin sääntö laskee alueella mallintamalla polynomi toiminnon. Tämä polynomi on tarkka muille polynomit aste 2 ja 3.
I = (((3-1) /8) /3) [f
(1) +4 {f (1,25) + f (1.75 ) + f (2,25) + f (2,75)} + 2 {f (1,5) + f
(2) + f (2,5)} + f
(3)]
= (0,25 /3) [0,841471 +4 (1.806062) +2 (1.359035) +0.047040)
= 0.902569
funktio I = puolisuunnikkaan (f_str, a, b, n) % puolisuunnikkaan puolisuunnikassääntöä yhdentymistä.% I = puolisuunnikkaan (F_STR, A, B, N) palauttaa puolisuunnikassääntöä lähentämisestä% varten kiinteä f (x) X = x = B, käyttäen N osaintervallien, jossa% F_STR on merkkijono edustus f.
= 0; g = inline (f_str), h = (ba) /n;
I = I + G ();
< p> ii = (+ h): h: (bh) I = I + 2 * g (ii); lopussa
I = I + G (b); I = I * h /2 ;
toiminto i = simprule (f_str, a, b, n)% SIMPRULE Simpsonin sääntö yhdentymistä.% I = SIMPRULE (F_STR, A, B, N) palaa Simpsonin sääntö lähentämisestä% varten kiinteä f (x) X = x = B, käyttäen N osaintervallien, jossa% F_STR on merkkijono esitys f.% virhe syntyy, jos N ei ole positiivinen, vaikka kokonaisluku .
= 0; g = inline (f_str), h = (ba) /n;
jos ((n> 0) && (rem (n, 2) == 0 )) I = I + G (); ja ii = (+ H): 2 * h: (bh) I = I + 4 * g (ii); end KK = (+ 2 * h): 2 * h: (b-2 * h) I = I + 2 * g (kk); end I = I + G (b);
I = I * h /3, muuten DISP ("väärä Vastinetta N ') päähän
>>> Romberg numeerinen integrointi - Matlab Script
>>> Lagrangen menetelmä