2392阅读
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 %%%%%%%%%% T:0#se
ex(1)=ex_low_m2; wvz_)bN~A
ex_low_m2=ex_low_m1; + 0 |d2_]E
ex_low_m1=ex(2); C^fUhLVSZ^
?TMrnR/d
hy(200)=ex_high_m2; u&Cu"-%=M
ex_high_m2=ex_high_m1; #xNXCBl]O
ex_high_m1=hy(199); q~6((pWi|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

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 ..

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