平行平面腔腔镜光场分布数值模拟.ppt

平行平面腔腔镜光场分布数值模拟 基本原理：菲涅耳-基尔霍夫衍射积分 基本方法：微元分析法(以矩形镜为例) 平行平面腔腔镜光场分布数值模拟 如图由M1及M2组成的矩形镜平行平面腔; 已知M1镜的初始场分布u1; 求经过腔长L0渡越后到达镜M2的场分布u2？ 平行平面腔腔镜光场分布数值模拟 微元分析法：将矩形镜分为M?N均匀网格， 当M和N足够大时，每一网格都可以看作一 个点光源，其传播r后的场可表为： 因此，镜M2上某点(网格)处的光场等于镜M1 上所有网格点光源传播到此处的场叠加。 平行平面腔腔镜光场分布数值模拟 网格划分： 平行平面腔腔镜光场分布数值模拟 镜M1上的任意点(i,j)的场u1(i,j)传播到镜 M2上的点(u,v)时的场为du2(u,v),则： 程序设计：界面 采用Delphi(Pascal) 编程软件：界面如图 所示(共13个对象)。 将初始光场分布分为 4个矩形区域。 程序设计：一次渡越后的结果 上图为强度-灰度图 下图为强度-条形图 M1镜初始场分布均匀 程序设计：基尔霍夫衍射积分 for u:=0 to M-1 do //基尔霍夫衍射积分 for v:=0 to N-1 do for i:=0 to M-1 do for j:=0 to N-1 do begin r:=sqrt((j-v)*(j-v)*dx2+(i-u)*(i-u)*dy2+L02); cosa:=L0/r; Cr:=(1+cosa)/r; ur:=u1[i][j][0]; ui:=u1[i][j][1]; t:=k*r; coskr:=Cos(t); sinkr:=Sin(t); u2[u][v][0]:=u2[u][v][0]+(ur*coskr+ui*sinkr)*Cr; u2[u][v][1]:=u2[u][v][1]+(ui*coskr-ur*sinkr)*Cr; end; 程序设计：完整的代码 数据定义： 基本常数赋初值： 光强分布初始化操作： 基尔霍夫衍射积分及图形绘制： * * 专题讨论 专题讨论 M1 M2 L0 专题讨论 专题讨论 M1 M2 M格 N格 i行 j列 M格 N格 u行 v列 专题讨论 专题讨论 专题讨论 专题讨论 专题讨论 *