Scribd
Importer
3K vues40 pages

14 Examens en Analyse Numérique + Corrigé PDF

Transféré par

Saad Ettahiri
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF ou lisez en ligne sur Scribd
3K vues40 pages

14 Examens en Analyse Numérique + Corrigé PDF

Transféré par

Saad Ettahiri
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF ou lisez en ligne sur Scribd
Ecole Nationale Préparatoire aux Etudes d’Ingéniorat 3eme Année Préparatoire ‘Année Scolaire 2001/2002 Partiel n° Module; Analyse Numérique Semestre: 2 Date : 09/04/2002 Durée : 02h Bardme : Ex! : 08pts, Ex 2 Odpts, Ex 3 : 04 pts , Exd - 04 pis Exercice t : Soit I’équation f(x) = x3 - 9x? + 18x 1 = 0, x © [0,7] 1) Montrer que fx) = 0 admet trois racines réelles 0 < ay avec une précision de 10°? Comparer N avec le nombre d'iterations effectif, 3) A Vaide de la méthode du point fixe, déterminer une fonction g telle que g(x) = x aflx) = 0 pour laquelle Pitération converge plus rapidement vers la pius grande racine a3 avec une précision de 10° 4) Par la méthode de Newton, trouver une approximation de la plus petite racine a avec la méme précision, Exercice 2 Soit le sytéme tx = bet dx’ = b! ou A, beth’ sont définis comme sui. _ {12159 posed) (31781 \ ae ) \ 196418 317811 \ sine J OTL 1) Caleuler deta et inverse B= A 2) Résoudre par Cramer Ax = 6 et Ax! = b 3) Comparer b' & beta’ ax. Exercice 3 fried) (2) | 2 Sone) eee a } 1233 |b Liasa) (ad 1) Résoudre par la méthode de Gauss le systéme Ax = b en fonction de b;,b2,3,bs oli les b, sont des paramétres réels quelconques. 2) En déduire A 3) Que se passe-t-il lorsque 6 est remplacé par 'un des vecteurs de Ia base canonique de 84 e1€2.01e8 Exercice 4: Soit le systéme électrique ci-dessous, Le probleme consiste & ealculer les courants des branches, Si Pon considére que le montage csi composé de trois mailles, les courants peuvent étre calculés par la résolution d'un systéme & trois equations. Les trois équations découtent directement de la deuxiéme loi de Kirchhoff. Résoudre le systéme par ia méthode d ECOLE NATIONALE PREPARATOIRE AUX ETUDES D'INGENIORAT 3éme Année Préparatoire Année Seolaire : 2002/2003 Partiel 1 Module : Analyse Numérique Semestre:2 Date: 12/04/2003 Durée : 02h a @ @ Observation Baréme 10 06 04 Exercice t Soient f(x) =x) +3x-1 = 0, g() 1) Montrer que les équations f(x) et que a, B €}0, If 2) Montrer que les suites xn.1 = ah et Yat = a convergent respectivement vers a et B pour tous te va B43r-3=0,0ER Oct g(x) = Oadmet chacune une racine réelle a (respectivement B) dans tout R xo.¥o € (0, 1]. Calculer x6, 3) Pour xo = yo = 0.5, déterminer les nombres d’itérations suffisantes no et pour avoir [ty ~a| < 10, Va > mo et yn Bl < 10°, Yn > m 4) Expliquer la différence entre no et m; et comparer avec le nombre nz d’itérations suffisantes par Dichotomie pour avoir la méme précision. 5) Caculer une approximation de a & 10 prés par Palgorithme de Newton zp.i = zn ~ £2). commengant en S Gn) 2 = 015 et comparer avec xs Exercice 2: 1 a 1 1 Soienta € Reta © Mx(R), be R°définispard= | -¢ 202 ~a-a? -a 1 -a-a? L+a?tat lea? 1) Pour a # 0, résoudre Ax = b par Gauss, 2) Déduire det(4) et les factorisations A = LU, Lu = let A = LU, Ly = 1, 3) Sid! = (b),b2,b3)" © R%, résouder Ax’ = b’ par la factorisation.A = LU, 4) En déduire A~! = 1,2,3, La = Exereice 3 Soient f(x) = xe* 2x, x ER. 1) Montrer que la fonction fadmet un minimum unique a dans Wintervalle [0,1]. 2) Montrer que l’algorithme de Newton permet de déterminer ce minimum a. Ecrire la suite de Newton xy 3) Calculer xs sachant que xp = 1 Corrigé Exercice 1: 10 pts Soient lx) = x2 +3x~-1 = 0, g(x) = x'+ 31-3 = Ox R 1) fet g sont continues strictement croissantes sur ear f(x) = gx) = 3x2 + 1) > 0, Vee B, De plus AOYU) = giO)g(1) = -3 < 0, done il existe a. fi © R uniques tels que, fla) = Oct g(A) = 0. de plus 2,8 © 0.1L 2) Les fonction g(x) = spy et ww a vérifient fox) = 0. x= Gla), Vx € [0,1], gO) =O. x= yO). vx € [ON] ety = 3e sur[0.1) ety sont continues et décroissantes sur (0. 1], done @((0.1)) = fo(1),0) = EEF] <0. Hetwdo.ID = (v0) D’autre part Lic (0.1), ‘ 2x py) - IEE) > ecineee 910) = Gap lod = Gpayr 27 MEE aT done sp ote WD 181 = x? 180-2) 59, ve ewe & 10=£) = 0, ve (0th done sup v= WD Diaprés le théoréme du point fixe, les suites xy: = OC») et yan) = WOR) convergent respectivement vers ae B pour tous xo,y0 € [0,1]. Pour le calcul des x, voir la table | 3) fen als ets nals eenizem > Loe th De )irogt Ws) = Gage We) = xo = ‘on trouve 19 = 6. 5.21 is Poure = 10°. L = e= 10% L= 5 om trouve ni = Poure = 10%, L= $, yo = 0.5, 91 = s7pgvon trouve m = 14 Par Dichotomie ixn—ais 5 < cexigen > tog(2=# )ilog2, done MR 10°6, on trouve m » 20. y = 0.Sete= 10° onal: = ¢ 0, f"(a) = (x + 2)et > 0, Va [0.1], done, aprés le théoréme des valeurs intermédiaires et la stricte monotonic,. _f (x) = 0 admet une racine et une seule @ e]0, [. De plus f est eroissante sur (0,1), donc f(2) < f (a) = 0, Vx € [0,al et f(x) > f(@) = 0. Vx ela, I]. done f(a) = 0, fest décroissante sur {0,a{ et eroissante sur Ja, 1], 11 s’enswit que.x = a est un minimum de f sur [0, Jet il est unique. ea € [0,1] est bien défini, done it permet de déterminer a pour xo bien choisi (convergence locale). Ie, i converge vers « pour tout xo © [0,1] car f © CHO. 1), FOC) < 0, f" > Osu [0.1]. ") = (e+ Bet > 0, Wx € [0.1] ete = Dest tet que YEO mingf O/C? verte OL Pen) _ xa xe— 14 20% 2) Comme f'(x) = (s+ 2)e" + 0, ¥x © [0.1], Palgorithme xn1 = Xe sb-aeyst 3) not = 89 — iC ve + xo = 1, donne les xy de la table 3. Xn n tn ex, 0's) o 05 0 1 1 0.30769230 | 4 0.33333333 1 0,57858629 2 0.32313575 2 0.32222222 2 0.40127775 3 0.32212170 3 0,32218535 3 0.37531258 4 0.32218961 4 0,32218535 4 0,37482269 5 0,32218506 5 0,32218535 5 0.37482252 6 0.32218537 6 0.32218535 6 0.37482252 ECOLE NATIONALE PREPARATOIRE AUX ETUDES D’INGENIORAT me Année Préparatoire Année Scolaire 2002/2003 SYNTHESE Module : Analyse Numérique Semestre : 2 Date : 26/05/2003 Durée : 02h Exercice 1 Exercice 2 Exercice 3 Baréme | 07 pts 05 pts 08 pts: Exercice 1: p-4qo00 ~ ~q 0 soita=| “4 ? “9 , ot p,q € Ri tels que :p > 1, q > lp >. eae 0 0 -g p 1) Montrer que les algorithmes de Jacobi et Gauss-Seidel appliqués a Ax = 6 € R* convergent 2)Calouler J, Li, |lull2. Lill. 3) Dans le cas p = 3, q = 1, x° = (0,0,0,0)" ete = 10°, calculer le nombre d’itérations sar Jacobi et Gauss-Seidel pour avoir respectivement [jx — xl]2 < cet |jx*—xl]o < pot Ps Exercice 2: I Soit a = et © Ré quelconque. rid fez es} 14 10 20 1) Montrer que I’algorithme de Gauss-Seidel appliqué a Ax = } converge. 2) Calculer {2) = det(J— AJ) ot Jest 1a matrice de Jacobi associée a A. 3) Calculer -2) et f(-1), en déduire que Palgorithme de Jacobi appliqué & Ax = b diverge. 1 2 3 4 Exercice 3: Soit lx) = cos(ax), x € [-1,+1] et P le polynéme d’interpolation de Lagrange de faux points ~1,~1/3, 13, L sur (-1, 1. 1) Jusifier sans calcul le degré exact de P. 2) Appliquer ’algorithme de Newton pour le calcul de P (différences divisées). En déduire Fexpression de P. 3) Donner I’estimation de I'erreur d’ interpolation f(x) — P(x), x € [-1,+1]. 4) Estimer l’erreur d’ interpolation en x = 0 et comparer avec l’erreur exacte. Corrigé Synthése Analyse Numérique Exercice 1: 07pts 1) Comme p > 24, 4 est a diagonale strictement dominante, done les méthodes de Jacobi et Gauss-Seidel appliquées & Ax = b © R¥ convergent. o ¢ 0 0 0100 An ava 4 oe 0G). P 9 oro flo Gy Gy g nee o 8)" FY BY det(J' — AD = a —32? +1, done piv’) = 18 a alors p) = Zp) 4). (4) (4) 2 4((47 4% gi) a fa Wits = ($)'+(F)'+ (4) = $((F) +4 +1) <4 ig PTE 0 ° eceoo n (t= la) pe = in il> n= Mallee ps Be. i Injili Exercice 2 : 05 pts NAL =1>0, a2 = bad =1>0,As=] 1 2 3 |=1>0, A= deta =1>0, 136 done A est définie positive, la méthode de Gauss-Seidel appliquée a Ax = b € R* converge. 0 -1 -1 -10 -12 0 -32 -2 “6-112 0-5/3 120 -1/5 -12 0 3) f-2) - 3B yy = -47- _ ) A-2) = ZS = 1.3042... > 0, f-I) = ~ 3.7292... < 0, comme fest continue sur [-2,—I}, i existe Ay €] ~2,If tel que f{4o) = 0, Ao est une valeur propre de J, done p(J) > [Aol > et la méthode de Jacobi appliquée a Ax = b & R* diverge. Exercice 1: 08 pts Qs , det(J— Al) = 14 - 1) Comme fest paire et les points dinterpolation ~1,~1/3, 1/3, sont symeétriques par rapport a Vorigine, alors P(x) est pair (P(-x) = P(x), Vx € R), puisque dU P < 3, alors dU P= Oou dUP = 2, done dU P = 2 car P change de signe sur [-1, +1] (P-L)P(-1/3) < Oat PULB)PQ) < 0). fe) fod floyd HL -13,13,0) ao | 2-18 12 | 918 ! i W312 0 276 [i | 1/94! ame | 0 PR) = 1+ ROH) - Bet Note) = Th = he 3) Va € [LH 38, © [-Lo41] 20) Pt) = (+ DO Dye ea - OED Comme f(x) = x*cos(ex), alors ft) ~ Pew] s Arie? = Noe = Ah wee E141] 4) (AO) s “x = 0.45097... AO) ~ P(O) = 1 a 03125 Ecole Nationale Préparatoire aux Etudes d’Ingéniorat Analyse Numérique. 38me Année. 2002-2003 Tattrapage 06 Juillet 2003. Durée 02h Baréme : Exl (07 pts), Ex2 (08 pts), Ex3 (05 pts) Exercice 1: Soit f € O(a, 0] et H(z) le polynome de degré < 2 vérifiant : H(a) = f(a), H'(a) = f'(a), Hb) = f(b) 1) Montrer, & aide dun systéme linéaire, que H(c) existe et est unique. Soit © > 0 et L,(z) le polyndme d'interpolation de Lagrange de f en a,a+¢,b sur [a,b] 2) Determiner par la formule de Newton expression de L,(z). 3) Montrer que limL(r) = H(2). Exercice 2: xercice a Soit f € C%[-1, +1] telle que |f| < 1 sur [-1, +1} et 7 “oe, re [-1,+1}. 2 16 2 puis t= s+ 4 dans T(z) = 0, déterminer les racines x, k = 1,6 de T(z) 1) En posant ¢ telles que -1 O et A€ My(R) définie par A=} | eT 1 1 1 olte 1) Calculer []Alloo- 2) Caleuler A~ (on pourra résoudre Az = y € R*). 3) Calculer Cond,.(A). 4) Pour quelles valeurs de ¢ > 0, A est-elle mal conditionnée ? bien conditionnée ? Exercice 4: riid ' eeea, 2 Sot A=] 19 3 3 1234 1) Montrer que A est. symétrique définie positive. 2) Ecrire l’algorithme de factorisation de Cholesky. 3) Déterminer la matrice L triangulaire inférieure telle que A = LL". Bonne Chance ! Corrigé Rattrapage Analyse Numérique 2002/2003 Exercice 1: 2) H(z) = e+ yx + o92?, co, C1502 € Ry H'(e) = 61 + 2a H(a) = f(a) + ea + ena? = f(a) H(a)= fa) #) qt2o0=f'@) H(e) = $0) co + cab + cnb? = f(b) Lad\/o ie) 0 1 2 G f(a) | @ Mc= F, det M = (b- a) #0. 1b 8 0) f(b) 2) L(x) = fla) + float — a) + fla,a+¢,0)(x —a)(z-a—e) = 16) +L IOg as fered fete \'e ae) 10) flee ) slate) + LEHI LO, a4 a) ua 0, Lele) — fla) + (ole) + P= te? = 10) + Foe ~ a) +10) Sle) = 6-ODGF : = (1-27) 19 ee m0 3) Le polynome (2) = (1 F=")? Ke) +(e N6G— 2 ot vérifie bien ls conditions H(a) = f(@), Hla) = £0), HQ) = 10, Tepes 1) on a: done limbe(z) = H(e). (8) est de degré Exercice 2: 3 1 5 5 1 = x8 — Lat + — = = lor Bry ae tt = ae 1) Tle) 3 ee 3 0, t= 2? donne t! ate ii wat ot ime 2, #35 =0,done Lae, 2) x = cos8, cos(66) ©08(68) = 32.cos? 0-48 cost 0418 cos® 4— naan topo cost O-+ Ban an = 32.7 (00s). 3) We € [1,4], 36, € [-1,41]: J(@) - ps2) = The- Pa 6) =r) 6 A Done siz = 00s8, fle) ~ pala) = Tle ~ 2) = 2 cao! fr) et [f] <1 sur a ! [1,41], i vent 1 1 4) fe) Pil hg = agigg YEE LH Exercice 3: lte 1 1 1 Soit.e > 0 et AE MuR) definie par A= |} *F* tte i} CS 1 1 olte 1) [Allo = 446. -bo-1 = GHe Bte 1 “1 34e -1 -1 Ata * “AA Thoo NA =| 1-1 ate -1 | Mle care -1 -1 -1 3+e 6+e 6 =1+e. é 3) Condse(A) = [Aah the = (A+ DT yoy jonnée, pour € loin de 0 (par exemple ¢ > 1) A est bien 4) Pour e —+0*, A estelle mal condi conditionnée. Exercice 4: 1) Ay = Ag = Ag = Ay=1> 0, done Ast définie positive. as = (0x9 ~ 13, - Lh)" Lag = 28 Lgha ~ Laakso rr " Ly = (044 ~ Ly ~ Ug — Dig)? = 1 ECOLE NATIONALE PREPARATOIRE AUX ETUDES D'INGENIORAT géme Année Préparatoire Année Scolaire : 2003/2004 Partiel 1 Module : Analyse Numérique Semestre: 2 Date : 06/04/2004 Durée : 02h Baréme : 02 points par question Soit f(z) = 2x3 — 3x? -1, r ER. 1) Montrer que P’équation f(x) = 0 admet une racine réelle « et une seule dans tout R et que cette racine est localisée dans ’intervalle (3/2, 2] 2) Quel est le nombre m; @’itérations suffisantes a effectuer par la méthode de Dichotomie pour obtenir ‘une approximation de a & 10~® prés. 302 4 1\¥8 ; Soit (2) = ® ft) . ve 1/22] 3) Montrer que la fonction y et l'intervalle (3/2, 2] vérifient les conditions a) f(z) =0 x= 9(2), Ve € [8/2,2], bd) y([3/2,2]) C /2.2], ©) sup [p(@J=L<1 se/2a) En déduire que Vitération 2,41 = p(t») converge vers a pour tout sto € (3/2, 2] 4) Quel est le nombre np d’itérations suffisantes A effectuer par Vitération 2,41 = g(t) commencant en 9 = 1.75 pour obtenir une approximation de a A 10°S pros. Syn) 5) Calculer ya, 2 = 1,5, of tat =n — eas Yo = 1.75 1 ‘ F(a) On pose x= t+ 5 dans f(z) =0, ensuite t= a +b dans ’équation obtenue, avec ab = 3 1 7 6) Montrer que a? et 6° sont les solutions de P'équation s? pte 7) En déduire a, 6,¢ et la valeur exacte de a = 30 + 73+ 2V2+ ¥/3 — 2V2) = 1.677650699 3a? +1)" Soient g(a) = 62? ~83?=1, re Ret ole) = ( 7 ) ,r€ [0,1] 8) Montrer que Vitération 2.41 = (x) converge vers l'unique racine f de l'équation g(r) = 0 dans {0,1] pour tout zo € [0, 1} 9) Déterminer le nombre ns d'itérations suffisantes par la méthode de Dichotomie pour avoir une approximation de @ & 10°® prés ainsi que le nombre ng d'itérations 2,41 = (zn) commencant en 2 = 0.5 suffisantes pour atteindre la méme précision: 10) Comparer n; avec np et ns avec ng et expliquer dans quel cas la méthode de Dichotomie converge vers la racine plus vite que celle du point fixe et pourquoi ? Bonne Chance ! Corrigé Partiel Analyse Numérique 2003/2004 1) f(z) = 23 — 30? - 1, f(z) = 62(z - 1), Ve ER, f0) = -1, fl) = -2 f(-00) = ~00, f(doo) = +00, f(3/2) = -1, f(2) = 3. Done f(z) < 0, Vz < Let Jha €]l, oof: f(a) = 0 et @ €]3/2, 2 ‘au tableau des variations de f. ba 4 rine, @ = 3/2, b= 2, € = 10°, done n > 18.931... et 3)a) =z) = ( - ) Fle) =0, Vr €[3/2,2) 'b) y est strictement croissante sur [3/2,2], donc 'p ((3/2,2)) = [v (3/2), 9 (2)] = [1.570 ..., 1.866 ...] ¢ [3/2, 2). a? $1) 3 $4) >0, 262 8 (e+1)(e-1) (22 +1)" (* *) <0,2€ Ba 3)_3/8\%" Par conséquent L= sup|y'| = (2) == (3) = 0.608 009 1585... < 1. La fonction v est done 29) 2 2\31 contractante sur [3/2, 2] de constante de contraction L. D’aprés Je théoréme du point fixe, Pitération Tava = la) converge vers Punique point fixe a de y dans [3/2, 2}, qui est "unique racine de f(x) = 0 dans [3/2,2], pour tout 29 € [3/2,2] 4) One fem —al $ pol ~ nol Se > n> nO pint, L = 0608009158..., € 5 1 = 0 1078, 2 = 1.75, 21 = pao) = 1.720597 188..., done n > 22.560 ... et mz = 23. L(Yn) _ Ayn = Su +1. = eT; wy = 1.75, done Filvn) ~ Galan ~ 1)’ ‘th = 1.682539683, yp = 1.677675 278, ys = 1.677 65070, 44 = 1.677 650.699, yy = 1.677650699 3,3 1 3 Oaatts dans f(z) =0 donne t*— 7 = Oct t =a+bdonne a? +b*+ (a+b)(3ab— 7) — 5) Yat = Ua 4 34 2v2 1 1 done ab = } = a8 + 8° =2 et a? otf sont les solutions de s? ~ 304 7 = 0, done s = 4 4 4 s 1 8 7) D’aprés 6), on aa’ i (3+ 2v2) et = ; (3—2v2) (ou le contraire), done a et b= 5/3 W, done t= at bet 1 lis my _ ans tarb= 5 (14 Y542v2 + YS— 272) = 1.677 050609. dont le développement décimal coincide avec ys de la question 5). ‘3a? 4.1\ 1% 8) Comme ¢(c), Ia fonction ¥(z) a) © = (2) 4 625 — 3x? —1 = g(x) =0, We € [0,1], b) o est strictement croissante sur [0,1], done ((0, 1) = [¥(0),¥(1)} = 0.550... ,0.873...] c [0, 1) 38 >0, Ve € (0,1), verifie + wz) ay ye erne-neet rr (4) >0, Ve € [0,1], 1/3)" done L = suplv'| = (I) = 5 G) = 0.436 7902324... < 1 joa 3\2, Le théoréme du point fixe implique que T'itération zn+1 = (4) converge vers unique point de y dans (0, 1], qui est Punique racine de g(x) = 0 dans (0, 1], pour tout 20 € [0,1] 9) Par la Dichotomie, on a |on— 0) < 28 cen zme—?, In2, a=0, b=1, €= 10°, done n > 19.93156857... et ny = 20. D’autre part, -sls An =a) n> aa fiat, 1 = 04367902824..., € = 1a x = (z0) 663176 2013..., donc n > 15.183... et mg = 16, 10) En résumé, on am = 19, nz = 23, nj = 20, ny = 16. Dans le cas dena on aL ce qui explique nm < nq pout le meme zp = (a+b)/2= 1.75. Danslecas dens ona L cco qui justifie ns > m4 pour le méme 29 = (a+0)/2 = 0.5 10", zo 0.608 0.436 fixe B > 05 <05 ECOLE NATIONALE PREPARATOIRE AUX ETUDES D*INGENIORAT 3éme Année Préparatoire Année Scolaire : 2003/2004 Synthése Module: Analyse Numérique Semestre: 2. Date : 26/05/2004 Durée ; 02h Baréme : Ex.1 (04 points), Ex.2 (06 points), Ex.3 (06 points), Ex.4 (04 points) Exercice 1: rid Soit A € Ma(R) définie par A=] 5 $y | etbeRt 1123 1) Résoudte le systdme linéaire Ax = b par la méthode de Gauss, En déduire A~? 2) Montrer que les factorisations Doolittle, Crout et Cholesky de A sont identiques. 3) Montrer que 2” Ax = 01, (y(2))*, Ve € RY, of les y(x), ¢ = 1,4 sont a déterminer. Exercice 2 : 1a -g 1 Sot A={ 0 1 0 |,b=(1],a8eR Bal 1 1) Pour quelles valeurs de a, 8 € R a-ton p(J) <1? (J est la matrice de Jacobi associée & A) 2) On suppose a = f = 2et 2° = (2,2,2)". 2.1) La méthode de Jacobi est-elle convergente ? 2.2) Montrer que x* = (1,1,1)?, Vk > 1. Conelusion 3) Soit a = 0 et § = 1/2. 3.1) Déterminer le nombre d'itérations ky a effectuer par la methode de Jacobi commengant en (2,2,2)" pour avoir |[x* — lly < 10-*, Wk > ho. 3.2) Peut-on assimiler 2? & la meilleure solution approchée de x a 10-® prés ? Justifier. 4) Soit £; la matrice de Gauss-Seidel associéo a A. 4.1) Ecrire Palgorithme de Gauss-Seidel appliqué a Ax 4.2) Pour quelles valeurs de a, € Raton p(L1) <1 En déduire £y Exercice 3 : Soient A € M,(R) une matrice A diagonale strictement dominante, i.c. [ay] > Dygs lass, Yi = 1) Montrer que A est inversible. 2) Bn déduire que pour toute valeur propre \ de A, il existe i € {1,2,...,n} tel que lay A] < Dijgi las] 3) Montrer que si A est symétrique & diagonale strictement dominante et & éléments diagonaux stricte- ment positifs, alors A est définie positive. 4) En dédire que A= (ayhrcijeni Oy = O87 =i, ay =O si jf =itlet ay =Osij # iit Lest definie positive Va, 8 € R tels que a > 2{5], @ > 0. Exercice 4: Soit f € C*(R) telle qu'il existe M > 0 vérifiant : [f(«)| < M*, Vk EN, Ve © R et poni(x) le polyndme d'interpolation de Lagrange de f aux points ~n,...,—1,1,...,n sur R. Montrer que si M <2, alors lim pan1(0) = £(0). Viel Glitek ! Corrigé Synthése Analyse Numérique 2003/2004 ri tik Lili qt 0 1 -1 0 -h ‘it O01 -10 b-b on 0 = 1 AM=1 9 yo 1h fA =] 00 1 1 thm |? 00 1 2 kh 00 1 2 &-h li iib rtm +tetts = jou] Ob 0 bah mm = bab =] 00 1 1 by +b —2h mtn = byt b2— 2 0.0 0 1 bb-by—b th m= by — by — ba ty Thy — dbz — 3by + by 4-3 1 pa { Hat Bba + 2b —be 2-1 Done == | ab, 4 by + hy —by | 2 A by bby tbe 1-1-1 ; eG AHL, by, 61,234 1.2)A=1U,b=} 5 »U=) 99 mil LW, Ug=1, i= 1,2,3,4 it = 12.84 1011 oo o1 A=LIT, Ly > 0, i= 1,2, ‘Les 3 factorisations sont identiques suite leurs unicité. 1.8) oT Ax = a LET x = (L?2)"(U"2) = ||L7 23 = Doh, (le)? Ve € RY, done a(t) = 2+ 424 454-24, Yo(t) = 22 — 2a, Yo(t) = 3 + 24, Ce) = te 0 -a B 2a) =] 0 0 0 },de(s—X) =A? - 6"), B-a 0 oJ) <1 |p] Let Jacobi diverge. 2.2.2) 2! = Ja®+D-1b = (1,1, 1)" et 2? = Ja + Db = 41,1,1)", done ak = (1,1,1)", Wk 2 1 Ici la méthode converge car 2° est choisi comme étant la solution exacte x de Ax = b > 2 = Jz+c- 2.3.1) Si a = Det B = 1/2, ald) = 1/2, 2 = (2,2,9)", 21 = (21,2) fle — 2% = 1 lle = pl) = 1/2,€= 10-3, k> w(§ =e) nla) = 10.965 78428, by = 1. [2 — 22 2.3.2) Oui. 2? = (2,1,2)" est la solution exacte. 0,46, (J) = A}, done okt! = ark + Bok +1 0 -c B Bay a 2A) } Tha = gett goth pak Li={o 0 0 aktt = pak! — ark! 41 = -afrht Pug +B-a+l a 2 = (22,2) 0 -ap FF 2.4.2) det(£, —M) = *(A— f), A= 0,8, (Lr 3.1) Soit ¢ ER tel que Ax = 04 D7, ait; =0, Vi max; |x| #0 et alors — Dyitigay > alla B\, plLr) < 14 |B] 0, al vient, » FA — 04 > fase] = Eales > 0 Done A >0, VA v.p. de A qui est alors définie positive. 3.4) a > 2/6] implique Jai| > Dj ye layl, Vi = Tm et ay > 0, Vi = Tm. D'aprés 3.3) A est S.D.P. 4) V2 ER, HER: f(a) — = Te -L2o #S(@) — Pn-a(2) = TI] (@ — Pe, done ere) “ 1f(0) — Pon—1(0)| = (aE Comme |f(")| < M?" sur R, alors an 110) ~ Pana) < ty 08 69 = (nt ens _ (0D) MP? nt (mt MEP ene aT GEM 7 Done Ty wave a <1 Pour M <2 2 0 c8F ENPEI - 3éme Année Préparatoire - 2003/2004 Rattrapage Analyse Numérique - 12/06/2004 - Durée : 02h Baréme : 07 + 07 + 06 pts Exercice 1 : Soit f la fonction définie pour x > 0 par f(z) =Inz-+2? +2r—5. 1) Montrer que Péquation f(z) = 0 admet une racine a et une seule que Von localisera entre deux entiers consécutifs a et @ + 1. 2) Montrer que 'itération ta41 = 2a — f(a); Z0 € [a,a-+ I] ne converge pas vers a. 3) Pour ques valeurs de 2, Vitération zy. = Fn — Af (Zn); Zo € [2,a-+ 1) converget-elle vers a? Déterminer Ia valeur Ap de A pour laquelle max, 1 ~ 24" (z)| est minimal. —Inz, ate, 1 € fa, a+ 1} converge-telle vers a? 5) Ecrire Palgorithme de Newton appliqué & f(z) = 0. Pour quelles valeurs de 1» converge-t-il vers a? 4) Litération tay = Exercice 2 : nia Bhay ZO. Soit (z,) la suite définie par t= 1, 21 = 1) Montrer que A, = (-1)", Vn > 0. 2) Calculer x, pour n= 13, 14,15, 16. 3) Déduire la solution des systémes ee =)(2)-(%)*(s 2) (2)=(%8) 233 377) \ m )~ \ 610 233 377) \ 2) = \ 609)” 4) La matrice A de ces systémes est-elle bien conditionnée ? 5) Caleuler Cond,.(A) = JAlleol|A"leo- ), Zaza = Tn + Fass et An Exercice 3 : 2al a 1 6 ),ota,f eR vérifient a =1et bE R?. 162 1) Calculer p(X) = det(J — 1), o J est la matrice de Jacobi associée & A. 2) Montrer que p(—1) <0, et étudier la convergence de la méthode de Jacobi appliquée & Az = b. 3) Calculer la matrice £, de Gauss-Seidel associée & A, en déduite les valeurs de a, 8 pour lesquelles la méthode de Gauss-Seidel appliquée & Ax = b converge. Soit A= Bonne chance ! Corrigé Rattrapage Analyse Numérique 2003/2004 1.1) f est continue pour z > 0 et f(1) = -2 <0, f(2) = In2+3 > 0, done Ja €}1,2[: f(a) = 0. D’autre part, f"(x) = 1/z+ 2x +2>0, Vz € [1,2], donc a est unique. 1.2) giz) a =5-2-2?-Ing, p(2) = —1-2e-1/z, |p''(z) = 2-1/2? > 0, Vz € (1,2), done min} = |9"(1)| = 4 car |p est croisante, done ane — al = len) — (@)] = I¥'Elktn ~ al = Alen — a. ee eee done 2, a, n> 00. 1.3) f(a) =1/2+2r+2>0, f">0 sur [1,2] et M=max f= f(Q) = 13/2, m= mins’ = f'Q) donc max|t Ay" = max{]1 — Am], | — AM|} si A> 0. 5 (pour <0, 1— Af’ > Let 24» @ comme en 2). En effet, m051-AM S1-Af"S1-dm, de plus |L-AM| < 149 -1 < (13/2)A-1<1@00, Vee (1,2), ce Cet algorithme ey pour tout 2» € [1,2] car . fecf, 2), - f(z) £0, Vz € [1,2 {"e) = 2-1/2? > 0, Ve € [1,2 FDFQ2) <0, f'(1) =5, f"(2) = 13/2, done c= 1, od (c)| = min{|f(DL FQ} et [F0)/F(O| < b- a 2/5 <1. Le théoréme de convergence globale de l’algorithme de Newton implique que im {1,2} Tap = a — a, Vo € 21) An = taltn + tat) - Thy = ThA to ao ont) = = th tate = —Ant et Ao = tot2 — 2} = 1, done Ay = (-1)"A Aas, vn EN. 2.2) Partons de z = 1 et ), on obtient : r13 = 144, ty = 233, 245 = 377, 16 = 610. 23) {2ursente g (1H 28)(!) = (3% ) a> (1) tiesto ay ta Jer systéme car det: = -1 7 0, La solution du 2eme systéme est af = AWW = AnYb+ Ab 1)", done car U = b+ db, od 5b = (1, 7 -377 233 = 1 610) _ ( -609 (ay ‘)qeneearo-(1)+($°)=(a) 2.4) 52 = 2! — x = (-610,377)" est trés grand devant 5b = Y — b= (1,1), done A est mal conditionnée. 2.5) |[Alloo = 377, [A-tMleo = 610, donc Condao(A) = |[Alooll A™ Noo = 229970 >> 1. 0 -a/2 -1/2\ PX 3.1) & 3.2) s= i 0-6 oe 1) <0aat+p >4/2 -1/2 -B/2 0 + 6? > 208 =2> 1/2 Done Bo €] — oa If: p(%o) = 0 et alors p(J) > el > Let Jacobi diverge. ~a/2 =1/2 0 -a/2 1/2 3.3) A= /2 a/2—B 0 a2 a/2—B a -a(of -1)/4 -08/4+ 6?/24+1/4 0-0 tie det (C; — M) =X? = (a? + 6?) 02/2 4/4 A= oes pt Voth + #2) /4, done pls) = (0 + f+ Var F=2) [A < Le 2S oP +H < 5/2 jet(J — AI) = — 3 (2a? +267 +1) AFG ECOLE NATIONALE PREPARATOIRE AUX ETUDES D’INGENIORAT 3ame Année Préparatoire . Année Scolaire : 2004/2005 Partiel Module: Analyse Numérique Semestre: 1 Date : 07/12/2004 Durée : 02h Baréme : Ex.1 (10 points), Ex.2 (05 points), Bx.3 (05 points) Exercice 1: Soit f la fonction définie pour « > 0 par f(x) = In +2? + 2x —5. 1) Montrer que Péquation f(x) = 0 admet une racine a et une seule que lon localisera entre deux entiers consécutifs a et a +1 2) Montrer que l'itération zm41 = zn ~ f(2n)i Zo € [a,a + 1] ne converge pas vers a. 8) Pour quelles valeurs de 4 > 0, la fonction p(x) = x ~ Af(z) est-elle contractante sur [a, a+ 1] ? 4) Déterminer la valeur 4y de pour laquelle max, ~AY"(2| et minimal Inzy 4 A F8; 29 € [a+ 1] converge-telle vers a? 6) Eorire l'algorithme de Newton appliqué & f(x) = 0 ct montrer qu'il converge vers a pour tout 1 € [aa +1] 5) Litération 2ne1 Exercice 2: Soit y la fonction définie sur Ry par ) (x) Vr > 0; p(0) =0et a = elt-V8/2 =Ine 1) Le Théoréme du point fixe est-il applicable a la fonetion ¢ sur (0, a] ? 2) Montrer que V'itération 2,41 = g(z,) converge vers € = 0 pour tout 9 € [0,4] Exercice 3 : Soit (2,) la suite définie par Ing = Ty + Tass Wr > Oet An = tatnye~ Thar, 220. La 1) Montrer que A, = (-1)*, ¥n 2 0. 2) Les valeurs de 2, pour n = 44,45, 46, 47 sont données par 44 = 433494437, 4s = 701408733, 245 = 1134903170, x47 = 1836311903, Déduire la solution du systéme 433494437 701408733 \ (2 \ _ ( 1134903170 + ¢ 701408733 1134903170 ) \ x2 ) ~ \ 1836311903 ~ < pour ¢ = 0 et ¢ #0, Conclusion ? Bonne chance ! Corrigé Partiel Analyse Numérique 2004/2005 1.1) f est continue pour x > 0 et f(1) =-2 <0, f(2)=In2+3> 0, done Ja 1,2: fla) = Diautre part, f'(z) = te 22+2>0, Yr >0, done a est unique. 1.2) 6 done min | -1-2-1, Wre)=2-3 >0, We € (1,2), ~ f(a) =5—2~2?-Inz, ¢(2) (9"(1)] = 4 car |p| est eroissante, il vient [tees =a] = lplen) — 9a) = |v Enza ~ @| 2 dlzn ~ a, done [rq ~ a] 2 4*|20— a] —> +00 sit #0, ie. ty a, n— 00. 1.8) /@ = Laare2>o, "> Osur (1, det M= mags = 12) =F me mips = #0) =5, done max t— Af = max{]t~ Am [DAM A> 0. En effet, m< J’ M = m0; to € [1,2] F'Gq) ~ 24 2t_+ Way = | 1.6) f(z)= fet >0, Ve € [1,2], test = te Cet algorithme converge Wzo € [1,2] car fe c7,2), J'(a) 40, Ve € [1,2], S"() =2~ 3 >0,¥ee (1,2, (DF(2) < 0, FQ) =5, #2) donc ¢= 1, 08 |f'(ol = min(|/"A, 1/1} et LA < ba eo} $1. Le Théor®me de convergence F(0) globale de l’algorithme de Newton implique que lim (z,) = a, Vo € [1,2]. 2.1) ye) = me? Ve < €2, p est strictement croissante pour z < e, s #eeta<1 done (0, <1) = (40), eC@)] = (0, 7 << (dal, 0, Vz € [0,a], ¢ est strictement croissante sur (0,a], done (1 v5)/2 © (I=Ina (= v5)/2) La fonction n’étant pas contractante, le théoréme du point fixe n’est pas applicable a y sur (0, a]. a= 2.2) On distingue deux cas : 1» € [0, af et 2p = a. Pour tont 29 € [0, a1, il existe 0 < a” Lt = g(a") < ya) = Let y est L*-contractante sur [0,a"), done zny1 = @(n) converge pour Ie zp choisi dans [0, a, Maintenant, pour z = a, 1 = (to) = Tone 0, od Rest une matrice triangulaire infériroure, Partie 2: (06 = 01 + 01 + 02 + 02 pts) 2.1. Kerire les algorithmes de Jacobi et Gauss-Seidel appliqués a Az = 2.2. Montrer que ces algorithmes convergent. 2.3. Calculer J, Je» pl) et £1, (Lillo, p(Ls). 24. Pour «° = (0,0,0)" et e = 10°, calculer les nombres d'itérations ky et hy & effectuer par ces deux algorithmes pour avoir respectivement |[x* — lho Se, Wk > ky par Jacobi et |Ix* — Joo <&, Vk > ky par Gauss-Seidel. Partie 3: (05 = 02 + 02 + 01 pts) 3.1. Calculer A? et déterminer sa factorisation de Cholesky A? = ATR, of Rest une matrice triangulaire supérieure. 3.2. Montrer que si A = QR, ott R ost la matrice de 3.1, alors Q"Q =I 3.3. Calculer @. Partie 4: (04 = 02 + 01 + 01 pts) 4.1. Calculer det(A—N). En déduire Ay, Xz, As ainsi que V?, V2, V2 tels que AVE = A,V#, k= 1,2,3. 7 4.2. Déterminer une matrice U telle que U-1AU = A = diag{y, Ag, As} et U-! 4,3, Bn déduite les formes yy : R° +B telles que 2” Ax = Sh, Auyel(z), Ve € BS Bonne chance ! Corrigé Synthese Analyse Numérique 2004/2005 2-1 0 2-1 0 Uae 0 LL a={-1 3-1}, 4M=[0 52 -1),4@=[0 52 -1 -1 2 Oo -1 2 0 0 3/5 1 00 2-1 0 A=; L={ -1/72 10],u=4aM=(0 52 -1 |. 0-2/5 1 0 0 8/5 h 1 [ Sut 20+ ye 12. Le=ys2=| n/2tn etUr=z32=— | +4m+2ys | = Ay. t4/5+ 2y0/5-+ ys 8 \ tent bu ge 000 2 0 0 13. 4=ipi7; b=[-12 10),p=[052 0), 0 2/5 1 0 0 8/5 0 0-2/5 1 0 Va 0 2/3 2/25 yt = (ck +1) /2 +1) /2 2.1. {¢ sieve . { (es ane > 2° € B® donné. ( ( 1 00 v2 0 oe () 0 R=LDV =| -1/2 2) 0 VR 0 (He ‘5/2 0 0 a 1 Bad aft = (xf + 1) /2 ht = (oh 41) /2 2.2. A est a diagonale strictement dominante (2>1+0, 3>1+1, 2>041). 01/2 0 0 1/2 0 - 2.3.72 (13 0 13),a=(0 1/6 1/3 ), Wle=2 Bale tao 0 afte 176 J | Wills = 1/3, (C1) = 1/ En effet, det(J — AZ) = X° — A/3 et det(L, — A) = d° — 12/3 (ou bien p(Lx) = p*(J)). 4, Jacobi (0,0,0)", 2° = (1/2,2/3, 1/2)", lla! — 2° jn = 2/3, € = 107%, k>In (S Li) /In ||J]Jeo = 35.782, hy = 36. Gauss-Seidel : 2° = (0,0,0)", x = (1/2,5/6,11/12)", ljz" ~ 20 = 11/12, ¢ = 108, (1 = [Latheo)e k>In IRE) /intsts = 20806 , = 2 -5 1 v5 -VB vB un 4) 10m no 0 ve ue ) 55 0 0 4/3/15 1? = RTQT, donc R™R = A? = ATA = RTQTQR qui donne Q7Q =I. UVv5 Wve 1/v30 3.3. Q= ART = ( -1/V5 Ve/3_ V30/15 |. 0 -1/V¥6 v30/6 4.1. det(A — 1) = 8-70? 414-8, = 1, W=2, =A, V1 = (1,1,1)?, V? = (-1,0,1)", V3 21)". 5 1, = oe | vi yz ys _ ( v3 -1/v2 1/6 ) IV']2 = v3 EWR IBM oe BOG Je Wik U est la matrice de passage de {e*}a {V*/||V*]2}, elle vérifie U-TAU = A A= UAUT et U-1 = UT car A est symétrique et les vecteurs colonnes de U sont orthonormaux deux a deux. 4.3, o? Ar = c"UAUT 2 = (UT2)A(UT2) = yl Ay = 23, ay avec y = UT x, done Ine) = (+ 25 +)/V5, ala) = (21+ 35)/V2, @) = (Ora +25) V6 ENPEL séme Année Préparatoire. 2004/2005 Epreuve de Rattrapage d’ Analyse Numérique 13/06/2005. Durée 02h Exercice 1: Soit f Ia fonction polynomiale définie par f(z) = 120 — 600z + 6002? — 200z* + 2524-25, ER. 1) Montrer que V'équation f(z) =0 admet cing racines réelles strictement positives a1, 43,42, 4) 25: 2) Localiser ces racines dans des intervalles {a,b} du type [p,p + 1], o8 p € N. 3) Calculer des approximations de ces racines & 10° prés par la méthode de Newton b 1 sonoma = 22 dec. "En) 1-1 0 1 Soient A={ -1 2 -1], 1 ],e>0. 0 -1 l+e 1) Ecrire les algorithmes de Jacobi, Gauss-Seidel appliqués & Ax = b. 2) Calculer les matrices de Jacobi J et Gauss-Seidel £, associées & A. 3) Calculer p(J) et déduire p(C;) pour e # Let ¢ = 1. 4) Montrer que ces algorithmes convergent. 5) Dans le cas ote = 1, 2 = (0,0,0)", € pour avoir respectivement \\x* — zll2 <¢, Vk > hy par ’algorithme de Jacobi, \j2* — alla < ¢, Vk > ky par V’algorithme de Gauss-Seidel, Exercice 2 : 10-5, calculer les nombres d'itérations fi, kz & effectuer (Choisir resp. L = p(J), (C1) dans la majoration |[z* — 22 $ 2p [= |e \P> = 0 [05 [0 | 1.500000 ||-0 | 3.500000 ||"0 | 7500000 || 0 | 12.500000 1 [Oar20181 |f 7 [1.413873 || 1 | 3.600204 || 7 | 7.140002 || T | 12.650750 1.3. [2 [02389630 || 2 | 1.413403 || 2 | 3.596431 || 2 | 7.087021 || 2 | 12.640850 [3 | 0.2627561 |] 3 | 7.413403 |['3 6 || 3 | 7.085811 |[ 3 | 12.640800 4 | 0.2635594 || 4 [1.413403 |[4 | 7.085810 || 4 | 12.640800 5 | 0.2635603 J 5 | 1.413403 5 [ 7.085810 || 5 | 12.640800 | a = 0.2635603..., a2 = 1.413.403 sy 4 = 7.085810..., 5 = 12.640800 tasted $41 QL.) cht = (kaha) /2 , dat = (ot 4 ob +1)/2 , 2? € RY donne. a= (b+) (tte) Cah = (er +1) /0 +2) 0 1 0 0 1 0 2.2. c= | 1/2 0 12),4=|0 1/2 2 0 1/(lte) 0 0 1/2(e+1) 1/2(€+1) a e+? + _v3 2.3, det(J — AI) =! arcs Via AJ) = Ter Sie=1, ad) =F AL) = Pd) = 5*? sie, ples) = car A est tridiagonal. e+2 < 1ou bien A est tridiagonale symétrique définie positive, car 24. p(6s) = PD = ay Ay = Ap = Ag =1>0, Ay =e > 0, Jacobi et Gauss-Seidel sont de méme nature, done ces deux algorithmes convergent. 2.5. Jacobi : 2° = (0,0,0)", 2! = (1,1/2,1/2)", ll" - k>In (Gee) /1n p(J) = 111, 4309839, ky = 112. Gauss-Seidel : 2° = (0,0,0)", at = (1,1,1)", lz! — 2°lla = V3, (C1) =3/4, ¢= 10°, k>In (Gee) pin p(es = 54.751 79920, hy = 55. = ¥372, pd) = ENPEL. 38me Année Préparatoire. 2005/2006 Partiel d’Analyse Numérique 27/03/2006. Durée 02h Baréme: 44+ 6+6+4 Exercice 1: Soient f(x) 1 e@) =14+ yoo Fe [hl 1) Montrer que Vitération tai1 = p(t») converge vers unique racine réelle a de (2 m9 € [1,3] 2) Quel est le nombre d’itérations suffisantes a effectuer pour atteindre la précision [rq — a < 10-* sachant que 9 = 1.5 3) Comparer ce nombre avec celui qu’on devrait effectuer par la ruéthocle de Dichotomie pour attein- dre la meme précision. <2 0 pour tout Exercice 2 : Soient F(x) = x ~ f(r) =0, f(2) 1) Montrer que |/"(z)| < 5, Vr € Ret os tal < g Yn € N, of C est une co In(x? +4), ER. Ondefinit tay = (zn), 70 € R déterminer. 2) Montrer que la fonction 2 f(r) admet un point fixe unique s dans (2,3) 3) Literation zy41 = f(t) est-elle convergente vers s pour tout 29 € [ 10~ 4) Six = 2, montrer que {r, — 8] $ 5-7, Wn EN 5) Pour ¢ = 10°, déterminer ny € N° tel que [tn —s| Se, Vn 2 my 6) Pour le méme ¢, déterminer np € N* tel que [rq — 5] <2, Vn > nz par la méthode de Dichotomie La valeur de ng est-elle prévisible par rapport a celle de nm ? Exercice 3: 1) Séparer les racines de Véquation f(x) = sin (95x) + VI — z? = 0 situées dans l'intervalle (0,0.2) 2) Combien de racines l’équation g(z) = x” ~ 2r° — 32° + dr? ~ 5x + 7.54 = 0 admet-elle dan Vintervalle [-2, +2] ? 3) Montrer que l’équation h(x) = 2* — 162° + 722? ~ 962 + 24 = 0 admet 4 racines réelles positives ay < a9 < a3 < a4 que I’on localisera dans des intervalles de la forme [k,k + 1), ke N. 4) Déterminer des approximations de 1,42, 03,a4 4 10~® prés par la méthode de Newton com- mengant respectivement aux milieux des intervalles trouvés précédemment. Exercice 4: Soit x un nombre réel et fl(x) la représentation de xs en base décimale par n chifftes significatifs exacts dont p est le premier chiffre de cette représentation (A gauche). Montrer que Verreur relative due & cette représentation est inférieure & 1/ (p x 10") N.B. Les n chiffres significatifs exacts de x dans fl(x) sont ceux qui vérifient |e ~ fl(x)] $0.5 x 10- Bonne Chance ! Corrigé Partiel Analyse Numérique 2005/2006 LA. z+ f(a) = 23-2? — 2-1 est continue sur R et f(1) = -2 <0, f(2) = 1 > 0, done Ba €]1,2[: f(a) = 0. De plus, f'(z) = 322 ~ 2x — 1 = (3r+1) (x= 1) > 0, Ve €]1,2], donc a est unique. En fait, d’apris le tableau de variations de f, ona: f(z) <0, Vx < Let f est strictement croissante pour > 1, donc a est unique dans R, Pour tout # € [1,3], ona meet 2 2 2 en g(a) ratty esis fo ees 4 @(r-1P(@@+l) s=2e2-2-2-1=0, re (1,2). (0.5 pt) g(t) =14+V2(e41"", Ye) = 2 (a+ 1)”, Vz € [1,2], done y est décroissante sur [1,2]. ) (1,2) = £42), 0) = [1+ ¥278,2] 11,2). (0.5 pt) ole@l= Peery”, t= 22 ery <0, done suplel= W001 Bysster (08 pb) Diaprés a), b), ¢) msi = (en) Converge vers a pour tout zo € [1,2] en vertu du théoréme du point fixe (0.5 pt). 1.2, L=1/4, €= 10, 29 = 15, 21 = y(t) = n= [mn (C =D) jint| +1= [9.503] +1 894427191..., 0. lex — aol 21: (0.25 pt), formule de ny : (0.25 pt), 9: (0.5 pt) 1.3. a=1, b=2, m= [in | [19.931] = 19. On voit que ny = 19 > ng = 10 car € L=1/4 < 1/2. En fait, ny ~ 20 = 2ng et L = (1/2). Calcul deny: 05 pt, Comparaison : 0.5 pt 2. f(x) =In(2?-+4), f(z) =2., 2 ER, done B+4) ols te ZL ct <5 Al] <2? 446 (lx]-27 20, VWeeER. (0.5 pt) 1 P+a>2 \F(@_) — S(n—1)1 = IF" En)(En ~ tn—1)] S Biter teal Sn $ Alar ~tol, YEN, done C= zr—zol- (0.5 pt) 2.2. F est continue sur [2,3], F(2) = 2—1n8 = —0.07944 < 0, F(3) = 3 —In13 = 0.435 > 0, il existe s €)2,3[: F(s) =0 4 s = f(s). P(z) =1— f'(2) > 0 car 1/2 < f(z) < 1/2, Vr ER, done sest unique. (0.5 pt) 2.3. On applique le théoréme du point fixe a f sur le fermé I = [2,3]. f est }-contractante sur I daprés (2.1) car f est derivable et |f"| < 1/2 sur I. D’autre part F(1) = [F(2), F(3)] = [In8, In 13] = (2.079... 2.564...) CT car f est strictement croissante sur I, done J est stable par f. La convergence de zn41 = {(@x) pour tout x9 € I découle du théoréme du point fixe. Stabilité: (0.5 pt), _ Contraction : (0.5 pt), Convergence : (0.5 pt) 2.4. Ona|zg—s| << btnsr — tal \tn — tna) z ler aol, Yn EN’, 08 = 1/2, 20 a , «ry = f(o) = In8, comme fer — 70] = [In8 —2| = 0.0794... < 0.08 = 2° x 10°, il vient I=L fey no) < CPO e210 ID lta — 3] < (01 pt) i-1p 2.5. Ona |r, — 3] < lade n> Pour L = 1/2, €= 10-5, 29 =2, x =In8, on trouve n > 13.955, done m= 14. (01 pts) (Avec la majoration = ‘10°, on trouve n > 4 +31n10/In2 = 13.965, donc le méme m1). 2.6. Par la dichotomie, on a [tq —s|< sy SeN41Z nina, Comme a = 2, b= 3, € = 10-5, on trouve n > 5In 10/In2— 1 = 15.609, done nz = 16. La valeur do nz est prévue proche de m sans calcul car L=1/2. (01 pt) Pour zp = 2.5, cette comparaison est encore plus prévisible, on trouve 2; = In (10.25) et n > 15. 076, done ny = nz = 16. 3.1. f(z) = sin(95z) + VIA, = € [0,0.2), On utilise les racines de sin( 0 dans (0,0.2] pour esquisser un tracé du graphe de z ++ sin(95uz) et son intersection avec x ++ —V1— 7? (le quart du cerdle).. sin(Q52) = 0 = 952 = ha => 24 = kn /95, KEN. 2, € (0,0.2] 420 OH 05 15 45 9.5 0.2764085 | 1.772059 | 4.537162 | 9.399695 0.3204643 | 1.745897 | 4.536621 | 9.395081 0.3225431 | 1.745761 | 4.536620 | 9.395071 0.322547 | 1.745761 | 4.536620 | 9.395071 0.322547 | 1.745761 | 4.536620 | 9.395071 oeoenros 0.25 pt par racine (1 pt = 0.25 x 4) 4. Soit m le nombre de chiffres aprés la virgule de la représentation de «, p le premier chiffre cette représentation, E, Verreur absolue et B, Perreur relative. On a = 009. .Onamsba-eBmevey —— fU(E) = 002 Onan bide + ins n>1 E, <0.5x 10, done z>pxiorn™? (Apt) z>pxio-"!—05x10-, (1 pt) donc i: im . EB 0.5 x 10" 0.5, 1 1 (1 pt) sat px 10" pour toutn>1 (1 pt) Ecole Nationale Préparatoire aux Etudes d’Ingéniorat Analyse Numérique. 3éme Année. 2005/2006 Epreuve de Synthése 22 Mai 2006. Durée 02h Exercice 1: 05 pts =2+1+2 1100 1520 Soit A=1 9 2 10 3 003 7 1) Déterminer la factorisation de Doolittle A = LU, Ly = 1 2) En déduire la factorisation A= LDL", L, 3) Montrer que A est définie positive et déduire de 2) sa factorisation de Cholesky A= RR. Exercice 2: 05 pts =2+1+2 1-1 1 0 u oo 4_{-l 3-3 2 _|® Soint A=] “1 3 5 5 | 9= 1 wy 02-59 wa 1) Déterminer le systéme triangulaire A®)x = y® obtenu a la 3¢me étape d’élimination de Gauss. 2) Déterminer la factorisation A = LU, Ly =. 3) Montrer que y®) = L~!y, en déduire L* & partir de 1) Exercice 3: 05 pts =1+1+2+1 3 01 3 1 6 42-4 0} oa Soient A=] ~3 9 9 Lo 5 | eer 3 -la 3 0 1) Triangulariser le systéme Ax = b par la méthode de Gauss ordinaire 2) Pour quelles valeurs de a la matrice A admet-elle une factorisation A = LU, Ly = 1. Déterminer cette factorisation. 3) Résoudre le systéme Ax = b par cette factorisation. 4) Que représente z pour la matrice A~. Justifier ! Exercice 4: 05 pts =2+1+1+1 22-11 4 _{ 43-12 6 Soint A=| 5 5 _3 4 a 3.3 -2.2 6 La triangularisation gaussienne de A est équivalente a la multiplication de A par un certain de nombre de matrices inversibles Ex, k = 1,2,3. 1) Determiner ces matrices Ex, k = 1,2, 3. 2) En déduire la factorisation A = LU de Doolittle. 3) Résoudre le systéme Ax = b par cette factorisation. 1 Corrigé Synthése Analyse Numérique 2005/2006 Exercice 1: 05 pts =2-+1+2 1100 110 04200 042 0) = = DAM=1 9 2 0 3 aaah 1 008 000 eee _{1 1 0 *=loip 1 00 31 Beoo 2) Asymétrique > A = LU = LDL", oa D = diag{1,4,9,16}, DY? = diag{1,2,3,4}. 3) A: =1>0, Ar=2>0, As=3>0, Ay=24>0, done Aest SDP 10 0 0\/1000 roo0 £4 Uol[ozas 1200 = RRT, od R= LD = = ASHE oh Ree 012 1 0}Jo030 0130 0 0 31)/\o004 oo14 Exercice 2: 05 pis =2+1+2 T-1 1 0m 14t10 aw ai i -1 3-3 2m) qy_][0 2-2 2 ntm i 0 2 zee Ke 1-3 6-55 |'4 0-2 5-5 -mty 0 2-5 9M Dako 1-1 1 0% 1-1 10m gou{0 2-2 2ntm jou|® 2-2 2+» 0 0 3-3 mtu 00 3-3 ytw 0 0-3 7 —w-mtm 0 0 0 4 -wtmtu i 0\(n n 40 = (AM|y rye | ° 2 || 2 nth AP) = (Ay), done AMZ = y «2 | s in 0 —nt wt Teva 11000) game 1-1 1 0 [eUcel= o1-d 3) Ar sy # Le =y Ura Ly et Ac =y ¥), done y = Ly n 1000 nim | 1100 nen do 0110 ~wt is ty 1011 Exercice 3: 05 pts =1+1+241 301 31 3 g,{-§ az4aal x 0 - = DA=| 3 20 20} 4 =I 0 Sein 80 0 1 0 0 0 301 3 _{[-2 1 00 O44 2 2) aR quelconque, A=LU, =| “1 yp 1 9 | 4=lo0-1 0 1-1/4 -a 1 00 0 12 3) Ly =b > y= (1,2,0,-1/2)" et Uz =y > « = (4/3,1,0,—1)". 4) x = (4/3,1,0,—1)” est la 1ére colonne de A~! car la tere colonne de A. Exercice 4: 05 pts =2+2+1 AMz = BO ¢ EArt = Exb AM = BA Eyb 1) § AM =10 & BAM = BLO >) AD = RAY = ERA , ExEyb A@Mz = b° > B,AMr = Fgh A® = BsA® = EE, EA 0 = ByEpE\b 22 -1 1 22 -1 1 22-11 o-1 1 0 o-1 1 «0 o-1 1 0 = (2) = 0 = mato 3 1 0 Pe A?=10 0 -2 0 }A"=l0 0 -2 0 0 0 -1/2 1/2 0 0 -1/2 1/2 0 6 0 11/2 1 000 1000 10 0 0 > 2 100 Ole pee ooo B=! 4 o10]'"=(0-310]'*=|00 1 0 3/2001 Oo ot 00 -1/4 1 Si on raisonne par coefficients, i.e. on considérant ine matrice de transformation Ti, pour chaque coefficient a éliminer, il fandra déterminer 6 matrice Te, k = 1,6 telles que T3727, = Ey, TsTs = Ey, Ts = Es. Ces matrices T;, sont données par 1000 Teor0. 0 1 000 2100 0 0 1 ip) 0 te} -4 0 0 0001 0 1 0 00 0 00 0 10 0 Sooronce oune roo 2) U = A® = BPE, A => A= LU; L = Ey'By' Ey = 0 0 0 1 3) Ly = b= y= (4,-2,2,-1/2)", Uz = y x= (1,1,-1,-1)" Partiel de Rattrapage d’Analyse Numérique 5 ENPEI 2005/2006 Exercice 1: 05 pts Soicnt les dows équations f(a) =a + —L= 0, g(r) = 04 2r-1=0, © € [OT 1) Montrer que f(r) = 0 admot ane racine et une senle-a € (0.1) et gle) = 0 ade w 3€ [0,1] 2) Montrer que les itérations ty = 1/(03 + 1) et rugs = 1/( +2) convergent respectivement vers a et 3 pour tont 9 € (0,1 3) Quel est le nombre Witérationss ny a effvctuer par aryey = 1/(03 + 1s a ty 0 1) Montver que f(a} = 0 adnet trois racianes a, 4,7 strieten raciues dans des intervalles di type [pp +t), PEN 2) Quel est Ie nombre np dit i cut positives tels que 0 << 4 <7. Localiser ces ‘ations i effectuer par dichotomie pour approcher ces racines a. 10% pid 8) Soit ple) = Fe (G+ Oe — 4) > 0. Mons appartenant & Mintervalle (p,p + 1) tron que Vitération souy1 = eGrn} converge vers a pour ton ey eu 1) contenant a 1) Quel est le nombre ny Witérations & effoctuer party ga = que ty est le milion de Pitervalle {p,p+ 1) rome en 1) 5) Comparer ng et. vy ot expliquer ponrquoi le méthode de dichotomie converge phis vite. 1a) pont avoir [ry —a] < 10“, Yi S11 sachant Bxer 10 pts Soit f(r) = 0! — 20t Gr? + e+ = 0, ER. 1) Montrer que Péquation fe) 2) Lacaliser y.a2.ergoetg pris en donner des valeurs approchées & 10-5 prés par Valgorithme de Newton, 3) Bu posant r= 1+ 1/2 dans f(r) = 0, montrer que cetic Gquation se réduit a f+ af? & b+ © = 0, af a,b,c € B sont a calenler 4) Bn utilisant la factorisation (4 af? + bt-be = (2 + att dU? + M+ ji), déterminer le systénie nonelingaire winili6 par 0%, 8,54 5) Montrer que A est solution «nue équation polynomiale de dogré 6 8 trouver 6) En posant © = A%, mourer quan se ( audiet quatre vacines eGelles onary, 0r4, 044 sue Amie équtation polynomiale duu 3edme degté. Montrer que les racies exaetes ay.oz,q.e4 sont données. 8 Forclte pres. par 14 V5 + 104275) = 3.520 47021 og = B+ v5 ~ V+ 2V5) = —0.284 0790438 = 0,557 5365158 -1. 793604 193 Bonue Chance ! Ecole Nationale Préparatoire aux Etudes d’Ingéniorat ‘Analyse Numérique. 3¢me Année. 2005/2006 Epreuve de Rattrapage Exercice 1: 05 pts =2-+1-+2 4222 2175 5 25 38 8 25 8 67 1) Déterminer la factorisation de Crout A= LU, Us =1 Soit A= 2) En déduire la factorisation A = UT DU, 1 3) Montrer a partir de 2) que 2" Az = Sct, al Yy2(2), Vz € RY, od af”, k = 1,2,3,4 sont les pivots de Gauss et x(x), k = 1,2,3,4 sont des formes linéaires indépendantes & déterminer. Exercice 2: 05 pts= 214141 1 1 ¥2 1/3 by aae{ 1 2 1 28] ,_[ Sot A=] jp 1 aja a/3 Jom] by | ob arPrbnbobv eR 1/3 2/3 af 1 by 1) Triangulariser le systéme Az = b par la méthode de Gauss ordinaire et indiquer pour quelles valeurs de a cette triangularisation est possible. 2) En déduire la factorisation A= LU, Lic = 1. 3) Pour quelles valeurs de ar la matrice A est-elle définie positive. 4) Donner alors sa décomposition de Cholesky A = RR? Exercice 3 : 04 pts +2 lte 1 1 Soit ees €M,(R), € > 0. 1 o 1 olte 1) Résoudre Az = y € R°. En déduire A. 2) Calculer Condgo(A) = |}Allool]A~"lfoo- Exercice 4: 06 pts=1+24+1+4+141 2 acl soit A= | 5 0 1) Montrer que A est symétrique définie positive. 2) Déterminer les valeurs propres 4;, 42,3, A4 de A et les vecteurs propres associés Vj, V2, Va, V4. Indication : det(A— I) = At — 12° + 47X* — 72 + 36 et Ay = 1, Ap = 2. 3) Si A = diag {X, 2,42, As}, déterminer P € M,(R) inversible telle que A = PAP. 4) En déduire Q € M,(R) inversible telle Qt = QF et A= QAQ™. 5) Quelles sont alors les formes linéaires indépendantes ¢,(z), k = 1,2,3,4 vérifiant a? Az = D4 dvi(z), Ve € Rt Bonne Chance ! Corrigé Rattrapage Analyse Numérique 2005/2006 Exerci pts= 24142 2 42°22 aon ae 5 oi 4 4 ow 4 4 = 0 = 0) = ae sa o4ar7 | 00 36 6 |? 67 04 7 66 00 6 65 22 1 0 00 4222 44 12 1 00 ow 44 = ae 1 36 | PATLUL=) a2 ay 1 of} Y=] 0 0 36 6 0 64 12 1/4 1/6 1 00 0 64 Ml vient : IP =UTLT = LU; UT, U=L", Us ) A= LU, Us = Let L = UD implique A = UT DU; D = diag {4, 16,36, 64) = diag {ars aS) 08), af) } nal (Uz) = ty. ale » ealz)=met tt Exercice 2: 05 pts =2+1+141 "Dg; g = Uz. Done Ay3(x) + 1693 (2) + 36p3(z) + 64y3(z), Ve € RS, » 9(2) = 25+ Les formes linéaires y,(z), k = 1,2,3, 4 sont indépendantes, car y,(2) = Be pati = te , k= 1,2,3,4 > 2 = 0. 1 1 12 1/3 by Be Piece leet 2/3 by Sot A=| 179 1 a/2 a/3 ma b » otk @, by, ba, ba, bg € R. 1/3 2/3 a/3 1 by igieeest/2 13 by o1 1/2 1/3 by + be 9p =) YAM = 09) 9g 12 (a-1)/3 |7=| wrth |°%* 00 0 —2a/9+1 ~2by/3 + ba gee teri 2 01/3 1 0 0 0\/i1 172 1/3 y{ 4 2 2 28 1 1 0 o}fo1 12 1/3 a 1/2 1 a/2 af3 y2 12 1 0f{ 00 @-/2 (@-HA |? 1/3. 2/3 a/3 1 173 1/3 2/31) \00 0 —2a/9+1 3) (1) /2 > et —2a/9-+1 > 0, done « €}1,9/2I- 4) DY? = diag {1,1, V@=N72, /@=2a5 73} 1 0 0 o raat 0 0 = RRM; R= LDV? = eee ee 12 1/2 J@-n-p 0 1/3 1/3 $V2(@—1) V@-2a)/3 Exercice 3 : 04 pts = 2 +2 (l4e)t tagta3t---+ n= (l+e)aptaet2st-+tn =” 1) Are yor| BEUH Dz bast tae wen +t = hh ay baat eytee + (146) n= te ey + £4 = Yn — (te ntmtyt tia =) Rant e-wd/e gi = (Lb e) 2+ rg bay teeta = (1+) zt (n= Dart (2 baste tn (= Dm) ley done (n+ e)ar = (21 +e) ays =~ ta) /E et alors wn) [e(r+2)) a) /(e(n+e)) ua) /(e(n+e)) - Done (n-14+e)m— me — tht (n—1+2)ye— ys — —n~wt(n-1l+e)ys— Beem — mamta na $2) mn) (el +2) pelbe “boos ad 1 n-lte i Av =y@2= Ay avec A) = way y ave AY = Tare) ’ 3 ay _2-V +e 2) [lAlleo =e, JATfoo = “Ba I*S = = ) [Allo =e, HAMIL (nt) => Condo(A) = 1+ Exercice 4: 06 pts =14+2+1+1+1 1) Aeest symétrique, A; = 2>0, Ar =7>0, Ay =24>0, Ay = 36> 0, done A est S.D.P 2) det(A— AM) = M— 1229-4470? — 72+ 36 = (A 1) (A—2) (97 - 92 + 18), donc Az = 3, A 2 1 0 1 -1 0 1 _|2 f=] 0} a)" 2 2 -1 0 1 1000 2101 _[o200 _ _{-1 0 1 2) \eipe a= Po DA=| 903.0 [Pa MeVeVeMd= 1 “1 9 1 2 | Mle A= PUAP 0006 2 01 7 ee yva a Bae -1/vi 0 1/v2_ 2/v10 ar Mtoe Ome /2) 2/00) cae avo -1/v2 0 1/vi0 14Q et Q? = Q7 car les vecteurs colonnes de Q sont orthonormaux 2 a 2. AQT => a? At = aT QAQT2 = (QP2)" A (Q?x) = VTAY = Dhar MUR), VEE RY A V2 Va Vs 00= i 0 O= [Taran Tae We Lam), Y= g(a H)» (22+ 29), a= L Wy = ee (tr + 2g + 2s + 44). Les formes linéaires jy(x), k= 1,2,3,4 sont indépendantes, car (2) =0, k= 1,2,3,4 > 7 =0.

Vous aimerez peut-être aussi