激光谐振腔的模式计算研究.docxVIP

激光谐振腔模式研究的MATLAB实现光信1001班 刘吉祥 U201013222摘要：谐振腔内的模式计算是分析激光器输出光束质量的前提和基础。本文在matlab环境下，采用Fox_Li数值迭代法计算了条形腔、矩形腔、圆形腔、倾斜腔的自再现模的振幅分布和相位分布，并比较了腔形、菲涅尔数、初始光强分布、倾斜扰动等因素对最终模式的影响，具有一定的实际应用价值。原理说明设初始时刻在镜I上有某一个场分布，则当波在腔中经第一次渡越而到达镜II时，将在镜II上形成一个新的场分布，场经第二次渡越后又将在镜I上形成一个新的场分布。每次渡越时，波都将因为衍射损失一部分能量，并引起能量分布变化，如此重复下去……由于衍射主要是发生在镜的边缘附近，因此在传播过程中，镜边缘附近的场将衰落得更快，经多次衍射后所形成的场分布，其边缘振幅往往都很小（与中心处比较），具有这种特征的场分布受衍射的影响也将比较小。可以预期：在经过足够多次渡越之后，能形成这样一种稳态场：分布不再受衍射的影响，在腔内往返一次后能够“再现”出发时的场分布，即实现了模的“自再现”。光学中的惠更斯—菲涅尔原理是从理论上分析衍射问题的基础，该原理的严格数学表示是菲涅尔—基尔霍夫衍射积分。设已知空间任意曲面S上光波场地振幅和相位分布函数为,由它所要考察的空间任一点P处场分布为，二者之间有以下关系式：式中，为与连线的长度，θ为S面上点处的法线和上述连线之间的夹角，为S面上的面积元，k为波矢的模。本文采用Fox—Li数值迭代法实现了条形腔、矩形腔、圆形腔、倾斜腔的自再现模的形成。实现方案2.1条形腔条形腔是一种理想模型，即一个方向有限长，而另一个方向上无限延伸的腔形，故只在长度有限的那个方向上发生衍射现象，迭代公式为一维的菲涅尔—基尔霍夫衍射积分：将条形腔的左镜面S上沿着之间划分N-1等分，则有N个点，每个区间为。右边镜面P上每一点的求解都需将左边镜面上的点逐点计算一遍并相加，如此循环迭代下去，最终会达到稳态分布。2.2矩形腔在矩形腔中，与连线的长度可以表示为，经过计算与推导可知：矩形腔的计算不需考虑整个面上的点的影响，可以按照、两个方向分离变量为，其中的计算与条形腔相同。2.3 圆形腔圆形腔的迭代思想与矩形腔相同，只是划分与矩形腔不同。圆形腔是按照径向和角向划分，在极坐标(r,Φ)下完成数值迭代，但在最后显示的时候，需要将极坐标还原成笛卡尔坐标系。具体思路是：由极坐标和直角坐标的转换关系，X=r COSΦ，Y=r sinΦ，其中，r、Φ为极坐标参量。将X、Y用相应的极坐标参量代换并代入衍射方程，得 （4）为了分离变量，对圆形镜谐振腔，其场分布函数经常采用如下形式：（5）式中：p表示场分布在径向的变化；f表示场分布按方位角以不同的正弦或余弦方式变化。将式(5)代入式(4)，可得：式中，右边积分可以分离为Φ和r的积分，方括号内Φ的积分可以仿照圆形镜共焦腔来进行，利用积分关系式中，Jl为l阶第一类贝塞尔函数。再将式(7)代入式(6)，可以将方程(6)化简为只含径向的本征方程：数值迭代开始前需要给定初始的场分布尺，对TEM00模，设初始场分布为均匀平面波，将0≤r≤a等分为N个点，令R1L(r)=1，即镜面上各点振幅均为1。第q次迭代后，r1,r2,…rN各点本征值为2.4 倾斜腔严格的平行平面腔只是一种理想情况，实际情况下出现一定的不平行性是不可避免的，这里主要考察倾斜条形腔对自再现模的影响，如图3所示：图3 倾斜平行平面腔的示意图两个镜面相对其理想位置(即两镜面与其公共轴线严格垂直的位置) 沿相反方向偏离同样大小的微小角度β,在镜的边缘处与理想位置的偏离线度δ。在δ甚小的情况下,且只考虑腔的旁轴光线,镜面上两点的距离M1′M2′与理想情况下相应两点的距离M1M2之差为:，故，于是衍射积分方程变为：，类似于条形腔，可以计算出倾斜条形腔的自再现模。实验结果与分析激光谐振腔模式的各种分析方法的比较特征向量矩阵法，Fox—Li数值迭代法、厄米一高斯展开法、快速傅立叶变换法（FFT）、有限差分法（FDM）和有限元法（FEM）。特别是Fox—Li数值迭代方法，它是一种模式数值求解中普遍适用的一种方法，只要取样点足够多，它可以用来计算任何形状开腔的自再现模，而且还可以计算诸如平行平面腔中腔镜的倾斜、镜面的不平整性等对模的扰动。其缺点是在菲涅耳数F很大时，计算工作量很大。 特征向量法是对腔镜进行有限元单元划分的，构造光束传输矩阵，通过求解特征矩阵的特征向量，即可获得腔镜光场分布的振幅和相位分布。矩阵运算时间与矩阵维数有着近似平方的关系，二维衍射积分方程的传输矩阵的维数高，计算需要数小时甚至数天的时间，对于大菲涅耳系数的谐振腔计算难度更大。3.2菲涅尔数对结果的影响对于条形腔，菲涅尔数F定义为