基本信息

---

下载

   友情提示：（无需注册或充值，赞助后即可获取资源下载链接）

---

     嘿，亲！知识可是无价之宝呢，但咱这精心整理的资料也耗费了不少心血呀。小小地破费一下，绝对物超所值哦！如有下载和支付问题，请联系我们QQ(微信同号)：813200300

本次赞助数额为： 2 元　

支付宝支付

微信支付

   源码介绍

---

FDTD自由空间实现平面波

% 计算PML系数

for i = 1:8    
    Xnum = NPML - i;
    Xd = NPML;
    Xxn = Xnum / Xd;
    Xn = 0.33 * Xxn^3;
    Gi2(i) = 1.0 / (1.0 Xn);
    Gi2(IG - 1 - i) = 1.0 / (1.0 Xn);
    Gi3(i) = (1.0 - Xn) / (1.0 Xn);
    Gi3(IG - 1 - i) = (1.0 - Xn) / (1.0 Xn);
    Xxn = (Xnum - 0.5) / Xd;
    Xn = 0.25 * Xxn^3;
    Fi1(i) = Xn;
    Fi1(IG - 2 - i) = Xn;
    Fi2(i) = 1.0 / (1.0 Xn);
    Fi2(IG - 2 - i) = 1.0 / (1.0 Xn);
    Fi3(i) = (1.0 - Xn) / (1.0 Xn);
    Fi3(IG - 2 - i) = (1.0 - Xn) / (1.0 Xn);    
end

for j = 1:8
    Xnum = NPML - j;
    Xd = NPML;
    Xxn = Xnum / Xd;
    Xn = 0.33 * Xxn^3;
    Gj2(j) = 1.0 / (1.0 Xn);
    Gj2(JG - 1 - j) = 1.0 / (1.0 Xn);
    Gj3(j) = (1.0 - Xn) / (1.0 Xn);
    Gj3(JG - 1 - j) = (1.0 - Xn) / (1.0 Xn);
    Xxn = (Xnum - 0.5) / Xd;
    Xn = 0.25 * Xxn^3;
    Fj1(j) = Xn;
    Fj1(JG - 2 - j) = Xn;
    Fj2(j) = 1.0 / (1.0 Xn);
    Fj2(JG - 2 - j) = 1.0 / (1.0 Xn);
    Fj3(j) = (1 - Xn) / (1.0 Xn);
    Fj3(JG - 2 - j) = (1.0 - Xn) / (1.0 Xn);    
end