[求助]为什么一有sigma，结果就非常非常小，正常吗？

这个程序是DENNIS M.SULLIVAN中的真空三维程序，我把它改为matlab语言写的。我在真空的基础上改了频率，给sigma赋了值，相应的波长网格划分也改了，程序如下，很简单的程序。 A&d_!u>  
?i\$U'2*z3  
clear all; "Z"`X3,-z  
clc; rm<`H(cT  
ie=60; sDwE,f0h  
je=60; /AIFgsaY  
ke=60; >?{> !#1  
R_IT${O  
ib=ie+1;     G&yF9s)Lvs  
jb=je+1;   j[=P3Z0q  
kb=ke+1;   ^
C8f(  
yg%T{hyzH  
cc=3*10^8; #b1/2=PA  
S+pm@~xe  
ex=zeros(ie,jb,kb); @ 3b-  
ey=zeros(ib,je,kb); InG<B,/
W?  
ez=zeros(ib,jb,ke); i;67<f}-  
hx=zeros(ib,je,ke); L @8[.  
hy=zeros(ie,jb,ke); !
/"y  
hz=zeros(ie,je,kb); rjHW  

xQ=L2pX  
dx=5.995849160000000e+006;  %2.5hz对应的波长的20分之一 i;#AW($+a  
dt=0.9/sqrt(3.0./dx/dx)/cc; !)OA7%3m  
&RlYw#*1.  
ic=ie/2; ([h
d  
jc=je/2; /j:-GJb*!u  
kc=ke/2; } ejc  
]Wc 2$  
epsz=8.85e-12; H
"rIOoxf  
eps=1.0; w.Kp[  
sig=1.0; "'c =(P  
freq=2.5; C2?p>S/q  
muz=4.0*pi*1.0e-7;          %permeability of free space rzKn5
Z  
d`KW]HJw  
%specify the dipole 0urM@/j+  
B)ibxM(n*  
ca=(1.0-(dt*sig)/(2.0*epsz*eps))/(1.0+(dt*sig)/(2.0*epsz*eps)); oX-h7;SD  
cb=(dt/epsz/eps/dx)/(1.0+(dt*sig)/(2.0*epsz*eps)); 3UGdXufw  
da=1.0; HT
P~5J  
db=dt/muz/dx; t<qXXQ&5  
<gSZ<T  
tau= 0.2;   %高斯波的脉宽，经过计算后的值（不是随便取的） N4,!b_1  
t0=4.5*tau/dt;   %t0也是计算后的值 `"|u NVn  
spread=tau/dt;  ePI)~  
T=0; E"%G@,|3*  
nsteps=150; ;4(FS  
/7b$C]@k  
for n=1:1:nsteps -KCQ!0\F  
    T=T+1
}/z\%Y  
     }x
:nhy`  
ex(1:ie,2:je,2:ke)=ca*ex(1:ie,2:je,2:ke)+... E
JO6k1  
                   cb*(hz(1:ie,2:je,2:ke)-hz(1:ie,1:je-1,2:ke)+... ;Z#DB$o\  
                       hy(1:ie,2:je,1:ke-1)-hy(1:ie,2:je,2:ke)); D^2yP~(  
7A>glZ/x  
ey(2:ie,1:je,2:ke)=ca*ey(2:ie,1:je,2:ke)+... ;gW?Fnry;  
                   cb*(hx(2:ie,1:je,2:ke)-hx(2:ie,1:je,1:ke-1)+... 0Yc
#fD  
                       hz(1:ie-1,1:je,2:ke)-hz(2:ie,1:je,2:ke)); mr`EcO0  
                     dRLvej,  
ez(2:ie,2:je,1:ke)=ca*ez(2:ie,2:je,1:ke)+... q*{Dy1Tj  
                   cb*(hx(2:ie,1:je-1,1:ke)-hx(2:ie,2:je,1:ke)+... xW'(]Z7_  
                       hy(2:ie,2:je,1:ke)-hy(1:ie-1,2:je,1:ke)); C9cQ} j:  
                     M{~KT3c  
   pulse=exp(-((T-t0)/spread)^2); t/:]\|]WB  
   ,m.IhnCV\  
   ez(ic,jc,kc)=pulse; Z+x`q#ZQr  
%*********************************************************************** 4+r26S,T  
%     Update magnetic fields p"/B3  
%*********************************************************************** z,=k F I  
wz5*?[4  
hx(2:ie,1:je,1:ke)=hx(2:ie,1:je,1:ke)+... 0t}&32lL&  
                   db*(ey(2:ie,1:je,2:kb)-ey(2:ie,1:je,1:ke)+... jiAN8t*P  
                       ez(2:ie,1:je,1:ke)-ez(2:ie,2:jb,1:ke)); ,46k8%
WW  
                 !4 G9`>n  
hy(1:ie,2:je,1:ke)=hy(1:ie,2:je,1:ke)+... eI*o9k$Qs  
                   db*(ex(1:ie,2:je,1:ke)-ex(1:ie,2:je,2:kb)+... qQUCK  
                       ez(2:ib,2:je,1:ke)-ez(1:ie,2:je,1:ke)); 38eeRo  
                 421ol  
hz(1:ie,1:je,2:ke)=hz(1:ie,1:je,2:ke)+... OO+QH 2j  
                   db*(ex(1:ie,2:jb,2:ke)-ex(1:ie,1:je,2:ke)+... )k|_ CW~  
                       ey(1:ie,1:je,2:ke)-ey(2:ib,1:je,2:ke)); n6 a=(T  
   qU,u(El  
smooth(ez(:,45,30))   ;   RgT|^|ZA  
figure(1); [!#<nY/C  
subplot(1,2,1) {QTnVS't 0  
surf(ez(:,:,30)); (lVHKg&U[  
colorbar; ln6Hr^@5  
pause(0.1); 2Fsv_t&*>  
disp('ez()='); Q'S"$^~{  
disp(ex(45,45,30)); o Kfm=TbY  
AA(T)=ex(45,45,30); )R'~{;z }  
')>&:~  
\`M8Mu9~w  
subplot(1,2,2) _}-Ed,.=  
imagesc(ez(:,:,30)); NJraol  
colorbar; W{(q7>g  
pT`oC&  
[q{[Avqf  
end s!j[Ovtx  
^p433  
T=(1:1:nsteps); ("B[P/  
figure(2); h|j
$Jy  
plot(T,AA .. h+7THMI  
oMV<Yn_<  

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

2.5hz的电磁波

请问理论上，当频率很低时，能这样用FDTD吗，这样有没有什么问题呢？还是说就是这样。。。。频率很低，因为有的地方只能用低频，比如海洋测井等等。所以用2.5hz。为什么一有sigma,值就超级小，正常吗？ ^xyU*A}D  
谢谢！！

要知道FDTD是是差分的麦克斯维方程得到的。你的2.5GHZ电磁波存在么，或者说你的2.5GHZ的波是不是符合麦克斯维方程的别变化规律，这个我不知道，不好乱讲。