[求助]程序求助啊

第一次做FDTD程序不对，求助 ]
x5(bnWx  
(mOqv9pn  
clc }R?v"6aBS  
clear 3> n2  
9Y!0>&o  
u(TgWp5WF  
c=2.99792458e8;                A
>ug'.  
muz=4.0*pi*1.0e-7;             9!Q $GE?vl  
epsz=1.0/(c*c*muz);         GA%"w=M\  
freq=1.0e+9;                   {d*qlztO  
lamdaz=c/freq;                 ~{3o(gzl  
omiga=2.0*pi*freq;             &+|jJ{93z  
dx=lamdaz/20.0;                >F LdI  
dt=dx/c;                       BIHHRCe:@n  
N=200;                       
rTm>8et
 
nmax=round(12.0e-9/dt);                      5Tcl<Y6l  
3.GdKP.%  
f<$K.i  
Q
 KDb  
"`8H:y  
eps=1.0; -=v/p*v0o  
sig=5.0e-3; 0&SrKn  
sigm=0; CguU+8]  
Ot4;,UZ  
Sl+jduc  
2s]]!{Z#  
:"y7Weh  
CA=(1-(sig*dt/(2*epsz*eps)))/(1+(sig*dt/(2*epsz*eps))); 2:6Y83  
CB=(dt/(eps*epsz*dx))*(1/1+(1+(sig*dt)/(2*epsz*eps))); +lk\oj$S+  
->RF`SQu  
&zg$H,@Qp  
8)ebXc  
q:D0$YY0  
Ex(1:N+1)=0.0; nf+"vr}1  
Hy(1:N)=0.0; /4f 5s#hR  
b'Piymx  
y:^>
(l#;  
x=linspace(dx,N*dx,N); lN=m$J  
subplot(2,1,1),plot(x,Ex(1:N)/2,'r'),axis([-3 3 -1 1]); GakmROZ@9  
ylabel('Ex'); "\R@lUx.Y  
subplot(2,1,2),plot(x,Hy,'b'),axis([0 3 -3.0e-3 3.0e-3]); f6aT[Nw<  
xlabel('x (meters)');ylabel('Hy'); H*:r>
Lm=  
rect=get(gcf,'Position'); <(6-9(zHa  
rect(1:2)=[0 0]; zwniS6R1  
M=moviein(nmax/2,gcf,rect); L`VQ{|&3V  
,*O{jc`(  
)ZuQ;p  
for n=1:nmax hD9'`SQ  
,sRrV $,"  
>8Yrmq  
Ex(1)=sin(omiga*dt*n); JE8p5WaR  
r=Ex(N); pyb}ha  
Ex(2:N)=CA*Ex(2:N)+CB*(Hy(2:N)-Hy(1:N-1)); BSq;R
G(  
3gfV0C\  
eA~_)-Z-  
r=Ex(N+1); ^r?sgJ  
x:&L?eOT  
,~=z_G`R  
Hy(1:N)=DA*Hy(1:N)+CB*(Ex(2:N+1)-Ex(1:N)); ,co9f.(w  
c?3F9w#  
y>pq*i  
if mod(n,2)==0; *9U4^lJjn  
time=num2str(round(n*dt/1.0e-9)); H(*=9  
IZ(CRKCGBl  
[gg7Z|Hu  
subplot(2,1,1),plot(x,Ex(1:N)/2,'r'),axis([0 3 -1 1]); b`={s  
title(['time = ',time,' ns']); | xErA  
ylabel('Ez'); =+sIX3
  
#mz,HK0|aC  
uN?O*h/(  
subplot(2,1,2),plot(x,Hy,'b'),axis([0 3 -3.0e-3 3.0 .. LBiv]3  
U}C#:Xi>$  

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

哪有，da,db的参数？

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

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