DSP技术05电子通信10.pptxVIP

第五章 数字信号处理应用程序设计;一、FFT运算的DSP实现; 时域 —— 频域 离散信号 周期谱 周期信号 离散谱 离散周期信号 离散周期谱;;DFT的定义;;DFT和FFT计算量对比;2、N点蝶2FFT;;上式的最后一步是因为;我们已经将一个N点的DFT分解成为两个N/2点的DFT。 但是，X(k)有N点，但X1(k)和X2(k)都只有N/2点，因此，前面计算的只是X(k)的前一半项的结果。;对于后一半X(k)，有; 这样，只要计算出（0，N/2-1）区间的X1(k)和X2(k)，也就可以很方便地计算整个（0，N-1）区间的全部X(k)，从而大大地节省了运算量。;将N点DFT分为两个N/2点DFT;将N/2点DFT分为两个N/4点DFT;2点DFT;8点蝶2FFT;按频率抽取（DIF）;;DIT与DIF对比; ;定点DSP计算1024点FFT的时间;浮点DSP计算1024点FFT的时间;顺序、混序与位倒序;3、FFT运算的DSP实现;1. 基二实数FFT运算的算法;第三步，分离复数FFT的输出为奇部分和偶部分 分离FFT输出为RP、RM、IP和IM四个序列，即偶实数、奇实数、偶虚数和奇虚数四部分。 第四步，产生2N点的复数FFT输出序列 产生2N=256个点的复数输出，它与原始的256个点的实输入序列的DFT一致。;2. FFT运算模拟信号的产生及输入;3. 实序列FFT汇编源程序及链接命令文件;DATA .space 1024 .copy mdata1.inc ; mdata1.inc 为模拟输入信号数据，起始地址标号INPUT N .set 128 ;复数点数 LOGN .set 7 ;蝶形级数 sav_grp .usect tempv,3 ;定义组变量值 sav_sin .set sav_grp+1 ;定义旋转因子表 sav_idx .set sav_grp+2 OUTPUT .usect OUTPUT,256 ;信号功率谱 BOS .usect stack,0Fh ;定义堆栈 TOS .usect stack,1 .copy “twiddle1.inc” ;正弦表系数由twiddle1.inc ;文件给出，起始地址标号TWI1;.copy “twiddle2.inc“ ; 余弦表系数由twiddle2.inc文件给 ；出，起始地址标号TWI2 .text _c_int00 b start nop nop start: STM #TOS,SP LD #0,DP SSBX FRCT …………;;二、数字滤波器的DSP实现;理想低通滤波器;低通、高通、带通与带阻滤波器;高通滤波器;带通滤波器;带阻滤波器;截止频率与带宽;通带、阻带与过渡带;滚 降（过渡带）;阻带衰减;???内平坦度;巴特沃尔斯滤波器（Butterworth);;契比协夫滤波器（Chebyshev）;契比协夫滤波器（Chebyshev）I型;契比协夫滤波器（Chebyshev）II型;椭圆滤波器; ;二阶低通滤波器;幅频特性和相频特性;; 2、数字滤波器;数字滤波器的差分方程;IIR滤波器的主要特点;y[k] = a0x[k] + a1x[k-1] + a2x[k-2] –b1y[k-1] –b2y[k-2];零极点子系统交换位置（传输函数不变，但差分方程不同）;标准型二阶节;二阶节的C54实现;级 连;多级二阶节级连成多级 IIR;并 联;IIR滤波器的设计;冲激响应不变法 h(nts)= ha[k];Z平面的单位圆映射为S平面的±fs/2的区间;双线性变换;有限冲击响应（FIR）滤波器;FIR滤波器的结构;频率响应;FIR滤波器的主要特点;系统非线性相移造成输出信号失真; ;相位特性;要求线性相位的例子;线性相位的FIR滤波器;线性相位设计;假设在 时H(?) = 0，有：;截取有限项得近似H(z)：;用DSP实现FIR滤波器的存储器安排;FIR滤波器编程实现;计算16阶的FIR滤波器输出256个点;C54的FIRS指令;矩形窗;Hamming窗;Blackman窗;矩形窗与blackman窗的比较; ;数字滤波器设计软件包DFDP4/plus; 1．FIR滤波器的基本结构 ;FIR滤波器的结