Nz=400;
Nt=200;
N_source=Nz/2;
spread=12;
E=zeros(Nz,1);
H=E;
for t=1:Nt
for m=2:Nz;
E(m)=E(m)+0.5*(H(m-1)-H(m));
end
pulse=exp(0.5*(t/spread)^2);
E(N_source)=pulse+E(N_source);
for m=1:(Nz-1)
H(m)=H(m)+0.5*(E(m)-E(m+1));
end
%display
figure(1)
plot(E,t);
xlabel('Ex')
ylabel('t')
axis([0 Nz -1 1])
pause(0.001);
end
(敬请帮忙指正,我的QQ328248640,希望不吝赐教) |