低通插值滤波器说明书.docxVIP

n TheUniversityofSouthChina 数字信号处理课程设计说明书学院名称指导教师一班级学号学生姓名一2010年6月 设计一个按因子1=5的内插器，要求镜像滤波器通带最大衰减为 0.1dB，阻带最小衰减为30dB，过渡带宽不大于兀/20，设计FIR滤 波器系数h（n） 一、初始设计（1）幅度指标可以两种方式给出。第一种，叫做绝对指标，它提出了对幅度回应函数IH（jw）l的要求。这些指标一般可直接用于FIR滤波器。第二种方法叫做相对指标，它以分贝（dB）值的形式提出要求，其定义为： dB=-20logHS）10 101H（ejw）l max 经过定义中所包含的归一化，所有滤波器的相对幅频特性最高处的值为0dB，由于定义式中有一个负号，幅频特性小的地方，其dB值反而是正的。 绝对指标： [0,wp]段叫通带，61是在理想通带中能接受的振幅波动或（容限） [ws,]段叫做阻带，62是阻带中能接受的振幅波动或（容限） [wp,ws]叫做过渡带，在此段上幅度回应通常没有限制，也可以给些弱限制。 低通滤波器的典型幅度指标 相对指标（dB）: 七是通带波动的dB值; 气是阻带衰减的dB值。 由于绝对指标中的|H（jw）|max=（1+81），因此e1-8八 夫广-20%*1〉0， 1A=-20log”么1s101+8 1 -R 1-1020逆向的关系为8=-R_1+1020cc-A-A 8=(1+8)1020牝10201低通FIR滤波器阶数的估计-20lg(前)-13 N牝Y_^l14.6(①一①)/2兀滤波器结构分析： 整数倍内插器的FIR直接实现 整数I倍内插是在已知的相邻两个原采样点之间等间隔插入I-1个新的采样值。对已知的采样序列x(nT)进行D/A转移，得道原来的模拟信号x(t)，然后再对x(t)进行较高采11aa 样率的采样得到y(nT)，这里T=IT2I为大于1的整数，称为内插因子。 整数倍内插是先在已知采样序列x(nT)的相邻两个样点之间等间隔插入I-1个0值点，11然后进行低通滤波器，即可获得I倍内插的结果。内插方案如图所示： 图中fI表示在X(njp相邻样点之间插入I-1个0值采样，称为零值内插器。 x(nT)、火nT)的傅里叶变换为：X(幻叫)、Y(e川?)，二者均为周期函数，若二1122者都用模拟频率Q表示，贝0X(e1=X(幻电)，周期为。皿=2兀/〈；Y(e2)=Y(。哄)，周期为。广2兀/T=2兀/(〈/I)=I。1。 T一一x(n—p)当n=0,±I,±210其它 V(2)的傅里叶变换k(E为V(ejw2)=黄v(nT)e-形％=X(幻屯)=X(幻叫)22勺=一8要想从V(ejw2)得到Y(e2),就必须滤除这些镜像频谱。所以，要求滤波器h(nT)的理想低通幅频特性如图所示。实际工作中，。凶2。。所以允许H(e^)有一定的过渡带，可用线性相位FIR滤波器实现。根据其功能，将h(叩°)称为镜像滤波器。 salsalsal一兀/10兀/1兀2兀镜像滤波器的理想幅频特性 将理想镜像滤波器的阻带截止频率换算成数字频率为。t—兀T—兀~2^~2—Tj—71 所以，理解情况下，镜像滤波器h(nT)的频率响应特性为冗 C0lw1H(ejw2)]——-oil七1兀 式中，C为定标系数，因此输出频谱为CX(ejiw2)olwl-Y(e2){^7-0-lw2l-而后得出：定标系数C=I。 内插器的输入、输出关系y(nT)=党v(mT)h(nT—mT)222222m=—s 频域输入、输出关系Y(ejw2)=V(ejw2)H(ejw2)镜像滤波器俱Q采用FIR结构时，】倍内插器的皿直接实现结构如图所示: x(nT)- x(nT)- TI w(nT) h(0) 11k 22~~ki 7-11 htr^\1 2 2 K) h(0)内插器的线性相位FIR h(0) 内插器的线性相位FIR直接实现结构 ▼ ——h[(N-2)T]— — 】F h[(N-1)T] ? Z-1 2 ZT 2 按整数因子I内插系统的直接型FIR滤波结构 当满足线性相位条件［二h(N-1-n2)T2)］时，可用线性相位结构实现，将乘法次数再减少一半。取N=9，画出内插旗的线性相位FIR直接高效实现结构，如图所示: x(nT)11 多相滤波器结构 按整数因子I内插系统的高效FIR滤波器结构可以用一组较短的多相滤波器组实现。如果FIR滤波器总长度为M=NI，则多相滤波器组由I个长度为N=M/I的短滤波器构成。且I个短滤波器轮流分时工作。 从m=0开始，整数因子I内插系统的输出序列y(m)计算如下： y(m)=劳p(n)x(j-n)=p(n)*x(n)kk n=0式中，m=jI+k;k=0,1,2,....,I-1;j=0,1,2,....。显然