5双光束干涉的仿真

  • 格式:doc
  • 大小:26.50 KB
  • 文档页数:2

西安邮电大学光学课程设计
5.图

-0.4-0.200.20.4-1-0.500.51x 10-301234
-1.5
-1
-0.5
0
0.5
1
1.5
x 10
-3

程序
% µ¥É«¹âµÄË«·ì¸ÉÉæ
clear;
Lamd=500e-9;
d=2*0.001;
Z=1;
yMax=5*Lamd*Z/d;
xs=yMax;
Ny=100;
ys=linspace(-yMax,yMax,Ny);
B=zeros(100,100);% y·½Ïò·Ö³É100µã
for i=1:Ny % ¶ÔÆÁÉÏÈ«²¿µã½øÐÐÑ-»·¼ÆËã
% ¼ÆËãµÚÒ»ºÍµÚ¶þ¸ö¹âÔ´µ½ÆÁÉϸ÷µãµÄ¾àÀë
L1 = sqrt((ys(i)-d/2).^2 + Z^2 );
L2 = sqrt((ys(i)+d/2).^2 + Z^2 );
Phi = 2*pi*(L2-L1)/Lamd; % ´Ó¾àÀë²î¼ÆËãÏàλ²î
I(i,:) = 4*cos(Phi/2).^2; % ¼ÆËã¸Ãµã¹âÇ¿£¨ÉèÁ½Êø¹âÇ¿Ïàͬ)
% Èô¿¼ÂǹâÆ×µÄ·Çµ¥É«ÐÔ,°ÑǰÁ½¾ä¸ÄΪºóËľä
%Nl=11; dL=linspace(-0.1,0.1,Nl); % Éè¹âÆ×Ïà¶Ô¿í¶ÈÕý¸º10%,
%Lamd1=Lamd*(1+dL'); % ·Ö11¸ùÆ×Ïß,²¨³¤ÎªÒ»¸öÊý×é
%Phi1 = 2*pi*(L2-L1)./Lamd1; % ´Ó¾àÀë²î¼ÆËã¸÷²¨³¤µÄÏàλ²î
%I(i,:) = sum(4*cos(Phi1/2).^2)/Nl; % µþ¼Ó¸÷²¨³¤Ó°Ïì,¼ÆËã¹âÇ¿
end
% ÔÚÆÁÉÏ»-³öͼÏó
subplot(1,2,1);
imagesc(xs,ys,I);
subplot(1,2,2);
plot(I(:),ys)