[PDF] [PDF] Méthodes dintégration numérique

Figure 5 3 – Méthode composite de Simpson (p = 2) pour m = 3 intervalles (c'est à dire n = 6 sous- intervalles et n +1=7 points au total) À nouveau, l'erreur est 



Previous PDF Next PDF





[PDF] Analyse Numérique

Corrigé du TD 4 EXERCICE 1 Estimation de l'erreur d'intégration élémentaire La fonction q La convergence de la méthode des trap`ezes composée est



[PDF] Série dexercices no2/5 Calculs approchés dintégrales - Université

(b) (1 point) Rappeler la formule de quadrature par la méthode de Simpson pour une soit une formule d'intégration numérique exacte pour les polynômes de 



[PDF] CORRIGÉS DES TRAVAUX DIRIGÉS DE lUE - Jerome BASTIEN

Intégration 21 Exercices Ce polycopié constitue les corrigés de TD de Méthodes Numériques de Base du département Mécanique (1) L' approximation de l'intégrale par la méthode des trapèzes composite (voir tableau A 4) est IT N = h



[PDF] Méthodes dintégration numérique

Figure 5 3 – Méthode composite de Simpson (p = 2) pour m = 3 intervalles (c'est à dire n = 6 sous- intervalles et n +1=7 points au total) À nouveau, l'erreur est 



[PDF] Analyse numérique Exercices corrigés

Intégration numérique Exercice 6 Déterminer par la méthode des trap`ezes puis par celle de Simpson ∫ π 2 0 f(x)dx sur la base du tableau suivant : x 0 π 8



[PDF] Fichier pdf

Intégration numérique Le but de ce Proposition 1 - La valeur approchée de l' intégrale de f sur I par la méthode des rectangles MÉTHODE DES TRAPÈZES



[PDF] Correction - Feuille de TD 2 : Méthodes dintégration numérique

On définit, pour toute fonction f continue sur [−1,1], la méthode d'intégration numérique T de la façon suivante : T(f) = λ1f(x1) + λ2f(x2) 1 Montrer que T est exacte 



[PDF] Exercices corrigés

Les formules de Newton-Cotes ont un degré de précision n + 1 si n est pair Exercice 23 (Formule composite des trapèzes) La méthode pour obtenir une formule 



[PDF] Méthodes numériques Intégration L3 2017/2018 Exercice 1 On veut

cos(x)dx = sin(2) par la méthode des trap`ezes et de Simpson Calculer l'erreur EN pour N ∈ 12, 22, , 210l pour les méthodes des trap`ezes et de Simpson et



[PDF] الجمهوريــــــــــــــــــــــة الجزائريـــــــــــة الديمقراطيـ - DSpace - USTO

3 Dérivation et intégration numérique 33 3 3 2 Méthode des Trapèzes Ce document notes de cours d'analyse numérique avec exercices corrigés re-

[PDF] intégration numérique simpson

[PDF] intégration par changement de variable exercices corrigés

[PDF] intégration par parties exercices corrigés

[PDF] intégrer antidote dans word mac

[PDF] intégrer dauphine en l3

[PDF] integrer ecole ingenieur apres licence

[PDF] intégrer une école d'ingénieur après la fac

[PDF] inter caisse

[PDF] interaction entre l'homme et l'environnement

[PDF] interaction gravitation exercices corrigés 3ème

[PDF] interaction microbienne pdf

[PDF] interaction onde matière bac math

[PDF] interactiunea corpurilor wikipedia

[PDF] interactiunea dintre corpuri

[PDF] interclub badminton departemental

Chapitre5

Méthodesd'intégrationnumérique

Lebut

Lebu tdecechapitre estd'abord erlecalc ulgénéraldel'intégraled 'unefonctionf(x)surundomai ne

finidéli mitépardesbornesfiniesaetb(lescasd esbornes infiniesn'est doncpascouvertic i): I= b a f(x)dx.(5.1)

Lesm otivations

Danscertain scastrèslimités,unetell eintégral epeutêtrecal culéeanalytiquement(àla main).Cep en-

dant,cen'estquetrès rar ementpossibl e,etlepl ussouventundes cassuivantsseprésente: longuesàévaluer -Cet teintégralen 'apasd'expressionanalytiqu e(parexem plelafon ctionerreur:Erf(x)= 2 x 0 e "x !2 dx Danstousce scas,on préfèreracalcu lernumériqu emen tlavaleurdel'intégraleI.

Lepri ncipe

L'idéeprincipalee stdetrouverdesmétho desquipermette ntdecalc ulerrapidementune valeurapproc hée

Idel' intégraleàcalculer:

I!I(5.2)

Commetoujour s,unprogrammenumériquen'in venter ien,etnefaitqueprocédertrè srapidementàuncalcul

quel'onpourrait enp rincipefaireàlamain .Uneméth odebienconnueconsis teparexem ple àdivis erl'aire

sousla courbeenu ngrandnombredepet its rectanglesd'ai re I k etdel essommer .Lerésult at I= k I k

estalor suneapproximati ondel 'intégraleI.Ce tteapproxima tionestd'autantmeilleurequelalargeurhdes

rectanglestendvers0,c'es tàdire:lim h$0 I=I.Ce tteméthodedi tedesrect anglesestunex emple parmi

d'autres.Nouslereverr ons,maisnousver ronsaussid' autresméthodes,plusgénéralesetpl usperformantes .

Pourpresqu etouteslesméthodes(s auflaméthodedeMonte -Carlo), l'intégralenumériqueestca lculéeà

partirdel'évaluatio ndel afonctionf(x)enunnom brede pointn+1distincts:f k =f(x k ),k"[0,n].Elle s'écritalors:

I=(b#a)

n k=0 w k f k (5.3) 33
UniversitéPaulSabatier2014-2 015CHAPITRE5.MÉ THODESD'INTÉGRATIONNUMÉRIQUE

Danscecas ,onparled eméthodesde quadrature.

Nousallonsvoir 4typesd eméthode sdi

érentes:

1.1-Lesmétho desdeNewto n-Cotessimples

1.2-Lesméth odesdeNewt on-Cotescomposites

1.3-Lesméth odes deGauss-Legendre

1.4-Lesmétho desdeMont e-Carlo

Performances

Lape rformanced'uneméthodesejuge encomparant

•lapré cisiondurésultat:Celle-cisecaractérise enesti mantl'erreur!entrel'approximati onetlavaleur

réelledel'intégrale : !=I#

I(5.4)

Lav aleurdel'erreurnep eutpasê trecalculéeexactemen tpuisqu'e ngénéral,onneconn aîtpasl'inté-

graleIquel'onche rch eàcalculer.Cependant,une majoratio npeutsouventêtree stiméeenétudiant

ledé veloppementensériedeTaylordelafonctionf(x).

•Larapid itéd'exécutionnécessairepouratteindr ecerésultat.D emanièregénérale,toutesles méthode s

peuventatteindredetrè sgrandesprécisions.Cependan t,let empsdecalcul augmenteaveclapréci- sion.Cetempsn'a ugm entepasdelamêm emanièrepourtouteslesméthodes sibienque certai nes s'avèrentpluse cacesque d'autres.Enpar ticulier,letempsde cal culdesméthodesdequadratureest proportionnelaunombredepointsoùlafo nctionf(x)estév aluée.

5.1LoisdeN ewton-Cotessi mples

Commenousall onslevoir,le sméthodesdeNewton- Cotessimplesneper mettentpas,àelles-seules, d'atteindredesprécisionssu santessurdesinterv all es[a,b]finisetneson tdonc jama isutiliséesda nsce

cas.Enr evanc he,ellesdeviennentprécisesl orsque|b#a|$0,etellesconstituentalorslabaseélémentaire

desmétho descompositesprésentéesdanslasectio nsuivant e.

5.1.1Princ ipe

Leprin cipegénéraledesméthodes deNewton-Cotessimplesestd 'approxim erlafonctionf(x)àinté-

grerpar unpolynômeP(x)!f(x).Sicetteapproximationestsu"sammentbonnealors,l'intégraledece polynôme I= b a

P(x)dx(5.5)

seraunebonneappr oxi mationdeI= b a valeurexactede I.Danscesméthodes,onchoisitdespolynômesdedegrépquicoïnc identavecf(x)enp+1

pointsdistincts ,espacésrégulièremententreles bornesaetb.Ces point ssontsituésauxposi tions:

{x k =a+kh,k"[0,p]}avech= b#a p (5.6)

Onaalors %k"[0,p]P(x

k )=f k =f(x k

Despolynôme sdedegrésdi

érentsdéfinisse ntdesméthodesdi

érentesauxperformances di

érentes.

Nousallonsvoir lesplusc ourantes,c 'estàd irelesméthodes d'ordreslesplusbas. 34
CHAPITRE5.MÉTHODES D'INTÉGR ATIONNUMÉRIQUEUniversit éPaulSabatier2014-2015

5.1.2Méth odedurectangle(p=0)

Cetteméthodeu tiliselepolynômedede gréleplusbas,àsavoir lepolyn ômeconstant: P 0 (x)=f(a)=f 0 .(5.7)

L'intégraleapproché e

I 0 b a P 0 (x)dxsecalc ulealorstrivial e- mentetdonne: I 0 =(b#a)f 0 (5.8)

Ils' agitdel'airedurec tangle .

Cetteintégral enumériquenécessiteuneuniq ueévaluationdelafonctionf(enx 0 =a)etreprésentedonc cequ' onpeutfairedepl usrapide.

L'erreurpeutêtrees timéeenutilisantlesd éveloppe mentsensériedeTa yloroulethéorè medesaccrois-

sementsfinisontrouv ealorspour h=b#a: &""[a,b]! 0 h 2 2 f (")c.a.d.|! 0 h 2 2 Sup [a,b] (|f |)(5.9)

Démonstration:Pourcalcul erl'erreur,onpeututilis erlethéorèmedesaccroisseme ntsfini s:!x"[a,b],#!"[a,b]

telque : f(x)=f(a)+(x$a)f Enremplaça ntdansl'expressiondel'intégrale etdel'erreur,ontrouve: "=I$ I b a (f(x)$P0(x))dx= b a (f(x)$f(a))dx b a (x$a)f (!)dx=f b"a 0 xdx (b$a) 2 2 f h 2 2 f L'erreur!n'estpasco nnuecarla valeurde""[a,b]resteindétermin ée.Cependant,onpeutlamajorer

parlaplusg randeval eurde ladérivéesurledoma ineconsidéré.Quelquesremar quessur cetteerreur :

-Cet teméthoded'i ntégrationestexactep ourtouteslesfonctionsfconstantes(danscecas! 0 =0 puisquequ'ellesv érifientf =0).Dans lecasplu sgén éralcettemé thodeestd'autant pluspréciseque lesvariationsd efsontfaibles( f petit). -Plusledomaine[a,b]estpetit ,plusl'erreurestfa ible.Cetteer reurdécroitenh 2 35
UniversitéPaulSabatier2014-2 015CHAPITRE5.MÉ THODESD'INTÉGRATIONNUMÉRIQUE

5.1.3Méth odedupointmilieu(p=0)

Cetteméthode utiliseégalementlepol ynômeconstantpour approximerlafonctionf.Cep endant,elleexploitemieuxle s symétriesduproblèmeenchoisissant lav aleurmilieu: P 0 !(x)=f a+b 2 =f 0 .(5.10)

L'intégraleapproch ée

I 0 b a P 0 (x)dxsecalcul ealorstrivi a- lementetdonne: I 0 !=(b#a)f 0 (5.11) Ils' agitdel'airedurect angle .Cetteméthodenécess iteun euniqueévaluationdelafonction f(enx 0

L'erreurpeutêtrees timéeenutilisantlesd éveloppe mentsensériedeT aylor,ouleth éorèmedesaccrois-

sementsfinis.Ontr ouvealorsp ourh=b#a: &""[a,b]! 0 h 3 24
f (")c.a.d.|! 0 h 3 24
Sup [a,b] (|f |)(5.12)

Démonstration:Pourcalcul erl'erreur,onpeututilis erlethéorèmedesaccroisseme ntsfini saudeuxième ordre:

!x"[a,b],#!"[a,b]telque: f(x)=f a+b 2 x$ a+b 2 f a+b 2 x$ a+b 2 2 f 2 Enremplaça ntdansl'expressiondel'intégrale etdel'erreur,ontrouve: "=I$ I b a (f(x)$P0(x))dx= b a f(x)$f a+b 2 dx b a x$ a+b 2 f a+b 2 x$ a+b 2 2 f 2 dx =f a+b 2 $"b"a 2 b"a 2 xdx+ f 2 "b"a 2 b"a 2 x 2 dx=0+ f 3 b$a 2 3 h 3 3 f L'erreur!n'estpasconn uecarlav aleurde""[a,b]resteindétermin ée.Cependant,onpeutlamajorer

parlaplus grandeva leurdeladér ivéesecondesurledomai neconsidéré.Quelquesr emarq uessurcette

erreur: pourlesfo nctionsa nes(dansc ecas! 0 !=0puisqu'ellesvérifientf =0). (f petit). -Plusledomaine[a,b]estpetit ,plusl'erreurestfa ible.Cetteer reurdécroitenh 3 ,c'estàdireplusvite quel'erreu rdelaméthodepréc édente: ! 0 0 petits,laméthodedupo intmil ieuesttoujourspluspr éciseq uelaméthode précédente. 36
CHAPITRE5.MÉTHODES D'INTÉGR ATIONNUMÉRIQUEUniversit éPaulSabatier2014-2015

5.1.4Méth odedutrapèze(p=1)

Pourapprox imerlafonctionf,cetteméthodeutiliselepolynôme d'ordre1(ladroite)q uipasse parf 0 =f(a)etf 1 =f(b): P 1 (x)= f 0 +f 1 2 f 1 #f 0 b#a x# a+b 2 (5.13)

L'intégraleapproché e

I 1 b a P 1 (x)dxsecalcul ealorsmathéma- tiquementougéométriquemente tdon ne: I 1 =(b#a) fquotesdbs_dbs1.pdfusesText_1