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

在加入一维入射波时，总场边界上任意点投影到一维波入射方向k上的垂足到源点的距离为d， [HF)d#A  
PNT.9 *d  
在离散化后，为什么是 @`}'P115@  
{eIE|  
d = dx*cos(phi)*cos(phi) + dy*sin(phi)*sin(phi)     ，其中phi为一维波入射角。 PZLWyp  
] 5
P{*  
我十分不明白为什么是这样？ i/O!bq[o  
zQ]IlMt  
参看了Allen Taflove 的computational electromagnetic：FDTD书中的公式, chapter 5, P197. t
l;?/  
-2> L*"^  
d = k dot r_comp, 则d至多也就是： d = dx*cos(phi)+ dy*sin(phi) Uo^s]H#:  
8i-?\VZD  
TW3
:Y\p  
@pytHN8( $  
1{o CMq/v  
"R*B~73  

HH^{,53%  
_?kf9.  

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

自己顶一个，以免沉了。

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

入射波搞定了没你

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

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

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

晕，你理论没跟上啊。

我是从别人的程序看来的。下面是平面波高斯正弦源的成员函数设置： 9.'h^#C  
+#U|skl  
int CFDTD_2D_TM::Init_plWave_GaussSin(double E0, double omega, double t0, double tw, &Z(K6U#.  
                            double PHI, int nxa, int nxb, int nya, int nyb) **9x?s  
{ =4V&*go*\  
    double phase = 0.0; ZkL8
e  
     pl@O N"=[  
    double eps_r = Mat[0][0];
NBl+_/2'w  
    double mu_r  = Mat[0][1]; j',W 64  
    double sigma = 0.0; glDcUCF3  
P-F)%T[  
    //the incidence angle of the plane wave Z;9>S=w!  
    phi = PHI;                                    //入射角 |4$M]Mf0  
]Chj T}  
    //the size of the total field zone &'cL%.  
    l_x = nxb - nxa + 1;                      //x方向总场的长度 fjvN$NgVs  
    l_y = nyb - nya + 1;                      //y方向总场的长度 T134ZXqqz  
    cos_phi = cos(phi); ah
x>q  
    sin_phi = sin(phi); L,y6^J!  
     Gjo
Im?  
    l_1D = (int) (l_x*abs(cos_phi) + l_y*abs(sin_phi)) + 10 + 2*n_PML;  //一维入射波计算的网格数 `E+Jnu,jC  
    //l_1D = nx; f}g\D#`]/  
     Lg8nj< TF  
    d = dx*cos_phi*cos_phi + dy*sin_phi*sin_phi;  //这就是我的问题所在？ s { #3r  
SJD@&m%?[  
    fdtd_1D.Init(l_1D, n_PML, dt, d, eps_r,  mu_r, sigma); w=b)({`M  
     4tL<q_  
    fdtd_1D.Init_PML_Param(); afjtn_IB  
U*90m~)  
    fdtd_1D.Init_GaussSin(E0, omega, phase, t0, tw); zvABU+{jD  
     720PjQ  
    fdtd_1D.Get_Data(Ez_1D,Hy_1D); F, U*
yj  
     "$s~SIUB  
    n_TS_xa = nxa; oFOnjK"|F  
    n_TS_xb = nxb; B8E'ddUw  
    n_TS_ya = nya; ?X@fKAj  
    n_TS_yb = nyb; UOQEk22  
     c/c$D;T  
    jel_plane_wave = 1; >WY#4  
~l;[@jsw F  
zdXkR]  
    if (Init_TS_Param()) [MuEoWrq(}  
    { dJ"iEb|4  
        // memory allocation problems /mo(_  
        return 1; oIUy-|  
    } x.Q&$#  
     U qG .:@T  
    return 0; ~y=T5wt  
} 3u%{dGa  
Gt9&)/#  
/cc\fw1+  
+P.+_7+:  
]Q[p@gLd  
GfM;saTz{  
X,9 M"E 2  
//////////////////////////////////////////////////////////////////////////////////////

按照allen书上的公式： =8U&[F  
Q:J^"  
d = k dot r_comp v4@Z(M  
n3J53| %v  
   = (i_comp - i_0)*cos(phi)  + (j_comp - j_0)*sin(phi) V5}B:SUB  
^eW}XRI  
可是在我上面贴出的部分代码中为何d = dx*cos_phi*cos_phi + dy*sin_phi*sin_phi; OY?y^45y  
'X shmZ0&  
这让我很迷惑，也不知道怎么推导？