一般的fft需要将点数补成2的整数次幂，MATLAB中有fft函数，输入N点序列则输出也是N点序列，其中N不一定为2的整数次幂。所以会疑惑MATLAB做的fft和N点序列对应的dft是否一样。经过验证，MATLAB中对N点序列做fft的结果与N点dft的结果是完全一样的，没有对N点序列进行补0后再做 ...
分类：
其他好文 时间：
2016-10-19 20:23:38
阅读次数：
235
实数DFT，复数DFT，FFTFFT是计算DFT的快速算法，但是它是基于复数的，所以计算实数DFT的时候需要将其转换为复数的格式，下图展示了实数DFT和虚数DFT的情况，实数DFT将时域中N点信号转换成2个(N/2+1)点的频域信号，其中1个(N/2+1)点的信号称之为实部，另一个(N/2+1)点的 ...
分类：
其他好文 时间：
2016-10-19 02:18:44
阅读次数：
129
计算2000次1000点的fft，分别使用CPU和GPU,得到的结果差4个数量级，是我程序优化的问题，不都是说GPU做计算更快一些，我得到的恰好相反，请高人指导，谢谢。 the running time of cpu is : 0.018837the running time of gpu is : ...
分类：
其他好文 时间：
2016-10-08 02:26:17
阅读次数：
426
i=0 f= 0.0Hz RealAmp=-0.0 ImagAmp=+0.0ii=1 f= 10000.0Hz RealAmp=-0.0 ImagAmp=+0.0ii=2 f= 20000.0Hz RealAmp=-0.0 ImagAmp=+0.0ii=3 f= 30000.0Hz RealAmp= ...
分类：
其他好文 时间：
2016-10-07 23:28:13
阅读次数：
314
题目链接 题意：给定两个长度为n的数组a与长度为m的数组b, 给定一个操作符op满足 x op y = x < y ? x+y : x-y. 有q个询问，每次给出询问c，问：有多少对(i, j)满足a[i] op b[j] = c ？ 0 <= c <= 100000, 其余数据范围在[0, 500 ...
分类：
其他好文 时间：
2016-10-07 01:19:14
阅读次数：
126
题意：dp[n] = ∑ ( dp[n-i]*a[i] )+a[n], ( 1 <= i < n) cdq分治。 计算出dp[l ~ mid]后，dp[l ~ mid]与a[1 ~ r-l]做卷积运算。 1 #include <bits/stdc++.h> 2 using namespace std ...
分类：
其他好文 时间：
2016-10-06 00:11:47
阅读次数：
186
花了一两个小时研究了一番powergui的fft。在powergui中，FFT Analysis Tool如图所示 那么如果我想提取出Simulink模型中的数据并且做fft处理，得到自己想要的分析图像，应该怎么办呢？ 方法一： 其实MATLAB在这个过程中调用了power_fftscope()这个 ...
分类：
其他好文 时间：
2016-10-05 21:47:03
阅读次数：
873
在数字信号处理中常常需要用到离散傅立叶变换(DFT)，以获取信号的频域特征。尽管传统的DFT算法能够获取信号频域特征，但是算法计算量大，耗时长，不利于计算机实时对信号进行处理。因此至DFT被发现以来，在很长的一段时间内都不能被应用到实际的工程项目中，直到一种快速的离散傅立叶计算方法——FFT，被发现 ...
分类：
编程语言 时间：
2016-10-04 11:42:52
阅读次数：
280
FFT是离散傅立叶变换的快速算法，可以将一个信号变换到频域。有些信号在时域上是很难看出什么特征的，但是如果变换到频域之后，就很容易看出特征了。这就是很多信号分析采用FFT变换的原因。另外，FFT可以将一个信号的频谱提取出来，这在频谱分析方面也是经常用的。 虽然很多人都知道FFT是什么，可以用来做什么 ...
分类：
其他好文 时间：
2016-10-04 09:55:33
阅读次数：
279
