2

我正在尝试使用 scipy.fftp 在 python 中平滑和插入一些周期性数据。我已经设法获取数据的 fft，删除高于 wn 的高阶频率（通过执行 myfft[wn:-wn] = 0），然后使用 ifft(myfft) 重建数据的“平滑”版本。ifft 创建的数组与原始数据的点数相同。我如何使用该 fft 创建一个具有更多点的数组。

x = [i*2*np.pi/360 for i in range(0,360,30)]
data = np.sin(x)
#get fft
myfft = fftp.fft(data)
#kill feqs above wn
myfft[wn:-wn] = 0
#make new series
newdata = fftp.ifft(myfft)

我还能够以与此处演示的分辨率相同的分辨率手动重新创建 系列使用 FFT 结果重新创建时间序列数据而不使用 ifft

但是当我尝试提高 x 值数组的分辨率时，它也没有给我正确的答案。

提前致谢

尼尔

4

1 回答 1

5

返回的是np.fft.fft位置 0 处的 DC 分量，然后是所有正频率，然后是 Nyquist 频率（仅当元素数量为偶数时），然后是相反顺序的负频率。因此，为了增加分辨率，您可以在奈奎斯特频率的两侧添加零：

import numpy as np
import matplotlib.pyplot as plt

y = np.sin(np.linspace(0, 2*np.pi, 32, endpoint=False))

f = np.fft.fft(y)
n = len(f)
f_ = np.concatenate((f[0:(n+1)//2],
                     np.zeros(n//2),
                     [] if n%2 != 0 else f[(n+1)//2:(n+3)//2],
                     np.zeros(n//2),
                     f[(n+3)//2:]))
y_ = np.fft.ifft(f_)
plt.plot(y, 'ro')
plt.plot(y_, 'bo')
plt.show()

在此处输入图像描述

于 2013-06-07T20:34:34.790 回答