Analyse Numérique

Calculs Numériques Approchés

Pr. LAKHBAB

EHTP
Printemps 2025

Pr. LAKHBAB Calculs Numériques Approchés 1 / 23

 Introduction

L'analyse numérique est une discipline des mathématiques

appliquées

L'objectif de l'analyse numérique est de concevoir et d'étudier des méthodes de

résolution de certains problèmes mathématiques, en général issus de la physiques,
l'ingénierie, la biologie et l'économie, et dont on cherche à calculer la solution à
l'aide d'un ordinateur.

L'accroissement de la disponibilité et de la puissance des ordinateurs depuis la

seconde moitié du 20ème siècle a grandement facilité l'utilisation de l'analyse
numérique.

Pr. LAKHBAB Calculs Numériques Approchés 3 / 23

 Introduction

Enjeux de l'analyse numérique :

Calcul approché des solutions, que l'on ne sait pas chercher autrement, et
"mieux" qu'on ne le faisait avant.
Exemple :
1 Evaluer la solution d'un système linéaires de grande taille.

2 Approcher le graphe d'une fonction connue à partir de ses valeurs en des

points.
3 Approcher une fonction inconnue qui est solution d'une équation diérentielle.

Pr. LAKHBAB Calculs Numériques Approchés 4 / 23

 Introduction

La physique numérique

La physique numérique (ou parfois physique informatique) est l'étude et

l'implémentation d'algorithmes numériques dans le but de résoudre des
problèmes physiques pour lesquels une théorie existe déjà.
Les sciences physiques se xent l'objectif de prédire les phénomènes à partir
de la connaissance d'un nombre ni de grandeurs microscopiques intervenant
dans la modélisation.
La modélisation consiste à construire un modèle mathématique permettant la
représentation du phénomène physique.

Pr. LAKHBAB Calculs Numériques Approchés 5 / 23

 Introduction

Algorithme
La recherche de résultats numériques approchés conduit à mettre au points des
algorithmes.
Dénition
Un algorithme est un ensemble de règles opératoires ou de procédés, dénis en
vue d'obtenir un résultat déterminé au moyen d'un nombre ni d'opérations.
La mise en ÷uvre pratique d'un algorithme consiste à codé celui-ci dans un
langage compréhensible par ordinateur, à l'aide d'un compilateur (traducteur), ce
qui donne un programme.

Pr. LAKHBAB Calculs Numériques Approchés 6 / 23

 Introduction

Logiciels

1 MatLab : un logiciel commercial permettant de résoudre numériquement des

problèmes mathématiques et physiques variés. La dernière version est R2016b.
[Link]
2 Scilab : est un logiciel libre développé par l'INRIA, il est très similaire au plus
connu Matlab.
[Link]
3 Maple et Mathematica : deux logiciels commerciaux qui dominent le domaine
de la résolution analytique. Les possibilités de chaque logiciel sont dans une
large mesure trés comparables.
[Link]
[Link]

Pr. LAKHBAB Calculs Numériques Approchés 7 / 23

 Calculs Numériques Approchés

Sources d'erreur

Un ordinateur ne peut représenter qu'un nombre ni de rationnels avec un nombre

ni de chires.
Donc il est nécessaire de représenter les nombres réels sous forme approchée.
Il en découle que les calculs fournis par l'ordinateur sont entachés d'erreurs.
Les erreurs d'arrondi
qui proviennent du fait qu'il n'est pas possible de représenter (tous) les réels
exactement dans un ordinateur.
Exemple :

1
= 0.142857 · · ·
x=
7
Lors de la manipulation de ce nombre, la machine devra choisir entre la
troncature ou l'arrondi à la décimale la plus proche.

Pr. LAKHBAB Calculs Numériques Approchés 9 / 23

 Calculs Numériques Approchés

Sources d'erreur

Les érreurs de propagation et génération

sont les erreurs qui apparaissent dans le résultat d'une opération comme
conséquence des erreurs des opérandes (dûs aux opérations arithmétiques
eectuées).

Exemple : Soient a = 0.3761 et b = 0.4687

Avec une représentation à trois chires signicatifs
ã = 0.376 et
b̃ = 0.468 ou b̃ = 0.469

a + b ≃ 0.844 ou 0.845

Pr. LAKHBAB Calculs Numériques Approchés 10 / 23

 Calculs Numériques Approchés

Sources d'erreur
Les érreurs de discrétisation ou de troncature
qui proviennent de simplications du modèle mathématique comme par exemple le
remplacement d'une dérivée par une diérence nie

Exemple : Evaluer de manière numérique l'exponentielle z = e x .

Une méthode consiste à utiliser les séries de Taylor
+∞ k
x
X x
z =e =
k!
k=0

En considérant des ressources de calcul limitées,

N
X xk
z̃ =
k!
k=0

On commet donc une erreur, pour le calcul de l'exponentielle !

Pr. LAKHBAB Calculs Numériques Approchés 11 / 23
 Calculs Numériques Approchés

Evaluation de l'erreur : Erreur absolue

Dénitions
Soit x , un nombre, et x̃ , une approximation de ce nombre. L'erreur absolue est
denie par :
∆x = |x − x̃|

Exemple
1 Si x = 3 et x̃ = 2, alors ∆x = 1

2 Si x = 109 et x̃ = 109 + 1, alors ∆x = 1

Dénir l'erreur relative

Pr. LAKHBAB Calculs Numériques Approchés 12 / 23

 Calculs Numériques Approchés

Evaluation de l'erreur : Erreur relative

Dénitions
Soit x , un nombre non nul, et x̃ , une approximation de ce nombre. L'erreur
relative est denie par :

x − x̃ ∆x ∆x
Er = | |= ≈
x |x| |x̃|

Exemple
1 Si x = 3 et x̃ = 2, alors Er = 50%

2 Si x = 109 et x̃ = 109 + 1, alors Er = 10−7 %

L'erreur absolue donne une mesure quantitative de l'erreur commise et l'erreur

relative en mesure l'importance.

Pr. LAKHBAB Calculs Numériques Approchés 13 / 23

 Calculs Numériques Approchés

N chires signicatifs
Dénition
⋆ Si ∆x ≤ 0.5 × 10m , alors le chire correspondant à la me puissance de 10 est
dit signicatif et tous ceux à sa gauche, correspondant aux puissances de 10
supérieures à m, le sont aussi.
x ∈ [x̃ − 0.5 × 10m , x̃ + 0.5 × 10m ]

Exemple
On donne une approximation de π au moyen de la quantité
22
x̃ = = 3.142857 · · · , on a
7
22
=⇒ ∆x = |π − | = 0.00126 · · · ≈ 0.126 × 10−2 ≤ 0.5 × 10−2
7
le chire des centièmes est signicatif et on a en tout 3 chires signicatifs
(3.14).
Considérer l'approximation x̃ = 3.1416.
Pr. LAKHBAB Calculs Numériques Approchés 14 / 23
 Calculs Numériques Approchés

Représentation approchée des nombres réels

Un ordinateur stocke généralement un nombre réel de la manière suivante :
x = (−1)s mb p (1)

0 si x est positif

s=
1 sinon .
p ∈ Z est l'exposant,
m est la mantisse, (La mantisse possède un nombre maximum N de chires
signicatif imposé par la mémoire d'un ordinateur)
N
m = 0, a1 a2 · · · aN = 0 ≤ ai ≤ b − 1, a1 ̸= 0
X
ai b −i ,
i=1

Dans la base décimale (b = 10), les ai sont les chires décimaux de m (ou les
digits)

Pr. LAKHBAB Calculs Numériques Approchés 15 / 23

 Calculs Numériques Approchés

Remarque
La condition a1 = 0 assure qu'un nombre ne peut pas avoir plusieurs
représentations. Par exemple, sans cette restriction, le nombre 101 pourrait être
représenté (dans le système décimal) par 0.1 × 100 , mais aussi par 0.01 × 101 , etc.

Exemple
1
x= 15= 0.066666666 · · · . Dans le cas d'une représentation tronquée nous
aurons, pour N = 5, x ≃ 0.66666 × 10−1 . Donc nous avons modié l'exposant
an de respecter la règle qui veut que le premier chire de la mantisse ne soit pas
nul.

Pr. LAKHBAB Calculs Numériques Approchés 16 / 23

 Arithmétiques ottante

Opérations machine

Les ressources d'une machine étant limitées ! Seul un sous ensemble F de

cardinal ni de R peut être représenté.
F = {les nombres à virgule flottante}.
Un nombre réel x est en général tronqué par la machine, dénissant ainsi un
nouveau nombre (le nombre à virgule ottante), noté fl(x), qui ne coïncide
pas nécessairement avec le nombre x original.
Pour eectuer une opération sur deux nombres réels, on eectue l'opération sur
leurs représentations ottantes et on prend ensuite la représentation ottante du
résultat.
˜
fl(x)⊕fl(y ) = (fl(x) ⊕ fl(y ))(1 + r )
⊕ ∈ {+, −, ×, /} et ⊕ ˜ représente l'opération de l'ordinateur (op ). où |r | < eps ,
et eps est la précision machine
la précision machine est dénie comme le plus petit nombre positif eps tel que
fl(1 + eps) > 1.

Pr. LAKHBAB Calculs Numériques Approchés 18 / 23

 Arithmétiques ottante

Opérations machine

Eectuer les calculs en arithmétique ottante à 3 chires signicatifs :

Perte d'associativité
Soient x = 7, 34; y = 0, 00365; z = 0, 00498. (Les calculs sont eectués avec
3 chires signicatifs).
Calcule la somme x + y + z de deux façons
1 (x + y ) + z
2 x + (y + z)
Perte de la distributivité de la multiplication par rapport à l'addition
1 (
122 333 + 695).
2 (122 × 333) + (122 × 695).

Pr. LAKHBAB Calculs Numériques Approchés 19 / 23

 Arithmétiques ottante

Opérations machine

Remarque
On ne peut additionner ou soustraire que des nombres ayant le même exposant
égal au plus grand en violant temporairement la règle que le 1er chire de la
mantisse ne soit pas nul.

Exemple (En arithmétique ottante à 3 chires)

Calculer 0.1 + 0.9 × 10−2

Pr. LAKHBAB Calculs Numériques Approchés 20 / 23

 Conditionnement et stabilité

Conditionnement et stabilité

⋄ On dit qu'un problème est bien conditionné si une petite variation de ces
données entraîne une petite variation sur les résultats. Ainsi, si un problème est
mal conditionné, alors la moindre erreur dans les données provoquera une erreur
très importante dans la solution trouvée.

⋄ Un problème peut être bien conditionné et la méthode utilisée pour le résoudre

peut être sujette à une propagation importante des erreurs numériques, on dira
que le calcul est numériquement instable. Si les erreurs introduites dans les
étapes intermédiaires ont un eet négligeable sur le résultat nal, on dira que le
calcul est numériquement stable.

Pr. LAKHBAB Calculs Numériques Approchés 22 / 23

 Conditionnement et stabilité

Conditionnement et stabilité

Soit le système linéaire :

x − 1.99995y =3

(2)
2x − 4y =1

On considère le système perturbé suivant

x − 1.9999y =3

(3)
2x − 4y =1

Pr. LAKHBAB Calculs Numériques Approchés 23 / 23