P ROGRAMMATION AVEC MATLAB

Par

Dr Kounhinir SOME
Maître de conférence, CAMES

UFR-ST, UNZ
Burkina Faso

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 1 / 103

Plan
1 Introduction
2 Généralités et prise en main
3 Calculs élémentaires
4 Les vecteurs
5 Graphiques
6 Résolution des systèmes
7 Fonctions numériques
8 Compléments
9 Programmation avec MATLAB
10 Optimisation
11 Résolution numérique des EDO et des EDP
12 Exercices
Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 2 / 103
 Introduction

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 3 / 103

Introduction
Types d’outils, Société, Domaine

 Type d’outil
MATLAB est un langage de programmation et un
environnement de développement.
 Société
MATLAB est commercialisé par The MathWords une
société américaine.
 Domaine
MATLAB est utilisé aujourd’hui dans tous les domaines :
Calcul numérique, éducation, recherche, industrie,
développement de projet, ...

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 4 / 103

Introduction
Historique, Capacité

 Historique MATLAB veut dire Matrix Laboratory. Il fut conçu

par le mathématicien Cleve Moler de l’Université du
Nouveau Mexique à la fin des années 70 à partir des
bibliothèques FORTRAN, LINPACK, et EISPACK. Par suite,
Jack Little, un ingénieur, comprend rapidement les capicités
et entreprend avec Steve Bangert, de le recoder en C. En
1984 Jack Little, Cleve Moller et Steve Bangert créèrent la
société The MathWorks afin de commercialiser la version
1.0 de MATLAB.
 Capacité
MATLAB est un langage interactif, de haut niveau exécution
de tâches plus rapide que C, C++ et FORTRAN.

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 5 / 103

Introduction
Outils et Modules associés, Alternatives

 Outils et Modules associés

Plusieurs outils sont incorporés dans MATLAB pour la
résolution de problèmes spécifiques : Communication
Toolbox, Control System Toolbox, Excel Link, MATLAB
Compiler, Neural Network Toolbox, Optimisation Toolbox,
Real-Time Workshopr, Robust Control Toolbox,
SimMechanics, Simulink, Statistics Toolbox, System
Identification Toolbox, Virtual Reality Toolbox.
 Alternatives
D’autres outils ont été développés, certains gratuits, et sont
des alternatives à MATLAB tels que : FERMAT, GNU,
OCTAVE, IDL, JMathLib, LMS [Link], AMEsim,
Scilab, Sysquake.

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 6 / 103

 Généralités et prise en main

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 7 / 103

Généralités et prise en main
Aide, documents en ligne, espace de travail

Pour démarrer le MATLAB, double cliquer sur l’icône de

démarrage

F IGURE: MATLAB 7.1

On obtient
Master de Physiqueune fenêtre subdivisées
(UFR-ST) en 4 partir :
Dr Kounhinir SOME Avril 2022 8 / 103
Généralités et prise en main
Aide, documents en ligne, espace de travail

 Tout en haut la barre des menus qui se présente dans la

même configuration que WORD ;
 En dessous de la petite fenêtre qui affiche le répertoire
courant de travail et les document qui s’y trouvent ;
Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 9 / 103
Généralités et prise en main
Aide, documents en ligne, espace de travail

 En dessous de cette fenêtre une autre fenêtre qui donne

l’historique des commandes utilisées ;
 A droite la fenêtre d’exécution des commandes MALAB ou
workspace.
Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 10 / 103
Généralités et prise en main
Aide, documents en ligne, espace de travail

Remarque
Il existe plusieurs fonctions intégrées dans MATLAB qu’on peut
utiliser sans forcement toutes les retenir. Il est donc intéressant
de savoir utiliser l’aide en ligne qui donne toute les fonctions et
les syntaxes correspondantes. L’aide est uniquement disponible
en langue anglaise. Elle permet de trouver toutes les fonctions
utiles.

Aller dans le menu le plus en haut cliquer sur Help puis sur
MATLAB help ou Product help une nouvelle fenêtre s’ouvre
partagée en deux :

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 11 / 103

Généralités et prise en main
Aide, documents en ligne, espace de travail

A gauche pour la recherche et à droite les résultats de la

recherche
Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 12 / 103
Généralités et prise en main
Aide, documents en ligne, espace de travail

Nous avons les espaces suivants dans help :

 Contents :Il est dans la fenêtre de gauche. Permet de
présenter la table des matière en ligne. En cliquant sur un
élément présente dans cette liste, son contenu est affiché
dans la fenêtre de droite.
 Index : L’index de l’aide en ligne, qui répertorie toutes le
commandes MATLAB et permet d’accéder à leur mode
d’emploi spécifique. Pour une commande connue, on peut
taper le nom dans la fenêtre de saisie Search index. Si on
ne connaît pas la commande on peut la rechercher dans la
liste alphabétique proposée. Pour avancer dans l’alphabet
taper la première lettre de la commande recherchée dans la
fenêtre de saisie.

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 13 / 103

Généralités et prise en main
Aide, documents en ligne, espace de travail

Remarque
Si on connait la commande et on ne maîtrises pas le mode
d’emploi on peut aussi taper dans le Workspace >lookfor nom
de la commande.

 Favorite : On peut stocker dans favoris un chapitre que l’on

pense l’utiliser couramment. Pour cela on le sélection et on
clique sur favorite dans la barre des menus puis on clique
sur add to favorites en enfin sur save
 Demos : C’est une aide très importante. Elle se présente
sous la forme de programme de démonstration auquels on
peut accéder ? partir de help puis demos

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 14 / 103

 Calculs élémentaires

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 15 / 103

Calculs élémentaires
Déclaration des variables

Remarque
Chaque commande passée à MATLAB est interprétée, le
résultat est affiché, et l’utilisateur reçoit de nouveau le prompt 
pour la commande suivante : c’est le mode interactif

MATLAB permet de créer et d’initialiser des variables. La

déclaration des variables en matlab suit les règles suivantes :
 toutes les variables sont des matrices ;
 pas de déclaration de type. Le type est établi
automatiquement à partir des valeurs affectées à la
variable.

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 16 / 103

Calculs élémentaires
Déclaration des variables

Les types de variables utilisés par MATLAB sont les suivants :

 réel ;
 complexe ;
 chaîne de caractères
 logique

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 17 / 103

Calculs élémentaires
Exemples

Example
  a = 2, ’définit la variable a de type réel et lui donne la
valeur 2 ; (taper entrer)
  b = 3, ’définit la variable b de type réel et lui donne la
valeur 3 ; (taper entrer)

Remarque
Ces valeurs sont affectées et garder en mémoire à chaque fois
que tu tape a ou b il va t’afficher la même jusqu’à ce qu’on
rentre une nouvelle valeur pour la variable. par exemple en
faisant  b = a

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 18 / 103

Calculs élémentaires
Exemples

  z = 2 + i et  z1 = 3 + j, ’Définit la variable z de type

complexe et lui donne la valeur 2 + i’
  rep =0 oui 0 , ’Définit la variable rep de variable de type
chaîne de caractères et lui attribue la valeur oui.’
  whos Permet d’afficher toutes les variables avec
indication sur leur taille, leur codage et leur type.

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 19 / 103

Calculs élémentaires
Exemples

En écrivant Whos et en tapant sur entrer on obtient le tableau

suivant :
Name Size Bytes Class
a 1×1 8 double array
b 1×1 6 double array
rep 1×3 6 char array
z 1×1 16 double array
z1 1×1 16 double array

On peut aussi retrouver le type d’une variable grâce aux

commandes : ischar, islogical et isreal :

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 20 / 103

Calculs élémentaires
Exemples

 ischar(x) retourne 1 si x est de type chaîne de caractères et

0 sinon ;
 islogical(x) retourne 1 si x est de type logique et 0 sinon.
 isreal(x) est à utiliser avec discernement : elle retourne 1 si
x est réel ou de type chaîne de caractères et 0 sinon (x est
complexe à partie imaginaire non nulle ou n’est pas un
tableau de valeurs réelles ou de caractères).
Par exemple :
  ischar(rep), donne ans=1
  ischar(x), donne ans=0
  isreal(z), donne ans=0
  isreal(x), donne ans=1
  isreal(rep), donne ans=1

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 21 / 103

Calculs élémentaires
Opérations élémentaires

  1 + 9, donne ans=10.
  6/2, donne ans=3 (division à droite).
  2 \ 6, donne ans=3 (division à gauche).
  1/3, donne ans=0.3333.
  format long, donne ans=0.333333333333333 (permet
d’afficher une variable avec plus de chiffre après la virgule ;
15 chiffres).
  format short est le format par défaut.
  format rat est le format rationnel, il permet d’afficher les
nombres à virgule sous forme de fraction.
  25 (2 ˆ 5) donne ans=32.

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 22 / 103

Calculs élémentaires
Variables spéciales

  pi, donne ans=3.1416

  realmin, donne ans = 2.2251e-308 (Plus petit entier
positif que le logiciel peut enregistrer.)
  realmax, ans =1.7977e+308 (Plus grand entier positif
que le logiciel peut enregistrer.)
  eps, donne ans =2.2204e-016 (Valeur assez petite et
positive)

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 23 / 103

Calculs élémentaires
Fonctions mathématiques

  sin(pi/2), donne ans=1

 il existe bien d’autres fonctions telles :
 cos et tan
 sinh, cosh et tanh
 exp et log
 sqrt
 ...

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 24 / 103

Calculs élémentaires
Gestion des variables

  who pour afficher toutes les variables.

  whos pour afficher toutes les variables avec indication
sur leur taille.
  size(a) pour afficher les dimensions de la matrice a.
  clear a pour effacer la variable a.
  clear all pour effacer toutes les variables.

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 25 / 103

 Vecteurs

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 26 / 103

Les vecteurs
Création par énumération

  v 1 = [1254], donne v = 1254

  v = [1 2 5 4], donne v = 1 2 5 4 ou bien
  v = [1, 2, 5, 4], donne la même chose
 
1
 2 
  w = [1; 2; 5; 4], donne w =  5  ou bien on peut


4
écrire
  w = v 0 cela donne le même w

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 27 / 103

Les vecteurs
Création par description

  x = 0 : pi/11 : pi, donne :

X= Columns 1 through 7
0 0.2856 0.5712 0.8568 1.1424 1.4280 1.7136
Columns 8 through 12
1.9992 2.2848 2.5704 2.8560 3.1416
  y = linspace (0, pi, 12)
y = Columns 1 through 7
0 0.2856 0.5712 0.8568 1.1424 1.4280 1.7136
Columns 8 through 12
1.9992 2.2848 2.5704 2.8560 3.1416
  y = linspace (0, pi ), donne 100 éléments par défaut.

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 28 / 103

Les vecteurs
Indexation et adressage

Pour afficher des éléments d’un vecteur :

  v (3), donne ans=5
  v (2 : 1 : 4), donne ans=2 5 4
  v (2 : 4), donne la même reponse
  v ([4, 3, 2, 3]), donne ans= 4 5 2 5
  V ([2, 3]) = [0, 10], donne v=1 0 10 4

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 29 / 103

Les vecteurs
Opérations scalaires

  a = 2 ∗ v , donne a= 2 0 20 8
  b = a − 5, donne b= -3 -5 15 3
  c = a. ˆ 2, donne c= 4 0 400 64
  d = a ∗ a0 , donne d=468

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 30 / 103

Les vecteurs
Combinaison de vecteurs

  K = [a, b ], donne K=2 0 20 8 -3 -5 15 3

  L = [K ([3, 7, 4]), K ([5 : −1 : 1])], donne L=20 15 8 -3
8 20 0 2

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 31 / 103

Les vecteurs
Récapitulatif : utilisation des symboles

 [] énumération d’éléments
 : description d’éléments
 () ensemble d’arguments
 , séparateur d’arguments
 ; séparateur des lignes ; non affichage du résultat de
l’exécution d’une instruction
 . force l’opérateur à s’appliquer sur chaque éléments.

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 32 / 103

 Matrices

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 33 / 103

Les matrices
Création

Une matrice est un ensemble de lignes comportant toutes le

même nombre de colonnes.
 M1 = [12 3; 4 5 6; 
7 8 9]  M2 = [1 2 3; 11 13 15]
1 2 3  
1 2 3
M1 =  4 5 6  M2 =
11 13 15
7 8 9

Remarque
On pouvait saisir les matrices M1 et M2 par :
M1 = [1 : 3; 4 : 6; 7 : 9] et M2 = [1 : 1 : 3; 11 : 2 : 15]

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 34 / 103

Les matrices
Opérations scalaires

 M3
 = 10 ∗ M1    M4 = M1/2 
10 20 30 0.5000 1.000 1.5000
M3 =  40 50 60  M4 =  2.0000 2.5000 3.0000 
70 80 90 3.5000 4.000 4.5000

 M5 = M1 ˆ 2  M6 = M1. ˆ 2 
30 36 42 1 4 9
M5 =  66 81 96  M6 =  16 25 36 
102 126 150 49 64 81

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 35 / 103

Les matrices
Opérations entre matrices

N1 = M2 ∗ M1   N2 = M1 ∗ M2
30 36 42
N1 = N2 impossible
168 207 246

 P = [4 8 6; 7 9 4; 7 2 1]  N4 = M2 ∗ inv (P )
N3 = M2/P
 N3 =   N4 = 
0.8525 −0.5902 0.2459 0.8525 −0.5902 0.2459
3.6967 −2.2131 1.6721 3.6967 −2.2131 1.6721

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 36 / 103

Les matrices
Matrices particulières


 eye (4)   zeros (2, 3)
1 0 0 0  
 0 1 0 0  0 0 0
ans = 
 0 0 1
 ans =
0  0 0 0
0 0 0 1

 ones (2, 4)  diag

 ([4, 5, 3])
  4 0 0
1 1 1 1
ans = ans =  0 5 0 
1 1 1 1
0 0 3

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 37 / 103

Les matrices
Exercices d’application

On considère la matrice M1 précédente.

1 Renommer cette matrice A.
2 Mettre à 0 l’élément de la 2e ligne 3e colonne.
3 Mettre à 4 tous les éléments de la 1? ligne.
4 Mettre à 5 tous les éléments de la 2e colonne.
5 Appeler cette matrice B.
6 Créer la matrice C en accolant à droite de M1 la matrice B
7 Créer la matrice D en extrayant les 2 premières lignes et les
2 dernières colonnes de A.
8 Créer le matrice diagonale E dont la diagonale principale
est égale à la diagonale principale de A.
Remarque
Faire toutes ces opérations sans boucle.
Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 38 / 103
Les matrices
Compléments

Déterminant : Le déterminant des matrices s’obtient en

utilisant la fonction det. La syntaxe est la suivante pour
toute matrice A :
 det (A)

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 39 / 103

Les matrices
Compléments

Déterminant : Le déterminant des matrices s’obtient en

utilisant la fonction det. La syntaxe est la suivante pour
toute matrice A :
 det (A)

Example
Calculer le déterminant des matrices A, B et C suivantes

A = [2, 0; 12, 6] et B = [1, 5, −2; 6, 4, 11; −7, 3, 9]

C = [2, 6, 9; −1, −5, 13]

.

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 39 / 103

Les matrices
Compléments

Diagonalisation : L’utilisation de la fonction eig permet de

déterminer pour toute matrice carrée A, la matrice de
passage (P ) et la matrice diagonale (D ). La syntaxe est la
suivante :
 [P, D ] = eig (A)

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 40 / 103

Les matrices
Compléments

Diagonalisation : L’utilisation de la fonction eig permet de

déterminer pour toute matrice carrée A, la matrice de
passage (P ) et la matrice diagonale (D ). La syntaxe est la
suivante :
 [P, D ] = eig (A)

Example
Diagonaliser les matrices A et B suivantes

A = [2, 0; 12, 6] et B = [1, 5, −2; 6, 4, 11; −7, 3, 9]

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 40 / 103

Les matrices
Compléments

Décomposition : Une matrice carrée A peut être

décomposée de sorte à lui extraire sa matrice triangulaire
supérieure et/ou sa ma matrice triangulaire inférieure en se
servant respectivement des fonctions triu et tril. La syntaxe
est la suivante
 triu (A) ou tril (A)

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 41 / 103

Les matrices
Compléments

Décomposition : Une matrice carrée A peut être

décomposée de sorte à lui extraire sa matrice triangulaire
supérieure et/ou sa ma matrice triangulaire inférieure en se
servant respectivement des fonctions triu et tril. La syntaxe
est la suivante
 triu (A) ou tril (A)

Example
Extraire des matrices A et B les matrices triangulaires

A = [2, 0; 12, 6] et B = [1, 5, −2; 6, 4, 11; −7, 3, 9]

.
Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 41 / 103
Les matrices
Compléments

Normes : Avec matlab on peut calculer la norme des

vecteurs et aussi celle des matrices pour les normes
classiques que nous connaissons déjà. Pour un vecteur ou
une matrice A est syntaxe est :

 norm (A, t )
avec t=1, 2, inf ou ’fro’ (Frobenius).

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 42 / 103

Les matrices
Compléments

Normes : Avec matlab on peut calculer la norme des

vecteurs et aussi celle des matrices pour les normes
classiques que nous connaissons déjà. Pour un vecteur ou
une matrice A est syntaxe est :

 norm (A, t )
avec t=1, 2, inf ou ’fro’ (Frobenius).
Example
Calculer pour les A et B les quatre normes

A = [2, 0, −8, 12, 6] et B = [1, 5, −2; 6, 4, 11; −7, 3, 9]

.
Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 42 / 103
 Graphiques

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 43 / 103

Graphiques
Dimension 2

Pour faire une représentation graphique on utilise la fonction

plot. Par exemple representons f (x ) = sin(x ) :
  x = linspace (0, 2 ∗ π, 30); y = sin(x ) ;

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 44 / 103

Graphiques
Dimension 2

Pour faire une représentation graphique on utilise la fonction

plot. Par exemple representons f (x ) = sin(x ) :
  x = linspace (0, 2 ∗ π, 30); y = sin(x ) ;

  plot(x,y)

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 44 / 103

Graphiques
Dimension 2

Pour faire une représentation graphique on utilise la fonction

plot. Par exemple representons f (x ) = sin(x ) :
  x = linspace (0, 2 ∗ π, 30); y = sin(x ) ;

  plot(x,y)

  plot(x,cos(x))

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 44 / 103

Graphiques
Dimension 2

Pour faire une représentation graphique on utilise la fonction

plot. Par exemple representons f (x ) = sin(x ) :
  x = linspace (0, 2 ∗ π, 30); y = sin(x ) ;

  plot(x,y)

  plot(x,cos(x))

  hold on

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 44 / 103

Graphiques
Dimension 2

Pour faire une représentation graphique on utilise la fonction

plot. Par exemple representons f (x ) = sin(x ) :
  x = linspace (0, 2 ∗ π, 30); y = sin(x ) ;

  plot(x,y)

  plot(x,cos(x))

  hold on

  plot(x,sin(x))

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 44 / 103

Graphiques
Dimension 2

  hold off

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 45 / 103

Graphiques
Dimension 2

  hold off
  clf

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 45 / 103

Graphiques
Dimension 2

  hold off
  clf
  plot (x, sin(x ), x, cos (x ))

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 45 / 103

Graphiques
Dimension 2

  hold off
  clf
  plot (x, sin(x ), x, cos (x ))
  legend(’courbe de sinus’,’courbe de cosinus’)

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 45 / 103

Graphiques
Dimension 2

  hold off
  clf
  plot (x, sin(x ), x, cos (x ))
  legend(’courbe de sinus’,’courbe de cosinus’)
  title(’courbes de fonctions circulaire’)

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 45 / 103

Graphiques
Dimension 2

  hold off
  clf
  plot (x, sin(x ), x, cos (x ))
  legend(’courbe de sinus’,’courbe de cosinus’)
  title(’courbes de fonctions circulaire’)
  xlabel(’abscisse’),ylabel(’ordonnée’)

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 45 / 103

Graphiques
Dimension 2

  hold off
  clf
  plot (x, sin(x ), x, cos (x ))
  legend(’courbe de sinus’,’courbe de cosinus’)
  title(’courbes de fonctions circulaire’)
  xlabel(’abscisse’),ylabel(’ordonnée’)
  plot (x, sin(x ),0 r +0 , x, cos(x ),0 b <0 )

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 45 / 103

Graphiques
Dimension 2

  hold off
  clf
  plot (x, sin(x ), x, cos (x ))
  legend(’courbe de sinus’,’courbe de cosinus’)
  title(’courbes de fonctions circulaire’)
  xlabel(’abscisse’),ylabel(’ordonnée’)
  plot (x, sin(x ),0 r +0 , x, cos(x ),0 b <0 )
  help plot

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 45 / 103

Gaphiques
Dimension 3

On utilse la fonction plot3 pour les courbes paramétrées et la

fonction mesh pour les courbes de fonctions de deux variables.

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 46 / 103

Gaphiques
Dimension 3

On utilse la fonction plot3 pour les courbes paramétrées et la

fonction mesh pour les courbes de fonctions de deux variables.
Example
représenter la courbe paramétrée définie par :

x = sin(t )
f (x, y ) = ∀t ∈ [0, 10 ∗ π ]
y = cos(t )

  t = linspace (0, 10 ∗ π ) ;
  plot3(sin(t ), cos(t ), t )

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 46 / 103

Graphiques
Dimension 3

Pour les fonctions à deux(2) variables on définit un maillage des

avec les deux variables x et y en utilisant la fonction meshgrid.
Par exemple considérons les points x et y définie
respectivement par :
  x = [1, 2, 3] ; y = [4, 5, 6];

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 47 / 103

Graphiques
Dimension 3

Pour les fonctions à deux(2) variables on définit un maillage des

avec les deux variables x et y en utilisant la fonction meshgrid.
Par exemple considérons les points x et y définie
respectivement par :
  x = [1, 2, 3] ; y = [4, 5, 6];
  [X , Y ] = meshgrid (x, y ) ;

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 47 / 103

Graphiques
Dimension 3

Pour les fonctions à deux(2) variables on définit un maillage des

avec les deux variables x et y en utilisant la fonction meshgrid.
Par exemple considérons les points x et y définie
respectivement par :
  x = [1, 2, 3] ; y = [4, 5, 6];
  [X , Y ] = meshgrid (x, y ) ;

1 2 3 4 4 4
X = 1 2 3 et Y = 5 5 5
1 2 3 6 6 6

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 47 / 103

Graphiques
Dimension 3

Pour les fonctions à deux(2) variables on définit un maillage des

avec les deux variables x et y en utilisant la fonction meshgrid.
Par exemple considérons les points x et y définie
respectivement par :
  x = [1, 2, 3] ; y = [4, 5, 6];
  [X , Y ] = meshgrid (x, y ) ;

1 2 3 4 4 4
X = 1 2 3 et Y = 5 5 5
1 2 3 6 6 6
1 2 3
4 (1,4) (2,4) (3,4)
5 (1,5) (2,5) (3,5)
6 (1,6) (2,6) (3,6)
Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 47 / 103
Graphiques
Dimension 3

Example
Représenter la fonction

f (x, y ) = cos(x + y ), ∀(x, y ) ∈ [0, 10] × [−10, 10]

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 48 / 103

Graphiques
Dimension 3

Example
Représenter la fonction

f (x, y ) = cos(x + y ), ∀(x, y ) ∈ [0, 10] × [−10, 10]

  x = 0 : 0.1 : 10 ; y = −10 : 0.1 : 10;

  [X , Y ] = meshgrid (x, y ) ;
  Z = cos(X + Y ) ;
  mesh(X , Y , Z )

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 48 / 103

Graphiques
Application

Exercice
Représenter la fonction numérique sur le domaine D = [−π, π ]2
en utilisant successivement les fonctions mesh et surf
sin(x 2 + y 2 )
f (x, y ) =
x2 + y2
et
2 −y 2
g (x, y ) = xe −x

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 49 / 103

Gaphiques
Application

Exercice
Le but de cet exercice est d’interpoler ces donner, c’est dire
tracer une courbe qui passe le plus possible par l’ensemble des
point d’abscisse le temps et l’ordonnée la mesure
correspondante.
Temps 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Mesures 0 0.5 1 1.3 2 1.9 2 3 1.5 1 0.5

1 Utiliser d’abord la fonction plot

2 utiliser les fonctions : polyfit, polyval, spline (Cherche ces
fonctions dans help).

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 50 / 103

 Résolution des systèmes

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 51 / 103

Résolution des systèmes
Décomposition de matrice et valeurs et vecteurs propres et déterminant

Matlab permet de séparer des matrices selon plusieurs

décompositions tels que LU, QR, SVD..
  A = [1, 2, 3; 4, 5, 6; 7, 8, 9] ;
  [L, U, P ] = lu (A) ; avec P ∗ A = L ∗ U
  [Q, R ] = qr (A)
  [U, S, V ] = svd (A)

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 52 / 103

Résolution des systèmes
Décomposition de matrice et valeurs et vecteurs propres et déterminant

Matlab permet de séparer des matrices selon plusieurs

décompositions tels que LU, QR, SVD..
  A = [1, 2, 3; 4, 5, 6; 7, 8, 9] ;
  [L, U, P ] = lu (A) ; avec P ∗ A = L ∗ U
  [Q, R ] = qr (A)
  [U, S, V ] = svd (A)

Pour les valeurs et vecteurs propres l’instruction est :

  [P, D ] = eig (A)
  det(A).

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 52 / 103

Résolution des systèmes
Décomposition de matrice et valeurs et vecteurs propres et déterminant

Exercice
Résoudre le système d’équations linéaires suivant :


 x − 2y + 2z − t = 9
−2x + 5y − 3z + 5t = −27


 2x − 3y + 4z + 2t = 5
−x + 5y + 2z + 10t = −35


Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 53 / 103

 Fonctions numériques

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 54 / 103

Fonctions numériques
Simplification d’équations

Matlab dispose des fonctions qui permet de simplifier les

équations symboliques. Pour ce faire, il suffit de déclarer les
différents symboles via la commande syms et ensuite donner ?
la fonction simple/simplify l’équation
  syms x ;
  f = cos2 (x ) + sin2 (x ) ;
  f = simple (f ) ; il affiche f = 1

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 55 / 103

Fonctions numériques
Simplification d’équations

Matlab dispose des fonctions qui permet de simplifier les

équations symboliques. Pour ce faire, il suffit de déclarer les
différents symboles via la commande syms et ensuite donner ?
la fonction simple/simplify l’équation
  syms x ;
  f = cos2 (x ) + sin2 (x ) ;
  f = simple (f ) ; il affiche f = 1
  syms x ;
  g = cos(3 ∗ acos (x )) ;
  simple (g ) il affiche 4 ∗ x 3 − 3 ∗ x

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 55 / 103

Fonctions numériques
Dérivées et intégrales

Matlab est capable de dériver ou intégrer une des expressions

symboliques.
  syms x ;
  diff (sin(x 2 )) ; il affiche 2 ∗ x ∗ cos(x 2 )
  diff (x 6 , 6) ; il affiche 720
  int (−2 ∗ x /(1 + x 2 )2 ) ; il affiche 1/(x 2 + 1)
1
  int (x ∗ log (1 + x ), 0, 1) ; il affiche
4

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 56 / 103

Fonctions numériques
Dérivées et intégrales

Matlab est capable de dériver ou intégrer une des expressions

symboliques.
  syms x ;
  diff (sin(x 2 )) ; il affiche 2 ∗ x ∗ cos(x 2 )
  diff (x 6 , 6) ; il affiche 720
  int (−2 ∗ x /(1 + x 2 )2 ) ; il affiche 1/(x 2 + 1)
1
  int (x ∗ log (1 + x ), 0, 1) ; il affiche
4
  syms x t ;
  diff (sin(x ∗ t 2 ), t ) ; il affiche 2 ∗ t ∗ x ∗ cos(x ∗ t 2 )
  int (x /(1 + t 2 ), t ) ; il affiche x ∗ atan(t )
  int (2 ∗ x, sin(t ), 1) ; il affiche cos2 (t )
Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 56 / 103
Fonctions numériques
Équations et systèmes linéaires

Matlab est capable de résoudre des équations avec des

paramètres inconnus en utilisant la fonction solve.
  syms a b c x ;
  solve (a ∗ x 2 + b ∗ x + c )
  solve (sin(t + m )) il affiche −m
  solve (a ∗ x 2 + b ∗ x + c, b ) ; il affiche −(a ∗ x 2 + c )/x

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 57 / 103

Fonctions numériques
Équations et systèmes linéaires

Matlab est capable de résoudre des équations avec des

paramètres inconnus en utilisant la fonction solve.
  syms a b c x ;
  solve (a ∗ x 2 + b ∗ x + c )
  solve (sin(t + m )) il affiche −m
  solve (a ∗ x 2 + b ∗ x + c, b ) ; il affiche −(a ∗ x 2 + c )/x
  syms x y ;
  S = solve (x + 2y == 1, x − 1 ∗ y == 5) ;
  S = [S.x, S.y ] ; il affiche [11/3, −4/3]

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 57 / 103

Fonctions numériques
Équations et systèmes linéaires

Matlab est capable de résoudre des équations avec des

paramètres inconnus en utilisant la fonction solve.
  syms a b c x ;
  solve (a ∗ x 2 + b ∗ x + c )
  solve (sin(t + m )) il affiche −m
  solve (a ∗ x 2 + b ∗ x + c, b ) ; il affiche −(a ∗ x 2 + c )/x
  syms x y ;
  S = solve (x + 2y == 1, x − 1 ∗ y == 5) ;
  S = [S.x, S.y ] ; il affiche [11/3, −4/3]
  syms x, y, a ;
  S = solve (x + a ∗ y == 1, x − y == 5) ;
  S = [S.x, S.y ] ; il affiche
[(1 + 5 ∗ a)/(1 + a), −4/(1 + a)]
Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 57 / 103
Fonctions numériques
Équations différentielles

Matlab est capable de dériver ou intégrer une des expressions

symboliques.
  dsolve (0 Dy == ay 0 ,0 y (0) == b 0 ) il affiche b ∗ exp (a ∗ t )
 
dsolve (0 D2y == −2 ∗ a ∗ y 0 ,0 y (0) = 10 ,0 Dy (π/a) == 00 )

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 58 / 103

Fonctions numériques
Équations différentielles

Matlab est capable de dériver ou intégrer une des expressions

symboliques.
  dsolve (0 Dy == ay 0 ,0 y (0) == b 0 ) il affiche b ∗ exp (a ∗ t )
 
dsolve (0 D2y == −2 ∗ a ∗ y 0 ,0 y (0) = 10 ,0 Dy (π/a) == 00 )
  syms x, y ;
  Z = dsolve (0 Dx == y 0 ,0 Dy == −x 0 ) ;
  Z = [Z .x, Z .y ] ;

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 58 / 103

Fonctions numériques
Équations différentielles

Matlab est capable de dériver ou intégrer une des expressions

symboliques.
  dsolve (0 Dy == ay 0 ,0 y (0) == b 0 ) il affiche b ∗ exp (a ∗ t )
 
dsolve (0 D2y == −2 ∗ a ∗ y 0 ,0 y (0) = 10 ,0 Dy (π/a) == 00 )
  syms x, y ;
  Z = dsolve (0 Dx == y 0 ,0 Dy == −x 0 ) ;
  Z = [Z .x, Z .y ] ;

  syms x, y ;
  Z = dsolve (0 Dx = y 0 ,0 Dy = −x 0 ,0 y (0) = 2, x (0) = 10 ) ;
  Z = [Z .x, Z .y ] ;
Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 58 / 103
 Compléments

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 59 / 103

Compléments
Vecteurs : Opérations

Matlab dispose de fonctions qui réalisent des opérations

propres aux vecteurs. Il s’agit des fonctions suivantes :
Considérons le vecteur x défini par x = [2, 4, 5]

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 60 / 103

Compléments
Vecteurs : Opérations

Matlab dispose de fonctions qui réalisent des opérations

propres aux vecteurs. Il s’agit des fonctions suivantes :
Considérons le vecteur x défini par x = [2, 4, 5]
 sum. Permet de calculer la somme des éléments du vecteur
 sum (x ) affiche 11

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 60 / 103

Compléments
Vecteurs : Opérations

Matlab dispose de fonctions qui réalisent des opérations

propres aux vecteurs. Il s’agit des fonctions suivantes :
Considérons le vecteur x défini par x = [2, 4, 5]
 sum. Permet de calculer la somme des éléments du vecteur
 sum (x ) affiche 11
 prod. Permet de calculer le produit des éléments du vecteur
 prod (x ) affiche 40

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 60 / 103

Compléments
Vecteurs : Opérations

Matlab dispose de fonctions qui réalisent des opérations

propres aux vecteurs. Il s’agit des fonctions suivantes :
Considérons le vecteur x défini par x = [2, 4, 5]
 sum. Permet de calculer la somme des éléments du vecteur
 sum (x ) affiche 11
 prod. Permet de calculer le produit des éléments du vecteur
 prod (x ) affiche 40
 diff. Permet de calculer les différences entre les éléments
consécutifs du vecteur de la droite vers la gauche.
 diff (x ) affiche 2 1

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 60 / 103

Compléments
Vecteurs : Opérations

Matlab dispose de fonctions qui réalisent des opérations

propres aux vecteurs. Il s’agit des fonctions suivantes :
Considérons le vecteur x défini par x = [2, 4, 5]
 sum. Permet de calculer la somme des éléments du vecteur
 sum (x ) affiche 11
 prod. Permet de calculer le produit des éléments du vecteur
 prod (x ) affiche 40
 diff. Permet de calculer les différences entre les éléments
consécutifs du vecteur de la droite vers la gauche.
 diff (x ) affiche 2 1
 mean et std. Permet de calculer respectivement la moyenne
et l’écart type des éléments du vecteur
 mean(x ); std (x ) affiche 3, 6667 et 1.5275
Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 60 / 103
Compléments
Vecteurs : Sous forme polynôme

Matlab représente un polynôme sous forme d’un tableau de ses

coefficients classés dans l’ordre des puissances décroissantes.
Considérons le polynôme P d’expression P (x ) = x 2 − 6x + 9.
Alors il est représenté sous Matlab par :
 P = [1 − 6 9]

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 61 / 103

Compléments
Vecteurs : Sous forme polynôme

Matlab représente un polynôme sous forme d’un tableau de ses

coefficients classés dans l’ordre des puissances décroissantes.
Considérons le polynôme P d’expression P (x ) = x 2 − 6x + 9.
Alors il est représenté sous Matlab par :
 P = [1 − 6 9]
On peut calculer les racines d’un polynôme par la fonction roots.
Par exemple
 roots (P ) affiche 3 3

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 61 / 103

Compléments
Vecteurs : Sous forme polynôme

Matlab représente un polynôme sous forme d’un tableau de ses

coefficients classés dans l’ordre des puissances décroissantes.
Considérons le polynôme P d’expression P (x ) = x 2 − 6x + 9.
Alors il est représenté sous Matlab par :
 P = [1 − 6 9]
On peut calculer les racines d’un polynôme par la fonction roots.
Par exemple
 roots (P ) affiche 3 3
Example
Déterminer les racines du polynôme Q (x ) = x 3 + 2x − 3

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 61 / 103

Compléments
Vecteurs : Sous forme polynôme

Matlab représente un polynôme sous forme d’un tableau de ses

coefficients classés dans l’ordre des puissances décroissantes.
Considérons le polynôme P d’expression P (x ) = x 2 − 6x + 9.
Alors il est représenté sous Matlab par :
 P = [1 − 6 9]
On peut calculer les racines d’un polynôme par la fonction roots.
Par exemple
 roots (P ) affiche 3 3
Example
Déterminer les racines du polynôme Q (x ) = x 3 + 2x − 3

On peut évaluer la valeur d’un polynôme en un point

quelconque a ∈ R par la fonction polyval. Par exemple :
 polyval (P, a) calcule la valeur de P (a).
Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 61 / 103
Compléments
Vecteurs : Sous forme polynôme

Example
Trouver la valeurs P (1) et Q (0) pour les fonctions P et Q
précédente.

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 62 / 103

Compléments
Vecteurs : Sous forme polynôme

Example
Trouver la valeurs P (1) et Q (0) pour les fonctions P et Q
précédente.

On peut aussi obtenir le polynôme "minimal" à partir de ses

racines. Cela se fait par la fonction poly. Par exemple, le
polynôme dont les racines sont 1, 2 et 3 s’obtient par :
 poly ([1, 2, 3]).

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 62 / 103

Compléments
Vecteurs : Sous forme polynôme

Example
Trouver la valeurs P (1) et Q (0) pour les fonctions P et Q
précédente.

On peut aussi obtenir le polynôme "minimal" à partir de ses

racines. Cela se fait par la fonction poly. Par exemple, le
polynôme dont les racines sont 1, 2 et 3 s’obtient par :
 poly ([1, 2, 3]). Il correspond à x 3 − 6x 2 + 11x − 6.

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 62 / 103

Compléments
Vecteurs : Sous forme polynôme

Example
Trouver la valeurs P (1) et Q (0) pour les fonctions P et Q
précédente.

On peut aussi obtenir le polynôme "minimal" à partir de ses

racines. Cela se fait par la fonction poly. Par exemple, le
polynôme dont les racines sont 1, 2 et 3 s’obtient par :
 poly ([1, 2, 3]). Il correspond à x 3 − 6x 2 + 11x − 6.
Le produit et la division euclidienne de deux polynômes
d’obtiennent respectivement en utilisant la fonction conv et
deconv . Par exemple le produit des fonctions P et Q s’obtient :
 conv (P, Q ) ce donne 1 − 6 11 − 15 36 − 27

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 62 / 103

Compléments
Complexes

Un nombre complexe défini dans matlab par z = a + b ∗ i.

Toutes les opérations algébriques sont possibles sur Matlab.
Considérons le complexe z = 2 + 3 ∗ i
 abs. Permet de calculer le module du complexe
 abs (z ) affiche 3, 6056

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 63 / 103

Compléments
Complexes

Un nombre complexe défini dans matlab par z = a + b ∗ i.

Toutes les opérations algébriques sont possibles sur Matlab.
Considérons le complexe z = 2 + 3 ∗ i
 abs. Permet de calculer le module du complexe
 abs (z ) affiche 3, 6056
 imag. Permet d’afficher la partie imaginaire du complexe
 imag (z ) affiche 3

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 63 / 103

Compléments
Complexes

Un nombre complexe défini dans matlab par z = a + b ∗ i.

Toutes les opérations algébriques sont possibles sur Matlab.
Considérons le complexe z = 2 + 3 ∗ i
 abs. Permet de calculer le module du complexe
 abs (z ) affiche 3, 6056
 imag. Permet d’afficher la partie imaginaire du complexe
 imag (z ) affiche 3
 real. Permet d’afficher la partie réelle du complexe.
 real (z ) affiche 2

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 63 / 103

Compléments
Complexes

Un nombre complexe défini dans matlab par z = a + b ∗ i.

Toutes les opérations algébriques sont possibles sur Matlab.
Considérons le complexe z = 2 + 3 ∗ i
 abs. Permet de calculer le module du complexe
 abs (z ) affiche 3, 6056
 imag. Permet d’afficher la partie imaginaire du complexe
 imag (z ) affiche 3
 real. Permet d’afficher la partie réelle du complexe.
 real (z ) affiche 2
 angle. Permet de calculer l’argument du complexe.
 angle (z ) affiche 0, 9828 (en radian)

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 63 / 103

Compléments
Figure

Avec Matlab on peut représenter plusieurs figures sur une seule

page en subdivisant la page par la fonction matlab subplot. Elle
s’exécute sous la forme subplot(a,b,c) où :
 a désigne le nombre de ligne souhaitée
 b désigne le nombre de colonne souhaitée
 c précise la position de la figure courante

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 64 / 103

Compléments
Figure

Avec Matlab on peut représenter plusieurs figures sur une seule

page en subdivisant la page par la fonction matlab subplot. Elle
s’exécute sous la forme subplot(a,b,c) où :
 a désigne le nombre de ligne souhaitée
 b désigne le nombre de colonne souhaitée
 c précise la position de la figure courante
 x = −π : 0.1 : π ;
 subplot (2, 1, 1); plot (x, sin(x ),0 b 0 )
 subplot (2, 1, 2); plot (x, cos(x ),0 r 0 )
Example
Représenter les fonctions suivantes : ln(x 2 ), cos (x ), sin(x ) et
ex sur la même page.
Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 64 / 103
 Programmation avec MATLAB

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 65 / 103

Programmation avec MATLAB
M-files ou scripts
Lorsqu’on a un programme comportant plusieurs instructions, il
faut le saisir dans un M-file. Cela permet de l’enregistrer et de
l’exécuter par la suite. Un script (ou M-file) est un fichier d’exten-
sion .m contenant des instructions MATLAB.
Il est important de bien ranger les fichiers MATLAB afin de ne pas
avoir de problème à l’exécution. Pour cela il faut créer des dos-
siers et des sous dossiers sur votre machine pour les différents
programmes.

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 66 / 103

Programmation avec MATLAB
M-files ou scripts
Lorsqu’on a un programme comportant plusieurs instructions, il
faut le saisir dans un M-file. Cela permet de l’enregistrer et de
l’exécuter par la suite. Un script (ou M-file) est un fichier d’exten-
sion .m contenant des instructions MATLAB.
Il est important de bien ranger les fichiers MATLAB afin de ne pas
avoir de problème à l’exécution. Pour cela il faut créer des dos-
siers et des sous dossiers sur votre machine pour les différents
programmes.
Example
Créer un M-file qui permet la représentation des quatre
fonctions suivantes :
1 1+x
f (x ) = x 2 , g (x ) = ln( ), h(x ) = sh(x ), T (x ) = th(x ) sur
2 1−x
une même page.
Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 66 / 103
Programmation avec MATLAB
Opérateurs

Dans la programmation on a besoin des opérateurs relationnels

et des opérateurs logiques. Les plus courants sont :
 Opérations relationnels :
 <,(strictement plus petit) ;
 >,(strictement plus grand) ;
 <=,(plus petit ou égal) ;
 >=, (plus grand ou égal) ;
 ==(égalité) ;
 ∼=(non égalité).
 Opérations logiques :
 &, (et) ;
 |, (ou ) ;
 ∼, (non ).

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 67 / 103

Programmation avec MATLAB
Instructions de contrôle

Comme les langages de programmation classiques, Matlab

dispose des commandes du type for, while et if.

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 68 / 103

Programmation avec MATLAB
Instructions de contrôle

Comme les langages de programmation classiques, Matlab

dispose des commandes du type for, while et if.
L’instruction if à deux syntaxes que l’on retrouve dans la
programmation :

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 68 / 103

Programmation avec MATLAB
Instructions de contrôle

Comme les langages de programmation classiques, Matlab

dispose des commandes du type for, while et if.
L’instruction if à deux syntaxes que l’on retrouve dans la
programmation :
Cas 1 :
if (expression)
commande 1 (si l’expression est vraie) ;
end

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 68 / 103

Programmation avec MATLAB
Instructions de contrôle

Comme les langages de programmation classiques, Matlab

dispose des commandes du type for, while et if.
L’instruction if à deux syntaxes que l’on retrouve dans la
programmation :
Cas 1 :
if (expression)
commande 1 (si l’expression est vraie) ;
end
Cas 2 :
if (expression)
commande 1 (si l’expression est vraie) ;
else
Commande 2 (si l’expression n’est pas vraie) ;
end.
Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 68 / 103
Programmation avec MATLAB
Instructions de contrôle

Example
Créer un M-file qui, à partir d’une moyenne donné, affiche
Ajourné ou Admis.

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 69 / 103

Programmation avec MATLAB
Instructions de contrôle

Example
Créer un M-file qui, à partir d’une moyenne donné, affiche
Ajourné ou Admis.

Code Matlab
m= 14 % moyenne compris entre 0 et 20
if (m>=10)
disp (’Admis en classe supérieur’) ;
end

Remarque
Comme il n’y a pas de commande pour m < 10 il y aura rien
d’affiché.
Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 69 / 103
Programmation avec MATLAB
Instructions de contrôle

Code Matlab
m= 14 % moyenne compris entre 0 et 20
if (m>=10)
disp (’Admis en classe supérieur’) ;
else
disp (’Ajourné’)
end

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 70 / 103

Programmation avec MATLAB
Instructions de contrôle

La boucle for
for (expression)
commande ;
end.

Example
Générer le vecteur formé des dix (10) premières valeurs
multiples de cinq (5).

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 71 / 103

Programmation avec MATLAB
Instructions de contrôle

Code Matlab
% initialisation du vecteur à l’ensemble vide
V=[ ] ;
for i=1 :inf
% test sur le reste de la division (i/5) doit être nul
if (rem(i,5)==0)
V=[V i] ;
end
if (length (V)==10)
break
end
end

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 72 / 103

Programmation avec MATLAB
Instructions de contrôle

L’instruction While
while (expression)
commande ;
end.

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 73 / 103

Programmation avec MATLAB
Instructions de contrôle

L’instruction While
while (expression)
commande ;
end.

Example
reprendre l’exemple précédent avec while.

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 73 / 103

Programmation avec MATLAB
Instructions de contrôle

Code Matlab
V=[ ] ;
i=0 ;
while (length (V)<10)
i=i+1 ;
if (rem(i,5)==0)
V=[V i] ;
end
end

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 74 / 103

Programmation avec MATLAB
Création de fonction

Une fonction est un fichier qui possède des paramètres d’appel

et des paramètres de retour. Beaucoup de commande Matlab
sont des fichiers fonctions. On peut créer sa propre fonction
(commande) sous Matlab. La syntaxe et la suivante :

Code Matlab
function ouput=nom(input)
les instructions de calcul
les instructions d’affichage
end

Example
Créer une fonction qui permet de générer le factoriel d’un
nombre.
Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 75 / 103
Programmation avec MATLAB
Création de fonction

Code Matlab
function fact=factorielle (n)
% tester si le nombre est bien définie pour le factoriel
if ∼ ((fix (n) == n)&(n >= 0))
error(’le nombre doit être un entier positif’)
end
% calculer la valeur du factoriel
fact=prod(1 :n) ;
end

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 76 / 103

Programmation avec MATLAB
Fonctions en ligne

Une fonction ne comportant qu’un petit nombre d’instructions

peut être définie directement dans une fenêtre de commandes
de la manière suivante :
  f=inline(’sin(x)+cos(y)’)
f =
inline function :
f (x, y ) = sin(x ) + cos(y )
π
  f ( , 0)
2
ans =
2

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 77 / 103

Programmation avec MATLAB
Fonctions en ligne

Les arguments de la fonction f sont normalement fournis à

l’appel dans l’ordre d’apparition dans la définition de la fonction.
On peut aussi spécifier les arguments d’appel explicitement.
  f=inline(’sin(x)+cos(y)’,’x’,’y’)
f =
inline function :
f (x, y ) = sin(x ) + cos(y )
  f (90, 0)
ans =
2
On peut aussi ne pas utiliser l’instruction inline.
  f=’sin(x)+cos(y)’
dans ce cas à l’exécution on fait :
 x = 90; y = 0; eval (f )
ans =
Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 78 / 103
 Optimisation

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 79 / 103

Optimisation
Création des fonctions numériques

On peut utiliser le caractère @ Pour créer les fonctions

numériques

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 80 / 103

Optimisation
Création des fonctions numériques

On peut utiliser le caractère @ Pour créer les fonctions

numériques
Example
  f = @(x )2x 2 − 5x + 3
  f (2) on obtient f (2) = 1

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 80 / 103

Optimisation
Création des fonctions numériques

On peut utiliser le caractère @ Pour créer les fonctions

numériques
Example
  f = @(x )2x 2 − 5x + 3
  f (2) on obtient f (2) = 1

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 80 / 103

Optimisation
Création des fonctions numériques

On peut utiliser le caractère @ Pour créer les fonctions

numériques
Example
  f = @(x )2x 2 − 5x + 3
  f (2) on obtient f (2) = 1

En dimension n (n ≥ 2)la variable x est considéré comme un

vecteur.
Exercice
  g = @(x, y )(x )2 + (y )2 + 2
  g (1, 2) on obtient g (1, 2) = 7

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 80 / 103

Optimisation
Recherche d’un extremum : minimum

Nous utilisons ici la fonction fminsearch est dans la bibliothèque

MATLAB. Cette fonction correspond la méthode Simplexe de
Nelder Meald qui n’est pas à confondre à celle de Dantzig

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 81 / 103

Optimisation
Recherche d’un extremum : minimum

Nous utilisons ici la fonction fminsearch est dans la bibliothèque

MATLAB. Cette fonction correspond la méthode Simplexe de
Nelder Meald qui n’est pas à confondre à celle de Dantzig
Example
  f = @(x )2x 2 − 5x + 3
  [x, fval ] = fminsearch(f , 0)

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 81 / 103

Optimisation
Recherche d’un extremum : minimum

Nous utilisons ici la fonction fminsearch est dans la bibliothèque

MATLAB. Cette fonction correspond la méthode Simplexe de
Nelder Meald qui n’est pas à confondre à celle de Dantzig
Example
  f = @(x )2x 2 − 5x + 3
  [x, fval ] = fminsearch(f , 0)

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 81 / 103

Optimisation
Recherche d’un extremum : minimum

Nous utilisons ici la fonction fminsearch est dans la bibliothèque

MATLAB. Cette fonction correspond la méthode Simplexe de
Nelder Meald qui n’est pas à confondre à celle de Dantzig
Example
  f = @(x )2x 2 − 5x + 3
  [x, fval ] = fminsearch(f , 0)

En dimension n (n ≥ 2).

Exercice
  g = @(x )(x (1))2 + (x (2))2 + 2
  [x, fval ] = fminsearch(g, [1, 1])

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 81 / 103

Optimisation
Projet

Exercice
On considère une fonction numérique f et on définit un
opérateur :
n1 o
J (x ) = argmin f (x ) + | (x − x0 | , avec λ et x0 quelconque
2λ
On construit l’algorithme numérique suivant :
n 1 o
J (x k +1 ) = argmin f (x k ) + | (x k − x0 | , avec λ et x0 quelcon
2λk

Travail à faire :
pour x0 fixer et pour quelques valeurs de λk .

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 82 / 103

 Résolution numérique des EDO et des EDP

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 83 / 103

Résolution numérique des EDO
Formulation

Definition
On appelle Équation Différentielle Ordinaire (EDO) toute
équation pouvant se mettre sous la forme : y 0 = f (t, y ) où f est
une fonction continue.
 0 Elle dite de Cauchy si on peut mettre
y = f (t, y ), t ∈ [a, b ]
sous la forme :
y (a) = y0

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 84 / 103

Résolution numérique des EDO
Formulation

Definition
On appelle Équation Différentielle Ordinaire (EDO) toute
équation pouvant se mettre sous la forme : y 0 = f (t, y ) où f est
une fonction continue.
 0 Elle dite de Cauchy si on peut mettre
y = f (t, y ), t ∈ [a, b ]
sous la forme :
y (a) = y0

Pour obtenir une approximation numérique de la solution y (t )

sur l’intervalle [a, b ], on estime la valeur de la fonction en un
nombre fini de points ti , i = 0, n constituant les nœuds d’un
maillage de l’intervalle [a, b ].

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 84 / 103

Résolution numérique des EDO
Résolution

L’écart entre deux points de la discrétisation est noté h et est

appelé pas de discrétisation. Dans les méthodes dites simples ce
pas h est constant et on a : hn = tn − tn−1 , n = 1, N. La solution
numérique discrète obtenue aux points tn est notée yn = y (tn ).

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 85 / 103

Résolution numérique des EDO
Résolution

L’écart entre deux points de la discrétisation est noté h et est

appelé pas de discrétisation. Dans les méthodes dites simples ce
pas h est constant et on a : hn = tn − tn−1 , n = 1, N. La solution
numérique discrète obtenue aux points tn est notée yn = y (tn ).
En réécrivant l’EDO en chaque point de discrétisation tn et on
obtient :  0
y (tn ) = f (tn , y (tn )), tn ∈ [a, b ]
y (a) = y (t0 ) = y0
En utilisant la formule décentrée à droite (Euler progressif) de la
y ( tn + 1 ) − y ( tn )
dérivée qui est : y 0 (tn ) =
tn+1 − tn

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 85 / 103

Résolution numérique des EDO
Résolution

L’écart entre deux points de la discrétisation est noté h et est

appelé pas de discrétisation. Dans les méthodes dites simples ce
pas h est constant et on a : hn = tn − tn−1 , n = 1, N. La solution
numérique discrète obtenue aux points tn est notée yn = y (tn ).
En réécrivant l’EDO en chaque point de discrétisation tn et on
obtient :  0
y (tn ) = f (tn , y (tn )), tn ∈ [a, b ]
y (a) = y (t0 ) = y0
En utilisant la formule décentrée à droite (Euler progressif) de la
y ( tn + 1 ) − y ( tn )
dérivée qui est : y 0 (tn ) =
tn+1 − tn
0 yn+1 − yn
En posant : y (tn ) =
h

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 85 / 103

Résolution numérique des EDO
Résolution

On obtient le schéma numérique dit d’Euler progressif ou

explicite de la forme :

yn+1 = yn + h.f (tn , yn ), tn ∈ [a, b ]
y (a) = y0

Remarque
On obtient ainsi une suite récurrente d’ordre 1 dont on peut
calculer tous les termes de façon consécutive.

Example
y 0 = 3t − 3y , t ∈ [0, 100]

Pour l’EDO suivante : Trouver le
y (0) = 1
schéma numérique puis représenter sa solution.
Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 86 / 103
Résolution numérique des EDO
Code Matlab

Code Matlab

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 87 / 103

Résolution numérique des EDO
Résolution

Example
Considérons l’équation différentielle

 y (0) = 0
y 0 (0) = 5
y − y 0 + 6y = xe −2x , x ∈ [0, +∞[
 00

1 Déterminer le schémas numérique de cette équation.

2 Représenter cette solution numérique.

Remarque
yn+1 − 2yn + yn−1
La dérivée centrée d’ordre 2 s’écrit :
h2
Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 88 / 103
Résolution numérique des EDO
Code Matlab

Code Matlab

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 89 / 103

Résolution numérique des EDP
Formulation

Definition
Les équations aux dérivées partielles (EDP) du second ordre
sont de la forme :
∂2 U ∂2 U ∂2 U ∂U ∂U

a(x, y ) 2
+ b ( x, y ) + c ( x, y ) 2
+ f x, y , , =0
∂x ∂x ∂y ∂y ∂x ∂y

En posant ∆ = b2 − 4ac on a :
∆ > 0, nous avons une EDP hyperbolique (équation des
∂2 U ∂2 U
ondes 2 = c 2 2 ).
∂t ∂x

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 90 / 103

Résolution numérique des EDP
Formulation

Definition
Les équations aux dérivées partielles (EDP) du second ordre
sont de la forme :
∂2 U ∂2 U ∂2 U ∂U ∂U

a(x, y ) 2
+ b ( x, y ) + c ( x, y ) 2
+ f x, y , , =0
∂x ∂x ∂y ∂y ∂x ∂y

En posant ∆ = b2 − 4ac on a :
∆ > 0, nous avons une EDP hyperbolique (équation des
∂2 U ∂2 U
ondes 2 = c 2 2 ).
∂t ∂x
∆ = 0, nous avons une EDP parabolique (diffusion).

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 90 / 103

Résolution numérique des EDP
Formulation

Definition
Les équations aux dérivées partielles (EDP) du second ordre
sont de la forme :
∂2 U ∂2 U ∂2 U ∂U ∂U

a(x, y ) 2
+ b ( x, y ) + c ( x, y ) 2
+ f x, y , , =0
∂x ∂x ∂y ∂y ∂x ∂y

En posant ∆ = b2 − 4ac on a :
∆ > 0, nous avons une EDP hyperbolique (équation des
∂2 U ∂2 U
ondes 2 = c 2 2 ).
∂t ∂x
∆ = 0, nous avons une EDP parabolique (diffusion).
∆ < 0, nous avons une EDP elliptique (Laplace, Poisson).

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 90 / 103

Résolution numérique des EDP
Exemple didactique
Considérons l’équation de la chaleur suivante :

∂U ∂2 U
, x ∈]0, π [, t > 0


 =
∂t ∂x 2

 U ( x, 0 ) = x −1
 U (0, t ) = U (π, t ) = 0

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 91 / 103

Résolution numérique des EDP
Exemple didactique
Considérons l’équation de la chaleur suivante :

∂U ∂2 U
, x ∈]0, π [, t > 0


 =
∂t ∂x 2

 U ( x, 0 ) = x −1
 U (0, t ) = U (π, t ) = 0

Pour obtenir le Schéma numérique de cette EDP on discrétise

sont domaine (espace, temps) par la méthode des différences
finies dont le principe est le suivant :

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 91 / 103

Résolution numérique des EDP
Exemple didactique
Considérons l’équation de la chaleur suivante :

∂U ∂2 U
, x ∈]0, π [, t > 0


 =
∂t ∂x 2

 U ( x, 0 ) = x −1
 U (0, t ) = U (π, t ) = 0

Pour obtenir le Schéma numérique de cette EDP on discrétise

sont domaine (espace, temps) par la méthode des différences
finies dont le principe est le suivant : On se propose d’approcher
U en xj = j∆x, j = 0, K et tn = n∆t, n = 0, N ou ∆x et ∆t sont
positifs et appelés à tendre vers 0.

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 91 / 103

Résolution numérique des EDP
Exemple didactique
Considérons l’équation de la chaleur suivante :

∂U ∂2 U
, x ∈]0, π [, t > 0


 =
∂t ∂x 2

 U ( x, 0 ) = x −1
 U (0, t ) = U (π, t ) = 0

Pour obtenir le Schéma numérique de cette EDP on discrétise

sont domaine (espace, temps) par la méthode des différences
finies dont le principe est le suivant : On se propose d’approcher
U en xj = j∆x, j = 0, K et tn = n∆t, n = 0, N ou ∆x et ∆t sont
positifs et appelés à tendre vers 0. Posons

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 91 / 103

Résolution numérique des EDP
Exemple didactique
Considérons l’équation de la chaleur suivante :

∂U ∂2 U
, x ∈]0, π [, t > 0


 =
∂t ∂x 2

 U ( x, 0 ) = x −1
 U (0, t ) = U (π, t ) = 0

Pour obtenir le Schéma numérique de cette EDP on discrétise

sont domaine (espace, temps) par la méthode des différences
finies dont le principe est le suivant : On se propose d’approcher
U en xj = j∆x, j = 0, K et tn = n∆t, n = 0, N ou ∆x et ∆t sont
positifs et appelés à tendre vers 0. Posons
∂U U (xj , tn+1 ) − U (xj , tn )
( x , tn ) =
∂t j ∆t
∂2 U U (xj +1 , tn ) − 2U (xj , tn ) + U (xj −1 , tn )
et (xj , tn ) =
∂x 2 (∆x )2
Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 91 / 103
Résolution numérique des EDP
Exemple didactique

Ainsi l’équation de la chaleur peut s’écrire :

∆t h i
U (xj , tn+1 ) = U (xj , tn ) + ∗ U ( xj + 1 , t n ) − 2U ( xj , t n ) + U ( x j − 1 , tn )
(∆x )2

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 92 / 103

Résolution numérique des EDP
Exemple didactique

Ainsi l’équation de la chaleur peut s’écrire :

∆t h i
U (xj , tn+1 ) = U (xj , tn ) + ∗ U ( xj + 1 , t n ) − 2U ( xj , t n ) + U ( x j − 1 , tn )
(∆x )2

En posant U (xj , tn ) = Ujn on obtient :

∆t h i
Ujn+1 = Ujn + ∗ U n
j + 1 − 2U n
j + U n
j − 1
(∆x )2

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 92 / 103

Résolution numérique des EDP
Exemple didactique

Ainsi l’équation de la chaleur peut s’écrire :

∆t h i
U (xj , tn+1 ) = U (xj , tn ) + ∗ U ( xj + 1 , t n ) − 2U ( xj , t n ) + U ( x j − 1 , tn )
(∆x )2

En posant U (xj , tn ) = Ujn on obtient :

∆t h i
Ujn+1 = Ujn + ∗ U n
j + 1 − 2U n
j + U n
j − 1
(∆x )2

Les conditions aux bords donnent :

Uj0 = xj − 1
 
U ( xj , 0 ) = xj − 1
⇒
U (0, tn ) = U (π, tn ) = 0 U0n = UKn = 0

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 92 / 103

Résolution numérique des EDP
Résolution
D’où le schéma numérique suivant :

n + ∆t
 h i
n +1 n − 2U n + U n
 U = U ∗ U j −1
 j j
(∆x )2 j +1 j

0
U = xj − 1
 jn


U0 = UKn = 0

En faisant varier le j de 1 à K et en posant X = (U1 , U2 , · · · , UK )0 l’EDP

s’écrit :
 X n = I + ∆t ∗ A n ∗ X 0
  
K
(∆x 2 )
X = (U10 , U20 , · · · , UK0 )0
 0

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 93 / 103

Résolution numérique des EDP
Résolution
D’où le schéma numérique suivant :

n + ∆t
 h i
n +1 n − 2U n + U n
 U = U ∗ U j −1
 j j
(∆x )2 j +1 j

0
U = xj − 1
 jn


U0 = UKn = 0

En faisant varier le j de 1 à K et en posant X = (U1 , U2 , · · · , UK )0 l’EDP

s’écrit :
 X n = I + ∆t ∗ A n ∗ X 0
  
K
(∆x 2 )
X = (U10 , U20 , · · · , UK0 )0
 0

Remarque
A est une matrice symétrique d’ordre K . C’est une matrice
bande de trois parallèles dont la diagonale est -2 et le autres 1.
Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 93 / 103
Résolution numérique des EDP
Code Matlab

Code Matlab

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 94 / 103

Résolution numérique des EDP
Code Matlab

Code Matlab

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 95 / 103

Exercices
Polynômes et systèmes

Exercice
Déterminer les racines des polynômes suivants :
1 P (x ) = x 2 + 3x + 8.
2 Q (x ) = x 5 + 2x 4 − 2x 3 + 4x 2 + 3x + 5.

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 96 / 103

Exercices
Polynômes et systèmes

Exercice
Déterminer les racines des polynômes suivants :
1 P (x ) = x 2 + 3x + 8.
2 Q (x ) = x 5 + 2x 4 − 2x 3 + 4x 2 + 3x + 5.

Exercice
Construire les polynômes à partir des racines suivantes :
1 2 et 1.
2 2, 2, 3 et −5.

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 96 / 103

Exercices
Polynômes et systèmes

Exercice
Déterminer les racines des polynômes suivants :
1 P (x ) = x 2 + 3x + 8.
2 Q (x ) = x 5 + 2x 4 − 2x 3 + 4x 2 + 3x + 5.

Exercice
Construire les polynômes à partir des racines suivantes :
1 2 et 1.
2 2, 2, 3 et −5.

Exercice
6
Décomposer en éléments simples f (x ) =
x 4 + 6x 3 + 11x 2 + 6x
Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 96 / 103
Exercices
Polynômes et systèmes

Exercice
Résoudre les systèmes linéaires suivantes :
 2  p
x =1 x |x | + cos(xy ) = 2
et
x − y = −1 x + sin(y ) = 1

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 97 / 103

Exercices
Polynômes et systèmes

Exercice
Résoudre les systèmes linéaires suivantes :
 2  p
x =1 x |x | + cos(xy ) = 2
et
x − y = −1 x + sin(y ) = 1

Exercice
Résoudre
 les équations suivantes :
 2x1 + x2 = −5
4x − 3x2 + 2x3 = 0
 1
x1 + 2x2 − x3 = 1

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 97 / 103

Exercices
Vecteurs

Exercice
On considère les nombres complexes u = 11 − 7i et
v = −1 + 3i
1 Calculer les modules de u et de v , les produits uv et v u
2 Déterminer la partie imaginaire de u 3 + v 2

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 98 / 103

Exercices
Vecteurs

Exercice
On considère les nombres complexes u = 11 − 7i et
v = −1 + 3i
1 Calculer les modules de u et de v , les produits uv et v u
2 Déterminer la partie imaginaire de u 3 + v 2

Exercice
On note u, v et w trois vecteurs de R3 par :
u = (1, −1, 2)t , v = (10, −1, 3)t et w = (5, −1, 4)t
1 Calculer : 3u, ||u ||2 , 2u − v + 5w, ||2u − v + 5w ||1 et
||w − 4v ||∞
2 déterminer l’angle formé par les vecteurs v et w.
Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 98 / 103
Exercices
Matrices

Exercice
On note
 A, B et C les
 trois matrices suivantes :
1 3 2    
 −5 3 1  1 −2 5 10 −5
A=   , B= et C =
−10 0 3  6 1 −1 3 1
1 0 −2
1 Calculer les matrices AB, BA, AB t et D = I2 − BB t
2 Calculer les déterminants de : A, B, C, D, et E = AAt .
Trouver leurs inverses.
3 Déterminer les valeurs propres de E et préciser le rayon
spectral.
4 Déterminer les valeurs propres de E.
5 Déterminer D et P telle que E = PDP −1 .
Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 99 / 103
Exercices
Graphiques

Exercice
Considérons
 les fonctions suivantes
x4

f (x1 , x2 ) = 4 − 2.1x12 + 1 x12 + x1 x2 + − 4 + 4x22 x22


3
x = (x1 , x2 ) ∈ [−5, 5]2

x12 x22

x
− cos (x1 ) + 1 cos ( √2 ) + 1



g (x1 , x2 ) = +

50 50 2
x = (x1 , x2 ) ∈ [−15, 15]2


1 Représenter les fonctions f et g.

2 Resoudre les équations f (x1 , x2 ) = 0, g (x1 , x2 ) = 0 et
f ( x1 , x2 ) = g ( x1 , x2 ) = 0

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 100 / 103

Exercices
Résolution de système et d’EDO

Exercice
Résoudre
 les équations suivantes :
 2x1 + x2 = −5
4x − 3x2 + 2x3 = 0
 1
x1 + 2x2 − x3 = 1

Exercice
Résoudre les équations différentielles suivantes :
1 y 0 = sh(y ), y (1) = 1 ;
2 y 00 + 2y 0 + y = tet , y (0) = 0 ;
 0
3
x = −y + sin(at )
y 0 = x − cos(at ), a ∈ R

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 101 / 103

Exercices
Résolution d’EDP

Exercice
Considérons l’EDP suivante :

∂2 U ∂2 U
(x, t ) = 0, (x, t ) ∈ [0, π ] × [0, 1];


 ( x, t ) +
∂t 2 ∂x 2

 U (x, 0) = 0, U (x, 1) = 1
 U (0, t ) = 0, U (π, t ) = 1, t ≥ 0.

1 Donner un schéma numérique de cette EDP.

2 Faire la simulation numérique de la solution de l’EDP
ci-dessus.

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 102 / 103

 Merci de votre aimable attention

Master de Physique (UFR-ST) Dr Kounhinir SOME Avril 2022 103 / 103