Analyse numérique élémentaire 2007 Génie Mécanique et Conception Université de Technologie de Belfort Montbéliard

Jmartin

Description

Examen du Supérieur Université de Technologie de Belfort Montbéliard. Sujet de Analyse numérique élémentaire 2007. Retrouvez le corrigé Analyse numérique élémentaire 2007 sur Bankexam.fr.

Sujets

• 2007

Informations

• Publié par: bankexam
• Publié le: 29 juin 2008
• Nombre de lectures: 32
• Langue: Français

Extrait

Utbm mt40Examen finalAutomne 2007 Nb: Onrédigera les exercices sur des feuilles séparées. On pourra admettre tout résultat intermédiaire afin de poursuivre la résolution d’un exercice.

Exercice 1 Le but de cet exercice est d’approximer la fonction logarithme népérien, notéeln, par différentes méthodes numériques. 1.1Polynôme d’interpolation a)Donner la forme de Newton du polynôme d’interpolation de la fonctionlnsur le support{1,2}. b)On considère la fonction erreure 1 , définie pare 1 (x) = ln(x)−p(x). Étudier les variations dee 1 sur[1,2]. En déduire que l’erreur est maximale, en valeur absolue, pourx= 1/ln(2); donner la valeur maximale dee 1 sur[1,2]. 1.2Méthode du point milieu Pour la suite du problème on rappelle que la fonctionlnest définie par:  x 1 ln(x) =dt t 1 Nb :On remarquera ici que la variable d’intégration esttet nonx; la variablex, elle, définit le segment d’intégration ! a)Montrer en utilisant la méthode d’intégration du point milieu, que :  x 1x−1 dt≃2 t x+ 1 1

b)Étudier les variations de la fonction erreure M , définie sur[1,2]pour cette méthode par : x−1 e M (x) = ln(x)−2. x+ 1

c)En déduire la valeur maximale dee M sur[1,2]. 1.3Intégration gaussienne  x 1 Maintenant on cherche à évaluerln(x) =dtpar une intégration de Gauss-Legendre à deux points. 1 t a)Montrer en précisant le changement de variable à réaliser que :   x1 1x−1 dt=du t(x−1)u+x+ 1 1−1

b)En déduire, en majorant l’erreur de méthodee G commise dans l’évaluation deI, que cette méthode −2 donne une approximation deln(x)à10près sur[1,2].

1

.../...

1.4Discussion

a)Remarques générales :

•La différence de précision entre les méthodes1.2et1.3est-elle surprenante ? •Au vu des calculs effectués, par quel type d’objets mathématiques a-t-on approchéln(x)dans les parties1.2et1.3? Même question pour la partie1.1.

b)Généralisation

•Si on veut une précision d’ordreεen utilisant la méthode de la partie1.3, comment opérer ? Fournir le plan de mise en oeuvre.

Exercice 2 L’objet de l’exercice est l’étude de la méthode de Müller, extension de la méthode de la sécante utilisée dans le cadre de la résolution d’équations non linéaires. Soitfune fonction numérique qu’on suppose définie surRSoitpour des raisons de simplicité.iun entier naturel supérieur ou égal à2; on considère trois réels notésx i ,x i−1 etx i−2 . Onfournit en page 4 une représentation graphique générique de la situation où figurentfet son polynôme d’interpolation sur le support{x i , x i−1 , x i−2 }considère l’équation. On(E) :f(x) = 0et on notella solution cherchée.

2.1Etude mathématique préalable a)Interpolation Fournir l’écriture de Newton du polynôme interpolateurpdefsur le support{x i , x i−1 , x i−2 }en fonction des différences divisées defet des points du support. b)Transformation de l’écriture dep •Montrer que pour toutxréel on a: 2 (x−x) (x−x) = (x−x() +x−x) (x−x). i i−1i ii i−1

•En déduire que pour toutxréel on peut écrire:

sachant que

2 p(x) =f(x i ) + (x−x i )c i +f[x i , x i−1 , x i−2 ] (x−x i )

c i =f[x i , x i−1 ] +f[x i , x i−1 , x i−2 ] (x i −x i−1 ). .../...

2

′ ′ c)Soit(E)l’équation approchée de(E), définie par:(E) :p(x) = 0. ′ •Combien l’équation(E)admet-elle de solutions dansCen général ?Existe-t-il des exceptions ? •On poseX=x−x i considère l’équation. On(E”) :P(X) = 0, oùP(X)est donnée par: 2 P(X) =f(x i ) +c i X+f[x i , x i−1 , x i−2 ]X Montrer que si l’on souhaite construire par ce procédé une suite(x i )dont on espère qu’elle i∈N converge versl, on devra : ′ —définirx i+1 comme la racine convenable de l’équation(E); —et par conséquent poser x i+1 =x i +X 1 , oùX 1 désigne la solution de(E”)de plus petit module. 2.2Algorithme de Müller Un certain nombre de sous modules de l’algorithme proposé ci-dessous, sont seulement nommés, sans être précisément décritsUne fois complété, cet algorithme pourra être implémenté sous langage et opérationnel! On demande, grâce à l’étude mathématique préalable, de fournir le détail des modules suivants: Détermination de(E”),Calcul deX 1 etDéfinition devect(i+ 1) en identifiant nettement les outils et l’environnement nécessaires. *********************** mu¨ller(x 2 , x 1 , x 0 , nmax→vect) 1. En-tête •Entrées: —x 2 , x 1 , x 0 valeurs initiales de la racine cherchée; —nmaxd’itérés de Müller demandés (: nombrenmax≥3) •desSorties: Vecteurnmaxpremiers itérés de Müller. 2. Corpsd’algorithme (a)Initialisations: •{déclaration devectcomme vecteur denmaxcomplexes; •vect(1)←x 0 ;vect(2)←x 1 ;vect(3)←x 2 } (b)faire pouri←3jusqu’ànmax−1 1.Détermination de(E”); 2.Calcul deX 1 ; 3.Définition devect(i+ 1). fin de faire eni 3. Finde fonctionmu¨ller( )

***********************

2.3On suppose ici que la fonctionfest polynôme.

a)La méthode de Newton permet-elle de déterminer des racines multiples, des racines complexes? Pourquoi?

b)Mêmes questions pour l’algorithme de Müller.

3

Figure 1:

4