Code de diffraction

[tpscience]

Bonjour à tous,

J'ai un problème dans un code de diffraction écrit en Matlab, concernant apparemment la fonction 'imagesc', mais je ne vois pas trop ce qui cloche...?

Voici la partie de code concernée :

clear
N=300;
format compact
clf
hold off
start=cputime;
x=[0:.1:pi];  
z=zeros(size(x));
StartPLD=6.2;
EndPLD=6.35;
Increment=.002;
intensity=zeros(1,2000);
OPL=zeros(1,2000);
k=1;
figure(2);clf;
for pld=StartPLD:Increment:EndPLD,  
  z=zeros(size(x));
  a=0;
  for j=1:N, 
     y=sin(3.*x+a);
     z=z+y;   
     a=a+pld;
  end
  intensity(k)=sum(z.*z);  
  OPL(k)=pld./(2*pi);
  imagesc(OPL(1:k-1),intensity(1:k-1))
  ylabel('Observed irradiance  (Mean-square of sum of all reflections)')
  xlabel('Pathlength difference between adjacent grooves, in wavelengths')
  title(['First-order diffraction pattern for grating with '  num2str(N)  ' grooves.'])
  
  k=k+1;
end
hold off
figure(2)
ElapsedTime=cputime-start
figure(3)
plot(OPL(end-1),intensity(end-1))

Merci pour vos suggestions.
-----

[invite1c6b0acc]

Bonjour,

J'ai un problème dans un code de diffraction écrit en Matlab, concernant apparemment la fonction 'imagesc', mais je ne vois pas trop ce qui cloche...?

Et si tu nous disais quel problème tu as ? Parce qu'en ne sachant ni ce que ton programme doit faire, ni quelle erreur se produit, on va mettre plus de temps à trouver ...

[tpscience]

Eh bien comme énoncé dans le premier post c'est un programme qui doit me calculer le profil d'intensité d'un faisceau lumineux diffracté par un objet.
Je cherche donc à représenter cette figure de diffraction...

[tpscience]

Petit up...

[A voir en vidéo sur Futura]

[invite1c6b0acc]

En ce qui me concerne, j'attends toujours de savoir quel erreur tu as ...