基于GPU的光纖非線性薛定諤方程的仿真與計算
發(fā)布時間:2020-10-02 06:37
光纖中信號傳播的仿真在光傳輸系統的研究發(fā)展中起著重要作用。研究信號在光纖中的傳輸特性主要是轉化為數學模型去數值求解光纖非線性薛定諤方程。數值方法的性能主要考慮誤差精度和和計算效率,目前求解單模非線性傳輸方程的數值方法已較為完善,定步長方法有全局誤差四階精確的龍格庫塔交互圖片法,為進—步優(yōu)化方法性能,衍生出大量變步長計算準則,使用較為廣泛的是局部誤差法。多模非線性傳輸方程的計算相對單模更為復雜,目前只采用了精度較低的定步長方法,因此很有必要優(yōu)化多模變步長數值方法。此外,在長距離多模光纖通信中,尤其非線性作用占主導時,仿真時間往往長達幾天,提高計算效率已成為重中之重。對數值方法的優(yōu)化更多是在提高數值精度,要大規(guī)模加速計算,很有必要利用GPU并行計算來加速求解過程。本次研究工作的主要成果如下:(1)提出三種求解多模光纖非線性傳輸方程的誤差估計準則—max,sum,ave準則,將多模誤差向量轉換為誤差標量,基于對稱分步傅里葉的局部誤差法實現多模傳輸自適應步長統一變化。通過仿真單個高斯脈沖在漸變折射率多模光纖中的傳輸,驗證了定變步長方法在不同準則下局部誤差與全局誤差的性能。實驗結果表明三種準則的變步長算法都具有收斂性,且利用sum準則計算局部誤差控制步長變化,在相同計算量的情況下能得到更高的數值精度,相同全局誤差的情況下計算量相對更少,對進一步提高多模非線性傳輸方程的計算效率有參考意義。(2)基于GPU并行計算多模非線性傳輸方程取得顯著加速效果。當GPU利用率接近最大值時,加速效果最為明顯。實驗仿真多個高斯脈沖入射到多模光纖中,通過調節(jié)脈沖傳輸個數來使得GPU利用率達到90%,研究表明單精度條件下,傳輸15個模式100個脈沖,相比CPU的計算時間,GPU的加速比可達一百多倍,計算效率得到大幅提升。此外,分別在單雙精度條件下,仿真了多模多脈沖數值算法的誤差精度與傳輸步數及模式數的關系,結果表明,在不同模式數下,隨著傳輸步數的增大,雙精度條件下的誤差精度會更高。
【學位單位】:北京交通大學
【學位級別】:碩士
【學位年份】:2019
【中圖分類】:O241.82;TN929.11
【部分圖文】:
圖2-1邋RK4IP計算步驟圖解逡逑Fig邋2-1邋Calculation邋step邋diagram邋ofRK4IP逡逑
jV邋=邋211,時域窗口為7\vin=l0A,則采樣頻率力=Twm邋/邋#。設置步長A=l00W,分別用逡逑定步長方法S3FM、RK4IP及變步長RK43IP方法進行數值計算,將傳輸前后脈沖逡逑的時域波形及頻域頻譜放在一起進行對比,如圖2-2所示,從二維角度觀察,一段逡逑時間后,傳輸前后的波形及頻譜均重疊在一起,即滿足時域及頻譜無畸變穩(wěn)定傳逡逑輸。從直觀角度表明了所用數值方法的正確性。逡逑3|邐邐邐■——邐—[邐0.08.邐邐邋—[逡逑傳輸前頻ii"!邐+傳輸前波形逡逑2.5.邐:邐1-^M-邐0邋07邋■邐::卜傳輸后逡逑0.06'邋:逡逑2邋■逡逑^邐^邋0.05邋■逡逑M1-5'邐if004邋■逡逑w9邐,丨邐15邋0.03-邐-逡逑1邋.邋■逡逑0_I__邐4_I_逡逑-1邐-0.5邐0邐0.5邐1邐-5邐0邐5逡逑頻率邐x1013邐時間/ps邐x10-1。逡逑a)邐b)逡逑圖2-2邋—階光孤子傳輸前后時域及頻域對比圖a)頻譜對比b)時域對比逡逑Fig邋2-2邋Time-domain邋and邋frequency-domain邋comparison邋diagram邋before邋and邋after邋first-order逡逑optical邋soliton邋transmission邋a)邋Spectrum邋comparison邋b)邋Time邋domain邋comparison逡逑2.4.2二階光孤子脈沖傳輸仿真逡逑為進一步驗證數值算法的性能
同距離且初始步長一樣的情況下傳輸步數也一致,可以更方便對比局部誤差與傳逡逑輸距離的關系。經計算,RK43IP法傳輸一個孤子周期需傳輸390步,所以初始步逡逑長設置為0.506/390m,如圖2-4所示,a)圖為一個孤子周期內步長與傳輸距離的關逡逑系,h-RK43IP為RK43IP算法計算的步長,h-average為平均傳輸步長,實際也是逡逑固定采用初始步長傳輸的情況?梢园l(fā)現,在開始傳輸的極短距離內步長迅速增逡逑大,變化趨勢接近一條垂直線段,然后步長變化呈現為一正向拋物線形狀。同時逡逑因兩種步長對應傳輸步數是一致的,所以兩條步長曲線與兩條坐標軸圍成的面積逡逑是相同的。b)圖展示的是三個孤子周期內步長的變化情況,可以發(fā)現變化趨勢同第逡逑一個孤子周期內形狀,呈現穩(wěn)定的正向拋物線變化趨勢,由此可見初始條件固定逡逑后
本文編號:2832053
【學位單位】:北京交通大學
【學位級別】:碩士
【學位年份】:2019
【中圖分類】:O241.82;TN929.11
【部分圖文】:
圖2-1邋RK4IP計算步驟圖解逡逑Fig邋2-1邋Calculation邋step邋diagram邋ofRK4IP逡逑
jV邋=邋211,時域窗口為7\vin=l0A,則采樣頻率力=Twm邋/邋#。設置步長A=l00W,分別用逡逑定步長方法S3FM、RK4IP及變步長RK43IP方法進行數值計算,將傳輸前后脈沖逡逑的時域波形及頻域頻譜放在一起進行對比,如圖2-2所示,從二維角度觀察,一段逡逑時間后,傳輸前后的波形及頻譜均重疊在一起,即滿足時域及頻譜無畸變穩(wěn)定傳逡逑輸。從直觀角度表明了所用數值方法的正確性。逡逑3|邐邐邐■——邐—[邐0.08.邐邐邋—[逡逑傳輸前頻ii"!邐+傳輸前波形逡逑2.5.邐:邐1-^M-邐0邋07邋■邐::卜傳輸后逡逑0.06'邋:逡逑2邋■逡逑^邐^邋0.05邋■逡逑M1-5'邐if004邋■逡逑w9邐,丨邐15邋0.03-邐-逡逑1邋.邋■逡逑0_I__邐4_I_逡逑-1邐-0.5邐0邐0.5邐1邐-5邐0邐5逡逑頻率邐x1013邐時間/ps邐x10-1。逡逑a)邐b)逡逑圖2-2邋—階光孤子傳輸前后時域及頻域對比圖a)頻譜對比b)時域對比逡逑Fig邋2-2邋Time-domain邋and邋frequency-domain邋comparison邋diagram邋before邋and邋after邋first-order逡逑optical邋soliton邋transmission邋a)邋Spectrum邋comparison邋b)邋Time邋domain邋comparison逡逑2.4.2二階光孤子脈沖傳輸仿真逡逑為進一步驗證數值算法的性能
同距離且初始步長一樣的情況下傳輸步數也一致,可以更方便對比局部誤差與傳逡逑輸距離的關系。經計算,RK43IP法傳輸一個孤子周期需傳輸390步,所以初始步逡逑長設置為0.506/390m,如圖2-4所示,a)圖為一個孤子周期內步長與傳輸距離的關逡逑系,h-RK43IP為RK43IP算法計算的步長,h-average為平均傳輸步長,實際也是逡逑固定采用初始步長傳輸的情況?梢园l(fā)現,在開始傳輸的極短距離內步長迅速增逡逑大,變化趨勢接近一條垂直線段,然后步長變化呈現為一正向拋物線形狀。同時逡逑因兩種步長對應傳輸步數是一致的,所以兩條步長曲線與兩條坐標軸圍成的面積逡逑是相同的。b)圖展示的是三個孤子周期內步長的變化情況,可以發(fā)現變化趨勢同第逡逑一個孤子周期內形狀,呈現穩(wěn)定的正向拋物線變化趨勢,由此可見初始條件固定逡逑后
【參考文獻】
相關期刊論文 前3條
1 鐘鳴宇;劉東風;胡長俊;;高階非線性薛定諤方程的分步小波方法[J];光子學報;2012年08期
2 李瑩;崔慶豐;;基于分步傅里葉變換法對非線性薛定諤方程的數值仿真[J];長春理工大學學報(自然科學版);2011年01期
3 邱橋飛;張冠茂;張曉萍;;快速數值差分遞推改進算法在超寬連續(xù)譜模擬中的應用研究[J];光子學報;2011年01期
相關博士學位論文 前3條
1 吳云虎;光纖非線性薛定諤方程數值仿真方法優(yōu)化研究[D];華中師范大學;2015年
2 周勤;光孤子傳輸特性的解析研究[D];武漢大學;2014年
3 韓佳巍;大容量模分復用光傳輸系統的若干關鍵技術研究[D];北京郵電大學;2013年
相關碩士學位論文 前5條
1 鄧春雨;光纖中基于受激布里淵散射的模間色散補償及傳感[D];蘭州理工大學;2017年
2 蔡躍進;耦合飛秒光孤子的相互作用研究[D];聊城大學;2017年
3 劉志飛;多模光纖高速傳輸關鍵技術研究[D];北京交通大學;2016年
4 林睿哲;光孤子相互作用的解析研究與數值模擬[D];哈爾濱工業(yè)大學;2015年
5 孟佳;一類非線性薛定諤方程的數值解法[D];河南大學;2015年
本文編號:2832053
本文鏈接:http://www.sikaile.net/kejilunwen/xinxigongchenglunwen/2832053.html