第六讲 多采样率信号处理.ppt

6.1、采样率降低——整数Ｍ倍抽取 6.2 、 采样率提高——整数Ｌ倍内插 采样率变换一个有理因数 6.3 抽取与内插的FIR结构一、抽取的FIR结构 抽取器的多相（Polyphase）滤波器结构 二、内插的FIR结构 * 第六章 多采样率信号处理 在实际应用中，各系统之间的采样率往往是不同的，例如，在音频范围内，广播系统的采样率为32kHz，CD唱机的采样率为44.1kHz，而数字录音带（DAT）的采样率为48kHz。这就需要设计一类数字滤波器，用于将被处理信号的采样率转换成与相应系统所要求一致的采样率。多采样率技术还有许多其他的应用，过采样技术是其中的一个典型应用实例。 0 5 10 15 20 25 30 35 40 -1 -0.5 0 0.5 1 x(n) 0 2 4 6 8 10 12 14 16 18 20 -1 -0.5 0 0.5 1 y(m) 根据采样定理,为防止混叠,应滤除高频分量 直接抽取 滤波抽取 2倍抽取产生的频谱混迭 X (e j ? ) 1 p - p ? - 4p 0 - 2p 2p 4p X (e j ?’ /2 ) 1 p - p ?’ - 4p 0 - 2p 2p 4p 1 p - p ?’ - 4p 0 - 2p 2p 4p X ( e j( ?’- 2p ) /2 ) p/2 - p/2 ?’ - 2p 0 - p p 2p Y ( e j ?? ) 1/2 例: 对采样频率为100kHz的信号进行M=50倍的抽取，抽取滤波器的截止频率为fs/2M=1000Hz 。通带和阻带边缘频率为950Hz和1000Hz,对应 假定通带和阻带的波纹峰值为0.005 如果基于Kaiser窗设计符合上述要求的FIR滤波器，根据教材(4.56)式，则滤波器系数个数N=5431 现采用多级级联的抽取结构，设计相应的FIR滤波器，并比较两种不同方法的滤波器系数个数。 滤波器的设计指标如图（a）所示，采用图（b）所示的两级级联，现在使 H1(z) 25 H2(z) 2 0.019 0 1 H(f) π 0.02 (a) (b) 0.019 0 1 H(f) π 0.06 0.475 0 1 H(f) π 0.5 通带边界依然是950Hz，则 由于抽取因子是25，阻带边界可以定为 间出现混叠不会影响1kHz以下的信号 所以 通带波纹峰值为 约0.2dB。阻带最小衰减为0.005，大于47dB。 此时的滤波器长度为 见教材(4.56)式，式中 为过渡带宽。 的滤波器 ，此时的采样频率为 则滤波器的通带边界950Hz对应于 而阻带边界定为 通带波纹峰值为 约0.2dB。阻带最小衰减为0.005，大于47dB。 此时的滤波器长度为 式中 是 的十五分之一 工作于极低的采样频率上，计算量将进一步下降。 事实上 0 2 4 6 8 10 12 14 16 18 20 -1 -0.5 0 0.5 1 0 10 20 30 40 50 60 -1 -0.5 0 0.5 1 插零的结果 L=4 插零并且滤波的结果, L=4 p - p 2 p 2 - p p - p L=2 如要求y(0)=x(0)，则应有G=L，即对理想的内插器要求能恢复抽取前的信号，增益Ｇ必须等于Ｌ。 内插器 低通滤波 抽取器 问题:低通滤波的指标如何确定? 例: 数字录音带（DAT）驱动器的采样频率为48kHz,而激光唱盘（CD）播放机则以44.1kHz的采样频率工作。为了直接把声音从CD录制到DAT，需要把采样频率从44.1kHz转换到48kHz。为此，考虑完成图6—16的采样率转换系统。求L和M的最小可能值以及适当的滤波器 完成这个转换。 解：已知 及 为改变采样频率，需 所以，如果我们选择内插因子L＝160而抽取因子M＝147，便得所求采样率转换。我们所求的滤波器截至频率为 其增益应等于160，采样频率为44.1×160=7056kHz h(0) h(1) y(m) h(N-1) z-1 z-1 x(n) M h(0) h(1) y(m) h(N-1) x(n) M M M z-1 z-1 (a) FIR滤波与抽取器级联(b) 提高运算效率的结构 抽取的FIR结构 k=Mq+i，i=0,1,…M-1, q=0,1,…,N/M-1 , N取M的整数倍 令 如 y1 (m) yM-1(m) y0(m) xM-1(m) x0(m) y(m) x(n) M M M hM-1 (m) h1(m) h0(m) x1 (m) z-1 z-1 抽取器的多相滤波器结构 h(0) h(1) y(m) h(N-1) x(n) L z-1 z-1 h