用MATLAB开发的FDTD系统研究电磁散射远场特性.docVIP

用MATLAB开发的FDTD系统研究电磁散射远场特性 ? ?中国电磁兼容网 ?? 摘要： 利用LOVE场等效原理，实现了从近场到远场的变换，有助于时域有限差分法对电磁散射远场问题的解决，并在此基础上开发了建立在MATLAB平台上的时域有限差分法软件系统，以方便对电磁散射远场问题进行研究。利用该系统解决了三个电磁散射实例中的远场问题，并将其计算结果与理论或其它计算方法所得的结果进行了比较，结果吻合较好。　　关键词： 远场特性； 时域有限差分法； MATLAB Study on Far-field Characteristics of Electromagnetic ScatteringProblems with FDTD System Based on MATLAB Plane DOU Wei-ping， ZHANG Lin-chang(EMC Research Section, Northern Jiaotong University, Beijing 100044, China) 　　Abstract: LOVEs Equivalence Theorem is presented to solve far-field characteristics of electromagnetic scattering problems in finite-difference time-domain method. Correspondingly, a finite-difference time-domain method system based on MATLAB plane is developed by the authors. Some examples are studied by the system, and results are compared with ones calculated with other methods for validation.　　Key words: Far-field characteristics； Finite-difference time-domain method； MATLAB 1　前言 　　由于所需计算内存少、能模拟复杂的介质结构等诸多优点，时域有限差分法(FDTD)已经成为研究电磁散射的主要数值方法，但由于计算机存储量的局限性，很难实现将计算的空间扩大到令人满意的远场范围，故如何实现从近场到远场的变换成为求解电磁散射远场的一个焦点。实现从近场到远场的变换一般有两种方案：一是在频域范围内进行变换［1］，二是在时域范围内进行变换［2］。分别利用相应的等效原理，求解任意空间点的远场。鉴于已知被研究电磁散射问题的源的频率，本文采用的是第一种方案，从而可以避免考虑第二种方案中存在的较为复杂的场值收敛问题。近远场变换采用LOVE场等效原理，通过FDTD法计算所得的近区内的电磁场分量，求解远场的分布。　　为实现上述思想，我们专门开发了FDTD应用软件系统，并以MATLAB5.0作为开发工具。MATLAB是一种功能强、效率高的便于进行科学和工程计算的交互式软件包。它集丰富的数学函数库和强大的绘图功能于一体，并可在WINDOWS95下直接运行，具有漂亮方便的图形操作界面和多任务运行等优点，因此它在数值计算等领域体现的优势是诸如Fort-ran 5.0、NDP Fortran 、Borland C++等一般高级语言无法比拟的。 2　算法介绍 　　到目前为止，FDTD法的基本理论已经相当成熟，由于篇幅限制，这里不对FDTD法本身做具体介绍，仅给出作者用计算机编程语言给出的FDTD法中电磁场六个分量的基本迭代公式：　　Ex(i,j,k)=Cx(i,j,k){-Dx(i,j,k)Ex(i,j,k)+［Hz(i,j+1,k)-Hz(i,j,k)］／dy　　　　　　　-［Hy(i,j,k+1)-Hy(i,j,k)］／dz}　　Ey(i,j,k)=Cy(i,j,k){-Dy(i,j,k)Ey(i,j,k)+［Hx(i,j,k+1)-Hx(i,j,k)］／dz　　　　　　　-［Hz(i+1,j,k)-Hz(i,j,k)］／dx}　　Ez(i,j,k)=Cz(i,j,k){-Dz(i,j,k)Ez(i,j,k)+［Hy(i+1,j,k)-Hy(i,j,k)］／dx　　　　　　　-［Hx(i,j+1,k)-Hx(i,j,k)］／dy}　　Hx(i,j,k)=C′x(i,j,k){-D′x(i,j,k)Hx(i,j,k)-［Ez(i,j,k)-Ez(i,j-1,k)］／dy　　　　　　　+［Ey(i,j,k)-Ey(i,j,k-1)］／dz}　　Hy(i,j,k)=C′y