UNIVERSITE DE MAROUA Ect ‘OLE NATIONALE SUPERIEURE POLYTECHNIQUE DE MAROUA DEPARTEMENT DE GENIE CIVIL ET ARCHITECTURE & THE UNIVERSITY OF MAROUA RATTRAPAGE DES METHODES NUMERIQUES IC3 ; semestre : 5 ; session EXERCICE 1 : 04 points (questions de cours) 1 Matlab. Solve the linear system of equations by using Matlab: durée: 02 Heures Dites ce que font les commandes : quit ; clear all ; ele ; help dans environnement Ipt xty+2z=1 x-y-3z=0 Ipt -2x-Syt2=4 3. Donner les relations qui permettent d’approximer numériquement une intégrale 1=[f(x}de et d’évaluer Perreur r—frcra| par la méthode des Trapézes. Ipt 4. Ecrire les schémas numériques d’Euler expli et Runge-Kutta d’ordre 4 qui permet de résoudre numériquement une équation différentielle du type y'= f (ty). Ipt EXERCICE 2: 14,5 points (Dosage colorimétrique du manganése dans un acier. On se propose de tracer la courbe d'absorbance A=f(A)a partir des résultats expérimentaux consignés dans le tableau ci-dessous. (A(nm)) A [400 [410 [420 [430 440 450] 460 [470 [480 [490 [500 [ | | 0,073 | 0,053 | 0,046 | 0,048 | 0,059 A Ae 0,216 | 0,313 | 0,397 [510 [520 [530 [540 [550 [560 [570 580 [£0,507 | 0,569 [ 0,638 | 0,581 | 0,606 | 0,408 | 0,358 | 0,224 I- Créer les tableaux Aet A (0,75x2)=1,5pt 2- Tracer A=/f (A) (afficher les points expérimentaux uniquement) 1,5pt 3. Tracer la courbe par interpolation avec un pas de Inm. Pour cela utiliser la commande «spline cubic » pour l’interpolation et « hold on » pour afficher les deux courbes sur le méme graphique. (0,5x3)=1,5pt Scanné avec CamScanner4- Déterminer A pour A quelconque (a aide de la spline cubic ). EXERCICE 3 1,5 points I~ Soit "équation différentielle; “42y=e" avec la condition initiale y(0)=1 et dt OSx<5 (h=0.1). 1.1- Déterminer la solution exacte de cette équation différentielle. 0,5pt 1.2. Eerire un code Matlab qui résout cette équation par le schéma numérique d’Euler explic 1,5pt 1.3- Tracer sur le méme graphe Ia solution exacte et la solution numérique donnée par Euler explicite, (Metre en valeur le graphe a I’aide d'un titre et d’informations portées sur l'abscisse et ’ordonnée ainsi que la légende). 1,5pt 2+ Un projectile de masse m, propulsé dans I'air a une vitesse initiale ¥, sous un angle de tir &. Ses coordonnées (x(0), »(@)) au cours du temps vérifient les équations différentielles x()+kxe()=0 yOrtky@=-g x) =0 >} yO=0 et re [0:10] x=, cosa yO=nsine Les solutions exactes sont : j vo sin ar -t 10=(2E 4 Ee fr 21- Eerire deux sous-programmes foncl.m et fonc2.m qui définissent ces équations différentielles. 2pts 2.2. Berire des codes Matlab qui pour résoudre ces équations par le schéma numérique @°Buler explicite et le schéma numérique de Runge-Kutta d’ordre 4. (1,5x2) =3pts 23- Tracer sur le méme graphe les solutions exactes et les solutions obtenues numériquement par les méthodes d’Buler explicite et Runge-Kutta @ordre 4, (Metre en valeur le graphe a l'aide d'un titre et d'informations portées sur Pabscisse et 'ordonnée ainsi que la Iégende) 3pts Prendre : v=Sm/s ; g=10m/s* ; k=20 et a=30° Scanné avec CamScannerUNIVERSITE DE MAROUA THE UNIVERSITY OF MAROUA ECOLE NATIONALE sup; ‘OLE NATIONALE SUPERIEURE POLYTECHNIQUE DE MAROUA DEPARTEMENT DE GENIE CIVIL ET ARCHITECTURE CONTROLE CONTINU DES METHODES NUMERIQUES 1C3 ; semestre : 8 ; session : 1 ; durée : 02 Heures ERCICE 1 : 05 points (questions de cours) 1. Dites ce que font les commandes : quit ; clear all ; cle ; help dans l'environnement Matlab Ipt 2. Use the command roots to calculate the solutions of the equation: -2x+2° +47 =4 0,5pt x+y+2 3. Solve the linear system of equations by using Matlab: } x y'—3z Ipt -2x-Sy+ 4. Donner les relations qui permettent d’approximer numériquement une intégrale J = J f(x} et d’évaluer Perveur {7 ~f (shi par la méthode des Trapézes 1,25pt : | cite ct Runge-Kutta d'ordre 4 qui permet de 5. Berire les sehémas numérique: xpi résoudre numériquement unc ¢quation différentielle du type y'= f (ty) 1,25pt EXERCICE 2 : 07 points (Mise en valeur des résultats expérimentaux) (On considére ici un ressort que I’on tend plus ou moins tout en mesurant la force qui lui est appliquée 4 Paide d’un dynamoméue. Les résultats que l'on a obtenus sont les suivants : Longueurfen Cm] [42 150 [60 [7.0 [80 [9.0 [100 {11.0 | 12.0 13.0 [14.0 Forcefen N] 00 11.1120 [32139 146 [58 [70 |83 [90 [95] T- Créer les tableaux Long (correspondant a Ja longueur) et Force. 1,5pt ¢ eu fonction de 1a longueur (afficher les points expérimentaux 2pts 2+ Tracer le graphe de la fore: uniquement). 3+ Mettre le graphe en vi Vordonnée. nomiale modélisant au mieux cette caractéristique ; celle d’un ressort est 4- Rechercher une loi poly: généralement lingaire, éventuellement cubique. Afficher les courbes sur le meme woke as Spts. mettre la légende, aleur a ’aide d'un titre et d’imformations portées sur 'abscisse et Ipt (cs GCA wt a aes Scanné avec CamScannerEXERCICE 3 : 08 points E oe x 1" Soit l'équation différentielle : art ay set avec la condition initiale y(0)=1 et 0<1<10 ( h=0.1). Déterminer la solution exacte de cette équation différentielle. 2pts Ecrire un code Matlab qui résout cette équation par le schéma numérique d’Euler explicite (1,5x2) =3pts et le schéma numérique de Runge-Kutta d’ordre 2. Tracer sur le méme graphe la solution exacte et les solutions obtenues numériquement par lee méthodes d’Euler explicite et Runge-Kutta d'ordre 2. (Mettre en valeur le graphe & Vordonnée ainsi que la légende) Laide d'un titre et d’informations portées sur I’abscisse et 3pts Scanné avec CamScanner— as UNIVERSITE DE MAROUA fod THE UNIVERSITY OF MAROUA ECOLE NATIONALE SUPERIEURE POLYTECHNIQUE DE MAROUA DEPARTEMENT DE GENIE CIVIL ET ARCHITECTURE CONTROLE CONTINU DE METHODES NUMERI QUES EXERCICE 2: 05 points Soit Méquation uifférentictle It) y avec Ia condition initiate y(0)=1 et O et d’évaluer lerreur js -Jr( “94 par la méthode de Simpson ‘tpt EXERCICE 3: 05 points 1. Créer le fichier f1.m correspondant la fonction X> fi(x)=exp-x'+1) ape Scanné avec CamScannerDonner le graphe sur Mintervalle [0:1] et dire quette est ta méthode d ‘atégration A humeérique adaptée pour le calcul approché de: | = Jexp(-x! +1)dx. (140,5)=1,5pt Ecrire ta formule de Simpson qui permet dapproximer | d 1s un intervalle [a;b] puis calculer Ta valeur approchée 1, de 1 par la méthode de Simpson lorsqu’on utilise 1a subdivision(0 1/10 2/10 3/10 4/10 S/O 6/10 7/10 sO 9/0 1). (141,5)=2,5pts. EXERCICE 2: 05 points Dosage colorimétrique du manganése dans un acier. On se propose de tracer Ia courbe dabsorbance A= f (2) a partir des résultats expérimentaux consignés dans le tableau ci-dessous. ( A(nm)) A474 ]4]4]4 ]4]4 ]4]4]4]5 15 ]2 > 11> 11> o}1 ]2/3 ]4]s Jo }7]8 }9 Jo |i 2 13 Ja |s Jo |7 [8 o |0 |o |o {o Jo Jo Jo Jo Jo Jo jo Jo Jo Jo Jo Jo Jo Jo [0 [0.70.70 70.0. fo. Jo. fo. Jo. Jo. fo fo. 70 fo Jo fo fo. fo. Ao }o jo Jo jo jo {1 J 1 ]2 ]3 43 ]5 |s Jo Js Jo J4 ]3 2 7 |5 |4]4 {5 |7]o /s5 |i ]1 [9 Jo 6 ]3 [8 Jo Jo {5 2 3 13 16 |8 |9 [5 [4 |9 J6 [3 |7 |7 19 |8 Jt Jo [8 18 14 1- Créer les tableaux Aet A. (0.5x2)=1pt 2- Tracer A= f (A) (afficher les points expérimentaux uniquement) 1,5pt 3- Tracer la courbe par interpolation avec un pas de Inm. Pour cela utiliser la commande « spline cubic » pour interpolation et « hold on » pour afficher les deux courbes sur Je méme graphique. (0,25+0,25+0,5)=1pt Déterminer Apour A quelconque (a l'aide de la spline cubic ) O,5pt - Scanné avec CamScannerUNIVERSITE DE MAROUA. he see THE UNIVERSITY OF MAROUA ECOLE NATIONALE SUPERTEURE POLYTECHNIQUE DE MAROUA DEPARTEMENT DE GENIE CIVIL ET ARCHITECTURE EXAMEN DES METHODES NUMERIQUES IC3 ; semestre : 5 ; session : 1 ; durée : 03 Heures EXERCICE 1 : 04 points (questions de cours) “1, Dites ce que font les commandes : quit ; clear all ; ele ; help dans environnement Matlab. Ipt xty+2z=1 0 Apt 2. Solve the linear system of equations by using Matlab: )_x-y—3: -2x-Sy +z 5 3. Donner les relations qui permettent d’approximer numériquement une intégrale refs :jds et d?évaluer erreur rf f(x) | narta méthode des Trapézes. Ipt 4. Bette ls schémas mumériques Euler explicite et Runge-Kutta ordre 4 qui permet de résoudre numériquement une équation différentielle du type y L(y). Ipe EXERCICE 2 : 04,5 points (Dosage colorimétrique du mang dans un acier. On se propose de tracer la courbe d’absorbance A=f(A)a partir des résultats expérimentaux consignés dans le tableau ci-dessous. (A (nn) ) A \|400_[410_[420 [430 [430 [450 [460 [470 [480 | 490 | 500 A 9,073 | 0,053 | 0,046 | 0,048 | 0,059 | 0,075 | 0,104 | 0,159 | 0,216 | 0,313 | 0,397 510_| 520 [530 [540 [550 [560 [570 [580 9,507 | 0,569 | 0,638 | 0,581 | 0,606 [0,408 | 0,358 | 0,224 1+ Créer les tableaux det 4, 2 Tracer A=f(2) eee (alficher les points expérimentaux uniquement) 1,5pt % Tracer la courbe par interpolation avec un pas de Tum. Pour eela utiliser commande « spline cubic »; i ® pour Vinterpolati eee polation et « hold on » pour afficher 1 la les (0,5x3)=1,5pt ——— ———__ rr, EXAMEN Méthodes Nunériques les Gan ‘2on0-2021 Paget 4 ( Mf Scanné avec CamScannerUNIVERSITE DE MAROUA Saf TIE UNIVERSITY OF MAROUA ECOLE NATIONALE SUPERIEURE POLYTECHNIQUE DE MAROUA DEPARTEMENT DE GENIE CIVIL ET ARCHITECTURE EXAMEN DES METHODES NUMERIQUES 1C3 ; semestre : 5; session : 1 ; durée : 03 Heures EXERCICE 1 : 04 points (questions de cours) 1. Dites ce que font les commandes : quit ; clear all ; ele ; help dans l'environnement Matlab, Ipt xtyt22= 2. Solve the linear system of equations by using Matlab: | x—y~3z Ipt -2x-Sytz=4 + 3. Donner les relations qui permettent dapproximer numériquement une intégrale . free} par la méthode des Trapézes. Ipt 4, Ecrire les schémas numériques d’Buler explicite et Runge-Kutta d’ordre 4 qui permet S (ty). Ipt de résoudre numériquement une équation différentielle du type y'= EXERCICE 2 : 04,5 points (Dosage colorimétrique du manganase dans un acier. On se propose de tracer la courbe dabsorbance A=/(A)a partir des résultats ‘expérimentaux consignés dans Ie tableau ci-dessous. (A(nnt)) 400 [410 [420 [430 [440 [450 [460 [470 [480 [490 [500 0,073 | 0,053 | 0,046 | 0,048 | 0,059 | 0,075 | 0,104 | 0,159 | 0,216 | 0,313 | 0,397 A 3510 [520 [530 [540 [550 [560 [570 [580 0,507 | 0,569 | 0,638 | 0,581 | 0,606 | 0,408 | 0,358 _| 0,224 1- Créer les tableaux Aet A. (0,75x2)=1,5pt 2- Tracer A= f(A)(afficher les points expérimentaux uniquement) L Spt 3- Tracer la courbe par interpolation avec un pas de Inm, Pour cela utiliser Ia commande « spline cubic » pour l’interpolation et « hold on » pour afficher les deux courbes sur le méme graphique. (0,5x3)=1,5pt Scanné avec CamScanner