[PDF] Exercices corrigés PDF exo_corriges

Si vous avez des questions concernant ces exercices, n'hésitez pas à envoyer un mail à votre enseignant d'analyse numérique pour lui poser une question Si

Previous PDF

Next PDF

[PDF] Analyse Numérique

Analyse Numérique Corrigé du TD 4 EXERCICE 1 Formule des trap` une estimation de l'erreur d'intégration Estimation de l'erreur d'intégration élémentaire

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

Exercice 2 (De l'interpolation à l'intégration numérique) Soit f : [−1,1] −→ R une fonction de classe C2

[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] Exercices corrigés

Si vous avez des questions concernant ces exercices, n'hésitez pas à envoyer un mail à votre enseignant d'analyse numérique pour lui poser une question Si

[PDF] Analyse numérique, Matmeca 1ere année Corrigé de la feuille 5 I

Analyse numérique, Matmeca ‚ere année Corrigé de la feuille 5I ˜™ defgilE nEp chée et l'intégration exacte est égale à l'intégrale de l'erreur d'interpolation D' apr s le cours On a vu à l'exercice que la formule de Simpson sur [ −1;1] est

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

soit une formule d'intégration numérique exacte pour les polynômes de degré inférieur ou égal à 4 sur [-1,1] (b) (2 points) Cette méthode est-elle d'ordre 5? D'

[PDF] Recueil dexercices pour les cours MTH2210x

Exercices pour les cours de calcul scientifique pour ingénieurs MTH2210B Édition du Proposer une méthode d'intégration numérique pour estimer le temps

[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] CORRIGÉS DES TRAVAUX DIRIGÉS DE lUE - Jerome BASTIEN

http://utbmjb chez-alice fr/Polytech/MNB/TDcorMNB pdf Intégration 21 Ce polycopié constitue les corrigés de TD de Méthodes Numériques de Base du Dans cet exercice, on a montré que l'erreur d'interpolation entre f et son polynôme

[PDF] TD no 3 : Intégration numérique

sur un intervalle [a, b] et donner la valeur de l'erreur Exercice 2 Interpolation/ intégration numérique utilisant les polynômes de Hermite a) Soit αi,i = 1 4

Grenoble INP - Pagora Analyse numérique

1ère année

Exercices corrigés

NB : Les exercices corrigés ici sont les exercices proposés durant les séances de cours. Les corrections données

sont des corrections plus détaillées que celles fournies durant le cours (si le temps a permis de donner ces

corrections). Si vous avez des questions concernant ces exercices, n"hésitez pas à envoyer un mail à votre

enseignant d"analyse numérique pour lui poser une question. Si vous trouver des coquilles, des erreurs dans

le présent document, n"hésitez pas à le signaler à votre enseignant par un mail.

Chapitre 1 : Introduction au calcul approché

Exercice 1Montrer que9325s"écrit bien(10010001101101)2en base2puis reconvertir(10010001101101)2

en base10.Pour convertir un entier de la base10à la base2(on verra que la méthode diffère légèrement pour un

nombre décimal un peu plus tard), on divise l"entier par2(division euclidienne) et le reste correspond au

dernier chiffre de l"entier en base2. Pour9325, cela donne

9325 = 24662 +1

et on itère le processus sur le quotient obtenu (jusqu"à ce qu"il vaille1). Ainsi puisque

4662 = 22331 +0

On peut réécrire9325sous la forme

9325 = 2(22331 + 0) + 1 = 222331 + 210+ 201

et01sont les 2 derniers chiffres de9325écrit en base2. Pour enfoncer le clou, on détaille encore l"itération

9325 = 2

2(21165 + 1) + 210+ 201= 231165 + 221+ 210+ 201

et101sont les 3 derniers chiffres de9325écrit en base2. On affiche ensuite le processus itératif dans son

entier :9325 = 24662 +1

4662 = 22331 +0

2331 = 21165 +1

1165 = 2582 +1

582 = 2291 +0

291 = 2145 +1

145 = 272 +1

72 = 236 +0

36 = 218 +0

18 = 29 +0

9 = 24 +1

4 = 22 +0

2 = 21 +0

1 = 20 +1

1 Donc, on peut décomposer9325de la manière suivante

9325 = 2

131+ 2120+ 2110+ 2101+ 290+ 280+ 270

+ 2

61+ 251+ 240+ 231+ 221+ 210+ 201

On a bien montré que(9325)10= (10010001101101)2. Il ne reste plus qu"? reconvertir ce nombre binaire en

base10. Pour ce faire, on va procéder de manière itérative. On commence par cette première étape,1est

le chiffre le plus fort (le plus à gauche) de(10010001101101)2et on construit le résultat intermédiaire de la

manière suivante, on multiplie par 2 le résultat intermédiaire précédent (au départ 0) et on ajoute le chiffre

le plus fort restant à traiter. On commence donc par

20 +1= 1 = 201

Le résultat intermédiaire est donc1et il ne reste plus qu"à traiter0010001101101du binaire (10010001101101)

2. On itère le processus. On multiplie par2le résultat intermédiaire (ici 1) puis on ajoute

le chiffre le plus fort restant à traiter (soit ici0). D"où

21 +0= 2 = 211+ 200

Il ne reste plus qu"à traiter010001101101du binaire(10010001101101)2. Si on détaille l"étape suivante, on a

22 +0= 4 = 221+ 210+ 200

Il ne reste plus qu"à traiter10001101101du binaire(10010001101101)2. On affiche ensuite le processus itératif

dans son entier :20 +1= 1

21 +0= 2

22 +0= 4

24 +1= 9

29 +0= 18

218 +0= 36

236 +0= 72

272 +1= 145

2145 +1= 291

2291 +0= 582

2582 +1= 1165

21165 +1= 2331

22331 +0= 4662

24662 +1= 9325

On vérifie donc bien que(9325)10= (10010001101101)2. Notons bien que de processus décrit ici est juste le

premier processus mais pris en sens inverse.

Exercice 2Écrire(34)10et(27)10en binaire puis effectuer l"opération en binaire(34)10+(27)10et vérifier

que le résultat obtenu soit le bon. 2 Convertissons tout d"abord34en binaire. Cela donne

34 = 217 +0

17 = 28 +1

8 = 24 +0

4 = 22 +0

2 = 21 +0

1 = 20 +1

On a donc(34)10= (100010)2. Convertissons maintenant27en binaire. On a

27 = 213 +1

13 = 26 +1

6 = 23 +0

2 = 21 +1

1 = 20 +1

et(27)10= (11011)2. On effectue maintenant l"addition de(100010)2et(11011)2. Pour rappel, l"addition en

binaire fonctionne de la manière suivante+01 001 1110

D"où l"opération suivante

1 0 0 0 1 0

+ 1 1 0 1 1= 1 1 1 1 10 1 On a(100010)2+ (11011)2= (111101)2. Or(34)10+ (27)10= (61)10, vérifions si(61)10= (111101)2.

20 +1= 1

21 +1= 3

23 +1= 7

27 +1= 15

215 +0= 30

230 +1= 61

On a bien(61)10= (111101)2, le résultat obtenu en binaire est bien conforme au résultat obtenu en base10.

Exercice 3Écrire(90)10et(97)10en binaire puis effectuer l"opération en binaire(90)10(97)10et vérifier

que le résultat obtenu est le bon.Convertissons tout d"abord90en binaire. Cela donne

90 = 245 +0

45 = 222 +1

22 = 211 +0

11 = 25 +1

5 = 22 +1

2 = 21 +0

1 = 20 +1

3 On a donc(90)10= (1011010)2. Convertissons maintenant97en binaire. On a

97 = 248 +1

48 = 224 +0

24 = 212 +0

12 = 26 +0

6 = 23 +0

3 = 21 +1

1 = 20 +1

et(97)10= (1100001)2. On effectue maintenant la multiplication de(1011010)2par(1100001)2. Pour rappel,

la multiplication en binaire fonctionne de la manière suivante01 000 101

D"où l"opération suivante

1 0 1 1 0 1 0

1 1 0 0 0 0 11 0 1 1 0 1 0

+ 0 0 0 0 0 0 00 + 0 0 0 0 0 0 00 0 + 0 0 0 0 0 0 00 0 0 + 0 0 0 0 0 0 00 0 0 0 + 1 0 1 1 0 1 00 0 0 0 0 + 1 0 1 1 0 1 00 0 0 0 0 0= 1

10101011101010 0 1 1 0 1 0

On a(1011010)2(1100001)2= (10001000011010)2. Or(90)10(97)10= (8730)10, vérifions si(8730)10= (10001000011010) 2.

20 +1= 1

21 +0= 2

22 +0= 4

24 +0= 8

28 +1= 17

217 +0= 34

234 +0= 68

268 +0= 136

2136 +0= 272

2272 +1= 545

2545 +1= 1091

21090 +0= 2182

22182 +1= 4365

24365 +0= 8730

On a bien(8730)10= (10001000011010)2, le résultat obtenu en binaire est bien conforme au résultat obtenu

en base10. 4

Exercice 4Si on dispose de4bits (bit de signe compris), quelles valeurs peuvent prendre les entiers codés

sur ces4bits?Si on dispose de4bits dont1de signe, il ne reste plus que3bits pour coder les entiers naturels (ceux plus

grand que0). Ils ne peuvent donc prendre que23valeurs distinctes dont la valeur0. Les entiers naturels

codés sont ainsi0,1,2,3,4,5,6, et7 = 231. Maintenant, si on tient compte du bit de signe, les entiers

codés devraient pouvoir varier entre7et7.

Cependant deux combinaisons auraient la même valeur0:1000et0000, le chiffre en gras désigne ici le bit

de signe. Pour éviter cette redondance, on pose1000 =8(classiquement, le signe de bit lorsqu"il vaut1

indique un nombre négatif).

Finalement, si on dispose de4bits (bit de signe compris), on peut coder les entiers de valeurs comprises

entre8 =23et7 = 231.

Exercice 5Vérifier l"égalité entre(9;90625)10et(1001;11101)2.On distingue la partie entière et la partie décimale à traiter. On vérifier tout d"abord que(9)10= 10012, en

effet9 = 24 +1

4 = 22 +0

2 = 21 +0

1 = 20 +1

On a donc

9;90625 = 231+ 220+ 210+ 201+ 0;90625

Mais on a aussi

9;90625 = 231+ 220+ 210+ 201+ 21(20;90625)

9;90625 = 231+ 220+ 210+ 201+ 211;8125

9;90625 = 231+ 220+ 210+ 201+ 21(1+ 0;8125)

On vient donc de calculer le premier chiffre après la virgule de9;90625en binaire (soit ici1). On réitère le

même processus pour avoir le chiffre après la virgule suivant

9;90625 = 231+ 220+ 210+ 201+ 211+ 22(20;8125)

9;90625 = 231+ 220+ 210+ 201+ 211+ 221;625

9;90625 = 231+ 220+ 210+ 201+ 211+ 22(1+ 0;625)

Le deuxième chiffre après la virgule (en binaire) est donc1. Voici enfin directement, les traces des calculs

pour obtenir tous les chiffres nécessaires après la virgule

9;90625 = 231+ 220+ 210+ 201+ 211+ 221+ 231;25

9;90625 = 231+ 220+ 210+ 201+ 211+ 221+ 23(1+ 0;25)

9;90625 = 231+ 220+ 210+ 201+ 211+ 221+ 231+ 240;5

9;90625 = 231+ 220+ 210+ 201+ 211+ 221+ 231+ 240+ 25(20;5)

9;90625 = 231+ 220+ 210+ 201+ 211+ 221+ 231+ 240+ 251

On vérifie donc bien que(9;90625)10= (1001;11101)2. 5 Chapitre 2 : Résolution d"équations non-linéaires Exercice 6On définit la méthode du point fixe suivante x0fixé dans[a;b] x n+1=g(xn)

On suppose que cette suite admet une limite sur[a;b]notéel. Cette méthode est d"ordrepsijxn+1ljjxnljp

admet une limite réelle strictement positive lorsquentend vers l"infini.

On supposeg pfois dérivable sur[a;b]. En utilisant la formule de Taylor, montrer que la méthode est d"ordre

psi et seulement si g

0(l) =g00(l) =:::=g(p1)(l) = 0 etg(p)(l)6= 0On rappelle tout d"abord la formule de Taylor.

Soitk1un entier et soitfune fonction deRdansRkfois dérivable ena2R, alors il existe une fonction kdeRdansRtel que f(x) =f(a) +f0(a)(xa) +f00(a)(xa)22 +:::+f(k)(a)(xa)kk!+k(x)(xa)k et limx!ak(x) = 0 Maintenant, sigestpfois dérivable sur[a;b], on peut écrire que g(xn) =g(l) +g0(l)(xnl) +g00(l)(xnl)22 +:::+g(p)(l)(xnl)pp!+p(x)(xnl)p avec limx n!lp(xn) = 0

Posonsen=xnl, on axn+1=g(xn)et

e n+1=g(xn)g(l) =g0(l)(xnl) +g00(l)(xnl)22 +:::+g(p)(l)(xnl)pp!+p(xn)(xnl)p =g0(l)en+g00(l)e2n2 +:::+g(p)(l)epnp!+p(x)epn Si la métho deest d"ordre palors le quotienten+1e pntend vers un réel non-nul (carjxn+1ljjxnljpadmet une

limite réelle strictement positive et on enlève juste les valeurs absolues). Or d"après ce qui précède, on

a l"égalité suivante e n+1e pn=g0(l)e p1n+g00(l)2ep2n+:::+1e ng (p1)(l)(p1)!+g(p)(l)p!+p(xn)

A gauche de l"égalité, le terme tend vers un réel non-nul. Concernant les termes à droite, on a queen

tend vers0quandntend vers l"infini donc lim n!+1g 0(l)e p1n=0 sig0(l) = 0

1sinon

De même

lim n!+1g

00(l)2ep2n=0 sig00(l) = 0

1sinon

jusqu"à lim n!+11e ng (p1)(l)(p1)!=0 sig(p1)(l) = 0

1sinon

Par contre

g(p)(l)p!est un réel etp(xn)tend vers0. Par unicité de la limite des deux côtés de l"égalité,

il faut queg0(l) =g00(l) =:::=g(p1)(l) = 0pour que le côté droit tende vers un réel,g(p)(l)p!en

l"occurence. La limite réelle est en plus différente de0, donc on doit avoirg(p)(l)6= 0.

Si on a

0(l) =g00(l) =:::=g(p1)(l) = 0 etg(p)(l)6= 0

alors e n+1=epnp!g(p)(l) +epnp(xn) et en+1e pn=1p!g(p)(l) +p(xn) d"où commep(xntend vers0quandntend vers l"infini lim n!+1jen+1e pnj=1p!jg(p)(l)j 2R+ et la méthode est bien d"ordrep. Exercice 7 (ordre de convergence de la méthode de Newton)On rappelle ici la méthode de New-

ton, il s"agit d"une méthode du point fixe pour rśoudref(x) = 0sur[a;b]définissant la suite suivante8<

0fixé dans[a;b]

x n+1=g(xn) =xnf(xn)f 0(xn)

On suppose que cette suite admet une limite sur[a;b]notéel. Montrer que sifest3fois dérivable sur[a;b]

et quef0(l)6= 0, alors la méthode de Newton est d"ordre2au moins.La fonctiongvaut ici g(x) =xf(x)f 0(x)

Sa dérivée vaut

0(x) = 1f0(x)f0(x)f(x)f00(x)(f0(x))2=f(x)f00(x)(f0(x))2

et enl, cette dérivée vaut0carf(l) = 0,lest la solution au problème étudié. Donc la méthode est d"ordre2

au moins (d"ordre2sig00(l)6= 0, d"ordre supérieur à 2 sinon). Si on va plus loin, calculons la dérivée seconde

deg g

00(x) =(f0(x)f00(x) +f(x)f000(x))(f0(x))22f(x)f0(x)(f00(x))2(f0(x))4

Elle vaut enl

00(l) =f00(l)f

0(l)

Donc sif00(l)6= 0alors la méthode de Newton est d"ordre2, sinon elle est d"ordre supérieur à 2.

Chapitre 3 : Approximation de fonctions

Méthode des moindres carrés

Exercice 8On dispose denpoints(xi;yi),i= 1;:::;net on suppose que la fonction modèle est de la forme

f(x;0) =0

Trouver0minimisant la quantité

S(0) =nX

i=1(yif(xi;0))2On cherche le minimum0tel queS(0) = min

02RS(0). Pour cela, on a besoin de calculer sa dérivéeS0(0).

Celle-ci vaut

0(0) =nX

i=1[2(yif(xi;0))] =2nX i=1(yi0)

Pour trouver un extremum local (minimum ou maximum local), il faut résoudre le problèmeS0(0) = 0.

Dans notre cas, on en déduit que

0=1n n X i=1y i

Il reste cependant à vérifier que0est bien un minimum et que celui-ci est global. Pour cela, étudions le

comportement deS(0). On a le tableau des variations suivant 0S

0(0)S(0)1

0 00+ +1+1S(0)0S(0)0+1+1Grâce à ce tableau, on vérifie que0=1n n X i=1y iminimise bien la quantitéS(0).

Exercice 9 (Régression linéaire)On dispose denpoints(xi;yi),i= 1;:::;net on suppose que la fonc-

tion modèle est de la forme f(x;a;b) =ax+b 8

On cherche le minimum de

S(a;b) =nX

i=1(yiaxib)2

Pour cela, il faut chercher où le gradient deSvaut0. On doit donc calculer les dérivées partielles deSpar

rapport àaetb. Celles-ci valent

8>>>><

>>>:@S@a =2nX i=1x i(yiaxib) @S@b =2nX i=1(yiaxib) Le minimum deSest donc atteint pour(a;b)solution du système suivant

8>>>><

>>>:n X i=1x i(yiaxib) = 0 n X i=1(yiaxib) = 0

On notexetyles moyennes desxietyi:x=1n

n X i=1x iy=1n n X i=1y i Exprimerben fonction dea,xetypuis exprimera.La seconde ligne du système à résoudre donne n X i=1(yiaxib) = 0()nX i=1y ianX i=1x ibnX i=11 |{z} =n= 0

D"où

b=1n nX i=1y ianX i=1x i! =yax Maintenant réécrivons le système à résoudre :

8>>>><

>>>:n X i=1x i(yiaxib) =nX i=1x i(yiya(xix)) = 0 n X i=1(yiaxib) =nX i=1(yiya(xix)) = 0 Si on multiplie la deuxième ligne parxet qu"on la soustrait à la première ligne, on obtient n X i=1x i(yiya(xix))x nX i=1(yiya(xix)) = 0x0 = 0 9 Or 0 = nX i=1x i(yiya(xix))x nX i=1(yiya(xix)) =nX i=1(xix)(yiya(xix)) nX i=1(xix)(yiy)anX i=1(xix)2

Et finalement, on obtient que

a=n X i=1(xix)(yiy)n X i=1(xix)2

Remarque : si vous obtenez le résultat suivant

a=n X i=1(xiyixy)n X i=1(x2ix)=n X i=1x iyinxy n X i=1x 2inx celui-ci est aussi valide. En effet, on l"égalité suivante pour le numérateur n X i=1(xix)(yiy) =nX i=1(xiyixy ixiy+xy) =nX i=1x iyix nX i=1y i |{z} =ny y nX i=1x i |{z} =nx +nxy=nX i=1x iyinxy et l"égalité suivante pour le dénominateur n X i=1(xix)2=nX i=1(x2i2xx i+x

2) =nX

i=1x 2i2x nX i=1x i |{z} =nx +nx 2=nX i=1x 2inx 2 Exercice 10 (Moindres carrés linéaires)On dispose denpoints(xi;yi),i= 1;:::;net on suppose que la fonction modèle est de la forme f(x;) =mX k=1 kk(x) avecle vecteur contenant tous les paramètres que nous allons chercher =0 B BB@ 1 2... m1 C CCA

On cherche le minimum de

S() =nX

i=1(yif(xi;))2 10

Pour cela, il faut chercher où le gradient deSvaut0. On doit donc calculer les dérivées partielles deSpar

rapport auxk. Celles-ci valent @S()@ k=2nX i=1(yif(xi;))@f(xi;) k= 0quotesdbs_dbs1.pdfusesText_1

[PDF] [PDF] Exercices corrigés

Grenoble INP - Pagora Analyse numérique

1ère année

Exercices corrigés

Chapitre 1 : Introduction au calcul approché

9325 = 24662 +1

4662 = 22331 +0

On peut réécrire9325sous la forme

9325 = 2(22331 + 0) + 1 = 222331 + 210+ 201

2331 = 21165 +1

9325 = 2

2(21165 + 1) + 210+ 201= 231165 + 221+ 210+ 201

4662 = 22331 +0

2331 = 21165 +1

1165 = 2582 +1

582 = 2291 +0

291 = 2145 +1

145 = 272 +1

72 = 236 +0

36 = 218 +0

18 = 29 +0

9 = 24 +1

4 = 22 +0

2 = 21 +0

1 = 20 +1

9325 = 2

131+ 2120+ 2110+ 2101+ 290+ 280+ 270

61+ 251+ 240+ 231+ 221+ 210+ 201

20 +1= 1 = 201

2. On itère le processus. On multiplie par2le résultat intermédiaire (ici 1) puis on ajoute

21 +0= 2 = 211+ 200

22 +0= 4 = 221+ 210+ 200

21 +0= 2

22 +0= 4

24 +1= 9

29 +0= 18

218 +0= 36

236 +0= 72

272 +1= 145

2145 +1= 291

2291 +0= 582

2582 +1= 1165

21165 +1= 2331

22331 +0= 4662

24662 +1= 9325

34 = 217 +0

17 = 28 +1

8 = 24 +0

4 = 22 +0

2 = 21 +0

1 = 20 +1

27 = 213 +1

13 = 26 +1

6 = 23 +0

2 = 21 +1

1 = 20 +1

D"où l"opération suivante

1 0 0 0 1 0

20 +1= 1

21 +1= 3

23 +1= 7

27 +1= 15

215 +0= 30

230 +1= 61

90 = 245 +0

45 = 222 +1

22 = 211 +0

11 = 25 +1

5 = 22 +1

2 = 21 +0

1 = 20 +1

97 = 248 +1

48 = 224 +0

24 = 212 +0

12 = 26 +0

6 = 23 +0

3 = 21 +1

1 = 20 +1

D"où l"opération suivante

1 0 1 1 0 1 0