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

View - Download Correction - Feuille de TD 2 : Méthodes d’intégration numérique

Previous PDF

Next PDF

Faculté des sciences et ingénierie (Toulouse III) Année universitaire

Département de mathématiques 2019-2020

L2 Maths, UE d"Analyse numérique

Correction - Feuille de TD 2 : Méthodes d"intégration numériqueExercice 1.(Une méthode sur[-1,1]) Soientx1,x2?[-1,1],x1< x2, etλ1,λ2?R.On définit, pour toute fonctionfcontinue sur[-1,1], la méthode d"intégration numériqueTde la façon suivante :

T(f) =λ1f(x1) +λ2f(x2).

1. Mon trerque Test exacte d"ordre au moins 1 sur[-1,1]si et seulement siλ1= 2x2x

2-x1etλ2=2x1x

1-x2. Correction :D"après le cours,Test exacte d"ordre au moins1si et seulement si elle est exacte pour les polynômes d"une base deR1[X], par exemple les polynômes de la base canonique :1etX,i.e

T(1) =?

1 -1dx= 2etT(X) =? 1 -1xdx= 0. Or,

T(1) =λ1+λ2etT(X) =λ1x1+λ2x2

ce qui conduit équivaut au système suivant :

1+λ2= 2

1x1+λ2x2= 0

Ce qui donne finalementλ1=2x2x

2-x1etλ2=2x1x

1-x2. 2. P ourquelles v aleursde λ1,λ2,x1etx2,Test-elle au moins exacte d"ordre 3? Quel est alors l"ordre d"exactitude de la méthode? Correction :De même que précédemment,Test au moins d"ordre3si et seule- ment si elle est exacte pour les polynômes de la base canonique deR3[X], soit, si et seulement si elle est exacte d"ordre1et pourX2etX3, ce qui, d"après la question précédente équivaut à

1+λ2= 2

1x1+λ2x2= 0

T(X2) =λ1x21+λ2x22=?1

-1x2dx=23

T(X3) =λ1x31+λ2x32=?1

-1x3dx= 0. 1

1=2x2x

2-x1, λ2=2x1x

1-x2λ1x21+λ2x22=23

1x31+λ2x32= 0.

1=2x2x

2-x1, λ2=2x1x

1-x2x21x2-x1x22=x2-x13

x

1x2(x1-x2)(x1+x2) = 0.

Or, on sait par hypothèse quex1< x2. De plus, les casx1= 0etx2= 0mènent à une contradiction avecx21x2-x1x22=x2-x13 . On a doncx2=-x1et, par un calcul simple, la2-ème équation donnex1=-⎷3 3 x2=⎷3 3 , ce qui conduit finalement à

1=λ2= 1.

La méthodeTest donc de degré d"exactitude au moins3, reste à savoir siT(X4) =?1 -1x4dx. Un calcul élémentaire donne :

T(X4) = 2⎷3

43
4=29 ?=25 1 -1x4dx ce qui finit de prouver que la méthode est exacte d"ordre 3. 3. Déduire des qu estionsprécéden tesune métho ded"in tégrationd"ordre 3 sur un segment[a,b]quelconque. Correction :Pour cette question, on utilise le changement de variable affine vu en cours : ?b af(x)dx=b-a2 1 -1f?b-a2 t+a+b2 dt.

On approche alors

?b af(x)dxparb-a2 T? x?→f(b-a2 x+a+b2 , c"est-à-dire T a,b(f) =b-a2 f(-(b-a)⎷3 6 +a+b2 ) +f((b-a)⎷3 6 +a+b2 Reste à montrer que cette méthode est exacte pour les polynômes de degré63.

SoitPun tel polynôme. Il est clair que

Q(X) =b-a2

P?b-a2

X+a+b2

est aussi un polynôme de degré inférieur ou égal à3. CommeTest une méthode de degré d"exactitude3, alors

T(Q) =?

1 -1Q(x)dx=b-a2 1 -1P?b-a2 x+a+b2 dx=? b aP(x)dx. Or, par construction,T(Q) =Ta,b(P), ce qu"il fallait démontrer. 2 Exercice 2.(De l"interpolation à l"intégration numérique) Soitf: [-1,1]-→Rune fonction de classeC2et soientx0,x1?[-1,1]avecx0?=x1. 1.

In terpolationde Lagrange aux no eudsx0, x1:

a. Donner l"expression du p olynômeP1d"interpolation de Lagrange defassocié aux noeudsx0,x1dans la base de Lagrange. Correction :Comme on l"a vu au chapitre pércédent, le polynôme d"interpo- lation de Lagrange associé à ces noeuds dans la base de Lagrange est : P

1(X) =f(x0)X-x1x

0-x1+f(x1)X-x0x

1-x0. b. Donner la form uled"erre urd"appro ximationde Lagrange sup x?[-1,1]|f(x)-P1(x)| en fonction deM= sup z?[-1,1]|f??(z)|. Correction :Le résultat principal du chapitre suivant donne ?x?[-1,1],|f(x)-p1(x)|6M2 |x-x0||x-x1|. 2. On considère la métho ded"in tégrationn umériques ur[-1,1]suivante pour appro- cherI(f) =?1 -1f(t)dt:

J(f) =?

1 -1P1(x)dx. a. Mon terqu"il existe λ0,λ1?Rtels queJ(f) =λ0f(x0) +λ1f(x1).

Correction :On se sert de la question 1)a) :

J(f) =?

1 -1p1(x)dx=? 1 -1? f(x0)x-x1x

0-x1+f(x1)x-x0x

1-x0? dxJ(f) =f(x0)? 1 -1x-x1x

0-x1dx+f(x1)?

1 -1x-x0x

1-x0dx

On a doncλ0=?

1 -1x-x1x

0-x1dxetλ1=?

1 -1x-x0x

1-x0dx

b. Donner une ma jorationde l"erreur |I(f)-J(f)|en fonction deM.

Correction :Par définition deJ(f)on a :

|I(f)-J(f)|=????? 1 -1f(x)dx-? 1 -1p1(x)dx????=????? 1 -1(f(x)-p1(x))dx????6? 1 -1|f(x)-p1(x)|dx En utilisant la majoration obtenue à la question 1-b) on a alors : |I(f)-J(f)|6? 1 -1M2 |x-x0||x-x1|dx

Finalement, si on majore|x-xi|par 2, on obtient

|I(f)-J(f)|64M 3

Exercice 3.(Exam 2016)

Soitf: [0,1]?→Rune application de classeC1.

1. A l"aide d"un dév eloppementde T aylorde la fonction

F(x) =?

x

0f(t)dt

montrer qu"il existec?]0,1[tel que?1

0f(t)dt=f(0) +f?(c)2

Correction :La première chose à rappeler est que la définition fait deFLA PRIMITIVE DEfQUI S"ANNULE EN0. En particulier,Fvérifie Fest de classeC1,?x?[0,1], F?(x) =f(x), F(0) = 0.(1) Comme suggéré dans l"énoncé on peut donc écrire de développement de Taylor

Lagrange deFen0à l"ordre2:

?x?[0,1],?αx?[0,1], F(x) =F(0) + (x-0)F?(0) +(x-0)22

F??(αx)

ce qui donne immédiatement d"après ?x?[0,1],?αx?[0,1], F(x) =f(0)x+x22 f?(αx). En particulier enx= 1, en reprenant la définition deF(1), on a doncc?[0,1]tel que

F(1) =?

1

0f(x)dx=f(0) +12

f?(c) (précisément ce qu"il fallait démontrer). 2.

On prop osed "approcherl"in tégraleI(f) =?1

0f(t)dtpar une formule du type

J(f) =f(0) +f?(α)2

pourα?]0,1[à déterminer. a. Mon trerque Jest exacte pour les polynômes de degré inférieur ou égal à1 quelque soit le choix deα?]0,1[. Correction :Pour montrer queJest exacte pour les polynômes de degré inférieur ou égal à1, il faut (et il suffit) de montrer qu"elle est exacte pour la base canoniqu eueR1[X], c"est-à-dire

J(1) = 1etJ(X) =?

1

0xdx=12

Or, par définition deJ:

-J(1) = 1 + 0 = 1 4 -J(X) = 0 +12 =12 b. Déterminer αpour que l"approximationJ(f)soit exacte pour les polynômes de degré au plus égal à deux. Correction :Pour queJsoit de plus exacte pour les polynômes de degré2, d"après ce qui précède il suffit d"avoirJ(X2) =?1

0x2dx=13

Or,

J(X2) = 0 +2α2

La condition nécessaire et suffisante est donc queα=13 c. P ourle c hoixde αde la question 2-b), quel est l"ordre d"exactitude de la méthode? Correction :On a vu à la question précédente que le degré d"exactitude est au moins2, il faut donc tester pour3,4etc jusqu"à trouver quel est le plus au degré pour lequelJest exacte. On va vérifier siJest exacte à l"ordre 3 :

J(X3) = 03+3×α22

=3×13 22
=16 et 1

0t3dt=14

. Ainsi,J(X3)?=? 1

0t3dtdoncJest de degré d"exactitude 2.

3. On fixe p ourla suite αcomme à la question 2-b). On suppose quef?C3([0,1]). A l"aide d"un développement de Taylor deF(x), montrer qu"il existed?]0,1[tel que

I(f) =f(0) +12

f?(0) +16 f??(0) +124 f(3)(d). A l"aide d"un développement de Taylor def?(α), montrer alors que |I(f)-J(f)|6572 sup x?[0,1]|f(3)(x)|. Correction :La fonctionfest maintenant de classeC3, ce qui fait deFune fonction de classeC4. On peut donc effectuer un développement de Taylor Lagrange

à l"ordre4en0et on obtient :

?x?[0,1],?αx?[0,1], F(x) =F(0)+xF?(0)+x22

F??(0)+x36

F???(0)+x424

F(4)(αx).

En particulier, enx= 1et en utilisant les propriétés (1) deF, on obtient qu"il existed=α1?[0,1]tel que

I(f) =F(1) =f(0) +12

f?(0) +16 f??(0) +124 f???(d) ce qu"il fallait démontrer dans un premier temps. On veut maintenant évaluerI(f)-J(f). CommeJ(f)fait intervenirf?(α), il 5 paraît raisonnable de faire un développement de Taylor Lagrange def?(α)en0à l"ordre maximal (i.e2puisquef?est de classeC2) : fquotesdbs_dbs11.pdfusesText_17

[PDF] intégration globale

[PDF] intégration globale définition

[PDF] integration globale en anglais

[PDF] intégration globale ifrs

[PDF] intégration horizontale et verticale management

[PDF] intégration horizontale exemple d'entreprise

[PDF] intégration latérale définition

[PDF] intégration nouveau salarié dans l'entreprise

[PDF] intégration numérique méthode de simpson

[PDF] integration numerique methode de trapeze pdf

[PDF] intégration ou inclusion éléments pour contribuer au débat

[PDF] intégration proportionnelle 99 02

[PDF] intégration proportionnelle exemple

[PDF] intégration proportionnelle ifrs

[PDF] intégration proportionnelle normes françaises