[PDF] Analyse Numérique PDF CTD3.pdf

Proposition de corrigé du TD 3 EXERCICE 1 Interpolation de Lagrange Soit x0, x1, , xn, n + 1 points distincts a Soit (Li)i=0,n n + 1 fonctions de Pn vérifiant

L2 Maths, UE d'Analyse numérique Feuille de TD 1 - Correction : Interpolation de Lagrange Exercice 1 (Identification) On considère x, y ∈ R4 donnés par : x

[PDF] Série dexercices no5/6 Interpolation polynomiale

Autrement dit, connaissant Pn1 , il suffit de calculer an pour connaître Pn a) Montrer que le polynôme d'interpolation de Lagrange de la fonction f aux points

[PDF] Analyse Numérique

Proposition de corrigé du TD 3 EXERCICE 1 Interpolation de Lagrange Soit x0, x1, , xn, n + 1 points distincts a Soit (Li)i=0,n n + 1 fonctions de Pn vérifiant

[PDF] TD Interpolation

Montrez par récurrence que Pk,j avec k ≥ j est le polynôme d'interpolation de Lagrange pour les points x0, x1, , xj−1, xk 3 Qu'en concluez-vous pour Pk,k ?

[PDF] Analyse numérique Exercices corrigés

Exercice 3 Avec quelle précision peut-on calculer √ 115 `a l'aide de l' interpolation de Lagrange, si on prend les points : x0 = 100, x1 = 121, x2 = 144 Corrigé :

[PDF] Réponses aux exercices du chapitre 5

Le polynôme d'interpolation par Lagrange est donné par : pn(x) = n ∑ i=0 f(xi)Li( x) où les (n + 1) fonctions Li(x) sont définies par : Li(x) = (x - x0)···(x - xi-1)(x

[PDF] Exercices avec corrigé succinct du chapitre 5 - UTC - Moodle

les points ti = i/n, i = 0,1, ,n, `a l'aide du polynôme d'interpolation de Lagrange de degré n Expliquer le résultat 2 Même question pour la fonction g(t) = tn+1

[PDF] Analyse - Gloria FACCANONI

2 juil 2010 · [3 pt] Construire le polynôme de Lagrange P qui interpole les points (0, 2), (1, 1), (2, 2) et (3, 3) Exercice ƒ : interpolation polynomiale [2 pt] [2 pt]

[PDF] Exercices corrigés

Exercices corrigés NB : Les exercices corrigés ici sont les exercices proposés durant les et li(x), polynôme de Lagrange nécessaire pour l'interpolation

Universit´e de Nice Sophia-Antipolis

Licence L3 Math´ematiques Ann´ee 2008/2009Analyse Num´erique

Proposition de corrig´e du TD 3EXERCICE 1

Interpolation de Lagrange

Soitx0, x1, ..., xn,n+ 1points distincts.

a. Soit(Li)i=0,nn+ 1fonctions dePnv´erifiantLi(xj) =δij. Montrer que (Li)i=0,nest une base dePn(ensemble des polynˆomes de degr´e inf´erieur ou

´egal `an). Construire cette base.

(Li)i=0,nest une base dePn •On a (Li)i=0,n? Pn. •On aCard? (Li)i=0,n? =n+ 1 =dimPn. •Soit (ai)i=0,n?Rtels que n? i=0a iLi(x) = 0.

Donc 0 =

?n i=0aiLi(xj) =ajpour chacun desj? {0,...,n}, ou encoreaj= 0 pour chaquej. D"o`u la famille (Li)i=0,nest libre.

Par suite la famille (Li)i=0,nest une base dePn.

Construction de la base(Li)i=0,nSoiti? {0,...,n}. Pour toutj? {0,...,n}j?=i, Li(xj) = 0. Donc L i(x) =ciΠnj=0,j?=i(x-xj).

DeLi(xi) = 1, on d´eduit

c i=1Π nj=0,j?=i(xi-xj). D"o`u L i(x) =Πnj=0,j?=i(x-xj)Π nj=0,j?=i(xi-xj)= Πnj=0,j?=i?x-xjx i-xj? b. Soitpn? Pnv´erifiant :pn(xi) =f(xi)?i= 0,...,n. D´ecomposerpnsur la base des(Li)i=0,n. Un telpnest-il unique? 1

Universit´e de Nice Sophia-Antipolis

Licence L3 Math´ematiques Ann´ee 2008/2009D´ecomposition depnsur la base(Li)i=0,nOn a p n(x) =n? i=0a iLi(x) avecaj?R.

Depn(xj) =f(xj)?j= 0,...,n, on obtient

p n(x) =n? i=0f(xi)Li(x).

Unicit´e depn

Soientpn, qn? Pntels quepn(xi) =f(xi)?i= 0,...,netqn(xi) =f(xi)?i= 0,...,n. n´ecessairementr= 0. c. ´Ecrire le polynˆome d"interpolation associ´e aux points donn´es dans le tableau suivant :x i-1-1/201/21 f(xi)-3/201/400

Tab.1 - Tableau pour l"interpolation.

On a p

4(x) =f(x0)L0(x) +f(x1)L1(x) +f(x2)L2(x) +f(x3)L3(x) +f(x4)L4(x),

=-32

L0(x) +14

f(x2)L2(x). o`u L

0(x) =(x+12

)x(x-12 )(x-1)(-1 +12 )(-1-0)(-1-12 )(-1-1)=x4-x3-14 x2+14 x3 2 L

2(x) =(x+ 1)(x+12

)(x-12 )(x-1)(0 + 1)(0 + 12 )(0-12 )(0-1)=x4-54 x2+141 4 D"o`u p

4(x) =-32

4-x3-14

x2+14 x3 2 14 x 4-54 x2+141 4 =x3-x2-14 x+14 2

Universit´e de Nice Sophia-Antipolis

Licence L3 Math´ematiques Ann´ee 2008/2009d. ´Etablir la majoration d"interpolation de Lagrangei.e.sif? Cn+1([a,b]), alors il existeξ?]a,b[tel que f(x)-pn(x) =Πnj=0(x-xj)(n+ 1)!f(n+1)(ξ).(1.1) •Six=xj?j= 0,...,n, alorsf(x)p(x) = 0, et toutξ?]a,b[ convient. •Six?=xj?j= 0,...,n, alors d´efinissons

φ(t) =f(t)-p(t)-k(x) Πnj=0(t-xj),?t?[a,b],

o`uk(x) est choisi de telle sorte queφ(x) = 0.

D"une part, on en d´eduit que

k(x) =f(x)-p(x)Π nj=0(x-xj).(1.2) D"autre part, la fonctiont?→φ(t) est de classeCn+1([a,b]) admet (n+2) racinesdistinctes x, x

0, x1, ..., xnsur ]a,b[. D"apr`es le th´eor`eme de Rolle :

t?→φ?(t) est de classeCn([a,b]) admet (n+ 1) racinesdistinctessur ]a,b[, appartenant chacune entre les intervalles ouverts d"extr´emit´es dex, x0, x1, ..., xncontenus dans ]a,b[.

Par application du th´eor`eme de Rolle,

t?→φ(2)(t) est de classeCn-1([a,b]) admetnracinesdistinctessur ]a,b[. Par application du th´eor`eme de Rolle une nouvelle fois, t?→φ(3)(t) est de classeCn-2([a,b]) admetn-1 racinesdistinctessur ]a,b[.

Ainsi de suite, par application du th´eor`eme de Rolle,t?→φ(n+1)(t) est de classeC0([a,b])

admet une racineξ?]a,b[,φ(n+1)(ξ) = 0.

Puisquepn? Pn, on ap(n+1)n= 0 et

0 =φ(n+1)(ξ) =f(n+1)(ξ)-(n+ 1)!k(x).(1.3)

Donc de (1.2) et (1.3) on tire,

k(x) =f(x)-p(x)Π nj=0(x-xj)=f(n+1)(ξ)(n+ 1)!.

D"o`u le r´esultat.

e. Soientf(x) = cos(x)etg(x) =e3xd´efinies sur[0,1]. Estimer le nombre minimum de points pour que l"erreur entre la fonction et son polynˆome d"in- terpolation de Lagrange soit inf´erieure `a0.1,0.01et0.001. 3

Universit´e de Nice Sophia-Antipolis

Licence L3 Math´ematiques Ann´ee 2008/2009Nombre de points mimimum pour satisfaire une tol´eranceεdonn´ee

Premier cas:f(x) = cos(x) sur [0,1].

On a cos

(n+1)(x) = cos? x+ (n+ 1)π2 donne (n+ 1)!≥ε-1. D"o`u •ε= 0.1?ε-1= 10 :3! = 6,

4! = 24?n+ 1≥4?n≥3.

•ε= 0.01?ε-1= 100 :

5! = 120?n+ 1≥5?n≥4.

•ε= 0.001?ε-1= 1000 :

6! = 720,

7! = 5040?n+ 1≥7?n≥6.

Deuxi`eme cas:g(x) = exp(3x) sur [0,1].

D"o`u •ε= 0.1 :

3(9+1)exp(3)(9 + 1)!

?0.3268383, 3 (10+1)exp(3)(10 + 1)! ?0.0891377?n≥10. •ε= 0.01 :

3(11+1)exp(3)(11 + 1)!

?0.0222844, 3 (12+1)exp(3)(12 + 1)! ?0.0051426?n≥12. •ε= 0.001 :

3(13+1)exp(3)(13 + 1)!

?0.0011020, 3 (14+1)exp(3)(14 + 1)! ?0.0002204?n≥14. 4

Universit´e de Nice Sophia-Antipolis

Licence L3 Math´ematiques Ann´ee 2008/2009EXERCICE 2

Interpolation de Hermite

Soitf? C1([a,b])etx1, x2deux points distincts. Soitpun polynˆome de a. Montrer qu"un tel polynˆome existe et est unique.

Existence

On posep(x) =a3x3+a2x2+a1x+a0, doncp?(x) = 3a3x2+ 2a2x+a1.

Les conditions surpetp?s"´ecrivent :

3x31+a2x21+a1x1+a0=f(x1),

3x32+a2x22+a1x2+a0=f(x2),

3a3x21+ 2a2x1+a1+ 0 =f?(x1),

3a3x22+ 2a2x2+a1+ 0 =f?(x2).

ou encore ((((((((1x1x21x31

1x2x22x32

0 1 2x13x21

0 1 2x23x22)

((((((((a 0 a 1 a 2 a 1) ((((((((f(x1) f(x2) f ?(x1) f ?(x2)) ou bien encore AX=B avec A=( ((((((((1x1x21x31

1x2x22x32

0 1 2x13x21

0 1 2x23x22)

)))))))), X=( ((((((((a 0 a 1 a 2 a 1) )))))))), C=( ((((((((f(x1) f(x2) f ?(x1) f ?(x2)) 5

Universit´e de Nice Sophia-Antipolis

Licence L3 Math´ematiques Ann´ee 2008/2009On adet(A) =-? x 2-x1?

4?= 0 puisquex1?=x2, d"o`u l"existence dep.

Unicit´e dep

Soientp, q? P3tels que

p(xi) =f(xi), i= 1,2, p ?(xi) =f?(xi), i= 1,2.

Soitr=p-q. Alors

r(xi) =f(xi), i= 1,2, r ?(xi) =f?(xi), i= 1,2,?r(x) =C(x)(x-x1)2(x-x2)2. Commer? P3, on ac(x) = 0, doncr= 0, puisp=q. D"o`u l"unicit´e. b. ´Etablir la majoration d"interpolation suivante : sif? C4([a,b]), alors il existeξ?]a,b[tel que f(x)-p(x) =(x-x1)2(x-x2)24! f(4)(ξ). •Six=x1oux2, alorsf(x)-p(x) = 0, et toutξ?]a,b[ convient. •Six?=x1,x2, on pose

φ(t) =f(t)-p(t)-(t-x1)(t-x2)k(x),?t?[a,b],

o`uk(x) est choisi de telle sorte queφ(x) = 0. Donc k(x) =f(x)-p(x)(x-x1)(x-x2),(2.1)

On a ´egalement :

la fonctiont?→φ(t) est de classeC4([a,b]) et admet 3 racinesdistinctesx, x1, x2sur [a,b]. D"apr`es le th´eor`eme de Rolle : t?→φ?(t) est de classeC3([a,b]) et s"annule en 2 points distinctsc1,c2?=x, x1, x2, c

1,c2?]min(x,x1,x2),max(x,x1,x2)[.

De plus

?(t) =f?(x)-p?(x)-k(x)?

2(t-x1)(t-x2)2+ 2(t-x1)2(t-x2)?

ce qui entraˆıneφ?(x1) = 0,φ?(x2) = 0. t?→φ?(t) s"annule en4 points distinctsc1,c2,x1, x2. t?→φ??(t) est de classeC2([a,b]) et admet 3 racinesdistinctesd1,d2,d3chacune appar- tenant `a l"intervalle ]zi,zk[ o`uzi,zk? {c1,c2,x1, x2}. t?→φ(3)(t) est de classeC1([a,b]) et admet 2 racinesdistinctese1,e2chacune appartenant `a l"intervalle ]yi,yk[ o`uyi,yk? {d1,d2,d3}. t?→φ(4)(t) est de classeC0([a,b]) et admet 1 racineξ?]e1,e2[?]a,b[. 6

Universit´e de Nice Sophia-Antipolis

Licence L3 Math´ematiques Ann´ee 2008/2009On a

0 =φ(4)(ξ) =f(4)(ξ)-p(4)(ξ)-(4!)k(x).(2.2)

Commep? P3, on ap(4)= 0. Des ´equations (2.1) et (2.2), on d´eduit f(x)-p(x)(x-x1)2(x-x2)2=k(x) =f(4)(ξ)4! ,pourx?=x1,x2.(2.3)

D"o`u le r´esultat.

c. Trouver une base(A1, A2, B1, B2)deP3telle que p(x) =f(x1)A1(x) +f(x2)A2(x) +f?(x1)B1(x) +f?(x2)B2(x). et exprimer cette base en fonction des polynˆomes d"interpolation de Lagrange L

1etL2.

La conditionp(x1) =f(x1) s"´ecrit :

f(x1)A1(x1) +f(x2)A2(x1) +f?(x1)B1(x1) +f?(x2)B2(x1) =f(x1), ou encore f(x1)? A

1(x1)-1?

+f(x2)A2(x1) +f?(x1)B1(x1) +f?(x2)B2(x1) = 0.

Ce qui s"´ecrit encore :

??A

1(x1) = 1,

quotesdbs_dbs8.pdfusesText_14

[PDF] [PDF] Analyse Numérique

Universit´e de Nice Sophia-Antipolis

Proposition de corrig´e du TD 3EXERCICE 1

Interpolation de Lagrange

Soitx0, x1, ..., xn,n+ 1points distincts.

´egal `an). Construire cette base.

Donc 0 =

Par suite la famille (Li)i=0,nest une base dePn.

DeLi(xi) = 1, on d´eduit

Universit´e de Nice Sophia-Antipolis

Depn(xj) =f(xj)?j= 0,...,n, on obtient

Unicit´e depn

Tab.1 - Tableau pour l"interpolation.

4(x) =f(x0)L0(x) +f(x1)L1(x) +f(x2)L2(x) +f(x3)L3(x) +f(x4)L4(x),

L0(x) +14

0(x) =(x+12

2(x) =(x+ 1)(x+12

4(x) =-32

4-x3-14

Universit´e de Nice Sophia-Antipolis

φ(t) =f(t)-p(t)-k(x) Πnj=0(t-xj),?t?[a,b],

D"une part, on en d´eduit que

0, x1, ..., xnsur ]a,b[. D"apr`es le th´eor`eme de Rolle :

Par application du th´eor`eme de Rolle,

Puisquepn? Pn, on ap(n+1)n= 0 et

0 =φ(n+1)(ξ) =f(n+1)(ξ)-(n+ 1)!k(x).(1.3)

Donc de (1.2) et (1.3) on tire,

D"o`u le r´esultat.

Universit´e de Nice Sophia-Antipolis

Premier cas:f(x) = cos(x) sur [0,1].

On a cos

4! = 24?n+ 1≥4?n≥3.

5! = 120?n+ 1≥5?n≥4.

6! = 720,

7! = 5040?n+ 1≥7?n≥6.

Deuxi`eme cas:g(x) = exp(3x) sur [0,1].

3(9+1)exp(3)(9 + 1)!

3(11+1)exp(3)(11 + 1)!

3(13+1)exp(3)(13 + 1)!

Universit´e de Nice Sophia-Antipolis

Interpolation de Hermite

Existence

Les conditions surpetp?s"´ecrivent :

3x31+a2x21+a1x1+a0=f(x1),

3x32+a2x22+a1x2+a0=f(x2),

3a3x21+ 2a2x1+a1+ 0 =f?(x1),

3a3x22+ 2a2x2+a1+ 0 =f?(x2).

1x2x22x32

0 1 2x13x21

0 1 2x23x22)

1x2x22x32

0 1 2x13x21

0 1 2x23x22)

Universit´e de Nice Sophia-Antipolis

4?= 0 puisquex1?=x2, d"o`u l"existence dep.

Unicit´e dep

Soientp, q? P3tels que

Soitr=p-q. Alors

φ(t) =f(t)-p(t)-(t-x1)(t-x2)k(x),?t?[a,b],

On a ´egalement :

1,c2?]min(x,x1,x2),max(x,x1,x2)[.

De plus

2(t-x1)(t-x2)2+ 2(t-x1)2(t-x2)?

Universit´e de Nice Sophia-Antipolis

0 =φ(4)(ξ) =f(4)(ξ)-p(4)(ξ)-(4!)k(x).(2.2)

D"o`u le r´esultat.

1etL2.

La conditionp(x1) =f(x1) s"´ecrit :

1(x1)-1?

Ce qui s"´ecrit encore :

1(x1) = 1,