DSP-应用实例-课件教程.ppt

第8章 应用实例分析;教学内容及要求;本章主要内容;8.1 FIR滤波器的C54X实现方法;的冲击响应系数。因此可使用Matlab语言实现这一过程，得到滤波器的冲击相应系数；再将冲击响应系数应用到DSP汇编语言程序中，以实现上面的计算公式，就可以实现FIR滤波器，达到滤波效果。图8.1是N阶FIR滤波器的方框图。图中表明了实现滤波器的整个过程。 图8.1 N阶FIR滤波器的原理图 ;8.1.2 FIR滤波器程序设计;npoints=length(k)： t=0:npoints-1： des(k)=linspace(aa(2*i-1)，aa(2*i)，npoints)： if i=1 wt(k)=wtx(i)*(1.5+cos((t)*pi/(npoints-1)))： elseif i=nbands wt(k)=wtx(i)*(1.5+cos(pi+(t)*pi/(npoints-1)))： else wt(k)=wtx(i)*(1.5-cos((t)*2*pi/(npoints-1)))： end end ;;;8.2 IIR滤波器的C54X实现方法;;;8.2.2 IIR滤波器实现举例;;;;;;;;;;;;;;;;;;;;;8.3 快速傅里叶变换（FFT）的C54X实现方法; 对于有限长离散数字信号它的离散频谱可由离散傅里叶变??求得。DFT定义如下 ; 式中（有时简写为W）代表。不难看出，是周期性 的，且周期为N，即 ; 来说，由DFT的定义可以看出，在为复数序列的情况下，直接运算N点DFT需要次复数乘法和次复数加法。因此，对于一些相当大的N值（如1024点）来说，直接计算它的DFT的计算量是很大的。一个优化的实数FFT算法是一个组合以后的算法。原始的2N个点的实输入序列组合成一个N点的复序列，然后对复序列进行N 点的FFT运算，最后再由N点复数输出拆散成2N点的复数序列，这2N点的复数序列与原始的2N点的实数输入序列的DFT输出一致。FFT的基本思想在于：将原来的N点序列分成两个较短的序列，这些序列的DFT可很简单地组合起来得到原序列的DFT。 ; 例如，若N为偶数，将原有的N点序列分成两个（N/2）点序列，那么计算N点DFT将只需要约[（N/2）2*2]=N2/2次复数乘法。即比直接计算少作一半乘法。因子（N/2）2表示直接计算（N/2）点DFT所需要的乘法次数，而乘数2代表必须完成两个DFT。上述处理方法可以反复使用，即（N/2）点的DFT计算也可以化成两个（N/4）点DFT（假定N/2为偶数），从而又少作一半乘法。使用这种方法，在组合输入和拆散输出的操作中，FFT的运算量减半。这样，利用实数FFT算法来计算实输入序列的DFT的速度几乎是一般复FFT算法的两倍。假定序列的点数是2的整数次幂，按上述处理方法， 定义为; 的偶数项和奇数项的（N/2）点序列和;因考虑到 可写成： ; 上式表明，N点DFT可分解为按上式的规则加以组合的两个（N/2）点DFT。如果直接计算（N/2）点的DFT，则不难看出，计算N点DFT需要（N2/2+N）次复数乘法。在N值很大的情况下，这表示可以节约一半左右的计算时间。下面用流程图来说明两个4点DFT变换计算8点DFT的处理方法。计算流程如图8.19所示。输入序列先进行位倒序成偶数项和奇数项得到;图8.19 8点FFT蝶形图; 对于64点的FFT运算，基于上述的思想，可以将64个点分为32个点，32个点也可继续分为16个点，如此继续分下去，一直分到最基本的两点傅里叶变换，如图8.20所示。; 由此可得出8点基2FFT的流程图，如图8.21所示。由图8.21可推出任意2的整数次幂的FFT的流程图。;8.3.2 FFT的C语言实现;struct compx EE(struct compx b1，struct compx b2) { struct compx b3： b3.real=b1.real*b2.real－b1.imag*b2.imag： b3.imag=b1.real*b2.imag+b1.imag*b2.real： return (b3)： } 4. 定义位反函数。 ;void displace(void)： { …… nv2=N/2： for(i=1：i=nm1：i++) { if(ij) {t=xin[j]：xin[j]=xin[i]：xin[i]=t：} k=nv2： while(kj){j=j－k：k=k/2： {;if(ij) {t=xin[j]：xin[j]=xin[i]：xin[i]=t：} k=nv2： while(kj){j=j－k：k=