Tracer un contour plot

Version imprimable

1 pièce(s) jointe(s)

Bonjour,

J'essaie de faire un contour plot pour avoir un diagramme comme ceci :

Pièce jointe 158236

J'ai mes données sous un fichier Excel que j'ai réussi à exporter sans trop de souci.

Donc voila le code que j'ai écrit :
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 clear all; close all; clc; EngineSpeed=xlsread('Data.xls','a2:a219'); Torque=xlsread('Data.xls','c2:c219'); be=xlsread('Data.xls','al2:al219'); figure contour(EngineSpeed,Torque,be) xlabel('Regime moteur [1/min]'); ylabel('Couple [Nm]'); shading interp colorbar
Lorsque j’exécute le code il me met cette erreur :
Code:

1 2 3 4 5 ??? Error using ==> contour at 73 Z must be size 2x2 or greater. Error in ==> CartoMoteur at 13 contour(EngineSpeed,Torque,be)
Que faut-il que je modifie dans mon code ?

Je ne comprends pas l'erreur car j'ai le même nombre de variable dans mes colonnes excel.

Je précise, cette fois ci je suis bien sous matlab :P

Merci :zoubi:

08/10/2014, 13h44
Invité
Bonjour,
Citation:
Envoyé par Documentation de la fonction contour

contour(X,Y,Z), contour(X,Y,Z,n), and contour(X,Y,Z,v) draw contour plots of Z using X and Y to determine the x and y values.

If X and Y are vectors, then length(X) must equal size(Z,2) and length(Y) must equal size(Z,1). The vectors must be strictly increasing or strictly decreasing and cannot contain any repeated values.
08/10/2014, 13h44
Jerome Briot

Les données importées sont sous forme de vecteurs alors qu'il faut des tableaux comme ceux générés par la fonction meshgrid (par exemple)

Pourrais-tu nous fournir un fichier xls ?

1 pièce(s) jointe(s)

Je t'ai mis le fichier un pièces-jointes. Celui-ci est différent, car l'original est trop lourd. Donc j'ai également ré-ecrit quelques détails dans le programme notamment pour l'ouverture :
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 clear all; close all; clc; EngineSpeed=xlsread('Data.xlsx','a2:a219'); Torque=xlsread('Data.xlsx','b2:b219'); be=xlsread('Data.xlsx','c2:c219'); figure contour(EngineSpeed,Torque,be) xlabel('Regime moteur [1/min]'); ylabel('Couple [Nm]'); shading interp colorbar
L'erreur est la même sauf le nom du fichier que j'ai changer :oops:
Code:

1 2 3 4 5 ??? Error using ==> contour at 73 Z must be size 2x2 or greater. Error in ==> Essais at 11 contour(EngineSpeed,Torque,be)
Donc d'après ce que je comprends il faut que j'utilise des données sous forme de tableaux avec contour et non sous forme de vecteur ?

1 pièce(s) jointe(s)

Le problème est que le domaine est irrégulier.

Essaie avec cette contribution : Tricontf by Rich Pawlowicz
Code:

1 2 3 4 5 6 7 8 9 tri = delaunay(EngineSpeed, Torque); figure subplot(211) scatter(EngineSpeed, Torque, 36, be, 'fill') title('scatter') subplot(212) tricontf(EngineSpeed, Torque, tri, be, 5); title('tricontf')

09/10/2014, 09h19
amelianna

Qu'est ce que sa veut dire domaine irrégulier ? Je suis comprends plus rien...:(

J'ai téléchargé le fichier tricontf et j'ai 3 fichier en .mat. Dans quel dossier de matlab dois-je les placer ?
09/10/2014, 09h46
Jerome Briot

Citation:

Envoyé par amelianna

Qu'est ce que sa veut dire domaine irrégulier ?

Les points de mesure ne sont pas répartis régulièrement selon l'axe des ordonnées et les bords du domaine ne forme par un rectangle.

Citation:

Envoyé par amelianna

J'ai téléchargé le fichier tricontf et j'ai 3 fichier en .mat. Dans quel dossier de matlab dois-je les placer ?

Ce sont des fichier .m, pas des fichiers .mat. Place-les dans le même dossier que celui qui contient le code que tu as écrit.
09/10/2014, 10h43
amelianna

Peux tu m'expliquer le code que tu as écrit ?

a quoi sert

Code:

subplot(212)

?

Code:

tricontf(EngineSpeed, Torque, tri, be, 5);

est la fonction comme contour.

A quoi correspondent les paramètres tri et 5 ?
09/10/2014, 11h02
Jerome Briot
Citation:

Envoyé par amelianna

a quoi sert

Code:

subplot(212)

?

Je te renvoi vers la documentation :

Code:

doc subplot

J'ai simplement ajouté le graphique du haut pour te montrer que le domaine n'est pas une grille rectangulaire.

Citation:

Envoyé par amelianna

Code:

tricontf(EngineSpeed, Torque, tri, be, 5);

est la fonction comme contour.

Oui, sauf qu'elle a besoin d'un maillage triangulaire généré à partir des données initiales (EngineSpeed et Torque)

Citation:

Envoyé par amelianna

A quoi correspondent les paramètres tri et 5 ?

La variable tri contient les indices de la triangulation générée avec la fonction delaunay :

Code:

doc delaunay

Tu peux visualiser ce maillage en faisant :
Code:

1 2 figure trimesh(tri, EngineSpeed, Torque);
Le dernier argument de tricont fonctionne comme celui de la fonction contour :
Code:

1 2 help tricont doc contour
09/10/2014, 11h24
amelianna

Je comprends mieux le premier graphique maintenant :mrgreen:

Comment puis-je interpoler les valeurs non existante avec matlab ? (je vais chercher un peu :roll:)

Le dernier argument de tricont sert juste à la précision ? J'ai changé plusieurs fois de valeurs et selon la valeur, il y a plus ou moins de courbes sur le graphique du bas.
09/10/2014, 11h31
Jerome Briot

Citation:

Envoyé par amelianna

Comment puis-je interpoler les valeurs non existante avec matlab ? (je vais chercher un peu :roll:)

Fais une recherche avec "griddata" sur ce forum, tu trouveras des discussions en rapport avec ce que tu cherches à faire.

Citation:

Envoyé par amelianna

Le dernier argument de tricont sert juste à la précision ? J'ai changé plusieurs fois de valeurs et selon la valeur, il y a plus ou moins de courbes sur le graphique du bas.

Si le dernier argument est un scalaire, alors il donne le nombre de contour à tracer.
09/10/2014, 13h04
amelianna

1 pièce(s) jointe(s)
J'ai essayer avec griddata. J'ai trouvé sa sur le forum.
Code:

1 2 3 4 5 [xg,yg] = meshgrid(x,y); zg = griddata(x,y,z,xg,yg,'cubic'); figure(3); mesh(xg,yg,zg); hold on; plot3(x,y,z);
avec
x=EngineSpeed
y=Torque
Z=be

Sauf que quand j'affiche le domaine pour voir s'il est rectangulaire cela me met la même chose que précédemment. Donc je crois que la fonction ne fait rien.

J'ai mis un tableau en lien. (ce sont les mêmes valeurs que le excel précédent, c'est juste plus lisible et plus facile à comprendre )

En fait ce que je voudrais c'est faire une interpolation des valeurs que je n'ai pas dans le tableau. Faut-il que j'utilise griddata ?

Voici comment générer la grille :

Code:

1
2
3
4
5
6
7
8
9
xi = linspace(min(EngineSpeed(:)), max(EngineSpeed(:)), 18);
yi = linspace(min(Torque(:)), max(Torque(:)), 10);
[xi, yi] = meshgrid(xi, yi);
 
figure
subplot(211)
scatter(EngineSpeed, Torque, 36, be, 'fill')
subplot(212)
plot(xi, yi, 'r*')

Sa marche bien sa rempli tout le domaine d’après ce que je vois.

D'après ce que j'ai compris, matlab aurait interpoler toutes les valeurs que je n'ai pas dans le tableau excel ? (case vide c'est a dire le be, je ne cherche pas a extrapoler EngineSpeed et Torque mais seulement be sur les différents points de fonctionnement que je n'ai pas)

Pour la suite j'utilise :
x=EngineSpeed
y=Torque
z=be

Voila l'erreur qu'il me met quand je met :

Code:

plot3(xi,yi,z)
Code:

1 2 3 4 5 ??? Error using ==> plot3 Vectors must be the same lengths. Error in ==> Essais at 58 plot3(xi,yi,z)
Et quand j'utilise :

Code:

tricontf(xi, yi, tri2, z, 8);

J'ai :
Code:

1 2 3 4 5 6 7 8 9 10 ??? Index exceeds matrix dimensions. Error in ==> tricont at 93 xx=[ Xp(M(:,[1 2 3 1])')]; Error in ==> tricontf at 81 [Xp,Yp,M,Zp,nv,CS,xx,yy,zz]=tricont(args{1:nin}); Error in ==> Essais at 40 tricontf(xi, yi, tri2, z, 8);
Je ne comprends pas pourquoi il y a des erreurs. Si on vient d'interpoler les valeurs manquantes de be, le domaine devrait être rectangulaire et cela ne devrait pas poser de problème non ?

09/10/2014, 16h15
Jerome Briot

Citation:

Envoyé par amelianna

Comment puis-je interpoler les valeurs non existante avec matlab ?

Pour éviter tout malentendu, pourrais-tu préciser ce que tu entends exactement par "valeurs manquantes" ? Où devraient-elles être situées ?
10/10/2014, 09h17
amelianna

1 pièce(s) jointe(s)

Bonjour,

Dans le tableau en pièce jointe, ce sont les cases en rouges dont je n'ai pas les valeurs que je voudrais interpoler.

Est-ce possible ? Si oui avec l'aide de quelle fonction ?