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

在加入一维入射波时，总场边界上任意点投影到一维波入射方向k上的垂足到源点的距离为d， #Fzb8Yo  
FG#nap{  
在离散化后，为什么是 hS_.l}0yf  
vJThU$s-  
d = dx*cos(phi)*cos(phi) + dy*sin(phi)*sin(phi)     ，其中phi为一维波入射角。 8A4TAT4,  
PWG;&ma  
我十分不明白为什么是这样？ }!B.K^@)  
U8[Qw}T P  
参看了Allen Taflove 的computational electromagnetic：FDTD书中的公式, chapter 5, P197. ?(Bl~?zD  
{u4i*udG`)  
d = k dot r_comp, 则d至多也就是： d = dx*cos(phi)+ dy*sin(phi) '!@A}&]  
qc&jd  
?i%nMlcc  
FfSI n3  
A:)sg!Lt  
0_HJ.g!  

JMePI%#8  
UN(3i(d  

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

自己顶一个，以免沉了。

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

入射波搞定了没你

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

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

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

晕，你理论没跟上啊。

我是从别人的程序看来的。下面是平面波高斯正弦源的成员函数设置： k2a^gCBC  
mbK$Wp#  
int CFDTD_2D_TM::Init_plWave_GaussSin(double E0, double omega, double t0, double tw, 2 r)c?  
                            double PHI, int nxa, int nxb, int nya, int nyb) "
~6BC  
{ U
gJHSl  
    double phase = 0.0; 7;V5hul  
     BDg /pDnwg  
    double eps_r = Mat[0][0]; 2J5dZYW  
    double mu_r  = Mat[0][1]; J/ ~]A1fP6  
    double sigma = 0.0; 'Z[R*Ikzq  
Z9y:}:j"  
    //the incidence angle of the plane wave Kk#@8h>  
    phi = PHI;                                    //入射角 Hqk2W*UTl  
&-9D.'WzP  
    //the size of the total field zone Q*5d~Yr]R  
    l_x = nxb - nxa + 1;                      //x方向总场的长度 -_dgd:or  
    l_y = nyb - nya + 1;                      //y方向总场的长度 I&9B^fF6  
    cos_phi = cos(phi); Lj#6K@u@Z  
    sin_phi = sin(phi); VrG|/2  
     0%GWc}o  
    l_1D = (int) (l_x*abs(cos_phi) + l_y*abs(sin_phi)) + 10 + 2*n_PML;  //一维入射波计算的网格数 '_%Jw:4k  
    //l_1D = nx; {z7{ta  
     fr7/%{s  
    d = dx*cos_phi*cos_phi + dy*sin_phi*sin_phi;  //这就是我的问题所在？ -=Hr|AhE  
 u[u=:Y+  
    fdtd_1D.Init(l_1D, n_PML, dt, d, eps_r,  mu_r, sigma); :0l(Ll KD  
     X,b}d#\  
    fdtd_1D.Init_PML_Param(); l%O-c}X  
e# DAa  
    fdtd_1D.Init_GaussSin(E0, omega, phase, t0, tw); (`N/1}vk  
     n >^?BU  
    fdtd_1D.Get_Data(Ez_1D,Hy_1D); Ra5cfkH;  
     _nt%&f  
    n_TS_xa = nxa; 6r`g+Js/  
    n_TS_xb = nxb; \GN5Sy]r  
    n_TS_ya = nya; U7N<!6  
    n_TS_yb = nyb; @QofsWC  
     VI4d/2e  
    jel_plane_wave = 1; _<8y^y
mo  
:>5]A6Wi  
;-F#a+2]!  
    if (Init_TS_Param()) iT5%X  
    { i.]}ooI  
        // memory allocation problems bP[/  
        return 1; ??Lxb% 7R  
    } @NF8?>!  
     }`X$ '  
    return 0; w K+2;*bI  
} )8_0d)  
D'X'h}+2  
YY9q'x,w  
u'C4d6\wS  
VmW_,  
n.)-aRu[  
*w;f\zW  
//////////////////////////////////////////////////////////////////////////////////////

按照allen书上的公式： m-:k]9I  
cGDA0#r  
d = k dot r_comp (?Yz#Yf  
Hi9]M3Ub  
   = (i_comp - i_0)*cos(phi)  + (j_comp - j_0)*sin(phi) N\u-8nE5  
er_aol e  
可是在我上面贴出的部分代码中为何d = dx*cos_phi*cos_phi + dy*sin_phi*sin_phi; $t]DxMd  
SaNx;xgi  
这让我很迷惑，也不知道怎么推导？