- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
西安电子科技大学课程论文 数字图像处理
快速抽取傅里叶变换
班级:070821
作者:董文凯
学号 时间:2011-06-30
快速抽取傅里叶变换
实验要求:
编写基2时间抽取快速傅里叶变换算法程序,用快速傅里叶变换计算下列两个序列。其中的线性卷积,比较当时快速傅里叶变换方法与用卷积定义法求卷积的计算速度。
实验内容:
1.基2时间抽取FFT的推导过程:
1.1基本原理:
设序列点数为,若不满足,认为加0使之达到要求。设旋转因子为,快速抽取傅里叶变换将的序列按的奇偶分为两组:,其中。则有:
式中,和分别是和的的DFT,
另外,的取值范围是:
又利用
可得
考虑到
及前半部分,
可得后半部分
1.2复杂度分析:
对任一次迭代,共有对结点,因此共需次复数加法和次复数乘法;因此次迭代的总计算量为:复数加法次数,复数乘法次数为;算法复杂度为。
2.算法流程:
2.1 蝶形运算信号流图:
蝶形运算符号见图1,以2点、4点和8点为例表示出蝶形信号流程图如图2、图3和图4所示。
图1 蝶形运算符号
图2 两点蝶形信号流图
图3四点蝶形信号流图
图4八点蝶形信号流图
2.2 算法流程图:
图5 基2时间抽取FFT算法流程图
其中倒序部分算法流程见图6。
图6 倒序部分算法流程
实验结果及分析:
将序列和基2时间抽取FFT进行变换得
当N=4 时计算结果分别为:
X1 =
0.5839 - 0.4697i -0.9665 - 0.0000i 0.5839 + 0.4697i 6.1342 + 0.0000i
X2 =
1.4142 - 0.7321i -1.3178 - 0.0000i 1.4142 + 0.7321i 4.1463 + 0.0000i
当N=8 时计算结果分别为:
X1 =
Columns 1 through 5
12.0848 + 0.9318i 4.1977 - 2.9399i -1.3322 - 3.7609i -3.6975 -
0.0000i -1.3322 + 3.7609i
Columns 6 through 8
4.1977 + 2.9399i 12.0848 - 0.9318i 29.4783 + 0.0000i
X2 =
Columns 1 through 5
2.5554 + 2.2279i 1.8637 - 1.1417i -0.4848 - 1.7721i -1.7502 -
0.0000i -0.4848 + 1.7721i
Columns 6 through 8
当N=16 时计算结果分别为:
X1 =
1.0e+002 *
Columns 1 through 5
0.2394 + 0.3432i 0.2346 + 0.0775i 0.0883 + 0.0001i 0.0680 -
0.0619i 0.0178 - 0.0796i
Columns 6 through 10
-0.0276 - 0.0692i -0.0585 - 0.0396i -0.0694 - 0.0000i -0.0585 +
0.0396i -0.0276 + 0.0692i
Columns 11 through 15
0.0178 + 0.0796i 0.0680 + 0.0619i 0.0883 - 0.0001i 0.2346 -
0.0775i 0.2394 - 0.3432i
Column 16
1.2072 + 0.0000i
X2 =
Columns 1 through 5
2.5675 + 8.1495i 3.9564 + 2.8984i 3.6344 + 0.0469i 2.4751 -
1.7258i 0.8794 - 2.5275i
Columns 6 through 10
-0.7271 - 2.3729i -1.9126 - 1.4172i -2.3489 - 0.0000i -1.9126 +
1.4172i -0.7271 + 2.3729i
Columns 11 through 15
0.8794 + 2.5275i 2.4751 + 1.7258i 3.6344 - 0.0469i 3.9564 -
2.8984i 2.5675 - 8.1495i
Column 16
40.4692 + 0.0000i
当N=32 时计算结果分别为:
X1 =
1.0e+002 *
Columns 1 through 5
0.4796 + 1.5601i 0.4611 + 0.6776i 0.4309 +
文档评论(0)