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

[attachment=21043] LrT? ]o  
Qxfds`4V9i  
[attachment=21044]   -jzoGzC3  
u2OrH3E4E3  
Ez[0] = E0*exp( - (time - t0)*(time - t0)/(tw*tw) ); ,$Cr9R&/  
T
$sm}=  
高斯点源加载在Ey上，位置分别是nx=0, nx/2处。 m%nRHT0KAf  

"xNP"S  
当加载在x轴上nx=0处，自左向右传播。现在是波形是正确的了，可是Ey振幅达到1.0，而Hz振幅为何只达到0.0025。 < lUpvr  
-p"}K~lt:  
当加载在x轴上nx/2处，自左向右传播。现在是波形是正确的了，可是Ey振幅达到0.5，而Hz振幅为何只达到0.0013。  /9,y+"0SQz  
22*~CIh~x  
不知道这样的结果是否正确？ wq|7sk{  
'w=aLu5dY  
2FY]o~@  
在Dennis M. Sullivan：electromagnetic simulation using the FDTD method.书中第三3页中图1.2中， 当源加载在nx/2时，E值振幅为一， T7|=`~  
d@tf+_Ih  
而且H值振幅也为一。这就让我感到很奇怪了，不知道哪个是正确的了？ U
02  
o;21|[z  
ktkS$  
希望大家能帮助解释一下什么样的波形才是正确的，进而才能让我确定程序是正确的？ dsOt(yNo  
*zMt/d*<&  
0+k..l  
Y6T{/!  
< .. ?Yx2q_KZk  
/Ny#+$cfk  

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

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

说的有道理，是这样的。 R?,v:S&i7;  
!6sR|c"~j  
还是有两个问题： PR AP~P&^  
Sn 7h$  
    一、那为何加载在nx=0处振幅就为1.0， 加载在nx/2处振幅为0.5. 3RI6+Cgmn  
%KN2iNq  
    二、是不是说明Dennis M. Sullivan：electromagnetic simulation using the FDTD method.书中第三3页中图1.2是错的呢？ [} d39  
a+CJJ3T-  
      但是在测试他的程序时，结果如他的图一样。这让人糊涂了！！ 7
l09  
WX`wz>KK^  
t$p%UyVE  
-UAMHd}4  
1K)9fMr]  
|Q@(<'8=  
#6mw CA|  
Uqx@9z(  
(lq%4h  
qlg.\H:W~  
L_+Fin  
[#SO}'1n  
0"N4WH O  
============================================================================================

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

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

因为像两个方向传播，而能量分半，这样说也有道理。 A=\:b^\  
ZKoISuM  
我想起来了，Dennis M. Sullivan：electromagnetic simulation using the FDTD method.书中对场量是做了归一化的。 H.!\j&4j  
n~?n+\.&a  
所以导致EH的振幅都是1，它们是处于同一个数量级。而我现在是没有归一化，可能就是这个原因。 *ZV=4[#bT  
xC-&<s  
]Jx_bs
~g  
至于为什么加载在nx=0这点，没有特别的原因，只是这样去做了。 e<HHgC#J  
\@GKVssw  
]E}eM@xdD  
p!H'JNG  
1\%2@NR  
?9:~d#p  
3mLtnRX[m  
fzW!-  
&Q[|FO;[  
VK^m]??s_  
.v_-V?7  
5HW'nhE  
~cb7]^#u1l  
;Zw? tU  
================================================================

有些程序里用的是国际单位制，有些是用的高斯制，各有各的方便。 BY`vs+]XY  
高斯下面电场，磁场幅度都是一样的了。

你做之前没对你要做的东西进行分析吗？ VWrb`p@  
对于平面波在真空中传播这种简单例子而言，它是有解析解的，也就是说， bMB*9<c~  
在真空中如何传播，是已知的。 <RuLIu  
建议你先了解自己要做的物理模型，然后再动手编写程序。 $ D.*r*c6  
一般而言，数值模拟的流程是这样的： 5Zn:$?7  
1）分析要解决的问题，建立合理的模型，理论上可以得到什么结果，准备用数值模拟得到什么样的结果。。 [C{oj*"c]  
2）进行分析，选定方程，和初始以及边界条件，推导迭代公式。 OM?FpRVU8  
3）设置模型参数，可以代入特例，与已知的理论结果进行对比。 0( q:K6zI}  
4）代入实际参数，对模拟结果进行预先估计。 1|:'jK#gE  
5）分析模拟结果。

斑竹 高见！！ keStK8  
\HQ.Pwr 6  
又学习了。 J|@D @\?7  
-B(p8YH  

QL2 `X2  
我刚接触FDTD还有很多东西不了解，现在只是在尝试一般的现象。至于建模我还不怎么理解是怎么回事？ ,?`kYPZ  
Zv*uUe  
是不是像葛德彪书后的附录二：1.建模，那样。我现在还没能写出从文件读入数据的函数， 7xR:\FBa^  
B\qy:nr j  
x(p/9$.#  
=W>a~e]/  
vZs~=nfi#|  
KC%&or  
E,cQ9}/  
y=jZ8+M   
t:xTmK&vt  
P*OT&q  
V#R; -C  
====================================================================================

这样的话，那说明你对编程了解很少啊。有空还是先看看编程方面的书吧，很简单的，基本上一个星期就差不多搞定了。 -Rr !J37  
一般来说： dP>FXgY  
1）定义变量，包括整型，double,float，数组，自定义结构等。 LlVbY=EX7  
2）结构，if , for , while NS;LFeGD  
3）熟悉一些基本的函数，包括数学函数，文件操作等。 H ]!P[?  
4）字符串操作，数据类型间的转换等等。 hW6Ksn,*  
uD[T l  
大体上来说，你找本书看个一天就差不多了，然后再用两三天动手编写几个例子，最好规范化点，以后也好调用。 H\a\xCP3  
比如就如你说的读取文件，保存文件等等。 a4Z e!l(  
563ExibH  
不知道你用什么语言，c/c++ ,fortran还是matlab，不管哪个仅仅读取文件也就廖廖几行代码，五六分钟就可以弄好了。