牛顿环
其实还是比较简单的
看我发的
function niudunhuan(R,lamda)
% 演示牛顿环干涉实验
% R为牛顿环的半径
% lamda为光波长
% 光强分布图
% 要求R与lamda的比值不能太大，否则效果不好
clc
t=R/500;
x=-R:t:R;y=x;
[x,y]=meshgrid(x,y);
r=sqrt(x.^2+y.^2);
d=R-sqrt(R.^2-r.^2);
e=sin(4*pi*d./lamda);
E=sin(4*pi*d./lamda).^2;
imwrite(E,'e.jpg')
subplot(121)
imshow('e.jpg')
subplot(122)
plot(r,e)
axis([0 R/1.25 -1 1])
保存到work里后，直接运行就可以了
|