[zukaldl]

以下程序实现不了~~
不知道哪里出了问题
请指教：
Lambda=500e-9;
theta=1e-9;
X=0.01;
Ni=8000;
ds=linspace(0,0.02,Ni);
for k= 1:Ni
y(k)=ds(k)/sin(theta);
Delta=2*ds(k)+Lambda/2;
Phi=2*pi*Delta/Lambda;
B(k,=4*COS(Phi/2).^2;
end
figure(gcf);
NCLevels=250;
Br=(B/4.0)*NCLevels;
subplot(1,2,1);
image(x,Y,Br);
colormap(gray(NCLevels));
subplot(1,2,2);
plot(B(,Y);