Forum |  HardWare.fr | News | Articles | PC | S'identifier | S'inscrire | Shop Recherche
1788 connectés 

  FORUM HardWare.fr
  Programmation
  C++

  Pb code calcul de la diffraction fraunhofer par methode numerique

 
 


 Mot :   Pseudo :  
 
Bas de page
Auteur Sujet :

Pb code calcul de la diffraction fraunhofer par methode numerique

n°1667526
Anodo-kun
Posté le 08-01-2008 à 10:23:54  profilanswer
 

Bonjour à tous,
 
Voila je post ce message car j'ai un code en langage c à réaliser qui calcule par une méthode numérique la diffraction de fraunhofer. Voici un sujet similaire au mien : http://www.impmc.jussieu.fr/impmc/ [...] 5/td5.html( attention monsujet ne comporte que la partie I).
J'ai utilisé la méthode de simpson pour faire mon calcul d'intégrale. Mon problème ce que mon code me rend une intensité nulle et je ne comprend pas pourquoi ??? :??:. J'ai utilisé 4 fonctions ( F1, F2, F3, F4) pour séparer la partie cosinus et la partie sinus et aussi pour séparer les x des y.
 
Voila si quelqu'un pourrait m'aider un peu car je patauge depuis un petit moment ?!?! :pt1cable: ^^
Merci d'avance !  
 
Voici la code que j'ai réalisé :

Code :
  1. #include <stdio.h>
  2. #include <math.h>
  3.    
  4. #define pi M_PI
  5. /* Longueur d'onde de la source lumineuse */
  6. #define lambda 500e-9
  7. /* Focale de la lentille */
  8. #define f 0.15
  11.        
  12. double F1(double x, double X)
  13. {
  14.     return cos((2*pi/(lambda*f))*(x*X));
  15. }
  17. double F2(double x, double X)
  18. {
  19.     return sin((2*pi/(lambda*f))*(x*X));
  20. } 
  21. double F3(double y, double Y)
  22. {
  23.     return cos((2*pi/(lambda*f))*(y*Y));
  24. }
  26. double F4(double y, double Y)
  27. {
  28.     return sin((2*pi/(lambda*f))*(y*Y));
  29. } 
  30. double simpson(int N, double x1, double x2,double y1, double y2, double X, double Y)
  32. {
  33.     double xx,yy, h1, IappF10, IappF11, IappF12, IappF1, h2, IappF20, IappF21, IappF22, IappF2, IappF30, IappF31, IappF32, IappF3, IappF40, IappF41, IappF42, IappF4,Iapp;
  34.     int NN, i;
  35.    
  36.     // Etape 1
  37.     h1 = (x2 - x1) / N;
  38.     h2 = (y2 - y1) / N;
  40.     // Etape 2
  41.     IappF10 = F1(x1,X) + F1(x2,X);
  42.     IappF11 = 0.0;
  43.     IappF12 = 0.0;
  44.     IappF20 = F2(x1,X) + F2(x2,X);
  45.     IappF21 = 0.0;
  46.     IappF22 = 0.0;
  47.     IappF30 = F3(y1,Y) + F3(y2,Y);
  48.     IappF31 = 0.0;
  49.     IappF32 = 0.0;
  50.     IappF40 = F4(y1,Y) + F4(y2,Y);
  51.     IappF41 = 0.0;
  52.     IappF42 = 0.0;
  53.     // Etape 3
  54.     NN = N -1;
  55.     for (i=1; i<=NN; i++)
  56.     {
  57. // Etape 4
  58. xx = x1 + i*h1;
  59. yy = y1 + i*h2;
  60.    
  61. // Etape 5
  62. if ((i%2) == 0)
  63. {
  64.     IappF12 = IappF12 + F1(xx,X);
  65.     IappF22 = IappF22 + F2(xx,X);
  66.     IappF32 = IappF32 + F3(yy,Y);
  67.     IappF42 = IappF42 + F4(yy,Y);
  68. }
  69. else
  70. {
  71.     IappF11 = IappF11 + F1(xx,X);
  72.     IappF21 = IappF21 + F2(xx,X);
  73.     IappF31 = IappF31 + F3(yy,Y);
  74.     IappF41 = IappF41 + F4(yy,Y);
  75. }
  76.     }
  77.    
  78.     // Etape 6
  79.     IappF1 = (IappF10 + 2.0 * IappF12 + 4.0 * IappF11) * h1 / 3.0;
  80.     IappF2 = (IappF20 + 2.0 * IappF22 + 4.0 * IappF21) * h1 / 3.0;
  81.     IappF3 = (IappF30 + 2.0 * IappF32 + 4.0 * IappF31) * h2 / 3.0;
  82.     IappF4 = (IappF40 + 2.0 * IappF42 + 4.0 * IappF41) * h2 / 3.0;
  83.     // Etape 7
  85.     Iapp = (IappF1*IappF2-IappF3*IappF4)*(IappF1*IappF2+IappF3*IappF4)+(IappF3*IappF2+IappF1*IappF4)*(IappF3*IappF2+IappF1*IappF4);
  87.     // Etape 8
  88.     return (Iapp);
  89. }
  91. int main(int argc, char *argv[])
  92. { 
  94.     int i,n;
  95.     double X,Y,x1,x2,y1,y2;
  96.    
  97.     x1=-0.000025;       /* Borne inferieure sur x*/
  98.     x2=+0.000025;       /* Borne superieure sur x*/
  99.     y1=-0.0025;         /* Borne inferieure sur y*/
  100.     y2=-0.0025;         /* Borne superieure sur y*/
  101.     n=100;              /* Nombre d'iteration */
  102.    
  103.     double I;
  105. /* Creation d'un fichier de donnée contenant l'intensité I(X,Y) demandé */
  106.     FILE *fp;
  107.     fp=fopen("intensite.data","w" );
  108.     for(X=0 ; X<=49 ; X=X+0.1)
  109.     {
  110. for(Y=0; Y<=49; Y=Y+0.1)
  111. {
  112.      I=simpson(n,x1,x2,y1,y2,X,Y);
  113.     fprintf(fp,"%lf %lf %lf\n",I,X,Y);
  115. }
  116.     }
  117.     fclose(fp);
  119. }


mood
Publicité
Posté le 08-01-2008 à 10:23:54  profilanswer
 

n°1667532
Joel F
Real men use unique_ptr
Posté le 08-01-2008 à 10:31:21  profilanswer
 

ta fonction simpson me parait douteuse à vue de nez.
Tu l'as sort de ou ?
 
 
Et accessoirement, c'est du C pas du C++


Message édité par Joel F le 08-01-2008 à 10:54:46
n°1667663
Anodo-kun
Posté le 08-01-2008 à 12:42:43  profilanswer
 

voici le lien ou j'ai trouvé ma fonction simpson que j'ai ensuite adapté à mon pb :  
http://www.cppfrance.com/codes/INT [...] 27783.aspx

n°1667691
nargy
Posté le 08-01-2008 à 13:42:09  profilanswer
 

Code :
  1. #     y1=-0.0025;         /* Borne inferieure sur y*/
  2. #     y2=-0.0025;         /* Borne superieure sur y*/


si c'est une intégration sur un plan, avec ces bornes ça fera 0 pour toute fonction.

n°1667773
Anodo-kun
Posté le 08-01-2008 à 15:04:22  profilanswer
 

ah ouais merci erreur de frappe je suppose :p
je vais corriger ca !

n°1668117
neg'gwada
rafix 971
Posté le 09-01-2008 à 02:24:20  profilanswer
 

theoriquement ce genre de calcul est possible à coup de transformation de fourier, encore plus pour le cas fraunhofer car tu fais des hypothèses sur le fait que tu te trouves à r>>lambda. Passer par du simpson ne semble pas etre la meilleur idée! (apres j'ai pas vraiment tout lu en detail, je parle juste par experience)  [:spamafote]


---------------
--- WinSplit Revolution ---
n°1668696
Anodo-kun
Posté le 10-01-2008 à 09:54:33  profilanswer
 

ouais va falloir que je modifie qque chose car malgré la correction sur les bornes de y ca me rend toujours une intensité nulle ...

n°1668716
Anodo-kun
Posté le 10-01-2008 à 10:46:13  profilanswer
 

je viens d'installer la librairie fftw3. Quelqu'un aurait un site ou un document qui m'expliquerait comment l'utiliser ?? mci d'avance


Aller à :
Ajouter une réponse
  FORUM HardWare.fr
  Programmation
  C++

  Pb code calcul de la diffraction fraunhofer par methode numerique

 

Sujets relatifs
Methode template intelligenteProblème code VBA
calcul coordonnées pixels d'un droite[PHP] Du code dans un echo
Un menu html qui se répéte dans plusieur page mais un seule code ![C++] undefined reference to WinMain@16 sous Code::Blocks
[Delphi] Problème dans un code / N'affiche jamais la variableInsérer une image DANS le code html
Code windows media player 
Plus de sujets relatifs à : Pb code calcul de la diffraction fraunhofer par methode numerique

Forum MesDiscussions.Net, Version 2010.2
(c) 2000-2011 Doctissimo
Page générée en 0.067 secondes

Copyright © 1997-2022 Hardware.fr SARL (Signaler un contenu illicite / Données personnelles) / Groupe LDLC / Shop HFR