clc????????????????%%%?《信息光學數字實驗室?matlab版》-錢曉凡編著??實驗二??二維光場分析（球面波）
clear
close?all
lamda=6328e-10;??????????????????????%波長，單位：米
k=2*pi/lamda;????????????????????????%波數
x0=0.001?????????????????????????????%點光源的x坐標，單位：米
y0=0.001;????????????????????????????%點光源的y坐標，單位：米
z=0.3????????????????????????????????%觀察面到點光源的垂直距離，單位：米
L=0.005??????????????????????????????%觀察面的尺寸，單位：米
x=linspace（-L/2L/2512）;y=x;????????%構建x坐標和y坐標
[xy]=meshgrid（xy）;?????????????????%構建二維坐標網格
U1=exp（j*k*z）.*exp（j*k.*（（x-x0）.^2+（y-y0）.^2）/2/z）;??%發散球面光波
ph1=k.*（（x-x0）.^2+（y-y0）.^2）/2/z;????%發散球面波的實際相位
figuresurfl（ph1）shading?interpcolormap（gray）??????%%%?圖1
phyp1=angle（U1）;?????????????????????%發散球面波的包裹相位?????-pi到pi
figureimshow（phyp1[]）??????????????????????????????%%%??2
U2=exp（-j*k*z）.*exp（-j*k.*（（x-x0）.^2+（y-y0）.^2）/2/z）;?%會聚球面光波
ph2=-k.*（（x-x0）.^2+（y-y0）.^2）/2/z;???%會聚球面波的實際相位
figuresurfl（ph2）shading?interpcolormap（gray）???????%%??3
phyp2=angle（U2）;?????????????????????%會聚球面波的包裹相位
figureimshow（phyp2[]）?????????????????????????????????%%%４
figure?plot（ph2（257:）‘--‘）????????%實際相位的剖線　　　圖５
hold?on??????????????????????????????%保持當前圖像
plot（phyp2（257:）‘r‘）???????????????%包裹相位的剖線??????圖5
diff1=U1+1;??????????????????????????%觀察面上發散球面光與垂直照射平行光的干涉???
I1=diff1.*conj（diff1）;???????????????%觀察面上的光強????????conj（）取復共軛
figureimshow（I1[0max（max（I1））]）???%??圖6
diff2=U2+1;??????????????????????????%觀察面上會聚球面光與垂直照射平行光的干涉
I2=diff2.*conj（diff2）;??%觀察面上的光強?????或者?I2=abs（diff2）.*abs（diff2）;
figureimshow（I2[0max（max（I2））]）???????%??7

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

?????文件??????28781??2016-04-17?11:44??2D?light?field?analysis\1.jpg

?????文件?????115639??2016-04-17?11:44??2D?light?field?analysis\2.jpg

?????文件??????24362??2016-04-17?11:43??2D?light?field?analysis\3.jpg

?????文件?????115667??2016-04-17?11:43??2D?light?field?analysis\4.jpg

?????文件??????16824??2016-04-17?11:43??2D?light?field?analysis\5.jpg

?????文件??????75137??2016-04-17?11:43??2D?light?field?analysis\6.jpg

?????文件???????1810??2016-04-17?11:38??2D?light?field?analysis\CX2_1.m

?????文件??????75137??2016-04-17?11:43??2D?light?field?analysis\二維光場分析（球面波）.jpg

?????目錄??????????0??2016-04-17?11:46??2D?light?field?analysis

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

???????????????453357????????????????????9