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

最近在编写一个3DFDTD的C程序，在设置波源时遇到了问题，我的想法是设置一个正弦点波源（软源），程序是： ZV07;`I  
T$V8n_;  
void powerSource(int timestep)  9Lqz:4}  
{ gX* &RsF  
source = 10 * sin(2*4*atan(1)*6e9*dt*n); frWY8&W^H  
$% W.=a'5  
zS?DXE  
Hz[isource][jsource][ksource] += source; {|e7^_ke  
} @*$"6!3s5  
hiWs:Yq  
结果出来的结果发现并不是正弦波，如果把 += 改成 = 的话结果就是错的，算出来电场值有10几个数量级 .. 2%]#rZ  
zEeix,IU  

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

应该是＝号，表示硬源。 &s&Ha{(!w  
w{Y:p[}  
建议把正选波，变为渐近的正玄波试试。 5OC3:%g  
"3NE%1T  
source = sin(...) * ramp_function

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

ramp_function : 71ctjU`U2  
2nGQD{  
f(t)=0, t<0 > %U
  
8'B   
f(t) = \alpha t, when 0<t, t<C A12#v,  
m# ]VdO'f  
f(t) = 1 when t>C T};fy+iq  
,#PeK(  
\alpha 和 C 都是常量。通常C是两三个周期。 OK}+:Y  
EJrn4QOs  
使用硬源和误差无关。

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