数字信号处理DSP10-1.pptVIP

窄带FIR滤波器的高效实现 例：用Kaiser窗口实现一个窄带滤波器， 要求采样频率为40k； 通带截止频率[13.2k 13.8k]； 阻带截止频率[13k 14k]； 带内纹波 0.02dB ，带外衰减 60dB。 窄带FIR滤波器的高效实现 窄带FIR滤波器的高效实现 利用经验公式 可以得到 滤波器最低阶数为725阶； 该滤波器通带相对宽度小于1/20，可以进行20倍抽取；利用抽取滤波器进行设计，则该滤波器的实现成本只相当于一个36阶的滤波器！ 设计原理： 首先设计宽过渡带的低阶滤波器及其互补滤波器； 利用内插压缩频谱，使过渡带变窄； 采用频率响应掩蔽法选取通频带和阻带。 窄过渡带的宽带FIR滤波器设计 低阶滤波器 互补滤波器 目标滤波器 通过内插进行频谱压缩 窄过渡带的宽带FIR滤波器设计 设计要求： 在标准数字频率范围内，过渡带宽度为d的锐截止宽带低通滤波器； 设计步骤： 先设计截止频率为 /2的低通滤波器 ，设定其过渡带宽度为 ； 同时实现其互补滤波器 ； 窄过渡带的宽带FIR滤波器设计 通过L倍内插，使频谱压缩为1/L并周期化，各过渡带宽度变为d； 利用宽过渡带滤波器 和 作为内插滤波器（设其过渡带宽度为 /L），选出几个周期的镜像频谱进行组合，就可以实现具有锐截止过渡带的宽带滤波器； 窄过渡带的宽带FIR滤波器设计 为了减少设计难度，可将上述各滤波器相对过渡带宽设计为相等，即有： ，由此可以得出过渡带宽为d的滤波器设计时的内插倍数为 。 若采用相对过渡带宽度描述，则为： 窄过渡带的宽带FIR滤波器设计 例： 若设计滤波器要求的过渡带宽度d为0.01，则可以选取L=17，各滤波器的过渡带宽可以设置为0.17，放宽程度可以达到17倍，滤波器阶数下降为原来的1/17 。 窄过渡带的宽带FIR滤波器设计 第十章 变采样系统 抽取与内插 临界抽取滤波器组设计 数字信号处理 变采样系统设计 数字信号由采样值构成，采样频率应为信号带宽的2倍以上； 在信号处理过程中，信号的带宽会发生变化；及时调整采样率，不仅可以提高数据保存和传输效率，也可以提高数字系统的设计效率。 信号抽取（下采样） 从原始数字序列中进行等距抽取构成新的数字序列 抽取前后信号的关系 采样周期加大，数据量减少； 采样频率降低，标准频谱展宽，可能出现混叠失真！ 信号抽取的效果 抽取系统的构成 通常在抽取之前，需要先对信号进行抗混叠滤波，限制信号带宽；对于低通信号，抽取滤波器为低通滤波器，截止频率为 抽取滤波器 滤波器卷积方程 滤波器运算量为同阶普通滤波器的1/M 信号内插（上采样） 在原始数字序列中每2个数据间等距插入L-1个零，构成新的序列。 内插前后信号的关系 采样周期缩减，数据量增加； 采样频率提高，标准频谱压缩，数字频谱标准区间内出现镜像频谱！ 信号内插的效果 理想内插系统的构成 在插0之后，对信号进行抗镜像滤波，可以将插0点改为理想插值，该滤波器称为内插滤波器；对于低通信号，内插滤波器为低通滤波器，截止频率为 内插滤波器 内插滤波器 滤波器卷积方程 滤波器运算量为同阶普通滤波器的1/L 分数倍变采样系统 先对信号进行L 倍上采样，形成高采样率信号；再对信号进行M倍下采样，采样率变为L/M。 与抽取和内插有关的MATLAB函数 Decimate(x,r,n, filter)： 经过低通后进行下采样,得出输出信号； 输入参数： 输入信号序列x，下采样因子r， 低通滤波器阶数n， 滤波器类型：默认为切比雪夫滤波器 例1 用10阶切比雪夫滤波器进行3倍下采样 x=rand(100,1);r=3;n=10; y=decimate(x,r,10); n=1:100;n1=1:3:100; subplot(2,1,1),stem(n,x); title(初始信号); subplot(2,1,2),stem(n1,y); title(抽取信号); 与抽取和内插有关的MATLAB函数 与抽取和内插有关的MATLAB函数 例2 用50阶FIR低通滤波器进行5倍下取样 x=rand(200,1);r=5; y=decimate(x,r,50,fir); n=1:200;n1=1:5:200; subplot(2,1,1),stem(n,x); title(初始信号); subplot(2,1,2),stem(n1,y); title(抽取信号); 与抽取和内插有关的MATLAB函数 与抽取和内插有关的MATLAB函数 Interp(x,r,l,