资源描述：
《激光原理课程设计(自再现模的迭代法)》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、摘要激光器谐振腔内的模式计算是提高激光器输出光束质量和应用自适应光学系统校正腔内像差的前提和基础。在理论分析的基础上,着重采用数值迭代解法（Fox-Li方法即福克斯和厉鼎毅迭代法）计算平行平面腔（F-P腔）在初始光场三阶梯分布条件下,自再现模的光场振幅的分布。用数学软件MATLAB建模编写计算程序，计算结果表明,在经过260次渡越后,归一化的振幅分布实际上已不再发生变化,则已找到光腔的一个自再现模式或横模分布。关键词:数值迭代法；光场振幅分布；MATLAB数值模拟目录1引言12理论分析2.1原理22.2平行
2、平面腔迭代法的算法实现33MATLAB数值模拟53.1程序原代码53.2数值模拟结果63.3结果分析74总结8【参考文献】91引言谐振腔是激光器的主要构造之一，使激光通过增益物质，实现光的自激振荡。在激光器出光的过程中,谐振腔内存在许多扰动因素,如腔镜失调、增益介质不均匀、热效应、腔镜变形等,这些腔内扰动因素都会引起不同程度的腔内像差,带来光束质量的下降和光束能量的降低。目前采用自适应光学系统对腔内像差进行校正,定量得到腔内像差扰动对输出光束模式的影响,通过适当的控制算法对像差波前进行实时地校正,而腔内模式
3、计算是这一过程的前提。平行平面腔（又称为F—P腔）它由两块平行平面反射镜组成，在激光发展史上最先被采用。目前，在中等以上功率的固体激光器和气体激光器中仍常常采用它。其主要优点是，光束方向性极好、模体积较大、比较容易获得单横模振荡等。谐振腔的经典理论仅给出了部分简单腔型的模式解析解。对于平行平面腔以及在激光器的不断发展过程中所涌现的许多新型结构谐振腔通常是没有解析结果的,必须采用各种数值模拟方法进行求解。所以本文基于平行平面腔，研究初始光场三阶梯分布条件下,自再现模的光场振幅。由于平行平面腔振荡模所满足的自再
4、现积分方程：v(x,y)=γ至今尚得不到精确的解析解，因此本文致力于研究平面腔模的迭代解法（Fox-Li方法）。Fox-Li方法是一种模式数值求解中普遍适用的一种方法,只要取样点足够多,它原则上可以用来计算任何形状开腔的自再现模,并且,还可以计算诸如腔镜的倾斜、镜面的不平整性等因素对腔内模式造成的扰动.2理论分析2.1原理所谓迭代法：用变量的旧值不断递推出新值解决问题的方法，通常用于数值计算。对于开放式光腔，镜面上稳态场分布的形成可以看成是光在两个界面间往返传播的结果。因此，两个界面上的场必然是互相关联的：
5、一个镜面上的场可以视为由另一个镜面上的场所产生，于是求解镜面上稳态场的分布问题就归结为求解一个积分方程。本文基于初始光场三阶梯分布条件下，分析自再现模的光场振幅分布。利用迭代公式nθ（x'，y'）（x，y）U2U1菲涅尔—基尔霍夫积分公式中各量的意义(1.1)(1.2)其中，(1.2)式称为积分方程的核，可直接进行计算。首先，假设在某一个镜面上存在一个初始场分布，将它代入上式，计算在腔内经第一次渡越而在第二个镜面上生成的场，然后再将所得的代入(1.1)式，计算在腔内经第二次渡越而在第一个镜面上生成的场。如此
6、反复运算并注意经过足够多次以后，在腔面上能否形成一种稳态场分布。在对称开腔的情况下，当j足够大时，有数值计算得出的，，，能否满足下式关系，…式中，γ为复常数。如果直接数值计算得出了这种稳定的场分布，则可认为找到了腔的一个自再现模或横模。原理示意图如下：uq+1uqu1ⅠⅡLu1u3…u2u4…2au3uq+1u2uqu1u2u3初始入射波（a）（b）（c）开腔中自再现模的形成（a）理想开腔；（b）孔阑传输线；（c）自再现模的形成2.2平行平面腔迭代法的算法实现以对称条状腔为例，分析平行平面腔中自再现模的形成
7、。考察镜的宽度为2a，腔长为L的对称条状腔。矩形平面镜腔（x，y）ρ2aL（x'，y'）其中所以整理得：该条状腔的模式迭代方程应为这里以一列均匀平面波作为第一个镜面上的初始激发波。由于重要的只是振幅的相对分布，因此，可以取即认为整个镜面为等相位面，且镜面上各点波的振幅为1。代入迭代方程数值计算求出，然后将归一化，即取然后继续迭代运算。3MATLAB数值模拟3.1程序原代码clear,clcglobalstepsLkalamdalamda=input('波长lamda=');L=input('腔长L=');a
8、=input('镜长a=');N=input('渡越次数N=');k=2*pi/lamda;steps=500;x=linspace(-a,a,steps);u_=ones(1,steps);form=1:Nformm=1:stepsu0(mm)=QR(x(mm),u_);end;u_=u0/max(abs(u0));endsubplot(2,1,1)plot(x,abs(u0)/abs(u0(steps/2))
