OpenCV中的逆傅立叶变换
逆傅立叶变换（Inverse Discrete Fourier Transform，IDFT）是傅立叶变换（DFT）的逆过程，用于将频域信号转换回时域信号。在OpenCV中，逆傅立叶变换可以通过cv2.dft()函数实现。
在OpenCV中，cv2.dft()函数可以实现快速傅立叶变换（DFT）和逆傅立叶变换（IDFT）。其中，DFT将时域信号转换为频域信号，而IDFT将频域信号转换回时域信号。
逆傅立叶变换的公式如下：
X(u) = 1/N * sum(x(n) * e^(-2*pi*i*n*u/N))其中，X(u)表示频域信号，x(n)表示时域信号，N表示信号长度，u表示频率，i表示虚数单位，e表示自然对数的底。
在OpenCV中，可以使用以下代码实现逆傅立叶变换：
import cv2
import numpy as np
# 加载图像
img = cv2.imread('image.jpg', cv2.IMREAD_GRAYSCALE)
# 将图像填充为偶数尺寸
rows, cols = img.shape
if rows % 2 != 0:
img = cv2.copyMakeBorder(img, 0, 1, 0, 0, cv2.BORDER_REPLICATE)
if cols % 2 != 0:
img = cv2.copyMakeBorder(img, 0, 0, 0, 1, cv2.BORDER_REPLICATE)
# 快速傅立叶变换
dft = cv2.dft(np.float32(img), flags=cv2.DFT_COMPLEX_OUTPUT)
# 计算逆傅立叶变换
idft = cv2.dft(dft, flags=cv2.DFT_INVERSE | cv2.DFT_SCALE, nonzeroRows=rows // 2)
# 显示原始图像和逆傅立叶变换后的图像
cv2.imshow('Original Image', img)
cv2.imshow('IDFT Image', np.abs(idft))
cv2.waitKey(0)
cv2.destroyAllWindows()在上述代码中，首先加载图像并将其转换为灰度图像。然后，将图像填充为偶数尺寸，以便进行快速傅立叶变换。接下来，使用cv2.dft()函数进行快速傅立叶变换，并计算逆傅立叶变换。最后，显示原始图像和逆傅立叶变换后的图像。
总之，逆傅立叶变换是一种将频域信号转换回时域信号的方法，在OpenCV中可以使用cv2.dft()函数实现。
相关·内容
1回答
我想开始一些有趣的项目，并发现我对声音感兴趣。所以我想让一些东西随着声音的节奏移动。但我不知道，例如，winamp (在图像中)如何生成值来绘制这些条形图。我怎样才能在C#中实现呢？有人能给我一些建议吗？
编辑:我会弄清楚的：
浏览 1提问于2012-05-07得票数 1
回答已采纳
我需要改变一个音频文件的声音调制，尝试了一整天的SoundPool，但它只改变了音频的播放速度。我想改变调制方式。需要您的建议，提前谢谢。
浏览 0提问于2012-10-22得票数 0
回答已采纳
1回答
快速2D卷积实现？
、、、、
我已经做了一个用于2D卷积的CUDA程序，现在我想将它与一些非CUDA实现进行比较，以测量加速比。
我可以使用经典的多循环方法或matlab的conv2与我自己的纯C实现进行比较，但这不像是一个合法/公平的比较，因为它们不是最快的实现。
此外，我正在考虑尝试OpenCV，我一直在寻找一个SIMD优化版本，但没有运气。有什么建议吗，我应该用OpenCV吗？
注意:我已经阅读了其他问题，包括，但答案基本上与我的纯C代码或对各种可用方法的讨论相同。
浏览 1提问于2011-06-03得票数 3
回答已采纳
2回答
我读过一些关于如何使用信号的fft更有效地计算自相关的解释，将实数部分乘以复共轭(傅立叶域)，然后使用逆fft，但我在Matlab中实现这一点有困难，因为在细节层面上。
浏览 99提问于2010-10-16得票数 23
回答已采纳
我要逆狄拉克三角函数的傅里叶变换，δ(t)。所以我做了这些：
>> syms F w t
>> F(w)=dirac(w)
F(w) = (symfun) DiracDelta(w)
>> ifourier(F(w))
ans = (sym)
1
----
2*pi
虽然正确的结果应该是这样的：
为什么计算错误？
浏览 3提问于2020-06-12得票数 0
回答已采纳
1回答
opencv检测光频
、、、、
我目前正在做一个机器人汽车的爱好项目。在相机前面，我将放置两个光源，并使用计算机视觉检测光源，并将我的机器人汽车放在它前面。
在Raspberry Pi 3上使用Python和OpenCV，我做了以下工作:我使用了2个红色光源。使用以下OpenCV函数:cvtColor(帧到HSV)，inRange(红色滤镜)，findContours(查找2个光源)，Moments(获取光源的x，y坐标)，计算光源与帧边缘的距离，并生成汽车必须移动的方向。
现在的问题是，我无法确定哪个红色光源是左光源，哪个是右光源。例如，如果我用两个红色光源转动我的工作台，程序将不能告诉汽车必须在光源周围移动。
我有两个想
浏览 0提问于2016-05-11得票数 0
2回答
使用FFT实现2D卷积
、、、
对于使用大核(滤镜)卷积大图像，TensorFlow.conv2d()的速度非常慢。将1024x1024图像与相同大小的内核进行卷积需要几分钟时间。为了进行比较，cv2.filter2D()会立即返回结果。
我找到tf.fft2()和tf.rfft()了。
然而，我并不清楚如何使用这些功能执行简单的图像过滤。
如何使用快速傅立叶变换在TensorFlow中实现快速2D图像滤波？
浏览 2提问于2018-05-22得票数 3
1回答
我有两张大的灰度图像。PIL.Image或numpy数据结构。
如何沿一个轴对两个图像进行一维卷积？
我能想到的最好的是
def conv2(im1, im2, *args):
res = 0
for l1, l2 in zip(im1, im2):
res += np.convolve(l1, l2, *args)
return res
这是可行的，但速度不是很快。有没有更快的方法？
请注意，所有2D卷积函数可能都不相关，因为我对2D卷积不感兴趣。我以前在上面看到过这个问题，但我没有看到比我的代码更好的答案。所以我又把它撞上了。
浏览 12提问于2016-08-03得票数 0
我的目标是在DFT中获取OpenCV中的图像。
使用dft函数，我可以计算它，然后通过计算它的大小来绘制它(然后，应用日志，最后将它规范化，以便绘制0到1之间的值)。
我的结果是，对于下面的图像，我向您展示的结果(通过交换以使图像中心有较低的频率)：
然而，如果我将其与使用其他工具(如 )所获得的结果进行比较，我的结果似乎不正确，因为它似乎具有“很高”的值(我指的是OpenCV DFT值)：
我想可能是因为这些原因：
DFT (at OpenCV)与FFT (Halcon)的区别
我正在执行的操作是为了在OpenCV中显示震级。
第一个的问题是我很难分析，
浏览 1提问于2015-06-01得票数 4
回答已采纳
1回答
我在GTX690上使用opencv 。
我在内部检查了该函数，发现gpu::matchTemplate()没有在进程中使用任何FFT，而它的CPU对等项使用了。
我的问题是，在GPU中不使用FFT背后的原因是什么？
浏览 4提问于2014-01-10得票数 3
这是我第一次尝试在C#上分析音频流。在没有任何经验的情况下，感谢您分享您的解决方案。:)
使用NAudio(用于处理音频的C#开源)，我尝试将两个音频流与以下阶段进行比较：
在.wav文件中读取并以16位格式存储在float[]中:音频=新AudioFileReader("test.wav")；_buffer =新floatwave.length；audio.Read (_buffer，0，_buffer.Length)；
将float[]传递给SmbPitchShifter()，并设置PitchShift=1(不改变音高)，将osamp设置为跃点大小，并传递_buff
浏览 0提问于2018-04-17得票数 0
回答已采纳
我正在尝试实现一个过滤器，它看起来像
在Matlab中。我有：
omegas = (-length(t)/2:length(t)/2)*2*pi/tau/10;
SOS_freq_sum = zeros(1,length(omegas));
for i = 1:length(K_set)
k = K_set(i);
SOS_freq_sum = SOS_freq_sum + b_k(i)*sinc(omegas/(2*pi/tau)-k);
end
SOS_filter_in_frequency = (tau/sqrt(2*pi))*((SOS_freq_sum)'
浏览 5提问于2014-12-04得票数 0
1回答
有没有实现傅里叶逆变换的Java/Android库？我发现有几个实现了傅里叶变换，但没有一个实现了它的逆。
浏览 7提问于2011-09-14得票数 2
回答已采纳
1回答
java中的同态过滤
、、、、
有谁知道同态过滤的概念和工作原理，以及它的java实现。我需要比较两张图片。从摄像机到数据库。我使用的是朴素的相似性算法。它工作得很好。但是换光是个问题。我们如何解决照明问题？我被困在测试上了。android 2.2
浏览 0提问于2013-02-26得票数 1
3回答
我试图用matlab证明白噪声具有恒定的功率谱密度，但频谱的振幅看起来是随机振幅。
有人能告诉我为什么吗？
这是我的代码。
noise = randn(1,10000);
fft_noise=fft(noise);
plot(abs(fft_noise(1:5000)))
谢谢。
浏览 0提问于2012-04-23得票数 1
回答已采纳
1回答
从.m4a数据中过滤频率
、、、、
我试图分析一个声音文件(.m4a)，以获得随着时间推移的振幅，并绘制一个图表。我在网上发现了一些很好的代码(下面)。不过，我想另外过滤掉所有的声音，不是在一个目标的频率范围。我只想绘制1900年到2100赫兹之间的声音。我该怎么做？
var processingBuffer = [Float](repeating: 0.0, count: Int(readFile.arrayFloatValues.count))
let sampleCount = vDSP_Length(readFile.arrayFloatValues.count)
vDSP_vabs(readF
浏览 3提问于2017-12-03得票数 2
我正在用OpenCv用C++编写一个程序，它需要找到两个图像之间的矩阵，并找到存在于响应阈值之上的所有峰值。我确实找到了一个名为phaseCorrRes()的函数，但它使用minMaxLoc()只返回最高的峰值。如何修改此函数以检索所有峰值？
浏览 3提问于2013-11-25得票数 1
2回答
如何获得图像的低、中、高频分量？
、、、、
我目前正在寻找一些图像处理项目，只是想知道如何获得图像的低频，中和高频成分？例如，正如这张图片所显示的(我是从googling获得的，没有详细描述如何获得这张图片，但可能使用了一些过滤)。
另外，我遇到了使用离散余弦变换(DCT)的帖子，它可以帮助我们获得图像的低频和高频分量。只是想知道如何利用DCT来获得中频分量？
我对过滤也有非常基本的知识。我认为也有高斯高通/低通滤波器可用。以及基于小波的滤波。只是想知道高斯滤波、小波滤波和DCT滤波有什么区别？我应该用哪一种？
浏览 3提问于2015-06-17得票数 2
2回答
有没有办法在图像中找到莫尔条纹，我可以用Swift或者OpenCV在我的iOS应用中使用？
任何帮助都将不胜感激。
浏览 5提问于2018-08-30得票数 0
1回答
我需要比较两个音频信号，signal1和signal2。Signal1是一种白噪声。Signal2是与signal1相同的信号，但有一种特殊的均衡，用于切断或衰减某些频率。如何获取两个音频信号在频域中的比值？(例如:在100 to的频率下，与signal1相比，signal2衰减50% )。我需要这些信息来处理第三个信号，方法是应用与signal2中转换signal1相同的均衡。
我使用来处理我的数据，并将其从时域传递到频率域。这段代码对于signal1和signal2是相同的。
DoubleFFT_1D fft1 = new DoubleFFT_1D (FFT_SIZE);
double[]
浏览 0提问于2013-05-15得票数 1
回答已采纳
1回答
我应该如何在Lazarus中为Mac创建一个音调生成器？
我似乎在网上找不到任何关于这个的东西？
没有Wav文件相关的东西或TProcess相关的东西。
在这里，事情需要实时运行。
基本上，我希望能够播放SIN波并通过扬声器输出它们。
浏览 2提问于2014-03-15得票数 2
2回答
我有两个8个无符号字节的序列，我需要计算它们的循环卷积，这将产生8个无符号19位整数。当我重复这一百万次的时候，我想要优化。
简单的方法需要64次MAC操作。我已经知道如何通过SSE/AVX指令来加速，这不是我想要的。
有没有其他方法，可能是基于FFT或数论变换来减少运算次数或其他技术来获得一些加速？
实际上，我不需要8个值:最大的一个和相应的移位就足够了。
浏览 1提问于2016-06-25得票数 2
1回答
我想知道是否有人遇到过这个问题。
我可以从matlab中重建出与原始图像相似的图像，但是，实际值总是不同的。
例如，原始图像在矩阵中的值从0到1，而我的重建图像的范围从-0.2到0.4。
重建的图像看起来与原始图像相似，只是图像中的数据具有不同的尺度。
这是我所说的示例代码。
p=phantom(64);
theta=0:1:179;
r=radon(p,theta);
ir=iradon(r,theta);
figure
subplot(1,2,1);imagesc(p)
subplot(1,2,2);imagesc(ir)
浏览 3提问于2015-06-08得票数 2
1回答
我看到，FFTW计划需要几秒钟时间，所以最好只调用一次，然后在for循环中调用execute()。这是我的案子。但是当然，在每个循环中，输入数据是不同的，但是大小是一样的，那么我如何改进它呢？用布尔变量更好地解决这个问题吗？
fftw_plan my_plan;
in = (fftw_complex*) fftw_malloc(sizeof(fftw_complex)*N);
v = (double*) fftw_malloc(sizeof(double)*N);
out = (fftw_complex*) fftw_malloc(sizeof(fftw_complex)*N);
my_
浏览 6提问于2013-12-03得票数 0
回答已采纳
2回答
fft理解
、、
有没有人能很好地解释一下快速傅立叶变换图像变换是如何分析快速傅立叶变换图像和它的Re^2+Im^2图像的？我只是想在看图像的时候理解一些东西，以及它的频率。
浏览 2提问于2010-03-31得票数 7
回答已采纳
2回答
我有一些FFT数据，257维，每10毫秒，121帧，即1.21秒。我认为第一个维度可能是其他的，剩下的是FFT系数，我猜。它可能只是一些光谱数据。从关于快速傅立叶变换数据的评论来看，sqrt10和均值-方差-归一化可能已经应用于该数据。
从那里，我想计算一些44.1 Hz的PCM信号，这样我就可以播放声音了。我问过，但也许StackOverflow是一个更好的地方，因为我实际上想实现它。关于理论，我也问了同样的问题。
我该怎么做呢？也许我需要更多的信息(我必须以某种方式找出)--哪一个？也许这些缺失的信息可以通过某种方式智能地猜测出来？
这个问题既是一个理论问题，也是一个实际实施问题。我想这个
浏览 0提问于2015-08-16得票数 1
4回答
我已经读了很多关于快速傅立叶变换的书，并试图理解它的低级方面。不幸的是，谷歌和维基百科根本帮不上什么忙。我打开了5本不同的算法书，但也没有多大帮助。
我在试着找一些简单的东西的FFT，比如向量1，0，0。当然，我可以直接把它插入到Matlab中，但这并不能帮助我理解下面发生的事情。另外，当我说我想要找到一个向量的FFT时，这是否等同于我想要用一个更有效的算法来找到一个向量的DFT？
浏览 0提问于2010-07-12得票数 15
2回答
我正在努力用Android来模糊位图。
我看过很多关于使用简单内核的信息，比如
0 0 0 5 0 0 0
0 5 18 32 18 5 0
0 18 64 100 64 18 0
5 32 100 100 100 32 5
0 18 64 100 64 18 0
0 5 18 32 18 5 0
0 0 0 5 0 0 0
我的问题是，我真的不确定如何将它与我的位图以一种有效的方式相乘。
我是不
浏览 0提问于2010-09-09得票数 5
回答已采纳
1回答
OpenCV中的傅里叶反变换
、、、、
我是OpenCV和图像处理算法方面的新手。我需要在C++的OpenCV中做离散傅立叶逆变换，但我不知道怎么做。我在网上找了一下，也没找到答案。我正在用这个页面上的代码在我的程序中做傅立叶变换：。我尝试过与该代码相反的操作，但我不知道我在哪里做错了。我的代码如下(我认为整个代码都是错误的)：
void doFourierInverse(const Mat &src, Mat &dst) {
normalize(src, dst, 0, -1, CV_MINMAX);
int cx = dst.cols/2;
int cy = dst.rows/2;
Mat
浏览 0提问于2012-04-22得票数 4
回答已采纳
1回答
fft和小波
、、、
我可以使用快速傅立叶变换来获得加载的1秒音频文件的频率、相位和幅度，然后重新创建它。我现在要做的是找出每个频率在1秒音频文件中的开始和结束位置，并将数据放入一个数组中
例如:100 at从.23seconds到.34seconds，104.34 at从.35seconds开始，到.37seconds结束。
fft可以做到这一点吗?或者我需要将我的整个程序转换为使用小波？另外，在octave中有没有一些小波例子来展示如何完成我想要完成的任务？
我使用的是repo中的Ubuntu Linux 12.04和Octave 3.2.4
谢谢，瑞克
浏览 9提问于2012-08-23得票数 2
回答已采纳
1回答
我需要一个好的音调变化的解决方案为我的项目，以改变声音。很多的音调偏移js库都试过了，但是它们都没有提供想要的结果。主要的事情是不能控制结果的声音音色，我得到了米老鼠或地狱僵尸声音的东西，但不是真正的声音。然而在这里，如果用vega的声音测试，结果是非常出色的：不幸的是，我对音频处理一无所知，至少想知道它是如何完成的，这里使用了什么类型的移位算法节奏，以及我们如何实现对过程的音色/共振峰控制。如有任何提示，非常感谢。谢谢;)
浏览 7提问于2018-11-25得票数 0
我正在使用Tacotron2为一种新的语言训练TensorflowTTS模型。我成功地训练了模型(执行了预处理、规范化和解码生成的少数输出文件)--输出目录中的文件是.npy文件。这是有意义的，因为他们是梅尔-平面图。我正试图找到一种将上述文件转换为.wav文件的方法，以检查我的工作是否已经完成。
我用这个：
melspectrogram = librosa.feature.melspectrogram(
"/content/prediction/tacotron2-0/paol_wavpaol_8-norm-feats.npy", sr=22050,
wi
浏览 14提问于2022-02-18得票数 0
1回答
我在使用gnuplot和python来拟合这条曲线时遇到了问题。
我想要适合的文件是。我想用一个简单的正弦函数来拟合它，如下所示：
a+b*sin(c*x+d)
我在gnuplot和python的拟合方法中都有这样的结果：
从x-scale的微小差异中，python中的一个被重新缩放，两个数字都给了我一个不正确的拟合。
浏览 3提问于2015-10-26得票数 1
1回答
我使用MATLAB对图像应用了不同标准差的高斯低通滤波器，并记录了每种方法所需的时间。我看到在频域中实现滤波器要高效得多(更快)。有人对此有什么解释吗？
浏览 0提问于2011-06-23得票数 3
回答已采纳
4回答
我已经在Matlab中对时间序列波形执行了fft (快速傅立叶变换)，但我似乎在傅立叶变换图中实际上有一个奇怪的波形，尽管有尖峰，这个波形看起来只在时域中才能看到。有没有编程上的原因导致这种情况发生？
浏览 0提问于2009-12-17得票数 1
我正在尝试使用C++接口编写一个使用FFTW3.3和OpenCV 2.1的FFT/IFFT函数。我见过很多使用旧的OpenCV格式的例子，我做了一个直接转换，但有些东西不起作用。
我的函数的目标是返回一个带有FFT的实部和虚部的Mat对象，就像dft的默认OpenCV函数一样。下面是该函数的代码。在将im_data复制到data_in的代码行中，程序因内存问题而被阻塞。
有人知道我做错了什么吗？谢谢
Mat fft_sr(Mat& I)
{
double *im_data;
double *realP_data;
double *
浏览 0提问于2012-11-26得票数 1
回答已采纳
我一直在看OpenCV华丽的源代码来理解分数的计算。您可能已经注意到，该库可以在图像中找到具有置信度分数的模板。cvMatchTemplate函数给出了模板在图像中的概率和位置。
由于DFT很慢，我决定使用FFTW库进行相位相关。我得到了图像，将其扩展到2的幂的大小，然后对模板进行相同的处理，首先计算两个数组的FFT，将两个数组相乘，计算逆FFT，最后对结果进行归一化。
如果我在位图上显示这个结果，它会清楚地显示模板的位置，但是，我很难计算相关性分数。我试着在调试器下逐行运行cvMatchTemplate方法，但没有什么帮助，因为这个库使用了离散傅立叶变换和任意的图像大小。该方法中使用的所有表
浏览 0提问于2011-01-09得票数 4
回答已采纳
如何使用opencv (python)、或pil进行图像拼接？或者任何其他python库，我有一个看起来像的图像，它是从200个图像创建的，我需要按正确的顺序将它们拼接起来形成图像。
浏览 1提问于2012-03-30得票数 0
回答已采纳
1回答
有没有办法在这些图像上检测半圆形？
尝试了很多方法，结果都很糟糕
示例：
image = cv2.imread('path_here', 0)
output = image.copy()
ret, thresh = cv2.threshold(image, 0, 255, cv2.THRESH_OTSU)
circles = cv2.HoughCircles(thresh, cv2.HOUGH_GRADIENT, 1, 110,
param1=25, param2=5, minRadius=9, maxRadius=10)
dete
浏览 2提问于2021-09-12得票数 2
1回答
尝试用小波变换对歌曲进行频域表示，用离散小波变换进行分解，并用Matlab将演唱者的频率放在原处
问题是，dwt和分解mades仅在时域中表示它。
如果DWT不能表示它，我怎么能用频率表示它呢？
谢谢
浏览 1提问于2013-05-15得票数 0
快速傅立叶变换(fft；)将'a‘转换为它的傅立叶，频谱等效：
numpy.fft.fft(a, n=None, axis=-1, norm=None)
据我所知，参数n表示输出中有多少个样本，如果n小于a中的样本数，则会裁剪输出；如果n较大，则用零填充输出。
axis是做什么的？这到底是什么意思？我还没有找到任何明确的例子来说明它的用法。
浏览 0提问于2020-05-14得票数 0
3回答
如何生成包含wav/riff声音文件频谱的文件？我想使用linux命令行。
我知道生成png频谱图的sox函数很酷
sox sound.wav -n spectrogram
但我不需要频率谱的直观表示。我只想在数据文件中获得频率的频谱，这样我就可以处理它们了。我相信肯定有使用sox的选择。Sox需要在绘制数据之前生成数据。如何获得此功能？
不确定，也许第二种解决方案是将wav文件导出到dat文件中。来自dat文件的每个样本是膜在某一时刻的位置的测量。所以这不是频率谱。
sox sound.wav file.dat
如何将这些膜位置转换成我需要的光谱？
谢谢。
浏览 8提问于2014-02-13得票数 11
2回答
傅里叶逆变换FFT3W
、、、
我使用C++函数来求傅立叶逆变换。
int inYSize = 170; int inXSize = 2280;
float* outData = new float[inYSize*inXSize];
fftwf_plan mReverse = fftwf_plan_dft_c2r_2d(inYSize, inXSize,(fftwf_complex*)temp, outdata,
FFTW_ESTIMATE);
fftwf_execute(mReverse);
我的输入是具有复数的二维数组temp。所有元素的实值为1，虚值为0。
所以我期望这样的数组的InverseFFT应该是具有实
浏览 6提问于2013-04-25得票数 1
1回答
我需要在窦上使用FFT函数。
这是我的窦功能：
octave:1> t=(0:1/16000:1);
octave:2> A=2;
octave:3> x=A*sin(2*pi*80*t);
这是我尝试做fft的方法：
octave:60> l=length(x);
octave:61> k=fft(x,l);
它让我笑成三角形..我怎样才能正确地做到这一点？
@edit screen，sin的fft好看吗？我能以某种方式改进它吗？
@edit2我可以改进它，让它看起来更好吗？
浏览 2提问于2013-01-10得票数 1
回答已采纳
1回答
我想对图像进行傅里叶变换。但是如何将图片更改为数组呢？在此之后，我认为我应该使用numpy.fft.rfft2进行转换。以及如何从数组转换回图像？提前谢谢。
浏览 0提问于2010-04-16得票数 3
回答已采纳
我目前有自己的FFT代码，原型是
void doFFT(short* data, int nSamples);
参数data的大小为2*nSamples。输出在相同的变量中返回。
我想用libffmpeg中的代码替换我的FFT代码。我在上看过快速傅立叶变换的测试样本。然而，我并不清楚如何正确地使用它。我想知道是否有人可以分享一些伪代码或为我指出正确的方向。
浏览 1提问于2016-01-03得票数 1
3回答
我在试着处理音频数据。我正在使用Java。我已经将音频数据提取到一个数组中。现在，我应该将N个数据样本传递给一个函数，该函数计算离散傅立叶变换(或更有效的快速傅立叶变换)。我已经阅读了文档，但我越来越困惑。我要计算的是星等谱(|X(k)|)。有谁可以帮我？谢谢
浏览 1提问于2009-06-07得票数 13
回答已采纳
所以我写了这个：
function y = convolution(u,v)
[m,n] = size(u);
[w,z] = size(v);
y = zeros(m-w+1,n-z+1);
for i = 1:m-w+1
for j = 1:n-z+1
y(i,j) = sum(sum(u(i:i+w-1,j:j+z-1).*v));
end
end
end
然后，我将其与matlab中的conv2()进行了比较，结果如下：
function timedConv
a = im2double(rgb2gray(imread(&#
浏览 3提问于2016-11-08得票数 2
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
