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

我是学习FDTD方法的新手，我在用matlab编程“Drude模型下负折射率材料中的差分方程”时候，没有图像，不知道是公式推导不正确，还是程序编程出现了问题，请高手帮忙给指点下。（附件中为差分公式及对应的程序）  ojyIQk+  
 谢谢！ .Asv%p[W  
%LVm3e9  
.. u/cg|]x&T  
MjlP+; !  

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

没有图像是什么意思呢？ <sb~ ^B  
发散了？

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

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

“没有图像是什么意思呢？ 发散了？” LENq_@$  
[>I<#_
^~  
是发散了，我运行普通介质时候，可以得出TE、TM情况下的场分布图； (XTG8W sN  
k=$TGqQY?  
可是当我加入Drude模型后，考虑负介电常数的时候，根据推导出的差分公式，编程的时候就出现了发散。 q>_.[+6  
LVM%"sd?  
我感觉自己很有可能犯错的地方在于：在处理PML边界的时候，没有迭代好，所以出现很大的发散。 dl
h)gp;  
6Gl
J>r+n  
我想的是在边界迭代时候，要考虑主区域的有耗介质中的迭代，又考虑pml层的迭代，但是我不知道怎样在边界上同时考虑这两种迭代（单独在主区域中或者在pml层中的迭代要好处理些）。 ^CYl\.Y@  
v
4TQX<0s  
我附件中的程序没有同时考虑主区域和pml层中的迭代，所以可能是这个原因使得发散很大。 -m zIT4  
+HpA:]#Y  
由于自己是刚学习这方面的新手，我想尝试自己动手编写自己设计的结构，所以又很多还不知道，以上只是我自己的想法，我其实也不是很清楚，所以想请熟悉这方面的高手给指点下。  tU5zF.%  
#lo6c;*m5  
谢谢！

你把pml层和主区域的迭代分开处理不行么？ 5lmHotj#  
TER=*"!  
这样至少保证pml层的是对的，如果再有发散，那应该就是主区域的问题了。

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

引用第1楼vincentwl于2009-09-26 22:22发表的  : jKAEm  
没有图像是什么意思呢？ X"|['t  
发散了？ '6iEMg&3  

y*jp79G  
*SbMqASv4G  
你帮他调试一下

刚刚搜到了 xiaowei580  发的帖子： 左手（负折射）材料资料目录。 9IdA%RM~mH  
非常感谢xiaowei580 ！我要消化下。

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