大家好,欢迎来到IT知识分享网。
一维FDTD程序如下
https://wenku.baidu.com/view/699d2b5b4b35eefdc8d333cf.html
为了减小误差,采用二阶中心差分,所以看起来是十字状,呈时空交错的
function FDTDonedimensionpipei(L,d,T) %version1.0 终端匹配 %FDTDonedimensionpipei(6,0.18,0.5e-9) t0=3*T; c=3e8; u=4*pi*1e-7; e=8.8541878e-12; dz=T*c/10; Nz=fix(L/dz); dt=dz/2/c; Ex=zeros(1,Nz+1); B=zeros(1,Nz+1); Hy=zeros(1,Nz); Nt=2*Nz; for n=0:Nt t=n*dt; F=exp(-(t-t0).^2./T^2); Ex(1)=F; for k=1:Nz Hy(k)=Hy(k)+dt./u.*(Ex(k)-Ex(k+1))./dz; end for k=1:Nz-1 Ex(k+1)=Ex(k+1)+dt./e.*(Hy(k)-Hy(k+1))./dz; end Ex(1)=B(2)+(c*dt-dz)./(c*dt+dz).*(Ex(2)-B(1)); Ex(Nz+1)=B(Nz)+(c*dt-dz)./(c*dt+dz).*(Ex(Nz)-B(Nz+1)); Vref1=d.*Ex(Nz-300); Vref2=d.*Ex(Nz-100); plot(t,Vref1,'s'); hold on; plot(t,Vref2,'rx'); hold on; B=Ex; end
免责声明:本站所有文章内容,图片,视频等均是来源于用户投稿和互联网及文摘转载整编而成,不代表本站观点,不承担相关法律责任。其著作权各归其原作者或其出版社所有。如发现本站有涉嫌抄袭侵权/违法违规的内容,侵犯到您的权益,请在线联系站长,一经查实,本站将立刻删除。 本文来自网络,若有侵权,请联系删除,如若转载,请注明出处:https://yundeesoft.com/34392.html