[RFEDA原创]发一个二维TM FDTD 平面波入射到金属方柱上散射程序

谢谢分享，有fortran版的吗？ i!}6FBZ  

课设中要算RCS，学习学习

好东西，非常感谢

多谢楼主的分享，我要做一个三角形的，不过先看懂你这个，然后再根据你的来修改。谢谢哈！！！

谢谢分享劳动成果，雪中送炭啊，呵呵

谢谢你的程序，如果是matlab就更好了

//边界Mur ABC D 6
y,Q  
            //i=0; x的左边边界 0]
  
            for(j=1;j<NUM_of_YAXIS-1;j++) oS..y($TI  
            { PgK7CG7G  
             ez[0][j]=ez_x_low_m1[j]+ce*(ez[1][j]-ez[0][j])-cf*(hx[0][j]-hx[0][j-1]+hy[1][j]-hy[1][j-1]); 77)WNL/ x  
            //ez[0][j]=ez_x_low_m1[j]+ce*(ez[1][j]-ez[0][j]); 1,;X4/*  
            ez_x_low_m1[j]=ez[1][j]; A|+QUPD  
            } ; iia?f1  
            // x右边 的边界 /o m++DxV  
            for(j=1;j<NUM_of_YAXIS-1;j++) }}tbOD)t  
            { S$Zi{bU`G  
            ez[NUM_of_XAXIS-1][j]=ez_x_high_m1[j]+ce*(ez[NUM_of_XAXIS-2][j]-ez[NUM_of_XAXIS-1][j])-cf*(hx[NUM_of_XAXIS-1][j]-hx[NUM_of_XAXIS-1][j-1]+hy[NUM_of_XAXIS-2][j]-hy[NUM_of_XAXIS-2][j-1]);             1s.2z[B~  
            //ez[NUM_of_XAXIS-1][j]=ez_x_high_m1[j]+ce*(ez[NUM_of_XAXIS-2][j]-ez[NUM_of_XAXIS-1][j]); xH*X5?  
            ez_x_high_m1[j]=ez[NUM_of_XAXIS-2][j]; -%yrs6  
            } m!%aB{e  
thJ~* 0^  
            // y=0的左边边界 Yjjh}R#  
            for(i=1;i<NUM_of_XAXIS-1;i++) <R@,wzK  
            { kc^,V|Nbq6  
                ez[0]=ez_y_low_m1+ce*(ez[1]-ez[0]); ]$lt  
                ez_y_low_m1=ez[1];     rjK`t_(=  
            } u7[}pf$}  
             pCh v;  
            // y=NUM_of_YAXIS-1的右边界 ]9y\W}j  
            for(i=1;i<NUM_of_XAXIS-1;i++) ={ms@/e/T  
            { *|dr-e_j  
                ez[NUM_of_YAXIS-1]=ez_y_high_m1+ce*(ez[NUM_of_YAXIS-2]-ez[NUM_of_YAXIS-1]); Q&I`uS=F  
                ez_y_high_m1=ez[NUM_of_YAXIS-2]; 0Z8"f_GK  
             
,u2Qkw  
            } [rT.k5_  
87 s*lS  
            //四个边角 ciS,  
            ez[0][0]=ez_x_low_m1[0]+cg*(ez[1][1]-ez[0][0]); WrGnLE kiV  
            ez_x_low_m1[0]=ez[1][1]; 1Q>nS[  
|M?vFF]TN  
            ez[0][NUM_of_YAXIS-1]=ez_x_low_m1[NUM_of_YAXIS-1]+ce*(ez[0][NUM_of_YAXIS-2]-ez[0][NUM_of_YAXIS-1]); 'q)g,2B%  
            ez_x_low_m1[NUM_of_YAXIS-1]=ez[0][NUM_of_YAXIS-2]; _5-h\RB)  
Y"MHs0O5>  
            ez[NUM_of_XAXIS-1][0]=ez_x_high_m1[0]+ce*(ez[NUM_of_XAXIS-2][1]-ez[NUM_of_XAXIS-1][0]); =,(TP  
            ez_x_high_m1[0]=ez[NUM_of_XAXIS-2][1]; /qX?ca1_4^  
-U=Ci  
            ez[NUM_of_XAXIS-1][NUM_of_YAXIS-1]=ez_x_high_m1[NUM_of_YAXIS-1]+ce*(ez[NUM_of_XAXIS-2][NUM_of_YAXIS-2]-ez[NUM_of_XAXIS-1][NUM_of_YAXIS-1]); }LT&BNZj  
            ez_x_high_m1[NUM_of_YAXIS-1]=ez[NUM_of_XAXIS-2][NUM_of_YAXIS-2]; rR#Ditn^  
?FAI@
4  

//边界Mur ABC #l=yD]tPU  
            //i=0; x的左边边界 y )QLR<wf  
            for(j=1;j<NUM_of_YAXIS-1;j++) qSY\a\.<  
            { /2tA n  
             ez[0][j]=ez_x_low_m1[j]+ce*(ez[1][j]-ez[0][j])-cf*(hx[0][j]-hx[0][j-1]+hy[1][j]-hy[1][j-1]); "#1\uoH  
            //ez[0][j]=ez_x_low_m1[j]+ce*(ez[1][j]-ez[0][j]); )bUnk+_  
            ez_x_low_m1[j]=ez[1][j]; Rb6BY-/J  
            } (GJtTp~2C4  
            // x右边 的边界 YYWD\Y`8  
            for(j=1;j<NUM_of_YAXIS-1;j++) FZ"n6hWA  
            { @{Rb]d?&F?  
            ez[NUM_of_XAXIS-1][j]=ez_x_high_m1[j]+ce*(ez[NUM_of_XAXIS-2][j]-ez[NUM_of_XAXIS-1][j])-cf*(hx[NUM_of_XAXIS-1][j]-hx[NUM_of_XAXIS-1][j-1]+hy[NUM_of_XAXIS-2][j]-hy[NUM_of_XAXIS-2][j-1]);             F%w!I 9  
            //ez[NUM_of_XAXIS-1][j]=ez_x_high_m1[j]+ce*(ez[NUM_of_XAXIS-2][j]-ez[NUM_of_XAXIS-1][j]); j-<-!jTd  
            ez_x_high_m1[j]=ez[NUM_of_XAXIS-2][j]; TM)INo^  
            } Y%KowgP\  
$vs],C"pX  
            // y=0的左边边界 ;UQGi}?CD  
            for(i=1;i<NUM_of_XAXIS-1;i++) Hjm> I'9  
            { 3=O [Q:8  
                ez[0]=ez_y_low_m1+ce*(ez[1]-ez[0]); sscbf  
                ez_y_low_m1=ez[1];     c!EA>:;(<  
            } )|~K&qn`  
             %*D=ni#(sT  
            // y=NUM_of_YAXIS-1的右边界 6}"%>9  
            for(i=1;i<NUM_of_XAXIS-1;i++) Lq(=0U\"P  
            { Qi18q|l8v  
                ez[NUM_of_YAXIS-1]=ez_y_high_m1+ce*(ez[NUM_of_YAXIS-2]-ez[NUM_of_YAXIS-1]); m<CrkKfpG  
                ez_y_high_m1=ez[NUM_of_YAXIS-2]; (xfc_h*xA  
             *:%&z?<Fw  
            } qc8Ge\3s  
iLy^U*yK  
            //四个边角 p6DI7<C<H  
            ez[0][0]=ez_x_low_m1[0]+cg*(ez[1][1]-ez[0][0]); (?z?/4>7<  
            ez_x_low_m1[0]=ez[1][1]; E7X6RB b  
PCT&d)}  
            ez[0][NUM_of_YAXIS-1]=ez_x_low_m1[NUM_of_YAXIS-1]+ce*(ez[0][NUM_of_YAXIS-2]-ez[0][NUM_of_YAXIS-1]); A8% e_XA  
            ez_x_low_m1[NUM_of_YAXIS-1]=ez[0][NUM_of_YAXIS-2]; F2N"aQ&  
4Mt3<W5  
            ez[NUM_of_XAXIS-1][0]=ez_x_high_m1[0]+ce*(ez[NUM_of_XAXIS-2][1]-ez[NUM_of_XAXIS-1][0]); x-%O1frc  
            ez_x_high_m1[0]=ez[NUM_of_XAXIS-2][1]; f_~T  
\{h_i FU!  
            ez[NUM_of_XAXIS-1][NUM_of_YAXIS-1]=ez_x_high_m1[NUM_of_YAXIS-1]+ce*(ez[NUM_of_XAXIS-2][NUM_of_YAXIS-2]-ez[NUM_of_XAXIS-1][NUM_of_YAXIS-1]); 7-744wV}Z  
            ez_x_high_m1[NUM_of_YAXIS-1]=ez[NUM_of_XAXIS-2][NUM_of_YAXIS-2]; 2F4<3k!&  
1Lb)S@Q`*R  
请问楼主  为什么四条边的吸收边界条件所用公式不一样啊  ？

感谢分享！呵呵……

多谢，下了学习一下