4
3
2
1
0
0.8
0.6
0.4
0.2
0
-0.2
50
45
40
35
30
25
20
15
-0.4
-0.6
10
5
-0.8
0
辅助激光光学分析与应用
作者: 刘良清
二零一三年四月
二零一四年 第三版
Matlab 辅助激光光学分析与应用
Matlab 辅助激光光学分析与应用
作 者:刘良清
单 位:武汉凌云光电科技有限公司
毕业院校:华中科技大学激光技术与工程学院
学 历:硕士研究生
研究方向:自适应光学、非线性光学、激光光学、固体激光
器件、激光工业应用
联系方式:beamcom@gmail.com QQ:19205977
2014 年五月 第三版
II
Matlab 辅助激光光学分析与应用
前 言(序)
每个人都会有自己的理想和梦想,只不过我们常常走着走着就把它们弄丢了。
也许你和我一样,小时候也曾信誓旦旦地说过自己的将来的理想是要当科学家什么的。
那时候不是因为我们理想很远大,其实我们根本就不知道理想是什么。我们不知道真正科学
家到底是什么,也不知道除了科学家之外还有什么可以叫做理想!
当然,我现在没能成为科学家,将来也没有太大的可能会成为科学家。我既不是什么教
授,也不是什么专家,顶多算是一个普通的技术员,普通得不能再普通了。记得研二的时候,
导师还曾担心过我能不能正常毕业。整个大学过得似乎都有些浑浑噩噩,没什么明确的目标,
只是盼望着大学毕业了能否找到一个温饱不愁的工作。现在我知道,这世界上根本就没有哪
一样东西是真正会让你不愁的。
很早,我就有想写一本书的想法。开始想写小说,结果构思来、构思去,发现好的情节
似乎都被作家们写完了。突然有一个新奇的思路的时候,才发现自己文笔真的不怎地,开篇
几段就写不下去了。罢了,作家之梦就这么完蛋了!
研究生的最后一年,为了毕业论文和文章,慢慢自己学会了一点 MATLAB 编程,处理
一些数据和图形。但是,后来发现之前做过的一些研究,突然有一天再去看的时候,竟然怎
么也看不懂。于是,我就想,干脆把自己研究过的或者是纯属娱乐研究过的东西都总结在一
起,把整个思路都写下来,这样可以方便自己以后参考。既然要写,那就得像模像样地写,
有多少就写多少。我不知道什么时候会写完,照这个想法,也许一辈子也写不完。因为,总
会有新的东西值得我们去研究。
那么就这样吧,这就是我的第一本书了。书中肯定会有各种问题和错误,很多结论只代
表我个人观点,如有雷同请找专家。尽信书不是好书!
III
【样品】Matlab辅助激光光学分析与应用(2014年 第三版)淘宝发售(搜索书名可找到)
Matlab 辅助激光光学分析与应用
目 录
MATLAB辅助激光光学分析与应用....................................................................................................I
MATLAB辅助激光光学分析与应用...................................................................................................II
作 者:刘良清 ..............................................................................................................................II
第一章 光的波动性和衍射 ............................................................................................................... 1
1.1 MAXWELL方程组和电磁波........................................................................................................1
1.2 波动方程 .......................................................................................................................................4
1.3 衍射 ...............................................................................................................................................5
1.3.1 小孔衍射................................................................................................................................6
1.3.2 双缝衍射..............................................................................................................................14
1.4 波前畸变 .....................................................................................................................................17
1.4.1 Zernike多项式......................................................................................................................17
1.4.2 畸变光束的衍射..................................................................................................................20
1.5 光束通过光学元件的变换 .........................................................................................................24
1.5.1 平行光束通过透镜的聚焦..................................................................................................24
1.5.2 高斯光束通过透镜的聚焦..................................................................................................30
1.5.3 自聚焦透镜..........................................................................................................................32
1.6 高斯光束 .....................................................................................................................................36
1.6.1 高阶高斯光束......................................................................................................................36
1.6.2 高斯光束的传输变换..........................................................................................................41
1.6.3 高斯光束的干涉..................................................................................................................45
第二章 激光谐振腔 ......................................................................................................................... 60
2.1 激光谐振腔的本征方程 .............................................................................................................60
2.2 无源腔的FOX-LI迭代方法 .......................................................................................................60
2.3 无源腔的矩阵特征向量方法 .....................................................................................................66
2.3.1 平行平面腔..........................................................................................................................67
2.3.2 双凹腔..................................................................................................................................71
2.4 基模谐振腔 .................................................................................................................................78
2.4.1 高斯光束的ABCD定律与谐振腔.......................................................................................80
2.4.2 基模谐振腔的稳定性..........................................................................................................91
2.4.3 基模谐振腔设计实例..........................................................................................................97
2.4.4 基模谐振腔设计的逆向方法............................................................................................107
第三章 激光工程设计 ................................................................................................................... 118
IV
【样品】Matlab辅助激光光学分析与应用(2014年 第三版)淘宝发售(搜索书名可找到)
Matlab 辅助激光光学分析与应用
3.1 速率方程的求解 .......................................................................................................................118
3.1.1 主动调Q速率方程的求解.................................................................................................118
3.1.2 被动调Q速率方程的求解.................................................................................................125
3.1.3 脉冲泵浦速率方程的求解................................................................................................132
3.2 非线性激光光学 .......................................................................................................................136
3.2.1 非线性晶体........................................................................................................................136
3.2.2 非线性晶体相位匹配........................................................................................................161
第四章 MATLAB数据与图像处理............................................................................................... 171
4.1 MATLAB图像处理应用...........................................................................................................171
4.1.1 自制特殊图形....................................................................................................................171
4.1.2 图像数据处理....................................................................................................................179
4.2 MATLAB数据处理应用...........................................................................................................197
4.2.1 批量数据处理....................................................................................................................197
4.3 MATLAB方程求解应用...........................................................................................................211
4.3.1含参数非线性方程组的求解.............................................................................................211
V
【样品】Matlab辅助激光光学分析与应用(2014年 第三版)淘宝发售(搜索书名可找到)
r
ε
=·E
=·B
0
0
(源于毕奥-萨瓦尔定律的高斯定律)
B
¶
¶
t
(法拉第定律)
=
0
´ +
E
(1.2)
(1.3)
Matlab 辅助激光光学分析与应用
第一章 光的波动性和衍射
1.1 Maxwell方程组和电磁波
十八世纪中叶,James Maxwell 将已知的各种电磁作用关系用一组方程组合起来,形
成了一个方程组:
(源于库伦定律的高斯定律)
(1.1)
´ -
B
m
0
e
0
E
¶
¶
t
=
J
(Maxwell 修正的安培定律)
(1.4)
式中, E 和 B 分别代表了电场和磁场分量。电荷密度 r 描述路径空间单位体积内的电荷量
分布;电流 J 描述电荷的移动( 单位电荷乘以速度) 。 0e 表示真空介电常数,其值 为
7
-
e
·T m A
10
0
12
C N m 。 0m 表示真空磁导率常数,其值为
-
8.854 10
m
0
= ´
p
4
=
´
/
/
2
2
·
(或者
k
·g m C )。
/
2
在安培定律中引入了一个关键参数之后,Maxwell 意识到,方程组构成了一个完美的电
磁现象自洽理论。此外,方程组预言了电磁波的存在,并以光速传播。在 Maxwell 时代之
前就已经有人对光速进行了测量,因此一个显而易见的结果(当时还难以令人置信)便是,光
是一种高频振荡表现,类似并超越了支配电流和电荷的影响因素。而在此之前,光学还仍然
作为一种独立于电学和磁学的主体进行讨论的。
这里,我们不再对电磁学的基本知识进行详细的讨论,因为它们在普通物理课程中都有
讲述,并且有大量的文献和书籍对其进行了细致的分析。但我们要简要的从波动方程出发,
求解旁轴近似下的 Maxwell 方程组,得到激光传输与变换的基本方程,以方便我们后续的
讨论和应用。
为了体现 Matlab 在可视化方面的优势,我们先以一个简单的例子作为本书的开篇,以
达到抛砖引玉的效果。在电动力学中,我们会遇到真空电磁场波动方程的旁轴近似解,众所
周知,其解为具有高斯分布的电场复振幅:
Y
( ,
r z
)
= Y
0
2
p
w
0
( )
w z
exp
é
ê
-
ê
ë
j
2
2
kr
2 ( )
( )
R z w z
r
2
-
exp
{
-
[
j kz
-
f
( )
z
}
]
(1.5)
ù
ú
ú
û
作者:大胡子 http://my.oecr.com/35699
1
【样品】Matlab辅助激光光学分析与应用(2014年 第三版)淘宝发售(搜索书名可找到)
Matlab 辅助激光光学分析与应用
式中,
=k
2 /p l
为光波传播常数。 ( )w z 、 ( )R z 、 ( )f z 是与光束有关的传播参数。分别
表示为:
( )
w z
=
w
0
1
æ
ç
+ç
ç
çè
l
z
p
w
2
0
2
ö÷
÷
÷
÷
ø
( )
R z
=
z
æ
ç
+ç
ç
çè
é
ê
1
ê
ê
ë
p
0
w
l
z
ù
22
ö
ú÷
÷
ú÷
÷
ø
ú
û
p
=
f
( )
z
=
tan
1
- æ
ç
ç
ç
çè
z
Z
R
ö÷
÷
÷
÷
ø
,
Z
R
2
0
w
l
光束远场发散角为:
或者
q
0
=
lim
¥
z
( )
w z
z
=
l
p
w
0
q
0
=
æ
ç
ç
ç
ç
è
( )
w z
( )
R z
2
ö
÷
÷ +
÷
÷
ø
æ
ç
ç
ç
ç
è
l
p
( )
w z
2
ö
÷
÷
÷
÷
ø
(1.6)
(0.7)
(0.8)
(0.9)
(0.10)
我们可以用几行简单 matlab 程序就可以画出具有高斯分布的电场强度,如图 1.1 所示,图形
美观,方便对光强的分布有一个感性的视觉认识。程序代码为:
clear;
clc;
w0=0.5;
r=linspace(0,3*w0,200);
Iopt=exp(-2*rho.^2/w0^2);
eta=linspace(0,2*pi,200);
[x,y]=pol2cart(theta,rho);
[rho,theta]=meshgrid(r,eta);
1-
2-
3-
4-
5-
6-
7-
8-
9-
surf(x,y,Iopt);
10- shading interp;
11- xlabel('位置 /mm');
12- ylabel('位置 /mm');
13- zlabel('相对强度 /a.u.');
14- title('高斯强度分布');
15- axis([-3*w0,3*w0,-3*w0,3*w0,0,1]);
16- colorbar;
17- colormap('hot');
18- box on;
19- grid off;
作者:大胡子 http://my.oecr.com/35699
2
【样品】Matlab辅助激光光学分析与应用(2014年 第三版)淘宝发售(搜索书名可找到)
Matlab 辅助激光光学分析与应用
图 1.1 高斯光强分布
另外,我们还可以画出高斯光束在自由传输过程中的强度变化,如图 1.2 所示,程序代码如
下:
clear;
clc;
lambda=1.064e-3;
w0=0.5;
ZR=pi*w0^2/lambda;
[py,pz]=meshgrid(y,z);
y=linspace(-4*w0,4*w0,200);
z=linspace(-2*ZR,2*ZR,200);
wz=w0*sqrt(1+(lambda*pz/pi/w0^2).^2);
1-
2-
3-
4-
5-
6-
7-
8-
9-
10- Iopt=w0^2./wz.^2.*exp(-2*py.^2./wz.^2);
11- surf(pz,py,Iopt);
12- shading interp;
13- xlabel('位置 /mm');
14- ylabel('位置 /mm');
15- zlabel('相对强度 /a.u.');
16- title('高斯强度分布的传输');
17- colorbar;
18- colormap('hot');
19- box on;
20- grid off;
作者:大胡子 http://my.oecr.com/35699
3
【样品】Matlab辅助激光光学分析与应用(2014年 第三版)淘宝发售(搜索书名可找到)