研究生数学建模竞赛优秀论文选《面向下一代通信的VCSEL激光器仿真模型》457Word格式.docx

1、𝑃, 𝜀（）见表 7-2），绘制 VCSEL 带宽模型的带宽响应曲线并分析，得出温度和偏置电流对带宽曲线的影响规律。对拟合参数进行敏感性分析，提出改变参数以得到更宽带宽的方案。针对问题 4：通过文献调研的方式，建立了基于激光器速率方程的小信号等效电路模型，提出了该模型的参数拟合思路；结合光功率转换效率𝑘的物理意义， 计算出了一套增大带宽的参数设计方案。关键词：VCSEL L-I 模型 带宽模型 遗传算法 误差分析1.1问题背景1 问题重述随着经济技术的发展，对互联网技术的要求越来越高，我国家庭固定网络速度已经快速发展到了今天的百兆，甚至千兆的光纤

2、宽带入户阶段。“光纤宽带入户”就是采用光纤来传输信号。光纤中传输的信号有激光信号和电信号，其中激光信号具有远高于电信号传输速率的特点，更适合于未来高速率的传输网络。在光纤通信传输系统设计前，往往会通过计算机仿真的方式研究系统设计的指标， 以便快速找到最适合的解决方案。因此在进行系统仿真时，需要准确掌握系统中各个器件的特性以保证仿真模型的精度，其中激光器作为光纤通信系统的核心器件是系统仿真中需要考虑的一个重要因素。1.2问题的提出尽管与我们生活息息相关的激光器种类繁多，但由于垂直腔面发射激光器（VCSEL: Vertical Cavity Surface Emitting Laser）具有使用简

3、单，功耗较低等特点，一般 VCSEL 的工作电流在 6mA8mA，所以本题目需要在满足如下要求时，完成本文的主要任务得到能准确反映 VCSEL 激光器特性的数学模型。（1）问题的要求针对垂直腔面发射激光器（VCSEL）的特点它的输出的光功率强度与器件的温度相关，当器件温度（受激光器自身发热和环境温度的共同影响）改变后， 激光器输出的光功率强度也会相应发生变化，在进行建模时，需要满足如下要求：1）准确反映 VCSEL 激光器特性；2）考虑激光器输出的功率强度与温度的关系，即该激光器可以在多大的外界环境温度范围内使用；3）如何设计激光器的参数使激光器具有更大的传输带宽，即 S21 曲线上纵坐标-1

4、0dB 位置对应的横坐标频率值更大。（2）待研究问题问题 1：基于 VCSEL 的 L-I 模型的经验公式，根据实测数据，确定模型参数（𝜂4），画出 10，20，30，90等温度下的L-I 曲线，并根据建立的 L-I 模型推测出在保证用户可以正常使用网络的条件下， 电信机房里 VCSEL 激光器的最高的工作环境温度。问题 2：根据问题 1 中模型精度和误差产生的原因对模型提出改进，并根据改进的模型画出 10，20，30，90等温度下的 L-I 曲线与问题 1 中的 L-I 曲线进行比较，或者通过完全不同方法得到更好的 L-I 模型。问题 3：根据 VCSEL 激光器小信号幅频响

5、应曲线数据和相应的驱动电流、输出光功率数据，建立恰当的激光器小信号幅频响应参数模型，给出参数构成及其确定方法，画出不同环境温度和不同偏置电流下的带宽响应曲线，其中一条必须是 20下 7.5mA 偏置电流的仿真输出曲线；利用该模型分析激光器的温度和激光器的偏置电流对器件带宽曲线的影响；提出除了通过改变环境温度和偏置电流来获得更宽带宽激光器的方法；确定通过改变哪些参数会有激光器在 3dB 范围内的部分频率处幅度高于 0 频位置的现象。问题 4：探索更好的带宽模型建模方式，使得模型运算速度更快；或者在相同的温度和偏置电流下，可以获得更宽的 3dB（或 10dB）带宽。若有，给出建模方案，包括可能的数

6、学公式，不同温度和偏置电流下的带宽响应曲线，并与问题3 的模型进行比较。2 模型假设（1）转换效率𝜂（𝑇）受温度影响较小，即𝜂）近似于常数𝜂；（2）加载于激光器的电流为直流电，即𝐼 = 0；（3）二极管热电压（diodes thermal voltage）不随温度的变化而变化，即𝑉𝑇为一常数；（4）阀值电流等于热偏置电流与一个常数的和，即𝐼（𝑁, 𝑇） = 𝐼𝑜𝑓） + 𝐼

7、0。3 符号说明符号意义单位P0激光器输出光功率的测定值mW，𝑃𝑆L-I 模型光功率的计算值I注入到激光器的外部驱动电流mA ， （T ） ，i转换效率无因次Ith （ N,T ）阈值电流N载流子数Ioff （T ）与温度相关的经验热偏置电流Ith0常数RthVCSEL 热阻抗/mWth|热时间常数T0环境温度KV ，V0输入电压的计算值和测定值V𝑉二极管热电压𝑅串联电阻𝐼二极管饱和电流时间sN0透明载流子数Ndb振幅dB4.1基本名词含义4 模型准备（1）垂直腔面发射激光器（VCSEL）垂直腔面发射激光器（Ver

8、tical-Cavity Surface-Emitting Laser，简称 VCSEL， 又译垂直共振腔面射型雷射）是一种半导体激光器，能将电信号转化为激光信号。与传统的激光器相比，VCSEL 具有使用简单，功耗较低等特点。根据题目要求， 本文将以这种激光器作为研究对象。（2）信号在物理或工程上可以将信号理解成变化着的物理量，最常用的是电压、电流、磁通量、光通量等。例如电信号可以通过电流的幅度、频率、相位变化来表示不同的消息。本文中研究的是光信号，而研究的重点是光信号的幅度和频率。（3）系统一般输入的信号和最终得到的信号不完全相同，而造成这种变化的中间环节的统称为系统。若系统的功能不随时间变

9、化，那么可以将该系统建模为线性常系数微分方程。通常该模型是通过机理分析法建立的，然后用输入和其对应的输出数据，通过数学方法求解参数，得出模型。本文模型的确定主要就是使用这种方法。（4）带宽在模拟信号系统中，带宽其实就是信号占用的频谱的度量，能标识传输信号所占有的频率宽度。其中频率宽度是指传输信号的最高频率和最低频率的差值。4.2数据预处理题中给出的数据可能与模型中的参数的单位不同，或者为了方便编程，需要对测试数据进行参数的转化。转化后的各个参数的单位见符号说明。需要说明的是振幅的单位 dB 与功率的单位 mW 之间的关系。根据调研文献1，它们之间的有如下换算关系：0dBm = 10 log （

10、1）mW=1mW x xxdBm = 10 log1010 mW=1010 mW ydBm-zdBm= （ y - z ）dB（4-1）5 问题 1：VCSEL 的 L-I 模型的参数确定激光器是将电能转换成光能的半导体器件，能量转换的过程，也是电子的电能转换为光子的光能的过程，在转换过程中，伴随着电子的运动，半导体器件会产生一定的热量。从能量守恒的角度看，转化为热能的能量越多（发热导致能量浪费了），器件温度越高，那么转化为光能的能量越少（输出光功率越低），可以利用的能量就越少。因此需要对激光器输出功率与工作电流的关系进行研究。5.1问题描述及分析本题的目的是找出一种方法，确定 L-I 模型的

11、参数𝜂、𝐼0、𝑅、𝑎0、𝑎1、𝑎2、𝑎3和𝑎4，再根据确定了参数的 L-I 模型绘制 L-I 曲线，并且计算在保证用户正常使用的情况下，VCSEL 激光器的最高工作温度。此题实际上为非线性多参数函数的拟合问题。题中已经给出了 1401 组数据，每组数据包括测得的驱动电流、电压和对应的光功率，其中温度为固定值 20。本题的关键是如何拟合 L-I 模型的 8 个参数，所以考虑应用拟合参数的一般方法最小二乘法进行拟合。由于需要拟合的函数是非线性的，而且需要拟合的参数一共有

12、8 个，所以找到合适的算法将是参数能否成功拟合的关键。题中要求绘制在不同环境温度下，光功率和电流强度的关系，而根据题中附录一，L-I 模型最后是光功率𝑃0和驱动电流𝐼、电压𝑉和环境温度𝑇0的函数关系，关系式中多了一个变量𝑉，所以找到一种方法消除这个变量是关键。最后一个要求是求取在保证用户正常使用的情况下，VCSEL 激光器的最高工作温度。可以在 L-I 模型确定后，通过改变环境温度，看不同电流强度下的最大光功率是否大于 2mW 即可找出最高温度。综上所述本题的关键是：（1）找到合适的算法拟合参数；（2）如何消除电压

13、这一变量。5.2L-I 模型的建立L-I 模型，即激光器的工作电流与输出光功率强度关系模型，国际上很多研究机构都对此模型做了大量研究，目前已经有一个经验公式获得了大多数人的认可P = （T ）（I - Ith （N,T ）（5-1）假设转换效率𝜂同时令𝐼） =0 + 𝐼），那么式（5-1）可以简化为将𝐼）表示为P = （I - Ith0 - Ioff （T ）（5-2）（T ） = a Tnoff nn=0（5-3）式（5-3）中的温度𝑇受外界环境温度𝑇0和自身的温度影响，自身的温度与器件产生的瞬时

14、功率𝐼相关T = T + （IV - P） R - dT（5-4）0 th th dt所以 VCSEL 的经验模型为P = （I - Ith0 - Ioff （T ）（T ） = a T n off n n=0 dT（5-5）T = T0 + （ IV - P） Rth - thdt5.3基于 L-I 模型的参数拟合模型的建立根据已经建立的 L-I 模型，将式（5-3）和式（5-4）代入式（5-2）可得P = I - I-aT +IV - PR n th00 （ ） th （5-6）因为在直流电中，式（5-4）中的𝑑𝑑 = 0，所以省略。一般情况

15、下，级数取前四项即可满足精度要求，所以式（5-6）可写为 4 n P = I - Ith0 - an T0 +（IV - P）Rth （5-7） n=0 分析式（5-7）可知，上述经验公式是一个𝑃关于𝐼，𝑇0 ，𝑉的隐函数𝑃 =（𝐼, 𝑉0, 𝑃）。在实际中我们知道，电流和功率不可能小于 0，所以有非负约束P 0 0 th0（5-8）以光功率的拟合值与实测值的离差平方和最小为目标函数，参数𝜂、1、𝑎4为决策变量，式（5-5）和式（5-8

16、）为约束条件建立光功率0和电流𝐼、环境温度𝑇0的函数关系参数的拟合模型Min（ P - P ）2I （T ） = a T n off n（5-9）s.t. T = T0 + （ IV - P） Rth - th dtP 0 Ith0 05.4L-I 模型的参数拟合模型的求解（1）求解思路在建立的光功率𝑃0的函数关系参数的拟合模型中，𝜂、4均为待定常系数。根据给定的测试数据，通过经验模型参数的一般确定方法，即参数拟合，运用遗传算法确定各个参数的值。（2）算法设计与模型参数确定对光功率和电流、环境温度的函数关系参数的拟合模型，以离差平

17、方和最小为优化目标，获得拟合参数。具体计算步骤如下：Step1：读取𝐼、𝑈、𝑃0、𝑇数据文件； Step2：设置遗传算法参数集 gaoptionset；Step3：调用 MATLAB 软件的 ga 函数和 Eval 评价函数；Step4：令𝑉𝑠𝑢𝑚 = 0，𝑖 = 0；Step5： = 𝑖 + 1，𝑇 = 𝑇0 + （𝐼 𝑃0）𝑅；𝑗=1St

18、ep6：计算𝐼 = 4 𝑎j 𝑇Step7：计算𝑃 = （𝐼 𝐼 𝐼），如果 P4?YP=0?i=1401Vsum=Vsum+（P-P0）2，i=i+1P=0P=*（I-Ith0-Ioff）Ioff=Ioff+aj*Tj，j=j+1i=i+1，T=T0+（IV-P0）*RthVsum=0，i=0调用Eval评价函数调用Matlab的ga函数设置遗传算法参数集gaoptionset读取I、U、P0、Ta数据文件输出结果图 5-1 程序计算流程图图 5-2 遗传算法收敛图从实测光功率与拟

19、合光功率的对比图（图 5-3）可以看出，实测曲线和拟合曲线基本重合，拟合效果较好。图 5-3 实测光功率与拟合光功率对比图5.5电压和电流函数关系的确定从 5.3 节中可以得出 VCSEL 的经验模型是一个𝑃，𝑉0的隐函数𝑃），该等式中的变量有四个，分别是光功率𝑃、电流𝐼0。为了求得在不同环境温度𝑇0下，光功率𝑃和电流𝐼的关系，必须消除其中的一个变量。根据调研，电压𝑉存在函数关系2。（1）电压和电流函数关系的参数拟合模型建立据研究，𝑉为

20、（𝐼）的函数2，即𝑉 = 𝑓）V = IR I + V ln 1 +（5-10）S T IS 根据假设（3）𝑉为常数和实际情况，电阻、电流和电压不能小于 0，所以有非负约束RS 0VT 0（5-11） S以电压的拟合值与实测值的离差平方和最小为目标函数，𝑅、𝑉和𝐼为决策变量，式（5-10）和式（5-11）为约束条件建立电压𝑉的函数关系参数的拟合模型（V - V0 ） I V = IRS + VT ln 1 + I R SV 0 S （5-12） T（2）求解思路IS

21、0在建立电压𝑉的函数关系参数的拟合模型中，𝑅均为待定常系数。根据给定的测试数据，通过经验模型参数的一般确定方法，即参数拟合， 运用 MATLAB 软件工具箱 curve fitting 拟合确定各个参数的值。（3）模型参数确定运用 MATLAB 软件工具箱 curve fitting 进行面板操作即可拟合得到参数𝑅的值分别为 0.05788、0.1501 和 9.710-6，其中判定系数R2 = 0.9946。参数拟合对比如图 5-4 所示，拟合参数如表 5-2 所示，从对比图和拟合优度的值可以看出拟合值与实测值误差较小，拟合效果较好。表 5-2

22、 拟合值及拟合优度值0.057880.15019.710-6SSE1.044R-square0.9946Adjusted R-squareRMSE0.02732图 5-4 实测电压与拟合电压对比图5.6基于电压和电流函数关系的 L-I 曲线绘制（1）改进的 L-I 模型及求解思路从 5.3 节可得到𝑃0与𝐼0的函数关系，从 5.5 节中可以得到电压𝑉和电流的函数关系，即电压可以用电流来表示。通过化简可以得到𝑃）的隐函数关系；然后给定不同的环境温度的值，即可得到仅仅是光功率和电流的函数关系；最后通过 MATLAB 软件 fsolve

23、 函数进行数值求解，得到不同电流下光功率的值。（2）算法设计与模型参数确定 L-I 曲线绘制已经可以得出函数关系𝑃），其中变量有𝑃，𝐼给定一个环境温度，再给出一组电流的值，即可根据 fsolve 函数计算出一组光功率的值；再给定另一个环境温度的值，即可算出这个环境温度下，一组电流与光功率的值，如此循环，最后得到不同环境温度下，不同电流下的光功率，根据这些数据绘制不同温度下的 L-I 曲线。具体的计算步骤如下：给定一组环境温度𝑇，和电流初值𝐼Step2：令𝐼 = 𝐼 + 0.01；给定

24、光功率初值𝑃 = 1，调用 fsolve 函数求𝑃, 𝑃）的值，如果𝐼 28，则进入 Step2；否则进入 Step4；程序计算流程如图 5-5 所示。输入环境温度T和电流初值I=0令I=I+0.01给定光功率初值P=1，调用fsolve函数求PI 28图 5-5 程序计算流程图通过 Matlab 编程（代码见附录 B）得出不同温度下的 L-I 曲线。从不同温度下的 L-I 曲线图（图 5-6）中可以看出：在同一温度下，随着电流的不断增大， 光功率的数值是先变大，后变小的；当电流增大到 22mA 后，所有温度下的光功率都等于 0；

25、在同一电流强度下，随着温度的升高，光功率不断减小；随着温度的降低，每个温度下最大光功率对应的电流逐渐增大。图 5-6 不同温度下和实测数据 L-I 曲线图5.7保证用户正常使用最高环境温度的确定（1）模型建立根据得到的 L-I 模型和电压与电流的函数关系，可以得到𝑃0的函数关系为𝑃），换一种表达方式即为𝑇0 = 𝑓1（𝐼结合题中所给条件， 当电信机房里 VCSEL 激光器在直流输入时输出的平均光功率低于 2mW 时，用户的光猫无法检测到信号。所以激光器在所求的最高工作温度下，其光功率不能小于 2mW，所以有约束P0 2在实际中电流大小不可能小于 0，