[求助]各位牛人看一下三维的迭代程序哪出错了 计算的值偏小一个数量级

这个是EX的迭代程序 计算的是填充部分介质波导 ，我想问一下 这个分裂场的写法是否正确呢 ？因为之前做散射时 这样的写法 结果也偏小，但是简单易懂 lG\lu'<C  
void computex(double ***ex,double ***exy,double ***exz,double ***hy,double ***hz,double ***sy,double *sz1,double ***epsilon,double dt,double dy,double dz,int meshnum_x,int meshnum_y,int meshnum_z,int ls,double **lag,double d) 4yRT!k}o  
{ rZJp>Q)s  
 int i,j,k,l; YI%7#L7C  
 double pd1,pd2; C!qW:H  
    for(i = ls-1;i<=meshnum_x-ls;i++) F{l,Tl"Jw  
 { |8+<qgQ  
  for(j = ls;j<=meshnum_y-ls;j++) \hi{r@
k>}  
  { c0Q`S"o+  
   for(k = ls;k<=meshnum_z-ls;k++) v#-%_V>ph  
   { yaR|d3ef?4  
    pd1=0; {zVJlJKxs  
    pd2=0; aMv  
    for(l = 0;l<= 2*ls-1;l++) *ZN"+wf\  
    { QREIr |q'  
     pd1 += lag[1][l]*hz[j-ls+l][k]; // 高阶部分  计算的一阶导数 R0. `2=  
     pd2 += lag[1][l]*hy[j][k-ls+l]; YXV![gw0  
    } *ArzXhs[  
    pd1=pd1/dy; >#!n"i;  
    pd2=pd2/dz; ),
M8W15  
    i"JF~6c<  
      exy[j][k]=exy[j][k] + d*dt/epsilon[j][k]*pd1;    //分裂的部分 ,{'~J @  
    exz[j][k]=func_r1(d*dt*sz1[k]/epsilon[j][k])*exz[j][k] - d*dt/epsilon[j][k]*func_s1(d*dt*sz1[k]/epsilon[j][k])*pd2; O,^,G<`  
                ex[j][k] = exy[j][k] +exz[j][k]; 1-w1k^e  
   } z W+wtYV4  
  .. '?3Hy|}  
zw5EaY  

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

这个按照葛德彪书上的 分裂写的 但最后得到波导的S11参数 值小了1个数量级(10),我检查了看到反射波小了,估摸着是这里合起来写出现的问题（即PML区 与空间传播区和在一起写程序),请各位帮帮忙！

S11曲线是否吻合？