fft

0热度

1回答

我已经采样了5kHz采样的传感器数据1分钟。因此，一个采样数据文件包含5,000 x 60 = 300,000个数据点。请注意，传感器会测量周期性数据，如60Hz交流电流。现在，我想对一个数据文件应用FFT（使用python numpy.rfft函数）。据我所知，FFT结果的数量是输入数据数量的一半，即在30万个数据点的情况下有150,000个FFT结果。但是，FFT结果的数量太大而无法分

0热度

1回答

将2个相位声码器帧混合在一起

我试图使用相位声码器来冻结声音。我通过存储光谱帧（幅度和相位）以及前一帧和当前帧之间的相位差来实现这一点。要播放冻结帧，我只需将频谱帧重复插入相位声码器的反相功能中，每次使用相位差值递增（并缠绕）相位。这里是我现在正在做的一些伪码（为了简洁），其中frameA和frameB是相位声码器的fft表示的幅度/相位表示。 void analyze(inputSignal) { // conv

1热度

1回答

不同的Python输出

我试图运行Python STFT（也称为短时傅立叶变换）上两个相同采样数据块。但是，由于我的数据的应用程序之一，我需要它的一些被收集到一个.txt文件。我很困惑，为什么STFT将致力于为从.csv进来的数据，而不是生产为从.txt文件进来的数据正确的输出。两者都是Panda.core.series.Series数据类型，并且都具有类似的幅度值。测试条件也相同。，我试图解决的代码是标记代码，C

1热度

1回答

傅立叶变换正常密度函数

我使用下面的MATLAB代码执行正常的密度函数的傅立叶变换： N=100; j=0:(N-1); a=-5; b=5; dx = (b-a)/N; x = a+j*dx; dt = 2*pi/(N*dx); f1 = -N/2*dt; f2 = N/2*dt; t= f1+ j*dt; GX = normpdf(x,0,1); fft_GX = real(fft(GX))';

1热度

1回答

寻找关于FFT模板的帮助

我目前正在研究需要FFT进行卷积的问题，但是当我从我的存档中引入FFT模板时，我发现输出存在问题。例如：输入：（0,0）（0,0）（4166667，0）（1，0）正确的输出：（4166668，0）（-4166667， 1）（4166666，0）（-4166667，-1）模板输出：（4166668，0）（-4166667，-1）（4166666，0）（-4166667，）代码： #defi

0热度

1回答

如何使用FFT实际实现多项式乘法？

我一直在算法教科书中研究这个话题。统一的复杂根源的巧妙使用似乎是数学上的工作。但是，我不明白如何在电脑中实际表示这一点。我能想到的两件事情：使用实/虚分解来表示复数。但是这意味着使用浮点数，这意味着我打开我的算法数值误差，即使我想用整数系数乘以两个多项式，我也会失去精度。将exp（i 2pi/n）表示为n。所以，我最终会得到一个欧米茄中的元组，如果我必须保持这种形式，我基本上会在欧米茄中进

0热度

1回答

在fft matlab中计算相位角

我正试图计算我在Matlab工作区中构建的合成信号的傅里叶谱的第二个分量的频域（在计算fft之后）的相位角。我相信这个阶段等于0（正如你在代码中看到的那样），但是我得到的结果是pi/2。代码如下： t = 0:pi/128:(2 * pi - pi/128); V = sin(t); L = length(V); n = 2^nextpow2(L); Y = fft(V, n); thr

1热度

1回答

数组太大而无法在matlab中划分但不是python

我遇到的问题是Matlab中的数组太大。阵列数据来自音频文件。我想得到冲动回应。我首先FFT了原始和录制的音频。然后由原来的记录划分。最后进行逆FFT以获得脉冲响应。这是我计划要做的事，但我被困在了部门。卡住使用Matlab，我发现python code可以做到这一点很好。我将代码重写到Matlab中，问题又回来了。代码不完整，但足以显示问题。希望得到很多建议和批评。由于计划这样做，但没有

4热度

1回答

matlab中fft，ifft和fftshift的使用

我试图实现split-step fourier方法来求解光学中的非线性薛定谔方程。它基本上分别对待线性部分和非线性部分。它利用傅里叶变换和时域非线性部分解决了线性部分。下面的代码从一本书复制： alpha = 0 beta_2 = 1 gamma = 1 T = linspace(-5,5,2^13); delta_T = T(2)-T(1); L = max(size(A));

0热度

1回答

得到一个像谱图在MATLAB中的阴谋

我想得到一个与spectrogram函数得到的相同类型的阴谋，我试图与contour但我没有得到相同的结果。我写了一个比较这两个情节的小函数。该功能记录音频一秒钟，然后绘制频谱图。积结果的 function graph_comparison a = audiorecorder(44100,16,1); recordblocking(a,1); y = getaudio