CHP快速傅立叶变换.pptVIP

实验5.1分析 实验5.1分析 实验5.2分析 xn1=[[0:25].*exp(-0.3*[0:25]),zeros(1,24)]; xn2=[[0:25].*exp(-0.3*[0:25]),zeros(1,74)]; %input(请输入序列x(n):);N=length(xn); subplot(411),stem(xn1); subplot(412),stem(xn2); X1=fft(xn1); X2=fft(xn2); subplot(413);stem([0:49],abs(X1)); subplot(414);stem([0:99],abs(X2)); 实验5.2分析 实验5.3分析 N=64;f1=10.8;f2=11.75;f3=12.55;fs=40;stepf=fs/N; n=0:N-1; t=2*pi*n/fs;n1=0:stepf:fs/2-stepf; x=sin(f1*t)+sin(f2*t)+sin(f3*t); M=N;W=exp(-j*2*pi/M);A=1; Y1=czt(x,M,W,A); subplot(311);plot(n1,abs(Y1(1:N/2)));grid on; Y2=abs(fft(x)); subplot(312);plot(n1,abs(Y2(1:N/2)));grid on; M=60;f0=8;DELf=0.12;A=exp(j*2*pi*f0/fs); W=exp(-j*2*pi*DELf/fs); Y3=czt(x,M,W,A); n2=f0:DELf:f0+(M-1)*DELf; subplot(313);plot(n2,abs(Y3));grid on; 实验5.3分析 ENDThanks! * 8点DIF基2FFT算法流图 4.3 频率抽取（DIF）基2FFT算法 4.3 频率抽取（DIF）基2FFT算法 4.3 频率抽取（DIF）基2FFT算法 DIT与DIF的相同之处： （1）DIF与DIT两种算法均为原位运算。 （2）DIF与DIT运算量相同。 DIT与DIF的不同之处： （1）DIF与DIT两种算法结构倒过来。 DIF为输入顺序，输出乱序。运算完毕再运行“二进制倒读”程序。 DIT为输入乱序，输出顺序。先运行“二进制倒读”程序，再进行求DFT。 （2）DIF与DIT根本区别：在于蝶形结不同。 DIT的复数相乘出现在减法之前。 DIF的复数相乘出现在减法之后。 4.5 分裂基算法 自从基2快速算法出现以来，人们仍在不断寻求更快的算法。1984年，法国的杜梅尔(P.Dohamel)和霍尔曼(H.Hollmann)将基2分解和基4分解糅合在一起，提出了分裂基FFT算法。其运算量比前几种算法都有所减少，运算流图却与基2FFT很接近，运算程序也很短。它是目前一种实用的高效新算法。 4.5 分裂基算法 分裂基算法又称基2/4算法,算法的基本思路是:对偶号输出使用基2算法,对奇号输出使用基4算法。 下面首先介绍基4算法： 令N=4M，对N点DFT可按下面方法进行频率抽取 分别令k=4r， k=4r+2， k=4r+1， k=4r+3，其中，r=0，1，2，…N/4-1，有 4.5 分裂基算法 4.5 分裂基算法 4.5 分裂基算法 4.5 分裂基算法 算法分析 对于N=4M点DFT，当输出项k为偶数时，采用基2算法，即 当输出项k为奇数时，采用基4算法，即 4.5 分裂基算法 4.5 分裂基算法 4.5 分裂基算法 分析： 一个N点DFT可以分解为一个N/2点DFT和两个N/4点DFT。由x(n) x(n+N/4) x(n+N/2)和x(n+3N/4)求N/2点DFT和N/4的DFT只需要两次乘法，可以减少运算量。 N/2点DFT可进一步分解为一个N/4点DFT和两个N/8的DFT。N/4的点DFT进一步分解为一个N/8点DFT和两个N/16的DFT。 这样一直下去，直到分解为两点或4点DFT为止。 4.5 分裂基算法 结论： 分裂基FFT算法结构同基2FFT算法结构相似，适用于N=2M的场合，并由M级运算实现。运算流图输入为顺序，输出为倒序。 分裂基FFT算法的计算量 以上提出FFT算法，可以很快地求出全部DFT值。即求出有限长序列x(n)的z变换X(z)在单位园上N个等间隔抽样点zk处的抽样值。它要求N为高度复合数。即N可以分解成一些因子的乘积。例N=2L 实际上： （1）也许对其它围线上z变换取样发生兴趣。 （2）只需要计算单位圆上某一段的频谱,即M不等于N。如窄带信号，希望在窄带频率内频率抽样能够非常密集，提高分辨率，带外则不考虑。 （