信号的抽取与插值许多数字信号处理应用中,一个系统中会存在不同采样频率的信号,称这样的系统为多速率信号处理系统数字音频:广播、数字压缩光盘和数字音频磁带使用的采样率不同,音频系统中经常需要在几种不同的采样率之间转换多媒体通信,所传输的信号包括音频、视频和数据等不同类型的信号,这些信号的带宽相差很大,因此它们的采样率也不同增加信号采样率的运算叫内插(interpolation)或上采样(up-sampling)降低采样率的运算叫抽取(decimation)或下采样(down-sampling)多速率信号处理系统的核心内容是信号采样率转换与滤波器组整数因子抽取整数因子插值采样率的分数倍转换抽取和插值的多相滤波器结构主要内容1整数因子抽取)(nx1D)(mxD把原序列每隔个点抽取一个点组成一个新序列DmxmxD其中,D为正整数↓D)(nx)(mxD051015202530-1-0.500.51InputSequenceTimeindexnAmplitude0123456789-1-0.500.51Outputsequencedown-sampledby3AmplitudeTimeindexn抽取后是否会丢失信息?)(mxD)(nx为此分析序列与原序列之间的频谱关系)(nxsf)(mxDDfs/)2/(Dfs以因子D进行抽取后不产生混叠的频率范围为2/sf如果要采样不产生混叠,则的频率限制在以内)(nx的频率分量如果x(n)中含有大于)2/(Dfs必然产生频谱混叠,导致不能从中恢复)(mxD)(nx0)(nxn...123456789-1-2-3-4-5-6-7-8-9...0)(nDn...123456789-1-2-3-4-5-6-7-8-9...0)(nxintn...123456789-1-2-3-4-5-6-7-8-9...0)(Dmxm...123-1-2-3...nDDnnxnx其它0,......2,,0)()(int)()()()()(intnnxiDnnxnxDi其中)()(iDnniD)()()(intDmxDmxmxDnnmDmDzDnxzmxzX)()()(DnnnnznxzDnx/intDint)()(的整数倍是由于n不是D的整数倍时)(intnx=0DnDnDnnDznnxznxzX)()()()(/int)(01110niDniDnWDDnrDDrDnnrDDrnDzWDnxzX101)()()(1])[(1110110DrDDrnDrDnDrzWXDzWnxD)e(1)e(2j10jDrDrDXDX)e(jDX可见,抽取序列的频谱是原序列频谱D倍展宽后按(2)的整数倍位移并叠加而成0)e(jXππ2ππ2π41π40)e(2/jXππ2ππ2π41π40)e(2/)2j(Xππ2ππ2π41π40)e(jDXππ2ππ2π41π42抽取产生频谱混叠为保证抽取后不丢失信息,应该保证抽取后信号的频谱不发生混叠,所以在抽取前对信号的频谱应该加以限制。如果在抽取前先用一数字低通滤波器(滤波器的带宽为/D,对x(n)进行滤波,使得DX/π||0)e(j则按因子D抽取后信号频谱不会发生混叠序列抽取不发生混叠的Nyquist条件)(nx)(mxD↓D)e(jDH其他0/π|||)e(|jDDHD0)e(jXππ2ππ2π41π4π3π30)e(3/jXπππ61π60)e(3/π)2j(Xπ21π4π6π60)e(3/π)4j(X1π6π6π2π40)e(3/jDXπππ61π632整数因子插值nIImImxmxI其它0,......2,,0)/()()(nxsf)(mxIsIf)(nx)(mxI↑InnmImIzInxzmxzX)/()()(0)(nxn1230m123567894)(mXIm0123567894)(mXI)(a)(b)(c)()(InInzXznx)e()e(jjIIXX内插后信号频谱被压缩了I倍0)e(jXππ2ππ20)e(jIXπ2ππ2ππ0)e(jIXππ2ππ2π/Iπ/I/I/I(a(((((b((((c(((()e(jIX)e(jX中不仅包含的高频成分还含有大于/I的基带部分为了能从中得到基带信号,需要通过低通滤波消除内插带来的镜像。)e(jIX低通滤波器的幅度响应为其他0/||1|)e(|jIHI)(nxI)e(jIH()()ymxm0)(nxn1230m123567894()Ixmm0123567894()Ixm)(a)(b)(c插值可以提高时域分辨率3采样率的分数倍转换假设输入输出信号采样率的转换因子为DIR/可以通过将D抽取与I插值结合来实现为了保证不丢失信息,分数倍采样率转换应先插值后抽取)(nxI)e(jIH)e(jDHD)(my)(kssSIffsSfDIf)/(Sf)(a)(b)(nxI)e(jHD)(my其他0]/π,/πmin[|||)e(|jDIDH)e(jH[例]在图示系统中,已知输入信号x(n)的采样率为4kHz,其频谱幅度如图所示。试确定该系统中的I、D和幅度特性,使输出信号y(m)的抽样频率为5kHz)(nxI)e(jHD)(my解:由于输出信号比输入信号的抽样频率提高了5/4倍,所以R=5/4取I=5,D=40|)(|jeX18.08.0024.0116.056.064.096.016.096.0|)(|jIeX0|)(|jeY164.064.0(a)(b)(c)π||π16.00π16.0||4|)(|jeH7.4多抽样率系统的多相滤波结构)(nx)(mxD↓D)e(jDH)(nxI)e(jIH)(my先滤波后抽取先插值后滤波对运算速度的要求相当高7.4.1抽取器与插值器的恒等变换恒等变换1↓D1()xn2()xn1()ym↓D↓D1()xn2()xn()ym2112恒等变换2两个信号分别定标以后再相加后的抽取等于它们各自抽取后再定标和相加12↑I↑I()xn↑I2()xn1()ym2()ym2()ym1()ym1信号分别定标以后的插值等于它们各自插值后再分别定标恒等变换3↓D()xn()ymDz2()ym↓D()xn1z信号延迟D个样本后作D倍抽取和先作D倍抽取再延迟一个样本是等效的'()()xnxnD'()()DXzzXz证明:设则x’(n)y’(m)11'101()()()DrDDDrYzXzXWzD11101()()DrDrDDDDrWzXWzD11101()DrDDrzXWzD'()()()DymxmxDm2()'(1)ymym12()()YzzYz1101()()DrDDrYzXWzD1111201()()()DrDDrYzzYzzXWzD12()()YzYz恒等变换4↑I()xn1z()ym↑IIz()ym()xn信号经过单位延迟后作I倍插值和先作I倍插值再延迟L个样本是等效的)(1)(110DrDDrDzWXDzX)(1)(1101DrDDrzWXDzV)()(1)()()(11011DrDDrzWXzHDzVzHzY)()()(2zXzHzVD)(])[(1)(11102DrDDDrDDrzWXzWHDzY)(nx)(1my↓D)(zH)(nx)(2my↓D)(DzH)(2nv)(1mv1rDDW)()(1)()(1)(1101102DrDDrDrDDrzWXzHDzWXzHDzY因为所以)()(21zYzY将信号按D因子抽取后再经过滤波器H(z),等价于将信号先经过滤波器H(zD),再按D因子抽取恒等变换5)(nx)(3my↑I)(IzH)(nx)(4my↑I)(zH)(4mv)(3nv)()(IIzXzX因为)()()(3zXzHzV)()()()(33IIIzXzHzVzY)()(4IzXzV)()()()()(44IIIzXzHzVzHzY)()(43zYzY将信号先经过滤波器H(z)后再按I因子内插,等价于将信号按I因子内插后再经过滤波器H(zI)恒等变换67.4.2抽取和插值的多相滤波器结构nnznhzH1111DDllDllDllzDDlhzDlhzDlhDllDDllDllzDDlhzzDlhzzDlh11)1(1令1,,1,0,DjzjDlhzEllj则DjDjjzEzzH10多相分量I型多相分解jenhDnjDjDjjzEzzH10)(nx)(my)(DzE0D1z1z)(DzE1......)(DDzE11z)(nx)(1my↓D)(zH)(nx)(2my↓D)(DzH)(2nv)(1mv)(nx)(my)(zE0D1z1z)(zE1......)(zED11zDD先滤波后抽取先抽取后滤波降低了对处理速度的要求,提高了实时处理能力减小运算误差的积累,提高运算精度降,低对字长的要求,每一支路上滤波器的阶数只有H(z)阶数的1/D)(nx)(mxD↓D)e(jDHnnznhzH1111IIllIllIllzIIlhzIlhzIlhIjIjjzEz10jIj1令IjIIjjIzEz110)1(定义zEzRjIj1IjIjjIzRzzH101II型多相分解适合于插值IjIjjIzRzzH101)(nx)(my)(0IzRI1z1z)(1IzR.........)(1IzRI1z)(nxI)(jIeH)(my)(nx)(3my↑I)(IzH)(nx)(4my↑I)(zH)(4mv)(3nv)(nx)(my)(zR01z1z)(zR1......)(1zRIIII1z...先插值后滤波先滤波后插值每一支路上滤波器的阶数只有H(z)阶数的1/I降低了对处理速度的要求,提高了实时处理能力减小运算误差的积累,提高运算精度,降低对字长的要求