单缝衍射

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('光强I/I0');   
title('单缝衍射强度分布'); 
lgray=zeros(256,3);  
for i=0:255   
lgray(i+1,:)=(255-i)/255;   end   
figure;   
imagesc(p1)   
title('单缝衍射模拟图');  colormap(lgray) ;