[求助]FDTD的问题求助

在进行Ez的迭代的时候需要用到周围的上一步长的磁场，那在迭代之前是不是要先将磁场备份以便在推得Ez的时候使用 L&gEQDPgq|  
<GFB'`L  
但是这样是不是会额外的消耗内存？ Bdk{.oh6  
KxTYc  
大家是怎 .. 5@&i:vs5y  
8Jy1=R*S  

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

如果你的迭代公式是高阶的  那就要根据需要保存上一次的场值  或者还要保存上上次的。 Q#bo!]H{t  
内存的消耗是必须的。

我的不是高阶的，是TMz假设下的二维问题，因该是最简单的，难道还需要保存吗？

这……你如果不保存能算出结果，那就不保存。如果必须保存了才能算出结果，那就保存。 Z'~FZRF  
而且二维TM模式的FDTD  除了Mur2阶边界，没有地方会用到上一次的值，你是不是对算法理解有出入？

在TMz模式下的电流源叠加时推导得到如下公式： }\F>z  
,$Mw/fA  
Ez=Ez-dt/(er0*er)*Jz +}!eAMQ  
$n) w4p_  
之后的这个Jz怎样求？ aap:~F{]X  
?bEYvHAzg  
我推导得到的是 <5 +?&i  
Jz=I/(dx*dy) &nss[w$%C  
其中I是想要注入的电流源，dx,dy分别是x与y的网格的大小 aNz%vbh\  
pJIH_H  
但是结果好像不对

回复三楼，是的，我的时间推进的时候错了，我的是将所有的HX，HY，Ez都与前一步的HX，HY，Ez有关，看了看书更正过来了。

为什么我在两次设置不同的dx 和dy后迭代的结果不一样？ CV`  I.  
如下图所示？ XW1
9hG  

图片:1.jpg

<%!@cE+y  
其中黑色的是疏网格的，绿色的是密网格的 F4k`x/ak  
{u:DC4eut  
有经验的帮忙指点一下。

如果是1D的，就把code贴上来吧。

将代码删除了

想用电流源1mA注入算端口上的总电场，总电流等， fsb_*sh&  
不知道怎样注入这个电流源，这个大小如何计算