2390阅读
4回复

[RFEDA原创]关于Sullivan Dennis中的一个例子

离线langren2007

UID ：6360

注册：2007-12-03
登录：2020-04-07
发帖：194
等级：仿真二级

0楼发表于: 2008-03-27 11:25:06

Sullivan Dennis那本书上第一章的第二个例子，一维的完全吸收边界。但为什么右行的波到达边界后发生了反射？
这是做到二维后才发现的，反射相当厉害。
希望大家给看看，不知道有没有同样的问题？

clear
clc
ke=200;
ex=rand(1,200);
hy=rand(1,200);
ex_low_m1=0;
ex_low_m2=0;
ex_high_m1=0;
ex_high_m2=0;
syms n k kc ke nsteps T t0 spread pulse;
%initialize
for k=1:200
ex(k)=0;
hy(k)=0;
end
kc=100;
t0=40;
spread=12;
T=0;
nsteps=1;

while nsteps>0
nsteps=input('nsteps-->');

n=0;
for n=1:nsteps
T=T+1;
for k=2:200
ex(k)=ex(k)+0.5*(hy(k-1)-hy(k));
end
pulse=exp(-0.5*(power((t0-T)/spread,2)));
ex(kc)=pulse;
%%%%%%% absorbing boundary %%%%%%%%%% ^?-:'<4q$
ex(1)=ex_low_m2; Qk0R a_
ex_low_m2=ex_low_m1; Sj-n;F|=X
ex_low_m1=ex(2); t=(!\:[D
"LHcB]^<
hy(200)=ex_high_m2; .S!mf
ex_high_m2=ex_high_m1; HD@$t)mn
ex_high_m1=hy(199); dV*]f$wQ
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

for k=1:199
hy(k)=hy(k)+0.5*(ex(k)-ex(k+1));
end
end

%end of main FDTD loop

fprintf('T=%5.0f\n',T);

%end
x=1:200;
subplot(2,1,1);
plot(x,ex(x));
axis([0,200,-2,2]);
grid on
title('FDTD cells ..

未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册