- 248
- 0
- 约7.42千字
- 约 20页
- 2018-12-28 发布于湖北
- 举报
快速傅里叶变换
如果一台通用计算机的速度为平均每次复乘5μs,每次复加0.5μs,用它来计算512点的的DFT[x(n)],问直接计算需要多少时间,用FFT运算需要多少时间。
解:
(1)直接计算
复乘所需时间
复加所需时间
所以
(2)用FFT计算
复乘所需时间
复加所需时间
所以
已知,是两个N点实序列,的DFT值,仅需要从,求,的值,为了提高运算效率,试用一个N点IFFT运算一次完成。
解:
依据题意 ,
取序列
对作N点IFFT可得序列。
又根据DFT性质
由原题可知,,都是实序列。再根据,可得
N=16时,画出基-2按时间抽取法及按频率抽取法的FFT流图(时间抽取采用输入倒位序,输出自然数顺序,频率抽取采用输入自然顺序,输出倒位序)。
解:
(1)按时间抽取,见图P4-3(a)。
图 P4-3(a)
(2)按频率抽取,见图P4-3(b)。
图 P4-3(b)
N=16时,导出基-4FFT公式、画出流图,并就运算量与基-2的FFT相比较(不计乘±1及乘±j的运算量)。
解:
依题意 N=4×4=r1r2
对于nN,有
同样令N=r1r2,对于频率变量k(kN)有
可得
根据上式
令
则
而
所以
计算量比较:
基-4:只在乘旋转因子时有复乘,复乘8次;复加64次。
基-2:复乘10次,复加64次。
基-4流图如图P4-4所示。
图 P4-4
试用N的组合数时的FFT算法求N=12的结果(采用基-3×4),并画出流图。
解:依题意
对于,有
同样 令,对于频率变量有
可得
根据上式,得
流图如图P4-5所示。最后输出为,是倒位序的,按可算出其相应的k值,在整序后,即可得正常顺序的输出。
图 P4-5
同上题。导出的结果,并画出流图。
解:依题意
对于nN,有
同样,令,对于频率变量有
可得
根据上式得
令
则
所以
流图如图P4-6所示。
图 P4-6
研究一个长度为M点的有限长序列。
我们希望计算求z变换在单位圆上N个等间隔点上的抽样,即在上的抽样。试对下列情况,找出只用一个N点DFT就能计算的N个抽样的方法,并证明之。
(1) (2)
解:
(1)依题意
设,则
因为
且令
所以
由此可见,对于,可先计算,然后对它求一次N点DFT,即可计算在单位圆上的N点抽样。
(2)若,可将补零到N点,即
则
已知一个8点序列
试用CZT法求其前10点的复频率。已知z平面路径为,,画出zk的路径及CZT实现过程示意图。
解:依题意
则
所以
而
把代入,可得
令
则
由(1)式可得的路径(模和相角),如表4-8所示。其CZT实现过程和的路径见图P4-8(a),P4-8(b)。
表 P4-8
k
0
1
2
3
4
5
6
7
8
9
|zk|
0.8
0.67
0.56
0.46
0.39
0.32
0.27
0.22
0.19
0.16
Agr[zk]
π/3
13π/30
16π/30
19π/30
22π/30
25π/30
28π/30
31π/30
34π/30
37π/30
图 P4-8(a)
图 P4-8(b)
在下列说法中选择正确的结论。线性调频z变换(CZT)可以用来计算一个M点有
限长序列在z平面的实轴上各点的z变换,使
(1)。
(2)。
(3)(1)和(2)两者都行。
(4)(1)和(2)两者都不行。即线性调频z变换不能计算在z为实数时的
原创力文档

文档评论(0)