激光原理课程设计--平行平面腔自再现模fox-li数值迭代解法及matlab实现

《激光原理课程设计--平行平面腔自再现模fox-li数值迭代解法及matlab实现》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、华南理工大学广州学院电气工程系生产实习报告激光原理课程设计题目：方形镜平行平面腔自再现模Fox-Li数值迭代解法及MATLAB实现院系理学院专业班级学生姓名指导教师提交日期2012年4月15日16华南理工大学广州学院电气工程系生产实习报告目录一、设计目的1二、设计要求和设计指标1三、设计内容13.1Fox-Li平行平面腔的迭代解法13.2matlab实现33.2.1迭代解法的过程33.2.2程序实现43.2.3自再现模形成的判断63.3GUI界面的制作8四、本设计改进建议9五、设计感想9六、主要参考文献9附录16激光

2、原理课程设计一、设计目的为了加深对激光原理中Fox-Li平行平面腔的迭代解法的理解，学习matlab的使用，锻炼运用数值方法解决专业问题的能力。二、设计要求和设计指标在matlab用Fox-Li平行平面腔的迭代解法求得激光器腔镜面上的光场分布。求出距离镜面中点为x处的光场的振幅A和相位P，并作出二维图像。在得出的一维图像的基础上作出镜面上光强的二维分布。三、设计内容3.1Fox-Li平行平面腔的迭代解法谐振腔是激光器必备条件之一，它使激光反复通过增益物质，从而实现光的自激振荡。在激光的发展史上最早提出的是平行平面腔，

3、又称为F—P腔，它由两块平行平面反射镜组成，第一台红宝石激光器的谐振腔就是用它来做成的。对于开放式光腔，镜面上稳态场分布的形成可以看成是光在两个界面间往返传播的结果。因此，两个界面上的场必然是互相关联的：一个镜面上的场可以视为由另一个镜面上的场所产生，于是求解镜面上稳态场的分布问题就归结为求解一个积分方程。考虑在开腔中往返传播的一列波。设初始时刻在镜I上有某一个场分布，则当波在腔中经第一次渡越而到达镜II时，将在镜II上形成一个新的场分布，场经第二次渡越后又将在镜I上形成一个新的场分布。每次渡越时，波都将因为衍射损失

4、一部分能量，并引起能量分布变化，如此重复下去……由于衍射主要是发生在镜的边缘附近，因此在传播过程中，镜边缘附近的场将衰落得更快，经多次衍射后所形成的场分布，其边缘振幅往往都很小（与中心处比较），具有这种特征的场分布受衍射的影响也将比较小。可以预期：在经过足够多次渡越之后，能形成这样一种稳态场：分布不再受衍射的影响，在腔内往返一次后能够“再现”出发时的场分布，即实现了模的“自再现”，具体过程图1所示：16激光原理课程设计图1开腔中自再现模的形成光学中的惠更斯—菲涅尔原理是从理论上分析衍射问题的基础，该原理的严格数学表示

5、是菲涅尔—基尔霍夫衍射积分。设已知空间任意曲面S上光波场地振幅和相位分布函数为,由它所要考察的空间任一点P处场分布为，二者之间有以下关系式：式中，为与连线的长度，θ为S面上点处的法线和上述连线之间的夹角，为S面上的面积元，k为波矢的模。而对于方形镜平行平面镜将按，的幂级数展开，当满足和时16激光原理课程设计从而得到将上式分离变量。令得到方形镜中。3.2matlab实现3.2.1迭代解法的过程本文采用Fox—Li数值迭代法得到了了镜面上自再现模在x方向的分布并推广到整个镜面，最终动态显示每次渡越镜面上光场分布。虽然是复

6、数积分，但其和实数积分实现方法相同，即取一定步长，用矩形面积的和代替函数的定积分。16激光原理课程设计下面是程序框图：3.2.2程序实现源程序：clear,clcglobalstepsLkalamda=input('波长lamda=');L=input('腔长L=');a=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

7、(mm)=QU(x(mm),u_);end;u_=u0/max(abs(u0));16激光原理课程设计endsubplot(2,1,1)plot(x,abs(u0)/abs(u0(steps/2)))xlabel('x');ylabel('相对振幅');angle_u0=angle(u0)/pi*180;angle_u0=angle_u0-angle_u0(steps/2);subplot(2,1,2)plot(x,angle_u0)xlabel('x');ylabel('相对相位');functiony=QU(x,

8、u)globalstepsLkax_=linspace(-a,a,steps);step_length=2*a/(steps-1);y=sqrt(1i/L*exp(-1i*k*L))*sum(exp(-1i*k/2/L*(-x_+x).^2).*u)*step_length;在菲涅耳数为6.25，渡越次数为1的时候做出下图：16激光原理课程设计