我希望能够计算一个特征函数已知的分布的密度函数。以正态分布为简单例子。
norm.char<-function(t,mu,sigma) exp((0+1i)*t*mu-0.5*sigma^2*t^2)
然后我想使用R的fft函数。但是我无法正确获取乘法常数,而且我必须重新排序结果(取值的第二半部分,然后是第一半)。我尝试了类似于
的东西。 xmax = 5
xmin = -5
deltat = 2*pi/(xmax-xmin)
N=2^8
deltax = (xmax-xmin)/(N-1)
x = xmin + deltax*seq(0,N-1)
t = deltat*seq(0,N-1)
density = Re(fft(norm.char(t*2*pi,mu,sigma)))
density = c(density[(N/2+1):N],density[1:(N/2)])
但这仍然不正确。有人知道在R中与密度计算相关的fft的好参考资料吗?显然问题在于连续FFT和离散FFT的混合。有人可以推荐一个步骤吗? 谢谢
fft
的帮助页面,我相信其中给出了方程式。 - Carl Witthoft