fdtd数值计算在微带天线分析中的应用

《fdtd数值计算在微带天线分析中的应用》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、电磁场报告电磁场期末课题报告FDTD数值计算在微带天线分析中的应用合作成员：竺可桢学院：混合005班张丹（3003001102）混合007班邓坚（3003001152）分工情况：张丹（3003001102）：负责理论阐述，公式推导，理论计算邓坚（3003001152）：负责程序实现，matlab仿真日期：2003年1月8日86电磁场报告目录时域有限差分基本知识概述3标量麦克斯韦方程5FDTD方程6解的稳定条件9边界条件12激励源的类型14时域有限差分在微带天线中的应用16微带贴片天线数值模拟的若干问题

2、16微带天线FDTD数值计算模拟实例21微带天线数值计算的基本理论介绍21激励源的设置22边界的处理23散射效应25有效长度，有效宽度以及谐振频率27程序原代码29Main函数部分29Coefficient函数部分40FT函数部分51Impendence函数部分55output函数部分61parameter函数部分62pattern函数65图形显示程序71程序分析74信号源74参数计算74场量的计算75等效电阻的计算77傅立叶变换77空间的方向图78部分语句的诠释80程序结果81方向图824个面上的辐

3、射方向的截面图83心得体会84参考文献…………………………………………………………………………………………8586电磁场报告FDTD数值计算在微带天线分析中的应用时域有限差分基本知识概述时域有限差分（FDTD）方法是60年代由K.S.Yee首先提出来，不仅在三维空间对场量离散，而且在时间轴上也对场量离散，并且用差分代替微分，得到麦克斯韦方程在一定边值与初值条件下的近似解。作为麦克斯韦方程的基本方程取为：(1.1)(1.2)式中：了解E和D及H和B的关系很重要。在真空和自由空间中，它们存在着简单的比例关

4、系，即其中，称为真空中的电容率。在一般各向同性媒质中具有简单关系式：和分别是媒质的磁导率和电容率对于各向异性媒质，上述关系变为下式：式中和分别称为张量磁导率和张量电容率。此处B的每一个分量和H的三个分量都有关（D和E关系相似）,例如等离子体和铁氧体的张量磁导率和张量电容率在忽略损耗时可分别为：和一般的讲，均匀各项同性媒质中和是标量常数。对于非均匀各项同性媒质中，和是随空间位置而变化的标量。对于各向异性媒质，和86电磁场报告是张量，均匀还是不均匀取决于和是否随位置的变化而变化。麦克斯韦方程的另外一种解的

5、形式是它的时谐场形式。假定电磁场随时间变化具有以下形式：对于无源区麦克斯韦第一、第二方程式化为复数形式(1.3)引入复数介电常数形式同理，在有磁损耗时亦可以表示为复数形自然界中很多媒质如水、铁磁材料、等离子体、生物体、各种吸波材料，通常其复介电常数、复磁导串是频率的函数。对于这种媒质，存在下述关系式：上述方程都是麦克斯韦方程在频域内的形式。而FDTD算法的主旨在于把场的问题在时域内进行计算，其基本方程的得到是先将麦克斯韦的旋度方程在直角坐标系中展成六个标量场分量的方程，再将问题空间沿三个轴向分成很多网

6、格单位，每个单位长度作为空间变元，响应得出时间变元。用有限差分式表示关于场分量对空间和时间分量的微分。标量麦克斯韦方程假定研究的空间是无源的，并且媒质参数不随时间而变化，在直角坐标系中，麦克斯韦方程（1.1）和（1.2）将化为下列六个标量方程：86电磁场报告(1.4)当不考虑磁损耗时，令，则以上方程与麦克斯韦方程等效。上述六个方程构成了完整三维问题的情形。在实际中，往往会遇到被研究模型沿一个轴向或二个轴向不变化的特殊情形。对于前者，设沿z向不变化(二维情形)，亦即是全部场分量对的偏导数等于0，则上述方

7、程组简化为二组独立的标量方程组，一组称为横磁模(TM)，另一组为横电模(TE)。对于后者，设被研究模型沿z和x二个轴向(一维情形)上无变化，即对z和x的偏导数为0，此时上述方程级简化为仅有二个场分量的标量方程组，称为横电磁横(TEM)。各简化情形的方程组如下1．二维情形对于TM模对于TE模2．一维情形86电磁场报告此时为TEM模，分二种情形仅有，仅有，FDTD方程将麦克斯韦六个标量方程中场分量对坐标和时间的偏导数用有限差分式来表示。其基本做法是，将问题空间沿三个坐标轴向上分成很多网格单元，用、和分别表

8、示单元沿三个铀向的长度，用表示时间增量。网格元顶点的坐标(x，y，z)可记为任意—个空间和时间的函数可表示为这里和为整数。其次，用中心有限差分式来表示函数对空间和时间的偏导数，这种差分式具有二阶精度，其表示式为：(1.5)为了实现空间坐标的差分计算，并考虑到电磁场在空间相互正交和铰链的关系，在基本网格单元上六个场分量的怔旨如图1．1所示86电磁场报告考虑到在时间上E和H有半个时间步的变化，由(1.4)方程组的第一式得到：取磁场分量的规一化值，即，则由上式