clear;
lam=500e-9;%設定波長，以lam表示
N=2;%設置縫的個數
a=2e-4;%設定縫寬
D=5;%設定縫到屏的距離
b=5*a;%設定兩個縫之間的間隔
ym=2*lam*D/a;%設定y方向的范圍
xs=ym;%設定x方向的范圍
n=1001;
ys=linspace（-ymymn）;%用線性采樣法產生一個一維數組ys，n是此次采樣的總點數
??????????????????????%采樣的范圍從-ym到ym，采樣的數組命名為ys
??????????????????????%次數組裝的是屏幕上的采樣點的縱坐標
for?i=1:n?%對屏幕上的全部點進行循環計算，則要進行n次計算??
????sinphi=ys（i）/D;%利用光柵衍射的理論公式計算衍射角的正弦值
????alpha=pi*a*sinphi/lam;%利用光柵衍射的理論公式分別計算單縫衍射因子和
????beta=pi*b*sinphi/lam;%縫間干涉因子
????B（i:）=（sin（alpha）./alpha）.^2.*（sin（N*beta）./sin（beta））.^2;%計算每一點的光強
????B1=B/max（B）;%計算相對光強
end%結束循環
NC=255;%確定使用的灰度等級為255級
Br=（B/max（B））*NC;%定標使最大光強對應于最大的灰度級（白色）
subplot（121）??????%用subplot創建和控制多坐標軸，subplot（1，2，1）表示把單前
????????????????????%窗口對象分成1*2塊矩形區域并在第一塊區域創建一個新的坐標軸
????????????????????%這個新的坐標軸設定為當前坐標軸
image（xsysBr）;?????%然后用image繪圖函數創建圖像
colormap（gray（NC））;%用灰度級顏色圖設置圖和明暗
subplot（122）??????%畫出沿y向的光強變化曲線，先用subpot函數創建和控制多坐標軸
????????????????????%把窗口對象分成1*2塊矩形區域并在第2塊區域創建一個新的坐標軸
????????????????????%把這個坐標軸設定成當前的坐標軸
plot（B1ys）;?????????%然后繪制以（B1ys）為坐標軸相連的折線

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????643??2011-01-08?20:14??牛頓環.m

?????文件???????1695??2011-01-08?20:14??楊氏雙縫干涉.m

?????文件???????1347??2011-01-08?20:14??圓孔衍射.m

?????文件????????360??2011-01-08?20:14??正弦光柵.m

?????文件????????924??2011-01-08?20:14??單縫衍射.m

?????文件???????1437??2011-01-08?20:14??光柵衍射.m

?????文件????????993??2011-01-08?20:14??矩孔衍射.m

?????文件????????581??2011-01-08?20:14??朗伯余弦體.m

?????文件????????516??2011-01-08?20:14??邁克爾遜干涉.m

-----------?---------??----------?-----??----

?????????????????8496????????????????????9