Université frères Mentouri Constantine 1
Faculté des Sciences Exactes
Département de Physique
Laboratoire de Physique Energétique
II.1 CALCUL POLYNOMIAL
Représentation d'un polynôme en MATLAB
Un polynôme de degré n s'écrit en mathématique sous la forme générale:
−1
( )= + −1 +⋯⋯⋯+ 1 + 0
En Matlab, P(x) est représenté par le vecteur des coefficients:
=[ −1 −2 ⋯⋯⋯ 1 0 ]
(du plus haut degré vers le plus bas degré).
Les fonctions MATLAB utilisées pour le calcul polynômial
Constitution du polynôme p(x) à partir de ses racines
P = poly (V)
V: vecteur dont les composantes sont les racines de P(x) = 0.
Evaluation des valeurs d'un polynôme p(x)
y= polyval (p,w)
w: scalaire ou vecteur
Racines d'une équation polynômiale, p(x)=0
V=roots(p)
V: vecteur contenant les racines de l'équation
Dérivée du polynôme p(x)
pr = polyder(p)
Intégrale d'un polynôme avec une constante d'intégration k
g = polyint(p,k)
Le polynôme de degré m qui ajuste n points (xi,yi): (Cours IV. Interpolation et ajustement)
p = polyfit(x,y,m)
Décomposition p(x)/q(x), en éléments simples
[ , , ]= ( , )
D'où la décomposition suivante:
( ) 1 2
= + +⋯⋯⋯⋯⋯+[ ]
( ) − 1 − 2
1
Produit de deux polynômes (Convolution)
h=conv(p,q)
Rapport de deux polynômes (Dé-convolution)
[d , r]=deconv (p , q)
d: le quotient et r: le reste de la division
Exemple
Résoudre les équations suivantes:
2
− −2=0
7 4 3
−3 + − −1=0
4 2
− +3=0
Exemple
Déterminer le polynôme P(x) d'ordre 4, qui permet d'ajuster les données suivantes:
x 1.15 1.75 2.35 4.00 4.55 5.25 6.00
y 3.25 5.35 7.25 10.00 11.75 10.50 9.00
Calculer les valeurs de P, en 10 points régulièrement espacés entre 0 et 6
II. VISUALISATION GRAPHIQUE
Matlab est équipé d'outils de visualisation et d'analyse graphique très puissants.
III.1 COURBES EN 2 DIMENSIONS
Courbe à partir de deux vecteurs de données x, y
plot(x,y)
Exemple: Courbe de la vitesse en fonction du temps.
temps = [ 0 1 2 3 4 5 6 7] ;
vitesse = [ 0 0.75 2.8 3.32 4.7 3.1 1.1 0.3];
plot (temps , vitesse)
Courbe d'une fonction ( )
fplot ( , [ ])
La fonction est définie par la syntaxe suivante:
= @( , , , … … ) expression de
Exemple: Tracer la courbe de la fonction:
2
− 2
( )=
2
√ +1
avec: − ≤ ≤
Options du graphisme
Plusieurs fonctions peuvent-être appliquées sur un graphe:
xlabel , ylabel, title, grid on, grid off, axis, gtext, legend, hold on, hold off ......
On peut également tracer une courbe sous forme de symboles (o , * , + , - , .... ) et utiliser une
couleur (r, g, b, y, w, k, m).
3
Exemple
Déterminer le polynôme de degré 3 qui permet d'ajuster les données suivantes:
x 0.00 0.25 0.61 0.87 1.00
y 0.00 0.85 1.36 3.85 5.62
Tracer les données et la courbe du polynôme d'ajustement sur le même graphe.
Programme
clc
clear all
x= [0 0.25 0.61 0.87 1];
y= [0 0.85 1.36 3.85 5.62];
p= polyfit(x,y,3)
yy=polyval(p,x);
plot(x,y,'ko',x,yy,'r')
grid on
xlabel('x')
ylabel('y')
title('Ajustement des données par un polynôme')
legend('Données','Polynôme d''ajustement')
III.2 COURBES EN 3 DIMENSIONS
La commande utilisée est plot3
Exemple(1)
0≤ ≤ 10 avec ∆ = , x=cos(t) , y=sin(t) , z=t ;
10
Programme
t = 0:pi/10:10*pi;
x=cos(t);y=sin(t);z=t;
plot3(x,y,z)
Exemple(2)
La sphère de Rayon r = 1.
Programme
clc
clear all
r=1;n=30;
phi=0:pi/10:2*pi*n;
teta=phi/(2*n);
x=r*sin(teta).*cos(phi);
y=r*sin(teta).*sin(phi);
z=r*cos(teta);
plot3(x,y,z,'r')
axis equal
4
III.3 COURBE D'UNE FONCTION A DEUX VARIABLES: = ( , )
Démarches à suivre:
Définir les vecteurs de données x , y
Etablir le maillage du domaine de calcul, par la construction des matrices X et Y,
contenants les coordonnées des nœuds:
[x,y]=meshgrid(x,y)
Calculer la valeur de z en chaque noeud.
Pour tracer la courbe = ( , ), utiliser l'une des fonctions graphique:
mesh , meshc , meshz, surface et contour
Exemple
−2 ≤ ≤2 , −2 ≤ ≤ 2 , ∆ = ∆ = 0.1
− 2− 2
Tracer la courbe: =
Programme
clc
clear all
x=-2:0.1:2;y=-2:0.1:2;
[X,Y]=meshgrid(x,y);
z=Y.*exp(-X.^2-Y.^2)
mesh(X,Y,z)
5
III.4 ORGANISATION DES FIGURES
Numérotation des figures : figure (n°)
Fermeture d'une figure: close (n°)
Partage d'une figure en zones graphiques: subplot(m,n,l)
m: nombre de lignes
n: nombre de colonnes
l: numéro de la partie de la figure
