[求助]程序求助啊

第一次做FDTD程序不对，求助 b3H;Ea?^^<  
{Bx\Z0+'&  
clc E J1:N*BA  
clear UpbzH(?#  
(WC<XKf  
:L,]<n  
c=2.99792458e8;                TocqoYX{{  
muz=4.0*pi*1.0e-7;             4(NI-|q0  
epsz=1.0/(c*c*muz);         w}iflAnjq  
freq=1.0e+9;                   ocvBKsfhE`  
lamdaz=c/freq;                 A<&9  
omiga=2.0*pi*freq;             8wOr`ho B  
dx=lamdaz/20.0;                `?:'_Ki  
dt=dx/c;                       <7XdT  
N=200;                        pR$c<p  
nmax=round(12.0e-9/dt);                      zI(Pti  
eUl[gHP  
Z(s} #-  
Q]\xO/  
LxbVRw  
eps=1.0; u[% #/  
sig=5.0e-3; shD$,! k  
sigm=0; Pdv&X*
KA  
E?-K_p  
qQb8K+t  
V
:lDR20*\  
}v:h EMO  
CA=(1-(sig*dt/(2*epsz*eps)))/(1+(sig*dt/(2*epsz*eps))); x-"7{@lz  
CB=(dt/(eps*epsz*dx))*(1/1+(1+(sig*dt)/(2*epsz*eps))); FUQT,7CA  
C]k\GlhB  
\%K6T)9  
L:31t
oGK  
8w\&QX  
Ex(1:N+1)=0.0; ;J4_8N-  
Hy(1:N)=0.0; ,{%[/#~6  
?lTQjw{  
hX^XtIC=  
x=linspace(dx,N*dx,N); up`6IWlLE  
subplot(2,1,1),plot(x,Ex(1:N)/2,'r'),axis([-3 3 -1 1]); 3cnsJV]  
ylabel('Ex'); TnC'<zm9!  
subplot(2,1,2),plot(x,Hy,'b'),axis([0 3 -3.0e-3 3.0e-3]); hAY_dM  
xlabel('x (meters)');ylabel('Hy'); N7NK1<vw2  
rect=get(gcf,'Position'); vK$W)(Z  
rect(1:2)=[0 0]; d"V^^I)yx&  
M=moviein(nmax/2,gcf,rect); u`ZnxD>  
 z\\MLyS  
%T&kK2d;  
for n=1:nmax H;v*/~zl  
G#csN&|,  
Dd*C?6  
Ex(1)=sin(omiga*dt*n); ].$N@tC  
r=Ex(N); 'rSM6j  
Ex(2:N)=CA*Ex(2:N)+CB*(Hy(2:N)-Hy(1:N-1)); C 9DRVkjj  
!#'*@a  
5!w
a\)wY  
r=Ex(N+1); rtPQ:CaA)?  
ZqSczS7uf  
Alrk3I3{  
Hy(1:N)=DA*Hy(1:N)+CB*(Ex(2:N+1)-Ex(1:N)); P:
Bg
()  
VH8,!#Q;  
RxAZ<8T_  
if mod(n,2)==0; 4kIy4x'*  
time=num2str(round(n*dt/1.0e-9)); Tfj%Sb,zM  
Cjw|.c`  
DRal{?CH  
subplot(2,1,1),plot(x,Ex(1:N)/2,'r'),axis([0 3 -1 1]); "zv

?qS  
title(['time = ',time,' ns']); H>2)R7h  
ylabel('Ez'); <"6}C)G  
c!b4Y4eJ  
Uh{|@D  
subplot(2,1,2),plot(x,Hy,'b'),axis([0 3 -3.0e-3 3.0 .. kid@*.I  
u1rT:\G1  

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

哪有，da,db的参数？

不好意思复制是删掉了，问题貌似不在这。DA=1 DB=CB

db=cb不对吧。你还是检查检查吧。一般系数对了没有什么问题的。