标签[fft] - 堆栈内存溢出

[英]Why does my own implementation of Fast Fourier Transform in N dimensions give different results than NumPy's?

我为 N 维快速傅立叶变换编写了以下代码，但它没有给我与 numpy 的 function 相同的结果。我在 1d 中计算 FFT 的代码如下：它的基本原理很难解释，它基于位反转，但我已经检查它是否正常工作，所以问题出在 N 维代码上。 ...

使用 python（类似 audacity）从音频文件中获取最大频谱

[英]Get maximum of spectrum from audio file with python (audacity-like)

我是数字信号处理的新手，我试图找到音频文件频谱的峰值，我通常使用 Audacity 和 plot 频谱打开音频文件。通过可视化上面的频谱，我可以找到 120HZ 的峰值，但这需要一些手动工作。我想用 Python 以更编程的方式找到峰值。我不确定 Audacity 中绘制了哪个频谱，但我假设它 ...

归一化复数空间中的向量：mean=0 and std-dev=1/sqrt(n)

[英]Normalize vectors in the complex space : mean=0 and std-dev=1/sqrt(n)

我有一个向量或一堆向量（按行存储在二维数组中）向量生成如下：平均值 = 0，标准偏差 = 1/SQRT（vec_len）操作之前或之后必须以相同的形式标准化我想在复杂的空间中对它们进行归一化。这是 function 的包装器：帮我填点。这是我使用的实值归一化。我不想发布这个，以免影响 ...

FFT 是如何工作的？ Rust 中每一帧的频谱似乎都看不到

[英]How does FFT work? Can't seem to the frequency spectrum for each frame in Rust

我每 0.03 秒只能获得 735 个频率的音频数据。我想弄清楚如何获得每一帧中的频谱。下面的代码每帧仅返回 735 个数据（因为samples_in_frame是 735，这是每帧中的样本数），但我想要整个 ~20,000hz 用于 0.03 秒的样本。我将如何 go 这样做？ path ...

计算多维数组中每个数组的每一列的傅里叶变换

[英]Computing the Fourier-transform of each column for each array in a multidimensional array

在下面的 4 维数组中，每一列代表机器学习4D开发的一个属性。import numpy as np from scipy.fft import fft, fftfreq, fftshift A = np.array([ [[[0, 1, 2, 3], ...

在一维矢量上使用一维 FFT 的二维 FFT

[英]2D FFT using 1D FFT on 1D vector

基于这篇文章，我可以在每一行上应用 1DFFT，然后在二维向量的每一列上应用。如果我有一个一维向量，我可以使用像这样的单元格将其视为二维向量v[rowIndex * columnCount + columnIndex] 我的 FFT1D 算法将 0 填充到矢量，直到下一个 2 的幂 positi ...

在 flutter 中检测扫描文档方向

[英]Detect scanned document orientation in flutter

有没有办法检测 flutter 中的扫描文档方向，python 中有示例根据文本方向检测图像方向角度和如何检测图像方向（文本），但有什么方法可以在 flutter 中实现这一点。 ...

在频域中过滤信号振幅最高的两个频率

[英]Filtering the two frequencies with highest amplitudes of a signal in the frequency domain

我试图过滤振幅最高的两个频率。我想知道结果是否正确，因为过滤后的信号似乎不如原始信号平滑？ FFT 函数的 output 包含基频 A0/C0 是否正确，将其包含在最大幅度的搜索中是否正确（它确实是最高的？）？我的代码（基于我的教授和同事的代码，到目前为止我还不了解每个细节）：我不太了解 F ...

Python / cvxpy 中快速矩阵向量乘法的优化问题

[英]Optimization Problem with fast matrix-vector multiplication in Python / cvxpy

我想解决以下（凸）最小化问题： min ||x||_1 在约束 sgn(A[x,R]=y) 和 ||x||_2 = 1 下其中A是 mx(N+1) 矩阵， x in R^N 是向量，而\[x,R\]是通过附加给定数字R创建的向量。目标是找到x的最优值。 A是傅里叶矩阵，有快速矩阵向量、求逆等算 ...

如何获得PSD的幅度

[英]How to get the amplitude of PSD

我分享了从此处提取的 Python 代码：在 Python 中绘制快速傅里叶变换。我正在制作傅里叶谱，然后我想“选择”以从图上的特定点获取值。在下一张图片中，您可以看到 50 Hz 的频率，但我想知道下一个峰值的频率值。我怎么能找到那个？ ...

如何使用 FFT 对振幅数据获得每分钟的节拍数？

[英]How can I get the beats per minute using FFT on the amplitude data?

我检查了关于这个主题的所有关于 stackoverflow 的答案，但任何解决方案都不适用于我的数据集。该数据集包含以每秒 60 个样本采样的 600 个振幅。在这个数据集中，10 秒内有 13 次心跳。所以结果应该是每分钟 78 次左右。我知道我可以通过多种方式使用峰值检测，但我想知道是 ...

时间序列的自举方法

[英]The bootstrapping method with timeseries

我收到有关我的代码的错误消息：Error in statistic(data, original, ...) : unused argument (original) 我需要运行引导程序方法来估计周期图中的 pvalue。我的时间序列是由包含 365 天数据的ts()函数创建的。首先，我定义了 ...

如何使用 np.fft.ifft2 对图像进行反向傅立叶变换 - Python

[英]How to reverse Fourier Transform on an image with np.fft.ifft2 - Python

我是 Python 的傅里叶变换新手。由于傅立叶变换，我想在图像上隔离一个场。这是我的照片：这是我应该获得的：到目前为止，这是我的代码：我的问题是我获得了这张图片而不是我想要的图片（见上文）：谁能帮帮我吗？ ...

使用 numpy fft 在频域中求导并转换回时间的缩放问题

[英]Scaling issue taking derivative in frequency domain and transforming back to time using numpy fft

我正在尝试使用 numpy 在频域中对信号进行导数。对于连续信号，这只是频谱乘以 2*pi*i*f。对于离散信号，考虑到 numpy FFT 实现，我认为它会是这样的：即 FFT 乘以 2*pi*i*k/n 但事实并非如此：不需要除以 n。这在下面的 plot 中显示：比较绿色曲线（除以 n ...

不能 plot matlab 二维快速傅立叶变换的结果

[英]Can't plot results of matlab 2D fast fourier transform

我使用 matlab fft2 命令（这个命令： https://www.mathworks.com/help/matlab/ref/fft2.html ）对 4436413x3 矩阵进行了二维傅立叶变换。乘积是一个包含复数值的 4436413x3 矩阵。当我尝试 plot 数据时，我最终得到一 ...

对任意频域信号进行IFFT

[英]Taking IFFT of Arbitrary Frequency Domain Signal

我习惯于从时域移动到频域，然后通过以下方式回到时域：对实值时域信号（即没有虚部的信号）进行 fft 从步骤 1 中获取频域信号的 ifft 忽略步骤 2 中 ifft 的虚部。但是，我想跳过第 1 步，而是在频域中定义一个任意复数（即实部和虚部）信号，然后使用 ifft 获得时域等效信号。要 ...

将 Mathematica 傅里叶级数代码转换为 Python

[英]Converting Mathematica Fourier series code to Python

我有一些简单的 Mathematica 代码，我正在努力将其转换为 Python 并且可以使用一些帮助： plot 看起来像这样：对于上下文，我有一个 function f(t)：我需要 plot 前 10 项的总和。在 Mathematica 中，这非常简单，但出于某种原因，我似乎无 ...

如何使用 np.fft.fft function 进行傅里叶变换？

[英]How can I take the Fourier Transform using np.fft.fft function?

我想生成两个正弦信号，其中 =5 Hz 和 f2=3 Hz，信号的持续时间为 =1。 x1(t)=sin(2) X2(t)=sin(2(1+f2)) X()=1()+2()因此，我将尝试分析该信号的频域表示。如何使用 np.fft.fft function 通过使用 DFT 大小 =64、采样 ...

计算 Julia 中的离散傅立叶变换矩阵？

[英]Compute a discrete Fourier Transform matrix in Julia?

在 Julia 中，我想随机生成一个大小为 n x n 的离散傅立叶变换矩阵。我目前不确定该怎么做。有没有人知道在 Julia 中执行此操作的方法？ ...

二维高斯 kernel 的傅里叶变换与空间域中的对应项不匹配

[英]Fourier transform of 2D Gaussian kernel is not matching up with its counterpart in the spatial domain

我们知道高斯滤波器的傅里叶变换在频域中也是高斯的，我编写了以下方法来构建高斯 kernel：def get_gaussian(size, sigma): g_kernel = np.zeros((size,size)) x_center = size // 2 y_center = si ...