matlab计算单模光纤耦合效率的积分,《高等光学仿真（MATLAB）版》——光波导、激光：前言...

前言

21世纪是信息时代，信息学科和信息产业的迅猛发展促使传统的光学仪器科学向光电信

息学科扩展。现代光电信息学科及其产业的发展要求新一代的科学研究人员与工程技术人

员除了具有扎实的理论基础，还应具有应用所学理论建模并仿真求解光电信息学科及现代

光学中各种问题的能力，成为知识结构新和创新能力强的高层次人才。光波导和激光器是现

代光电信息科学中两个举足轻重研究方向。本书以这两个方向的光学内容为基础，从基本的

物理概念出发，建立相应的理论模型，并将这些光学问题归纳为特征方程求根、积分求解、

常微分方程求解等几类数值求解问题，在对相应的数值分析方法进行简要介绍的基础上结

合MATLAB强大的数值计算和图形显示功能，完成光学问题的仿真计算并给出图形化的显示

结果。也就是通过光学仿真计算，利用MATLAB编程来完成现代光学典型问题的模型求解，

通过数据和图形来展示现代光学问题的本质，力求形成理工结合、经典理论与现代数值方法

紧密结合的新体系。长期以来，由于光学课程中的概念繁多、物理规律较为抽象，理论教学

对实验的依赖性较强，特别是其中的一些光学现象和规律缺乏细致的数学推导，再加上授课

教师一直沿袭传统的口授笔演的教学方式，这些都给学生学习该门课程带来了诸多困难，因

而该门课程的教学效果也总是不尽如人意。在这种情况下，笔者认为，光学课程需要形象生

动的教学，需要现代化的教学手段，千方百计地为学生提供观察光学物理现象的机会，提高

学生学习光学课程的兴趣，培养学生的思维水平和创新能力。

国内光学方面的著作内容多集中在理论方面，鲜有介绍如何对理论模型进行数值求解的

案例。而在科学研究和工程应用中不仅要求能够根据实际情况选择适当的理论建立模型，更

为重要的是要能够结合实际情况仿真求解理论模型，并在此基础上对模型的某些关键参数进

行优化，最终用于指导科学研究和工程应用。本书将MATLAB用于光学仿真中，将具体介绍

如何利用MATLAB来仿真光学中两个举足轻重的研究方向—光波导和激光器中的一系列理论

模型，通过这些仿真过程和结果能够进一步加深对光波导和激光器的理解和应用。本书的大

部分内容已经在北京航空航天大学等高校相关专业的研究生选修课程中讲授过，受到了学生

的普遍欢迎。由于MATLAB语言在很多理工专业的后续课程中有很大的作用，建议有条件的

学校也开设相应的课程。使学生能认识和掌握该语言，提高对光学问题仿真求解的水平，为

更深入的科学研究打下扎实的基本功。

本书作者在光学仿真与计算教学中，通过下列方式将MATLAB与光学课程有机地结合起

来：

一是以MATLAB为平台，开发制作了各种环境的导波光学和激光原理等高等光学现象仿

真程序，并运用于计算机所支持得课堂教学中，以其作为演示实验配合光学理论的讲授，很

好地解决了真实实验因环境限制而不能进入课堂的难题。

二是利用MATLAB的仿真与计算功能，鼓励学生通过自主探索，去研究高等光学课程中

的一些更深入的问题。在掌握理论知识的前提下，让学生建立相应的物理模型和数学模型，

然后利用MATLAB编写程序，去完成对知识的巩固与拓宽。这是一种探索过程，也是为学生

以后的研究工作奠定基础。因此，在本书每一章的最后一节都提供了一定数量的习题，供学

生巩固并加深理解该章的相关光学模型及其仿真计算。

三是利用MATLAB的计算与优化功能，启发学生对数学模型中的参数进行改变，根据

实际物理条件选择符合要求的最优值，并获得最优条件下的参数值，最终通过理论仿真来

指导实践。完成实践(参数获取)—理论(物理模型建立)—仿真(MATLAB计算)—优化

(MATLAB优化)—实践(最优参数选取)的过程，让学生真切感受科学技术是第一生产力。

基于以上指导思想，全书共分为6章。

第1章主要介绍光的电磁理论基础。本章从麦克斯韦方程出发，根据麦克斯韦电磁理论，

利用电矢量和磁矢量来分析光波在两介质表面的反射特性，并结合MATLAB仿真计算光波从

光疏介质进入光密介质，以及光波从光密介质进入光疏介质时的反射率、透射率、相位等随

入射角度的变换关系，得到布鲁斯特角、全反射、倏逝波等基本概念及特性。

第2章采用的是本征模方法，利用有限空间的波动光学理论分析光在介质光波导(理想平

板介质光波导)中传播特性。由于受到介质边界条件的限制，根据不同的边界条件，对麦克

斯韦方程或相应的波动方程求解后，可以得到其特征方程。在推导出理想平板介质光波导最

基本的TE模和TM模的特征方程后，利用MATLAB的方程求根函数即可对特征方程进行数值

求解，从而得到其中TE模(或TM模)的传输特性。

第3章首先介绍了光纤的基本概念；然后从光波在光纤(圆柱光波导)中传输的圆柱坐标

系下亥姆霍兹方程出发，得到阶跃折射率光纤中光波传输的Bessel方程，结合光波在光纤中

传输的边界条件，推导出弱导近似下的特征方程，并利用MATLAB中的Bessel函数以及数值

求根函数，对弱导近似下的特征方程进行数值求解，获得光纤的归一化工作频率、归一化横

向相位参数、归一化横向衰减参数等数据；再通过MATLAB的三维作图功能，将不同参数光

纤中的电场分布特性展示出来。

第4章着重介绍了高斯光束和光纤的耦合。由于高斯光束是光场亥姆霍兹方程在缓变振幅

近似下的一个特解，该特解不仅可以用来描述基模激光在空间中的传输，而且可以用来近似

描述从光纤端面出来的LP01模的辐射场，具有重要的物理意义。通过介绍高斯光束的基本性

质、复参数表示、ABCD定律，得到高斯光束通过复杂光学系统以及薄透镜的变换，并得到

其聚焦特性；在此基础上给出了光纤端面辐射场，并结合MATLAB仿真辐射场的函数特性；

最后给出了光纤的光功率发射和耦合，并用MATLAB的积分求解函数计算了LED与单模光纤

的耦合效率以及光纤与光纤连接的耦合效率。

第5章从激光的基本原理出发，介绍了辐射与物质的相互作用，关于自发辐射、受激辐射

和受激吸收的爱因斯坦关系式，吸收与光学增益，激光器的基本构成，激光速率方程和激光

调Q技术等内容。以一种典型的被动调Q的微晶片激光器为例，给出了其被动调Q的速率方程

组，该速率方程组可以简化为一个具有3个自变量的常微分方程组。利用MATLAB的常微分

方程初值问题求解函数即可对该被动调Q速率方程组进行仿真求解，得到被动调Q的微晶片激

光器的脉冲时域特性以及被动调Q过程中光子数密度和反转粒子数密度随时间的变化关系。

第6章重点介绍了一种新型的激光器—高功率双包层光纤激光器。光纤激光器是以掺杂

光纤作为增益介质的一类激光器。光纤激光器和其他类型的激光器一样，由能产生光子增益

的工作介质、使光子得到反馈并在工作介质中进行谐振放大的光学谐振腔和激励光子跃迁的

抽运源3部分组成。只不过光纤激光器的工作介质是同时起着波导作用的掺杂光纤。因此，光

纤激光器是一种波导型的谐振装置。在本章中分别给出了端面抽运和侧面抽运两种情况下的

双包层光纤激光器的理论模型，利用MATLAB的常微分方程边值问题求解函数，并结合双包

层光纤激光器的边值条件对其速率方程组进行仿真求解，得到抽运光、激光以及反转粒子数

密度沿光纤长度的变化。

书中大量运用到求解各类模型的数值计算方法，其中主要有方程求根的数值解法、数值

积分方法、常微分方程的初值问题数值求解、常微分方程的边值问题数值求解。考虑到学习

本书的读者数学方面的知识背景不尽相同，为了能够让读者在尽量少的时间内能够理解数值

求解的基本思路和方法，并运用MATLAB相应的模型进行数值求解，在每章的习题前都有一

节用于介绍数学和MATLAB相关方面的补充知识，从而减轻读者在学习过程中的负担。

本书全部源程序均在MATLAB R2008a中测试运行通过。程序源代码请到北京航空航天大

学出版社(http://www.buaapress.com.cn)下载中心下载。

同时，北京航空航天大学出版社联合MATLAB中文论坛(http://www.iLoveMatlab.cn)为本书

设立了在线交流版块，网址为https://www.ilovematlab.cn/forum-203-1.html，有问必答！作者会第

一时间在MATLAB中文论坛勘误，也会尽量根据读者要求陆续上传更多案例和相关知识。希

望这本不断“成长”的书能最大限度地解决您在学习、研究、工作中遇到的光学方面相关问

题。

由于作者水平有限，加之时间仓促，书中不免存在不足和疏漏，敬请读者批评指正。

编者

2011年1月