Répondre à la discussion
Affichage des résultats 1 à 15 sur 15

Problème représentation phase



  1. #1
    tpscience

    Problème représentation phase


    ------

    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.

    -----

  2. Publicité
  3. #2
    tpscience

    Re : Problème représentation phase

    Le fichier fait 12Ko mais ne veut pas passer. Je vais retenter...

  4. #3
    tpscience

    Re : Problème représentation phase

    Voici l'image :

    phi(x).jpg

  5. #4
    Thouxify

    Re : Problème représentation phase

    Bonjour,
    Soit vous utilisez des lentilles de fresnel, sinon c'est un repliement de la phase. Dans ce cas il y a la fonction unwrap qui permet de dérouler la phase. A+
    Il y a trois sortes de personnes : ceux qui savent compter et ceux qui ne savent pas.

  6. #5
    tpscience

    Re : Problème représentation phase

    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...

  7. A voir en vidéo sur Futura
  8. #6
    Thouxify

    Re : Problème représentation phase

    A priori c'est la bonne méthode puisque les rayons sont perpendiculaires au front d'onde, quelle est la focale de la lentille? Et quelle est la divergence que vous obtenez?
    A+
    Il y a trois sortes de personnes : ceux qui savent compter et ceux qui ne savent pas.

  9. Publicité
  10. #7
    tpscience

    Re : Problème représentation phase

    J'ai une focale de 0.1mm, pour des rayons de lentilles de 50 microns.

    Mais mon profil de phase est bien courbe, donc les rayons devraient avoir une direction différente...

    Voici ma partie de code :

    Code:
    for k=1:number_lens
        phi(:,k) = unwrap(angle(field(:,256/2,k)));
    end
    
    pas = 2*Rayon_lens/256;
    
    gradient_phi(:,:) = gradient(phi(:,:),pas);
    
    angle_i = acos(gradient_phi/k);
    ...?

    Merci encore.

  11. #8
    Thouxify

    Re : Problème représentation phase

    Ce ne serait pas plutôt
    gradient_phi(:,k) = gradient(phi(:,k),pas);
    ?
    Sinon il me semble que le gradient va se faire sur toutes les directions de la matrice non?
    Je ne comprend pas la dernière ligne, pourquoi diviser par le nombre de lentilles?
    A+
    Il y a trois sortes de personnes : ceux qui savent compter et ceux qui ne savent pas.

  12. #9
    obi76

    Re : Problème représentation phase

    Bonjour,

    d'autant plus que diviser par k en dehors de la boucle, c'est à éviter. Soit vous divisez par number_lens, soit vous ne faites rien, mais vous n'utilisez pas la variable d'une boucle en dehors de la boucle (surtout que là je crois que votre résultat sera tout simplement faux car k vaudra number_lens+1).
    Dernière modification par obi76 ; 21/07/2013 à 10h52.
    Paradoxalement, ce sont les débats stériles qui se reproduisent le plus.

  13. #10
    tpscience

    Re : Problème représentation phase

    Bonjour,

    C'est un pb de déclaration de ma part dans le script que je vous ai présenté. Ici j'ai bien déjà fini ma boucle sur la dernière ligne, et k ne représente que le nombre d'onde ...

  14. #11
    Thouxify

    Re : Problème représentation phase

    Je ne comprend pas la dernière ligne, j'aurais mis gradient=tan(angle) non?
    Enfin angle=atan(gradient) non?
    Il y a trois sortes de personnes : ceux qui savent compter et ceux qui ne savent pas.

  15. #12
    tpscience

    Re : Problème représentation phase

    Eh bien, je suis parti de l'équation eikonale, qui me dit que , où est la normale à la surface d'onde en un point donné, n l'indice du milieu (n=1 ici pour simplifier) et je considère mon axe de propagation suivant .
    Ainsi,



    A partir de quoi, je multiplie de part et d'autre part , et tenant compte des propriétés du produit scalaire de deux vecteurs :



    D'où mon arccos...

  16. Publicité
  17. #13
    Thouxify

    Re : Problème représentation phase

    J'avais trouvé la tangente par construction géométrique mais j'ai du me tromper quelque part... sinon il me semble que le gradient doit etre calculé dans la direction perpendiculaire à la direction de propagation, c'est peut être l'origine de vos résultats en pi/2
    A+
    Il y a trois sortes de personnes : ceux qui savent compter et ceux qui ne savent pas.

  18. #14
    tpscience

    Re : Problème représentation phase

    En fait, la figure obtenue est du style :

    theta.jpg

    Pour la partie de code liée au tracé :

    Code:
    x = Rayon_lens*ones(size(angle_i));
    [u,v] = pol2cart(angle_i,x);
    feather(u,v);
    
    figure
    feather(angle_i)
    Je ne vois pas trop comment avancer cela...?

  19. #15
    Thouxify

    Re : Problème représentation phase

    Re,
    pourquoi avoir multiplié x par rayon_lens? du coup votre x est un vecteur ou chaque point vaut Rayon_lens
    ce ne serait pas plutôt
    x=Rayon_lens.*linspace(0,1,siz e(angle));
    Il y a trois sortes de personnes : ceux qui savent compter et ceux qui ne savent pas.

Sur le même thème :

Discussions similaires

  1. Réponses: 3
    Dernier message: 21/03/2013, 22h00
  2. [Biologie végétale] Application pratique pour phase sombre et phase claire de la photosynthèse
    Par tnarol dans le forum Biologie
    Réponses: 0
    Dernier message: 25/09/2011, 13h22
  3. Réponses: 11
    Dernier message: 26/08/2009, 19h08
  4. Problème de représentation
    Par tomass' dans le forum Physique
    Réponses: 2
    Dernier message: 08/10/2008, 19h00
  5. représentation de la modulation de phase octovalente
    Par simbad2 dans le forum Internet - Réseau - Sécurité générale
    Réponses: 2
    Dernier message: 24/03/2008, 21h50