北京交通大学(数字信号处理生)ch7_1信号的抽取与内插探究.ppt

内插等式 例: L=M=2 基本单元的连接 例: L=2, M=3 如M和L互素，即M和L无公因子，则上述两种级联等价。 基本单元的连接 基本单元的连接 如M和L互素，即M和L无公因子，则上述两种级联等价。 基本单元的连接 当M和L互素时，有S1=S2 * * * * * * 基本单元 基本单元 近代数字信号处理(Advanced Digital Signal Processing) 电子信息工程学院 信号与图像处理研究室 信号的内插与抽取 抽取滤波器和内插滤波器 多相分解 半带滤波器 两通道滤波器组 多速率信号处理基础 常见的音频信号抽样频率 48kHz DAT (Digital Audio Tape) 44.1kHz CD 32kHz 数字广播 24kHz DAT 抽样频率的一半 22.1kHz CD 抽样频率的一半 16kHz 数字广播抽样频率的一半 问题的提出 利用抽样频率为16 kHz的播放系统，播放抽样频率分别为32kHz、16kH和8kHz的音频信号。 播放抽样频率为 32 kHz的信号 播放抽样频率为 16 kHz的信号 播放抽样频率为 8 kHz的信号 问题 32 kHz信号播放出的歌曲速率比正常情况慢 8 kHz信号播放出的歌曲速率比正常情况快 问题分析: 16 kHz 系统播放抽样频率 32 kHz信号 抽样频率为32kHz的离散信号 播放系统输出的连续信号 y(t)=x(t/2) 连续信号 抽样频率为8kHz的离散信号 播放系统输出的连续信号y(t)=x(2t) 连续信号 问题分析: 16 kHz 系统播放抽样频率 8 kHz信号 问题分析: 16kHz 系统播放抽样频率 16 kHz信号 抽样频率为16kHz的离散信号 播放系统输出的连续信号 连续信号 y(t)=x(t) 问题延伸 : 16kHz 系统播放抽样频率 24 kHz信号 播放系统输出的连续信号 抽样频率为24kHz的离散信号 连续信号 y(t)=x(2t/3) 问题延伸 : 16kHz 系统播放抽样频率 24 kHz信号 播放系统输出的连续信号 抽样频率为24kHz的离散信号 连续信号 y(t)=x(2t/3) 问题的提出 N=481100 问题的提出 N=288660 学 习 要 求 掌握序列抽取运算与内插运算的频谱变化规律。 掌握确定抽取滤波器与内插滤波器的频率指标。 掌握有理数倍抽样率转换的原理及方法。 了解抽取滤波器和内插滤波器的多级实现的基本思想。 了解半带滤波器的基本特性以及设计方法。 掌握两通道滤波器组的基本概念，以及两通道FIR PR滤波器组的设计方法和主要步骤。 抽取和内插的时域描述 抽取和内插的变换域描述 基本单元的连接 抽取等式 内插等式 信号的内插与抽取 (a) M倍抽取(down-sampler, M-fold decimator) 抽取和内插的时域描述 例：3倍抽取 例： M倍抽取是时变系统。 N=40; w0=0.6*pi; M=2 k=0:N-1; x=sin(w0*k); y=downsample(x,M); subplot(2,1,1); stem(k,x); title(x[k]); subplot(2,1,2); stem(0:length(y)-1,y); title(y[k]); 利用MATLAB实现序列抽取 0 5 10 15 20 25 30 35 40 -1 -0.5 0 0.5 1 x[k] 0 2 4 6 8 10 12 14 16 18 20 -1 -0.5 0 0.5 1 y[k] 利用MATLAB实现序列抽取 2倍抽取后的序列 (b) L倍内插(interpolation, up-sampler, L-fold expander) 抽取和内插的时域描述 例：3倍内插 N=20; w0=0.1*pi; L=3 k=0:N-1; x=sin(w0*k); y=upsample(x,L); subplot(2,1,1); stem(k,x); subplot(2,1,2); k=0:L*N-1; stem(k,y); 利用MATLAB实现序列内插 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 利用MATLAB实现序列内插 (a) M倍抽取 抽取和内插的变换域描述 M倍抽取后频谱的变换规律 证明 序列抽取不混叠的条件 3倍抽取后信号的频谱 X