How to find out the position matrix from Intensity matrix?

1 次查看(过去 30 天)
Hi sir,
I sending following program by using this I am finding intensity in the focal plane for an X-polarized (linearly polarized light) beam. By using this program I got intensity in the focal plane. But I facing the problem to locate position of this intensity matrix. How we can get position of intensity, where it is located in the focal plane. Can anyone help me in this regard.
Here kx,ky position vectors in pupil plane.
dim=32;
lamda=531e-9;
k1=2*pi/lamda;
[kx,ky]=meshgrid(-1:2/(dim-1):1);
circ=sqrt(kx.^2+ky.^2)<1;
alp=asin(0.6);
k0=1/sin(alp);
k=512;
kz=sqrt(k0^2-(kx.^2+ky.^2));
Gx=((k0*ky.^2+kz.*kx.^2)./(k0*(kx.^2+ky.^2)));
Gy=sqrt(k0./kz).*((kz-k0).*kx.*ky)./(k0*(kx.^2+ky.^2));
Gz=sqrt(k0./kz).*(kx./k0);
Ex=fftshift(fft2(Gx.*circ,k,k));
Ix=Ex.*conj(Ex);
Ey=fftshift(fft2(Gy.*circ,k,k));
Iy=Ey.*conj(Ey);
Ez=fftshift(fft2(Gz.*circ,k,k));
Iz=Ez.*conj(Ez);
Itot=Ix+Iy+Iz;
figure(1)
imagesc(Itot);
colormap gray;
axis image;
figure(2)
plot(Itot(k/2,:))
  7 个评论
Walter Roberson
Walter Roberson 2013-1-5
For example, once you have calculated Itot, what does Itot(7, 19) represent? For example does it represent the strength of Zeta Beam that Sardath would use to transport Vath Sarn 7 hours and 19 minutes into the past?

请先登录,再进行评论。

回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Lighting, Transparency, and Shading 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by