1和0组成的数组如何进行傅立叶变换
傅立叶变换（Fourier Transform）是一种将信号从时域转换到频域的数学方法。它可以将一个由连续时间或离散时间的信号组成的数组表示为一系列正弦和余弦波的加权和。
对于一个由1和0组成的数组进行傅立叶变换，首先需要将这个离散的时域信号进行零填充（zero-padding），使其长度变为一个较大的2的幂次方，以获得更精确的频域表示。
然后，可以使用离散傅立叶变换（Discrete Fourier Transform，DFT）或快速傅立叶变换（Fast Fourier Transform，FFT）算法来进行变换。DFT是一种精确计算傅立叶变换的方法，但计算复杂度较高，而FFT是一种基于DFT的高效计算方法，广泛应用于信号处理和频谱分析。
在傅立叶变换的结果中，频域表示的数组的每个元素代表了原始时域信号中对应频率的幅度和相位信息。通常通过计算每个频率对应的幅度谱（magnitude spectrum）或功率谱（power spectrum）来分析信号的频率成分。
傅立叶变换在许多领域中都有广泛的应用，包括信号处理、图像处理、音频处理、视频压缩、通信系统等。在云计算领域中，傅立叶变换常用于音频、视频和图像处理，例如音频频谱分析、图像滤波、视频编码等。
对于腾讯云相关产品，腾讯云提供了一系列与云计算相关的产品和服务，可以支持傅立叶变换的应用和计算需求。例如：
- 腾讯云云服务器（Cloud Virtual Machine，CVM）：提供了高性能的虚拟机实例，可以用于进行傅立叶变换的计算任务。
- 产品介绍链接：https://cloud.tencent.com/product/cvm
- 腾讯云函数计算（Serverless Cloud Function，SCF）：无需管理服务器的计算服务，可用于执行傅立叶变换等短时、低并发的计算任务。
- 产品介绍链接：https://cloud.tencent.com/product/scf
- 腾讯云弹性MapReduce（Elastic MapReduce，EMR）：提供了大数据分析和处理的服务，可以支持大规模数据集上的傅立叶变换计算。
- 产品介绍链接：https://cloud.tencent.com/product/emr
请注意，以上仅为腾讯云提供的部分相关产品，具体选择和使用哪种产品需要根据实际需求和场景进行评估和决策。
相关·内容
我正在尝试扩展fft代码，它可以很好地用于python中的一维图像数组。实际上，我知道问题出在扩展的逻辑上。我对FFT了解不多，我必须提交图像处理作业。我将对任何提示或解决方案表示感谢。
这是代码，实际上，我正在尝试用python为FFT创建一个模块，在rosetta code网站的帮助下，它已经在1D上工作得很好了。
from cmath import exp, pi
from math import log, ceil
def fft(f):
N = len(f)
if N <= 1: return f
even = fft(f[0::2])
odd
浏览 0提问于2013-03-03得票数 1
回答已采纳
2回答
我有一个n浮点流的列表，每个流都有不同的大小。
可以使用以下规则将流组合在一起：
您可以将流放在任何时间点(它在开始之前为零)。您可以多次使用同一个流(它可以自己重叠，甚至几次处于同一位置)，并且允许您根本不使用某个流。
例如：
输入流：
1 2 3 4
2 4 5 6 7
1 5 6
可以像这样组成：
1 2 3 4
1 5 6
1 5 6
在放置之后，输出流由这样的规则组成，即每个输出浮点数等于每个项的平方和的平方根。
例如：
如果某个位置的流是：
1
2
3
输出为：
sqrt(1*1 + 2*2 + 3*3) = sqrt(14) = 3.74...
因此，对于示例
浏览 1提问于2011-09-27得票数 12
1回答
我感兴趣的信号是风速矢量和时间矢量，它们形成了10分钟的风速剖面。我如何提取组成完整十分钟信号的各个波的信息(例如，10个波，平均速度为9m/s，频率为1 1hz等等)，所以我最终得到了所有组成波及其频率、最大/最小值等的列表。我尝试将几秒的“测试”波与10分钟向量进行互相关比较，但似乎不起作用
浏览 19提问于2020-02-05得票数 0
我需要将具有整数值的数组傅里叶变换到频域(以便稍后将其与另一数组相乘)。输出数组的大小必须为44100，但输入数组会有所不同。
我认为fftw3是一个很好的工具。但是我如何为输入和输出创建一个具有不同数组大小的“计划”呢？
下面是我写的函数：
fftw_complex* fourier( int* samples, int numberOfSamples ){
fftw_complex* input;
fftw_complex* output;
fftw_plan plan;
input = (fftw_complex*)fftw_malloc( sizeof
浏览 0提问于2016-06-30得票数 1
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
3回答
我正在使用FMOD库从MP3中提取PCM。我得到了整个2声道16比特的东西，我还得到了44100 16的采样率是1秒内44100个“声音”样本。我不明白的是，16位的值到底代表什么。我知道如何在xy轴上绘制坐标，但我要绘制的是什么？Y轴代表时间，x轴代表什么？音量？这和振幅一样吗？如何确定组成此值的不同声音。我的意思是，我如何从一个16位的数字中得到一个谱。
这可能是一个单独的问题，但它实际上是我真正需要回答的问题:如何获得每25毫秒的振幅？我要取44，0.025个值，除以40 (40 *100秒=1秒)吗？这给出了1102.5个样本；那么我会把1102值输入到黑盒中，给出那个时刻的振幅吗？
浏览 1提问于2012-05-01得票数 3
回答已采纳
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
回答已采纳
3回答
利用fft进行频率提取
、、、
我有一个信号，它是两个频率的线性组合，我正在尝试使用MATLAB来确定(主)信号由哪些频率组成。作为输入，我有采样率rate和带有信号数据segment的行向量。
现在我有了以下代码：
N=length(segment);
freq = rate*(0:N/2)/N;
X = fft(segment);
X=X(1:N/2+1);
plot(freq*2, abs(X))
这确实返回了一个具有两个峰值的图，每个峰值位于信号组成的每个频率处。现在我想提取这些频率。所以我想返回一个变量，signal，它是两个频率的向量。
因此，如果我的信号是由30 Hz的信号和60 Hz的信号的线性组合创建的
浏览 20提问于2019-11-08得票数 0
1回答
opencv检测光频
、、、、
我目前正在做一个机器人汽车的爱好项目。在相机前面，我将放置两个光源，并使用计算机视觉检测光源，并将我的机器人汽车放在它前面。
在Raspberry Pi 3上使用Python和OpenCV，我做了以下工作:我使用了2个红色光源。使用以下OpenCV函数:cvtColor(帧到HSV)，inRange(红色滤镜)，findContours(查找2个光源)，Moments(获取光源的x，y坐标)，计算光源与帧边缘的距离，并生成汽车必须移动的方向。
现在的问题是，我无法确定哪个红色光源是左光源，哪个是右光源。例如，如果我用两个红色光源转动我的工作台，程序将不能告诉汽车必须在光源周围移动。
我有两个想
浏览 0提问于2016-05-11得票数 0
2回答
我想在C++程序中用一个由16位整数组成的2-D数组尽可能高效地执行可变块大小的计算。我感兴趣的是一个实时块匹配代码。我想知道是否有可用的软件库来做这件事？代码在windows XP上运行，我不得不使用Visual Studio2010进行编译。中央处理器是一个双核AMD Athlon 64 x2 4850e。
对于可变块大小的绝对差和(SAD)计算，我指的是以下内容。
我有一个较小的二维数组，我将称之为template_grid，另一个较大的二维数组，我将称为image。我希望在图像中找到使模板中的像素和图像中区域中的像素之间的绝对差值之和最小的区域。
在C++ if中计算SAD的最简单方法
浏览 1提问于2011-11-12得票数 8
2回答
我有一些FFT数据，257维，每10毫秒，121帧，即1.21秒。我认为第一个维度可能是其他的，剩下的是FFT系数，我猜。它可能只是一些光谱数据。从关于快速傅立叶变换数据的评论来看，sqrt10和均值-方差-归一化可能已经应用于该数据。
从那里，我想计算一些44.1 Hz的PCM信号，这样我就可以播放声音了。我问过，但也许StackOverflow是一个更好的地方，因为我实际上想实现它。关于理论，我也问了同样的问题。
我该怎么做呢？也许我需要更多的信息(我必须以某种方式找出)--哪一个？也许这些缺失的信息可以通过某种方式智能地猜测出来？
这个问题既是一个理论问题，也是一个实际实施问题。我想这个
浏览 0提问于2015-08-16得票数 1
我正在尝试实现一个过滤器，它看起来像
在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
2回答
我有两个8个无符号字节的序列，我需要计算它们的循环卷积，这将产生8个无符号19位整数。当我重复这一百万次的时候，我想要优化。
简单的方法需要64次MAC操作。我已经知道如何通过SSE/AVX指令来加速，这不是我想要的。
有没有其他方法，可能是基于FFT或数论变换来减少运算次数或其他技术来获得一些加速？
实际上，我不需要8个值:最大的一个和相应的移位就足够了。
浏览 1提问于2016-06-25得票数 2
4回答
我已经读了很多关于快速傅立叶变换的书，并试图理解它的低级方面。不幸的是，谷歌和维基百科根本帮不上什么忙。我打开了5本不同的算法书，但也没有多大帮助。
我在试着找一些简单的东西的FFT，比如向量1，0，0。当然，我可以直接把它插入到Matlab中，但这并不能帮助我理解下面发生的事情。另外，当我说我想要找到一个向量的FFT时，这是否等同于我想要用一个更有效的算法来找到一个向量的DFT？
浏览 0提问于2010-07-12得票数 15
我想用python在一个浮动表示中量化声音的深度。例如，低沉的声音可能由低于0的分数表示，而高音(与深度相反)可能由高于0的分数表示。低沉的声音可以与阿诺德·施瓦辛格的声音相提并论，而高音则可以与小孩子的声音相提并论。
我还没有从谷歌搜索中找到任何明确的方法来解决这个问题。我已经尝试使用来获得音色分数，但似乎计算没有用。也就是说，低音和高音的得分没有太大的区别。
基于这个，我可以使用快速傅立叶变换分析。然而，我不明白如何将FFT图简化为语音深度的单一表示。
那么，在使用FFT获得分析后，我该如何进行呢？在更广泛的问题中，我如何将声音的深度表示为浮动表示，特别是使用Python？
编辑:为了获得
浏览 3提问于2019-12-02得票数 1
如果我在嘈杂的背景下录制音频，然后将录制的音频文件与本地保存的音频文件进行比较，以检查两者是不是相同的歌曲？在iPhone中
Shazam和Sound Hound应用程序如何在iPhone中生成录制音频的指纹？
任何人都可以探索在ios中从音频文件生成指纹的算法知识。
谢谢
浏览 2提问于2013-06-20得票数 0
回答已采纳
所以我写了这个：
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
回答已采纳
1回答
我正在尝试实现一个函数，它取一个wav文件，通过AForge的快速傅立叶变换运行100秒的音频。当我改变偏移量来改变我通过FFT在音频中计算的位置时，有时我会得到可以在我的图形中显示的结果，但大多数时候我会得到一个复杂的NaN数组。为什么会这样？
这是我的代码。
public double[] test()
{
OpenFileDialog file = new OpenFileDialog();
file.ShowDialog();
WaveFileReader reader = new WaveFileReader(file
浏览 2提问于2014-11-25得票数 0
1回答
我有2组x，y坐标，我试图确定它们的相关系数。一种方法是分别计算x值和y值的相关系数。那么用这些个体相关系数的算术平均值来估计x，y对的相关系数是否公平呢？
浏览 2提问于2015-07-09得票数 3
1回答
我使用Python中的statsmodels库对我的时间序列进行了季节性分解。如何获得组成观测信号的不同子信号的周期/频率？
有给出季节信号的傅里叶分解的函数吗？
浏览 1提问于2015-12-29得票数 2
1回答
尝试从加速度计数据中获取主频
、、、、
我正在开发一个应用程序，它将返回使用加速度计的手机正在经历的振动频率。
这是我个人用的..只是为了获得DSP和在我的手机中使用传感器的经验。
因此，我抓取了一些传感器数据，并将绝对值写入一个数组。我的下一步是使用快速傅立叶变换来尝试获得主要频率。这就是我迷路的地方。
实际上，我的最终目标只是确定数据是周期性的还是非实时的。所以，如果我说的都错了，请给我指出正确的方向。
如果我的方向是正确的，那么我在代码中哪里出错了呢？我是android和java的新手，可能我也犯了一些严重的逻辑错误。如果我可以获取一个数组，对其执行FFT，实时获得主频率，并显示结果……我在认真地听呢。
下面的代码片段：(数组
浏览 0提问于2015-02-28得票数 0
2回答
使用FFT实现2D卷积
、、、
对于使用大核(滤镜)卷积大图像，TensorFlow.conv2d()的速度非常慢。将1024x1024图像与相同大小的内核进行卷积需要几分钟时间。为了进行比较，cv2.filter2D()会立即返回结果。
我找到tf.fft2()和tf.rfft()了。
然而，我并不清楚如何使用这些功能执行简单的图像过滤。
如何使用快速傅立叶变换在TensorFlow中实现快速2D图像滤波？
浏览 2提问于2018-05-22得票数 3
我是Python的新手，但是有一个我没有准备好的家庭作业。请帮帮忙。我写道： import numpy as np
X,Y = np.loadtxt("wavedata.csv", usecols=(0, 1), delimiter=",", unpack=True) 以及为a_0、a_n、b_n和梯形规则定义的函数。但是CSV文件只有0到10的值。傅立叶系数不需要从-L积分到L吗？我如何在这里使用梯形规则，或者生成计算系数的函数？ 有什么建议吗？你可以和我分享程序的任何部分吗？非常感谢。
浏览 19提问于2019-10-09得票数 0
1回答
下面的代码只显示了一个图形，但是我需要声音的频率。我正在尝试录制语音并获得实时频率，这样我就可以弹奏钢琴或吉他的声音并找到频率。
public class AudioProcessing extends Activity implements OnClickListener {
int frequency = 8000;
int channelConfiguration = AudioFormat.CHANNEL_CONFIGURATION_MONO;
int audioEncoding = AudioFormat.ENCODING_PCM_16BIT;
private RealDoubl
浏览 0提问于2012-02-17得票数 3
我正在寻找一些关于如何在MATLAB中压缩时间序列数据的建议。
我有一些瞳孔大小的数据集，这些数据是在1秒内收集的，每次试验有25,000个点(我仍然不确定将这些数据称为“时间序列”是否合适)。从现在开始，我想要做的是将它们与另一个数据进行比较，我需要将点数压缩到大约10,000或更少，以最小化信息的损失。有什么方法可以做到这一点吗？
我试着搜索如何做到这一点，但我所能找到的是平滑数据或压缩数字图像的方法，这些方法已经完成或对我没有用处。
·数据集只由瞳孔直径组成，随着时间的推移而变化。对于每个试验，在1秒内收集了25,000个点的数据，这意味着1个点表示测量了0.04毫秒的瞳孔直径。我想要做
浏览 2提问于2011-11-19得票数 2
快速傅立叶变换(fft；)将'a‘转换为它的傅立叶，频谱等效：
numpy.fft.fft(a, n=None, axis=-1, norm=None)
据我所知，参数n表示输出中有多少个样本，如果n小于a中的样本数，则会裁剪输出；如果n较大，则用零填充输出。
axis是做什么的？这到底是什么意思？我还没有找到任何明确的例子来说明它的用法。
浏览 0提问于2020-05-14得票数 0
我目前有自己的FFT代码，原型是
void doFFT(short* data, int nSamples);
参数data的大小为2*nSamples。输出在相同的变量中返回。
我想用libffmpeg中的代码替换我的FFT代码。我在上看过快速傅立叶变换的测试样本。然而，我并不清楚如何正确地使用它。我想知道是否有人可以分享一些伪代码或为我指出正确的方向。
浏览 1提问于2016-01-03得票数 1
5回答
我被问到一个面试问题，我需要在哪里使用它，但我不知道它是什么。
那么，在简单的英语中，什么是，我如何使用它来求函数的导数，给定它的(x，y)值作为输入？
您将如何实现它？
编辑：
我之所以这样问，是因为给定一系列(x，y)值，我需要计算函数的外观，推导它并找到它不断变化的次数(即，(0，1)，(1，2)被算作1)或根本不变(0，5)，(1，5)也被算作一次变化)。
浏览 3提问于2011-02-14得票数 3
回答已采纳
1回答
我正在尝试使用NAN在2D数组上应用FFT。我不想用任何东西替换NaN值。有没有一种方法可以应用FFT而忽略NaN值？
arr=np.array([[1, 2, 3], [4, NAN, 6], [1, 2, NAN]])
f=np.fft.fft2(arr)
浏览 4提问于2018-11-23得票数 0
1回答
原始数据在上。它是两列数据，t和x。我做了以下离散fft变换。我不太明白主峰(尖峰)的高度比侧峰低，第二次子图显示，尖峰(最接近2.0)确实是主频。代码和数字如下：
import numpy as np
import math
import matplotlib.pyplot as plt
from scipy.fftpack import fft,fftfreq
freqUnit=0.012/(2*np.pi)
data = np.loadtxt(fname='data.txt')
t = data[:,0]
x = data[:,1]
n=len(t)
d=t[1]-t[
浏览 1提问于2015-11-18得票数 1
回答已采纳
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
回答已采纳
在matlab中，我一直使用函数"fft(x)“，其中"x”是复数的向量。我正在寻找一个易于使用的函数在C++中，将返回复数。
我找到了这个代码：
如果它是等效的，我如何使用它？感谢您的宝贵时间！
/*
This computes an in-place complex-to-complex FFT
x and y are the real and imaginary arrays of 2^m points.
dir = 1 gives forward transform
dir = -1 gives reverse transform
*/
浏览 0提问于2016-05-05得票数 0
三维空间(x，y，t)中的时间序列(X，Y，T)满足：
x(t) = f1(t), y(t) = f2(t),
where t = 1, 2, 3,....
换言之，坐标(x，y)随时间戳t而变化。计算x(t)或y(t)的FFT很容易，但如何计算(x，y)的FFT呢？我假设它不应该被计算为2D-FFT，因为这是针对图像的，而(x，y)只是一个序列。有什么建议吗？谢谢。
浏览 19提问于2020-03-25得票数 0
作为我们开发的一部分，我们使用Python 3.8 (使用NumPy 1.19)开发算法。然后我们使用neon NE10库在C中实现它。在一个阶段中，我们使用内置的IFFT函数计算傅立叶逆变换。当我们比较neon (NE10) (最新版本) IFFT和Python的NumPy 1.19 IFFT时，我们得到了不同的结果(它们差别不大，但差别远远大于数值精度，请参见下面的输出)。
Python代码：
import numpy as np
NFFT = 4096
FREQUENCY_BINS = 2049
NUMBER_OF_OUTPUTS = 4
in02_IFFT = np.zeros((
浏览 5提问于2020-11-05得票数 1
1回答
请随时指出现有代码中的任何错误/改进 所以这是一个非常基本的问题，我对信号处理只有一个初级的理解。我有一个1.02秒的加速度计数据，采样频率为32000 Hz。在python中执行FFT后，我希望提取以下频域特征- 平均频率，中值频率，功率谱变形，谱能量，谱峰度，谱偏度，谱熵，RMSF (均方根频率)，RVF (根方差频率)，功率倒谱。 更具体地说，我正在寻找这些功能的绘图作为最终输出。 包含数据的csv文件有四列:时间、X轴值、Y轴值、Z轴值(加速度计为三轴加速度计)。到目前为止，在python上，我已经能够可视化时域数据，对其应用卷积滤波器，应用快速傅立叶变换，并生成一个显示了一个有趣的
浏览 98提问于2019-04-11得票数 0
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
回答已采纳
2回答
fft理解
、、
有没有人能很好地解释一下快速傅立叶变换图像变换是如何分析快速傅立叶变换图像和它的Re^2+Im^2图像的？我只是想在看图像的时候理解一些东西，以及它的频率。
浏览 2提问于2010-03-31得票数 7
回答已采纳
1回答
我使用我的应用程序捕获声音。假设这个声音是一个正弦的1 KHz声音，并且存在一个背景声音。如何识别声音中是否存在此1 KHz声音？
我的意思是，我可以想象如何在图像中找到元素，例如，如果你在图像上寻找一个黄色正方形，你所要做的就是指定你想要的颜色，给它一定的容差，然后寻找一组具有这种颜色并形成正方形的像素。但是声音呢？如何识别波形和频率，当您得到的只是一个幅值，代表了1/44.000的波形在一秒？
我不想要代码，因为这对这篇文章来说太复杂了，但如果你们能给我指出如何做到这一点的正确方向，可以例证其背后的技术或数学的免费源代码，我将不胜感激。谢谢
浏览 3提问于2011-08-28得票数 0
回答已采纳
1回答
我有一个可以工作的FFT，但我的问题是如何将它转换为IFFT？我被告知IFFT应该就像你正在使用的FFT一样。那么如何从fft I c#生成ifft呢？有人告诉我应该只做几个改动就可以得到ifft。
我试着自己去做，但我没有得到我输入的相同的值。
所以我做了一个值数组，并把它放到fft中，然后ifft和我无法得到我输入的相同的值…
所以我不认为我以正确的方式改变了它。
这是我的FFT：
public Complex[] FFT(Complex[] x )
{
int N2 = x.Length;
Complex[] X = new Complex[N2
浏览 12提问于2014-07-22得票数 0
1回答
如何在FFT程序中表示一个使用声波(正弦波，1000 in，3秒，-3dBFS，44.1 FFT)？程序的输入是振幅和采样率的列表。
我的意思是如何将声音文件(例如: XYZ.wav文件)转换为快速傅立叶变换的输入，其中一个输入参数需要获取由振幅和其他输入参数组成的.dat文件，还需要获取采样率，如果需要的话。
浏览 2提问于2018-10-18得票数 1
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
1回答
我有一些带有重复模式的图像。我想根据相似的模式找到相似的图像。
这些图案是由十字架、三角形、正方形组成的，它们组合在一起，形成了由这些“原始形状”组成的更复杂的结构。例如，想象一个由三角形或六边形等组成的十字架。
这些装饰品出现在壁纸、地毯、短裙、羊毛衫甚至一些油画中。
维基百科的一个例子是
我的典型应用是，比方说，找到其中有相似图案的毛衣。暂时忽略颜色。
我尝试提取SIFT描述符(使用C++和OpenCV)，并在两个图像之间进行匹配。但是，它们匹配很小的区域，例如六边形和三角形的顶点，但理想情况下，我希望匹配三角形和矩形的实际形状。
如果我缩小图像，效果会更好一些，但我仍然觉得我需要一种不同
浏览 0提问于2018-08-22得票数 2
1回答
我已经编写了一个简单且非常小的Matlab代码，用于计算给定数组(或向量)的离散傅立叶变换。
我手动算出了答案，我的Matlab代码也给出了相同的答案。但fft通过掉期指数给出了一个与此不同的答案。下面是我做过的人工计算：
这是第二张图片：
这是第三张图片：
根据这些计算，很明显，我的答案是{12, -3-3j, -2, -3+3j}
下面是我使用的Matlab代码：
clc;
clear all;
close all;
inp=[1,2,3,4];
j=sqrt(-1);
op=zeros(1,length(inp));
for k=1:length(inp)
浏览 0提问于2017-07-27得票数 3
回答已采纳
我得到了连续傅立叶变换的基本概念，我对离散傅立叶变换有点模糊，对快速傅立叶变换更模糊。我想知道如何在maxima中计算和解释某些函数的快速傅立叶变换(FFT)。例如，我如何计算sin(x)的FFT？我知道如何在MATLAB中执行此操作，但我不确定如何在maxima中执行此操作。
浏览 1提问于2012-12-15得票数 0
1回答
这是我关于二维DFT的问题。我知道如何做二维离散傅立叶变换、卷积和滤波。然而，我仍然有一个基本的问题。也许我还没有理解这个概念，这就是为什么我不知道在这里做什么。
因此，基本上，当我们想通过将原始图像的傅里叶变换乘以高斯滤波器来过滤场景的空间频率内容时。如何删除每度0.5、1和2个周期以上的空间频率内容(cpd；每个图像12、24、49个周期的低通截止)。
我们还知道，对于3 dB的标准衰减，高斯滤波器的on是空间频率截止的函数。
非常提前感谢您。
浏览 1提问于2015-10-15得票数 0
此代码从时域计算DFT。有人能看到下面的代码并帮助我得到正确的答案吗？我的问题是:当我将N值更改为4、5、10或其他值时。X(1)会随之改变。但我认为X(1)对于N的每个值都必须是相同的，就像下面的形状:N值改变，但垂直值是相同的。如果你能帮我，我将不胜感激。谢谢。
clear; clc;
% %% Analytical
N=4;
k=0:N-1;
X=zeros(N,1);
t=k/N;
x=(5+2*cos(2*pi*t-pi/2)+3*cos(4*pi*t))
%x=abs((1-(0.012.*(pi.*52.*(t-0.3721)).^2)).*exp(-(pi.*52.*(t-
浏览 8提问于2020-11-10得票数 0
1回答
我有一个.wav文件，它是我在弹吉他音符的时候录制的。然后我使用下面的程序来读取我的.wav文件数据。我使用的是Naudio库。
AudioFileReader readertest = new AudioFileReader(@"E:\song\music.wav");
int bytesnumber = (int)readertest.Length;
var buffer = new float[bytesnumber];
readertest.Read(buffer, 0, bytesnumber);
for (int i = 0; i < buffer.Leng
浏览 0提问于2015-10-21得票数 1
2回答
我希望有人能审阅我下面的代码，并提供一些提示，如何加速tic和toc之间的部分。下面的函数试图比Matlab的内置函数更快地执行IFFT，因为(1)几乎所有的fft系数bin都是零(即，从10M到300M的10到1000 bin都是非零的)，并且(2)只保留IFFT结果的中心三分之一(第一和最后三分之一被丢弃-因此一开始就不需要计算它们)。
输入变量为：
fftcoef = complex fft-coef 1D array (10 to 1000 pts long)
bins = index of fft coefficients corresponding to fftcoef (10 t
浏览 5提问于2011-01-25得票数 2
回答已采纳
我有一个很大的数据集(大约17,000个值)，它由一个datetime对象数组和另一个关联的float值数组组成。我正在尝试将一个粗略的余弦函数映射到所有数据点，以便帮助分析睡眠模式。我一直在尝试找出将余弦曲线映射到所有这些数据的最好方法。我研究了一些方法，发现了一些类似的情况，但我遇到的一个挑战是如何在使用datetimes而不是其他任意值的情况下进行计算。 到目前为止，我编写的代码将必要的数据提取到两个数组中。一个是浮点数，另一个是关联的日期时间值。 示例数据： 下面是我正在处理的数据的示例，它是从调试器复制过来的。第一个数组包含日期时间，第二个数组包含整数值。我一直在寻找一些方法来拟合
浏览 8提问于2019-02-26得票数 0
我在Matlab中找到了一些代码：
y = fft(y,nfft);
其中y是512×443二维阵列，并且nfft =512。
我认为快速傅立叶变换是针对一维数组的，对于二维数组应该有fft2，但是快速傅立叶变换是有效的。这怎么可能，有人能给我解释一下吗？
浏览 4提问于2011-12-08得票数 2
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
即时通信 IM
实时音视频活动推荐
腾讯云开发者
