1144阅读
9回复

[求助]FDTD的离散方程写为程序时出现的问题

离线learnerlcy

UID ：43237

注册：2009-10-07
登录：2014-01-15
发帖：55
等级：仿真一级

0楼发表于: 2010-08-19 09:35:14

请问，我在编写一个三维的FDTD程序。
但是我看到两个程序，其中的一个将方程ex的计算方程中的磁场离散为：
ex[j][k] += (dtepsdy*(hz[j][k-1] - hz[j-1][k-1]) - dtepsdz*(hy[j-1][k] - hy[j-1][k-1])); 7vI ROK~
,Z7Ky*<j
另外的一个程序将ex的计算方程中的磁场离散为： ZZfi,0R
Ey[j][k]+= tmp1*(Hx(i,j,k)-Hx(i,j,k-1)-Hz(i,j,k)+Hz(i-1,j,k)) VD =f 'D
那么到底哪个是正确的，我写出来之后就是后面的这个，但是计算就是不正确。请问这两个方程哪个是正 .. frV_5yK'
n3 Rf:j^R
未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册

离线learnerlcy

UID ：43237

注册：2009-10-07
登录：2014-01-15
发帖：55
等级：仿真一级

1楼发表于: 2010-08-19 09:36:16

请问，我在编写一个三维的FDTD程序。

但是我看到两个程序，其中的一个将方程ex的计算方程中的磁场离散为：

ex[j][k] += (dtepsdy*(hz[j][k-1] - hz[j-1][k-1]) - dtepsdz*(hy[j-1][k] - hy[j-1][k-1])); -^rdB6O6j
4k6:
D2\EpL/
另外的一个程序将ex的计算方程中的磁场离散为： }F3}"Ik'L
9HlM0qE5b
Ey[j][k]+= tmp1*(Hx(i,j,k)-Hx(i,j,k-1)-Hz(i,j,k)+Hz(i-1,j,k)) t*c_70|@k
hm1.UE
那么到底哪个是正确的，我写出来之后就是后面的这个，但是计算就是不正确。请问这两个方程哪个是正确的呢？ KXu1%`x=%Z
8PG&/"K
谢谢诸位； %XXjQ5p

离线learnerlcy

UID ：43237

注册：2009-10-07
登录：2014-01-15
发帖：55
等级：仿真一级

2楼发表于: 2010-08-19 09:50:54

补充一下：ex,hy,hz分配的空间分别为：nx*（ny+1）*（nz+1）、nx*（ny-1）*（nz）、nx*（ny）*（nz-1）、

离线funnyhaha

UID ：16429

注册：2008-08-05
登录：2015-03-17
发帖：173
等级：积极交流四级

3楼发表于: 2010-08-21 04:41:37

你写的两个式子，一个是EX 一个是 EY

离线learnerlcy

UID ：43237

注册：2009-10-07
登录：2014-01-15
发帖：55
等级：仿真一级

4楼发表于: 2010-08-21 12:03:33

回 3楼(funnyhaha) 的帖子

哦，这个问题我已经解决了，其实都是正确的，至于具体原因好像语言不是很好表述，只可意会不可言传！

离线xingangahu

UID ：22197

注册：2008-11-28
登录：2016-05-31
发帖：900
等级：积极交流五级

5楼发表于: 2010-08-21 13:02:42

呵呵是程序的可读性和一般性的问题应该是space grid： dy dz是否相等的问题

离线learnerlcy

UID ：43237

注册：2009-10-07
登录：2014-01-15
发帖：55
等级：仿真一级

6楼发表于: 2010-08-21 14:00:25

回 5楼(xingangahu) 的帖子

其实就是说不能用程序中的i,j,k来恢复Yee原胞，但是我也不知道具体如何表述！不是dy, dz是否相等的问题！呵呵

离线xingangahu

UID ：22197

注册：2008-11-28
登录：2016-05-31
发帖：900
等级：积极交流五级

7楼发表于: 2010-08-21 15:26:20

恩差分格式的标准是不一样上面那个i,j,k整数点普遍的都减了1

离线learnerlcy

UID ：43237

注册：2009-10-07
登录：2014-01-15
发帖：55
等级：仿真一级

8楼发表于: 2010-08-21 21:16:49

回 7楼(xingangahu) 的帖子

呵呵，差不对就是这个意思；其实这里之所以减一是因为将边界条件考虑在了nx, ny, nz之内！

离线fuminglin85

UID ：40395

注册：2009-08-29
登录：2017-05-24
发帖：252
等级：仿真三级

9楼发表于: 2010-09-25 21:09:47

说白了还是边界条件没整对

尽人事，听天命！

发帖回复