发帖回复

1 2 3 4 567

12611阅读
68回复

[RFEDA原创]Electromagnetic Simulation Using The FDTD Method Sullivan书中PML部分中译加Matlab程序

离线yumo145

UID ：56254

注册：2010-04-01
登录：2012-12-27
发帖：25
等级：仿真新人

50楼发表于: 2010-04-26 08:50:54

感谢楼主，正是所需

离线inter

UID ：42650

注册：2009-09-28
登录：2020-04-23
发帖：120
等级：仿真二级

51楼发表于: 2010-05-03 12:45:11

非常不错！

我觉得我什么都不懂！

离线rain_1216

UID ：60058

注册：2010-05-20
登录：2010-06-12
发帖：1
等级：旁观者

52楼发表于: 2010-05-20 18:24:25

我怎么下载不了啊？请问如何下载？

离线lwfz

UID ：7619

注册：2008-01-15
登录：2022-08-21
发帖：1081
等级：积极交流六级

53楼发表于: 2010-05-28 08:11:30

强烈感谢楼主

离线h19870101l

UID ：62933

注册：2010-07-05
登录：2010-07-06
发帖：2
等级：旁观者

54楼发表于: 2010-07-05 16:51:52

好东西，学习学习！

离线雪山

UID ：63052

注册：2010-07-08
登录：2011-12-20
发帖：25
等级：仿真新人

55楼发表于: 2010-07-08 11:38:26

好书
很有用的一本学习FDTD的书

离线i2000s

UID ：50202

注册：2010-01-01
登录：2010-07-27
发帖：32
等级：仿真新人

56楼发表于: 2010-07-27 06:17:20

牛人阿
感谢一下

离线samuelp

UID ：63797

注册：2010-07-23
登录：2010-10-19
发帖：28
等级：仿真二级

57楼发表于: 2010-07-29 18:56:10

来一个python版的：

from pylab import *
from enthought.mayavi import mlab

IE = 60;   ic = IE/2; NSTEPS = 250; NPML = 8
JE = 60;   jc = JE/2
ddx = 0.01; dt = ddx/6.0e8; epsz = 8.8e-12

dz = zeros((IE,JE)); hx = zeros((IE,JE)); ihx = zeros((IE,JE))
ez = zeros((IE,JE)); hy = zeros((IE,JE)); ihy = zeros((IE,JE))
ga = ones((IE,JE))
gi2 = ones(IE); gi3 = ones(IE); fi1= zeros(IE); fi2 = ones(IE); fi3 = ones(IE)
gj2 = ones(JE); gj3 = ones(JE); fj1= zeros(JE); fj2 = ones(JE); fj3 = ones(JE)

t0 = 20.0; spread = 6.0; T = 0

xnum = NPML-arange(NPML)

xn = 0.33*(xnum/float(NPML))**3
gi2[0:NPML] =      1.0/(1.0+xn); gi2[-1:-1-NPML:-1] =      1.0/(1.0+xn)
gi3[0:NPML] = (1.0-xn)/(1.0+xn); gi3[-1:-1-NPML:-1] = (1.0-xn)/(1.0+xn)

gj2[0:NPML] =      1.0/(1.0+xn); gj2[-1:-1-NPML:-1] =      1.0/(1.0+xn)
gj3[0:NPML] = (1.0-xn)/(1.0+xn); gj3[-1:-1-NPML:-1] = (1.0-xn)/(1.0+xn)

xn = 0.25*((xnum-0.5)/float(NPML))**3
fi1[0:NPML] =               xn; fi1[-2:-2-NPML:-1] =               xn
fi2[0:NPML] =      1.0/(1.0+xn); fi2[-2:-2-NPML:-1] =      1.0/(1.0+xn)
fi3[0:NPML] = (1.0-xn)/(1.0+xn); fi3[-2:-2-NPML:-1] = (1.0-xn)/(1.0+xn)

fj1[0:NPML] =               xn; fj1[-2:-2-NPML:-1] =               xn
fj2[0:NPML] =      1.0/(1.0+xn); fj2[-2:-2-NPML:-1] =      1.0/(1.0+xn)
fj3[0:NPML] = (1.0-xn)/(1.0+xn); fj3[-2:-2-NPML:-1] = (1.0-xn)/(1.0+xn)

for n in range(NSTEPS):
    T += 1
    for i in range(1,IE):
        for j in range(1,JE):
            dz[j] = gi3*gj3[j]*dz[j]+gi2*gj2[j]*0.5*(hy[j]-hy[i-1][j]-hx[j]+hx[j-1]) V|}9d:&O
    pulse = sin(2*pi*1500*1e6*dt*T); dz[ic][jc] = pulse :tdx:
    ez = ga*dz Q'JEDH\
    for i in range(IE): ez[0] = 0.0; ez[-1] = 0.0 V-A^9AAPm
    for j in range(JE): ez[0][j] = 0.0; ez[-1][j] = 0.0 NUH;GMj,,
    for i in range(IE): sxgR;gf6
        for j in range(JE-1): nen6!bw4
            curl_e = ez[j]-ez[j+1] 7t5X
            ihx[j] = ihx[j]+fi1*curl_e oF.Fg<p(
            hx[j] = fj3[j]*hx[j]+fj2[j]*0.5*(curl_e+ihx[j]) k(=\&T
    for i in range(IE-1): uA C:&
        for j in range(JE-1): 0o At=S
            curl_e = ez[i+1][j]-ez[j] p$f#W
            ihy[j] = ihy[j]+fj1[j]*curl_e 10}oaL S
            hy[j] = fi3*hy[j]+fi2*0.5*(curl_e+ihy[j]) ,p1 (0i
x, y = np.ogrid[0:IE,0:JE] pPo?5s
pl = mlab.surf(x, y, ez, warp_scale="auto") rRcfZZ~` M
vc6UA%/f
Balanis书上用的是Berenger PML看来FDTD主要是在边界条件上变花样

图片:fd2d_3.2.JPG

离线yige2b

UID ：60536

注册：2010-05-27
登录：2010-08-19
发帖：9
等级：旁观者

58楼发表于: 2010-08-16 07:43:13

离线applenic21

UID ：66488

注册：2010-09-16
登录：2014-09-20
发帖：6
等级：旁观者

59楼发表于: 2010-10-06 22:31:56

很有用的東西
謝謝分享啦

发帖回复