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

[attachment=21043] =:|fN3nJ2  
,>DaS(  
[attachment=21044]   SM<kR1bo  
r*n_#&-7  
Ez[0] = E0*exp( - (time - t0)*(time - t0)/(tw*tw) ); z%F68f73  
LC!ZeW35  
高斯点源加载在Ey上，位置分别是nx=0, nx/2处。 |o:[*2-  
C*S%aR  
当加载在x轴上nx=0处，自左向右传播。现在是波形是正确的了，可是Ey振幅达到1.0，而Hz振幅为何只达到0.0025。 v,\2$q/  
5i 6*$#OM_  
当加载在x轴上nx/2处，自左向右传播。现在是波形是正确的了，可是Ey振幅达到0.5，而Hz振幅为何只达到0.0013。  6X@]<R  
<])w@QOA#  
不知道这样的结果是否正确？ nA!Xb'y&  
_li\b-  
@iW^OVpp<8  
在Dennis M. Sullivan：electromagnetic simulation using the FDTD method.书中第三3页中图1.2中， 当源加载在nx/2时，E值振幅为一， ;[)t*yAh  
"8bxb  
而且H值振幅也为一。这就让我感到很奇怪了，不知道哪个是正确的了？ ;+) M~2 =  
QO'=O}e  
Bn[5M[  
希望大家能帮助解释一下什么样的波形才是正确的，进而才能让我确定程序是正确的？ ?89ZnH2/  
a^9-9*  
aCL_cVOMR  
\Oku<5  
< .. Fu K(SP3  
~7N>tjB  

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

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

说的有道理，是这样的。 CVsc#=w0  
Dm':D  
还是有两个问题： 75`*aAZ3  
!$'s?rnh  
    一、那为何加载在nx=0处振幅就为1.0， 加载在nx/2处振幅为0.5. v%|()Z0  
Az2HlKF"L  
    二、是不是说明Dennis M. Sullivan：electromagnetic simulation using the FDTD method.书中第三3页中图1.2是错的呢？ 28L3"c  
3AWNoXh  
      但是在测试他的程序时，结果如他的图一样。这让人糊涂了！！ RHo|&.B;+  
n237%LH[  
|qS<{WZ!h  
L3M
]06y  
_DChNX   
oI9Jp`  
h(hb?f@1:  
PgZ~of&  
U!sv6=(y@  
fL7ym,?  
kI974:e42  
B9;-Blh  
s}Q%]W  
============================================================================================

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

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

因为像两个方向传播，而能量分半，这样说也有道理。 D!7-(3R  
>s dT=6v  
我想起来了，Dennis M. Sullivan：electromagnetic simulation using the FDTD method.书中对场量是做了归一化的。 K(jo[S  
u7||
]|2  
所以导致EH的振幅都是1，它们是处于同一个数量级。而我现在是没有归一化，可能就是这个原因。
bmCp:
6  
(U@Ks )  
~]f6@n  
至于为什么加载在nx=0这点，没有特别的原因，只是这样去做了。 |ppG*ee  
\jA#RF.W  
|m{u]9  
Iu"7  
{:=sCY!  
U]Y</>xGI  
7)Cn 4{B6  
suKr//_  
aF8'^xF  
,bTpD!  
7:/gO~gI  
2F7(Y)  

lV%oIf[OB  
}ejZk bP  
================================================================

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

你做之前没对你要做的东西进行分析吗？ 6lCpf1>6@  
对于平面波在真空中传播这种简单例子而言，它是有解析解的，也就是说， aoZ|@x  
在真空中如何传播，是已知的。 :{N*Z}]  
建议你先了解自己要做的物理模型，然后再动手编写程序。 {ZrlbDQX  
一般而言，数值模拟的流程是这样的： "b~C/-W I  
1）分析要解决的问题，建立合理的模型，理论上可以得到什么结果，准备用数值模拟得到什么样的结果。。 6`7tTn?n  
2）进行分析，选定方程，和初始以及边界条件，推导迭代公式。 Pc*lHoVL  
3）设置模型参数，可以代入特例，与已知的理论结果进行对比。 aPzn4}~/_  
4）代入实际参数，对模拟结果进行预先估计。 .Sn{a}XP4  
5）分析模拟结果。

斑竹 高见！！ i7dDklj4  
E]q>ggeNH  
又学习了。 xiW}P% bf  
GIlaJ!/  
oN}\bK  
我刚接触FDTD还有很多东西不了解，现在只是在尝试一般的现象。至于建模我还不怎么理解是怎么回事？ R V#w0 r  
MQ =x:p{  
是不是像葛德彪书后的附录二：1.建模，那样。我现在还没能写出从文件读入数据的函数， }P2*MrkcHB  
w(HVC  
}zA|M9%E  
-, +o*BP  
p7+{xXf  
m !:F/?B  
VY/r2o#  
^fKKsfIf  
hEjvtfM9\-  
pRfKlTU\  
O[Z$~  
====================================================================================

这样的话，那说明你对编程了解很少啊。有空还是先看看编程方面的书吧，很简单的，基本上一个星期就差不多搞定了。
;f3))x  
一般来说： >48)@sS  
1）定义变量，包括整型，double,float，数组，自定义结构等。 &)Wm rF  
2）结构，if , for , while wW/wvC-  
3）熟悉一些基本的函数，包括数学函数，文件操作等。 +Y~,1ai 5^  
4）字符串操作，数据类型间的转换等等。 9OF(UFgS  
)Z?\9'6e4  
大体上来说，你找本书看个一天就差不多了，然后再用两三天动手编写几个例子，最好规范化点，以后也好调用。 @(st![i+  
比如就如你说的读取文件，保存文件等等。 >m9ge`!9  
>*w(YB]/$V  
不知道你用什么语言，c/c++ ,fortran还是matlab，不管哪个仅仅读取文件也就廖廖几行代码，五六分钟就可以弄好了。