Discussion :

[tpscience]

Bonjour à tous,

Je cherche à représenter le profil d'amplitude et de phase d'un faisceau gaussien en propagation dans un ensemble de lentilles minces convergentes.
Lorsque je trace le profil d'amplitude au niveau de quelques lentilles, pas de pb, j'ai de belles gaussiennes. En revanche, quand je trace la phase, soucis (cf pièce jointe) !
Ma représentation est faite avec Matlab.

Auriez-vous une idée du pb ? Numérique...?

Merci.

[Invité]

Bonjour,

Le "problème" vient du fait que la phase est confinée dans l'intervalle [-2pi 2pi] donc toute valeur en sortant y sera ramenée en modulo 2pi.
Tu as la fonction unwrap pour contrer cet effet.

[tpscience]

Bonjour,

Merci beaucoup, en effet, le pb venait de là !

Un nouveau problème est survenu. J'aimerai tracer pour chaque élément de surface composant mon profil transverse, la direction du rayon incident associé. Je les supposent orthogonaux à chaque tangente respective.

Donc pour cela, je calcule mon profil de champ à une lentille, j'en sors le profil de phase associé (discrétisé sur 256 points) avec la fonction 'angle'. Ensuite, je cherche à calculer le gradient de cette phase en utilisant un pas défini par la largeur "physique" de mon axe transverse divisé par 256. Pour cela j'utilise la fonction 'gradient(phi,pas)'. Et enfin, je cherche à sortir l'angle à partir du gradient calculé (soit un calcul de pente).

Le pb vient du fait qu'en cherchant à représenter ça toutes mes directions sont quasiment liées à un angle de pi/2...