一维FDTD问题请教

我用Matlab编了一个程序，来仿真这样一个过程：一维条件，中部有个高斯脉冲激励源，其幅度最大值为1。程序运行后，有两个对称的高斯脉冲向两边传播。但我有两个问题： m<LzB_G\  
一、生成的两个脉冲幅度的最大值明显比1小，但我觉得应该和1差不多，不知道正确的应该是怎样。 =<uz'\Ytv%  
二、刚开始的一段时间内，源点位置有很大的负值分量，当两个脉冲分开后 .. kT=|tQ@  
' g!_Flk  

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

1)你直接加就源就是1 了，比如Ex(center) = f(n)。 (~o+pp!  
估计你写的是Ex(center) = Ex(center) + f(n) +T;qvx6  
}Ec"&  
2)我觉得你看错了，不应该是负的很大值，而是exp(-number)，是这个很大吧，无限接近于零。

（1）可以看看边界设置，这是衰减吧

把程序贴上来看看

是啊，把程序贴出来，大家帮你分析一下

不好意思，我当时比较急，后来就忘了来这里看了，先谢谢给位的关心，并顺祝元旦快乐！程序我后来稍微修改了一下。 9H<:\-:  
q&+GpR  
%初始化光速、介电常数、磁导率等常数 HTC7fS  
c=3e8; 4DNZ y2`  
ep0=8.854e-12; Og["X
0j  
ep=9*ep0; F"q3p4-<>  
u=pi*4e-7; lIf Our  
T=2e-11; a'|0e]  
t0=4*T; dZm{?\^_  
%初始化时间步长、空间步长和总的执行步数 =\3*;59\  
dt=3e-13; QMIQy  
dz=2e-4; 3|A"CU/z@  
tmax=2e-9; mHP1
.Z`  
step=round(tmax/dt); &I70veNY  
%设置介质分界面的坐标 XDn$=`2  
m1=[0.45 0.45]; ?TXFOr]g]2  
m2=[0.5 0.5]; GV9"8MZ6  
m=[-3 3]; c/s'&gG33z  
6x{B  
N=round(0.6/dz); 2`z+_DA  
n1=round(0.45/dz+1); f7`y*9^  
n2=round(0.5/dz+1); 2XE4w# [j  
ns=round(0.3/dz+1); w9 NUm  
fun=@(t) exp(-((t-t0)./T).^2); HdGy$m`  
Ex=zeros(1,N+1); ev; &$Hc  
Hy=zeros(1,N); 9(C Ke,  
t=0:dt:dt*step; NzEuiI}  
Z=0:dz:0.6; l6O2B/2
j  
Et=fun(t); 27MgwX NQ  
ca=dt/(dz*u); `1P &  
cb=zeros(1,N+1); 3 3V/<v  
cb(:)=dt/(dz*ep0); :( `Q4D~l  
cb([n1 n2])=dt/(dz*(ep0+ep)/2); U{ Y)\hR-  
cb(n1+1:n2-1)=dt/(dz*ep); {\%x{  
ct=(c*dt-dz)/(c*dt+dz); Z/0M9 Q%  
%电场和磁场相互迭代 jj^CW"IB  
for i=1:step (<f[$ |%  
    Ex(ns)=Ex(ns)+Et(i); (YR1ML3N  
    Hy(1:N)=Hy(1:N)-ca*(Ex(2:N+1)-Ex(1:N)); F2u{Wzr_@  
    temp1=Ex(2); -0W;b"]+A  
    temp2=Ex(N); +n0y/
0Au  
    Ex(2:N)=Ex(2:N)-cb(2:N).*(Hy(2:N)-Hy(1:N-1)); SZgH0W("L  
    Ex(1)=temp1+ct*(Ex(2)-Ex(1)); |h3YL!  
    Ex(N+1)=temp2+ct*(Ex(N)-Ex(N+1)); | o?@Eh  
     ^Ab|\5^3  
    %画图 V7&L+]!  
    plot(Z,Ex); ,0~'#x>  
    line(m1,m); 0$UE|yDs>  
    line(m2,m); 3 -5^$-7_  
    ylabel('Ex');xlabel('z'); 67#;.}4a  
    axis([0 0.6 -0.8 1.2]); 55#H A?cR  
    pause(0.0001); 8T%z{A1T  
end ci$o~b6V  
q H+~rj  

这个就如上面说的