1889阅读
4回复

[已解决]MUR边界发散是什么原因呢？

UID ：33555

注册：2009-05-25
登录：2010-04-25
发帖：56
等级：仿真一级

0楼发表于: 2010-03-03 18:30:15

        编写TM波真空中传播程序中加入一阶MUR边界，加入右边界或上边界后电场值E很大，这就是传说的发散吧，会是什么原因引起的呢？大家帮看看，谢谢了。
        边界代码如下：
aa = (c*delta_t-delta_s)/(c*delta_t+delta_s);
%左边界x=1
     for j=2:max_space_z-1
         Ez(1,,j) = ez(2,,j)+aa*( Ez(2,,j)-ez(1,,j));
     end
%右边界x=max_space_z
    for j=2:max_space_z-1
         Ez(max_space_x,,j) = ez(max_space_x-1,,j)+aa*( Ez(max_space_x-1,,j)-ez(max_space_x,,j));
    end
%下边界y=1
    for i=2:max_space_x-1
        Ez(i,1) = ez(i,2)+aa*(Ez(i,2)-ez(i,1));
    end
%上边界y=max_space_x
    for i=2:max_space_x-1
        Ez(i,max_space_z) = ez(i,max_space_z-1)+aa*(Ez(i,max_space_z-1)-ez(i,max_space_z));
    end
[backcolo ..

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