clear;
lam=500e-9;%設(shè)定波長，以lam表示
N=2;%設(shè)置縫的個數(shù)
a=2e-4;%設(shè)定縫寬
D=5;%設(shè)定縫到屏的距離
b=5*a;%設(shè)定兩個縫之間的間隔
ym=2*lam*D/a;%設(shè)定y方向的范圍
xs=ym;%設(shè)定x方向的范圍
n=1001;
ys=linspace（-ymymn）;%用線性采樣法產(chǎn)生一個一維數(shù)組ys，n是此次采樣的總點數(shù)
??????????????????????%采樣的范圍從-ym到y(tǒng)m，采樣的數(shù)組命名為ys
??????????????????????%次數(shù)組裝的是屏幕上的采樣點的縱坐標
for?i=1:n?%對屏幕上的全部點進行循環(huá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%結(jié)束循環(huán)
NC=255;%確定使用的灰度等級為255級
Br=（B/max（B））*NC;%定標使最大光強對應(yīng)于最大的灰度級（白色）
subplot（121）??????%用subplot創(chuàng)建和控制多坐標軸，subplot（1，2，1）表示把單前
????????????????????%窗口對象分成1*2塊矩形區(qū)域并在第一塊區(qū)域創(chuàng)建一個新的坐標軸
????????????????????%這個新的坐標軸設(shè)定為當前坐標軸
image（xsysBr）;?????%然后用image繪圖函數(shù)創(chuàng)建圖像
colormap（gray（NC））;%用灰度級顏色圖設(shè)置圖和明暗
subplot（122）??????%畫出沿y向的光強變化曲線，先用subpot函數(shù)創(chuàng)建和控制多坐標軸
????????????????????%把窗口對象分成1*2塊矩形區(qū)域并在第2塊區(qū)域創(chuàng)建一個新的坐標軸
????????????????????%把這個坐標軸設(shè)定成當前的坐標軸
plot（B1ys）;?????????%然后繪制以（B1ys）為坐標軸相連的折線

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

?????文件????????643??2011-01-08?20:14??牛頓環(huán).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