TP 04

I.N.S.S.E.T. – U. P. J.V.
Programmation avec Matlab

Licence Maths et SPI – 2ème année M. Marcoux
TP 4 : Polynômes et courbes
1. Introduction
Le but de ce TP est de voir les fonctions permettant de manipuler les polynômes, de tracer des
courbes et de réaliser des approximations de données.
Vous travaillerez dans le répertoire "c:/TPMatlab/votrenom/TP4".
2. Polynômes
a) Généralités
Trouver les racines d’un polynôme P (x) consiste à chercher les valeurs de x qui annulent ce
polynôme
Matlab représente un polynôme comme un vecteur ligne contenant les coefficients du
polynôme rangés par ordre décroissant d’exposant.
Par exemple, en lieu et place de P( x) = x 4 − 12 x 3 + 25 x + 116 ,
on écrit le vecteur P = [1 − 12 0 25 116]
On obtient alors instantanément les racines d’un polynôme avec la commande roots.
Regardez les possibilités de cette commande avec la fenêtre d’aide.
Dans Matlab, un polynôme et ses racines sont des vecteurs Le polynôme étant un vecteur
uniligne, et les racines étant un vecteur unicolonne.
b) Application
Valeur en un point
Dans un script intitulé "poly1.m" :
Ø Créez le polynôme P (x) précédent
• A l’aide de la fenêtre d’aide, déterminez les fonctionnalités de la commande polyval

Ø Déterminez la valeur de P (x) en x = 1 et affichez-la
Ø Créez un vecteur V dont les composantes sont les valeurs de 0 à 20 par pas de 0.5
Ø Calculez les valeurs de P (x) aux points contenus dans le vecteur V et affichez les
Ø Affichez ces valeurs dans un graphique pour connaître l’allure de ce polynôme
Recherche de racines
• Utilisez la fenêtre d’aide pour analyser la commande roots

Soient les deux polynômes suivants :
A( x) = x 3 + 2 x 2 + 3 x + 4 et B( x) = x 3 + 4 x 2 + 9 x + 16
Dans un script intitulé "poly2.m" :

Ø Créez les 2 polynômes A(x) et B (x) précédents
Ø Tracez leurs courbes dans une figure simple, sur l’intervalle [−10 ,10]
Ø Trouvez leurs racines et affichez les
• Vérifiez ces valeurs à l’aide de polyval
Addition
Ø Ajoutez au script précédent le calcul de la somme S (x) de A(x) et de B (x)
Ø Tracez le résultat obtenu
Multiplication
• Utilisez la fenêtre d’aide pour analyser la commande conv

Ø Ajoutez au script précédent le produit C (x) de A(x) et de B (x)
• Vérifiez à l’aide du calcul théorique
Division
• Servez vous de la fenêtre d’aide pour découvrir la commande deconv

Ø Ajoutez au script précédent le calcul de la division de C (x) par B (x)
3. Tracé de courbes
a) Introduction
La visualisation de données sous forme de courbes est un point important.
Lorsqu’une fonction est étudiée, comme cela est le cas en analyse numérique, il n’est pas très
intéressant de lister l’ensemble des valeurs prises par la fonction. Il est plus pratique et plus
lisible de tracer la courbe, comme nous venons de le voir.
Sous Matlab, l’affichage d’un graphique nécessite d’ouvrir une fenêtre graphique, ceci est
réalisé par la commande figure, que nous avons déjà rencontrée au TP n°1.
Il est possible de redimensionner et de déplacer cette fenêtre de façon classique.
A chaque fenêtre est affecté un numéro identificateur (1 pour le première, 2 pour la seconde,
etc.). La fermeture des fenêtres se fait alors avec la commande close.
Exemple :
• Ouvrir trois fenêtres graphiques (commande figure)

• Les déplacer et les redimensionner à loisir
• Fermer la seconde en utilisant la commande close
• Fermer tous les fenêtres graphiques restantes
b) Graphique
Comme nous l’avons déjà vu avec le TP n°1, la commande plot permet de tracer un vecteur
sous forme graphique.
La syntaxe la plus simple, que vous avez utilisé dans la première partie du TP, est la
suivante : plot(V), où V est un vecteur ligne ou colonne, peu importe).
Le graphique s’affiche alors sous forme de lignes brisées reliant les points définis par une
abscisse et une ordonnée. Pour chaque point, l’abscisse est l’indice de l’élément du vecteur V
et l’ordonnée est la valeur de l’élément.
Une fois l’affichage effectué, des commentaires peuvent s’ajouter au graphique : un titre
(title), la légende de l'axe des abscisses (xlabel) et de l’axe des ordonnées (ylabel).
Le graphique complet peut alors être éventuellement sauvegardé en allant dans le menu
File/save_as.
Une grille peut aussi être ajoutée à l’aide de la commande grid.
Application
Ø Créez un script "graphesin.m" réalisant les actions suivantes :
o Demande à l’utilisateur de trois valeurs : min, max et pas
o Génération d’un vecteur X contenant les points compris entre min et max avec le
pas spécifié
o Calcul du vecteur Y associé à la fonction f ( x) = sin( x)
o Tracé du graphique correspondant à la fonction f (x)
o Ajout d’un titre et des légendes des axes des x et des y par demandes à l’utilisateur
o Ajout ou pas d’une grille par questionnement
o Sauvegarde sur demande du graphique dans le répertoire sous un nom à demander
o Fermeture de toutes les fenêtres graphiques
Après avoir créez plusieurs graphiques, testez la commande open
c) Application aux polynômes
Ø Créez un script "graphepoly.m" réalisant les actions suivantes :

o Création d’un vecteur X contenant 100 points compris entre -1 et 3
o Calcul des valeurs du polynôme P (x) sur l’intervalle [− 1, 3] dans un vecteur V
o Affichage de la courbe ainsi obtenue.
o Ajout du titre et des légendes aux axes
o Création d’un fenêtre graphique pouvant contenir 2 sous fenêtres
o Tracé de la courbe de P (x) dans la première sous fenêtre.
o Calcul de la dérivée dP du polynôme P (x) (commande polyder )
o Tracé dans la seconde sous fenêtre graphique de la courbe dP sur [− 1, 3]
Ø Déterminez les extremums du polynôme P (x) sur l’intervalle [− 1, 3] .

Ø Affichez leurs valeurs
4. Approximation de données
a) Fonction polyfit
Une fonction très utile dans les laboratoires est la fonction polyfit(x,y,n) qui donne un
polynôme qui approxime les données.
Dans l’instruction, on donne les valeurs respectives de x et de y, puis on indique par n le degré
du polynôme qui approximera les données.
Polyfit retourne alors le polynôme d’approximation par le biais d’un vecteur contenant les
coefficients des puissances descendantes.
Ainsi, en appelant :
Polyfit(x,y,1) : on obtient la pente de la droite et la valeur à l’origine
Polyfit(x ;y,2) : donne de plus le coefficient de l’équation quadratique
Pour tracer la courbe par la suite, il suffit d’utiliser les coefficients du polynôme pour le
recréer et le tracer à l’aide de plot ou de fplot
b) Basic fitting
Si on connaît préalablement le degré du polynôme que sont supposées former les données,
cette méthode est excellente. Par contre, si on ne connaît pas exactement le degré à employer,
une seconde manière d’approximer s’impose.
Cette manière consiste à tracer les données désirées.

Une fois la figure produite, il faut aller dans le menu Tools/Basic Fitting.
Par la fenêtre qui apparaît, il est possible de sélectionner le groupe de données désiré, de
sélectionner le ou les approximations désirées selon le degré du polynôme.
De plus, il est possible de montrer l’équation de chaque approximation.
On peut aussi indiquer sur un graphique la différence entre chaque donnée et sa valeur
calculée avec l’approximation.
En cliquant en bas à gauche de la fenêtre, les résultats numériques sont disponibles.
Dans ces résultats, on trouve la fonction théorique du polynôme, ses coefficients, ainsi que la
norme des différences entre chaque donnée et sa valeur calculée avec l’approximation.
De plus, on peut enregistrer ces donnés dans des variables disponibles dans le Workspace.
En étendant la fenêtre une seconde fois avec la flèche, on peut calculer f(x) d’une valeur ou
d’un groupe de valeurs, selon l’équation sélectionnée.
Enfin, si cette méthode a été employée et que vous désirez conserver le graphique, n’oubliez
pas de l’enregistrer afin de conserver tous les changements que vous venez d’apporter.
c) Application
Ø A partir des données suivantes correspondant à des mesures de courant pour différentes
valeurs de tension, créez un script "ohm.m" produisant :
o un graphique correspondant aux données expérimentales
o une courbe associée au calcul de la résistance à chaque point
o la courbe approximant ces résistances avec un polynôme d’ordre 7, sur le même
graphique, avec les légendes
% mesures tension, courant

DATA = [ 0 0
20.5 2.03
30.2 2.98
39.8 4.02
50.3 5.05
60.0 6.01
70.5 7.09
80.0 8.00
89.3 8.99
100.5 9.97 ];
• Testez aussi ces données avec les différentes approximations du " basic fitting"

TP 04

Transféré par

Droits d'auteur :

Formats disponibles

TP 04

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

TP 04

Transféré par

Droits d'auteur :

Formats disponibles

I.N.S.S.E.T. – U. P. J.V.

Programmation avec Matlab

Par exemple, en lieu et place de P( x) = x 4 − 12 x 3 + 25 x + 116 ,

on écrit le vecteur P = [1 − 12 0 25 116]

• A l’aide de la fenêtre d’aide, déterminez les fonctionnalités de la commande polyval

• Utilisez la fenêtre d’aide pour analyser la commande roots

Dans un script intitulé "poly2.m" :

• Utilisez la fenêtre d’aide pour analyser la commande conv

• Servez vous de la fenêtre d’aide pour découvrir la commande deconv

• Ouvrir trois fenêtres graphiques (commande figure)

Après avoir créez plusieurs graphiques, testez la commande open

c) Application aux polynômes

Ø Créez un script "graphepoly.m" réalisant les actions suivantes :

Ø Déterminez les extremums du polynôme P (x) sur l’intervalle [− 1, 3] .

Cette manière consiste à tracer les données désirées.

% mesures tension, courant

Vous aimerez peut-être aussi