exn672.m

MATLAB及在电子信息课程中的应用.感觉很有用,本来也不花钱,有想用的来拿吧

%《MATLAB及其在理工课程中的应用指南》(十一五规划版)第六章例6－7－2程序
% 单色光的单缝衍射
% 西安电子科技大学出版社出版  陈怀琛编著 2007年7月
%
clear; % 清工作空间，输入波长，缝宽，和屏幕距离
Lambda = input('输入光的波长(单位为 nm): (书上取500) ');
Lambda = Lambda * 1e-9;   		% 将nm换 变为 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;          % 屏幕上的点数（沿y向）
ys = linspace(-ymax,ymax,Ny);
NPoints = 51;     % 缝上的点数（沿y向）
yPoint = linspace(-aWidth/2,aWidth/2,NPoints);   % 把缝上的点数设成数组
for j=1:Ny			% 对屏幕上y向各点作循环
% 对光缝中各点作循环,计算缝点到屏幕位置的距离
   L = sqrt((ys(j)-yPoint).^2 + Z^2 );		% L是一个数组
   Phi = 2*pi.*(L-Z)./Lambda;   % 计算相对于屏幕中心的相位差,也是一个数组
   %@ 求每个分量的累加和
   SumCos = sum(cos(Phi)); 	% 数组求和
   SumSin = sum(sin(Phi)); 
 % 求屏幕上的归一化光强;
 B(j) = (SumCos^2 + SumSin^2)/NPoints^2;  
 end
plot(ys,B,'*',ys,B);grid	% 屏幕上光强与位置的关系曲线
axis([-ymax, ymax, 0.0, 1.0]);
%xlabel('屏幕上的位置 (m)'); 
%ylabel('屏上的归一化光强');
%title('归一化的屏上光强');
set(gcf,'color','w')