%《MATLAB及其在理工課程中的應(yīng)用指南》(十一五規(guī)劃版)第六章例6－7－2程序
% 單色光的單縫衍射
% 西安電子科技大學(xué)出版社出版  陳懷琛編著 2007年7月
%
clear; % 清工作空間，輸入波長(zhǎng)，縫寬，和屏幕距離
Lambda = input('輸入光的波長(zhǎng)(單位為 nm): (書上取500) ');
Lambda = Lambda * 1e-9;   		% 將nm換 變?yōu)?m
aWidth = input('輸入縫寬 (單位為 mm): (書上取0.2或1或2) ');
aWidth = aWidth  * 0.001;  	% 將mm 變換為 m
Z = input('輸入縫到屏的距離 (單位為 m): (書上取1) ');
ymax = 3*Lambda*Z/aWidth;	%	屏幕范圍（沿y向）
Ny = 51;          % 屏幕上的點(diǎn)數(shù)（沿y向）
ys = linspace(-ymax,ymax,Ny);
NPoints = 51;     % 縫上的點(diǎn)數(shù)（沿y向）
yPoint = linspace(-aWidth/2,aWidth/2,NPoints);   % 把縫上的點(diǎn)數(shù)設(shè)成數(shù)組
for j=1:Ny			% 對(duì)屏幕上y向各點(diǎn)作循環(huán)
% 對(duì)光縫中各點(diǎn)作循環(huán),計(jì)算縫點(diǎn)到屏幕位置的距離
   L = sqrt((ys(j)-yPoint).^2 + Z^2 );		% L是一個(gè)數(shù)組
   Phi = 2*pi.*(L-Z)./Lambda;   % 計(jì)算相對(duì)于屏幕中心的相位差,也是一個(gè)數(shù)組
   %@ 求每個(gè)分量的累加和
   SumCos = sum(cos(Phi)); 	% 數(shù)組求和
   SumSin = sum(sin(Phi)); 
 % 求屏幕上的歸一化光強(qiáng);
 B(j) = (SumCos^2 + SumSin^2)/NPoints^2;  
 end
plot(ys,B,'*',ys,B);grid	% 屏幕上光強(qiáng)與位置的關(guān)系曲線
axis([-ymax, ymax, 0.0, 1.0]);
%xlabel('屏幕上的位置 (m)'); 
%ylabel('屏上的歸一化光強(qiáng)');
%title('歸一化的屏上光強(qiáng)');
set(gcf,'color','w')