[求助]请教一个FDTD的EH计算数量差别的问题？

在FDTD计算中中， <?!#QA  
Gh.02  
一维情况下，电场比磁场多一个，为什么要多一个呢？ JyV"jL   
>:.Bn8-  
|1EM )zh6  
二维情况下，电场Exy还是比磁场Hz多一个，或者Ez比Hxy多一个吗？ Z+4J4Ka^!(  
i4l?q#X  
三维时，Exyz也是 .. 6w'^,V  
BX?Si1c  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册

自己动手编一下就知道为什么了。 Y'wQ(6ok  
电场（磁场）计算都要涉及到前后磁场（电场）的值。 G7b>r  
比如1D， re:=fC:t5A  
qh{hpX)\D  
Ex (1:n) = Ex(1:n) + factor * (Hy(2:n+1) - Hy(1:n)) EHmw(%a|+  
Hy(1:n) = Hy(1:n) + factor * (Ex(2:n+1) - Ex(1:n)) /g@!#Dt  
En
yx+]9  
所以Ex要计算n个点的话，Hy就需要n+1个值，同样如果Hy要计算n个点的话，Ex就需要n+1个值。 }9=VhC%J  
至于到底哪个多一些，就看你自己编程了，只不过大家习惯E多一个。

哦，感谢斑竹回复！ O8qA2@,  
eh`n?C  
那么一维时我让Ex多一个。 7Rwn{]r  
!/2uO5  
当转为二维TE，是不是让Hz，比Ex和Ey多一个？或者Ex和Ey比Hz多一个？

恩，2D的时候也会这么处理的。 ):<9j"Z;At  
#z$g1\v  
你可以参考别人的代码，但看完一定要自己写。 :
9?y-X  
5|:t$  
写之前把推导公式写好，想一下可能遇到什么问题，怎么解决。 <gfkbDP2  
Ga,+  
如果出来的图形不对，再想一下，哪里不对，为什么不对。 C[c^zn  
MbbKo-7F$  
这样才能说掌握了这段代码。

     呵呵，楼主的帖子把我差点弄迷糊了，看了gwzhao的回复才知道楼主原来是问这个…… a@?2T,$  
FDTD的独特之处就在于电、磁场在空间和时间上都是交错排列的——其实你要想让电、磁场场量放置点一样多也是可以的，比如一维平面电磁波传播问题（不失一般性，假设电磁波从左向右传），左边界上的电场为给定的随时间变化的激励，而右边边界设在半点上，然后给磁场的辐射条件(如果给PML条件也是可以的，即只将右边界的半点放到导体里边并给定为0），这样就可保证存储电、磁场的点数一样多——多一个少一个其实不重要，是看需要和方便来定的。 8n2MZ9p]  
    另外，时域伪谱法（PSTD）以及我正在研究的时域有限体积法（FVTD）都是将电磁场放在同一点上，电、磁场点的数量自然也是一样多的。

谢谢！版主和hawkwhite

关键你要弄清楚Yee元胞到底是怎么回事  这样一维TEM模式 二维TE模式 TM模式  才能根据需要来分配网格 进而知道怎么分配数组大小  硬记E或H多一个意义不大