*   >> lukeminen koulutus Artikkelit >> education >> college and university

Simpsons Rule ja trapetsoidisäännöllä Numerical Integration - Matlab Scripts

Simpsonin sääntö ja trapetsoidisäännöllä Numerical Integration laskee viljelyala toiminto rikkomalla toiminto moniksi pienempiä alueita, jotka on helpompi laskea.

puolisuunnikkaan Sääntö

>> 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.

Simpsonin sääntö

>> 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

trapetsoidisäännöllä Matlab-Script

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 ;

Simpsonin sääntö Matlab Script

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

Muut Matlab-skriptit Luomani Sisällytä:

>>> Romberg numeerinen integrointi - Matlab Script

>>> Lagrangen menetelmä

Page   <<       [1] [2] >>

Copyright © 2008 - 2016 lukeminen koulutus Artikkelit,https://koulutus.nmjjxx.com All rights reserved.