- android - 多次调用 OnPrimaryClipChangedListener
- android - 无法更新 RecyclerView 中的 TextView 字段
- android.database.CursorIndexOutOfBoundsException : Index 0 requested, 光标大小为 0
- android - 使用 AppCompat 时,我们是否需要明确指定其 UI 组件(Spinner、EditText)颜色
我有一个复数数组 (waves
),其排序方式与 fft
返回的方式相同。如果我对该数组调用 ifft
,它会返回原始样本的近似值。
我想在 python 中自己实现 ifft
。我找到了 the formula of IFFT .我实现了它,但它看起来与 ifft
结果有点不同。我试图通过查看 ifft source 来修复它,但这是一个高度优化的版本,我无法找出它是如何工作的
这是我目前所拥有的:
import numpy as np
import matplotlib.pyplot as plt
from scipy.fftpack import fft
# create samples
data = np.zeros(1024)
for k in range(0, 1024):
data[k] = np.sin(k/20 * np.pi*2)
# plot original samples
#plt.plot(data)
#plt.show()
# apply FFT on samples
waves = fft(data)
# plot FFT result
#plt.plot(np.imag(waves),'r')
#plt.plot(np.real(waves),'b')
#plt.ylabel('value')
#plt.xlabel('period')
#plt.show()
#
res = np.zeros(1024)
for k in range(0, 1024):
val = 0.0
for n in range(0, len(waves)-1):
# https://dsp.stackexchange.com/a/510/25943
val += waves[n]*np.exp(-1.j * 2*np.pi * n * k / len(waves)) / len(waves)
res[k] = val.real
#np implementation
res2 = np.fft.ifft(waves)
plt.plot(data, 'b') # original
plt.plot(res,'g') # my implementation
plt.plot(res2,'y') # np implementation
plt.show()
也许必须以不同方式处理零频率项和负频率项。我不确定,因为在任何傅里叶变换的描述中都没有提到
最佳答案
这里只有两个错误:
for n in range(0, len(waves)-1):
应该是
for n in range(0, len(waves)):
因为 range
不包括其上限(与基于 0 的索引一起,与 Matlab 相比,这使得 FFT 类型算法的实现更容易一些)。
此外,
val += waves[n]*np.exp(-1.j * 2*np.pi * n * k / len(waves))
应该是
val += waves[n]*np.exp(1.j * 2*np.pi * n * k / len(waves))
符号约定不同;在 NumPy 中,直接变换有 -1j,逆变换有 1j。
当然,整个过程效率很低,但你可能想自己详细写出来。 NumPy 的矢量化操作将正常使用,以
开头data = np.zeros(1024)
for k in range(0, 1024):
data[k] = np.sin(k/20 * np.pi*2)
被替换为
data = np.sin(np.arange(1024)/20 * np.pi*2)
其他循环也是如此。
关于python - numpy中的一维逆离散傅里叶变换是如何计算的?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/47001893/
我正在尝试从底层叶进程写入命名管道并从顶层进程的管道读取。 为此,我首先在顶级进程中创建 FIFO,然后使用 for 循环派生更多进程。在 for 循环中,我正在检查叶子进程,如果它是叶子,我正在写入
我需要一个XPath表达式来验证XML列表元素是否只有某种类型的节点。 例: 我需要知道列表中是否只有图像。
我想从层次结构表中获取所有叶子及其根。叶子是最低层节点,根是最顶层节点。 给定一棵树: A --B --C --D E --F --G --H 叶子是节点:C, D根是:A、E 该表如下所
尝试从嵌套数组中删除相同的 sibling (叶子)/相同的数组。 例如 $data = [ 'test' => [ 'a' => [
我可以遍历目录并只打印文件夹/目录名称,但我想排除包含其他目录的目录的文件夹名称。出于某种原因,我称其为树结构中的“最终节点”,但我很可能是在自欺欺人,这不是第一次了。 =) 在审查其他答案列表时,也
我现在经常使用 Folium,在 Python 中拥有如此易于使用的东西真的很棒。但是他们的文档严重落后,我理解。所以我有2个问题。 我在尝试获得更多标记颜色时做错了什么?这是我尝试过的:map.si
我正在尝试使用 ape 包中的 plot.phylo 命令为 R 中系统发育类型图的边缘(线)添加颜色。这个例子是一个“粉丝”类型的图,虽然我希望这个方法与“系统图类型”或其他任何东西相同。 libr
我正在尝试呈现一个 JTree 叶以具有包含可编辑 JTextArea 的 JPanel,以便用户可以输入文本并能够在添加文本时向下滚动, 并且还使用包装器来限制文本的宽度,使其向下延伸以留出更多空间
这是要修改的正确 Linux 内核代码吗?我如何进行更改以模拟 CPUID 代码以及我需要更改哪个函数。谢谢 #include #include #include #include #incl
我是一名优秀的程序员,十分优秀!