光纤耦合器中光孤子传输的matlab仿真研究.doc
1计算物理实习(彭加福,0640502112,江苏科技大学,数理学院,应用物理)(李文凯,0640502109,江苏科技大学,数理学院,应用物理)(徐大程,0640502115,江苏科技大学,数理学院,应用物理)导师周青春题目光纤耦合器中光孤子传输的仿真研究2009年4月2日一、MATLAB仿真理论来源文献施娟,侯韶华光纤耦合器中光孤子传输的仿真研究J电子元器件应用2008,VOL10,NO46567光纤耦合器中光孤子传输的仿真研究施娟,侯韶华南京邮电大学光电工程学院,江苏南京210003摘要光纤耦合器因其在光纤通信中的广泛应用而得到深入研究。文中在分析了光脉冲耦合器中光孤子传榆特性的基础上,给出了求解光脉冲在N芯光纤耦合器中传榆信号的耦合模方程组的对称分步傅里叶解法,同时给出了采用此方法将光脉冲在双芯和三芯耦合器中进行传输演化的仿真结果。关键词光纤耦合器;分步傅里叶法;耦合模仿真;光孤子详细内容请查看该文章二、公式及MATLAB仿真1、描述光脉冲在N芯耦合器中传输的耦合方程组如下2,162,1231TAIKAITAITAZJNJNNJNNNNNN上式中,是第N个纤芯中模场的慢变振幅,表示线性失配,、和NN1 23分别为二、三阶色散和损耗,为非线性系数,是纤芯N和J之间的线性耦合系数,为模间NNJKNJ色散。2、计算脉冲在传输了步长H后的表达式为},2,2{EXP,,11NNJJJNNNTZAFCHITZAFGTZA,2HIZN},,{E,,11NNJJJNNNHZIZHZ其中,一式表示光脉冲传输前H/2时只受色散影响,二式表示非线性对光脉冲在步长H内的影响,三式表示光脉冲在传输后H/2只受色散影响。经过这三步运算,就可最终得到光脉冲在传输步长H后的数值解。二式中用到未知数,它在开始计算时可用代替,然后采用迭代算法。2,THZAN2,TZAN23、上三式中均含,它是关于频率的函数,正确计算圆频率很关键,因此要深入了解MATLAB中的NG快速傅立叶变换思想。MATLAB仿真时采用快速傅立叶变换FFT和IFFT两个函数,变换时根据正负3频选用正逆傅立叶变换函数,防止得到错误的对称解。5,4三、MATLAB仿真结果1、当光脉冲在单根光纤中传输时,取SEC,0THTA2IGN1初始脉冲ASECHT光脉冲在单根光纤中传输由图可见,在不考虑损耗的情况下,光孤子在单根光纤中传输时,保持原始形状不变,向前持续6传输,这是孤子的最主要特征和研究价值所在。2、当光脉冲在双芯耦合器中传输时,取SEC,01THTA2IGN1211KC也就是说,在光纤反常色散区只考虑二阶色散,而忽略损耗、高阶色散和模间色散。3光脉冲在双芯耦合器的第一根光纤中传输光脉冲在双芯耦合器的第二根光纤中传输3、当脉冲在三芯耦合器中传输时,取SEC,01THTA0,32IGN132213221KKC013也就是说,色散中只考虑二阶色散,耦合只存在相邻纤芯间。光脉冲在三芯耦合器的第一根光纤中传输4光脉冲在三芯耦合器的第二根光纤中传输光脉冲在三芯耦合器的第三根光纤中传输4、结果分析光脉冲在光纤耦合器中传输时,随着传输距离的增加,脉冲宽度也不断地展宽,这是由于光纤中色散和非线性效应共同作用的结果。需要注意的是,计算精度不够时会将这一结果放大,甚至会很7快扩散开,并因MATLAB中快速傅立叶变换的周期性边界条件,致使能量从时间窗口的一边跑到另一边,而导致错误结果,因此选择适当的精度也是很关键的。8从图中可以看出,输入的能量在耦合器的光纤中传递,保持总能量不变,光脉冲形状由于色散和非线性效应而展宽,但仍可持续传递。在三芯耦合器中,中间纤芯的能量变化周期约是边上两根纤芯能量变化周期的2倍。四、MATLAB仿真程序PARTONEWSPACEMBYTHOMASEMURPHYFUNCTIONWWSPACET,NTTHISFUNCTIONCONSTRUCTSALINEARLYSPACEDVECTOROFANGULARFREQUENCIESTHATCORRESPONDTOTHEPOINTSINANFFTSPECTRUMTHESECONDHALFOFTHEVECTORISALIASEDTONEGATIVEFREQUENCIESUSAGEWWSPACETVWWSPACET,NT5TVVECTOROFLINEARLYSPACEDTIMUESTSCALARREPRESENTINGTHEPERIODICITYOFTHETIMESEQUENCENTNUMBEROFPOINTSINTIMESEQUENCESHOULDONLYBEPROVIDEDIFFIRSTARGUMENTISSCALAROUTPUTWVECTOROFANGULARFREQUENCIESAUTHORTHOMASEMURPHYTEMUMDEDUIFNARGINPI/DTWKVWKV2PI/DTPARTTWOJ_F_PENGMBYJFPENGUSESIMULATIONOFOPTICALSOLITONTRANSMISSIONINOPTICALFIBERCOUPLERAUTHORJFPENGJINGYUJIAFU163COMDATE20080331COPYRIGHTWSPACEMTHOMASEMURPHY,J_F_PENGMJFPENG理论来源施娟,侯韶华光纤耦合器中光孤子传输仿真研究J详见参考文献CLEARCLCH迭代时选择的步长Z位置坐标T时间窗口1010NT时间区域个数(2M个)DT选择的时间步长T时间坐标W傅里叶域频率GN关于W的函数A单根光纤中模场的慢变振幅AT1双芯耦合器第一个纤芯中模场的归一化慢变振幅AT2双芯耦合器第二个纤芯中模场的归一化慢变振幅ATH1三芯耦合器第一个纤芯中模场的归一化慢变振幅ATH2三芯耦合器第二个纤芯中模场的归一化慢变振幅ATH3三芯耦合器第三个纤芯中模场的归一化慢变振幅H0001Z0H106T20NT212DTT/NTT1NT NT1/2DTWWSPACET,NT上式引用了THOMASEMURPHY的程序WSPACEM计算圆频率GNIW 2/2输入初始脉冲A1,SECHTAT11,SECHTAT21,ZEROSLENGTHT,1ATH11,SECHTATH21,ZEROSLENGTHT,1ATH31,ZEROSLENGTHT,1对称分步傅里叶变换法求解光脉冲在N芯耦合器中的耦合模方程组的数值解的迭代过程由于这里只考虑二阶色散,故可以不考虑FFT和IFFT的交换问题(本程序有考虑到)FORK11LENGTHZ1孤子(光脉冲)在单根光纤中的传输演化过程ANFFTEXPGNH/2IFFTAK,APANEXPIHABSAK,2AK1,FFTEXPGNH/2IFFTAPFORL12APANEXPIH/2ABSAK,2ABSAK1,2AK1,FFTEXPGNH/2IFFTAPEND孤子(光脉冲)在双芯耦合器中的传输演化过程AT1NFFTEXPGNH/2IFFTAT1K,IH/2IFFTAT2K,AT1PAT1NEXPIHABSAT1K,2AT2NFFTEXP