[求助]请教一维平面波加入的问题？

在加入一维入射波时，总场边界上任意点投影到一维波入射方向k上的垂足到源点的距离为d， W)Y:2P<.  
r5+MjR  
在离散化后，为什么是 woN d7`C}7  
QUF1_Sa  
d = dx*cos(phi)*cos(phi) + dy*sin(phi)*sin(phi)     ，其中phi为一维波入射角。 ''k}3o.K[  
;!OME*?m<  
我十分不明白为什么是这样？ ZmI#-[/  
buWF6LFC  
参看了Allen Taflove 的computational electromagnetic：FDTD书中的公式, chapter 5, P197. +eop4
|Z  
\lyHQ-gWhc  
d = k dot r_comp, 则d至多也就是： d = dx*cos(phi)+ dy*sin(phi) <l>L8{-3  
?ZkVk=t?  
w;
J#+ik  
'C;KNc  
RLGIST`  
%WYveY  

GkKoc v  
QqcAmp  

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

自己顶一个，以免沉了。

不知道你一开始的公式从哪理看来的，allen taflove书的公式没问题啊，这个就是矢量相乘的运算啊，应该不难理解

入射波搞定了没你

～～汗，没有，我检查了程序，确实发现了加减号有错误，可是改了以后，还是不行，我把入射波延z轴引入，把有泄漏的角点手动计算一下跟程序运行的值是一样的，现在一团浆糊

手动计算的和程序运行的一样？什么意思啊，是理论上就该这个值？

我的入射波是延z轴入射，在计算j=je_min=30左侧xz面的ex时，这个连接边界上（90，30，90）与（90，30，30）这两个点是有泄漏，ex场分量计算时根据程序这两个点投影到一维入射波来计算，我把这两个点投影到一维入射波手动计算的值与程序运行的值是相符的，但是我自己也糊涂，这么算到底验证了什么，所以是一塌糊涂，感觉自己还是没搞清

晕，你理论没跟上啊。

我是从别人的程序看来的。下面是平面波高斯正弦源的成员函数设置： Gg!))I+  
xv*mK1e  
int CFDTD_2D_TM::Init_plWave_GaussSin(double E0, double omega, double t0, double tw, -0]aOT--  
                            double PHI, int nxa, int nxb, int nya, int nyb) RhJ<<T.2  
{ T+[e6/|  
    double phase = 0.0; hD,-!R  
     x$Dv&4  
    double eps_r = Mat[0][0]; n;=A'g|Q  
    double mu_r  = Mat[0][1];  jPs+i  
    double sigma = 0.0; -?]ltn9!  
YL.z|{\e  
    //the incidence angle of the plane wave {c1qC zM4  
    phi = PHI;                                    //入射角 U;nC)'~YW9  
Qnw$=L:  
    //the size of the total field zone |HEw~x<=  
    l_x = nxb - nxa + 1;                      //x方向总场的长度 g\l;>  
    l_y = nyb - nya + 1;                      //y方向总场的长度 xL BG}C  
    cos_phi = cos(phi); "5YdmBy  
    sin_phi = sin(phi); N:<O  
     c]]F`B  
    l_1D = (int) (l_x*abs(cos_phi) + l_y*abs(sin_phi)) + 10 + 2*n_PML;  //一维入射波计算的网格数 6Te}"t>  
    //l_1D = nx; n=
&c5!  
     }-vBRY  
    d = dx*cos_phi*cos_phi + dy*sin_phi*sin_phi;  //这就是我的问题所在？ &8pXkD#A  
_R<HC  
    fdtd_1D.Init(l_1D, n_PML, dt, d, eps_r,  mu_r, sigma); ]M2<I#hF.  
     !lF^~x  
    fdtd_1D.Init_PML_Param(); z(WpOD  
Q*8-d9C  
    fdtd_1D.Init_GaussSin(E0, omega, phase, t0, tw); hG@ys5  
     `[KhG)Y7t  
    fdtd_1D.Get_Data(Ez_1D,Hy_1D); 2tv40(M:<  
     ]Rye AJ3  
    n_TS_xa = nxa; =l(JJ  
    n_TS_xb = nxb; Kzb@JBIF  
    n_TS_ya = nya; ["F,|e{y$  
    n_TS_yb = nyb; (*ng$zZ$  
     v2{O67j} o  
    jel_plane_wave = 1; NM@An2  
FNuu',:  
mig3.is  
    if (Init_TS_Param()) k)s 7Ev*  
    { ddP,_.0  
        // memory allocation problems 2FzS_\":I  
        return 1; (Pd>*G\  
    } s]c$]&IGG  
     j 7URg>i0  
    return 0; f 99PwE(=  
} m'tk#C  
0I((UA/7Zs  
1 =<|h  
{9".o,  
'EV *-_k  
.LhmYbQ2WE  
?U08A{ c  
//////////////////////////////////////////////////////////////////////////////////////

按照allen书上的公式： FH8mK)  
lM^!^6=v0l  
d = k dot r_comp h(Ed%  
%uVJLz  
   = (i_comp - i_0)*cos(phi)  + (j_comp - j_0)*sin(phi) -5 /v`  
i8_x1=A  
可是在我上面贴出的部分代码中为何d = dx*cos_phi*cos_phi + dy*sin_phi*sin_phi; Gr$*t,ZW  
M#]|$\v(  
这让我很迷惑，也不知道怎么推导？