
clear all;
%默认光强为1
lambda = 6*(1e-7);
theta = 0.005; %夹角
length = 0.001; %劈尖长度
x1 = length*cos(theta);
[x,y] = meshgrid(0:0.00001:x1,0:0.00001:0.001);
phi = 2*pi*(2*tan(theta)*x/lambda+1/2);
I = 4*cos(phi/2).^2; %干涉光强
z = tan(theta)*x;
surf(x,y,z,I)
shading interp;
colormap gray;
xlabel('x/m');
ylabel('y/m');
zlabel('z/m');
title(['劈尖夹角为:',num2str(theta/pi*180,4),'°']);
运行结果:
参考:屈玉福,陈沛戎 MATLAB光学仿真实用教程