Discussion :

[jimi567]

Bonjour à tous,

j'ai séché pendant longtemps et n'ai toujours pas trouvé la manière de produire un résultat identique à la fonction "surf" mais en coordonnées polaires (2D) (et donc créer de même un mesh dans les mêmes coordonnées).

En particulier je voudrais représenter la distribution de températures (avec interpolation, comme par exemple une interpolation de hermite) sur un disque et chaque noeud de température correspond à une maille, laquelle n'est évidemment pas rectangulaire. Je sais que je pourrais prendre un maillage rectangulaire très fin pour approcher la géométrie du disque mais je veux pouvoir le faire avec quelques mailles seulement...

J'avais trouvé des débuts de piste comme la fonction "patch" mais sans plus de certitude.

Pour plus de clarté, je joints une image du disque en question avec son maillage et une numérotation possible des noeuds.


D'avance un grand merci

[Jean Dumoncel]

Salut,

est-ce que ce lien correspond à ce que tu cherches :
3D Polar Plot?

[jimi567]

Salut,

j'avais déjà trouvé ce lien mais ce n'était pas vraiment ce que je cherchais

Par contre en faisant des recherches pour une autre partie du code j'ai trouvé finalement ce qu'il me fallait, il suffit de passer dans le plan cplx :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
%%%%%%%%%%%%%%%%%
[x,y]=meshgrid(-1:0.01:1,-1:0.01:1)
z=x+i*y;
[ch,ch]=contourf(x,y,-abs(z),0:-0.01:-1);
set(ch,'edgecolor','none');
axis equal
axis([-1 1 -1 1])
%%%%%%%%%%%%%%%%%

par exemple ; reste juste à évaluer correctement la fonction z en fonction de ma distribution de température...Un autre moyen que j'ai également trouvé (sur un forum du site de mathworks de mémoire) :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
%%%%%%%%%%%%%%%%%
r = [0 1 2 3];
T = [150 120 100 70];
maxr = max(r(:))./ 2;
 
for K = 1:length(r)
  cols(K,:) = hsv2rgb([T(K)./255 1 .8]);   
end
 
for K = length(r):-1:1
  rectangle('Curvature', [1 1], 'Position', [maxr-r(K)/2 maxr-r(K)/2 r(K)+eps r(K)+eps], 'FaceColor',cols(K,:))
end
axis equal
%%%%%%%%%%%%%%%%%

en passant par la fonction "rectangle"...

voilà merci quand même pour la réponse tschüss