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

在进行FDTD计算时，介质的四个参数的设置问题： '`>%RZ]  
1、真空的绝对参数 s L^+$Mq6  
double Epsilon0   = 8.854e-12；                // ε0  真空的介电系数，法拉/米  (F/m) 6"&cQ>$xh  
double Mu0            =  1.25663706e-6;        // μ0  真空的磁导率  : PI*4.0e-7，亨利/米  (H/m) YX$(Sc3.6  
double Sigma0     = 2.0e-4;                                           // σ0  真空的电导率，西门子/米(S/m)，空气为 0 yHeL&H  
double SigmaM0  = 0;                                                      // σm0 真空的导磁率，欧姆/米  (Ω/m) $ev+0m_  
2、介质的相对参数 Bqf(6\)F  
   Epsilon[0]   = 1                //真空的相对介电常数 w*F[[*j@.  
   Epsilon[1]   = 10              //  理想金属的相对介电常数 1G>Ud6(3<  
   6{h\CU}"  
   Mu[0]        = 1                       // 真空的相对磁导率    μ  磁导率 GG%b"d-  
   Mu[1]        = 1                       // 铜为抗磁体，相对磁导率为0.999999 [:8\F#KW  
  @<\oM]jX  
   Sigma[0]     = 1                          // 真空中 电导率：2.0e-4，空气中为0               ^O07GYF  
   Sigma[1]     = 5.8e+7               // .. y]db]pP5  
HpnF,4A>  

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

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

现在好一些了吗？

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

   Epsilon[0]   = 1 n7i~^nf>  
   Epsilon[1]   = 10 ] K&ca  
   Mu[0]        = 1                   H.M: cD:  
   Mu[1]        = 1                     T$mT;k  
   Sigma[0]     = 1               Fep@VkN  
   Sigma[1]     = 5.8e+7     _S[@?]=`b  
   SigmaM[0]    = 0           Zz (qc5o,F  
   SigmaM[1]    = 0

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

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

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

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