[讨论]再讨论介质参数设置问题，请高手指点^_^

在进行FDTD计算时，介质的四个参数的设置问题： n5z|@I`S_  
1、真空的绝对参数 F]fXS-@ c  
double Epsilon0   = 8.854e-12；                // ε0  真空的介电系数，法拉/米  (F/m) z,bK.KFSs  
double Mu0            =  1.25663706e-6;        // μ0  真空的磁导率  : PI*4.0e-7，亨利/米  (H/m) 4)8e0L*[B?  
double Sigma0     = 2.0e-4;                                           // σ0  真空的电导率，西门子/米(S/m)，空气为 0 |AT`(71  
double SigmaM0  = 0;                                                      // σm0 真空的导磁率，欧姆/米  (Ω/m) 4;Vi@(G)  
2、介质的相对参数 FmhAUe  
   Epsilon[0]   = 1                //真空的相对介电常数 10)RLh|+  
   Epsilon[1]   = 10              //  理想金属的相对介电常数 T0~~0G)k  
   =sAU5Ag68  
   Mu[0]        = 1                       // 真空的相对磁导率    μ  磁导率 Fy
0sn|  
   Mu[1]        = 1                       // 铜为抗磁体，相对磁导率为0.999999 t1G__5wp  
  hxMV?\MYj  
   Sigma[0]     = 1                          // 真空中 电导率：2.0e-4，空气中为0               56
MY@  
   Sigma[1]     = 5.8e+7               // .. m41%?uC/  
r-*j"1 e  

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

计算代码最好改得简洁一些，这样比较好找错误 GlTpK^
.  
把代码中系数用变量代换，这样代码清晰

现在好一些了吗？

没有改变啊，还是以前代码

   Epsilon[0]   = 1 v5\ALWy+p  
   Epsilon[1]   = 10 1aS:bFi`  
   Mu[0]        = 1                   nlhv  
   Mu[1]        = 1                     WgR%mm^  
   Sigma[0]     = 1               2HBey  
   Sigma[1]     = 5.8e+7     de2G"'F  
   SigmaM[0]    = 0           f/8&-L  
   SigmaM[1]    = 0

你的数据到底是放在变量Epsilon0,Epsilon1里面，还是放在数组Epsilon[0],Epsilon[1]里面了？ ndDF(qHr  
不解，感觉很乱啊。

是的，介质的绝对参数是放在变量里的，而介质的相对参数是放在数组里的，现在假设只有两种介质，真空和理想导体，不好意思 <VstnJo`Z  
不知道这样设置还有这些值及CA CB CP CQ的计算对不对？

用Fortran吧，葛老师书后面有现成的程序可以抄。。

书上的例子看过了，是二维的，我要做的是三维的，只能参考一下了，呵呵