上一篇随笔，简要写了一下FFT中数组重新排序的算法。现在把完整的FFT代码分享给大家(有比较详细的注释)。/*2015年11月10日于河北工业大学*/#include #include #include #include const int N=8;//数组的长度const double PI=3....
分类：
其他好文 时间：
2015-11-10 22:27:42
阅读次数：
302
最近一直在看傅里叶变换，看到FFT算法，其实算法的关键之一，蝶形运算，只要看懂了，编码实现并不难。反倒是其中位码倒序的环节，看很容易看懂，但是编码实现不是那么容易的。在网上参考了很多资料后，决定把下面这个算法分享给大家，在这里要感谢百度文库用户letsgotoyy123提供的《快速傅里叶变换FF.....
分类：
编程语言 时间：
2015-11-10 14:05:30
阅读次数：
275
网上看了一些大牛的关于FFT的见解，加上自己的一点儿理解，针对以下这几个问题来加深对FFT的理解。 不知道大家有没有类似以下几点的困惑：问题的提出对于1秒钟输出的连续信号，使用采样率Fs不同，就会得到不同的采样序列。在计算N点FFT之后得到N点的复数，幅频图中横坐标如何与我实际的频率相对应？ ...
分类：
其他好文 时间：
2015-11-05 21:58:10
阅读次数：
229
最近项目压力比较大，为了赶时间很多代码都得图简便，然而碰到的问题还是需要重新整理一下，即便当时不懂事后也得弄清楚。项目的主要任务是一个C6678的PCI板卡驱动，用于FFT计算，一个图形界面显示程序显示处理前后结果。设备操作上，需要实时从C6678的内存中读取两个数据，一个是64KB的unsig.....
分类：
其他好文 时间：
2015-11-04 02:05:21
阅读次数：
345
题目就是求一个大数的乘法这里数字的位数有50000的长度，按平时的乘法方式计算，每一位相乘是要n^2的复杂度的，这肯定不行我们可以将每一位分解后作为系数，如153 = 1*x^2 + 5*x^1 + 3*x^0 (在这里x可以理解成10)那么两个数字相乘就相当于系数相乘后得到新的系数组合如153 *...
分类：
其他好文 时间：
2015-10-26 13:41:53
阅读次数：
243
1. FFTW介绍FFTW由麻省理工学院计算机科学实验室超级计算技术组开发的一套离散傅立叶变换(DFT)的计算库，开源、高效和标准C语言编写的代码使其得到了非常广泛的应用，Intel的数学库和Scilib(类似于matlab的科学计算软件)都使用FFTW做FFT计算。FFTW是计算离散Fourier...
分类：
其他好文 时间：
2015-10-21 19:08:18
阅读次数：
319
图像的频率是表征图像中灰度变化剧烈程度的指标，是灰度在平面空间上的梯度。如：大面积的沙漠在图像中是一片灰度变化缓慢的区域，对应的频率值很低；而对 于地表属性变换剧烈的边缘区域在图像中是一片灰度变化剧烈的区域，对应的频率值较高。傅立叶变换在实际中有非常明显的物理意义，设f是一个能量有限的模拟信号，则....
分类：
其他好文 时间：
2015-10-19 23:51:47
阅读次数：
303
%fft and pwelch方法求取功率谱load x.matFs = 1; t = (0:1/Fs:1-1/Fs).'; Nx = length(x); % Window data w = hanning(Nx); xw = x.*w; % Calculate power nfft = ...
分类：
其他好文 时间：
2015-10-19 17:09:58
阅读次数：
159
题目连接：http://acm.hdu.edu.cn/showproblem.php?pid=1402hdu_1402:A * B Problem PlusTime Limit: 2000/1000 MS (Java/Others)Memory Limit: 65536/32768 K (Java/...
分类：
其他好文 时间：
2015-10-11 10:12:11
阅读次数：
212
NTT:快速数论变化，对于FFT精度减少的情况，NTT可以避免但是会慢一点，毕竟是数论有Mod,和快速米 引用：http://blog.csdn.net/zz_1215/article/details/40430041 周边介绍。利用原根，在ZP整数域（后悔没学好《信息安全数学基础》 ...
分类：
其他好文 时间：
2015-09-18 02:00:36
阅读次数：
223
