模拟杨氏双缝干涉实验的M程序_用MATLAB进行仿真_杨氏干涉_matlab

% 杨氏干涉的模拟 close all figure('position',[217 266 694 244]); % Author's email:[email protected] set(gcf,'DoubleBuffer','on'); axes('position',[0.12,0.08,0.6,0.8]);hold on; rectangle('position',[-1,1.1,0.1,3],'FaceColor',[0.1,0.3,0.4]); rectangle('position',[-1,-4,0.1,3],'FaceColor',[0.1,0.3,0.4]); rectangle('position',[-1,-0.8,0.1,1.6],'FaceColor',[0.1,0.3,0.4]); axis([-4,12,-5,5]); ha=plot([-3.8,-3.8],[0.9,0.9],'r','linewidth',3); hb=plot([-3.8,-3.8],[-0.9,-0.9],'r','linewidth',3); title('Double slit interference \copyright zjliu') L=-3.8; rectangle('position',[11,-4,0.1,8],'FaceColor',[0,0,0]); plot([-1,4.1],[-3.5,-3.5],'k'); plot([5.9,11],[-3.5,-3.5],'k'); plot([-0.7,-1,-0.7],[-3.3,-3.5,-3.7],'k'); plot([10.7,11,10.7],[-3.3,-3.5,-3.7],'k'); text(4.8,-3.5,'3m'); plot([-1.5,-1.5],[0.9,0.35],'k'); plot([-1.5,-1.5],[-0.9,-0.35],'k'); plot([-1.65,-1.5,-1.3],[0.6,0.9,0.6],'k'); plot([-1.6,-1.5,-1.3],-[0.6,0.9,0.6],'k'); text(-2.1,0,'2mm'); while L<=-0.95; L=L+0.05; set(ha,'xdata',[-3.8,L]); set(hb,'xdata',[-3.8,L]); pause(0.05); end xas=-0.95;yas=0.9; xbs=-0.95;ybs=-0.9; po=[-3.5:3.5]; ka=(po-yas)*i+(11+0.95);ka=ka./abs(ka); kb=(po-ybs)*i+(11+0.95);kb=kb./abs(kb); for n=1:8; ah(n)=plot(xas,yas,'r'); bh(n)=plot(xbs,ybs,'r'); end r=0; za=xas+yas*i;zb=xbs+ybs*i; Le=sqrt(11.95^2+3.5^2); while r<=Le; for n=1:8; Ta=za+r*ka(n); Ya=imag(Ta); if abs(imag(Ta))>=abs(po(n)); Ya=po(n); end set(ah(n),'xdata',[xas,min(11,real(Ta))],'ydata',[yas,Ya]); Tb=zb+r*kb(n); Yb=imag(Tb); if abs(imag(Tb))>=abs(po(n)); Yb=po(n); end set(bh(n),'xdata',[xbs,min(11,real(Tb))],'ydata',[ybs,Yb]); end r=r+0.05; pause(0.05); end s=meshgrid(linspace(4,-4,300))'; De=abs((s-1)*i+3000)-abs((s+1)*i+3000); lambda=0.6328e-3; de_A=De/lambda*pi*2; It=1-cos(de_A); axes('position',[0.72,0.16,0.26,0.64]); cc=cat(3,It/2,zeros(size(It)),zeros(size(It))); imshow(cc,[]); xlabel('\copyright zjliu')