張曉光

(北京郵電大學信息光子學與光通信國家重點實驗室，北京 100876)

自聚焦多模光纖模式色散的幾何光學與波動光學處理方法

張曉光

(北京郵電大學信息光子學與光通信國家重點實驗室，北京 100876)

自聚焦多模光纖相比階躍型多模光纖，模式色散要小兩個數(shù)量級，現(xiàn)在廣泛用于局域網(wǎng)。對于自聚焦光纖模式色散的處理方法主要是基于原用于量子力學的WKBJ近似方法。本文從幾何光學和波動光學兩個角度分析處理自聚焦光纖的模式色散，得到了與WKBJ方法完全相同的結果。在幾何光學的處理方法中，利用費馬原理作指導，引入遠軸光線與近軸光線概念的區(qū)別，用微擾的近似方法解光線方程，最后得到自聚焦光纖的模式色散。在波動光學的處理方法中，基于波動方程，將光纖模式近似成厄米高斯模式，依然采用近軸近似和遠軸近似的概念，通過計算傳播常數(shù)得到模式色散。兩種處理方法得到的模式色散結果一致，且處理過程物理意義清楚，非常適合教學時采用。

導波光學；光波導模式；自聚焦光纖；模式色散

0 導言

光波導理論是光纖通信系統(tǒng)的基礎理論之一，在光波導理論中，光波導模式理論和模式色散是重要的內(nèi)容。對于階躍折射率分布的多模光纖，其纖芯內(nèi)的折射率分布是均勻的，因此計算模式色散也是簡單的，不需要太多的數(shù)學理論。但是由于階躍折射率多模光纖的模式色散過大，在目前的光纖通信系統(tǒng)中已經(jīng)鮮有應用。在骨干網(wǎng)和城域網(wǎng)中，由于要求的通信容量很大，一般應用沒有模式色散的單模光纖。在局域網(wǎng)中，由于折射率漸變的自聚焦光纖模式色散相對來講比較小，應用得比較多。由于自聚焦光纖纖芯內(nèi)折射率分布不是均勻的，計算其模式色散，一般采用原用于量子力學的WKBJ近似方法(即所謂的相位積分法)求解[1-3]，數(shù)學比較復雜，計算過程物理意義也不十分明顯，一直是教學上的一個難點。

筆者在北京郵電大學講授光通信的物理基礎課時，開始研究自聚焦光纖模式色散的教學處理方法。由于講課對象是剛學完大學物理課程，還沒有學電磁場與電磁波或者電動力學課程的學生，用傳統(tǒng)全波動觀點講授光纖模式和利用WKBJ法處理模式色散的講授方式會讓學生難以接受。因此筆者試圖利用全幾何光學的處理方法來求解自聚焦光纖模式色散，發(fā)現(xiàn)結合光程、費馬原理的概念，利用光線方程求解自聚焦光纖模式色散，數(shù)學非常簡單，幾乎每一個處理過程其物理意義都很清楚，利于學生理解。其中引入了聚焦特性分析以及近軸光線與遠軸光線區(qū)別等概念，將聚焦特性、近軸光線、遠軸光線等物理概念與數(shù)學上所取近似程度緊密結合起來，建立了一幅清晰的物理圖像，說清楚了以下幾個物理問題：(1)近軸光線近似與費馬原理預計的等光程性等價，不能用近軸近似結果來計算自聚焦光纖的模式色散；(2)自聚焦光纖的模式色散來源于遠軸光線；(3)由于聚焦特性的不同，自聚焦光纖的模式色散比階躍折射率多模光纖小兩個數(shù)量級。另外從授課方式來看，從費馬原理的基礎理論出發(fā)，最后解決的是自聚焦光纖模式色散的工程問題，符合從理論到應用的典型學習與認知過程。這樣無論是宏觀上從理論到應用的安排，還是教學中物理概念與數(shù)學緊密結合的處理，構成一個較全面的教學模塊，有推廣價值。如果在教學實踐中，在講授完費馬原理和光線彎曲后，介紹一下太陽下山、海市蜃樓等光線彎曲的自然現(xiàn)象，最后過度為自聚焦光纖模式色散，構成一個從理論講授到自然現(xiàn)象的解釋，再到工程應用的完整架構，那就更加完美了。

為了給自聚焦光纖模式色散一個更完整的教學處理方案，筆者研究了文獻[4]，其中將自聚焦光纖模式處理成厄米高斯模式。在此基礎之上，引入近軸近似和遠軸近似方法，很方便地給出了自聚焦光纖模式色散的正確結果。這樣，本文針對自聚焦光纖模式色散給出了一個完整的幾何光學與波動光學的教學處理方案。

1 自聚焦光纖模式色散的幾何光學處理方法

1.1 費馬原理與光線彎曲

在中學物理的學習中，我們知道幾何光學建立在光線直線傳播、反射定律和折射定律這3大定律之上。大學物理的光學課，幾何光學是建立在費馬原理的基礎之上。也就是說，上述3大定律都可以看成是費馬原理的推論。

費馬原理表述成[5]：兩點之間的真實光線是所走光程相對于路徑的變分為平穩(wěn)的那一條，數(shù)學描述如式(1)和圖1。

(1)

圖1 路徑L與L′是從A點到B點的可能路徑，其中真實光線路徑ALB是光程相比相鄰路徑AL′B取平穩(wěn)值的那一條

這樣，光線的直線傳播只是費馬原理在均勻媒質(zhì)(折射率是均勻的)中的推論，對更普遍的情形，光學媒質(zhì)折射率是不均勻的，因此根據(jù)費馬原理，光線可以是彎曲的。以太陽落山為例，由于大氣層接近地表的空氣密度大，折射率高，遠離地表，折射率逐漸減小，太陽落山時，光線是彎曲傳輸?shù)模孛嫔系娜丝吹教栂律綍r，其實“太陽早就下山了”。海市蜃樓也是光線彎曲很好的例子。

光程變分取平穩(wěn)值，其中取極小值、極大值、拐點的情景，其真實光線是附近虛擬光線中取變分后的唯一一條，而取光程常數(shù)值的情景，得到的真實光線是一組光纖簇，光線簇中的每一根光線光程都相等，即所謂的等光程性。成像光學是光程取常數(shù)的典型例子。比如人眼成像(圖2)，物點A發(fā)出的傍軸光線經(jīng)眼球聚焦在視網(wǎng)膜上，形成像點，從物點A到像點B的所有傍軸光線的光程都應該相等。由光線所經(jīng)歷時間與所走光程之間的關系t=光程/c(c是真空中的光速)，光程相等等價于時間相等，亦即從物點A到像點B的所有傍軸光線所需時間都相等。設想一下，如果這個時間不等，人眼看到的圖像將是不穩(wěn)定的。

圖2 人眼成像

應用費馬原理，還可以得到彎曲光線尋跡的定量描述(即光線傳輸所遵循的方程)——光線方程(公式(2))。這里假定z軸為傍軸光線的主傳播方向，而傳輸光線的媒質(zhì)假定只在x方向變化[5]

(2)

(3)

從光線方程(2)可以看出，只要給出媒質(zhì)的折射率分布，就可以通過光線方程得到光線軌跡。

1.2 光纖模式與模式色散

如圖3所示，一般來說，光纖大致分階躍折射率分布光纖和漸變折射率分布光纖。階躍折射率光纖又根據(jù)纖芯的大小分單模光纖和階躍折射率多模光纖[6]，如圖3(a)和(b)所示。漸變折射率分布光纖也是多模的，如果將折射率分布做成如公式(4)的拋物線分布，則光線在其中可以周期性地聚焦，這種光纖也叫自聚焦光纖。

圖3 (a) 階躍折射率分布多模光纖； (b) 階躍折射率單模光纖； (c) 漸變折射率分布多模光纖

(4)

所謂光纖模式是光纖中能夠存在的穩(wěn)定光場形態(tài)。如果用幾何光學來描述光纖模式，可以粗略地用纖芯中不同角度全反射的不同光線來表示。如圖4(a)所示，接近臨界角θc入射的光線表示最高階模式，入射角大一些的光線代表低階模式，沿軸線傳播的光線代表基模。

圖4 (a) 光纖模式色散示意圖； (b) 模式色散導致的誤碼

一個光脈沖從光纖左端射入，由于不同模式光線沿軸線方向的傳播速度不同，造成在不同模式中的光脈沖到達光纖右端的時間不同，最后合成的光脈沖將是一個展寬的脈沖，這就是光纖中的模式色散現(xiàn)象。模式色散會導致光纖通信產(chǎn)生誤碼，如圖4(b)所示，假設待傳輸?shù)拇a流是“110101”，有光脈沖時記為“1”碼，沒有光脈沖時記為“0”碼。如果模式色散造成光脈沖展寬，展寬到一定程度，“0”碼處也有光功率，極易判斷為“1”碼，這就是誤碼。

對于階躍折射率分布光纖，其最大模式(接近臨界角θc的光線)傳輸L距離所需時間為

(5)

基模傳輸L距離所需時間為

(6)

模式色散定義為光傳輸單位距離產(chǎn)生的最大時間差(也叫時延差)為

(7)

將階躍折射率多模光纖的典型數(shù)值代入公式(7)，n1=1.500，n2=1.489，則Δ≈0.011，模式色散Δτ=55ns/km。如果對應于某傳輸碼率，傳輸信號的比特周期至少要求與模式色散相當，則傳輸20km，傳輸碼率被限制在1Mbit/s以下，這是一個非常低的傳輸容量，因此階躍折射率多模光纖不適合目前要求的高速率光纖通信系統(tǒng)。

1.3 自聚焦光纖模式色散的幾何光學處理方法

關于自聚焦光纖模式色散的幾何光學處理，是筆者的一項教學成果[7]。由于本文的目的是給出一個完整的教學處理方案，因此這里還是要給出文獻[7]的主要分析與結論。

公式(4)指出，自聚焦光纖具有拋物線型的折射率分布，中心的折射率高，折射率沿著半徑逐漸減少。因此進入自聚焦光纖的光線趨向于向軸線彎曲，即光線圍繞著軸線擺動。由于對稱性，我們只討論子午面內(nèi)光線的軌跡，在子午面z-x平面內(nèi)，折射率可以表示成

(8)

其中α2=2Δ/a2。

將上式代入光線方程(2)，并考慮傍軸光線近似，得

(9)

(10)

顯然，公式(10)是一個典型的簡諧振動方程，其通解為

(11)

其中C和φ是積分常數(shù)，分別代表光線擺動的幅度和相位。式(11)表明，近軸光線是周期性圍繞光軸正弦擺動的，α是擺動的空間頻率?？梢钥闯鏊薪S光線具有相同的空間周期2π/α，即近軸光線是周期性聚焦的，如圖5所示。根據(jù)費馬原理，從軸線A點發(fā)出的所有近軸光線都將聚焦于B點，A點和B點之間的近軸光線是等光程的，不會引起時間差。換句話說，只考慮近軸光線，無法計算出模式色散，模式色散來源于遠軸光線，這就是筆者將自聚焦光纖的傍軸光線分成近軸光線和遠軸光線的原因。

圖5 自聚焦光纖中的近軸光線與遠軸光線

可以預言，由于自聚焦光纖光線的這種聚焦特性，模式色散只來源于遠軸光線，相較于階躍折射率多模光纖光線根本沒有聚焦特性，自聚焦光纖的模式色散顯然比起階躍折射率多模光纖要小很多，這在隨后的討論中將被證實。

在討論遠軸光線時，利用泰勒展開，保留到αx的二階項，公式(9)變成

(12)

下面我們用微擾方法討論式(12)?？紤]遠軸光線仍具有正弦擺動特性，但是空間頻率有所不同

(13)

將式(13)代入式(12)，得到

(14)

其中遠軸光線的空間周期變成2π/σ，不是常數(shù)。因此在計及遠軸光線后，各光線不再可以周期性聚焦，存在著時間差。從式(14)可見，遠軸光線的空間頻率與近軸光線有差別，隨著擺動幅度C的不同而不同。最遠的遠軸光線幅度為a，則遠軸光線的最大空間頻率滿足

(15)

或者

(16)

令最大模式的對應光線軌跡函數(shù)為

x=asinσmaxz

(17)

計算最大模式傳輸L距離后的光程

(18)

將上式被積函數(shù)做泰勒展開，保留到α4a4項，則式(18)變?yōu)?/p>

(19)

式(19)的后一積分項在積分跨越2π/σmax點時為零，則

(20)

由光程差與時間差的關系，自聚焦光纖的模式色散為

(21)

這個結果與WKBJ近似法計算的結果完全一致[1,2]。將n1=1.500和n2=1.489代入式(21)，得Δτmax=0.25ns/km，比階躍折射率多模光纖小兩個數(shù)量級。

總結一下第1節(jié)的主要結論：(1)近軸光線近似與費馬原理預計的等光程性等價，不能用近軸近似結果來計算自聚焦光纖的模式色散； (2)自聚焦光纖的模式色散來源于遠軸光線，利用遠軸近似可以得出自聚焦光纖模式色散的正確結果； (3)由于聚焦特性的不同，自聚焦光纖的模式色散比階躍折射率多模光纖小兩個數(shù)量級。

1.4 階躍折射率多模光纖與自聚焦光纖的比較

前面我們預言，自聚焦光纖具有一定的聚焦特性，而階躍折射率多模光纖根本沒有聚焦特性，因此，自聚焦光纖模式色散小得多。下面我們還是用光程的泰勒展開來分析階躍折射率多模光纖的模式色散。

我們在圖6的階躍折射率多模光纖中選取光線OGF進行分析，假定它的空間周期也為2π/α，則有

(22)

其中n1L是軸線基模光線走過的光程?？梢?，由于階躍折射率多模光纖沒有聚焦特性，高階模與基模光線的相對光程差含有α2a2項。而自聚焦光纖有明顯的聚焦特性，這個差值只含有α4a4。顯然自聚焦光纖的模式色散比階躍折射率多模光纖小兩個數(shù)量級。

圖6 階躍折射率多模光纖的聚焦特性

為了計算模式色散，必須求得最大模式條件，最大模式顯然是對應臨界角的光線。因此

(23)

整理式(23)，得

(24)

將公式(24)代入式(22)，得到

(25)

則階躍折射率多模光纖的模式色散為

(26)

與式(7)是一致的。一般地Δ～10-2，則Δ2～10-4，因此自聚焦光纖的模式色散確實比階躍折射率多模光纖小兩個量級。

2 自聚焦光纖模式色散的波動光學處理方法

在通常的光纖通信與光波導文獻中，處理自聚焦光纖模式色散是利用原用于量子力學的WKBJ法來近似計算的。本文借鑒文獻[4]的處理方法，將自聚焦光纖中的模式按照厄米高斯函數(shù)這種正交函數(shù)展開來處理，也可以計算自聚焦光纖的模式色散。

2.1 自聚焦光纖模式的厄米高斯函數(shù)展開處理

從標量波動方程出發(fā)

(27)

令

E(x,y,z,t)=ψ(x,y,z)e-iωt

(28)

得到亥姆霍茲方程

(29)

對于自聚焦光纖，折射率n表示成

(30)

(31)

設ψ有分離變量解

ψ(x,y,z)=X(x)Y(y)Z(z)

(32)

代入方程(31)，整理得

(33)

如果上式對任意自變量(x,y,z)都成立，等號左邊的三項必然分別等于3個常數(shù)，即

(34)

和

(35)

其中

(36)

方程(35)的解為

Z(z)～e±iβz

(37)

代表沿z軸正向(正號)和負向(負號)的行波，這里取正號，即取沿z軸正向波。

如果令

(38)

則式(34)的兩個方程變?yōu)?/p>

(39)

這兩個方程與量子力學中的諧振子滿足的方程一樣，其有界解為厄米高斯函數(shù)。

方程(39)有有界解的條件為λ1-1與λ2-1取非負偶數(shù)，這樣才能使方程(39)的級數(shù)解中斷成多項式解(這樣才有界)[8]。

λ1-1=2l,l=0,1,2,…

λ2-1=2m,m=0,1,2,…

(40)

這樣

(41)

以及

Xl(x)=NlHl(κx)e

Ym(y)=NmHm(κy)e

(42)

其中Hl(κx)和Hm(κy)是厄米多項式，歸一化系數(shù)

(43)

這樣自聚焦光纖中的光場模式可以表示成

(44)

其中βl m是對應于(l,m)模式的傳播常數(shù)。由式(36)和式(41)有

(45)

2.2 自聚焦的光纖模式色散

對應于(l,m)模式的群速度為

(46)

相應模式的群延時(傳播單位長度所需時間)為

(47)

當考慮近軸近似時

(48)

顯然(l,m)模式的群時延

(49)

是常數(shù)，與基模群時延τ0一樣，沒有模式色散。

考慮遠軸近似

(52)

其中用到了式(45)，可得最大模式號

(53)

將式(53)代入式(51)，得到最大群延時

(54)

則自聚焦光纖的模式色散為

(55)

這結果與幾何光學的處理結果式(21)完全一樣，這里仍然用到了遠軸近似。

3 結論

本文給出了一個完整的解決自聚焦光纖模式色散的幾何光學與波動光學的教學處理方案。在幾何光學的處理方法中，以費馬原理和光程的概念為基礎，利用光線方程求解自聚焦光纖的模式色散。其中引入了近軸光線與遠軸光線的概念，并討論了它們的聚焦特性，給出了與WKBJ法一致的自聚焦光纖的模式色散公式。將聚焦特性、近軸光線、遠軸光線等物理概念與數(shù)學上的不同近似程度緊密結合，幾乎在教學的每一步驟，都給出了清晰的物理圖像，明確了下列物理問題：(1)近軸光線近似與費馬原理預計的等光程性等價，不能用近軸近似結果來計算自聚焦光纖的模式色散； (2)自聚焦光纖的模式色散來源于遠軸光線； (3)由于聚焦特性的不同，自聚焦光纖的模式色散比階躍折射率多模光纖小兩個數(shù)量級。在波動光學的處理方法中，基于波動方程，將自聚焦光纖模式近似處理成為厄米高斯模式，在此基礎上，結合近軸近似和遠軸近似的處理，也給出了自聚焦光纖模式色散的正確結果。

[1]MarcuseD,Theimpulseresponseofanopticalfiberwith

parabolic index profile[J]. Bell Syst. Tech. J., 1973, 52: 1169-1174.

[2] 顧婉儀，李國瑞.光纖通信系統(tǒng)[M]. (修訂版). 北京：北京郵電大學出版社，2006：45-50.

[3] Ghatak A, Lokanathan S. Quantum Mechanics: Theory and Applications[M]. Dordrecht: Kluwer Academic Publisher, 2004：423-459.

[4] 伽塔克A K, 謝伽拉揚K.近代光學[M].袁一方,譯.北京：高等教育出版社，1987:308-335.

[5] Ghatak A.光學[M].4版.張曉光，席麗霞，余和軍,譯.北京：清華大學出版社，2013：47.

[6] 原榮.光纖通信[M]. 北京：電子工業(yè)出版社，2002：52-78.

[7] 張曉光，黃瓊，萬仁浚.光程、費馬原理和光纖的模式色散[J].大學物理，1996, 15(4):4-6. Zhang X G, Huang Q, Wan R J. Optical path, Fermats’s priciple and intermodel dispersion[J]. College Physics, 1996, 15(4): 4-6. (in Chinese)

[8] 曾謹言.量子力學卷 Ⅰ[M].北京：科學出版社，2007：513-516.

■

THE TREATMENTS IN THE VIEWS OF GEOMETRICAL OPTICS AND WAVE OPTICS FOR TEACHING INTERMODAL DISPERSION OF A PARABOLIC INDEX PROFILE FIBER

Zhang Xiaoguang

(State Key Laboratory of Information Photonics and Optical Communications,Beijing University of Posts and Telecommunications, Beijing 100876)

The graded-index fiber with the parabolic profile has an intermodal dispersion which is about two order smaller than that of the step-index fiber, and is widely used in local area network. The intermodal dispersion of this kind of graded-index fiber was conventionally treated with WKBJ approximate method which is originally used in quantum mechanics. This paper proposes the teaching treatments from the views of geometrical optics and wave optics, and obtains the same result as that calculated with WKBJ method. In the geometrical optics treatment, based on Fermat Principle, we introduce the new concepts of near-axial and far-axial rays, and give the solution to the ray equation using a perturbation approximation. In the wave optics treatment, based on the wave equation, the fiber modes in the graded-index fiber are treated as the Hermite-Gaussian modes. The intermodal dispersion is obtained by analyzing the propagation constant with the aid of near-axial and far-axial approximations. With these treatments we obtain the same result, with clearer physical significance. We believe that the teaching treatments proposed in this paper are more appropriate in the undergraduate courses of optical fiber communications.

waveguide optics; modes in optical waveguide; graded-index fiber with parabolic profile; intermodal dispersion

2016-05-31

張曉光，男，教授，主要從事物理教學與光通信研究，研究方向為光纖中的偏振現(xiàn)象與控制，光纖通信信號處理等，xgzhang@bupt.edu.cn。

張曉光. 自聚焦多模光纖模式色散的幾何光學與波動光學處理方法[J]. 物理與工程，2017，27(1)：23-29，43.