[求助]程序求助啊

第一次做FDTD程序不对，求助 1/ a,7Hl  
AQ&;y&+QR  
clc Ugn"w E  
clear 9
}=Fdt  
$_y"P  
LakP'P6`E  
c=2.99792458e8;                zK&`&("4C  
muz=4.0*pi*1.0e-7;             c~V\,lcI  
epsz=1.0/(c*c*muz);         U*Q$:%72vO  
freq=1.0e+9;                   /#a$4 }2L  
lamdaz=c/freq;                 ^6=nL<L  
omiga=2.0*pi*freq;             y~\z_') <>  
dx=lamdaz/20.0;                C-Mop,w  
dt=dx/c;                       W<2-Q,>Y  
N=200;                        <K43f#%  
nmax=round(12.0e-9/dt);                      ;tHF$1!J  
!@Ox%vK  
<(v!Xj^yO  
8WvT0q>]  
r`=d4dK-  
eps=1.0; w/&#UsEIr  
sig=5.0e-3; /N$T[  
sigm=0; J-U}iU|  
*]LM2
J  
FY'f{gD^  
B>R6j}rh'k  
7-^d4P+|g  
CA=(1-(sig*dt/(2*epsz*eps)))/(1+(sig*dt/(2*epsz*eps))); 4x:fOhtP  
CB=(dt/(eps*epsz*dx))*(1/1+(1+(sig*dt)/(2*epsz*eps))); YRK4l\_`  
vXc<#X9  
<
X: 9y  
8WAg{lV
s  
t((0]j^
 
Ex(1:N+1)=0.0; )3;S;b  
Hy(1:N)=0.0; y,aASy!Q  
milU,!7J  
A9"ho}<  
x=linspace(dx,N*dx,N); js{ RaR=  
subplot(2,1,1),plot(x,Ex(1:N)/2,'r'),axis([-3 3 -1 1]); O_E[FE:+  
ylabel('Ex'); wDsEx!\#  
subplot(2,1,2),plot(x,Hy,'b'),axis([0 3 -3.0e-3 3.0e-3]); # RtrHm  
xlabel('x (meters)');ylabel('Hy'); fE(rDQI  
rect=get(gcf,'Position'); DV.m({?  
rect(1:2)=[0 0]; yEH30zSt  
M=moviein(nmax/2,gcf,rect); n7S[ F3  
EprgLZ1B  
<+6)E@Y  
for n=1:nmax yJ^
}uw  
0j*8|{|  
@=$;^}JS|  
Ex(1)=sin(omiga*dt*n); P63f0F-G  
r=Ex(N); Eq|_>f@@8  
Ex(2:N)=CA*Ex(2:N)+CB*(Hy(2:N)-Hy(1:N-1)); noacnQ_I$  
7&X^y+bMe6  
3+)i23[4=\  
r=Ex(N+1); $zJ.4NA  
t({:TQ  
NK#f Gz*,(  
Hy(1:N)=DA*Hy(1:N)+CB*(Ex(2:N+1)-Ex(1:N)); Uu G
;z5  
r!;NH3 *  
mF "ctxE  
if mod(n,2)==0; +;vfn>^!b  
time=num2str(round(n*dt/1.0e-9)); km *$;Nli  
-G{}8GM  
y'(;!5w  
subplot(2,1,1),plot(x,Ex(1:N)/2,'r'),axis([0 3 -1 1]); WKN\*N<  
title(['time = ',time,' ns']); MQhL>oQ  
ylabel('Ez'); ,ujoGSx}  
P?>p+dM  
1:8ZS  
subplot(2,1,2),plot(x,Hy,'b'),axis([0 3 -3.0e-3 3.0 .. Gv<K#@9T  
uoF9&j5E@Z  

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

哪有，da,db的参数？

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

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