第7章多采样率数字信号处理.ppt

第7章　多采样率数字信号处理 　　镜像频谱滤波器hI(n)的指标也有两项： 阻带截止频率和数字域阻带截止频率。 根据采样理论， 采样频率提高到Fy时， 采样信号序列y(m)的频谱以Fy为周期。 由于输入端信号x(n)的采样频率为Fx， 所以x(n)的频带宽度不会超过Fx/2（对应的数字频率为π）。 因此， 整数因子I内插器输出信号y(n)的频谱一定是带宽为Fx/2、 重复周期为Fy的周期谱。 我们将零值内插器输出的信号v(m)的频谱中除了y(m)频谱以外的其他频谱称为镜像频谱。 镜像频谱滤波器的作用就是让y(m)频谱尽量无失真通过， 滤除v(m)中的镜像频谱。 当I=3时， 图7.2.2中各点信号的频谱示意图如图7.2.3所示。 7.3 采样率转换系统的高效实现 实际上， 采样率转换系统的高效实现就是指其中的FIR数字滤波器的高效实现。 这里高效的含义有三个方面： 在满足滤波指标要求的同时， ① 滤波器的总长度最小； ② 使滤波处理计算复杂度最低； ③ 对滤波器的处理速度要求最低。  教材中介绍了采样率转换系统的两种实现方法： 直接型FIR滤波器结构、 多相滤波器实现。 各种实现方法的原理、 结构及其特点在教材中都有较详细的叙述， 本书不再重复。 题4解图 　　% 程序ex805.m　　% 调用remez函数设计按因子D=5的抽取器的抗混叠滤波器　f=［0.18， 1/5］； %对π归一化边界频率　　m=［1， 0］； 　　rp=0.1； rs=30； 　　dat1=(10∧(rp/20)－1)/(10∧(rp/20)+1)；dat2=10∧(－rs/20);　　　rip=［dat1， dat2］； 　　［M， fo， mo， w］=remezord(f， m， rip)； M=M+3； 　　hn=remez(M， fo， mo， w)；　　%以下为绘图部分（省略）　　运行程序， 绘图如题5*解图所示。 　　6*. 设计一个内插器， 要求内插因子I=2。 用remez函数设计镜像FIR滤波器， 图示滤波器的单位脉冲响应和损耗函数。 要求通带最大衰减为0.1 dB， 阻带最小衰减为30 dB，过渡带宽度为0.05π rad。 画出实现内插器的多相结构， 并求出多相实现时各子滤波器的单位脉冲响应。 　　滤波器单位脉冲响应h(n)数据为　　h(0)= －1.250515e－004 =h(73) 　　h(1)= －8.034005e－003 =h(72)　　h(2)= 8.680925e－003 =h(71)　　h(3)= －5.520558e－004 =h(70)　　h(4)= －4.859884e－003 =h(69)　　h(5)= －8.774395e－004 =h(68)　　h(6)= 4.896185e－003 =h(67)　　h(7)= 1.961462e－003 =h(66)　　h(8)= －5.290437e－003 =h(65)　　h(9)=－3.140232e－003 =h(64) 题6*解图 　　滤波器设计程序ex807.m如下： 　　% 程序ex807.m　　% 调用remez函数设计按因子2/5的采样率转换器中FIR滤波器　　f=［0.16， 0.2］； %对π归一化边界频率　　m=［1， 0］； rp=1； rs=30； 　　dat1=(10∧(rp/20)-1)/(10∧(rp/20)+1);dat2=10∧(－rs/20); 　　rip=［dat1， dat2］； 　　［M，fo，mo，W］=remezord(f， m， rip)；%M=M+2； 　hn=remez(M， fo， mo， W)； 　　n=0： M； y=［n； hn； M-n］； 　　fprintf(′h(%d)= %9e =h(%d)＼n′， y)　　%以下为绘图检验部分(省略) 　　h(9)=－1.272172e － 002 =h(48)　　h(10)= － 1.139736e － 002 =h(47)　　h(11)= － 5.945673e － 003 =h(46)　　h(12)= 2.482502e － 003 =h(45)　　h(13)= 1.155711e － 002 =h(44)　　h(14)= 1.836975e － 002 =h(43)　　h(15)= 2.020830e － 002 =h(42)　　h(16)= 1.549946e － 002 =h(41)　　h(17)= 4.44