diffraction of many particle.m

均匀分布多粒子引起的夫琅和费衍射且输出为三维图形

clear,clf
r=0.0003;           %小孔直径(m)
R=0.5;              %镜头半径(m)
%s=6*10e-6         %撒布总面积(m)
lamda=587.6e-009;  %可见光中心波长(m)
l=0.6;             %接收屏与衍射屏的距离(m)
DX=0.1;DY=0.1;
b=-0.5;d=0.5;c=200;  %b-区域最小值(m)，d-区域最大值(m)，c-区域格数
x=linspace(b,d,c);y=x;
[X,Y]=meshgrid(x,y);
strength0=zeros(200,200);
H0=zeros(200,200);
strength=zeros(200,200);
seta1=sin(atan(sqrt(X.*X+Y.*Y-R*R)/l));
z1=2*pi*R*seta1/lamda;
strength=strength+((2*besselj(1,z1))./z1).^2;
while ((X.*X+Y.*Y)<=R*R)
    any(real(strength)<=-Inf);
    strength=strength+ans;
end
%strength=1-strength0   %/I
surf(x,y,strength)
for number2=1:9              %(d-b)/DY-1
   for number1=1:9           %(d-b)/DX-1
      seta=atan(sqrt((X-(number1-(d-b)/2/DX)*DX).^2+(Y-(number2-(d-b)/2/DY)*DY).^2)/l);
      z=2*pi*r*sin(seta)/lamda;     %rou=ρ=r/(lamda*l)
      J1=besselj(1,z)
      H1=(2*J1)./z
      H0=H1+H0;
      strength1=((2*J1)./z).^2
      strength0=strength0+strength1
   end
end