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

在加入一维入射波时，总场边界上任意点投影到一维波入射方向k上的垂足到源点的距离为d， iXRt9)MT{  
x^2/jUc#B  
在离散化后，为什么是 v6\2mc.  
dRa<,@1"  
d = dx*cos(phi)*cos(phi) + dy*sin(phi)*sin(phi)     ，其中phi为一维波入射角。 gDNW~?/  
66^t[[  
我十分不明白为什么是这样？ ^)l@7XxD  

q"<-  
参看了Allen Taflove 的computational electromagnetic：FDTD书中的公式, chapter 5, P197. nk 9 K\I  
!\O,dq  
d = k dot r_comp, 则d至多也就是： d = dx*cos(phi)+ dy*sin(phi) y03a\K5[KQ  
U&$]?3?  

{~g  
=Gz>ZWF  
b'5pQ2Mq  
}{oBKm9_p  


9z/_`Xd_  
86F+N_>Z  

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

自己顶一个，以免沉了。

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

入射波搞定了没你

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

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

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

晕，你理论没跟上啊。

我是从别人的程序看来的。下面是平面波高斯正弦源的成员函数设置： )"7iJb<E  
Pd_U7&w,5  
int CFDTD_2D_TM::Init_plWave_GaussSin(double E0, double omega, double t0, double tw, 8}O lL,fP  
                            double PHI, int nxa, int nxb, int nya, int nyb) 0AL=S$B)
 
{ ivJ@=pd)B  
    double phase = 0.0; +~ P2C6@G  
     lR6@ xJd:@  
    double eps_r = Mat[0][0]; qm/22:&v5  
    double mu_r  = Mat[0][1]; ek*rp`y]  
    double sigma = 0.0; @niHl  
rlOAo`hd  
    //the incidence angle of the plane wave Rl?_^dPx  
    phi = PHI;                                    //入射角 D2Kp|F;  
tEvut=k'  
    //the size of the total field zone *0Skd  
    l_x = nxb - nxa + 1;                      //x方向总场的长度 12b(A+M   
    l_y = nyb - nya + 1;                      //y方向总场的长度 L [pBB  
    cos_phi = cos(phi); _-K2/6zy  
    sin_phi = sin(phi);  iu=7O  
     *k.G5>@  
    l_1D = (int) (l_x*abs(cos_phi) + l_y*abs(sin_phi)) + 10 + 2*n_PML;  //一维入射波计算的网格数 8;RUf
~q?  
    //l_1D = nx; 8V`WO6*  
     rZ}:Z'`  
    d = dx*cos_phi*cos_phi + dy*sin_phi*sin_phi;  //这就是我的问题所在？ KPKt^C  
qN9(S:_Px  
    fdtd_1D.Init(l_1D, n_PML, dt, d, eps_r,  mu_r, sigma); 3u+T~g0^  
     0z6R'Kjy A  
    fdtd_1D.Init_PML_Param(); f<d`B]$(  
:!WHFB o 8  
    fdtd_1D.Init_GaussSin(E0, omega, phase, t0, tw); I-]?"Q7Jz  
     p>v$FiV2N  
    fdtd_1D.Get_Data(Ez_1D,Hy_1D); 7x|9n  
     T $>&[f$6  
    n_TS_xa = nxa; $r@zs'N  
    n_TS_xb = nxb; Y!w`YYKP  
    n_TS_ya = nya; ; F"g$_D0  
    n_TS_yb = nyb; f@wquG'  
     -b9\=U[  
    jel_plane_wave = 1; t%/&c::(6  
84& $^lNV  
rr],DGg+B]  
    if (Init_TS_Param()) /T0F"e)Ci  
    { U]H#MiC!  
        // memory allocation problems wDal5GJp  
        return 1; l[0RgO*S  
    } PR#exm&  
     +>6iYUa  
    return 0; 7rc0yB  
} {HltvO%8  
_Xe>V0  
823Y\x~>  
Tztu}t]N  
KOk4^#h@  
LM<qT-/qs  
/K@XzwM  
//////////////////////////////////////////////////////////////////////////////////////

按照allen书上的公式： H(ARw'M  
X:"i4i[}{9  
d = k dot r_comp ^1I19q  
we//|fA<  
   = (i_comp - i_0)*cos(phi)  + (j_comp - j_0)*sin(phi) {T$9?`h~M  
]vUwG--*  
可是在我上面贴出的部分代码中为何d = dx*cos_phi*cos_phi + dy*sin_phi*sin_phi; v!~fs)cdE|  
{I%cxQ#y  
这让我很迷惑，也不知道怎么推导？