基本上,我想做以下几点:
t = np.linspace(-5, 5, 100)
y = np.exp(-(t**2)/2)
p = numpy.fft.ffti(y)
x = ?
我不熟悉数值傅里叶变换,所以我不知道计算出的概率密度值对应的x值是多少,我也不确定当我用x绘制p值时,我需要用哪个常数(如果有的话)乘以它们。我还没有找到有用的文档。
发布于 2019-11-08 23:04:10
这只是为了让你走。我使用了正变换,因为你的自变量看起来像时间(t)。
import numpy as np
import matplotlib.pyplot as p
%matplotlib inline
t = np.linspace(-5, 5, 100) #presumed to be time
y = np.exp(-(t**2)/2)
p.subplot(121)
p.plot(t,y)
p.subplot(122)
f = np.fft.fftshift(np.fft.fft(y))
freq = np.fft.fftshift(np.fft.fftfreq(100, d=t[1]-t[0]))
p.plot(freq,np.abs(f))
https://stackoverflow.com/questions/58756125
复制相似问题