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