[求助]请大伙帮忙见证一维TEM计算结果是否正确！！

[attachment=21043] b;jr;I  
ZWXA%u7V  
[attachment=21044]   }=$>w@mJ  
ni~1)"U.  
Ez[0] = E0*exp( - (time - t0)*(time - t0)/(tw*tw) ); v4Mn@e_#c  
Lm1 -  
高斯点源加载在Ey上，位置分别是nx=0, nx/2处。 !(nFq9~~Q  
Q|7l!YTzVu  
当加载在x轴上nx=0处，自左向右传播。现在是波形是正确的了，可是Ey振幅达到1.0，而Hz振幅为何只达到0.0025。 R?5v//[  
?M9?GodbP.  
当加载在x轴上nx/2处，自左向右传播。现在是波形是正确的了，可是Ey振幅达到0.5，而Hz振幅为何只达到0.0013。  5-2#H?:U  
}!*CyO*  
不知道这样的结果是否正确？ |{ TVW  
dh K<5E  
d<_#Q7]I4  
在Dennis M. Sullivan：electromagnetic simulation using the FDTD method.书中第三3页中图1.2中， 当源加载在nx/2时，E值振幅为一， SbK6o:[  
vugGMP;D(  
而且H值振幅也为一。这就让我感到很奇怪了，不知道哪个是正确的了？ /ei(Q'pc[  
f=paa/k0  
\#7@"~<  
希望大家能帮助解释一下什么样的波形才是正确的，进而才能让我确定程序是正确的？ D*}_

L   
n@_aTY  
uq\[^  
[5i}C K_=  
< .. }UX0 eI4  
Ju"K"  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册

E和H 相差一个ETA  波阻抗的倍数  377.*****

说的有道理，是这样的。 TRFza}4:i  
eKe[]/}e9  
还是有两个问题： o"g<Vz  
3+8{Y  
    一、那为何加载在nx=0处振幅就为1.0， 加载在nx/2处振幅为0.5. 0}-
MWbG  
e7n`fEpO  
    二、是不是说明Dennis M. Sullivan：electromagnetic simulation using the FDTD method.书中第三3页中图1.2是错的呢？ {4r }jH  
`gdk,L]  
      但是在测试他的程序时，结果如他的图一样。这让人糊涂了！！ ?3do-tTp  
Q#pgl
  
Vo9FlYj  
ZmS ]4WM<  
+oiuulA  
m0]Lc{  
UQq Qim  
u(~(+1W  
]Uul~T
 
4YA1~7R  
p*|Ct  
%Vhj<
gN  
LT:8/&\  
============================================================================================

没看过这本书 CwsC)]{/o  
你的源不是向一个方向传输的吧   点源在三维是球面波  波会向各个方向传播 b#:!b  
你加在nx/2的地方   波向+x -x两个方向传输的吧  不是从左向右传输的  那样能量会变成原先一半？  到Ey分量就是1/sqrt(2)?   我在一个类似的问题上也困惑过。希望大家讨论一下。

还有你为什么在nx=0的地方加源呢？   边界在哪里处理的？

因为像两个方向传播，而能量分半，这样说也有道理。 & _mp!&5XV  
_BBs{47{E  
我想起来了，Dennis M. Sullivan：electromagnetic simulation using the FDTD method.书中对场量是做了归一化的。 31~Rs?~f(  
OjfumZL#  
所以导致EH的振幅都是1，它们是处于同一个数量级。而我现在是没有归一化，可能就是这个原因。 =x}p>#o,J  
QZamf lk  
Gw?$.@L'I6  
至于为什么加载在nx=0这点，没有特别的原因，只是这样去做了。
uui3jZ:  
#lO~n.+P  
Z,5B(Xj  
8 J;\Z  
d@>1m:p  
,_RPy2N  
0'~Iv\s  
Pze{5!  
&,C;_3   
v}BXH4&Y  
Qs;MEt1  
PR~9*#"v..  
0p1~!X=I  
jt6,id)&  
================================================================

有些程序里用的是国际单位制，有些是用的高斯制，各有各的方便。 <U~at+M  
高斯下面电场，磁场幅度都是一样的了。

你做之前没对你要做的东西进行分析吗？ J/B`c(  
对于平面波在真空中传播这种简单例子而言，它是有解析解的，也就是说， +a0` ,Jc  
在真空中如何传播，是已知的。 M3Oqto
<8"  
建议你先了解自己要做的物理模型，然后再动手编写程序。 r>cN,C
 
一般而言，数值模拟的流程是这样的： W(\^6S)  
1）分析要解决的问题，建立合理的模型，理论上可以得到什么结果，准备用数值模拟得到什么样的结果。。 E_T!|Q.  
2）进行分析，选定方程，和初始以及边界条件，推导迭代公式。 ED^0t  
3）设置模型参数，可以代入特例，与已知的理论结果进行对比。 !Z<=PdI1Ys  
4）代入实际参数，对模拟结果进行预先估计。 9A ?)n<3d  
5）分析模拟结果。

斑竹 高见！！ vbn=ywz  
')0@J`  
又学习了。 #4./>}G  
U[02$gd0l  
-fXQ62:S  
我刚接触FDTD还有很多东西不了解，现在只是在尝试一般的现象。至于建模我还不怎么理解是怎么回事？ Vjs'|%P7  
qxDMDMN  
是不是像葛德彪书后的附录二：1.建模，那样。我现在还没能写出从文件读入数据的函数， S3l^h4  
:"MHmm=uU8  
^C1LQZ  
3}!u8,P  
e4,SR(O>  
8lzoiA_9  
kmX9)TMVO  
e5 L_<V^Jo  
|9=A"092{  
!Ve0:$  
G;flj}z  
====================================================================================

这样的话，那说明你对编程了解很少啊。有空还是先看看编程方面的书吧，很简单的，基本上一个星期就差不多搞定了。 qvSYrnpn  
一般来说： ;W|NG3_y  
1）定义变量，包括整型，double,float，数组，自定义结构等。 Qgj# k  
2）结构，if , for , while tRJ5IX##L  
3）熟悉一些基本的函数，包括数学函数，文件操作等。  e|!'  
4）字符串操作，数据类型间的转换等等。 O^`EuaL  
EN\cwa#FU  
大体上来说，你找本书看个一天就差不多了，然后再用两三天动手编写几个例子，最好规范化点，以后也好调用。 38p"lT  
比如就如你说的读取文件，保存文件等等。 TT/H"Ri}Jp  
NGSts\D'}  
不知道你用什么语言，c/c++ ,fortran还是matlab，不管哪个仅仅读取文件也就廖廖几行代码，五六分钟就可以弄好了。