|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
dx=4*Ra/(N-1);
dy=dx;
x=1:M;
x=(x-1)*dx-(M-1)/2*dx;
y=1:N;
y=(y-1)*dy-(N-1)/2*dy;
dy1=dy*d/f1;
dx1=dx*d/f1;
S3=zeros(M,N);
for nt=1:Nt
t(nt)=-3*tau+(nt-1)*dt;
for n=1:N
y1=(n-1)*dy1-(N-1)/2*dy1;
for m=1:M
x1=(m-1)*dx1-(M-1)/2*dx1;
S1=S0*exp(-t(nt)^2/2/tau^2);
%S1=fftshift(S1);
S2=S1.*exp(i*k*((x1-X).^2+(y1-Y).^2)/2/d);
S3(m,n)=S3(m,n)+trapz(trapz(S2))*dt;
end
end
end
SL=exp(i*k*d)/i/lamda/d*S3*dx*dy;
请大家帮忙看一下,这个菲涅尔衍射积分程序有什么问题,现在出来的结果能量不守恒。先谢谢啦! |
|