[求助]咨询一下有关波源的问题

最近在编写一个3DFDTD的C程序，在设置波源时遇到了问题，我的想法是设置一个正弦点波源（软源），程序是： 7^Uv7<pw  
sT.ss$HY9,  
void powerSource(int timestep)  TvM~y\s  
{ N[hG8f  
source = 10 * sin(2*4*atan(1)*6e9*dt*n); 0]L"H<W  
K:M8h{Ua  
^Pf WG*  
Hz[isource][jsource][ksource] += source; K}y f>'O  
} ?M2J wAK5  
0J|3kY-n>  
结果出来的结果发现并不是正弦波，如果把 += 改成 = 的话结果就是错的，算出来电场值有10几个数量级 .. ]|@^1we  
@iiT
<
 

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

应该是＝号，表示硬源。 >
j/w@Fj  
paK2xX
8E  
建议把正选波，变为渐近的正玄波试试。 Q?vlfZR`8  
TxD#9]Q`  
source = sin(...) * ramp_function

谢谢~ramp_function这种函数一般都是什么样的表达式啊，用硬源的话会不会使计算误差变得很大呢。。

ramp_function : xdPx{"C 3  
BtZyn7a  
f(t)=0, t<0 GgU/!@  
SbZ6t$"  
f(t) = \alpha t, when 0<t, t<C st*gs-8jJ;  
f);FoVa6  
f(t) = 1 when t>C \8tsDG(1 '  
z:O8Ls^\T  
\alpha 和 C 都是常量。通常C是两三个周期。 >_} I.\X  
!D6]J
PX  
使用硬源和误差无关。

我是设置软源计算的，结果都是好的。可以联系我详细讨论： 5;?yCWc  
sannxi@163.com