编程“Drude模型下负折射率材料中的差分方程”时没有图像，请高手帮忙给指点一下.

我是学习FDTD方法的新手，我在用matlab编程“Drude模型下负折射率材料中的差分方程”时候，没有图像，不知道是公式推导不正确，还是程序编程出现了问题，请高手帮忙给指点下。（附件中为差分公式及对应的程序）  ZHRMW'Ne  
 谢谢！ Bz{"K  
kJ* N`=  
.. mb'{@  
2x gk$E$7  

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

没有图像是什么意思呢？ H96|{q=  
发散了？

你可以参考一下FDTD板块里面色散介质的源代码。

“没有图像是什么意思呢？ 发散了？” mB"1QtD  
t$rla_rbY  
是发散了，我运行普通介质时候，可以得出TE、TM情况下的场分布图，可是当我加入Drude模型后，考虑负电解质常数的时候，推导出有耗迭代差分公式，再编程的时候，就出现了发散。我感到自己有可能犯错的地方在于，编程时候，场分量的迭代公式，在处理PML边界的时候，没有迭代号，所以出现很大的发散，可是我考虑边界时，出现了难题，就是我用的矩阵法编的程序，在边界迭代时候，既要考虑主区域的有耗介质，又要考虑pml层的迭代，所以在边界上处理的不好，肯能就是这样使得发散很大。这只是我自己的想法，我也不是很清楚，所以想请熟悉这方面的高手给指点下。谢谢！

“没有图像是什么意思呢？ 发散了？” l>KkAA
 
$lq.*UQ;0  
是发散了，我运行普通介质时候，可以得出TE、TM情况下的场分布图； $m A2AI  
m=SI *V  
可是当我加入Drude模型后，考虑负介电常数的时候，根据推导出的差分公式，编程的时候就出现了发散。 r-.>3J  
ck\W'Y*Q7  
我感觉自己很有可能犯错的地方在于：在处理PML边界的时候，没有迭代好，所以出现很大的发散。 /aIGq/;Y+a  
@_1cY#!  
我想的是在边界迭代时候，要考虑主区域的有耗介质中的迭代，又考虑pml层的迭代，但是我不知道怎样在边界上同时考虑这两种迭代（单独在主区域中或者在pml层中的迭代要好处理些）。 dFyGI?  
Qqs1%u;e8  
我附件中的程序没有同时考虑主区域和pml层中的迭代，所以可能是这个原因使得发散很大。 'by+hXk  
Ov^##E  
由于自己是刚学习这方面的新手，我想尝试自己动手编写自己设计的结构，所以又很多还不知道，以上只是我自己的想法，我其实也不是很清楚，所以想请熟悉这方面的高手给指点下。 jcb&h@T8kv  
u1=K#5^  
谢谢！

你把pml层和主区域的迭代分开处理不行么？ J=k=cFUX  
}ML2-k  
这样至少保证pml层的是对的，如果再有发散，那应该就是主区域的问题了。

我在推到公式时候分开处理，编程时主区域和边界也是分开处理的（在远离边界的地方），但就是在边界处，不知道怎样处理，看了葛德彪的书上说，可以采用平均值的方法，可是还是没看懂，具体怎样编程实现不知道处理，我准备先请教指点一下，再下去看。有可能是我的主区域的公式推导思路不正确，没有考虑到一些需要注意的因素。 YE~IO5   

_ERtL5^  
我刚刚在论坛上搜了下上面2楼gwzhao介绍的色散介质的FDTD推导过程和色散介质源代码，发现我推导公式的方法思路存在问题；下面引用的是有关的思路： G<n75!  
$3G^}A"  
“色散材料一般会给出D和E的频域函数关系D=eps.E（本构方程），为了使用时域方法求解，必须将这一公式转换到时域，而转换的基本思路：分子中的iw（其中w表示入射波圆频率，i为虚数单位）用关于时间的一次偏导数替换（注意，w**2=-(iw)**2），从而得到关于D和E及其时间导数的方程。 6]Hwr_/tk  
采用FDTD求解时，磁场的求解仍按照普通FDTD公式处理，由Maxwell关系式中关于D的时间导数的方程由n时刻递推得到n+1时刻的D，然后由时域本构方程求出n+1时刻的E（若时域本构方程含二次时间导数，可能要用到n-1时刻的D、E及其时间导数，注意在编程的时候特殊处理）。”——hawkwhite '<BLkr# @  
~TEn +  
请高手给我推荐一下有关色散介质FDTD这方面的公式推导和编程实现的文献和书。
E|"SMA,  
ZBK0`7#&EH  
谢谢！

引用第1楼vincentwl于2009-09-26 22:22发表的  : _ox+5?>  
没有图像是什么意思呢？ j7~FR{:j  
发散了？ h:?^0b!@  

YMX9Z||  
4@PA+(kvS  
你帮他调试一下

刚刚搜到了 xiaowei580  发的帖子： 左手（负折射）材料资料目录。 R;< q<i_l  
非常感谢xiaowei580 ！我要消化下。

你下资料方便的话，可以找找毕业论文，我记得好像有好几个硕士毕业论文里面都有色散介质的FDTD公式推导什么的。