Transformada de Fourier

Aqui vamos a analizar la transformada de Fourier. Comenzamos con la relacion de dispersion. La guia tiene nh=3.6

(set! geometry-lattice (make lattice (size 16 8 no-size)))

(set! geometry 
 (list
   (make block 
    (center 0 0) 
    (size infinity 1 infinity)
    (material (make medium (epsilon 12.96)))
   )
  )
 )

(set! sources 
 (list
  (make source
   (src (make continuous-src (frequency 0.1)))
   (component Ez)
   (center -7 0)
  )
 )
)

(set! pml-layers (list (make pml (thickness 1.0))))

(set! resolution 10)

(run-until 200
  (at-beginning output-epsilon)
  (at-end output-efield-z)

A continuacion planteamos el programa

load a00-ez-000200.00.h5
whos

Om    = 0.1;

nl    = 1.0;
nh    = 3.6;

Ql    = Om*nl;
Qh    = Om*nh;

Qi    = Ql;
Qf    = Qh;
NQ    = 100;
dQ    = (Qf-Qi)/NQ;

for iQ=1:NQ
Q  = Qi+iQ*dQ;
QV(iQ)=Q;

Er    = 0;
Ei    = 0;

for ix=20:150
Er = Er+ez(40,ix)*cos(2.0*pi*Q*(ix/10));
Ei = Ei+ez(40,ix)*sin(2.0*pi*Q*(ix/10));

end

E(iQ)  = Er*Er+Ei*Ei

end