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

[attachment=21043] IV0[!D  
m%pBXXfGYj  
[attachment=21044]   3L>d!qD  
tRteyNA  
Ez[0] = E0*exp( - (time - t0)*(time - t0)/(tw*tw) ); `fz,Lh*v  
P.cO6+jGR  
高斯点源加载在Ey上，位置分别是nx=0, nx/2处。 Txo@U  
ryg4hHspl  
当加载在x轴上nx=0处，自左向右传播。现在是波形是正确的了，可是Ey振幅达到1.0，而Hz振幅为何只达到0.0025。 j&WL*XP&5  
),=@q+{E{  
当加载在x轴上nx/2处，自左向右传播。现在是波形是正确的了，可是Ey振幅达到0.5，而Hz振幅为何只达到0.0013。  QQFf5^  
g5y;?fqJ  
不知道这样的结果是否正确？ 8V-,Xig;`  
fD{II+T  
%*kLEA*v  
在Dennis M. Sullivan：electromagnetic simulation using the FDTD method.书中第三3页中图1.2中， 当源加载在nx/2时，E值振幅为一， vIVr@1S  
r0\?WoF2C  
而且H值振幅也为一。这就让我感到很奇怪了，不知道哪个是正确的了？ A[`G^$  
[![ G7H%f  
J
*j5#V];  
希望大家能帮助解释一下什么样的波形才是正确的，进而才能让我确定程序是正确的？ [SkKz>rC  
cx{T  '1  
g`[`P@  
'.Z4 hHX  
< .. JS!*2*Wr  
D<t~e$H  

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

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

说的有道理，是这样的。 tCG76LH  
mLV[uhq   
还是有两个问题： )0 W`  
ikD1N  
    一、那为何加载在nx=0处振幅就为1.0， 加载在nx/2处振幅为0.5. Z6p5*+  
M$4[)6Y  
    二、是不是说明Dennis M. Sullivan：electromagnetic simulation using the FDTD method.书中第三3页中图1.2是错的呢？ :#~U<C@o  
u"1rF^j6k  
      但是在测试他的程序时，结果如他的图一样。这让人糊涂了！！ wIB`%V  
=3ioQZ^Vz  
b89a)k>^g  
Cqii}  
RwI[R)k  
7Z ;?b0W  
2dp>Z",  
?<Wb@6k
h`  
L"jY+{oLIJ  
Z/;Xl~  
0|&\'{  
uj>WgU  
FGWN}&K  
============================================================================================

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

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

因为像两个方向传播，而能量分半，这样说也有道理。 ; F"g$_D0  
h+g_rvIG*  
我想起来了，Dennis M. Sullivan：electromagnetic simulation using the FDTD method.书中对场量是做了归一化的。 /NI;P]s.  
84& $^lNV  
所以导致EH的振幅都是1，它们是处于同一个数量级。而我现在是没有归一化，可能就是这个原因。 rr],DGg+B]  
/T0F"e)Ci  
`EA\u]PwQ  
至于为什么加载在nx=0这点，没有特别的原因，只是这样去做了。
FrGgga$  
hF~n)oQ  
}HYbS8'  
2*;~S44  
"c
%0P"u  
nv|NQ Tk  
Fo5FNNiID  
P64PPbP  
_)3|f<E_t)  
:+^lJ&{U  
5H<m$K4z  
Qb-M6ihcc  
\kL3.W_  
'}53f2%gKa  
================================================================

有些程序里用的是国际单位制，有些是用的高斯制，各有各的方便。 0 1rK8jX  
高斯下面电场，磁场幅度都是一样的了。

你做之前没对你要做的东西进行分析吗？ gR;i(81U  
对于平面波在真空中传播这种简单例子而言，它是有解析解的，也就是说， wlqksG[B  
在真空中如何传播，是已知的。 ^6V[=!& H  
建议你先了解自己要做的物理模型，然后再动手编写程序。 :4/3q|cn  
一般而言，数值模拟的流程是这样的： [RhO$c$[\  
1）分析要解决的问题，建立合理的模型，理论上可以得到什么结果，准备用数值模拟得到什么样的结果。。 ea 'D td  
2）进行分析，选定方程，和初始以及边界条件，推导迭代公式。 YjKxb9  
3）设置模型参数，可以代入特例，与已知的理论结果进行对比。 g8% &RG  
4）代入实际参数，对模拟结果进行预先估计。 f.$af4 u  
5）分析模拟结果。

斑竹 高见！！ )M^ gT}M  
,fRq5"?  
又学习了。 `H+lPM66  
4&iCht =  
Z30A{6}  
我刚接触FDTD还有很多东西不了解，现在只是在尝试一般的现象。至于建模我还不怎么理解是怎么回事？ y_[vr:s5pG  
I`#JwMU;m  
是不是像葛德彪书后的附录二：1.建模，那样。我现在还没能写出从文件读入数据的函数， J~- 4C)  
 AOx[  
w2J<WC+_<  
6w77YTJ  
fB,_9K5i  
LH6vLuf  
:zR!/5  
D&zle~" J  
@o.I;}*
N  
;n},"&  

Mb=" Te>|  
====================================================================================

这样的话，那说明你对编程了解很少啊。有空还是先看看编程方面的书吧，很简单的，基本上一个星期就差不多搞定了。 V1JIht>Opo  
一般来说： ]
s748+  
1）定义变量，包括整型，double,float，数组，自定义结构等。 L_s:l9!r  
2）结构，if , for , while 9k=3u;$v  
3）熟悉一些基本的函数，包括数学函数，文件操作等。 bu"!jHPB  
4）字符串操作，数据类型间的转换等等。 yOKI*.}  
jA/w|\d!  
大体上来说，你找本书看个一天就差不多了，然后再用两三天动手编写几个例子，最好规范化点，以后也好调用。 Q5_o/wk  
比如就如你说的读取文件，保存文件等等。 1i] ^{;]  
[trwBZ^D~  
不知道你用什么语言，c/c++ ,fortran还是matlab，不管哪个仅仅读取文件也就廖廖几行代码，五六分钟就可以弄好了。