基于matlab的等光程差面

%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')