没有合适的资源?快使用搜索试试~ 我知道了~
Matlab辅助激光光学分析与应用(最终完整版20140512).pdf
需积分: 5 1 下载量 55 浏览量
2024-04-29
20:12:10
上传
评论
收藏 10.97MB PDF 举报
温馨提示
试读
213页
Matlab辅助激光光学分析与应用(最终完整版20140512).pdf
资源推荐
资源详情
资源评论
Matlab 辅助激光光学分析与应用
第一章 光的波动性和衍射
1.1 Maxwell方程组和电磁波
十八世纪中叶,James Maxwell 将已知的各种电磁作用关系用一组方程组合起来,形
成了一个方程组:
0
ε
r
=·E
(源于库伦定律的高斯定律) (1.1)
(源于毕奥-萨瓦尔定律的高斯定律) (1.2)
=·B 0
0
¶
´ + =
¶
B
E
t
(法拉第定律) (1.3)
0
0
e
m
¶
´ - =
¶
BE
J
t
(Maxwell 修正的安培定律) (1.4)
式中, 和 分别代表了电场和磁场分量。电荷密度
r
描述路径空间单位体积内的电荷量
分布;电流
J
描述电荷的移动(单位电荷乘以速度)。 表示真空介电常数,其值
为
。
表示真空磁导率常数,
其值为
(或者 )
。
E
8.854
B
10=´
0
e
12 2 2
0
/e
-
·CNm
2
k/·gm C
0
m
7
0
410 /mp
-
=´ ·Tm A
在安培定律中引入了一个关键参数之后,Maxwell 意识到,方程组构成了一个完美的电
磁现象自洽理论。此外,方程组预言了电磁波的存在,并以光速传播。在 Maxwell 时代之
前就已经有人对光速进行了测量,因此一个显而易见的结果(当时还难以令人置信)便是,光
是一种高频振荡表现,类似并超越了支配电流和电荷的影响因素。而在此之前,光学还仍然
作为一种独立于电学和磁学的主体进行讨论的。
这里,我们不再对电磁学的基本知识进行详细的讨论,因为它们在普通物理课程中都有
讲述,并且有大量的文献和书籍对其进行了细致的分析。但我们要简要的从波动方程出发,
求解旁轴近似下的 Maxwell 方程组
,得到激光传输与变换的基本方程,以方便我们后续的
讨论和应用。
为了体现 Matlab 在可视化方面的优势,我们先以一个简单的例子作为本书的开篇,以
达到抛砖引玉的效果。在电动力学中,我们会遇到真空电磁场波动方程的旁轴近似解,众所
周知,其解为具有高斯分布的电场复振幅:
[
{}
22
0
0
2
2
(,) exp exp ()
() 2 () ()
f
p
作者:刘良清
1
]
é ù
êú
Y=Y - - --
êú
ë û
w
kr r
rz j jkz z
wz Rz w z
(1.5)
Matlab 辅助激光光学分析与应用
式中, 为光波传播常数。 、 、 是与光束有关的传播参数。分别
表示为:
2/
pl
=k
()wz ()Rz ()
f
z
2
0
2
0
() 1
l
p
æö
÷
ç
÷
=+
ç
÷
ç
÷
ç
èø
z
wz w
w
(1.6)
2
2
0
() 1
p
l
é ù
æö
ê ú
÷
ç
÷
=+
ç
ê ú
÷
ç
÷
ç
èø
ê ú
ë û
w
Rz z
z
(0.7)
2
1
0
() tan ,
p
f
l
-
æö
÷
ç
÷
=
ç
÷
ç
÷
ç
èø
R
R
w
z
zZ
Z
=
(0.8)
光束远场发散角为:
0
0
()
lim
l
q
p
¥
==
z
wz
zw
(0.9)
或者
22
0
()
() ()
l
q
p
æöæ
÷÷
çç
=÷+
çç
÷
çç
÷
çç
èøè
wz
Rz wz
ö
÷
÷
÷
ø
(0.10)
我们可以用几行简单 matlab 程序就可以画出具有高斯分布的电场强度,如图 1.1 所示,图形
美观,方便对光强的分布有一个感性的视觉认识。程序代码为:
1- clear;
2- clc;
3- w0=0.5;
4- r=linspace(0,3*w0,200);
5- eta=linspace(0,2*pi,200);
6- [rho,theta]=meshgrid(r,eta);
7- [x,y]=pol2cart(theta,rho);
8- Iopt=exp(-2*rho.^2/w0^2);
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
-
colo
rm
ap
(
'ho
t');
18- box on;
19- grid off;
作者:刘良清
2
Matlab 辅助激光光学分析与应用
图 1.1 高斯光强分布
另外,我们还可以画出高斯光束在自由传输过程中的强度变化,如图 1.2 所示,程序代码如
下:
1- clear;
2- clc;
3- lambda=1.064e-3;
4- w0=0.5;
5- ZR=pi*w0^2/lambda;
6- z=linspace(-2*ZR,2*ZR,200);
7- y=linspace(-4*w0,4*w0,200);
8- [py,pz]=meshgrid(y,z);
9- wz=w0*sqrt(1+(lambda*pz/pi/w0^2).^
2);
10
-
Iop
t=
w0
^2
./w
z.
^2
.*
ex
p(-
2*
py
.^
2./wz.^2);
11- surf(pz,py,Iopt);
12
-
shad
in
g i
nte
rp;
13- xlabel('位置 /mm');
14- ylabel('位置 /mm');
15- zlabel('相对强度 /a.u.');
16- title('高斯强度分布的传输');
17- colorbar;
18- colormap('hot');
19- box on;
20- grid off;
作者:刘良清
3
Matlab 辅助激光光学分析与应用
图 1.2 高斯光束自由传输强度变化
以上我们以简单的例子展示了 Matlab 在可视化方面的强大功能,但本文不再对 Matlab
的基本功能和语法常识进行介绍,我们认为本书的读者已经具备了基本的 Matlab 编程技巧。
或者说,我们所做的只是将我们的实际运用跟读者进行交流讨论,促进大家共同进步。当然,
我们会在一些比较关键的地方指出编程过程中需要注意的问题。
1.2 波动方程
当 Maxwell 统一了电磁理论以后,他马上意识到,波动可能是该方程组的解的形式。
事实上,他希望找到一组满足波动形式的方程组,以辅助他完成找到真正的波动方程。既然
已经知道了光是以波动方式传播的,基尔霍夫首先注意到了
00
1/
em
正好给出了精确的光
速 (之前就已经被测量过),并且法拉第和克尔已经观测到强磁场和强
电场会影响光在晶体中的传播。对初始接触 Maxwell 方程组的人来说,并不能一眼就看出
它的解具有波动形式。但是经过适当的数学操作,我们就可以将它变为波动方程的形式。
8
3.00 10 /
=´
cms
我们来推到电场E的波动方程,磁场B的波动方程的推到过程是类似的。我们将方程
(1.3)进行卷积,可得:
()()
0
¶
´ ´ + ´ =
¶
EB
t
(1.11)
该方程可以由矢量微分恒等式简化:
(1.12)
()()
2
•´ ´ = -EEE
卷积 可由(1.4)
´B
式代换,由此得到:
()
2
00 0
•
em m
æö
¶¶
÷
ç
- + + =
÷
ç
÷
ç
èø
¶¶
E
EE
J
tt
0
(1.13)
作者:刘良清
4
再由(1.1)式代入上式,经过整理就可得到:
Matlab 辅助激光光学分析与应用
2
2
00 0
2
0
r
em m
e
æö
¶¶
÷
ç
÷
- = +
ç
÷
ç
÷
ç
¶¶
èø
E
E
J
tt
(1.14)
需要指出的是,上式中没有考虑到介质的极化。若考虑到介质的极化和实际一般光学问题中
自由电荷为零的条件,上式修正为:
(
22
2
00 0 0
22
0
1
•
em m m
e
¶
¶¶
- = + -
¶¶¶
E
E
free
J
P
P
ttt
)
(1.15)
式中, 为极化强度矢量。这样我们得到了一般的电场传播方程,该方程在非线性光
学中有很重要的地位。当光在真空中传播时,式(1.15)
P
中的右边所有项均为零,方程简化
为:
2
2
00
2
0
em
¶
- =
¶
E
E
t
(1.16)
这样我们就得到了电场传播的波动方程形式。当然在有些实际问题中,式(1.15)中右
边的项并不是都为零,至少会有一项不为零,这与介质的性质有关。
1.3 衍射
考虑一个振动频率为 的光场,其复振幅可以表述为 ,则它也必须满足波
动方程:
w
()
w-
E
jt
re
() ()
22
2
22
0
w
w
-
-
¶
-
¶
EE
jt
jt
ne
re r
ct
=
(1.17)
由于(是假设)电场振幅的含时部分是显式给出的,则方程(1.17)可以简化为:
(1.18)
() ()
22
0- =EErk r
式中 是波矢量的大小。(1.18)
/wºkn c
式就是所谓的赫姆霍兹方程。如果我们忽略
波动的矢量特性,而只考虑它的振幅(这里不再详细讨论其过程),那么在标量近似下,就得
到了标量赫姆霍兹方程:
(1.19)
() ()
22
0- =Er kEr
然后,我们考虑一束沿z轴传播的光束,它的电场复振幅写成 的形式。
我们将它代入标量赫姆霍兹方程(1.19)
()
,,
jk
z
Exyze
式,得到:
22 2
22 2
2
æö
¶¶ ¶¶
÷
ç
÷
++ + =
ç
÷
ç
÷
ç
¶¶ ¶¶
èø
jkz
EE EE
jk e
xy zz
0
(1.20)
作者:刘良清
5
在旁轴近似下,有
2
2
2
¶¶
¶¶
EE
k
zz
。即是说,我们假设了电场的复振幅沿 z 轴传播方
向是缓慢变化的,与平面波类似。但是我们允许振幅沿 z 轴在远大于波长量级的范围上有明
剩余212页未读,继续阅读
资源评论
天天Matlab科研工作室
- 粉丝: 2w+
- 资源: 7257
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功