打印

大虾指教一下关于FDTD一维MATLAB编程问题

[复制链接]
1126|0
手机看帖
扫描二维码
随时随地手机跟帖
跳转到指定楼层
楼主
immortaldeity|  楼主 | 2012-4-30 20:35 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
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,希望不吝赐教)

相关帖子

发新帖 我要提问
您需要登录后才可以回帖 登录 | 注册

本版积分规则

0

主题

2

帖子

1

粉丝