Plane Wave Propagation (Young's Experiment)

11 views (last 30 days)
Bektas
Bektas on 28 Nov 2013
Commented: Nitesh Anerao on 22 Feb 2015
Hello everyone,
I have a question regarding 2slit Young's Experiment. I am trying to simulate this by using plane wave propagation. I assume you are familiar with the concept. What I do is I create a vector of size N with 1s on the slits and zeroes otherwise at z=0. This is my source. Now I want to propagate this to z=z0 by using e^(-1i*distance*kz) and kz=sqrt(k0-kx) as I am not propagating on y axis and assuming that ky=0. The reasoning is like this, taking the FFT of the source at z=0 and multiplying with my propagation exponential. And than taking IFFT at z=z0 to obtain the angular spectrum. The problem is my source is not correctly propagated as I have a problem with creating correct kx values and so the multiplication results with somewhat wrong. I know I have to rescale coordinates of FFT somehow but I dont know how. To put it more understandable
A(kx,z0)=A(kx,0)e^(j*z0*kz) is what I want to achieve A being the fourier transformed source at 0 and z0. Any help would be appreciated, I can provide more information if needed but I believe this is understandable and I assume someone else had tried to perform this experiment on Matlab.
Best Regards Bektaş
  1 Comment
Nitesh Anerao
Nitesh Anerao on 22 Feb 2015
So as I understand that you need propogation in only x-direction. I think kx is defined as kx =(w/c)= (2*pi*fx)/(c) =(2*pi)/lam ;
Now w=2*pi*fx ;fx=spatial frequency of wave (samples per unit length)in x=direction ;c=speed of sound or light ;lam = wavelength of the wave
You can somewhat refer answer by Rick Rosson in link below http://www.mathworks.com/matlabcentral/answers/13896-fftshift-of-an-image

Sign in to comment.

Answers (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!