- 1、本文档共3页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
FFT递推运算框图
外层的循环控制变量M=,其中N输入数据的个数,一般取N为2的幂次方的整数,如果输入的数据的个数N≠2M,则在序列后面补零,让其达到N≠2M个,N,M均可以认为是事先给定的常数
对上述的递推运算框图有如下几点说明:
1) 该递推框图实现的程序是一个子程序块,调用时A,M,N由主程序赋值,因此在框图中,可以认为A,M,N是事先给定的常数。
2) 框图中,共有级(列)蝶形,每级有N/2个蝶形。
3) 框图中,用L表示运算的级数,则第L(L=1,2,……M)级中系数有个,它们是(可=0,1,……),U中存有,其起始值U=。
4)每一级各个系数对应的蝶形有个,它们依次相距LE=点。用循环变量I来控制同一种类型的蝶算,显然步进为LE=。这就是框图中的最内层循环。同一类蝶形中参加运算的两点相距为LE1=点。用J来控制各个类型的蝶形都做到,J的步进为1,这对应于框图中的次内层循环。也可以认为框图中的最内层循环完成每级蝶式运算,次内层循环完成因子的运算。
5)由L控制运算的级数。L由1变到M步进为1,对应于框图的最外层循环,每改变一个L值,则LE1、LE、系数U都会发生变化。
框图中,U、W、T为存放复数单元,相应的乘法为复数乘法,系数用递推公式:
(为乘号)
上述FFT递推运算框图的FORTRAN程序实现如下:
SUBROUTINE FFT(A,M,N)
COMPLEX A(N),U,W,T
N=
NV2=N/2
NM1=N-1
J=1
a DO 7 I=1,NM1
IF (I·GE·J)GO TO 5
T=A(J)
A(J)=A(I)
A(I)=T
5 K=KV2
6 IF=(K·GE·J)GO TO 7
J=J-K
K=K/2
GO TO 6
7 J=J+K
PI=3.14159
b DO 20 L=1,M
LE=
LE1=LE/2
U=(1.0,0.0)
W=CMPLEX(cos(PI/FLOAT(LE1)),-sin(PI/FLOAT(LE1)))
c DO 20 J=1,LE1
d DO 10 I=J,N,LE
IP=I+LE1
T=A(IP) * U
A(IP)=A(I)-T
10
20 U=U * W
RETURN
END
上述程序中,标号为a的DO循环来完成的是整序运算,然后整个FFT运算用三个嵌套循环完成,其中标号为d的DO循环完成每级蝶式运算,标号为c的DO循环完成因子的运算,标号为b的DO循环完成M级运算。
您可能关注的文档
最近下载
- GB_T 42588-2023系统与软件工程 功能规模测量NESMA方法.docx VIP
- 祝由十三科 收藏珍本.pdf
- 高级英语第三课Ships-in-the-Desert名师优质课获奖市赛课一等奖课件.ppt VIP
- (新教材)人教版高中物理必修1第三章第5节《共点力的平衡》优质说课稿.doc
- 新冀教版数学一年级上册全册课件(2024年秋新教材).pptx
- 《化妆品配方与制备技术》课件——第10章 唇部美容化妆品.pptx VIP
- JBT 10491-2022 额定电压 450750V 及以下交联聚烯烃绝缘电线和电缆.docx
- 《化妆品配方与制备技术》课件——第1章 护肤化妆品.pptx VIP
- 山东服装职业学院综合类招聘考试真题.pdf
- 合伙经营分公司协议书(7篇).docx
文档评论(0)