Calculul Integralei definite

Metoda DreptunghiurilorMetoda Trapezelor

Integrala definită. Interpretarea geometrică

Fie dată funcţia f(x) continuă pe segmentul [a,b]. Integrala definită

( )b

a

f x dxeste aria trapezului curbiliniu, determinat de axa 0X, dreptele x = a şi x = b, şi graficul funcţiei f(x) pe segmentul [a,b]

© Sergiu Corlat, 2004

Metoda dreptunghiurilor

Determinarea ariei unei figuri curbilinii este destul de dificilă, de aceea se utilizează procedura de aproximare a figurii iniţiale prin un set de figuri geometrice, ariile cărora se determină după formule standard.

© Sergiu Corlat, 2004

© Sergiu Corlat, 2004

11

0

1

0

;

, 0,..., .

( )

2

(2 1)2

i

b

a

ni i

i

n

i

b ah

nx a ih i n

f x dx

x xh f

hh f a i

Varianta dreptunghiurilor de mijloc

1

0

1

0

;

, 0,..., .

( )

i

b

a

n

ii

n

i

b ah

nx a ih i n

f x dx

h f x

h f a ih

Varianta dreptunghiurilor de stânga

Varianta dreptunghiurilor de dreapta

1

1

;

, 0,..., .

( )

i

b

a

n

ii

n

i

b ah

nx a ih i n

f x dx

h f x

h f a ih

Algoritmul general (număr fix de divizări)

1. Se introduc limitele de integrare a,b şi numărul de divizări n.

2. Se calculează pasul de deplasare h

3. Pornind de la a se calculează mijlocul fiecărui segment elementar zi f(zi), şi ariile

dreptunghiurilor elementare.

4. Se sumează ariile elementare.

© Sergiu Corlat, 2004

Estimarea erorii y

f(x)g(x)= f((xi,+ xi+1)/2)

0 xi xi+1 …

f(x) se aproximează prin g(x).

1

( ) ( )

2i i

f x g x

x xM x

Eroarea la integrare pe un segment elementar este integrala erorii de aproximare:

1 1

1

2

1

[ ]

( )2

sup ( )

i i

i i

i i

x x

i i

x x

x x

x xf x dx g x dx M

M f x

© Sergiu Corlat, 2004

Algoritmul general (pentru o eroare fixată)

1 1

1

2

1

[ ]( ) , sup ( )

2

i i

i ii i

x x

i i

x xx x

x xf x dx g x dx M M f x

Eroarea de calcul al integralei pe un segment elementar nu depăşeşte

Prin urmare eroarea de calcul al integralei pe [a,b] nu depăşeşte sume erorilor pe segmentele elementare

2

[ , ]( ) ( ) , sup ( )

2 4

b

a ba

h hf x dx S nM b a M M f x

Pentru o eroare fixată numărul de divizări se calculează apriori:

22 ( )

( ) ( ) 4 1;4 4

h b a Mb a M b a M n n

© Sergiu Corlat, 2004

Exemplu program:var a,b,h,S :real; j,k,i,n :integer; function f(x:real):real; begin f:=5-(x*x-sin(5*x)); end;begin for j:=1 to 3 do begin a:=-2; b:=2; n:=0; for k:=1 to 10 do begin n:=n+10; S:=0; h:=(b-a)/n; for i:=0 to n-1 do case j of 1: s:=s+h*f(a+i*h); 2: s:=s+h*f(a+i*h+h); 3: s:=s+h*f(a+i*h+h/2); end; writeln('n=',n:3,' I=',s:0:6); end; end;end.

Rezultate:

Dreptunghiuri stanga:

n= 10 I=14.777608

n= 20 I=14.748804

n= 30 I=14.727351

n= 40 I=14.714402

n= 50 I=14.705922

n= 60 I=14.699972

n= 70 I=14.695577

n= 80 I=14.692201

n= 90 I=14.689529

n=100 I=14.687361

Dreptunghiuri dreapta:

n= 10 I=14.342392

n= 20 I=14.531196

n= 30 I=14.582279

n= 40 I=14.605598

n= 50 I=14.618878

n= 60 I=14.627436

n= 70 I=14.633403

n= 80 I=14.637799

n= 90 I=14.641171

n=100 I=14.643839

Dreptunghiuri mijloc:

n= 10 I=14.720000

n= 20 I=14.680000

n= 30 I=14.672593

n= 40 I=14.670000

n= 50 I=14.668800

n= 60 I=14.668148

n= 70 I=14.667755

n= 80 I=14.667500

n= 90 I=14.667325

n=100 I=14.667278

Metoda trapezelor

Aproximarea ariei unui trapez curbiliniu este mult mai eficientă în cazul cînd pe fiecare din segmentele elementare este aproximată prin un trapez, şi nu prin dreptunghi.

Pe segmentul elementar [xi, xi+1] trapezul este determinat de extremităţile segmentului pe axa 0X (xi,0) (xi+1 0) şi de valoarea funcţiei f(x) în extremităţi: (xi,f(xi)) (xi+1 ,f(xi+1))

© Sergiu Corlat, 2004

Aparatul matematic

;

,

0,..., .i

b ah

nx a ih

i n

1 11 1

0 0

( ) ( ) ( ) ( )( )

2 2

b n ni i i i

i ia

f x f x f x f xf x dx h h

Estimarea eroriiFie f(x) - de două ori derivabilă pe intervalul [a,b].

Pe un interval elementar [xi, xi+1] g(x) aproximează funcţia f(x) şi coincide cu ea în extremităţi. Eroarea aproximării este determinată de formula:

1

1

[ ]( ) ( ) sup ( )

2 i i

i i

x x

x x x xf x g x M M f x

© Sergiu Corlat, 2004

Eroarea la integrarea pe segmentul elementar este integrala erorii de aproximare:

1 1 1

1

1

31

[ ]

( )2

( ) sup ( )12

i i i

i i i

i i

x x xi i

x x x

i ix x

x x x xf x dx g x dx M dx

Mx x M f x

Eroarea la integrarea pe segmentul [a,b] este suma erorilor de integrare pe segmentele elementare:

23

[ , ]

( )( ) sup ( )

12 12

b

a ba

M b a Mhf x dx S n h M f x

© Sergiu Corlat, 2004

Algoritmul general (număr fix de divizări)

1. Se introduc limitele de integrare a,b şi numărul de divizări n.

2. Se calculează pasul de deplasare h

3. Pornind de la a se calculează valoarea funcţiei în

extremităţile fiecărui segment elementar şi ariile trapezelor elementare.

4. Se sumează ariile calculate.© Sergiu Corlat, 2004

Exemplu programvar a,b,h,S :real; j,k,i,n :integer;

function f(x:real):real; begin f:=5-(x*x-sin(5*x)); end;begina:=-2; b:=2; n:=0; for k:=1 to 10 do begin n:=n+10; s:=0; h:=(b-a)/n; for i:=0 to n-1 do s:=s+h*(f(a+i*h)+f(a+(i+1)*h))/2; writeln('n=',n:3,' I=',s:0:6); end;end.

© Sergiu Corlat, 2004

Rezultate:

Trapeze

n=100 I=14.665600n=200 I=14.666400n=300 I=14.666548n=400 I=14.666600n=500 I=14.666624n=600 I=14.666637n=700 I=14.666645n=800 I=14.666650n=900 I=14.666653n=1000 I=14.666656

Dreptunghiuri de stanga:n=100 I=14.687361n=200 I=14.677280n=300 I=14.673802n=400 I=14.672040n=500 I=14.670976n=600 I=14.670264n=700 I=14.669754n=800 I=14.669370n=900 I=14.669071n=1000 I=14.668832

Dreptunghiuri de dreapta:

n=100 I=14.643839n=200 I=14.655520n=300 I=14.659295n=400 I=14.661160n=500 I=14.662272n=600 I=14.663010n=700 I=14.663536n=800 I=14.663930n=900 I=14.664236n=1000 I=14.664480