1

Université Mohammed Premier Année universitaire 2019-2020

Faculté des Sciences Filières Master «ME – PMR – OM »
Département de Physique Module « Calcul scientifique »
Oujda

Devoir à rendre

- Le devoir doit être rendu avant le lundi 27 juillet 2020 à 12h sur l’adresse e-mail suivante :
[email protected]
- Le devoir doit être rédigé à la main avec votre écriture personnelle, je n’accepte pas des
copies saisies sur ordinateur.

I. Analyse numérique linéaire :

1. Soit 𝐷 = (𝑑 ) ,…, une matrice diagonale d’ordre n > 0. Donner une condition nécessaire et
suffisante pour que D soit inversible.
Donner un algorithme qui permet de résoudre l’équation d’inconnue 𝑥 : 𝐷. 𝑥 = 𝑏, 𝑥 étant une
matrice uni-colonne donnée.
2. Soit 𝐿 = (𝑙 ) , ,…, une matrice triangulaire inférieure d’ordre n > 0. Donner une condition
nécessaire et suffisante pour que L soit inversible.
Donner un algorithme qui permet de résoudre l’équation d’inconnue 𝑦 : 𝐷. 𝑦 = 𝑏, 𝑦 étant une
matrice uni-colonne donnée.
3. Soit 𝑈 = (𝑢 ) , ,…, une matrice triangulaire supérieure d’ordre n > 0. Donner une condition
nécessaire et suffisante pour que U soit inversible.
Donner un algorithme qui permet de résoudre l’équation d’inconnue 𝑥 : 𝑈. 𝑥 = 𝑏, 𝑥 étant une
matrice uni-colonne donnée.
4. Normes vectorielles :
Soit un entier 𝑛 > 0 .
i. Montrer que les applications suivantes définies sur ℝ sont des normes de ℝ ∶
𝑥 ⟼ ||𝑥|| = ∑ |𝑥 |
/
𝑥 ⟼ ||𝑥|| = (∑ |𝑥 | )
𝑥 ⟼ ||𝑥||∞ = max |𝑥 |
,…,

ii. Montrer les relations suivantes sur ℝ :

||𝑥||∞ ≤ ||𝑥|| ≤ 𝑛 ||𝑥||∞
||𝑥||∞ ≤ ||𝑥|| ≤ √𝑛 ||𝑥||∞
||𝑥|| ≤ ||𝑥|| ≤ √𝑛 ||𝑥||
On dit que les trois normes sont équivalentes.
5. Conditionnement d’une matrice :
 2

Soit 𝐴 une matrice carrée d’ordre 𝑛. On définit les normes matricielles induites par les normes
vectorielles ci-dessus par :
||𝐴|| = max ,…, ∑ |𝑎 |
||𝐴|| = max ,…, ∑ |𝑎 |

||𝐴|| = 𝜌(𝐴 𝐴)
Pour 𝑝 = 1, 2, ∞ , on définit le conditionnement de 𝐴 par Cond (𝐴) = ||𝐴 || ||𝐴|| .
Montrer les propriétés suivantes :
i. Cond (𝐴) ≥ 1
ii. Cond (𝛼𝐴) = ||𝐴||
iii. Si 𝐴 est symétrique réelle alors :
| ( )|
Cond (𝐴) = | ( )|
, où 𝜆 (𝐴) est une valeur propre de 𝐴.

iv. Si 𝑈 est une matrice orthogonale (c.-à-d. 𝑈 = 𝑈 ) alors :

Cond (𝑈) =1 et Cond (𝐴𝑈) = Cond (𝑈𝐴) = Cond (𝐴)
6. Analyse des perturbations :
Soient 𝑥 et 𝑥 + 𝛿𝑥 les solutions des systèmes linéaires : 𝐴𝑥 = 𝑏 et (𝐴 + 𝛿𝐴)(𝑥 + 𝛿𝑥) = 𝑏
|| || || ||
i. Montrer que : || ||
≤ Cond (𝐴) || ||

ii. En déduire que si ||𝐴 𝛿𝐴|| < 1 alors :

||𝛿𝑥|| ||𝛿𝐴|| 1
≤ Cond (𝐴)
||𝑥|| ||𝐴|| 1 − ||𝐴 𝛿𝐴||
7. Méthodes itératives :

Soit 𝐴 ∈ ℳ (ℝ) une matrice carrées 𝑛 x 𝑛 inversible telle que ses éléments diagonaux soient tous
non nuls et soit 𝑏 un vecteur de ℝ . On souhaite résoudre le système linéaire 𝐴𝑥 = 𝑏 en utilisant la
méthode itérative suivante :

𝛼 étant un réel non nul et le vecteur 𝑥 ∈ ℝ étant donné, on construit la suite (𝑥 ) ∈ℕ par la
formule de récurrence :

𝑥 = (𝐼 − 𝛼𝐷 𝐴)𝑥 + 𝛼𝐷 𝑏 (1)

où 𝐼 est la matrice identité et 𝐷la matrice diagonale constituée de la diagonale de 𝐴 (𝐷 = 𝐴 ).

i. Montrer que si la suite (𝑥 ) ∈ℕ converge vers 𝑥̅ ∈ ℝ alors 𝑥̅ est solution du système 𝐴𝑥 = 𝑏 .

ii. Exprimer les coefficients de la matrice (𝐼 − 𝛼𝐷 𝐴) en fonction des coefficients de 𝐴 .
iii. On suppose que 𝐴 est à diagonale strictement dominante et que 0 < 𝛼 ≤ 1 . Montrer que :

||𝐼 − 𝛼𝐷 𝐴|| < 1

 3

iv. Montrer que, sous les hypothèses de la question précédente, la méthode itérative (1) converge.
v. Quelle méthode retrouve-t-on quand 𝛼 = 1 .

II. Interpolation :

Question de cours.

Soit 𝑛 + 1 points d’interpolation 𝑥 , 𝑥 , … , 𝑥 d’une fonction 𝑓(𝑥). Quelle est la condition nécessaire
et suffisante pour que le polynôme d’interpolation 𝑃 (𝑥) existe et soit unique ?

Exercice 1.

Soit 𝑓(𝑥) = . Déterminer le polynôme d’interpolation de Lagrange pour les points

d’interpolation d’abscisses : -2, -1, 0, 1, 2 .

( )
Discuter l’erreur d’interpolation, pour cela on donne : 𝑓 ( ) (𝑥) = ( )

Exercice 2.

Soit 𝑓(𝑥) donnée par ses valeurs aux points d’interpolation suivants :

𝑥 0.40 0.50 0.70 0.80

𝑓(𝑥) -0.916291 -0.6993147 -0.356675 -0.223144

Donner une estimation de 𝑓(0.6).

En fait on a 𝑓(𝑥) = ln (𝑥), comparer 𝑓(0.60) et ln (0.60).

( )
( )
Donner une estimation de l’erreur théorique donnée par : 𝐸(𝑥) = 𝛾 (𝑥) ( )!
(voir le cours)

Exercice 3.

Utiliser la formule d’interpolation de Lagrange pour trouver le polynôme d’interpolation passant par
0.4, 0.5, 0.7, 0.8 pour la fonction 𝑓(𝑥) = sin𝑥. Donner une estimation de l’erreur théorique.

III. Intégration numérique :

Question de cours. Soit 𝐼 = ∫ 𝑓(𝑥)𝑑𝑥 et soient (𝑥 , 𝑥 , … , 𝑥 ) une suite de points équidistants de

pas ℎ = .

1- Donner la formule de quadrature du point milieu simple approchant 𝐼.

Donner la formule de quadrature du point milieu composée approchant 𝐼.
2- Donner la formule de quadrature des trapèzes simple approchant 𝐼.
 4

Donner la formule de quadrature des trapèzes composée approchant 𝐼.

3- Donner la formule de quadrature de Simpson simple approchant 𝐼.
Donner la formule de quadrature de Simpson composée approchant 𝐼.
4- Quel est l’ordre des formules composées ci-dessus ?

Exercice 1.

/
Déterminer par la méthode des trapèzes puis par celle de Simpson 𝐼 = ∫ 𝑓(𝑥)𝑑𝑥 sur
la base du tableau suivant :

𝑥 0 𝜋/8 2𝜋/8 3𝜋/8 𝜋/2

𝑓(𝑥) 0 0.382683 0.707107 0.923880 1

Ces points sont ceux de la fonction sin 𝑥 , comparer le résultat obtenu avec le résultat exact.

Exercice 2.

On lance une fusée verticalement du sol et l’on mesure pendant les premières 80 secondes
l’accélération 𝛾(𝑡) :

𝑡(𝑠𝑒𝑐𝑜𝑛𝑑𝑒𝑠) 0 10 20 30 40 50 60 70 80
𝛾(𝑚/𝑠) 30 31.63 33.44 35.47 37.75 40.33 43.29 46.70 50.67

Sachant que ∫ 𝛾(𝑡)𝑑𝑡 = 𝑉(𝑡) − 𝑉(0), calculer 𝑉(80 𝑠)par la méthode des trapèzes puis par la
méthode de Simpson.

Exercice 3.

Calculer par la méthode des trapèzes l’intégrale ∫ sin 𝑥 𝑑𝑥 avec un pas ℎ = puis ℎ = .

Exercice 4.

Trouver le nombre 𝑁 de subdivisions nécessaires de l’intervalle d’intégration – 𝜋, 𝜋 pour évaluer à

0.5x10-3 près, l’intégrale 𝐼 = ∫ cos 𝑥 𝑑𝑥 par la méthode de Simpson.

Pour 𝐼 = ∫ 𝑓(𝑥)𝑑𝑥 , on donne l’erreur pour la méthode de Simpson :

(𝑏 − 𝑎)
𝐸(ℎ) = ℎ max |𝑓 ( ) (𝜉)|
180 ∈[ , ]
 5

Exercice 5.

Soit 𝐹(𝑥) = ∫ 𝑡𝑒 𝑑𝑡. Combien faut-il de subdivisions de [0,1] pour évaluer 𝐹(1) à 10-8 près en
utilisant la méthode des trapèzes.

Pour 𝐼 = ∫ 𝑓(𝑥)𝑑𝑥 , on donne l’erreur pour la méthode des trapèzes :

(𝑏 − 𝑎)
𝐸(ℎ) = ℎ max |𝑓 (𝜉)|
12 ∈[ , ]

IV. Equations différentielles ordinaires :

Exercice 1.

Soit l’équation différentielle à condition initiale 𝑦 (𝑡) = 𝑦(𝑡) + 𝑡 et 𝑦(0) = 1.

Approcher la solution de cette équation en 𝑡 = 1 à l’aide de la méthode d’Euler en subdivisant
l’intervalle de travail en 10 parties égales.
Comparer à la solution exacte, donnée par 𝑦(𝑡) = −1 − 𝑡 + 2𝑒 .
Donner la solution sous forme d’un tableau :
𝒏 0 1 2 3 4 5 6 7 8 9 10
𝒕𝒏 0 - - - - - - - - - -
𝒚𝒏 1 - - - - - - - - - -

Exercice 2.

On considère l’équation différentielle suivante :

2𝑡
𝑦 (𝑡) = 𝑦(𝑡) − et 𝑦(0) = 1
𝑦

Déterminer une solution approchée de cette équation différentielle à 𝑡 = 0.2 en utilisant le schéma
de Heun avec un pas ℎ = 0.2. Comparer avec la solution théorique « 𝑦(𝑡) = √2𝑡 + 1) ».

Exercice 3.

On considère l’équation différentielle suivante :

𝑦 (𝑡) = 2𝑡 − 𝑦(𝑡), 𝑡 ∈ [0,1] 
𝑦(0) = 1

Déterminer une solution approchée de cette équation différentielle à 𝑡 = 0.3 en utilisant le schéma
d’Euler explicite avec un pas ℎ = 0.1. Comparer avec la solution théorique « 𝑦(𝑡) = 2𝑡 − 2 + 3𝑒 ) ».
 6

Exercice 4 (avec solution).

Déterminer une solution approchée par le schéma d’Euler explicite de l’équation différentielle
modélisant le mouvement d’un pendule simple donnée par :

𝑎𝑦 (𝑡) + 𝑔 sin 𝑦(𝑡) = 0 avec 𝑦(0) = 𝑦 et 𝑦′(0) = 𝑧

Solution :

C’est une équation de 2ème degré, on va la transformer en un système de 2 équations différentielles

de 1er ordre. On pose 𝑦 (𝑡) = 𝑧(𝑡) d’où 𝑦 (𝑡) = 𝑧 (𝑡). On alors le système :

𝑦 =𝑧
𝑔  
𝑧 = − sin 𝑦
𝑎

Appliquons le schéma d’Euler :

𝑦 = 𝑦 + ℎ𝑧
𝑔  
𝑧 = 𝑧 − ℎ sin 𝑦
𝑎

Ainsi, connaissant 𝑦 et 𝑦′ = 𝑧 , on peut déterminer 𝑦 et 𝑧 pour tout 𝑛.