第五章IIR滤波器的设计-重庆大学电气工程学院.pptVIP

第四章  快速傅里叶变换（FFT） 主要内容 1. 基-2FFT算法； 2. 离散傅里叶反变换的快速算法； 3. 混合基FFT算法； 4. 线性调频Z变换； 5. 线性卷积与线性相关的FFT算法。 直接利用DFT进行谱分析时，存在一个突出矛盾，即当序列长度N较大时计算量大、计算时间长、数据占用内存多，难以利用DFT进行实时处理，其应用受到很大的限制。 1965年库利(Cooley)和图基(Tukey)提出了一种DFT的快速算法，这就是FFT。FFT算法使计算量大大降低，计算时间减少，特别是当序列长度N较大时，效果更为显著。 FFT并不是一种新的变换形式，它只是DFT的一 种快速算法。并且根据对序列分解与选取方法的不同而产生了FFT的多种算法。 图(a)为实现这一运算的一般方法，它需要两次乘法、两次加减法。考虑到-bW和bW两个乘法仅相差一负号，可将图(a)简化成图(b)，此时仅需一次乘法、两次加减法。图(b)的运算结构像一蝴蝶通常称作蝶形运算结构简称蝶形结，采用这种表示法，就可以将以上所讨论的分解过程用流图表示。 以N=8的频率抽取为例 4.2 IDFT的快速计算方法 IDFT与DFT的差别： 1）把DFT中的每一个系数 改为 。 2）再乘以常数 1/N ， 因此： 第一种方法：以上所讨论的时间抽取或频率抽取的FFT运算均可直接用于IDFT运算，当然，蝶形中的系数 应改为 。 4.3 N为复合数的FFT（混合基算法） 以P=3，Q=4, N=12为例 (1) 先将 x(n)通过 x(n1Q+n0)改写成 x(n1,n0)。因为 Q=4, n1=0,1,2, n0=0,1,2,3,故可得： x(0,0)=x(0) x(0,1)=x(1) x(0,2)=x(2) x(0,3)=x(3) x(1,0)=x(4) x(1,1)=x(5) x(1,2)=x(6) x(1,3)=x(7) x(2,0)=x(8) x(2,1)=x(9) x(2,2)=x(10) x(2,3)=x(11) (2) 求Ｑ个Ｐ点的DFT (3) X1(k0,n0)乘以 得到X1′(k0,n0)。 (4) 求P个Ｑ点的DFT，参变量是k0 计算量估算： (1)求Ｑ个P点DFT需要ＱP2次复数乘法和Ｑ·P·(P-1)次复数加法； (2)乘Ｎ个W因子需要Ｎ次复数乘法； (3)求P个Ｑ点DFT需要PQ2 次复数乘法和P·Ｑ·(Ｑ-1)次复数加法。 总的复数乘法量: QP2+N+PQ2=N(P+Q+1)； 总的复数加法量: Q·P(P-1)+P·Q·(Q-1)=N(P+Q-2) 当组合数 N=P1P2P3…Pm中所有的Pi均为4时，就是基-4FFT算法。 以N=43为例，第一级运算的一般形式为： 4.4 Chirp-z变换 一、算法原理 1、Chirp-z变换的计算步骤 3、Chirp-z变换的特点 用FFT计算线性卷积的方法如下: 1、计算X(k)=FFT[x(n)] 2、计算H(k)=FFT[h(n)] 3、计算Y(k)=H(k)X(k)，k=0～L-1 4、计算y(n)=IFFT[Y(k)]，n=0～L-1 可见，只要进行二次FFT，一次IFFT就可完成线性卷积计算。 计算表明，L32时，上述计算线性卷积的方法比直接计算线卷积有明显的优越性，因此，也称上述循环卷积方法为快速卷积法。 1、 重叠相加法 计算步骤： a.? 事先准备好滤波器参数 H(k)=DFT[h(n)]，N点; b.? 用N点FFT计算Xi(k)=DFT[xi(n)]; c.? Yi(k)=Xi(k)H(k); d.? 用N点IFFT求yi(n)=IDFT[Yi(k)]; e.? 将重叠部分相加。 2、重叠保留 则可以证明，rxy(n)的离散傅里叶变换为 其中: X1(k)=DFT[x1(n)], X2(k)=DFT[x2(n)], Rx1x2(k)=DFT[rxy(n)]，0≤k≤N-1 利用FFT求线性相关的步骤为： x=[1 3 -1 1 2 3 3 1];y=[2 1 -1 1 2 0 -1 3];k=length(x);xk=fft(x,2*k);yk=fft(y,2*k);rm=real(ifft(conj(xk).*yk));rm=[rm(k+2:2*k) rm(1:k)];m=(-k+1):(k-1);stem(m,rm);xlabel(m);ylabel(幅度);