大家帮忙看看这个3D_FDTD介质球的RCS有什么问题

吸收边界有问题，从新设置下。

11楼的兄台，如果吸收边界有问题，具体是哪里的问题，可否指出来？

这是我修改后的图

    for(i=0;i<PML_lx;i++) 2@pEuB3$?!  
    { OyTp^W`&  
        xnum=PML_lx-i; !*=+E%7  
        xd=PML_lx; s#V:! 7  
        xxn=xnum/xd; YnX6U1/^  
        PML_factor=(dt/2./Eps0)*((PML_m+1.)/(150.*Pi*ddx*sqrt(PML_epsr))); c+]5[6  
        xnpml=PML_factor*pow(xxn,PML_m);         YDiru  
        gi1=xnpml; hkR Jqta)  
        gi1[IE-2-i]=xnpml; t0o'_>*?A  
        gi2=1.0/(1.0+xnpml); RhumNP<M  
        gi2[IE-1-i]=1.0/(1.0+xnpml); 8Sg:HU\  
        gi3=(1.0-xnpml)/(1.0+xnpml); Kv6#WN~  
        gi3[IE-1-i]=(1.0-xnpml)/(1.0+xnpml); #W=H)6  
        xxn=(xnum-0.5)/xd; )R"deb=s  
        xnpml=PML_factor*pow(xxn,PML_m); Z$Mc{  
        fi1=xnpml; $4]"g}_  
        fi1[IE-1-i]=xnpml; 4(}V$#^+  
        fi2=1.0/(1.0+xnpml); ;NrU|g/ksX  
        fi2[IE-2-i]=1.0/(1.0+xnpml); n ,CMGe^:  
        fi3=(1.0-xnpml)/(1.0+xnpml); y ;T=u(}  
        fi3[IE-2-i]=(1.0-xnpml)/(1.0+xnpml); k[;(@e@c  
    } ~bSjZ1`  
PML扎按照红色部分修改吸收效果就好多了！ gX*i"Y#
 
这只是x方向的其他方向类似。

谢谢lx84学友的回答，我会按照你的方法进行调试！如有问题再请教，感谢！！！

有一个小问题：Eps0是真空中的介电常数吗？PML_epsr以及PML_m分别代表什么呢？上述PML公式来自什么资料呢？

Eps=08.85e-12;PML_factor=0.333

引用第14楼lx84于2010-12-07 16:52发表的  : /s
91[n(d  
    for(i=0;i<PML_lx;i++) A^/$|@  
    { Rj&V~or  
        xnum=PML_lx-i; JI{|8)S  
        xd=PML_lx; ?T+Uu  
        xxn=xnum/xd; RQzcsO  
....... n9.` 5BH7/  

y6.Q\=  
*82f{t]  
>Qk4AMIO  
楼主给出的程序是用于计算介质球双站RCS的，而lx84兄弟在14楼贴给出的图是计算双层球（外部为介质，内部为金属球）的结果。 9$-V/7@)  
我按lx84兄弟的方法修改了PML运行了一下后，感觉PML这样改过后跟原来的一样，并没起作用。当然也就得不到14楼那样的图。 oLJP@J  
所以我觉得这个程序的问题应该不是出在吸收边界上面。 J&%d(EJM  

pml影响RCS的幅值，但大体趋势是一样的。你要注意电场的极化方向，个人认为本程序的极化方向和葛德彪书上的介质球的散射极化方向不同。

楼住 ,,你有没有搞定这个问题?