计算物理实习-Read.DOCVIP

计算物理实习 （彭加福，0640502112，江苏科技大学，数理学院，应用物理） （李文凯，0640502109，江苏科技大学，数理学院，应用物理） （徐大程，0640502115，江苏科技大学，数理学院，应用物理） 导师：周青春 题目：光纤耦合器中光孤子传输的仿真研究 2009年4月2日 一、Matlab仿真理论来源 文献：施娟,侯韶华. 光纤耦合器中光孤子传输的仿真研究[J]. 电子元器件应用. 2008,Vol.10,No.4:65-67. 光纤耦合器中光孤子传输的仿真研究 施娟，侯韶华 (南京邮电大学光电工程学院，江苏 南京210003) 摘要：光纤耦合器因其在光纤通信中的广泛应用而得到深入研究。文中在分析了光脉冲耦合器中光孤子传榆特性的基础上，给出了求解光脉冲在N芯光纤耦合器中传榆信号的耦合模方程组的对称分步傅里叶解法， 同时给出了采用此方法将光脉冲在双芯和三芯耦合器中进行传输演化的仿真结果。 关键词：光纤耦合器；分步傅里叶法；耦合模仿真；光孤子 ——详细内容请查看该文章 二、公式及Matlab仿真 1、描述光脉冲在N芯耦合器中传输的耦合方程组如下： 上式中，是第n个纤芯中模场的慢变振幅，表示线性失配，、和分别为二、三阶色散和损耗，为非线性系数，是纤芯n和j之间的线性耦合系数，为模间色散。 2、计算脉冲在传输了步长h后的表达式为： 其中，一式表示光脉冲传输前h/2时只受色散影响，二式表示非线性对光脉冲在步长h内的影响，三式表示光脉冲在传输后h/2只受色散影响。经过这三步运算，就可最终得到光脉冲在传输步长h后的数值解。二式中用到未知数，它在开始计算时可用代替，然后采用迭代算法。 3、上三式中均含，它是关于频率的函数，正确计算圆频率很关键，因此要深入了解Matlab中的快速傅立叶变换思想。Matlab仿真时采用快速傅立叶变换fft()和ifft()两个函数，变换时根据正负频选用正逆傅立叶变换函数，防止得到错误的对称解。 三、Matlab仿真结果 1、当光脉冲在单根光纤中传输时，取： 初始脉冲A=sech (t) 光脉冲在单根光纤中传输 由图可见，在不考虑损耗的情况下，光孤子在单根光纤中传输时，保持原始形状不变，向前持续传输，这是孤子的最主要特征和研究价值所在。 2、当光脉冲在双芯耦合器中传输时，取： 也就是说，在光纤反常色散区只考虑二阶色散，而忽略损耗、高阶色散和模间色散。 光脉冲在双芯耦合器的第一根光纤中传输 光脉冲在双芯耦合器的第二根光纤中传输 3、当脉冲在三芯耦合器中传输时，取： 也就是说，色散中只考虑二阶色散，耦合只存在相邻纤芯间。 光脉冲在三芯耦合器的第一根光纤中传输 光脉冲在三芯耦合器的第二根光纤中传输 光脉冲在三芯耦合器的第三根光纤中传输 4、结果分析： 光脉冲在光纤耦合器中传输时，随着传输距离的增加，脉冲宽度也不断地展宽，这是由于光纤中色散和非线性效应共同作用的结果。需要注意的是，计算精度不够时会将这一结果放大，甚至会很快扩散开，并因Matlab中快速傅立叶变换的周期性边界条件,致使能量从时间窗口的一边跑到另一边，而导致错误结果，因此选择适当的精度也是很关键的。 从图中可以看出，输入的能量在耦合器的光纤中传递，保持总能量不变，光脉冲形状由于色散和非线性效应而展宽，但仍可持续传递。在三芯耦合器中，中间纤芯的能量变化周期约是边上两根纤芯能量变化周期的2倍。 四、Matlab仿真程序 *PART ONE : ****** Wspace.m by Thomas E. Murphy ****** function w = wspace(t,nt); % This