zoukankan      html  css  js  c++  java
  • fdtd二维解说(二)边界条件

    完美匹配层(PML) 

     

     

    figure();
    hold on;
    d=20;
    N=7;
    N_PML=2;
    for i=0:N
        plot([i*d,i*d],[0,N*d],'-k'); hold on;
        plot([0,N*d],[i*d,i*d],'-k'); hold on;
        plot(ones(N)*i*d,((1:N)-1/2)*d,'kx'); hold on;
        plot(((1:N)-1/2)*d,ones(N)*i*d,'kx'); hold on;
    end
    shift=5;
    plot([0,N*d,N*d,0,0],[0,0,N*d,N*d,0],'--r','Linewidth',2);
    plot([N_PML*d-shift,(N-N_PML)*d+shift,(N-N_PML)*d+shift,N_PML*d-shift,N_PML*d-shift], ...
        [N_PML*d-shift,N_PML*d-shift,(N-N_PML)*d+shift,(N-N_PML)*d+shift,N_PML*d-shift],'--r','Linewidth',2);
    for i=1:N
        plot(((1:N)-1/2)*d,ones(N)*(i-1/2)*d,'ko'); hold on;
    end
    axis([-20,160,-20,160]);
    

      

    figure();
    hold on;
    d=20;
    N=7;
    N_PML=2;
    shift=5;
    for i=0:N
        plot([i*d,i*d]-shift,[0,N*d]+shift,'-k'); hold on;
        plot([0,N*d]-shift,[i*d,i*d]+shift,'-k'); hold on;
        plot(ones(N)*i*d,((1:N)-1/2)*d,'kV','MarkerFaceColor',[.49 1 .63]); hold on;
        plot(((1:N)-1/2)*d,ones(N)*i*d,'k>','MarkerFaceColor',[ 1 .63 .49]); hold on;
    end
    for i=1:N
        plot(((1:N)-1/2)*d,ones(N)*(i-1/2)*d,'ko'); hold on;
    end
    
    
    plot([0,N*d,N*d,0,0]-shift,[0,0,N*d,N*d,0]+shift,'--r','Linewidth',3);
    plot([N_PML*d-shift, (N-N_PML)*d-shift, (N-N_PML)*d-shift, N_PML*d-shift, N_PML*d-shift], ...
        [N_PML*d+shift,N_PML*d+shift,(N-N_PML)*d+shift,(N-N_PML)*d+shift,N_PML*d+shift],'--r','Linewidth',3);
    
    axis([-20,160,-20,160]);
    

      

  • 相关阅读:
    Linux内核链表——看这一篇文章就够了
    2的幂和按位与&——效率
    fgets注意事项
    GDB TUI
    GDB TUI
    linux shell命令行选项与参数用法详解
    What are the differences between Perl, Python, AWK and sed
    What is the difference between sed and awk
    /proc/sysrq-trigger
    C++ Sqlite3的基本使用
  • 原文地址:https://www.cnblogs.com/Iknowyou/p/6930900.html
Copyright © 2011-2022 走看看