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

我是学习FDTD方法的新手，我在用matlab编程“Drude模型下负折射率材料中的差分方程”时候，没有图像，不知道是公式推导不正确，还是程序编程出现了问题，请高手帮忙给指点下。（附件中为差分公式及对应的程序）  {2F@OfuCF  
 谢谢！ Yhte&,D"  
5XoM)  
.. 5y8VA4L/o  
]p!)8[<  

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

没有图像是什么意思呢？ G %6P`:  
发散了？

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

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

“没有图像是什么意思呢？ 发散了？” L$ZjMJ  
CWj_K2=d  
是发散了，我运行普通介质时候，可以得出TE、TM情况下的场分布图； D tsZP (  
N'Gq9A  
可是当我加入Drude模型后，考虑负介电常数的时候，根据推导出的差分公式，编程的时候就出现了发散。 zTz}H*U  
hnnB4]c  
我感觉自己很有可能犯错的地方在于：在处理PML边界的时候，没有迭代好，所以出现很大的发散。 D!<F^mtl  
gD,&TW  
我想的是在边界迭代时候，要考虑主区域的有耗介质中的迭代，又考虑pml层的迭代，但是我不知道怎样在边界上同时考虑这两种迭代（单独在主区域中或者在pml层中的迭代要好处理些）。 ]-=L7a
 
54 lD+%E  
我附件中的程序没有同时考虑主区域和pml层中的迭代，所以可能是这个原因使得发散很大。 )`u17 {  
`+7F H  
由于自己是刚学习这方面的新手，我想尝试自己动手编写自己设计的结构，所以又很多还不知道，以上只是我自己的想法，我其实也不是很清楚，所以想请熟悉这方面的高手给指点下。 [);oj<  
#lYyL`B+~  
谢谢！

你把pml层和主区域的迭代分开处理不行么？ A*JOp8\)  
a:FU- ^B4~  
这样至少保证pml层的是对的，如果再有发散，那应该就是主区域的问题了。

我在推到公式时候分开处理，编程时主区域和边界也是分开处理的（在远离边界的地方），但就是在边界处，不知道怎样处理，看了葛德彪的书上说，可以采用平均值的方法，可是还是没看懂，具体怎样编程实现不知道处理，我准备先请教指点一下，再下去看。有可能是我的主区域的公式推导思路不正确，没有考虑到一些需要注意的因素。 B'Yx/c&n  
&RR;'wLoQT  
我刚刚在论坛上搜了下上面2楼gwzhao介绍的色散介质的FDTD推导过程和色散介质源代码，发现我推导公式的方法思路存在问题；下面引用的是有关的思路： /s?%ft#-9o  
(,At5T  
“色散材料一般会给出D和E的频域函数关系D=eps.E（本构方程），为了使用时域方法求解，必须将这一公式转换到时域，而转换的基本思路：分子中的iw（其中w表示入射波圆频率，i为虚数单位）用关于时间的一次偏导数替换（注意，w**2=-(iw)**2），从而得到关于D和E及其时间导数的方程。 cVi_#9u"  
采用FDTD求解时，磁场的求解仍按照普通FDTD公式处理，由Maxwell关系式中关于D的时间导数的方程由n时刻递推得到n+1时刻的D，然后由时域本构方程求出n+1时刻的E（若时域本构方程含二次时间导数，可能要用到n-1时刻的D、E及其时间导数，注意在编程的时候特殊处理）。”——hawkwhite l@`k:?  
Agcss20
.  
请高手给我推荐一下有关色散介质FDTD这方面的公式推导和编程实现的文献和书。 "4b{YWv  
35h|?eN_m!  
谢谢！

引用第1楼vincentwl于2009-09-26 22:22发表的  : }:?*n:g5  
没有图像是什么意思呢？ r@JMf)a]  
发散了？ Zzlt^#KLx  

5H'Iul<Os  
ZQ"dAR/y  
你帮他调试一下

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

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