CN101714085A - 光纤通信系统中光脉冲传输特性的matlab计算程序 - Google Patents
光纤通信系统中光脉冲传输特性的matlab计算程序 Download PDFInfo
- Publication number
- CN101714085A CN101714085A CN200910203182A CN200910203182A CN101714085A CN 101714085 A CN101714085 A CN 101714085A CN 200910203182 A CN200910203182 A CN 200910203182A CN 200910203182 A CN200910203182 A CN 200910203182A CN 101714085 A CN101714085 A CN 101714085A
- Authority
- CN
- China
- Prior art keywords
- optical fiber
- initial
- light pulse
- transmission
- fiber communication
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Optical Modulation, Optical Deflection, Nonlinear Optics, Optical Demodulation, Optical Logic Elements (AREA)
Abstract
本发明涉及光纤通信系统中光脉冲传输特性的MATLAB计算程序,属于光纤通信技术领域。本发明提供的计算程序,可通过调整不同参数(输入脉冲的初始时域波形、初始幅值、初始脉冲宽度和初始啁啾C,计算所用时域窗口大小、采样点数,不同光纤的色散参量、衰减系数、传输距离等等)计算光脉冲线性传输、非线性传输和高阶非线性传输等各种情况下的光脉冲传输特性。本发明提出的计算程序简洁易懂,操作方便,功能强大,计算速度快,作图简便,实用性强,可根据实际情况修改相关参数,为光纤通信系统优化及该领域的深入研究提供了有力支持,将推动本学科的快速发展。
Description
所属技术领域
本发明涉及一种光纤通信系统中光脉冲传输特性的MATLAB计算程序,属于光纤通信技术领域。
背景技术
美国科学家T.H.Maiman在1960年研制成功了世界上第一台激光器(红宝石激光器)以来,激光技术取得了迅速的发展,已经广泛应用于工农业生产、生命科学、信息科学、强场物理、“快点火”受控核聚变、探索极端物理条件下物质行为等重大科学研究领域,具有广泛的应用前景,在整个经济及社会发展中发挥出愈来愈重要的作用。得益于激光技术的发展,光纤技术得到了巨大的进步。至1970年,人们已经将光纤损耗降低到约20dB/km;到1979年,将1.55μm波长的光纤损耗降低到约0.2dB/km。低损耗光纤的获得,不仅导致了光纤通信领域的革命,而且也导致了非线性光纤光学新领域的出现。随后,掺铒光纤放大器得到了迅速发展,有效补偿了1.55μm波长区的光纤损耗,有力地推动了线性和非线性光纤光学的快速发展,激光脉冲在光纤通信系统中的线性和非线性传输特性研究得到了国内外科技工作者的普遍关注,特别是由广义非线性薛定谔方程
所描述的脉冲传输问题出现在物理学、等离子体、信息科学、生命科学等现代科学的各个分支中,并为相关学科中众多的类似问题的解决做出了重大贡献。然而,方程(1)是非线性偏微分方程,在一般情况下不适于解析求解,除非是在能使用逆散射方法的某些特殊情况下才有可能,这限制了该方程在相关学科的应用和阻碍了相应学科的发展。鉴于此,采用数值计算方法就不同学科中相关问题来求解方程(1)成为不同学科发展的重要研究课题。
发明内容
本发明针对方程(1)在光纤通信系统中难于解析求解的情况和计算机技术的快速发展,在Windows操作系统的MATLAB环境中按照分步傅里叶方法设计和撰写了光纤通信系统中光脉冲传输特性的MATLAB数值计算程序来解决上述问题。本发明提出的计算程序简洁易懂,操作方便,功能强大,计算速度快,作图简便,实用性强,可根据实际情况修改相关参数,为光纤通信系统优化及该领域的深入研究提供了有力支持,将推动本学科的快速发展。本发明解决其技术问题所采用的技术方案是:
为了方便利用方程(1)研究光纤通信系统中光脉冲的传输特性,本发明引入ξ=z/LD, u=NU,将方程(1)归一化,得
式中,u为光脉冲慢变包络振幅,ξ为归一化传输距离,τ为归一化时间参量;右边第一项为色散项,βn为各阶群速度色散系数,n表示色散阶数,第二项至第四项表示脉冲传输中的非线性项,第三项和第四项为高阶非线性项,其第二项为自相位调制效应;第三项为脉冲自变陡效应,对高阶孤子的影响显著,s=1/ω0T0,ω0为脉冲中心角频率;第四项为拉曼散射效应,τR=TR/T0,TR为脉冲内拉曼散射系数,与拉曼增益的斜率有关,通常为3fs;最后一项表示光纤损耗项,α为光纤损耗系数,LD是色散长度,T0是脉冲的半宽度。
本发明采用分步傅里叶方法(应用了快速傅里叶变换),按照方程(2)设计和撰写了光纤通信系统中光脉冲传输特性的MATLAB计算程序。一般来说,沿光纤的长度方向,色散和非线性是同时作用的。分步傅里叶方法通过假定在传输过程中,光场每通过一小段距离h,色散和非线性效应可分别作用,得到近似结果。更准确地说,从z到z+h的传输过程中分两步进行。第一步,仅有色散作用,方程(2)中的非线性作用为零;第二步,仅有非线性作用,方程(2)中的色散作用为零。分步傅里叶方法已广泛应用于各种光学领域,包括:大气中的光传输,折射率梯度光纤,半导体激光器,非稳腔及波导耦合器等。虽然用此方法运算速度快,相对简洁,但需要小心选择z和T的步长,以保证精度要求。该算法隐含着周期性边界条件,有时会造成到达窗口一边的能量从窗口的另一边自动进入的问题,通常需要使用一种“吸收窗口”,人为地吸收掉辐射到窗口边界上的能量。
本发明的有益效果是:本发明提出的计算程序简洁易懂,操作方便,功能强大,计算速度快,作图简便,实用性强,可根据实际情况快速修改相关参数,为光纤通信系统优化及该领域的深入研究提供了有力支持,将推动本学科的快速发展。
具体实施方式
光纤通信系统中光脉冲传输特性的MATLAB计算程序运行前,可根据实际情况修改相关参数,如:输入脉冲的初始时域波形、初始幅值、初始脉冲宽度和初始啁啾C,计算所用时域窗口大小、采样点数,不同光纤的色散参量、衰减系数、传输距离等等。
光脉冲线性传输情况:当输入脉冲的初始脉冲宽度在ps附近和光功率较小时(线性传输),可调节相关参数,使计算程序仅按照方程(2)左边第一项和右边第一项(色散项)运行,即按照运行,其中高阶色散可计算至六阶色散。当光脉冲传输距离较长时,可同时考虑方程(2)右边第五项(光纤损耗项)。
光脉冲非线性传输情况:当输入脉冲的初始脉冲宽度在ps附近和光功率较大时(非线性传输),可调节相关参数,使计算程序仅按照方程(2)左边第一项、右边第一项(色散项)和右边第二项(自相位调制项)运行,即按照运行。当光脉冲传输距离较长时,可同时考虑方程(2)右边第五项(光纤损耗项)。
光脉冲高阶非线性传输情况:当输入脉冲的初始脉冲宽度在fs附近和光功率很大时(高阶非线性传输),可调节相关参数,使计算程序仅按照方程(2)左边第一项、右边第一项(色散项)、右边第二项(自相位调制项)、右边第三项(高阶非线性项:自陡项)和右边第四项(高阶非线性项:脉冲内受激拉曼散射项)运行,即按照运行。当光脉冲传输距离较长时,可同时考虑方程(2)右边第五项(光纤损耗项)。
光脉冲分布拉曼放大传输情况:当输入脉冲在分布拉曼放大传输情况下,可调节相关参数,使计算程序仅按照方程运行,光纤损耗的拉曼增益系数α′R=αsexp[-αp(L′-ξ)LD]为归一到色散长度的参量,αp是在泵浦波长的光纤损耗系数,αs是在信号波长的光纤损耗系数,L′为归一化光纤长度。
总之,本发明提供了光纤通信系统中光脉冲传输特性的MATLAB计算程序,可通过调整不同参数计算各种情况下的光脉冲传输特性。我们应用该程序已经部分研究了光脉冲传输特性,取得了部分研究成果并发表在国内外重要刊物上。可见,本发明为光纤通信系统优化及该领域的深入研究提供了有力支持,将有力地推动本学科快速发展。
Claims (1)
1.本发明是光纤通信系统中光脉冲传输特性的MATLAB计算程序,属于光纤通信技术领域。其特征在于:
(1)本发明按照方程(2),在Windows操作系统的MATLAB环境中采用分步傅里叶方法设计和撰写了光纤通信系统中光脉冲传输特性的MATLAB数值计算程序。
(2)本发明提供的计算程序,可通过调整不同参数(输入脉冲的初始时域波形、初始幅值、初始脉冲宽度和初始啁啾C,计算所用时域窗口大小、采样点数,不同光纤的色散参量、衰减系数、传输距离等等)计算光脉冲线性传输、非线性传输和高阶非线性传输等各种情况下的光脉冲传输特性。
(3)本发明提供的计算程序简洁易懂,操作方便,功能强大,计算速度快,作图简便,实用性强,可根据实际情况修改相关参数,为光纤通信系统优化及该领域的深入研究提供了有力支持,将推动本学科的快速发展。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910203182A CN101714085A (zh) | 2009-06-03 | 2009-06-03 | 光纤通信系统中光脉冲传输特性的matlab计算程序 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910203182A CN101714085A (zh) | 2009-06-03 | 2009-06-03 | 光纤通信系统中光脉冲传输特性的matlab计算程序 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101714085A true CN101714085A (zh) | 2010-05-26 |
Family
ID=42417753
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200910203182A Pending CN101714085A (zh) | 2009-06-03 | 2009-06-03 | 光纤通信系统中光脉冲传输特性的matlab计算程序 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101714085A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102884736A (zh) * | 2012-06-28 | 2013-01-16 | 华为技术有限公司 | 一种光纤链路非线性传输损伤的计算方法和装置 |
CN103150429A (zh) * | 2013-02-28 | 2013-06-12 | 北京工业大学 | 全正色散全光纤激光器锁模动力学matlab计算方法 |
CN104142908A (zh) * | 2013-05-07 | 2014-11-12 | 中国人民解放军海军航空工程学院 | 电波传播抛物方程分步傅里叶变换解的上边界处理方法 |
US9698493B2 (en) | 2012-05-29 | 2017-07-04 | Huawei Technologies Co., Ltd. | Dual-polarized antenna radiating element and base station antenna |
-
2009
- 2009-06-03 CN CN200910203182A patent/CN101714085A/zh active Pending
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9698493B2 (en) | 2012-05-29 | 2017-07-04 | Huawei Technologies Co., Ltd. | Dual-polarized antenna radiating element and base station antenna |
CN102884736A (zh) * | 2012-06-28 | 2013-01-16 | 华为技术有限公司 | 一种光纤链路非线性传输损伤的计算方法和装置 |
WO2014000229A1 (zh) * | 2012-06-28 | 2014-01-03 | 华为技术有限公司 | 一种光纤链路非线性传输损伤的计算方法和装置 |
CN102884736B (zh) * | 2012-06-28 | 2015-04-08 | 华为技术有限公司 | 一种光纤链路非线性传输损伤的计算方法和装置 |
US9525481B2 (en) | 2012-06-28 | 2016-12-20 | Huawei Technologies Co., Ltd. | Method and apparatus for calculating nonlinear transmission impairment of optical fiber link |
CN103150429A (zh) * | 2013-02-28 | 2013-06-12 | 北京工业大学 | 全正色散全光纤激光器锁模动力学matlab计算方法 |
CN103150429B (zh) * | 2013-02-28 | 2016-03-30 | 北京工业大学 | 全正色散全光纤激光器锁模动力学matlab计算方法 |
CN104142908A (zh) * | 2013-05-07 | 2014-11-12 | 中国人民解放军海军航空工程学院 | 电波传播抛物方程分步傅里叶变换解的上边界处理方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Bai et al. | Recent advances in Brillouin optical time domain reflectometry | |
Mohammadi et al. | Simulation of a transmission system to compensate dispersion in an optical fiber by chirp gratings | |
Khan et al. | Determination of the best apodization function and grating length of linearly chirped fiber Bragg grating for dispersion compensation. | |
CN101714085A (zh) | 光纤通信系统中光脉冲传输特性的matlab计算程序 | |
CN101762290A (zh) | 基于分布式拉曼放大的布里渊光时域分析系统 | |
CN104111086A (zh) | 基于低布里渊散射阈值传感光纤的光时域反射仪的装置与方法 | |
CN104272072A (zh) | 波形重构装置、波形重构系统以及波形重构方法 | |
CN102539015B (zh) | 一种分布式温度传感光纤 | |
CN102564481A (zh) | 提高分布式光纤布里渊传感器信噪比的方法和装置 | |
CN202041323U (zh) | 一种分布式光纤多参量传感器 | |
Burdin et al. | Simulation of a long-haul fiber optic link with a two-mode optical fiber | |
CN204718706U (zh) | 一种高灵敏度光纤测温系统 | |
Chakkour et al. | Design and study of EDFA-WDM optical transmission system using FBG at 10 Gbits/s chromatic despersion compensation effects | |
Rani et al. | Exact solitary wave solution for higher order nonlinear Schrodinger equation using He’s variational iteration method | |
Sharma et al. | Optical performance of tellurite glass microstructured optical fiber for slow-light generation assisted by stimulated Brillouin scattering | |
CN206496737U (zh) | 一种温度增敏型分布式布里渊光纤传感器 | |
Azizan et al. | Simple distributed Brillouin scattering modeling for temperature and strain sensing | |
CN101666954B (zh) | 拉曼放大作用下的高速光孤子传输系统及孤子光谱的计算方法 | |
Andalib et al. | Analytical investigation and evaluation of pulse broadening factor propagating through nonlinear optical fibers (traditional and optimum dispersion compensated fibers) | |
EP2993849A1 (en) | Preprocessing and detection for nonlinear Fourier transform based optical receiver | |
RU221092U1 (ru) | Волоконный усилитель на основе неоднородного по длине активного световода для усиления и временного сжатия частотно-модулированных импульсов пико- и фемтосекундной длительности | |
Kah et al. | Fast Solution for the Nonlinear Schrodinger Equation in Optical Fibers by the Reduced Basis Method | |
Lachina et al. | FEMTOSECOND ERBIUM NONLINEAR FIBER LASER SYSTEM | |
Li et al. | Influence of doping on stimulated Brillouin scattering slow light in few-mode fibers | |
Šajgalíková et al. | Numerical modeling of degradation mechanisms in single mode optical fibre using finite difference method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20100526 |