DFT到FFT的理解
DFT简化计算理解(FFT)
DFT:
WN=e^(-j*2*pi/N)
DFT复杂度o(N^2)
降低与N^2的依赖 使N = LM (L^2+m^2 <= N^2)
N点DFT分解为M段L点DFT
一维的N点序列变为(L,M)二维序列,每一行分别进行DFT
举例两种一维到二维的映射关系
n = Ml+m
1 | 3 | 5 | 7 | 9 |
2 | 4 | 6 | 8 | 10 |
n = l+mL
1 | 2 | 3 | 4 | 5 |
6 | 7 | 8 | 9 | 10 |
与之所求的DFT 也可存入相对应的(q,p)矩阵中
以第一种(n = Ml+m)为例:k = Mp+q
找书麻烦这里给出推到:
重一维到二维
两种流程:
按列存入信号
计算每行M点DFT
乘以相位因子
计算每一列的L点DFT
按行读取所得数组
图示:
来看下基2_FFT算法:
上图的N/2点的dft可以分解为N/4的而N/4的DFT可以分解为N/8的……直到最后分解为2点的DFT
这儿的2点DFT其实是输出A+B,A-B两个值
为什么可以这样分解呢?其实他就是1个数学式子的’分解‘过程,来看下
N点的DFT是这玩意儿
将序列 奇偶分开
拆开成2km与1两项
X(2m) 可看为f(m)
而这个可看作一个新的N/2点DFT ------------可见N点的DFT已经分解为N/2点DFT
由采样定理,在频域上(N/2以为新的周期)F(k+N/2)=F(k),且, 。
所以得出下式:
复数乘法运算量 。而原始DFtT的量为N^2,当N够大时几乎减小了一半
注意 的确定,他是 分解出来的
这是举例看下8点DFT的奇偶分解
一级dft(抽第二级奇偶) | 二级dft(抽第三极的奇偶) | 第三极dft(最终) |
0 | 0 | 0 |
4 | 2 | 1 |
2 | 4 | 2 |
6 | 6 | 3 |
1 | 1 | 4 |
5 | 3 | 5 |
3 | 5 | 6 |
7 | 7 | 7 |
可用二进制倒序实现 即011100变为001110,感觉镜像啦下
对比下
可看作一个4点DFT,只需求k=0,1的DFT,即可根据关系得出k=2,3
第二级 等于W(4,0),W(1,4)
DFT到FFT的理解的更多相关文章
- FS,FT,DFS,DTFT,DFT,FFT的联系和区别
DCT变换的原理及算法 文库介绍 对于初学数字信号处理(DSP)的人来说,这几种变换是最为头疼的,它们是数字信号处理的理论基础,贯穿整个信号的处理. 学习过<高等数学>和<信号与系统 ...
- 一个蒟蒻对FFT的理解(蒟蒻也能看懂的FFT)
建议同学们先自学一下"复数(虚数)"的性质.运算等知识,不然看这篇文章有很大概率看不懂. 前言 作为一个典型的蒟蒻,别人的博客都看不懂,只好自己写一篇了. 膜拜机房大佬 HY 一. ...
- Matlab 之 FFT的理解和应用
网上看了一些大牛的关于FFT的见解,加上自己的一点儿理解,针对以下这几个问题来加深对FFT的理解. 不知道大家有没有类似以下几点的困惑: 问题的提出 对于1秒钟输出的连续信号,使用采样率Fs不同,就会 ...
- FS,FT,DFS,DTFT,DFT,FFT的联系和区别 数字信号处理
DCT变换的原理及算法 文库介绍 对于初学数字信号处理(DSP)的人来说,这几种变换是最为头疼的,它们是数字信号处理的理论基础,贯穿整个信号的处理. 学习过<高等数学>和<信号与系统 ...
- [信号与系统]傅里叶变换、DFT、FFT分析与理解
目录 一.前言 二.傅里叶变换 1.傅里叶级数 2.傅里叶级数系数求解 2.1.求解方法 2.2.三角函数的正交性 2.3.系数求解过程 2.4.关于傅里叶级数的个人感悟 3.引入复指数 4.总结 三 ...
- FFT算法理解与c语言的实现
完整内容迁移至 http://www.face2ai.com/DIP-2-3-FFT算法理解与c语言的实现/ http://www.tony4ai.com/DIP-2-3-FFT算法理解与c语言的实现 ...
- 【总结】对FFT的理解 / 【洛谷 P3803】 【模板】多项式乘法(FFT)
题目链接 \(\Huge\text{无图,慎入}\) \(FFT\)即快速傅里叶变换,用于加速多项式乘法. 如果暴力做卷积的话就是一个多项式的每个单项式去乘另一个多项式然后加起来,时间复杂度为\(O( ...
- DTFT、DFT、FFT
对于一般的周期信号可以用一系列(有限个或者无穷多了)正弦波的叠加来表示.这些正弦波的频率都是某一个特定频率的倍数如5hz.2*5hz.3*5hz--(其中的 5hz 叫基频).这是傅立叶级数的思想.所 ...
- 【转】离散傅里叶变换-DFT(FFT)基础
转:https://blog.csdn.net/zhangxz259/article/details/81627341 什么是离散傅里叶变换 matlab例子 本文是从最基础的知识开始讲解,力求用最通 ...
随机推荐
- 杂项-EMS:CRM
ylbtech-杂项-EMS:CRM CMS是"Content Management System"的缩写,意为"内容管理系统". 内容管理系统是企业信息化建设 ...
- redis实现消息发布/订阅
redis实现简单的消息发布/订阅模式. 消息订阅者: package org.common.component; import org.slf4j.Logger; import org.slf4j. ...
- IDA Pro 权威指南学习笔记(一) - 启动 IDA
启动 IDA 启动 IDA,有一个欢迎界面 之后有一个对话框 选择 New 将启动一个对话框来选择将要分析的文件 选择 Go 将使 IDA 打开一个空白的工作区 如果要选择分析的文件,可以直接拖到 I ...
- storm集群配置以及java编写拓扑例子
storm集群配置 安装 修改配置文件 使用java编写拓扑 storm集群配置 storm配置相当简单 安装 tar -zxvf apache-storm-1.2.2.tar.gz rm apach ...
- Oracle 10g RAC 如何配置 VIP IPMP
metalink note 283107.1介绍了如何设置VIP的IPMP,此处记录一下设置过程. o Existing 10g RAC installation ^^^^^^^^^^^^^^^^^^ ...
- wpa_supplicant移植与使用(转)
下载wpa_supplicant最新版和openssl(编译wpa_supplicant需要openssl的库) 我这里使用的是wpa_supplicant-0.7.3.tar.gz和openssl- ...
- Linux的基本指令--目录和文件和文件属性和文件用户组
目录和文件 一 . ls:列出目录的内容,未给出目录名或是文件名时,就显示当前目录的信息. -a 列出隐藏文件,文件中以”.”开头的均为隐藏文件,如:~/.bashrc -l 列出文件的详细信息 ...
- 中国大学MOOC 设置AutoCAD的绘图环境
- 回车换行0x0D和0x0A 小谈
在计算机还没有出现之前,有一种叫做电传打字机(Teletype Model 33)的玩意,每秒钟可以打10个字符.但是它有一个问题,就是打完一行换行的时候,要用去0.2秒,正好可以打两个字符.要是在这 ...
- PHP内核介绍及扩展开发指南—Extensions 的编写(下)
第一个参数是HashTable,在1.2.3节提到Zend使用HashTable来存储PHP函数,function_table用于指 定从哪个HashTable中获取函数.通常应该用CG(functi ...