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

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

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

这是我修改后的图

    for(i=0;i<PML_lx;i++) n 1MZHa,  
    { &C<K|F!j!  
        xnum=PML_lx-i; cHOtMPyQ  
        xd=PML_lx; XlUM~(7+v  
        xxn=xnum/xd; B*7Y5_N  
        PML_factor=(dt/2./Eps0)*((PML_m+1.)/(150.*Pi*ddx*sqrt(PML_epsr))); Z"PPXv-<jY  
        xnpml=PML_factor*pow(xxn,PML_m);         c7Sa|9*dR  
        gi1=xnpml; j78WPG  
        gi1[IE-2-i]=xnpml; ~$n4Yuu2[  
        gi2=1.0/(1.0+xnpml); L;zwqdI  
        gi2[IE-1-i]=1.0/(1.0+xnpml); k8
H@0p  
        gi3=(1.0-xnpml)/(1.0+xnpml); y\Kr@;q0w  
        gi3[IE-1-i]=(1.0-xnpml)/(1.0+xnpml); d4KTwn5g  
        xxn=(xnum-0.5)/xd; iSu7K&X9q  
        xnpml=PML_factor*pow(xxn,PML_m); lxb+
0fiN  
        fi1=xnpml; RUqN,C,m5I  
        fi1[IE-1-i]=xnpml; 's>  
        fi2=1.0/(1.0+xnpml); HE58A.Q&  
        fi2[IE-2-i]=1.0/(1.0+xnpml); 7S$Am84%  
        fi3=(1.0-xnpml)/(1.0+xnpml); 6yk=4l\  
        fi3[IE-2-i]=(1.0-xnpml)/(1.0+xnpml); VZo[\sWf  
    } 1O@ qpNm  
PML扎按照红色部分修改吸收效果就好多了！ 1=(jpy  
这只是x方向的其他方向类似。

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

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

Eps=08.85e-12;PML_factor=0.333

引用第14楼lx84于2010-12-07 16:52发表的  : )%6h9xyXt  
    for(i=0;i<PML_lx;i++) i .GJO +K  
    { v{44`tR  
        xnum=PML_lx-i; x,rK4L7U  
        xd=PML_lx; |H%[tkW6c  
        xxn=xnum/xd; mFa%d8Y  
....... 7J6D wh{  

cmu|d  
D_D76  
sPxDo?1x-  
楼主给出的程序是用于计算介质球双站RCS的，而lx84兄弟在14楼贴给出的图是计算双层球（外部为介质，内部为金属球）的结果。 +WTO_J7  
我按lx84兄弟的方法修改了PML运行了一下后，感觉PML这样改过后跟原来的一样，并没起作用。当然也就得不到14楼那样的图。 u6lcl}'  
所以我觉得这个程序的问题应该不是出在吸收边界上面。 QrSF1y'd  

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

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