利用变换域最小二乘法来改善FDTD算法中的数值色散特性研究.pdfVIP

2005’全国微波毫米波会议沦文集 利用变换域最牟帚委洼誊器蔫崔DTD算法中的数值色散特性 肖飞，唐小宏，张显静 (电子科技人学，电子T程系，电子系统和工程研究所，成都，610054) 摘要：本文讨论了一种对偏微分算r最优逼近的变换域最小二乘法。在保持原来时域有限差 分法(FDTD)的计算格式的前提下，仅仅选择不同的频率优化范围，就能利用这种方法确定 相应的优化系数，从而使FDTD算法的数值色散特性得到极大的改善。从而使这种方法其 有很大的灵活性，能满足不同问题的需要。 关键词：时域有限差分法，变换域最小二乘法，数值色散，稳定性条件。 时域有限差分法(FDTD)射t‘算电磁学中一种非常重要的方法。自从1966年出现以来， 它在电磁学中的很多方丽都得到了成功的应用。几十年过去了，仍然不断地有新的算法出现。 到目前为1I：，它已经发展成为一个庞大的算法类。具有代表性的FDTD算法有：基于Taylor finite 级数展开定理的T-FDTD算法f1．21j基于Nonstandard Fourier transform的PSTD算法【5】；基于 基于Wavelet的MRTD算法【4】；基于Fast method的US。FDTD算法[6]；等等。这些算法具有各自的优点 Alternating—directionimplicit 和缺点，它们的出现极大地丰富了FDTD算法的应用。然而，FDTD内在的数值色散是一个 很关键的误差来源。其结果导致数值相速成为频率和传播方向的函数。因此，改善FDTD 算法中的数值色散是非常重要的任务。一些研究者提出了一些方法【7．10]，但是仅仅针对某 一频率或者方向进行了改进。在本文中讨论了一种新的方法，它不象传统的方法在本域中来 考虑对偏微分算子的逼近，而是力图从变化域中来实现对偏微分算子的逼近。数值试验表明， 由新方法得到的FDTD并法的数值色散能够得到极大地改善。 1．变换域最小二乘法 为简单起见，本文忽略了边界的影响。以线性、各向均匀、无耗、无源媒质中的Il维 TM极化电磁波为倒。我们采用下面的中心差分格式来实现对时问偏微分算子的近似。 堂掣f，。“岳防；五叫，“]。一北训，=(。{)。] ㈣ 其中，d‘为待定系数。为了实现该中心差分格式在变换域中对时间偏微分算子进行最优逼 近，可以利用Fourier变换将(1)变换到频率域中。然后利用最小j乘法来确定系数∥。具体 做法如下： 先对(1)进行Fourier变换，得到： cz， 2 ， ∞z12dsin(10AAt r] 再利用最小二乘法来确定系数d 412 2005’全国微波毫米波会议论文集 q=卧划，sin(耕枥 ④ 墨：0 ’’(4) c?d‘ 其L卜l，石=coat为归一化频率。[-m．。，西。。]为感兴趣的归一化频率范围。 下面讨论对空间偏微分算子的近似处理。仪以x方向偏微分箕子为例，Y方向偏微分算 子可以做类似处理。仍然采用传统的四点中心筹分格式，如下： 望掣L。a石1蚤2靠弘札(1十阳嘲。川^批一z呦)小s， 对f5)进行Fourier变换，得到 颤z五2