登 录
註 冊
论坛
微波仿真网
注册
登录论坛可查看更多信息
微波仿真论坛
>
时域有限差分法 FDTD
>
发一个二维TM FDTD 平面波入射到金属方柱上散 ..
发帖
回复
1
2
3
4
5
6
7833
阅读
55
回复
[
RFEDA原创
]
发一个二维TM FDTD 平面波入射到金属方柱上散射程序
离线
jinli
UID ：36593
注册：
2009-07-04
登录：
2018-02-23
发帖：
597
等级：
积极交流五级
20楼
发表于: 2010-12-24 15:23:14
谢谢分享，有fortran版的吗？
A K/z6XGy
共
条评分
离线
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
i1b 4 J
//i=0; x的左边边界
S[W9G)KWp
for(j=1;j<NUM_of_YAXIS-1;j++)
LP5eFl`|T
{
](B& l{V
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]);
qPFG+~\c
//ez[0][j]=ez_x_low_m1[j]+ce*(ez[1][j]-ez[0][j]);
J{nyo1A
ez_x_low_m1[j]=ez[1][j];
'+ o :,6
}
Z]mM
// x右边 的边界
, aRJ!AZ
for(j=1;j<NUM_of_YAXIS-1;j++)
OoAr%
{
K^!e-Xi6
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]);
o9U0kI=W
//ez[NUM_of_XAXIS-1][j]=ez_x_high_m1[j]+ce*(ez[NUM_of_XAXIS-2][j]-ez[NUM_of_XAXIS-1][j]);
,om p F$%
ez_x_high_m1[j]=ez[NUM_of_XAXIS-2][j];
8\qCj.>S
}
Ee t+
OmT Z-*N
// y=0的左边边界
"C$!mdr7
for(i=1;i<NUM_of_XAXIS-1;i++)
d8Q_6(Ar|
{
% j[O&[s}
ez
[0]=ez_y_low_m1
+ce*(ez
[1]-ez
[0]);
R$!;J?SS
ez_y_low_m1
=ez
[1];
es.\e.HK
}
s=^r/Sz902
xG *lV|<7>
// y=NUM_of_YAXIS-1的右边界
xZ9} 8*Q&:
for(i=1;i<NUM_of_XAXIS-1;i++)
&=s|
{
]wkSAi5z*
ez
[NUM_of_YAXIS-1]=ez_y_high_m1
+ce*(ez
[NUM_of_YAXIS-2]-ez
[NUM_of_YAXIS-1]);
4i+PiD:H
ez_y_high_m1
=ez
[NUM_of_YAXIS-2];
s${T*)S@G
^D W#
}
O84v*=u A
!wLH&X$XT
//四个边角
]MqH13`)A
ez[0][0]=ez_x_low_m1[0]+cg*(ez[1][1]-ez[0][0]);
/%J&/2Wz
ez_x_low_m1[0]=ez[1][1];
Q|j@# @O 1
>gX0Ij#G
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]);
<F#*:Re_y
ez_x_low_m1[NUM_of_YAXIS-1]=ez[0][NUM_of_YAXIS-2];
]YisZE4s
LyJTK1]#
ez[NUM_of_XAXIS-1][0]=ez_x_high_m1[0]+ce*(ez[NUM_of_XAXIS-2][1]-ez[NUM_of_XAXIS-1][0]);
Uy$U8b-ov
ez_x_high_m1[0]=ez[NUM_of_XAXIS-2][1];
j61BP8E
8!o{W=m^4
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]);
1jUhG2y
ez_x_high_m1[NUM_of_YAXIS-1]=ez[NUM_of_XAXIS-2][NUM_of_YAXIS-2];
Vq\ . .!y
"Cyo<|
共
条评分
红色六月
离线
海聆羽
谁帮我改下程序啊 。。。万分感谢
UID ：81280
注册：
2011-08-09
登录：
2013-01-14
发帖：
13
等级：
仿真新人
27楼
发表于: 2011-08-10 18:18:57
//边界Mur ABC
7ck0S+N'b
//i=0; x的左边边界
,[ J '!NC1
for(j=1;j<NUM_of_YAXIS-1;j++)
vZ nO
{
#`vGg9
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]);
$MJDB
//ez[0][j]=ez_x_low_m1[j]+ce*(ez[1][j]-ez[0][j]);
t~<-4N$(
ez_x_low_m1[j]=ez[1][j];
Q 5Ghki
}
k,NU,^ &
// x右边 的边界
DO- K
for(j=1;j<NUM_of_YAXIS-1;j++)
2GD mZl
{
l>5]Wd{/
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]);
EkjK 9 2cF
//ez[NUM_of_XAXIS-1][j]=ez_x_high_m1[j]+ce*(ez[NUM_of_XAXIS-2][j]-ez[NUM_of_XAXIS-1][j]);
bJ,=yB+0
ez_x_high_m1[j]=ez[NUM_of_XAXIS-2][j];
rCSG@D.
}
H`6Jq?\
2`FsG/o\T~
// y=0的左边边界
aji~brq
for(i=1;i<NUM_of_XAXIS-1;i++)
?;KJ (@Va
{
WlQ&Yau
ez
[0]=ez_y_low_m1
+ce*(ez
[1]-ez
[0]);
h$ETH1Ue
ez_y_low_m1
=ez
[1];
K_lL\
}
H yX4ob[X
VrGb;L'[
// y=NUM_of_YAXIS-1的右边界
E]eqvT NH
for(i=1;i<NUM_of_XAXIS-1;i++)
b1 ,T!xL
{
<C.$Db&9
ez
[NUM_of_YAXIS-1]=ez_y_high_m1
+ce*(ez
[NUM_of_YAXIS-2]-ez
[NUM_of_YAXIS-1]);
Lx:9@3'7'
ez_y_high_m1
=ez
[NUM_of_YAXIS-2];
J|2H qd
v/TlXxfil
}
W'2-3J
^`>Ysc(@&
//四个边角
Q!+{MsZ
ez[0][0]=ez_x_low_m1[0]+cg*(ez[1][1]-ez[0][0]);
G4exk5
ez_x_low_m1[0]=ez[1][1];
l3 pW{p
m/F(h-?
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]);
u &d v[
ez_x_low_m1[NUM_of_YAXIS-1]=ez[0][NUM_of_YAXIS-2];
#L4Kwy
lRA Trp#T
ez[NUM_of_XAXIS-1][0]=ez_x_high_m1[0]+ce*(ez[NUM_of_XAXIS-2][1]-ez[NUM_of_XAXIS-1][0]);
!,JT9 1
ez_x_high_m1[0]=ez[NUM_of_XAXIS-2][1];
8apKp?~yW
k89gJ5B$
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]);
= a.n`3`Q
ez_x_high_m1[NUM_of_YAXIS-1]=ez[NUM_of_XAXIS-2][NUM_of_YAXIS-2];
K<7T}XzU$
= 3v 1]7 X
请问楼主 为什么四条边的吸收边界条件所用公式不一样啊 ？
共
条评分
红色六月
离线
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
多谢，下了学习一下
共
条评分
发帖
回复
