天津大学科学计算选讲第二次作业激光被动调Q速率方程组数值仿真学.docx

天津大学 科学计算选讲 第二次作业 激光被动调 Q速率方程组数值仿真 学生 学生: 徐世林 学号: 1011202008 学院: 精密仪器与光电子工程学院 专业: 光电子技术 ( 博士) 日期: 2012 年 5 月 22 日 PAGE PAGE 1 激光被动调 Q速率方程组数值仿真 物理背景 Q值是评定激光器中光学谐振腔质量好坏的指标 , 称为品质因数。Q值定义为在激光谐振腔内，储存的总能量与腔内单位时间损耗的能量之比 , 数学表达式为 W2 v0dW dtQ 。式中 W 2 v0 dW dt 即单位时间内损耗的能量， v0 为激光的中心频率。 Q值是评定激光器中光学谐振 腔质量好坏的指标，在激光谐振腔内储能不变情况下，损耗越小 Q值越大。 一般固体脉冲激光器的直接输出时， 由于激光器工作阈值附近， 输出脉冲是 系列尖峰振荡， 输出脉宽比较宽， 输出功率低。 激光调 Q是一种广泛用于产生巨脉冲功率激光的运转方式， 通过改变 Q值( 改变阈值 ) ，控制激光产生时间， 输出巨脉冲， 压缩脉冲宽度， 提高峰值功率。 激光调节 Q值的途径一般是采取改变腔内损耗的办法来调节腔内的 Q值。调 Q激光器输出的脉冲从时间上分为三阶段： 脉冲建立时间、 脉冲前沿时间和脉冲后沿时间。 从 Q值最小变到最大 Q值即损耗从最大变到最小需要的时间叫开关时间， 开关时间对激光脉冲的影响很大， 按开关时间的大小分为快、慢两种类型。 被动调 Q是利用材料一种可饱和吸收特性做为 Q开关，这种可饱和吸收染料是一种非线性吸收物质， 把它放在谐振腔内， 利用它对光的可饱和吸收特性来改变谐振腔内的吸收损耗， 起到 Q开关的作用。 可饱和吸收特性指材料的透射特性随着能量密度的增大而变透明，当能量密度达到某值时，材料的吸收达到饱和， 从而具有很高的透射率。吸收体对入射光强的改变应是这两个过程的综合过程。 高能态粒子越多，则吸收光强越少。当高能态粒子数 n2 和基态粒子数 n1 相等时， 这时从低能态到高能态的受激跃迁的粒子数和从高能态到基态放出光子的粒子 数相等， 即吸收的和放出的光子相等， 达到饱和。 常用的可饱和吸收材料有机染料，掺吸收性离子或含有色心（由点缺陷引起）的晶体。 激光器速率方程组是表征激光器腔内光子数和工作物质各有关能级上的粒 子数随时间变化的微分方程组。 对于一般的脉冲激光器， 脉冲形成时间长， 泵浦、自发辐射、受激跃迁过程都是存在的。采用 Q开关技术后，各参量之间关系发生 了很大变化，需要根据 Q开关的过程分析速率方程变化。 由于被动调 Q 脉冲通常 只有纳秒量级， 在脉冲期间可以不考虑 抽运源的影响， 研究连续脉冲的情况， 需要考虑抽运速 率、增益介质的上能级 寿命、可饱和吸收 体的恢复时间 等因素，得到下面考虑了可饱和吸收体激发态吸收的连续抽运被动调 Q 速率方程组为 d dt t r 2 nl 2 gsn gsl s 2 es(n 0s n gs )ls 1 (ln ) R dn n n Rp(1 ) c n （1） dt N T a dn gs dt n 0s n gs gs  gsc ngs 是激光腔中光子数密度； n是增益介质反转粒子数密度； n gs是可饱和吸收体基态粒子数密度， nes 是可饱和 吸收体激发态粒 子数密度 ； n 0s是可饱和 I 吸收体总粒子数密度 ; σ和l 分别是增益介质的受激发射截 面和长度；gs 和σes 分别是可饱和吸收体基态和激发态的吸收截面 ；l s是可饱和吸收体沿光腔轴线的长度； R 是输出镜的反射率； 是激光器腔体的耗散性损挺 ； c 是光速； r 为光在腔 中往返一用的时间； 是反转因子 ，对于四能级系统为 1 三能级系统为 2；RP(t ) 是抽运速率；NT 是增益介质的总粒子数密度 ； a是增益介质的上能级寿命 : gs是可饱和吸收体的恢复时间 ； 为了接近实际所添 加的因子（ 1-n/N T） 。 MATLAB中的 ode 求解器 一阶常微分方程的初值，其一般形式为 dy f (t , y) dt  ，有大量的常微分方程 y(t 0 ) y0 虽然从理论上讲，其解是存在的，但我们却无法求出其解析解，此时，我们需要寻求方程的数值解，在求常微分方程数值解方面， MATLAB具有丰富的函数，我 们将其统称为求解器（ solver ），其一般格式为： [T,Y]=solver(odefx,t,y0) 。 该函数表示在区间 t=[t0,tf] 上，用初始条件 y0 求解显式常微分方程。solver 为命令 ode45， ode23，ode113， ode15s， ode23s，ode23t ，ode23tb 之一，这些命令各有特点。 我们这里主要考虑