一维FDTD问题请教

我用Matlab编了一个程序，来仿真这样一个过程：一维条件，中部有个高斯脉冲激励源，其幅度最大值为1。程序运行后，有两个对称的高斯脉冲向两边传播。但我有两个问题： yN<fmi};c  
一、生成的两个脉冲幅度的最大值明显比1小，但我觉得应该和1差不多，不知道正确的应该是怎样。 }a1Sfl@`3  
二、刚开始的一段时间内，源点位置有很大的负值分量，当两个脉冲分开后 .. f}Tr$r  
Z:J.FI@  

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

1)你直接加就源就是1 了，比如Ex(center) = f(n)。 x1BOW  
估计你写的是Ex(center) = Ex(center) + f(n) ^5^}MB%  
Y <Znv%M  
2)我觉得你看错了，不应该是负的很大值，而是exp(-number)，是这个很大吧，无限接近于零。

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

把程序贴上来看看

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

不好意思，我当时比较急，后来就忘了来这里看了，先谢谢给位的关心，并顺祝元旦快乐！程序我后来稍微修改了一下。 .NJ|p=fy  
-Ep!- a  
%初始化光速、介电常数、磁导率等常数 )MZC>
:  
c=3e8; !}y1C
A  
ep0=8.854e-12; n`(~OO  
ep=9*ep0; m9$a"$c  
u=pi*4e-7; 8eluO ?p  
T=2e-11; x4m 5JDC  
t0=4*T; t9FDU  
%初始化时间步长、空间步长和总的执行步数 igDG}q3jG  
dt=3e-13; 
q-nER<  
dz=2e-4; q"akrI38  
tmax=2e-9; i9rS6<
V'  
step=round(tmax/dt); B**Nn!}0  
%设置介质分界面的坐标 !9;)N,  
m1=[0.45 0.45]; DG(%-w8p"  
m2=[0.5 0.5]; @x1%)1  
m=[-3 3]; [4*1}}gW%5  
!KcWH9  
N=round(0.6/dz); Qzqc .T  
n1=round(0.45/dz+1); V1B(|P  
n2=round(0.5/dz+1); s)zJT  
ns=round(0.3/dz+1); pMR,#[U<  
fun=@(t) exp(-((t-t0)./T).^2); =E(#YCx  
Ex=zeros(1,N+1); ~6`iY@)  
Hy=zeros(1,N); h >-'-Hx+  
t=0:dt:dt*step;
5RF4]$zT  
Z=0:dz:0.6; ACg;CTBb  
Et=fun(t); ExVDkt0  
ca=dt/(dz*u); AF}gSNX  
cb=zeros(1,N+1); %@?A_jS  
cb(:)=dt/(dz*ep0); VO=!8Yx[  
cb([n1 n2])=dt/(dz*(ep0+ep)/2); m:uPEpcU  
cb(n1+1:n2-1)=dt/(dz*ep); ?cKZ_c  
ct=(c*dt-dz)/(c*dt+dz); 3`*Kav>"  
%电场和磁场相互迭代 *6Q|}b[qcD  
for i=1:step =su]w2,Iy  
    Ex(ns)=Ex(ns)+Et(i); c 6"Ib)  
    Hy(1:N)=Hy(1:N)-ca*(Ex(2:N+1)-Ex(1:N)); /c`^iPb  
    temp1=Ex(2); \'Ssn(s  
    temp2=Ex(N); wN97_Y=`n  
    Ex(2:N)=Ex(2:N)-cb(2:N).*(Hy(2:N)-Hy(1:N-1)); 5\bJR0I@  
    Ex(1)=temp1+ct*(Ex(2)-Ex(1)); v4rW2F:X  
    Ex(N+1)=temp2+ct*(Ex(N)-Ex(N+1)); 6)qp*P$L  
     ba3*]01Yb  
    %画图 jV
O{$j  
    plot(Z,Ex); DkD
oA;m  
    line(m1,m); eQUe >*  
    line(m2,m); tIc 7:th  
    ylabel('Ex');xlabel('z'); w##^}nHOR  
    axis([0 0.6 -0.8 1.2]); Y-)xTn  
    pause(0.0001); :@eHX&  
end O^Q,-=tA\  
, sjh^-;  

这个就如上面说的