Cours Matlab - Partie 1

UNIVERSITÉ CHOUAIB DOUKKALI
Ecole Nationale des Sciences Appliquées d’El Jadida-ENSAJ
Cours : MATLAB
Filière : GEE/GI S5
Pr: M. SAYYOURI
Année universitaire 2017-2018

1
PLAN
Evaluation:
 Introduction
❖ Examen : 100%
 Calculs élémentaires
P  Programmation sous Matlab
L  Graphiques 2D et 3D
 Recherche de zéro d’une fonction
A  Résolution numériques des équations différentielles
N  Initiation à Simulink
 Initiation à Stateflow
2
Introduction
Introduction
Matlab (MATrix LABoratory) est l’une des principales plateformes de développement.
Calculs élémentaires
Graphiques 2D et 3D
Matlab est un système interactif de calcul numérique et de visualisation graphique.
Recherche de zéro Il possède un langage de programmation simple et efficace, constitué de :
d’une fonction
➢ de nombreux outils du calcul matriciel,
Résolution numérique
➢ des principaux algorithmes de l’analyse numérique (schémas d’intégration, différentiation...),
d’équations
➢ des outils pour le traitement du signal,
différentielles
Initiation à Simulink
➢ des outils de visualisation en 2D et 3D,
Initiation à Stateflow ➢ l’outil Simulink, pour modéliser et simuler des systèmes en utilisant des schémas-blocs,
➢ des outils pour concevoir des interfaces utilisateur graphiques (GUI pour « Graphical User Interface »)
personnalisées,
➢ la prise en charge d’échange de données et de signaux avec plusieurs type de matériels (oscilloscope,
carte d’acquisition, carte de contrôle, etc.),
➢ l’intégration d’applications externes écrites en C, C++, Fortran, Java, Excel, etc.
➢ la conversion d'applications MATLAB en C et C++ . 3
M.SAYYOURI
Introduction
Introduction
Graphiques 2D et 3D
Recherche de zéro d’une
fonction
d’équations
différentielles
La fenêtre de commande (ou ligne de
Initiation à Stateflow
commande):
c’est l’outil de base de Matlab. Elle permet entre
autre de faire n’importe quelle opération, de définir
et affecter les variables d’environnement, d’utiliser
les toolboxes via leurs fonctions, …
4
M.SAYYOURI
Introduction
Introduction
Graphiques 2D et 3D
fonction
d’ équations
différentielles
La fenêtre d’historique de commande :
Initiation à Stateflow permet de répéter les commandes précédentes en
double-cliquant sur la commande concernée dans
cette fenêtre. On peut obtenir le même résultat en
ligne de commande par l’appui répétitif sur la
touche « flèche haut » du pavé numérique.
5
M.SAYYOURI
Introduction
Introduction
Graphiques 2D et 3D
fonction
d’équations
différentielles
La fenêtre d’espace de travail :
Cette fenêtre permet de visualiser les différentes
variables existantes, en affichant leur nom, leur
contenu (ou leur dimension) ainsi que le min et le
max des données qu’elle contient. On peut obtenir
la liste des variables ainsi que leurs dimensions en
ligne de commande en tapant la commande
« whos ».
6
M.SAYYOURI
Introduction
Introduction
Graphiques 2D et 3D
fonction
d’équations
différentielles
Initiation à Simulink La fenêtre de répertoire courant :

Initiation à Stateflow Cette fenêtre permet de visualiser le répertoire
de travail dans lequel vous êtes. C’est ici que
sont chargés / enregistrés les scripts, les modèles
ou les données.
7
M.SAYYOURI
Introduction Matlab comme Calculatrice
➢ Le symbole [>>] indique à l’utilisateur où il faut rentrer la commande
Graphiques 2D et 3D
➢ Commençons par les operateurs les plus courants : +, -, *, /, ^. Le dernier signifie «puissance».
Les parenthèses s’utilisent de manière classique.
fonction
➢ Un dernier point sur les variables :
d’équations
▪ MATLAB fait la différence entre les minuscules et les majuscules.
différentielles
▪ Les noms de variables peuvent avoir une longueur quelconque.
▪ Les noms de variables doivent commencer par une lettre.
Initiation à Stateflow ➢ Matlab comporte un très grand nombre d’opérateurs, de commandes et de fonctions. On peut
taper à titre d’exemple les commandes suivantes :
• help permet d’obtenir l’aide de l’aide et donne une liste thématique ;
• help nom de fonction donne la définition de la fonction désignée;
• lookfor sujet donne une liste des rubriques de l’aide en ligne en relation avec le sujet
indiqué,
• clear permet d’effacer une partie ou toutes les variables définies jusqu’à Présent,
• who affiche les noms de toutes les variables en cours,
• whos affiche toutes les variables avec les détails..
M.SAYYOURI
Introduction
Dans Matlab, on trouve des variables prédéfinies :
• ans qui contient le dernier résultat de calcul
Graphiques 2D et 3D
• i et j qui représentent racine −1 .
• Pi qui représente 𝜋
fonction
• eps qui est un nombre prédéfini.
• inf qui représente plus infini
d’équations
• NaN : pas un nombre (not a number)
différentielles
Initiation à Simulink Attention : Ces variables ne sont pas protégées, donc si vous les affectez, elles ne gardent pas leur
Initiation à Stateflow valeur. C’est souvent le problème pour i et j que l’on utilise souvent spontanément comme indices
de boucles, de telle sorte qu’on ne peut plus ensuite définir de complexe ! !
 Toutes les fonctions courantes et moins courantes existent dans Matlab. On retiendra que pour
appliquer une fonction à une valeur, il faut mettre cette dernière entre parenthèses.
 Voici une liste non exhaustive :
• Fonctions trigonométriques et inverses : sin, cos, tan, asin, acos, atan
• Fonctions hyperboliques (on rajoute «h») : sinh, cosh, tanh, asinh, acosh, atanh
• Racine, logarithmes et exponentielles : sqrt, log, log10, exp 9
M.SAYYOURI
Introduction • La notion de fonction est plus générale dans Matlab, et certaines fonctions peuvent avoir
Calculs élémentaires plusieurs entrées (comme besselj par exemple) mais aussi plusieurs sorties.
Graphiques 2D et 3D

 Dans Matlab, toutes les variables représentent des matrices. On peut considérer donc qu’il y a un seul type
principal de données : le type matrice.
fonction
 Les éléments d’une matrice peuvent être :
➢ des nombres (réels ou complexes),
d’équations
➢ ou des caractères.
différentielles
Initiation à Simulink  On utilise les crochets [ et ] pour définir le début et la fin de la matrice. Ainsi pour définir une variable M
1 2 3
contenant la matrice 4 5 6 , on écrira: >> % 3 écriture d'une matrice
7 8 9 >> M=[1 2 3
456
7 8 9]
>> M=[1 2 3; 4 5 6 ; 7 8 9]
>> M=[1, 2 ,3 ; 4, 5, 6 ; 7, 8, 9]
 Matlab propose des commandes pour créer certaines matrices particulières très simplement. Pour plus
d’information, lire le help de chaque fonction.
10
M.SAYYOURI
Introduction
Commande Description
ones(n,m) Matrice de taille n × m ne contenant que des 1.
Graphiques 2D et 3D
zeros(n,m) Matrice de taille n × m ne contenant que des 0.
eye(n,m) Matrice de taille n × m contenant des 1 sur la première diagonale et des 0 ailleurs.
fonction
diag(v) Matrice diagonale où les éléments de la diagonale sont les composantes du vecteur v.
rand(n,m) Matrice de taille n × m contenant des nombres aléatoires
d’équations
différentielles
 Matlab dispose également de moyens très simples pour créer des listes. La commande [a:p:b] crée une
liste dont les éléments sont: a, a + p, a + 2*p, . . . , a + n*p, où n ∈ N,
a+n*p ≤ b < a+(n+1)*p
 Le cas particulier [a:b] est un raccourci pour [a:1:b]. Si les conditions initiales sont erronées,
Matlab renvoie un message d’erreur.
 linspace : Un autre cas particulier de [a:p:b] est la fonction linspace(a,b,n). Celle-ci crée une
liste de n éléments uniformément répartis entre a et b. Autrement dit, linspace(a,b,n) est la
même chose que [a: (b−a)/(n−1):b].
 Une dernière méthode pour créer des tableaux est la concaténation. Si A et B sont deux tableaux, alors
[A B], ou [A,B] est le tableau obtenu en collant B à la droite de A, et [A ;B] est le tableau obtenu en
collant B au-dessous de A. Comme d’habitude, il faut faire attention aux tailles de A et de11B.
M.SAYYOURI
Introduction
 Opérations avec les matrices
Calculs élémentaires  Matlab permet de faire certaines opérations avec des matrices. Dans ce qui suit, A et B sont des tableaux et
Graphiques 2D et 3D c est un scalaire.
Recherche de zéro d’une Commande Description
fonction A+B Addition terme à terme ; A et B doivent avoir le même format.
A+c = c+A Addition de c aux éléments de A.
A-B Soustraction terme à terme ; A et B doivent avoir le même format.
différentielles
A-c Soustraction de c aux éléments de A.
Initiation à Simulink A*B Produit matriciel standard ; nb. col. A doit être le même que nb de ligne de B
Initiation à Stateflow A.*B Multiplication terme à terme ; A et B doivent avoir le même format.
Aˆn (n ∈ 𝒁+ ) A * A * . . . * A (n fois) ; A doit être carrée.
𝑏𝑖𝑗
A.ˆB Tableau dont les éléments sont 𝑎𝑖𝑗
A’ Transposition et conjugaison.
B/A Le résultat est un tableau X tel que XA = B. Si A est inversible, alors X = B 𝐴−1 ;
nb. col. A doit être le même que nb. col. B.
A\B Le résultat est un tableau X tel que AX = B. Si A est inversible, alors X = 𝐴−1 B
; nb. Ligne de A doit être le même que nb. Ligne de B.
A./B Division terme à terme des éléments de A par ceux de B ; A et B doivent avoir le
même format.
A.\B Division terme à terme des éléments de B par ceux de A ; A et B doivent avoir le
même format. 12
M.SAYYOURI
Introduction
Fonctions sur les matrices
 Etant donnée une matrice A, il y a un certain nombre de choses que l’on peut calculer en rapport
Graphiques 2D et 3D
avec A. Nous présentons ici quelques fonctions définies dans Matlab prenant comme paramètre
des tableaux.
fonction Commande Description
Résolution numérique det(A) Renvoie le déterminant de A ; celle-ci doit être carrée.
d’équations trace(A) Renvoie la trace de A.
différentielles rank(A) Renvoie le rang de A
diag(A) Renvoie la première diagonale de A.
norm(v) Renvoie la norme euclidienne de v ; v est un vecteur. Il est aussi possible de
calculer d’autres normes ;
mean(A) Renvoie une liste contenant la moyenne des éléments de chaque colonne.
sum(A) Renvoie une liste contenant la somme des éléments de chaque colonne.
prod(A) Renvoie une liste contenant le produit des éléments de chaque colonne.
max(A) Renvoie une liste contenant la valeur maximale chaque colonne.
min(A) Renvoie une liste contenant la valeur minimale de chaque colonne.
length(A) Renvoie le maximum entre le nombre de lignes et de colonnes
 Finalement, on précise que toutes les fonctions mathématiques classiques (cos, sin, log, exp, etc)
s’appliquent également aux tableaux. Le résultat est un tableau où l’on a appliqué terme
13 à terme la
M.SAYYOURI fonction en question.
Exercices d’application
Introduction Exercice 1   2
3
11 , e , 9 , ln(5) , log (5) , sin( ) , 1   , 1  2i  (1  3i )

2 1.25

2
Calculs élémentaires 1. Calculer sous Matlab 10
4 2
 
Graphiques 2D et 3D
2. Afficher les résultats précédents sous format long , long e , short e, bank, et rational.
fonction Exercice 2
1. Créer le vecteur colonne
d’équations
différentielles 2. Saisir la matrice

a) Créer un vecteur colonne à 3 éléments, appelé ua; contenant les éléments de la troisième colonne de A:
b) Créer un vecteur colonne à 5 éléments, appelé ub; contenant les éléments de la deuxième ligne de A:
c) Créer un vecteur colonne à 9 éléments, appelé uc; formé par les éléments de la première, de la troisième et
de la cinquième colonne de A:
d) Créer un vecteur colonne à 10 éléments, appelé ud; formé par la première et la deuxième ligne de A:
14
M.SAYYOURI
Programmation sous Matlab
Introduction
 Un script Matlab est composé d'une suite d'instructions ou de commandes, toutes séparées par une virgule
Programmation sous
Matlab
(,), un point-virgule (;) ou un retour a la ligne ().
Graphiques 2D et 3D  Matlab facilite l’édition d’un nouveau script en utilisant le bouton (New Script) qui fait appel a l’éditeur de
Matlab.
fonction
 L’éditeur de Matlab crée par défaut un nouveau script sous le nom ‘untitled.m’. Le programmeur peut
d’équations spécifié ultérieurement un nom de son choix en utilisant le bouton (Save).
différentielles
 Généralement, le script est sauvegardé dans un fichier avec l’extension ‘.m’ (ex : script1.m) sous le dossier
Initiation à Stateflow courant. On peut aussi charger un fichier déjà édite en utilisant le bouton (Open).
 L’exécution du script se fait a travers l’invite de commande en spécifiant juste le nom du script.
 L’exécution peut être effectuée directement depuis l’éditeur en cliquant sur le bouton (Run). Néanmoins,
les sorties des scripts (affichage des résultats) se fait toujours sur la fenêtre des commandes.
15
M.SAYYOURI
Introduction
Affichage
Programmation sous
 L’affichage des variables se fait simplement par la spécification de leurs noms sans point virgule à la fin de
Matlab
Graphiques 2D et 3D la ligne.
Recherche de zéro d’une  La commande disp(X) affiche une variable X (chaine de caractères, entier, …) sur la fenêtre des
fonction
commandes. La commande disp n’accepte qu’un seul type de donnée à la fois.
 La commande fprintf() accepte tous genre de variable en utilisant leurs formats de types : %d pour les
d’équations
différentielles entiers, %f pour les réels, %s pour les chaines de caractères, %c pour le caractères, … . On peut ajouter un
Initiation à Simulink retour a la ligne via ‘\n’ et une tabulation via ‘\t’.
 La commande sprintf() sauvegarde tous genre de variables de la même manière que fprintf, mais sans
affichage sur la fenêtre des commandes. sprintf() retourne une chaine de caractères qui peut être affichée
ultérieurement en utilisant la commande disp().
Introduction des données à partir du clavier :

 La commande input(msg) fait introduire une valeur à partir du clavier en affichant d’abord un message msg.
16
M.SAYYOURI
Introduction
clear all
clc;
Programmation sous % lecture de deux valeurs
Matlab a = input ('Entrer la valeur de a : ');
b = input ('Entrer la valeur de b : ');
Graphiques 2D et 3D
% somme
Recherche de zéro d’une c = a + b;
% affichage
fonction
fprintf('%f + %f = %f \n', a, b, c);
Résolution numérique % introduction d'un vecteur
d’équations
vec = input('entrer un vecteur de 5 éléments : ');
% introduction d'une matrice
différentielles mat = input('entrer une matrice de 3x3 éléments : ');
Initiation à Simulink vec
mat
Entrer la valeur de a : -2.5

Entrer la valeur de b : 3.18
-2.500000 + 3.180000 = 0.680000
entrer un vecteur de 5 éléments : zeros(1,5)
entrer une matrice de 3x3 éléments : randi(3,3)
vec =
0 0 0 0 0
mat =
3 3 1
3 2 2 17
M.SAYYOURI 1 1 3
Introduction Les instructions d’alternative (IF et SWITCH):
Programmation sous  Comme dans la majorité des langages de programmation, Matlab permet d’utiliser l’instruction if selon
Matlab la syntaxe suivante :
Graphiques 2D et 3D
if condition
Commandes
elseif condition
fonction
Commandes
…
d’équations
else
différentielles Commandes
Initiation à Simulink end ;
Initiation à Stateflow  Les segments de elseif est else sont optionnels, par contre la parution de if et end est obligatoire
 Exercice
1. Écrire un script Matlab eq2degre.m qui permet de résoudre une équation de second degré à une seule
variable ,
2. Écrire un script Matlab qui permet de tester la parité d’un nombre entier naturel donné par l’utilisateur
18
M.SAYYOURI
Introduction %équation de second degré
a=input('entrer a: ');
Programmation sous
b=input(' entrer b: ');
Matlab c=input('entrer c : ');
Graphiques 2D et 3D delta=b^2-4*a*c;
if delta > 0
disp('l'‘équation admet deux solutions différentes');
fonction x1=(-b+sqrt(delta))/2*a
Résolution numérique x2=(-b-sqrt(delta))/2*a
elseif delta==0
d’équations
disp('l'‘équation admet une seule solution');
différentielles x2=-b/2*a
Initiation à Simulink else
disp('l'‘équation n''admet pas de solutions');
end
entrer a: 1
entrer b: -3
entrer c : 2
l‘équation admet deux solutions différentes
x1 =
2
x2 =
19
1
M.SAYYOURI
Introduction  On a aussi la possibilité d'utiliser l’instruction switch-case comme pour le langage C :
Programmation sous
Matlab switch expression % variable (comme a) ou % lecture d'un entier

Graphiques 2D et 3D expression de variables (a+b, f(a), …) a = input('Entrer un entier entre 1 et 5:');
case cas1 % valeur 1 de l’expression switch a
case 1
commandes
fonction disp('UN');
case cas2 % valeur 2 de l’expression case 2
commandes disp('DEUX');
d’équations
… case 3
différentielles
otherwise % autre valeurs de l’expression disp('TROIS');
commandes case 4
Initiation à Stateflow end ; disp('QUATRE');
case 5
disp('CINQ');
Exemple otherwise
Écrire un script Matlab qui demande un disp('cet entier n''est pas entre 1 et 5');
nombre compris entre 1 et 5 à l’utilisateur et end;
le convertit en lettre .
Entrer un entier entre 1 et 5 : 3
TROIS
Entrer un entier entre 1 et 5 : 8
20
M.SAYYOURI cet entier n'est pas entre 1 et 5
Introduction Les boucles:
Programmation sous
 En Matlab, on utilise deux formes de boucle:
Matlab
Graphiques 2D et 3D
❖ L’instruction for pour les boucles à nombre d’itération connu :
Recherche de zéro d’une for variable = valeur initiale : pas : valeur finale
fonction commandes
end ;
d’équations
différentielles
❖ et l’instruction while pour les boucles à condition de continuation.
Initiation à Simulink while condition de continuation
Initiation à Stateflow commandes
end ;
% affichage des nombres naturels % avec la boucle while

%impairs et inférieurs à 10 i = 1;
% avec la boucle for while i <= 10
for i = 1:2:10 disp(i);
disp(i); i = i + 2;
end; end; 21
M.SAYYOURI
Introduction
Les fonctions
Programmation sous
 Les fonctions peuvent être déclarées directement dans un fichier script qui doit obligatoirement porter le
Matlab
Graphiques 2D et 3D
même nom de la fonction à déclarer:
Recherche de zéro d’une  Une fonction est déclarée comme suit :
fonction
function nomdefonction(inpar1, inpar 2, …)
% une fonction sans paramètre de sortie.
d’équations
différentielles
function outpar = nomdefonction(inpar1, inpar 2, …)
Initiation à Simulink % une fonction avec un seul paramètre de sortie.
function [outpar1, outpar2, …] = nomdefonction(inpar1, inpar 2, …)
% une fonction avec plusieurs paramètres de sortie.
 Il faut noter que :
➢ La fonction ne se termine pas par end comme dans les autres langages évolués.
➢ Le point-virgule n’est pas nécessaire à la fin des lignes, sauf si un affichage du résultat de la ligne en
question est requis.
22
M.SAYYOURI
Introduction %le factoriel d’un nombre function seconeq(a,b,c); %équation de second degré
function fact(n) delta=b^2-4*a*c;
Programmation sous
s=1; if delta > 0
Matlab disp('l''équation admet deux solutions différentes');
for i=1:n x1=(-b+sqrt(delta))/2
Graphiques 2D et 3D
s=s*i; x2=(-b-sqrt(delta))/2
Recherche de zéro d’une end elseif delta==0
fonction S disp('l''équation admet une seule solution');
>> fact(5) x1=-b/2
s = elseif delta<0
d’équations
120 disp('l''équation n’admet pas de seule solution');
end
différentielles
>> seconeq(1,-3,2)
Initiation à Simulink function [s p]=somprod(n); l'équation admet deux solutions différentes
s=0; p=1; x1 =
for i=1:n 2
s=s+i; x2 =
p=p*i; 1
End
function [s p]=somprod(n);
>>[s p]=somprod(5)
s=0; p=1;
s = i=1;
15 while i<=n
p = s=s+i; p=p*i; i=i+1;
120 end;
23
M.SAYYOURI
Introduction  Opérateurs relationnels et logiques
Programmation sous
Des expressions relationnelles et logiques peuvent être utilisées dans MATLAB exactement comme dans les
Matlab
Graphiques 2D et 3D
autres langages de programmation.
Recherche de zéro d’une ➢ Opérateurs relationnels
fonction Les opérateurs relationnels sont : <, <=, >, >=, ==, ~=
La comparaison d’égalité se fait à l’aide de [==] et l’inégalité à l’aide de [~=]. Ces opérateurs peuvent être
d’équations
différentielles
utilisés avec des scalaires ou des matrices. Le résultat d'évaluation d'une expression relationnelle est 1 (vrai)
Initiation à Simulink ou 0 (faux).
Initiation à Stateflow ➢ Opérateurs logiques
Les expressions relationnelles peuvent être combinées en utilisant les opérateurs logiques suivants : &, | , ~
qui signifient respectivement "et" (AND), "ou" (OR) et "non" (NOT). Ces opérateurs sont appliqués sur les
matrices élément par élément. Les opérateurs logiques ont une priorité plus faible que les opérateurs
relationnels, qui à leur tour ont une priorité plus faible que les opérateurs arithmétiques.
!! Il est conseillé d'utiliser les parenthèses afin d'éviter toute ambiguïté.

M.SAYYOURI 24
Introduction >> a=(1<=5) >> A=[1 2 3 4]
Programmation sous a= A=
1 2 3 4
Matlab 1
>> B=[4 2 5 4]
Graphiques 2D et 3D >> b=(5==12/3) B=
Recherche de zéro d’une b= 4 2 5 4
fonction 0 >> A==B
Résolution numérique >> c=(a | b) ans =
0 1 0 1
d’équations c=
>> A<B
différentielles 1 ans =
Initiation à Simulink >> d= (a & b) | c 1 0 1 0
Initiation à Stateflow d= >> A~=B
1 ans =
1 0 1 0
>> f=(~a)&(~b)
>> ~(A<B)
f= ans =
0 0 1 0 1
x=(a~=b) >> A~=B
x= ans =
0 1 1 0
1
M.SAYYOURI 25
?
Calcul polynômial
Introduction
Les polynômes
Calcul polynômial
 Définition
Graphiques 2D et 3D • Un polynôme de la forme p ( x)  a x  a x  ...  a x  a x  a est représenté en Matlab par un
n
n
n 1
n 1
2
2
1 0
Recherche de zéro d’une vecteur ligne [a , a ,..., a , a , a ] (ou vecteur colonne) .

n n 1 2 1 0
fonction • Afin d’évaluer un polynôme pour une valeur fixe a, on utilise la commande polyval(P, a).
Résolution numérique  Opérations arithmétiques
d’équations
• Matlab fournit deux opérations arithmétiques de multiplication et de division à travers les commande conv
différentielles
et deconv. Il reste au programmeur de programmer les deux autres opérations d’addition et de soustraction.
• C = conv(A, B) est la convolution des tableaux A et B, c'est à dire les coefficients du produit des deux
polynômes.
• [Q, R] = deconv(A, B) est la déconvolution des tableaux A et B ou Q est le quotient de la division et R est le
reste (A = conv(B ,Q) + R).
 Exemple
Soient p et q deux polynômes définis par :
P(x)=2x^3-x^2+3x-5 et q(x)=3x^2-5x+3
Calculer p(2) et q(1.5)
Calculer s(x)=p(x)+q(x) ; r(x)=2p(x)-3q(x) ; t(x)=p(x)q(x) ; u(x)=p(x)/q(x)
M.SAYYOURI 26
?
Calcul polynômial
Introduction
Les polynômes
Calcul polynômial
>> p=[2 -1 3 -5], q=[3 -5 3] >> t=conv(p,q)
Graphiques 2D et 3D
p= t=
2 -1 3 -5 6 -13 20 -33 34 -15
fonction
q= >> [u,v]=deconv(p,q)
Résolution numérique 3 -5 3
d’équations >> a=polyval(p,-2) u=
différentielles a= 2/3 7/9
Initiation à Simulink -31
>> b=polyval(q,1.5) v=
b= 0 0 44/9 -22/3
9/4
>> s=p+[0 q] >> conv(q,u)+v
s= ans =
2 2 -2 -2 2 -1 3 -5
>> r=2*p-3*[0 q] >> p
r= p=
4 -11 21 -19 2 -1 3 -5
M.SAYYOURI 27
?
Calcul polynômial
Introduction  Racines et Interpolations
Calcul polynômial • La commande roots(P) fait extraire les racines d’un polynôme P.
Graphiques 2D et 3D • La commande poly(V) retourne un polynôme depuis ses racines stockées dans le vecteur V.
• La commande Polyfit(X, Y, n) permet d’approximer un polynôme de dégrée n qui passent
fonction
approximativement par les points (X, Y). Cette approximation P(X(i))=Y(i) est au sens des moindres carrées :
d’équations p1=[1 -5 6] >> x=[0 1 4]; y=[0 2 3];

différentielles p1 = >> p11=polyfit(x,y,2)
Initiation à Simulink 1 -5 6 p11 =
Initiation à Stateflow >> a=roots(p1) -0.4167 2.4167 -0.0000
a= >> p12=polyfit(x,y,1)
3 p12 =
2 0.6538 0.5769
>> poly(a) >> polyval(p11,4)
ans = ans =
1 -5 6 3
>> poly([1 1 1]) >> polyval(p12,4)
ans = ans =
1 -3 3 -1 3.1923
M.SAYYOURI 28
?
Calcul polynômial
Introduction  Dérivation et Intégration d’un polynôme
Calcul polynômial • La commande polyder(P) retourne la dérivée d’un polynôme P.
Graphiques 2D et 3D • La commande polyint(P) retourne l’intégral d’un polynôme P. Si on veut calculer l’intégral entre deux
point x1 et x2, on utilise la commande polyval pour évaluer l’intégral sur ces deux points puis on applique
fonction
la soustraction.
d’équations >> p=[1 1 1 1]; % p(x)=x^3+x^2+x+1

différentielles >> der_p=polyder(p)
Initiation à Simulink der_p =
Initiation à Stateflow 3 2 1 %p’(x)=3x^2+2x+1
>> der_p=polyder(polyder(p))
der_p =
6 2 %p’’(x)=6x+1
>> int_p=polyint(p)
int_p =
1/4 1/3 1/2 1 0 %intégral de p(x)=1/4x^4+1/3x^3+1/2x^2+x+0
>> int_p12=polyval(int_p,2)-polyval(int_p,1)
int_p12 =
103/12
M.SAYYOURI 29
Résolution des Systèmes Linéaires
?
Introduction
 Résolution des Systèmes Linéaires
Résolution des Systèmes
• Un système linéaire est un ensemble de n équations à p variables de la forme :
Linéaires
Graphiques 2D et 3D  a x  a x  ...  a x
11 1 12 2 1p p
b1

 a x  a x  ...  a x b
21 1 22 2 2p p 2

fonction
...
d’équations
 a x  a x  ...  a x b
différentielles
 p 11 1 p 12 2 p 1 p p p

 a x  a x  ...  a x
p1 1 p2 2 pp p
b p
Les sont aij appelés les coefficients du système ; les bk sont le second membre du système (ou la partie
gauche).
En calcul matriciel, un système linéaire s’écrit sous la forme: AX = B, ou
• A matrice des coefficients.
• X vecteur colonne des inconnues.
• B vecteur colonne comme second membre. 30
M.SAYYOURI
?
Introduction  Résolution d’un système linéaire

Résolution des Systèmes ➢ Matlab fournie plusieurs techniques afin de résoudre un système linéaire. Néanmoins, il reste a
Linéaires
l’utilisateur de vérifier la compatibilité du système linéaire (s’il existe solution).
Graphiques 2D et 3D
➢ Par exemple, si le déterminant de A est non-nul alors le système accepte une solution. Par contre, si le
Recherche de zéro d’une déterminant est nul le système peut accepter un nombre infini de solution ou un ensemble vide des
fonction solutions.
Résolution numérique ➢ La commande X = A\B fait résoudre le système linéaire AX = B (noter l’antislash).
d’équations ➢ La commande X = B/A fait résoudre le système linéaire XA = B (noter le slash).
différentielles On peut calculer directement l’inverse de A en le multipliant par B. Il faut noter que la multiplication
matricielle n’est pas commutative !
➢ La commande X = inv(A)*B fait résoudre le système linéaire AX = B.
➢ La commande X = B*inv(A) fait résoudre le système linéaire XA = B.
Exemple
Résoudre le système suivants:
x  2 y  z  1  1 2 1 x   1 
  1 1 1 y    0 
x  y  z  0     
2 x  3 y  z  1  2 3 1 z  1 
M.SAYYOURI       31
?
Introduction Pour les commande A\B et B/A, si la matrice A n’est pas singulière, le
Résolution des Systèmes
>> A=[1 2 1;1 -1 1;2 3 1],B=[1 ;0; 1] résultat est (0/0) pour un nombre infini de solutions, et (k/0, k≠0) pour
A= l’ensemble vide (). Pour les commande inv(A)*B et B*inv(A), si la
Linéaires
1 2 1 matrice A n’est pas singulière, le résultat est toujours (k/0, k≠0).
Graphiques 2D et 3D 1 -1 1
Recherche de zéro d’une 2 3 1 >> A1=[1 2;2 4],B1=[2 ;3] >> A1=[1 2;2 4],B1=[2 ;4]
fonction B= A1 = A1 =
Résolution numérique 1 1 2 1 2
d’équations
0 2 4 2 4
différentielles
1 B1 = B1 =
>> X=A\B 2 2
X= 3 4
-1/3 >> X1=A1\B1 >> X1=inv(A1)*B1
1/3 X1 = X1 =
2/3 1/0 1/0
>> X=inv(A)*B -1/0 1/0
X= >> X1=inv(A1)*B1 >> X1=A1\B1
-1/3 X1 = X1 =
1/3 1/0 0/0
2/3 1/0 0/0
M.SAYYOURI 32
?
Calcul symbolique Matlab

Introduction
 Résolutions des équations et des inéquations, linéaires et non-linéaires
Calcul symbolique
Matlab utilise la commande solve afin de résoudre les équations et les inéquations, les systèmes linéaires et
Matlab
non linéaires.
Graphiques 2D et 3D

➢ Tout d’abord, il est nécessaire de déclarer des variables symboliques des équations (ou inéquations, …)
fonction
en utilisant la commande syms.
Résolution numérique ➢ Il faut noter que solve retourne des solutions sous forme symbolique ou formelle au lieu de numérique
d’équations (racine de 2 a la forme numérique de 1.4142 et la forme formelle 2^(1/2)).
différentielles >> sqrt(sym(2)) >> syms x %équation de 2 >> syms a b x inéquation de 2 degré
Initiation à Simulink ans = degré >> solve(a*x+b,x)
Initiation à Stateflow 2^(1/2) >> solve(x^2-3*x+2) ans =
>> a=sym(1)/sym(3) ans =
a= 1 -b/a
1/3 2 >> syms a b c x
>> b=sym(2)/sym(5) >> syms x, solve(x^2+x+1) >> solve(a*x^2+b*x+c, x)
b =2/5 ans = ans =
>> c=a+b - 1/2 - (3^(1/2)*i)/2 -(b + (b^2 - 4*a*c)^(1/2))/(2*a)
c =11/15 - 1/2 + (3^(1/2)*i)/2 -(b - (b^2 - 4*a*c)^(1/2))/(2*a)
M.SAYYOURI 33
?
Calcul symbolique sous Matlab

Introduction La commande solve permet aussi de résoudre une inéquation ou un système linéaire ou non linéaire des
Calcul symbolique inéquations.
Matlab
>> syms x y % système à deux inconnues >> syms x % inéquation de 2 degré
Graphiques 2D et 3D
>> S = solve(x+y-5, x-y-1) >>solve('x^2-9 <=0')
S= ans =
fonction x: [1x1 sym]
Résolution numérique y: [1x1 sym] Dom::Interval([-3], [3])
d’équations >> [S.x S.y] >> syms x, solve('x^2-9 <0')
ans = ans =
différentielles
[ 3, 2]
>> syms x y z % système à trois inconnues Dom::Interval(-3, 3)
>>S=solve(x+y+z-5, x-y-z-1,2*x+y-2*z+1) >> syms x, solve('x^2-9 >=0')
S= ans =
x: [1x1 sym]
y: [1x1 sym] Dom::Interval([3], Inf)
z: [1x1 sym] Dom::Interval(-Inf, [-3])
>> [S.x S.y S.z]
ans =
[ 3, -1, -1]
M.SAYYOURI 34
?

Introduction La boîte à outils de Calcul symbolique permet de calculer la dérivé , l’intégral et les limites d’une fonction.
Calcul symbolique
 Calcul de la dérivée d’une fonction
Matlab
diff(f) calcule la dérivée de f par rapport à x :
Graphiques 2D et 3D

Pour calculer la dérivée seconde : diff(f,2) ou bien diff(f,x,2)
fonction La fonction diff peut également agir sur une matrice. Dans ce cas la dérivation est effectuée sur chaque
Résolution numérique élément de la matrice,
d’équations
>> syms a b x >> A=[cos(3*x) sin(2*x); tan(x) 1]
différentielles
>> f=a*cos(x)+b A=
[ cos(3*x), sin(2*x)]
f=
Initiation à Stateflow [ tan(x), 1]
b + a*cos(x) >> diff(A)
>> diff(f) ans =
ans = [ -3*sin(3*x), 2*cos(2*x)]
-a*sin(x) [ tan(x)^2 + 1, 0]
>> diff(A,2)
>> diff(x^3+2*x^2+3*x+1,x)
ans =
ans = [ -9*cos(3*x), -4*sin(2*x)]
3*x^2 + 4*x + 3 [ 2*tan(x)*(tan(x)^2 + 1), 0]
M.SAYYOURI 35
?

Introduction  Calcul de l’intégral d’une fonction
Calcul symbolique
Si f est une expression symbolique, alors int(f)tente de trouver une autre expression symbolique F telle que
Matlab
diff(F)=f , c’est-à-dire que int(f) donne l’intégrale indéfinie de f si elle existe.
Graphiques 2D et 3D
Recherche de zéro d’une >> syms x >> syms a b w t

fonction >> int(x^3+x+1) >> f=a*cos(w*t+b)
Résolution numérique ans = f=
d’équations
(x*(x^3 + 2*x + 4))/4 a*cos(b + t*w)
>> int(x^3+x+1,0,1) >> int(f,t)
différentielles
ans = ans =
7/4 (a*sin(b + t*w))/w
>> int(log(x)) >> int(f,t,0,pi)
ans = ans =
x*(log(x) – 1 (a*(sin(b + pi*w) - sin(b)))/w
 Exercice
La vitesse d’une particule est exprimée par la fonction v(t)=tsin(t/2)
1. Trouver la position de la particule 𝑥(𝑡) sachant que la position initiale est 𝑥(0) = 0
2. Trouver l’accélération de la particule 𝑎(𝑡) 36
M.SAYYOURI
?

Introduction  Calcul des limites d’une fonction
Calcul symbolique Les commandes syms a x, limit( f(x),x,a) permettent le calcul de la limite de f en a
Matlab
>> syms x >> syms n x >> syms x >> syms x
Graphiques 2D et 3D
>> limit(sin(x)/x,x,0) >> limit((1+x/n)^n,n,inf) >> limit(log(x)/x,x,inf) >> limit(x*log(x),x,0)
ans ans = ans = ans =
fonction 1 exp(x) 0 0
d’équations  Calcul des sommes

différentielles
On peut calculer des sommes symboliques quand elles existent en utilisant la commande symsum
Exemple
1 

2
• Ecrire la commande qui permet de calculer  

k
k 0
6 2
1 
x 
k
• Ecrire la commande qui permet de calculer pour |x| < 1
1 x
k 0
>> syms k n >> syms k x

>> symsum(1/k^2,1,inf) >> symsum(x^k,k,0,inf)
ans = ans =
pi^2/6 piecewise([1 <= x, Inf], [abs(x) < 1, -1/(x - 1)])
M.SAYYOURI 37
Graphique 2D
Introduction
 Les possibilités graphiques de Matlab sont innombrables. Pour créer le graphe de la fonction sinus par
exemple, nous commençons par définir deux vecteurs: x = -2*pi : pi/4 : 2*Pi ; y = sin(x);
Graphiques 2D et 3D
 Ensuite nous construisons le graphe au moyen de la commande: plot(x , y)
fonction
 On peut superposer un quadrillage (ou le faire disparaitre) par les commandes: grid (ou grid off)
Résolution numérique  On peut aussi mettre un titre, mettre des étiquettes aux axes: title ('graphe de sinus'), xlabel('x'),
d’équations ylabel('sinus(x)')
différentielles
Initiation à Stateflow >> x = -2*pi : 0.01 : 2*pi ;

>> y = sin(x);
>> plot(x, y)
>> grid % afficher la grille
>> % afficher le titre et les étiquettes
>> title ('Graphe de Sinus')
>> xlabel('x')
>> ylabel('y=sinus(x)')
38
M.SAYYOURI
Graphique 2D
Introduction  Afin d'afficher deux ou plusieurs graphes en même temps, une série des coordonnées (x, y) est
Calculs élémentaires introduise a la fonction plot: plot(x1, y1, x2, y2, …) .
Graphiques 2D et 3D  La legende (surtout en cas de plusieurs graphes) peut être affichée en utilisant la commande
Recherche de zéro d’une legend(legend1, legend2, …).
fonction  La commande legend('show') permet de montrer la légende, par contre la commande
Résolution numérique legend('hide') permet de la cacher.
d’équations
différentielles
>> x = -2*pi : 0.01 : 2*pi ; y = sin(x); z=cos(x);
>> plot(x, y, x, z)
>> legend('sin(x)','cos(x)')
>> legend('hide')
>> legend('show')
>> title ('Graphe de Sinus et cosnus')
>> xlabel('x')
>> ylabel('y')
39
M.SAYYOURI
Graphique 2D
Introduction
 La commande plot(x, y, LineSpec) affiche un graphe avec la spécification donnée par LineSpec qui est
formée au plus de trois valeurs: style de ligne, type de marquer, et la couleur.
Graphiques 2D et 3D Par exemple plot(x, y, '-+b') affiche un graphe bleu de la forme -+-+-+…
Recherche de zéro d’une  Noter bien que la spécification 'type de marqueur' est utilisée pour les points choisies (vecteurs x et y), par
fonction contre la spécification 'style de ligne' couvre toute la courbe.
Symbole Style de ligne Symbole Type de Marqueur
d’équations
- Ligne solide (par défaut) o Cercle
différentielles
-- La ligne non continue + Signe plus
: La ligne pointillée * Astérisque
Initiation à Stateflow -. La ligne tiret-point . Point
Symbole Couleur x Croix
y Jaune s Carré
m Magenta d Diamant
c Cyan ^ Triangle vers le haut
r Rouge V Triangle vers le bas
g Vert > Triangle pointant à droite
b Bleu < Triangle pointant à gauche
w Blanc p Pentagramme
k Noir h Hexagramme 40
M.SAYYOURI
Graphique 2D
Introduction
Exemple
Graphiques 2D et 3D >> x = -2*pi : 0.1 : 2*pi ;

>> y = sin(x);
fonction
>> z=cos(x);
d’équations
>> legend('sin(x)','cos(x)')
différentielles >> plot(x, y, '-ob', x, z, '-.dr')
Initiation à Simulink >> legend('sin(x)','cos(x)')
Initiation à Stateflow >> title('Graphe de Sinus et cosnus')
>> xlabel('x')
>> ylabel('y')
41
M.SAYYOURI
?
Graphique 2D
Introduction  Exemples 𝑓 𝑥 = sin 𝑥 𝑠𝑖 𝑥 > 0
Calculs élémentaires Imaginons que l'on veut définir la fonction suivante :ቊ
𝑓 𝑥 = sin 2𝑥 𝑠𝑖 𝑥 < 0
Graphiques 2D et 3D
C'est une fonction définie par morceaux puisque sa valeur dépend d'un test logique définissant, dans cet
exemple, deux domaines de définition.
fonction
Voilà comment définir y en fonction de x:
>> x=-4*pi : pi/20 : 4*pi;
d’équations
>> y=sin(x).*(x>0)+sin(2*x).*not(x>0);
différentielles
>> plot(x,y,'-b*')
Initiation à Simulink >> xlabel('x')
Initiation à Stateflow >> ylabel('y')
>> title ('sin(x) et sin(2x)')
x=-4*pi:pi/20:4*pi;
inds1=find(x>0);
x1=x(inds1);
inds2=find(x<=0);
x2=x(inds2);
y1=sin(x1);
y2=sin(2*x2);
y=[y2, y1];
plot(x,y,'-b*') 42
M.SAYYOURI
Graphique 3D
Introduction
 Courbes dans un espace 3D
Comme une courbe en 2D est définie par une liste de doublets (x,y) on peut définir une courbe en 3D par
Graphiques 2D et 3D
une liste de triplets (x,y,z). Puisque l'instruction plot attendait deux arguments, le vecteurs des abscisses, et
Recherche de zéro d’une celui des ordonnées, l'instruction plot3 en attend trois : le vecteur des abscisses Vx, le vecteur des
fonction ordonnées Vy et le vecteur des cotes Vz.
Résolution numérique Syntaxe: plot3(Vx, Vy, Vz) ;
d’équations Voici un exemple de courbe paramétrée tracée dans un espace à trois dimensions :
différentielles
Initiation à Simulink  x=sin(t)e-t/10


 y=cos(t)e
-t/10
t  [-2π,2π]
 z=e-t

t = -2*pi:2*pi/100:2*pi;
plot3(exp(-t/10).*sin(t), exp(-t/10).*cos(t), exp(-t));
grid;
43
M.SAYYOURI
Graphique 3D
Introduction
 Tracé d'un nuage de points 3D
Le principe est le même que pour un nuage 2D: les points possèdent 3 coordonnées (x , y , z) représentés
Graphiques 2D et 3D
par des marqueurs non liés.
fonction T1=ones(1,100);
Résolution numérique % Premier nuage de valeurs
d’équations
A = 0.2 * randn(3,100) + [1*T1 ; 2*T1 ; 1*T1]
différentielles
plot3( A(1, : ) , A(2 , : ) , A(3 , : ) , '.' )
hold on
% Deuxième nuage de valeurs
B = 0.6 * randn(3,100) + [2*T1 ; 3*T1 ; 2*T1] ;
plot3( B(1, : ) , B(2 , : ) , B(3 , : ) ,'+r')
grid on ; axis([0 3 0 4 0 4])
44
M.SAYYOURI
Graphique 3D
Introduction  Tracé d'une surface
Calculs élémentaires Une surface est une application de ( x , y) vers z de la forme z= f ( x , y) .
Graphiques 2D et 3D Exemple: le plan est décrit par l'équation z=a x+b y+z0
Recherche de zéro d’une La représentation est calculée pour un maillage des
fonction valeurs de l'espace de départ; le tracé de ce maillage
Résolution numérique forme une grille rectangulaire ou carrée. Chaque point du
d’équations
maillage est caractérisé par le couple des coordonnées
(x , y) .
différentielles
Initiation à Simulink Exemple de plan

a = 1.5 ; b = .7 ; z0 = 2 ;
x = [ -5 : .5 : 5 ]; y = [ -7 : .3 : 7 ];
[X Y]=meshgrid(x,y)
Z=a*X+b*Y+z0 ;
surf(Z)
Les fonctions les plus courantes sont :

mesh, qui trace une série de lignes entre les points de la surface en mode «lignes cachées» ;
meshc, qui fonctionne comme mesh mais en ajoutant les courbes de niveau dans le plan ;
surf, qui «peint» la surface avec une couleur variant selon la cote ;
surfl, qui «peint» la surface comme si elle était éclairée ;
surfc, qui fonctionne comme mesh mais en ajoutant les courbes de niveau dans le plan 45
M.SAYYOURI
Graphique 3D
Introduction Exemple : Tracer la représentation graphique de la fonction f(x,y)=x^2+y^2 sur [-1 1] x [-2 2]
Graphiques 2D et 3D x = -1:0.2:1;
Recherche de zéro d’une y = -2:0.2:2;
fonction
[X,Y] = meshgrid(x,y);
Z = Y .^2 + X.^2;
subplot(221)
d’équations
mesh(X, Y, Z);
différentielles
xlabel('x'); ylabel('y'); zlabel('z'); title('mesh’);
subplot(222)
Initiation à Stateflow meshc(X, Y, Z);
xlabel('x'); ylabel('y');zlabel('z'); title('meshc');
subplot(223)
surf(X, Y, Z);
xlabel('x'); ylabel('y');zlabel('z'); title('surf');
subplot(224)
surfl(X, Y, Z);
xlabel('x'); ylabel('y’); zlabel('z'); title('surfl');
46
M.SAYYOURI
Graphique 3D
Introduction  Tracé de courbes de contour
Calculs élémentaires Pour tracer les courbes z=cste d'une surface définie par z=f(x,y) , on peut utiliser la fonction contour.
Graphiques 2D et 3D Elle s'utilise comme les instructions précédentes, mais fournit un graphe 2D dans le plan . Plusieurs paramètres
Recherche de zéro d’une optionnels peuvent être spécifiés, notamment le nombre de courbes de contours à afficher.
fonction x = -1:0.05:1; x = -1:0.05:1; y = -2:0.05:2;
Résolution numérique y = -2:0.05:2; [X,Y] = meshgrid(x,y);
[X,Y] = meshgrid(x,y); Z = Y .^2 + X.^2;
d’équations
Z = Y .^2 + X.^2; ch= contour(X,Y,Z,5);
différentielles
contour(X,Y,Z,5); clabel(ch);
il est également
intéressant d'ajouter
sur le tracé les valeurs
de le long des
différentes lignes de
contour. Pour cela il
faut affecter à une
variable la valeur de
sortie de la commande
contour et utiliser la
fonction clabel
47
M.SAYYOURI

Cours Matlab - Partie 1

Transféré par

Droits d'auteur :

Formats disponibles

Cours Matlab - Partie 1

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Cours Matlab - Partie 1

Transféré par

Droits d'auteur :

Formats disponibles

UNIVERSITÉ CHOUAIB DOUKKALI

Ecole Nationale des Sciences Appliquées d’El Jadida-ENSAJ

Année universitaire 2017-2018

Recherche de zéro d’une

Recherche de zéro d’une

Recherche de zéro d’une

Recherche de zéro d’une

Initiation à Simulink La fenêtre de répertoire courant :

Recherche de zéro d’une

11 , e , 9 , ln(5) , log (5) , sin( ) , 1   , 1  2i  (1  3i )

différentielles 2. Saisir la matrice

Introduction des données à partir du clavier :

Entrer la valeur de a : -2.5

Matlab switch expression % variable (comme a) ou % lecture d'un entier

% affichage des nombres naturels % avec la boucle while

!! Il est conseillé d'utiliser les parenthèses afin d'éviter toute ambiguïté.

Recherche de zéro d’une vecteur ligne [a , a ,..., a , a , a ] (ou vecteur colonne) .

d’équations p1=[1 -5 6] >> x=[0 1 4]; y=[0 2 3];

d’équations >> p=[1 1 1 1]; % p(x)=x^3+x^2+x+1

Introduction  Résolution d’un système linéaire

Calcul symbolique Matlab

Recherche de zéro d’une

Calcul symbolique sous Matlab

Calcul symbolique sous Matlab

Recherche de zéro d’une

Calcul symbolique sous Matlab

Recherche de zéro d’une >> syms x >> syms a b w t

Calcul symbolique sous Matlab

d’équations  Calcul des sommes

• Ecrire la commande qui permet de calculer  

>> syms k n >> syms k x

Initiation à Stateflow >> x = -2*pi : 0.01 : 2*pi ;

Graphiques 2D et 3D >> x = -2*pi : 0.1 : 2*pi ;

Initiation à Simulink  x=sin(t)e-t/10

Initiation à Simulink Exemple de plan

Les fonctions les plus courantes sont :

Vous aimerez peut-être aussi

Initiation à Stateflow >> x = -2pi : 0.01 : 2pi ;

Graphiques 2D et 3D >> x = -2pi : 0.1 : 2pi ;