首页 > 代码库 > FDFD一维
FDFD一维
一维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
FDFD一维
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。