UH2C–FSJES Aı̈n Sebaâ 2019-2020
Département SMAEG Méthodes numériques S4–LF MASS
TP no 1
Exercice 1.
Ecrire en utilisant les instructions de boucles, un programme sous Matlab effectuant le pro-
duit de deux matrices construites grâces au génerateur de nombres aléatoires rand. Comparer
le programme obtenu et le produit matriciel de Matlab, en terme de temps de calcul.
Exercice 2.
On considère la fonction f définie par
0 si x ≤ a,
x−a
si a ≤ x ≤ b,
f (x) = b−a
c−x
si b ≤ x ≤ c,
c−b
0 si x ≥ c.
2.1 - On suppose que X est un vecteur de réels. Construire une fonction Matlab admettant
les réels a, b, c et le vecteur X en argument et renvoyant un vecteur contenant la valeur de
la fonction f pour chaque composante de X.
2.2 - Construire un vecteur X contenant les valeurs (−2, −1, 1, 2, 5). Taper l’instruction
f ind(X > 0).
En utilisant si besoin l’aide en ligne, expliquer ce que fait la fonction f ind.
2.3 - utiliser la fonction f ind pour optimiser le programme écrit en 2.1, en terme de temps
de calcul.
Exercice 3.
1 −1 1
3.1 - Créer sous Matlab la matrice A = 2 5 −7 et le vecteur (matrice ligne)
3 4 8
B = (2, 4, 7). Pouvez-vous faire le produit A ∗ B ?
3.2 - Construire At et B t . De même, construire une matrice diagonale admettant les valeurs
de B comme éléments de la diagonale ainsi qu’un vecteur contenant la diagonale de A.
3.3 - Construire une matrice de taille (3, 5) ne contenant que des 0, une matrice de taille
(5, 4) ne contenant que des 1, une matrice identité de taille (5, 5).
3.4 - Extraire les parties triangulaires supérieure et inférieure de la matrice A.
3.5 - Construire une matrice C aléatoire de taille (8, 6). En extraire des sous-matrices.
1
UH2C–FSJES Aı̈n Sebaâ 2019-2020
Département SMAEG Méthodes numériques S4–LF MASS
Exercice 4. Approximation de π et calcul d’intégrale par la
méthode de Monte Carlo
Principe de la méthode : On place N points au hasard dans un domaine d’aire connu Ad .
Soit n le nombre de points ”tombés” dans un sous domaine d’aire As . On a alors
n
As = lim Ad .
N →+∞ N
O 1
Objectif : Donner une approximation de de l’aire du quart de disque de centre (0, 0) et de
rayon 1 contenu dans le carré [0, 1]2 , c’est-à-dire de π/4.
4.1 - Ecrire une fonction matlab MonteC qui prend comme argument N le nombre de
points à tirer aléatoirement (fonction rand de matlab) dans le carré [0, 1]2 et renvoie une
approximation de π.
4.2 - Tester la méthode pour N = 100, N = 1000, N = 10000 et N = 100000.
4.3 - Tracer l’erreur, en fonction de N entre l’approximation calculée et la valeur exacte
de π.
4.4 - Répondre aux mêmes questions pour donner une approximation de l’intégrale
Z 1
I= x2 dx.
0
Exercice 5.
Coder et tester les programmes des exercices de TD1.
