%page 77
%等光程差面
%Equal optical path difference surface
lamda=630*1e-9;
d=0.5*1e-3;
m=3;
[x,y]=meshgrid([-10:0.15:10],[-10:0.15:10]);
u=(x.^2+y.^2)./( (d/2)^2-(m*lamda/2).^2 )+1;
z=(u.*(m*lamda/2).^2).^0.5;
mesh(x,y,z);
%colormap autumn
hold on
n=5;
[x1,y1]=meshgrid([-10:0.15:10],[-10:0.15:10]);
u1=(x1.^2+y1.^2)./( (d/2)^2-(n*lamda/2).^2 )+1;
z1=(u1.*(n*lamda/2).^2).^0.5;
mesh(x1,y1,z1);
hold on
q=1;
[x2,y2]=meshgrid([-10:0.15:10],[-10:0.15:10]);
u2=(x2.^2+y2.^2)./( (d/2)^2-(q*lamda/2).^2 )+1;
z2=(u2.*(q*lamda/2).^2).^0.5;
mesh(x2,y2,z2);
hold off
xlabel('x axis')
ylabel('y axis ')
title('Equal optical path difference surface')
基于matlab的等光程差面
最新推荐文章于 2023-08-30 19:06:28 发布