激光腔中模式模拟.docVIP

激光腔中模式模拟 The simulation of the mode in the resonant cavity of laser 光电子科学与工程学院0304班 张欢 摘要：对于激光腔通常没有解析解而只能有数值解，其研究方法有包括Fox_li方法在内的很多方法。此处是在Fox_li方法的基础上采用一种矩阵相乘的方法实现从一个面到另一个面的迭代的。同时还利用孔阑矩阵和二次位相弯曲矩阵实现了圆形腔和共焦腔的模场研究。 summary: The cavity of the laser usually has numerical value results but no analysis results, and there are many methods including the method of Fox-li to discuss about the problem .in this passage, I use a matrix-multiplying method which is based on the method of Fox-li to realise the alternation from one face to the other. Meanwhile, I use the matrix of aperture diaphragm and the secondary bending of phrase to realise the research on the round-cavity and the confocal cavity. 关键字： 平平腔 条形腔 共焦腔 自再现模式 FFT变换法 矩阵相乘Fox_li迭代方法 引言：对于激光腔的研究方法最常见的就是几何法和衍射方程法。在谐振腔菲涅尔数（）比较大的情况下考虑近轴光，镜面的衍射效应对模场的影响很小，因而可以采用几何法对谐振腔进行研究，然而当谐振腔的菲涅尔数N越小，则衍射效应的影响就越大，因而几何光学理论描述的局限性也表现得越明显。在N值足够小时，由真正的衍射效应所造成的光能溢出损耗变得不能忽略。而腔内波形场的振幅和相位分布可能远远偏离于 几何光学所预计的简单情况，此时必须采用衍射理论描述。 说明： 如下图：对于从一个面到第二个面的光场 由菲涅尔衍射其迭代公式为： 在L〉〉a〉〉的前提下可化为： 这样迭代公式化为： 若采用Fox_li迭代方法，则对于矩形腔可以将场分布利用变量分离分解为X向和Y向的场分别迭代然后再相乘即： 这个方法固然可以，但是在用matlab编程时是完全不用分解成两个分量的。这是因为matlab是基于矩阵的运算，而注意到矩阵运算的特点：对于两个矩阵相乘：其中均为的矩阵。若用表示A矩阵的第m行n列的元素，则的值就是的m行与CH的n列对应元素相乘之后的和，即：，将该式与迭代方程相比较可知：若将CH矩阵表示为的一系列值那么迭代方程（一维）就等价于：。对于二维的情况。易证明：迭代方程可表示为：下面用一个简单的示意图来说明这种算法：假设一个腔镜为的矩阵，矩阵中的每一个值包含了腔镜上对应点相对复振幅的信息。，为简单起见先假设将腔面划分为以后各点的间距为1（单位）则构造的CH矩阵（称为相位转换矩阵）为： 则上述关系为： 有矩阵运算法则知： 这正是二维迭代方程：的离散表达式， 由此可见：在用迭代法研究激光枪的模特性时用矩阵的相乘就可以很好的完成了，这正是本程序的核心算法。 在这个基础上，我们不难想到，对于圆形腔镜的平平腔由于每一次迭代的U都是一个呈圆分布的场，那么就可以在矩形平平腔的基础上乘上一个圆形孔阑矩阵将圆孔以外的场函数置为0再进行下一次迭代。 对于球面腔的情形，由信息光学的知识我们知道：球面的反射镜就相当于一凸透镜，而透镜对场的影响就想当于引入了一个二次位相弯曲（其中f为焦距）。于是若在每次迭代后在场函数上乘以一个二次位相弯曲矩阵，那么平平腔就顺理成章的成为了球面腔。 以上三种腔所提到的算法正是本程序所用到的算法 有于程序代码太长，在此处仅给出核心矩阵的算法： CH（转换矩阵）： CH=ones(N,N); for m=1:N CH(m,:)=(m-1)*D/N:(-1)*D/N:(m-N)*D/N; end CH=CH.^2; CH=exp(-i*k*CH/2/L); CH=sqrt(i*exp(-i*k*L)/lamda/L)*CH; CP矩阵（球面腔二次位相弯曲）： CP=ones(N,N); CPX=-(N-1)*D/2/N:1*D/N:(N-1)*D/2/N; CPX=exp(-i*k*CPX.^2/L); for n=1:N CP(n,:)=CPX*exp(-i*k*(((n-(N+1)