我模拟的程序输出图与文献给出的图不一样,请问哪位大神可以帮我看看。
lambda=532e-9;
L=3e-2;
w0=3e-6;
d=20e-6;
j=sqrt(-1);
x1=270e-6;
x2=290e-6;
x3=310e-6;
x4=330e-6;
y4=0;
y3=0;
y2=0;
y1=0;
A=0:10:600;
B=0:10:600;
[a,b]=meshgrid(A,B);
m1=(exp(-(sqrt(-1).*2.*pi.*(x1.*a+y1.*b)/(lambda.*L)))).^2;
m2=(exp(-(sqrt(-1).*2.*pi.*(x2.*a+y2.*b)/(lambda.*L)))).^2;
m3=(exp(-(sqrt(-1).*2.*pi.*(x3.*a+y3.*b)/(lambda.*L)))).^2;
m4=(exp(-(sqrt(-1).*2.*pi.*(x4.*a+y4.*b)/(lambda.*L)))).^2;
n=exp(-(2.*pi.*w0/(lambda.*L)).^2.*(a.^2+b.^2));
I=m1.*n+m2.*n+m3.*n+m4.*n;
mesh(a,b,I)
axis=([0,600,0,600,-0.5,1]);
colormap(hsv)
这是我写的程序。模拟的是图Fig. 5b.
文献是下面这篇。![哪位大神可以帮我看看程序怎么改]()
untitled20140410.jpg |