兩步分裂步長時域有限差分法

林智參

摘要：提出一種新的兩步分裂步長時域有限差分（TS-FDTD）法，該方法基于Split-Step方案和Crank-Nicolson方案，采用新的矩陣分解形式，與傳統的FDTD算法、傳統分裂步長時域有限差分法相比，減少計算復雜度，算法的推導簡單，提高了計算精度。本文還加入一階Mur吸收邊界條件，給出一階Mur吸收邊界差分方程。最后，通過實例仿真，比較TS-FDTD、傳統FDTD方法兩種算法的仿真結果，驗證了TS-FDTD算法的可行性及其高精度性。

關鍵詞：時域有限差分法；分裂步長；邊界條件；精度

中圖分類號：TM15 文獻標識碼：A 文章編號：1007-9416（2018）05-0144-02

時域有限差分法[1]（Finite-Difference Time-Domain-FDTD）是一種簡便的電磁波時域分析方法，此方法用Yee氏網格為基礎，把電磁場離散化，將麥克斯韋旋度方程差分化，建立差分方程，從而簡便有效的處理各種電磁場中復雜的問題，目前已經廣泛的應用于電磁場的各個方面。但是，傳統的FDTD算法也有不足之處，其推導公式較為復雜，運算過程負擔頗大大，因而，人們也從多個方向對FDTD進行改進[2]。本文提出了一種基于Split-Step[3]方案和Crank-Nicolson[4]方案新型FDTD算法，以TM波為例子，采用一種新的矩陣分解方法，來簡化運算公式，減輕計算負擔。

1 TS-FDTD算法理論推導

選擇一無源區域作為研究空間，其中介質均勻無耗并且各向同性，介電常數為ε，磁導率為μ，可將二維TM波麥克斯韋方程組以微分形式表示如下：

在分步2中，電場分量Ez在二維空間四個邊界上的一階Mur吸收邊界差分方程式可參考分步1，其形式近似，此處不再展開贅述。

3 實例仿真

本例將TS-FDTD算法用于運算一個二維自由空間TM波傳播及電場分布情況，空間的尺寸大小為100cm*100cm，并且采用一階Mur邊界條件，激勵源為sin（2*pi*f*t），放置于二維空間的中間位置，激勵源的頻率為f=1.5GHz。……