FDFD一维_IT分享知识网

FDFD一维

老牧童 • 2024-01-06 07:45 • 未分类 • 阅读 110

大家好，欢迎来到IT知识分享网。

一维FDTD程序如下

https://wenku.baidu.com/view/699d2b5b4b35eefdc8d333cf.html

FDFD一维

为了减小误差，采用二阶中心差分，所以看起来是十字状，呈时空交错的

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

FDFD一维

相关推荐

发表回复