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

这个程序是DENNIS M.SULLIVAN中的真空三维程序，我把它改为matlab语言写的。我在真空的基础上改了频率，给sigma赋了值，相应的波长网格划分也改了，程序如下，很简单的程序。 c0 WFlj9b  
(}:C+p 'I  
clear all; &gc`<kLu  
clc; Vdn.)ir~P  
ie=60; &D#+6M&LK{  
je=60; G*\h\
@  
ke=60; :o"8MZp  
~@8+hnE]  
ib=ie+1;     yW.COWL=)  
jb=je+1;   K<v:-TjQZ:  
kb=ke+1;   FXo2Y]K3`L  
P.sg
RsL  
cc=3*10^8; J`#`fX  
x:t<ZG&Xwg  
ex=zeros(ie,jb,kb); s T3p
>8n  
ey=zeros(ib,je,kb); 0W>9'Rw  
ez=zeros(ib,jb,ke); NfE.N&vI_c  
hx=zeros(ib,je,ke); dZ:r&Qa  
hy=zeros(ie,jb,ke); %M
cO6.M@  
hz=zeros(ie,je,kb); ^^(<c,NX#M  
Ak~4|w-  
dx=5.995849160000000e+006;  %2.5hz对应的波长的20分之一 W>spz~w%j  
dt=0.9/sqrt(3.0./dx/dx)/cc; fn9#>~vrD  
uG

>nV  
ic=ie/2; V)D-pV V  
jc=je/2; @<_`2eW'/R  
kc=ke/2; 83.E0@$  
w5]l1}rl  
epsz=8.85e-12; g$j6n{Yl  
eps=1.0; AH.9A_dG  
sig=1.0; KIL18$3J  
freq=2.5; wUi(3g|A  
muz=4.0*pi*1.0e-7;          %permeability of free space 3z{S}~  
p-GT`D  
%specify the dipole FoPginZ]J  
U '[?9/T  
ca=(1.0-(dt*sig)/(2.0*epsz*eps))/(1.0+(dt*sig)/(2.0*epsz*eps)); bO49GEUT _  
cb=(dt/epsz/eps/dx)/(1.0+(dt*sig)/(2.0*epsz*eps)); jiIST^Zq#t  
da=1.0; uC~g#[I QM  
db=dt/muz/dx; )%du@a8  
.F0]6#(  
tau= 0.2;   %高斯波的脉宽，经过计算后的值（不是随便取的） H]&!'\aUz  
t0=4.5*tau/dt;   %t0也是计算后的值 W'_/6_c$!  
spread=tau/dt; .V%*{eHLL  
T=0; hbuZaxo<  
nsteps=150; 
c!&Qj  
R V!o4"\]  
for n=1:1:nsteps C]`eH*z~8  
    T=T+1 DM3B]Yl  
     GH'O!}  
ex(1:ie,2:je,2:ke)=ca*ex(1:ie,2:je,2:ke)+... !v!N>f4S$  
                   cb*(hz(1:ie,2:je,2:ke)-hz(1:ie,1:je-1,2:ke)+... SZVV40w  
                       hy(1:ie,2:je,1:ke-1)-hy(1:ie,2:je,2:ke)); b2h":G|s  
LD+f'^>>Z  
ey(2:ie,1:je,2:ke)=ca*ey(2:ie,1:je,2:ke)+... |0{ i9.=  
                   cb*(hx(2:ie,1:je,2:ke)-hx(2:ie,1:je,1:ke-1)+... za,2r^  
                       hz(1:ie-1,1:je,2:ke)-hz(2:ie,1:je,2:ke)); M@?"t_e1  
                     pH9HK  
ez(2:ie,2:je,1:ke)=ca*ez(2:ie,2:je,1:ke)+... NMjnL&P`  
                   cb*(hx(2:ie,1:je-1,1:ke)-hx(2:ie,2:je,1:ke)+... ZHy><=2  
                       hy(2:ie,2:je,1:ke)-hy(1:ie-1,2:je,1:ke)); g=A$<k
 
                     zj]b&In6;  
   pulse=exp(-((T-t0)/spread)^2); b?kPN:U#N/  
   ID
8k/t!  
   ez(ic,jc,kc)=pulse; CKTrZxR"  
%*********************************************************************** \iMyo  
%     Update magnetic fields Sb^ b)q"  
%*********************************************************************** |*Ot/TvG  
L2Ux9_S  
hx(2:ie,1:je,1:ke)=hx(2:ie,1:je,1:ke)+... xaXV^ZM3  
                   db*(ey(2:ie,1:je,2:kb)-ey(2:ie,1:je,1:ke)+... MWq$AK]  
                       ez(2:ie,1:je,1:ke)-ez(2:ie,2:jb,1:ke)); @)sc6 *lnW  
                 PW_`qP:  
hy(1:ie,2:je,1:ke)=hy(1:ie,2:je,1:ke)+... 4`mO+.za1  
                   db*(ex(1:ie,2:je,1:ke)-ex(1:ie,2:je,2:kb)+... Sa]mm/G  
                       ez(2:ib,2:je,1:ke)-ez(1:ie,2:je,1:ke)); I$N7pobh  
                 R'EW7}&  
hz(1:ie,1:je,2:ke)=hz(1:ie,1:je,2:ke)+... ) Ypz!  
                   db*(ex(1:ie,2:jb,2:ke)-ex(1:ie,1:je,2:ke)+... KR3-Hb4  
                       ey(1:ie,1:je,2:ke)-ey(2:ib,1:je,2:ke)); E_[ONm=,  
   j%M @#  
smooth(ez(:,45,30))   ;   J5T=!wF (  
figure(1); B-Y+F  
subplot(1,2,1) 0&`}EXe<f  
surf(ez(:,:,30)); NyJnO
w(  
colorbar; 9YSVK\2$  
pause(0.1); 5fHYc0  
disp('ez()='); ZBj6KqfST%  
disp(ex(45,45,30)); IYNM
U\s  
AA(T)=ex(45,45,30); !cCg/  
-,>:DUN2  
+ x_wYv  
subplot(1,2,2) >=!AL,:  
imagesc(ez(:,:,30)); ci7~KewJ*  
colorbar; Qp&?L"U)2  
`?SGXXC  
v+tO$QZ`  
end 96W4c]NT  
*4#on>  
T=(1:1:nsteps); N^@ \tg=  
figure(2); y<r44a_!  
plot(T,AA .. K<,Y^3]6?  
*IlQ5+3I  

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

2.5hz的电磁波

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

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