pijian.m

光学模拟劈尖干涉的程序

% 关于光的等厚干涉(劈尖)的动画
close all;clc;clear;
set(gcf,'DoubleBuffer','on'); 
axes('position',[0.12,0.08,0.6,0.8]);hold on;
% Author: zjliu
% Author's email: [e-mail]
title('关于光的等厚干涉(劈尖)的动画','fontsize',16) 
axis manual;
rectangle('position',[0.15,0.2,0.7,0.01],'FaceColor',[0.1,0.3,0.4]);
plot([0.2,0.2],[0.2,0.4]);
plot([0.8,0.2],[0.2,0.4]);
style='none';
for k=1:10;
    h(k)=plot([0.21+(k-1)*0.06]*[1,1],[0.6,0.6],'r',...
        'linewidth',1,'EraseMode',style);
end
G=0.6;
while G>=0.2; 
   G=G-0.01; 
   set(h,'ydata',[0.6,G]); 
   pause(0.1); 
end
for k=1:10;
    hr(k)=plot([0.21+(k-1)*0.06]*[1,1],[0.2,0.2],'r',...
        'linewidth',2);
end
while G<0.4; 
   G=G+0.01; 
   set(hr,'ydata',[0.2,G]); 
   pause(0.1); 
end
z=[0.8,0.2]+i*[0.2,0.4];
N=400;
An=angle(diff(z))-pi/2;
An=exp(i*An);
zz=linspace(z(1),z(2),N);
Si=abs(sin(linspace(0,pi*20,N))).^2/20;
zz=zz+Si*An;
plot(zz,'b');