登 录
註 冊
论坛
微波仿真网
注册
登录论坛可查看更多信息
微波仿真论坛
>
时域有限差分法 FDTD
>
发一个二维TM FDTD 平面波入射到金属方柱上散 ..
发帖
回复
1
2
3
4
5
6
7830
阅读
55
回复
[
RFEDA原创
]
发一个二维TM FDTD 平面波入射到金属方柱上散射程序
离线
jinli
UID ：36593
注册：
2009-07-04
登录：
2018-02-23
发帖：
597
等级：
积极交流五级
20楼
发表于: 2010-12-24 15:23:14
谢谢分享，有fortran版的吗？
9^ mrsj
共
条评分
离线
ljhnuaa
Stay hungry， Stay foolish！
UID ：42803
注册：
2009-09-29
登录：
2023-01-03
发帖：
133
等级：
仿真二级
21楼
发表于: 2011-01-18 18:30:28
课设中要算RCS，学习学习
共
条评分
离线
xx240498420
UID ：78009
注册：
2011-05-24
登录：
2011-06-02
发帖：
5
等级：
旁观者
22楼
发表于: 2011-05-24 21:38:38
好东西，非常感谢
共
条评分
离线
heavenshu
一片空白
UID ：75800
注册：
2011-04-18
登录：
2012-05-16
发帖：
11
等级：
仿真新人
23楼
发表于: 2011-05-26 20:50:13
多谢楼主的分享，我要做一个三角形的，不过先看懂你这个，然后再根据你的来修改。谢谢哈！！！
共
条评分
离线
wwl1987
没有学不会的
UID ：78133
注册：
2011-05-26
登录：
2012-09-18
发帖：
24
等级：
仿真新人
24楼
发表于: 2011-05-26 22:57:44
谢谢分享劳动成果，雪中送炭啊，呵呵
共
条评分
离线
mywyk5522
UID ：48910
注册：
2009-12-15
登录：
2011-09-15
发帖：
78
等级：
仿真二级
25楼
发表于: 2011-06-30 12:03:58
谢谢你的程序，如果是matlab就更好了
共
条评分
离线
海聆羽
谁帮我改下程序啊 。。。万分感谢
UID ：81280
注册：
2011-08-09
登录：
2013-01-14
发帖：
13
等级：
仿真新人
26楼
发表于: 2011-08-10 18:18:10
//边界Mur ABC
}(MI}o}
//i=0; x的左边边界
nLvF^%P8
for(j=1;j<NUM_of_YAXIS-1;j++)
>`{B
{
OQ*BPmS-
ez[0][j]=ez_x_low_m1[j]+ce*(ez[1][j]-ez[0][j])-cf*(hx[0][j]-hx[0][j-1]+hy[1][j]-hy[1][j-1]);
IZ2(F,{o
//ez[0][j]=ez_x_low_m1[j]+ce*(ez[1][j]-ez[0][j]);
;YGCsLT<xt
ez_x_low_m1[j]=ez[1][j];
^qR2 !fwm<
}
d-%bRGo/
// x右边 的边界
-/)>DOgUq
for(j=1;j<NUM_of_YAXIS-1;j++)
1 >}x 9D
{
}Zuk}Og9+
ez[NUM_of_XAXIS-1][j]=ez_x_high_m1[j]+ce*(ez[NUM_of_XAXIS-2][j]-ez[NUM_of_XAXIS-1][j])-cf*(hx[NUM_of_XAXIS-1][j]-hx[NUM_of_XAXIS-1][j-1]+hy[NUM_of_XAXIS-2][j]-hy[NUM_of_XAXIS-2][j-1]);
9Su4nt `i
//ez[NUM_of_XAXIS-1][j]=ez_x_high_m1[j]+ce*(ez[NUM_of_XAXIS-2][j]-ez[NUM_of_XAXIS-1][j]);
`|rF^~6(dR
ez_x_high_m1[j]=ez[NUM_of_XAXIS-2][j];
OS - Xh-:z
}
."X}A t
|!Ryl}Oi
// y=0的左边边界
*h([ai"1-
for(i=1;i<NUM_of_XAXIS-1;i++)
@}{lp'8FYi
{
c2E*A+V#u
ez
[0]=ez_y_low_m1
+ce*(ez
[1]-ez
[0]);
YmwUl> @{
ez_y_low_m1
=ez
[1];
Gn?NY}.S
}
"/ 9EUbca
M_BG :P5
// y=NUM_of_YAXIS-1的右边界
IJ[r!&PY
for(i=1;i<NUM_of_XAXIS-1;i++)
ue8"_N
{
u$M, &Om
ez
[NUM_of_YAXIS-1]=ez_y_high_m1
+ce*(ez
[NUM_of_YAXIS-2]-ez
[NUM_of_YAXIS-1]);
qz_'v{uAj
ez_y_high_m1
=ez
[NUM_of_YAXIS-2];
MJ.K,e
a%g |E'\Jw
}
]0d j##5tJ
qp##>c31X
//四个边角
!v fbgK
ez[0][0]=ez_x_low_m1[0]+cg*(ez[1][1]-ez[0][0]);
arWP]%E0W
ez_x_low_m1[0]=ez[1][1];
./l^Iz&0
,; D$d#\"
ez[0][NUM_of_YAXIS-1]=ez_x_low_m1[NUM_of_YAXIS-1]+ce*(ez[0][NUM_of_YAXIS-2]-ez[0][NUM_of_YAXIS-1]);
(8Ptuh6\\2
ez_x_low_m1[NUM_of_YAXIS-1]=ez[0][NUM_of_YAXIS-2];
U gD|tuz]
&hnI0m=X
ez[NUM_of_XAXIS-1][0]=ez_x_high_m1[0]+ce*(ez[NUM_of_XAXIS-2][1]-ez[NUM_of_XAXIS-1][0]);
70Z#Ej
ez_x_high_m1[0]=ez[NUM_of_XAXIS-2][1];
a*kvU "]
V vFMpPi
ez[NUM_of_XAXIS-1][NUM_of_YAXIS-1]=ez_x_high_m1[NUM_of_YAXIS-1]+ce*(ez[NUM_of_XAXIS-2][NUM_of_YAXIS-2]-ez[NUM_of_XAXIS-1][NUM_of_YAXIS-1]);
Ez)hArxns
ez_x_high_m1[NUM_of_YAXIS-1]=ez[NUM_of_XAXIS-2][NUM_of_YAXIS-2];
W[VbFsI&b
'r 0kX||
共
条评分
红色六月
离线
海聆羽
谁帮我改下程序啊 。。。万分感谢
UID ：81280
注册：
2011-08-09
登录：
2013-01-14
发帖：
13
等级：
仿真新人
27楼
发表于: 2011-08-10 18:18:57
//边界Mur ABC
{YK7';_E*
//i=0; x的左边边界
e'A 1%g)
for(j=1;j<NUM_of_YAXIS-1;j++)
#b9V&/ln
{
G{74o8
ez[0][j]=ez_x_low_m1[j]+ce*(ez[1][j]-ez[0][j])-cf*(hx[0][j]-hx[0][j-1]+hy[1][j]-hy[1][j-1]);
(B0QBDj!
//ez[0][j]=ez_x_low_m1[j]+ce*(ez[1][j]-ez[0][j]);
}2hU7YWt
ez_x_low_m1[j]=ez[1][j];
B9dc *
}
] ^53Qbrv
// x右边 的边界
t{-*@8Ke
for(j=1;j<NUM_of_YAXIS-1;j++)
DX_?-jw})f
{
VDFs.;:s
ez[NUM_of_XAXIS-1][j]=ez_x_high_m1[j]+ce*(ez[NUM_of_XAXIS-2][j]-ez[NUM_of_XAXIS-1][j])-cf*(hx[NUM_of_XAXIS-1][j]-hx[NUM_of_XAXIS-1][j-1]+hy[NUM_of_XAXIS-2][j]-hy[NUM_of_XAXIS-2][j-1]);
{X10,
//ez[NUM_of_XAXIS-1][j]=ez_x_high_m1[j]+ce*(ez[NUM_of_XAXIS-2][j]-ez[NUM_of_XAXIS-1][j]);
T STkMlCG
ez_x_high_m1[j]=ez[NUM_of_XAXIS-2][j];
l4gZHMh'
}
d7G DIYH<
K=1prv2
// y=0的左边边界
@|:yK|6O
for(i=1;i<NUM_of_XAXIS-1;i++)
H =*lj.x
{
` >loleI
ez
[0]=ez_y_low_m1
+ce*(ez
[1]-ez
[0]);
BA8g[T A7K
ez_y_low_m1
=ez
[1];
9]vy#a#
}
4rLc] >
Y|6g g
// y=NUM_of_YAXIS-1的右边界
&$heW,
for(i=1;i<NUM_of_XAXIS-1;i++)
?G[=pY:=
{
`=l o. c
ez
[NUM_of_YAXIS-1]=ez_y_high_m1
+ce*(ez
[NUM_of_YAXIS-2]-ez
[NUM_of_YAXIS-1]);
[]G@l. ]W
ez_y_high_m1
=ez
[NUM_of_YAXIS-2];
b}&.IJ&40j
VCzb[.
}
z.Vf,<H
pQi |PQq
//四个边角
DQ@M?~1hp
ez[0][0]=ez_x_low_m1[0]+cg*(ez[1][1]-ez[0][0]);
f#McTC3C
ez_x_low_m1[0]=ez[1][1];
hmB`+?,z*
w7c0jIf{
ez[0][NUM_of_YAXIS-1]=ez_x_low_m1[NUM_of_YAXIS-1]+ce*(ez[0][NUM_of_YAXIS-2]-ez[0][NUM_of_YAXIS-1]);
sI MN""@Y^
ez_x_low_m1[NUM_of_YAXIS-1]=ez[0][NUM_of_YAXIS-2];
&2nICAN[
\ }J"`J\Q
ez[NUM_of_XAXIS-1][0]=ez_x_high_m1[0]+ce*(ez[NUM_of_XAXIS-2][1]-ez[NUM_of_XAXIS-1][0]);
; +1ooeU
ez_x_high_m1[0]=ez[NUM_of_XAXIS-2][1];
ZypK''&oc
/r8sL)D+
ez[NUM_of_XAXIS-1][NUM_of_YAXIS-1]=ez_x_high_m1[NUM_of_YAXIS-1]+ce*(ez[NUM_of_XAXIS-2][NUM_of_YAXIS-2]-ez[NUM_of_XAXIS-1][NUM_of_YAXIS-1]);
~^PNMZk
ez_x_high_m1[NUM_of_YAXIS-1]=ez[NUM_of_XAXIS-2][NUM_of_YAXIS-2];
lNz1|nS(Kd
NiYT%K%
请问楼主 为什么四条边的吸收边界条件所用公式不一样啊 ？
共
条评分
红色六月
离线
sujie198804
suliangjie
UID ：71338
注册：
2011-01-04
登录：
2013-04-06
发帖：
77
等级：
仿真一级
28楼
发表于: 2011-08-20 16:45:31
感谢分享！呵呵……
共
条评分
大家好
离线
xuanzehuiqu
UID ：73254
注册：
2011-03-06
登录：
2011-09-01
发帖：
8
等级：
旁观者
29楼
发表于: 2011-09-01 09:23:30
多谢，下了学习一下
共
条评分
发帖
回复
