[求助]FDTD谐振点位置计算误差较大

还有想问一下，有没有做频谱分析的好的代码，我不想自己编这个程序，想直接用别人的函数了 gp'n'K]  
uuhvd h=  
还有现在的很多FFT还要自己去编程序得到其频率分量，有没有和CST中的一样的直接对数据（时间，幅度）进行分析，得到（频率，幅度，相位）的函数？

贴上我用debye计算的结果 6;pREM+  

图片:未命名.jpg

\-0`%k"&  
都很多步长了，都感觉没有收敛

什么吸收边界？ zB6u-4^wT  
FFT自己写，很简单的，你设置频率点时设置大一些，少一些频率点计算，这样速度会快很多。

不是的，以前我做等离子体色散介质的时候，和其它介质一样，没什么特别要求步长的。 Kajkw>z  
你肯定程序中哪里有问题吧，就算通常的介质要几千步收敛也很夸张啊，这样万一哪里有错误，你怎么调试程序？

v1 ,v2 是什么？ 这个图是什么意思 ？ #G[t X6gU  
感觉你程序有问题，一般正常的图都很平滑的。 / $_M@>  
!~)90Z!  
matlab工具箱里不是有很多函数么，自己调用好了。 7{4w
2)  
还有个时频分析工具箱，里面就更加全面了。

图中的V1是时间（单位秒)，V2是电流（单位安） R&*@@F-dx  
bY)#v?  
X
h}&uZ`A  
MATLAB中只能得到FFT变换后的值，变换后还需要自己得到需要的频率以及幅度，而不是直接拿过来用。 oQ\&}@(V  
k"0;D-lTZ>  
这个程序我调了好几天了，都不知道哪错了

刚刚看了看文献，许多文献计算我说的模型复介电常数问题er=er1-jer2 nB] Ia?  
用的不是直接的debye模型而是用(FD)2TD模型 :kG)sw7  
不知道是什么原因？ D(\$i.,b2  
!%[S49s  
上传文献如下： gKi{Y1
 
i=rH7k  
时域有限差分法中有耗媒质的介电常数.pdf (271 K) 下载次数:3 /{*$JF  
zK;XFN#U^  
请知道这方面区别的给一些建议，谢谢！

到今天早上为止，这个问题我已经完全通过FDTD实现了，感谢大家的帮忙，帖上我的结果对比图 Me3dpF  
\?fIt?  

图片:未命名.jpg

gE7L L=x  
总结： (P|pRVO  
1、书上的所有FDTD公式必须自己从头到尾推导一遍，不论其多么简单和复杂 nHXX\i  
2、特别注意公式当中的系数，是相对介电常数还是绝对介电常数，是总时间还是步进时间，是步进次数还是时间等都要仔细 +0$/y]k  
3、注意FDTD中的物理量的求解速度 FY3IUG  
以我的二维为例： chI.{R
j  
Hx->Hy->Ez->source->Boundary :l u5Uu~  
TLa]O1=Bf.  
evuZY X@  
另外，通过结果比对，发现谐振点处的频率有点偏移，我想是不是因为网格的原因，如果采用自适应技术和subgridding技术是否可以让结果很接近CST的结果？

计算结果很好了。 nv0@xnbz  
从曲线看，高频的误差大，可能是网格步长取值大了造成了 e9pOisZ;8  
“采用自适应技术和subgridding技术”自己实现这些技术，那后面的问题就更多了，呵呵

感谢老大的不断关注。