[已解决]Sullivan 写的两个程序FDTD1D_1.5和FDTD_2.1模拟相同的情景为什么输出的结果不同呀

离线xiahuiz

UID ：20367

注册：2008-10-29
登录：2012-04-07
发帖：30
等级：仿真一级

10楼发表于: 2009-05-07 17:28:15

回 9楼(gwzhao) 的帖子

真是大好人呀，谢谢谢谢！

离线gwzhao

方恨少

UID ：17098

注册：2008-08-24
登录：2019-01-09
发帖：1374
等级：荣誉管理员

11楼发表于: 2009-05-07 18:28:13

回 10楼(xiahuiz) 的帖子

2.1里E(k) = dielectric_space_a(i)*(D(k)-I(k));
I(k) = I(k)+dielectric_space_b(i)*E(k);
这两个公式怎么推导出来的？

逆流而上

离线xiahuiz

UID ：20367

注册：2008-10-29
登录：2012-04-07
发帖：30
等级：仿真一级

12楼发表于: 2009-05-07 19:08:32

回 11楼(gwzhao) 的帖子

2.1里的公式推导放在附件里，非常感谢！

描述:公式推导
附件:

公式推导.doc (260 K) 下载次数:8

离线gwzhao

方恨少

UID ：17098

注册：2008-08-24
登录：2019-01-09
发帖：1374
等级：荣誉管理员

13楼发表于: 2009-05-07 22:15:27

回 12楼(xiahuiz) 的帖子

看了一下公式推导，电流卷积，显然公式没有问题，然后再看了一下参数，也没问题，只能调试程序了。
运行1。5，结果比较正常。
运行2。1，发现在0，100的区域之间，波形就开始衰减，而在这个区域内是没有加入损耗介质的，却出现这种现象，那肯定是程序的问题了。
同时显示电场和电流，发现电流开始就不为零，显然不对，找到电流赋值的地方，发现公式
I(k) = I(k)+dielectric_space_b(i)*E(k);
明显dielectric_space_b(i) 应该为dielectric_space_b(k)
因为 i 一直为200。

修改代码，重新运行，结果应该就对了。
你看一下。

共1条评分

cem-uestc

技术分 +1

积极参与讨论+技术分论坛感谢您的参与

2009-05-07

逆流而上

离线gwzhao

方恨少

UID ：17098