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

这个是EX的迭代程序 计算的是填充部分介质波导 ，我想问一下 这个分裂场的写法是否正确呢 ？因为之前做散射时 这样的写法 结果也偏小，但是简单易懂 g;Ugr8  
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) I'xc$f_+  
{ Rxdj}xy  
 int i,j,k,l; -W!M:8  
 double pd1,pd2; KTYjC\\G  
    for(i = ls-1;i<=meshnum_x-ls;i++) X>$Wf
3  
 { $6m@gW]N  
  for(j = ls;j<=meshnum_y-ls;j++) +E. D:  
  { bIm4s  
   for(k = ls;k<=meshnum_z-ls;k++) XW aa`q  
   { u^xnOVE  
    pd1=0; UG\2wH_  
    pd2=0; g7323m1=  
    for(l = 0;l<= 2*ls-1;l++) 0j8fU7~6S  
    { igL5nE=n  
     pd1 += lag[1][l]*hz[j-ls+l][k]; // 高阶部分  计算的一阶导数 9Qszr=C0  
     pd2 += lag[1][l]*hy[j][k-ls+l]; l^I?@{W  
    } YC;@^  
    pd1=pd1/dy; A$p&<#  
    pd2=pd2/dz; r)~?5d  
    xD*Zcw(vj~  
      exy[j][k]=exy[j][k] + d*dt/epsilon[j][k]*pd1;    //分裂的部分 {ccc[G?>.Q  
    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; 25Ee+&&%  
                ex[j][k] = exy[j][k] +exz[j][k]; ?5't1219  
   }
!{*yWpZ:  
  .. XOX$uLm  
dWzf C@]  

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

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

S11曲线是否吻合？