二维光场分析球面波。MATLAB代码。

二维光场分析（球面波）。编程实现单色球面波在给定二维平面上的构建，并用相位表达式直接计算其相位，与通过复光场得到的相位相比较，观察相位包裹（缠绕）现象。与平行光相干涉，观察干涉光强分布。

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