982阅读
6回复

[求助]关于一维 FDTD Mur 边界的问题（求助）

资源共享呗

UID ：71951

注册：2011-01-22
登录：2014-09-27
发帖：389
等级：仿真三级

0楼发表于: 2011-09-10 12:20:10

      有人发过这样的帖子·我刚学FDTD ，也遇到了这样的问题。就是“最近在看葛德彪的《电磁波时域有限差分方法》一书，在看到Mur边界时，自己做的一维Mur边界的吸收效果不好。用Gauss脉冲激励（最大值为1），在边界边会有0.2的反射波不能被吸收。
程序调了很久，就是不能消除反射波，不知道有没有人做过这个简单的一维FDTD程序，有没有出现这样的情况？”

           在论坛找了好久也没找到解决问题的办法和相关的代码希望知道的朋友能告诉一下。先谢谢了··
          并把我的错误程序放上，希望能给指出错误···
[backcolor=#c ..

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

附件:

fdtd_1.rar (1 K) 下载次数:16

离线tingo

临渊羡鱼，不如退而结网

UID ：28673

注册：2009-03-29
登录：2019-09-21
发帖：126
等级：仿真三级

1楼发表于: 2011-09-10 21:30:31

你还是把程序直接贴出来吧，我RF币不够，连你的附件都下不了，就没法帮你看程序了！

仿真新人！

离线da376805618

资源共享呗

UID ：71951

注册：2011-01-22
登录：2014-09-27
发帖：389
等级：仿真三级

2楼发表于: 2011-09-11 15:27:40

回 1楼(tingo) 的帖子

程序如下：

max_time = 300;
max_space = 200;

eta = .5;
delta_x = .01;
delta_t = delta_x / (2*3e8);
c = 3e8;

E = zeros(max_space,1);
H = E ;
Ez_n=zeros(1,4);

spread = 12;
center_problem_space = max_space / 2;
t0 = 40 ;

%定义边界处的系数
K_M1 = (c*delta_t-delta_x)/(c*delta_t+delta_x);

for n  = 1:max_time
     %左边界
     Ez_n_1(1) = Ez_n(1);
     Ez_n(1) = E(1);
     Ez_n_1(2) = Ez_n(2);
     Ez_n(2) = E(2);
     %右边界
     Ez_n_1(3) = Ez_n(3);
     Ez_n(3) = E(max_space-1);
     Ez_n_1(4) = Ez_n(4);
     Ez_n(4) = E(max_space);
    for k = 2:max_space
        E(k) = E(k) + eta*(H(k-1) - H(k));
    end

     %Mur ABC
     %at x=1;
     E(1) = Ez_n_1(2) + K_M1*(E(2) - Ez_n_1(1));
     %at x=max_space
     E(max_space) = Ez_n_1(3) + K_M1*(E(max_space-1) - Ez_n_1(4)) ;

    pulse = exp(-.5*((t0 - n)/spread)^2);
    E(center_problem_space) = pulse ;


    for j = 1:max_space-1
        H(j) = H(j) + eta*(E(j) - E(j + 1));
    end
    figure(1)
    plot(E)
    axis([1 max_space -1.1 1.1])
    pause(.05)
end
figure(2)
subplot(2,1,1)
plot(E,'b')
ylabel('Ex')
axis([1 max_space -1.1 1.1])
subplot(2,1,2)
plot(H,'r')
ylabel('Hy')
axis([1 max_space -1.1 1.1])

资源共享

离线da376805618

资源共享呗

UID ：71951

注册：2011-01-22
登录：2014-09-27
发帖：389
等级：仿真三级

3楼发表于: 2011-09-11 15:33:52

回 1楼(tingo) 的帖子

顺便再问一个问题··我看别人的程序中对于吸收边界的处理如下：
  E(1) = E_low_m2;
    E_low_m2 = E_low_m1;
    E_low_m1 = E(2);
    E(max_space) = E_high_m2;
    E_high_m2 = E_high_m1;
    E_high_m1 = E(max_space - 1);

它依据的原理是什么呢？还有它是《电磁波时域有限差分方法》一书中所说的行波延迟法么？我看他得书上这样写的:
  当c*delta_t  =  delta_x 时有 E n+1 (k) =  E n(k-1)    (n+1和n 是E 的上角标)
  还请大侠指点下··谢谢了  呵呵·

资源共享

离线tingo

临渊羡鱼，不如退而结网

UID ：28673

注册：2009-03-29
登录：2019-09-21
发帖：126
等级：仿真三级

4楼发表于: 2011-09-12 08:47:45

呵呵，你的问题出在蛙跳格式出错了！

仿真新人！

离线tingo

临渊羡鱼，不如退而结网

UID ：28673

注册：2009-03-29
登录：2019-09-21
发帖：126
等级：仿真三级

5楼发表于: 2011-09-12 08:55:23

回 2楼(da376805618) 的帖子

不过，如果是一维貌似差别不大，我再看看你的程序！

离线da376805618

资源共享呗

UID ：71951

注册：2011-01-22
登录：2014-09-27
发帖：389
等级：仿真三级

6楼发表于: 2011-09-12 16:34:50

回 4楼(tingo) 的帖子

大侠··· 谢谢你能帮我，不过我不咋懂蛙跳额。。嘿嘿···具体是哪一行呢？？我昨天把delta_t = delta_x / (2*3e8);
改为 delta_t = delta_x / (3e8);即：（ delta_t = delta_x / c ）后发现能够完全吸收了··
不知道这是为什么呢？难道说，传播一个网格要一个时间步长时，对于边界处的迭代只要一个时间步长即可？？我这样说对么？？

资源共享

发帖回复

[求助]关于一维 FDTD Mur 边界的问题 （求助）

[求助]关于一维 FDTD Mur 边界的问题（求助）