Discussion :

[tpscience]

Bonjour à tous,

J'ai écris un code MATLAB nécessitant l'expression d'une fonction discrétisée sur une grille (x,y) de points (512*512). J'aimerais passer cette représentation en polaire. Or, je ne vois pas trop quelle démarche suivre. Il faut, j'imagine, commencer par passer ma grille cartésienne en polaire...?


Merci

[tpscience]

Je vous transmet la partie de mon code me permettant de générer ma grille cartésienne :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
 
Np = 512;
max_grille = 2*G;
pas = max_grille/Np;
 
vec_np = 1:Np;
vec_axe = -max_grille/2 + pas/2 + (vec_np-1)*pas;
vec_axe_fft = -1/(2*pas) + (vec_np-1)*1/(Np*pas);
 
[X,Y] = meshgrid(vec_axe);
[FFT_X,FFT_Y] = meshgrid(vec_axe_fft);

...

[Invité]

Bonjour,

La fonction cart2pol te conviendrait-elle ?

[souviron34]

x = r cos a
y = r sin a

[tpscience]

Bonjour,

La fonction cart2pol te conviendrait-elle ?

Bonjour,

J'avais déjà effectivement essayé cette fonction en rajoutant :

Code :

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

1
2
3
 
[X,Y] = meshgrid(vec_axe);
[X,Y] = cart2pol(X,Y);

Mais la représentation graphique ensuite n'est pas bonne.

[Invité]

Dans ce cas, dis-nous plus précisément ce que tu souhaites faire.

[souviron34]

si (xc, yc) sont les coordonées du centre (si le pole est au centre)

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
pour i = 0  jusqu'à i < Ny   par pas de dy
   pour j = 0  jusqu'à j < Nx   par pas de dx
 
   x = j*dx - xc
   y = i*dy - yc
 
   r = sqrt( x^2 + y^2)
   a = atan2 ((y-yc), (x-xc)
 
   plot ( r, a )
 
 
   fin pour
fin pour

[tpscience]

Eh bien j'aimerais juste passer mon code cartésien du deuxième post en polaire, mais...

[Invité]

Eh bien j'aimerais juste passer mon code cartésien du deuxième post en polaire, mais...

Mais la représentation graphique ensuite n'est pas bonne.

En quoi ?

[tpscience]

Je ne retrouve pas de symétrie cylindrique, on dirait que la représentation se fait toujours en cartésien...

[Jerome Briot]

Pour tracer en coordonnée polaires, il faut utiliser la fonction polar

[tpscience]

Pour être plus explicite, et sans trop rentrer dans le détail du programme qui est très long : je cherche à représenter une somme de deux faisceaux, un gaussien auquel je somme un faisceau d'ordre supérieur, ici un laguerre-gauss d'ordre quelconque. Pour faire cela, j'avais déjà écris tout mon code en cartésien pour générer le faisceau gaussien. Or, maintenant j'aimerais rajouter cette partie en laguerre-gauss que je génère à partir de paramètres polaires. Je pense donc que sommer deux contributions dans deux systèmes de coordonnées différentes n'est pas adéquat !

D'où le fait que je cherche à passer dès le début de mon programme en polaire, et cela correspond aux lignes postées sur mon post #2.

Si je trace juste un mode laguerre-gauss d'ordre mn=10, et que je laisse tout en cartésien j'obtiens la figure jointe. Il y a sur celle-ci un pb qui est bien entendu que la symétrie circulaire n'est pas respectée, le profil semble étalé. Il me semble que cela est dû au fait que le profil correct ne soit enregistré que sur une des deux dimensions...

[Invité]

Si j'ai bien compris, tu peux t'inspirer de ce message.

[tpscience]

Merci beaucoup, je vais regarder cela !