Code de diffraction

Version imprimable

Code de diffraction

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 :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
 
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.

Les variables OPL et intensity sont des vecteurs :
Code:

1 2 3 4 5 >> whos OPL intensity Name Size Bytes Class Attributes OPL 1x2000 16000 double intensity 1x2000 16000 double
La ligne suivante pose donc problème :

Code:

imagesc(OPL(1:k-1),intensity(1:k-1))

Qu'espères-tu faire avec la fonction imagesc qui prend classiquement en entrée une seule matrice à 2 dimensions.

29/09/2013, 17h00
tpscience

Bonjour,

J'espérais juste faire une représentation 2D du vecteur d'intensité...
29/09/2013, 17h06
Jerome Briot

C'est assez vague... pourrais-tu être plus précis ou mieux nous montrer un exemple de ce que tu cherches à obtenir ?

Pourquoi ne pas utiliser plot comme un peu plus loin dans ton code ?

Une remarque au passage, les crochets sont inutiles dans la ligne suivante :

Code:

x = [0:.1:pi];
29/09/2013, 17h13
tpscience

1 pièce(s) jointe(s)

Effectivement, en fait j'essaie de représenter une figure de diffraction par un objet circulaire. La forme est du type de l'image jointe.
01/10/2013, 16h28
Jean Dumoncel

Bonjour,

tu ne peux pas représenter intensity sous cette forme puisque c'est un vecteur, il faut revoir le code si tu veux que le domaine de définition soit en 2D.