基于MATLAB的正三角形夫瑯禾費(fèi)衍射現(xiàn)象

基于MATLAB的正三角形夫瑯禾費(fèi)衍射現(xiàn)象夫瑯禾費(fèi)衍射振幅公式：

觀察屏上的光強(qiáng)表達(dá)式為：

其中指開孔平面上光的分布，其一般是均勻的，故通常為常數(shù)，這里用A表示。理論基礎(chǔ)單縫衍射模擬單縫衍射計(jì)算機(jī)模擬代碼clc;clear;a=-2*pi:0.0001*pi:2*pi;p1=(1-sinc(a)).^2;%方便下面著色p2=sinc(a).^2;figure;plot(a,p2);xlabel('kasinθ');ylabel('光強(qiáng)I/I0');title('單縫衍射強(qiáng)度分布');lgray=zeros(256,3);fori=0:255lgray(i+1,:)=(255-i)/255;endfigure;imagesc(p1)title('單縫衍射模擬圖');colormap(lgray)S*單縫衍射實(shí)驗(yàn)裝置屏幕單縫的夫瑯和費(fèi)衍射線光源在透鏡L1的物方焦平面接收屏在L2象方焦平面多縫衍射多縫衍射縫數(shù)2多縫衍射代碼clc;clear;n=input('多縫衍射縫數(shù)');a=-2*n*pi:0.0001*pi:2*pi*n;p1=1-(sin(n*a)./sin(a)).^2;%方便下面著色p2=(sin(n*a)./sin(a)).^2;figure;plot(a,p2);xlabel('kasinθ');ylabel('光強(qiáng)I/I0');title('多縫縫衍射強(qiáng)度分布');lgray=zeros(256,3);fori=0:255lgray(i+1,:)=(255-i)/255;endfigure;imagesc(p1)title('多縫衍射模擬圖');colormap(lgray);θLP0P多縫的夫瑯和費(fèi)衍射多縫衍射縫數(shù)3多縫衍射縫數(shù)3光強(qiáng)分布正三角形孔衍射代碼clearalla=0.00004;lmda=500e-9;f=10;H=a*sin(pi/3);x=-1:0.005:1;y=-1:0.005:1;fori=1:1:401forj=1:1:401X(i)=2*pi*x(i)/(tan(pi/3)*lmda*f);Y(j)=2*pi*y(j)/(lmda*f);I1(i,j)=(sin(0.5*H*(Y(j)-X(i))))^2/((Y(j)-X(i))^2+eps);I2(i,j)=(sin(0.5*H*(Y(j)+X(i))))^2/((Y(j)+X(i))^2+eps);I3(i,j)=2*cos(H*X(i))*(sin(0.5*H*(Y(j)-X(i))))*(sin(0.5*H*(Y(j)+X(i))))/(Y(j)*Y(j)-X(i)*X(i)+eps);I(i,j)=(4/(3*X(i)*X(i)+eps))*(I1(i,j)+I2(i,j)-I3(i,j));endendm=max(I(:));n=min(I(:));I0=(I-n)/(m-n);figure(1)imshow(I0)figure(