目的是要利用MATLAB模擬 做出一矩形狹縫配合一透鏡於屏幕上觀測其繞射分佈
而小弟於網路上有找到單縫繞射的MATLAB模擬
其程式如下
lamda=500e-9; %波長
N=1;%縫數
a=2e-4;D=5;d=5*a
ym=2*lamda*D/a;xs=ym;%屏幕上Y的範圍
n=1001;%屏幕上的點數
ys=linspace(-ym,ym,n);%定義區域
for i=1:n
sinphi=ys(i)/D;
alpha=pi*a*sinphi/lamda;
beta=pi*d*sinphi/lamda;
B(i,:)=(sin(alpha)./alpha).^2.*(sin(N*beta)./sin(beta)).^2;
B1=B/max(B);
end
NC=256;%確定灰度的等級
Br=(B/max(B))*NC;
subplot(1,2,1)
image(xs,ys,Br);
colormap(hot(NC));色調處理
subplot(1,2,2)
plot(B1,ys,'k');

手邊沒電腦，不能幫你跑模擬,我的猜測如下。第一，這是解析解，並不是模擬;第二，兩個sinc應為繞射項加上干涉項，所以a應為狹縫寬度，D為傳播距離，d為干涉參數，5代表slit有五個。