改变曲率半径
lamda = 0.000632;
[X,Y] = meshgrid(-2:0.01:2);
for R=3200:-500:200
I=cos(pi/lamda*(X.^2+Y.^2)/(R)+pi/2);
imshow(I) %将光强显示为灰度值
pause(0.8)
end
双缝干涉
d是双缝的距离
ym=1.25;
y=linspace(-ym,ym,101);
z=1000;
lambda=5e-4;
for d=0.2:0.2:3
L1=sqrt((y-d/2).^2+z^2);
L2=sqrt((y+d/2).^2+z^2);
phi = 2*pi*(L2-L1)/lambda;
I = 4*(cos(phi/2)).^2;
B = I*255/5; % 定标取255个级别,使I/5最亮
image(B) % 以图案表示干涉条纹
colormap(gray(255)); % 用灰度级别显示图案
pause(1)
end