[求助]FDTD的问题求助

在进行Ez的迭代的时候需要用到周围的上一步长的磁场，那在迭代之前是不是要先将磁场备份以便在推得Ez的时候使用 [FC7+ Ey^  
HPCgv?E3  
但是这样是不是会额外的消耗内存？ i?'HVx  
x-U:T.+{  
大家是怎 .. brNe13d3~"  
@|%t<{y^I  

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

如果你的迭代公式是高阶的  那就要根据需要保存上一次的场值  或者还要保存上上次的。
p#qQGJe  
内存的消耗是必须的。

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

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

在TMz模式下的电流源叠加时推导得到如下公式： 8^IV`P~2M  
u<L<o2  
Ez=Ez-dt/(er0*er)*Jz Sg%h}]~   
1O#
]qZS}]  
之后的这个Jz怎样求？ l7h6R$7; 0  
:1*q}R   
我推导得到的是 bG0 |+k3O  
Jz=I/(dx*dy) 5D]3I=kj  
其中I是想要注入的电流源，dx,dy分别是x与y的网格的大小 Wd`*<+t]  
TxJoN]Z.  
但是结果好像不对

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

为什么我在两次设置不同的dx 和dy后迭代的结果不一样？ ?@z/#3b  
如下图所示？ 9[`c"Pd  

图片:1.jpg

Lu~E5 ,  
其中黑色的是疏网格的，绿色的是密网格的 d-
C%R9  
UlWm). b;v  
有经验的帮忙指点一下。

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

将代码删除了

想用电流源1mA注入算端口上的总电场，总电流等， 5-$D<}Z  
不知道怎样注入这个电流源，这个大小如何计算