关于matlab的多缝衍射代码注释

求大神给代码注释
clear
lamda = 500e-9;
N = 2;
a = 2e-4; D = 5;d= 5*a;
ym = 2*lamda*D/a; xs = ym;
n=1001;
ys = linspace(-ym,ym,n);
for i = 1:n
sinphi = ys(i)/D;
alpha = pi*a*sinphi/ lamda;
beta = pi*d*sinphi/ lamda;
B(i,:) = (sin(alpha)./alpha).^2.*(sin(N*beta)./sin(beta)).^2;
B1 = B/max(B);
end
NC = 256; %确定灰度的等级
Br = (B/max(B))*NC;
subplot(1,2,1)
image(xs,ys,Br);
colormap(hot(NC)); %色调处理
subplot(1,2,2)
plot(B1,ys,'k');