[求助]3d-UPML程序,大家都来看看顺便给我解答一下疑问(有些疑问解决的地方我给出了注释)

我把时间步扩大了，但是没有出现数据上发散的现象，不知道你是什么情况，你可以再检查一下给我留贴咱们再讨论

这个你要看书了，书上有公式推导，论坛上有herosword师兄进行了翻译的中文版，你可以去下来看看

看得我投都晕了 ;aH3{TS  
怎么会有那么多的系数呢？ \b#`Ahf`  
我现在也在烦恼UPML的问题，但是我用的是葛得彪、闫玉波的书，发现没那么多的系数啊。 ?y,KN}s_  
不知道你研究过葛得彪的书没？ k<RaC=   
我在做这个程序，但出现了发散的现象 Ym2Ac>I4  
能不能多和你交流呢？
#;h> x  
我的QQ是24765428

哦，我现在再仿波在无限长的微带线上面传播，用的就是这个upml吸收边界，可是老有反射！有一部分介质层和upml重合了，想问一下在介质层里面电导率sigmam应该怎么取值呢？还是和上面的一样吗？

ki是怎么定义的啊？

我的qq243640120

不知道为啥加了个电场强度的激励元exp(-((n-ndelay)^2/tau^2))；就不对了呢? HE0m#  
也是遇到改变激励元就发散，谁能解释一下

你把参数tau调整一下，我以前也是改了激励后发现有反射，调整了tau参数就没事了，你不妨试试

基于楼主的FDTD程序，把激励源设置成ez(is,js,ks)=(n-ndelay)*sin(0.5*(n-ndelay))*exp(-((n-ndelay)^2/tau^2)); 5ka6=R(r  
然后画出了ez(is,js,ks)随时间变化的图，结果ez(is,js,ks)的最大值为0.3，但是sin(0.5*(n-ndelay))*exp(-((n-ndelay)^2/tau^2))的最大幅值为1； WT}xCni  
这是为什么呢？ V5gr-^E  
还有几个问题 ^-F#"i|Cn  
PML的吸收特性不是基本与频率无关么？为什么加上有的高斯脉冲激励后，在激励附近的场值会是负值， QY*F(S,\  
还有，仿真在真空中传播，为什么空间另外一点的场值衰减会比较大 ，加上软源时，即使在ez(is,js,ks)点，它的幅值衰减为加入的激励的1/3; 7R
n 4gT  
加上软源时，激励点的电场图形和原先加上的激励为啥差别这么大？  &w:0ad|  
我的微波基础较差，如果问题比较幼稚，请轻拍 bk]g}s  
9<o*aFgCa  

这个程序的PEC WALL为什么要这样设置，直接将最外面一层的电场设为0不就是PEC了么，不就完事了么？ dP"cm0  
这个程序在系数上用电导率设为无穷大代入系数计算得出-1 和0 ，本质上和直接将电场设0有区别么，两个PEC的方法哪个更好，编程时又有什么区别？