[求助]能帮忙看看我这easy的程序么？为什么出来的图总是发散的？

200步是达到10^57了  散的厉害 L)y}  
6X7r=w  
3P^eD:) w  
%定义几个常量 MR#j
I  
c=3*10^8; rZKv:x}{6  
e0=8.854*10^(-12); !`=r('l  
u0=4*pi*10^(-7); <saS2
.4  
h'MX{Wm.  
deta_x=0.001;%定义空间网格长度，为1/20波长 Hu1w/PLq  
deta_y=0.001; [9'5+RXw3  
deta_z=0.001; ;s9!ra:3  
x_max=50;%定义总的网格数 ;yBq'_e3  
y_max=50; k3sP,opacX  
z_max=50; &EE6<-B-  
%定义时间步长 %
nFZA)B[  
deta_t=1/c/sqrt((1/deta_x)^2+(1/deta_y)^2+(1/deta_z)^2); 0=t_a]+  
%定义迭代的次数 TcauCL  
step=100; xDRK^nmC  
%场量置零 IR5 S-vO  
hx1=zeros(x_max+1,y_max,z_max); UF,T  
hx2=zeros(x_max+1,y_max,z_max); 9oKRu6]D-  
hy1=zeros(x_max,y_max+1,z_max); !xj>~7  
hy2=zeros(x_max,y_max+1,z_max); AJCWp4,  
hz1=zeros(x_max,y_max,z_max+1); i:qc2#O:J  
hz2=zeros(x_max,y_max,z_max+1); A>R ^iu  
ez1=zeros(x_max+1,y_max+1,z_max); PM[6U#  
ez2=zeros(x_max+1,y_max+1,z_max); MCz+l0  
ex1=zeros(x_max,y_max+1,z_max+1); _YmYy\g  
ex2=zeros(x_max,y_max+1,z_max+1); v:JFUn}  
ey1=zeros(x_max+1,y_max,z_max+1); _~HGMC)  
ey2=zeros(x_max+1,y_max,z_max+1); EpOVrk  
%定义电流源的位置 YRW<n9=3  
source_x=25; WuSRA<{P  
source_y=25; [EB2o.EsO  
source_z=25; ]r'b(R; S  
%定义观察点的位置 ?9mWMf%t  
field_x=23; S4~^HvMG[Y  
field_y=23; _q3SR[k+`  
field_z=23; Wq&TbWR  
%定义激励源 '9$xOrv  
%时谐源 ]L}<Y9)t  
frequency=15*10^9;%波长为0.02 >Q[]i4*A  
omiga=2*pi*frequency;
0t/S_Q  
%高斯源 o"7,CQye  
t0=0.0667*10^(-9); Vxo3RwmR  
tt=0.0667*10^(-9); 'G@Npp)&^  
%定义保存结果的变量 Da?0B9'  
result=zeros(step,1,1); }eEF/o  
%迭代开始 {7Dc(gNS  
for t=1:1:step +&(sZFW5o  
    t }..}]J;To  
   % current(t)=sin(omiga*t*deta_t); Qy0bp;V/  
  current(t)=exp(-4*pi*(t*deta_t-t0)^2/tt^2); sBwkHsDD  
 %ex的计算； d!Y,i!l!  
for i=1:1:x_max 5C#&vYnq  
    for j=2:1:y_max /5U?4l(6[f  
        for k=2:1:z_max r}"Ty  
            if i==source_x&j==source_y&k==source_z%判断是不是源点 OFQsfW3O  
               ex2(i,j,k)=ex1(i,j,k)+deta_t/e0*((hz2(i,j,k)-hz2(i,j-1,k))/deta_y-(hy2(i,j,k)-hy2(i,j,k-1))/deta_z)-deta_t/e0*current(t); w2/%e$D!9  
            else l]oGhM;  
               ex2(i,j,k)=ex1(i,j,k)+deta_t/e0*((hz2(i,j,k)-hz2(i,j-1,k))/deta_y-(hy2(i,j,k)-hy2(i,j,k-1))/deta_z); n;v8Vc'  
            end _.?$~;7  
        end x*Z"
~'DI  
    end *5*d8;@>  
end hCW8(Zt  
  _Xsn1  
%ey的计算； 6l:CDPhR  
for i=2:1:x_max `.L8<-]W  
    for j=1:1:y_max J[VQ6fD%  
        for k=2:1:z_max X+P3a/T  
            if i==source_x&j==source_y&k==source_z%判断是不是源点 Z|'tw^0e5  
               ey2(i,j,k)=ey1(i,j,k)+deta_t/e0*((hx2(i,j,k)-hx2(i,j,k-1))/deta_z-(hz2(i,j,k)-hz2(i-1,j,k))/deta_x)-deta_t/e0*current(t); eHPGzNXb  
            else i+21tG$  
               ey2(i,j,k)=ey1(i,j,k)+deta_t/e0*((hx2(i,j,k)-hx2(i,j,k-1))/deta_z-(hz2(i,j,k)-hz2(i-1,j,k))/deta_x); ?y]3kU  
            end t
op3o{4  
        end wy:.  
    end MbYa6jrF  
end "=P@x|I  
%ez的计算； ${6 ;]ye  
for i=2:1:x_max Zt[ PkBi  
    for j=2:1:y_max k3S**&i!CR  
        for k=1:1:z_max +SUQRDF@i  
            if i==source_x&j==source_y&k==source_z%判断是不是源点 \-V  
                ez2(i,j,k)=ez1(i,j,k)+deta_t/e0*((hy2(i,j,k)-hy2(i-1,j,k))/deta_x-(hx2(i,j,k)-hx2(i,j-1,k))/deta_y)-deta_t/e0*current(t); u'Mq^8  
            else m"H9C-Y  
                ez2(i,j,k)=ez1(i,j,k)+deta_t/e0*((hy2(i,j,k)-hy2(i-1,j,k))/deta_x-(hx2(i,j,k)-hx2(i,j-1,k))/deta_y); AD*+?%hj  
            end XImb"7|  
        end .!
nFy`  
    end Q?W]g%:)  
end > {fX;l  
%hx的计算； xCoQ>.4p  
for i=1:1:x_max+1 n+Fl|4  
    for j=1:1:y_max owmA]f  
        for k=1:1:z_max #G)ZhgB^  
            hx2(i,j,k)=hx1(i,j,k)-deta_t/u0*((ez1(i,j+1,k)-ez1(i,j,k))/deta_y-(ey1(i,j,k+1)-ey1(i,j,k))/deta_z); VNXB7#ry  
        end %S$P+B?  
    end uh
:  
end MJ}VNv|S  
%hy的计算； rIS \#j  
for i=1:1:x_max NBUM* Z  
    for j=1:1:y_max+1 %z.u % %  
        for k=1:1:z_max 89'nbg  
            hy2(i,j,k)=hy1(i,j,k)-deta_t/u0*((ex1(i,j,k+1)-ex1(i,j,k))/deta_z-(ez1(i+1,j,k)-ez1(i,j,k))/deta_x); 9)l[$X  
        end ~l{CUQU  
    end ;<ed1%Le,  
end 47ir QK*  
%hz的计算； PS(LD4mD  
for i=1:1:x_max +c\fDVv  
    for j=1:1:y_max O23f\pm&  
        for k=1:1:z_max+1 ec"L*l"  
            hz2(i,j,k)=hz1(i,j,k)-deta_t/u0*((ey1(i+1,j,k)-ey1(i,j,k))/deta_x-(ex1(i,j+1,k)-ey1(i,j,k))/deta_y) .. zzq/%jki  
``>WFLWTn  

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