求完善程序matlablam=500e-9;N=2;a=2e-4;z=5;d=3*a;xm=5*lam*z/a;y0=xm;n=1001;x0=linspace(-xm,xm,n);for i=1:nsinphi=x0(i)/z;alpha=pi*a*sinphi/lam;beta=pi*d*sinphi/lam;B(i,:)=(sin(alpha)./alpha).^2.*(sin(N*beta)./sin(beta)).^2;B1=B/max(B);end取B1的