抽取与内插滤波器

合集下载

PCM／CVSD系统中内插和抽取高效滤波算法

ＳＨＥＮｉｍｉＨｕ．ｎ
（ｈ４ｈＲｓｒｎｔｒＥＣ，ｈｉｈａｇｔｂｉ．０１Ｃｉ）Ｔｅ５ｔｅａｃｓｔｅｆＣＴＳｉｚｕｎｔｅ００８，ｈａｅｈｌｉａｏｊａｅ５ｎ
ＡｂｔａｔＡｔｆｓ，ｈｐｌｃｔｏａｋｒｕｄｏｎｔｒｏａｉｎａｄｄｃｍａｉｎｆｔｒｔｅｒａｔｍｅｃｎｅｓｏｙｔｍｅｗｅｎＰＣＭｎｄｓｒｃｒｔｔｅａｐｉａｉｎｂｃｇｏｎｆｉｅｐｌｔｎｅｉｔｌｅｓ，ｈｅｌｉｏｖｒｉｎｓｓｅｂｔｅｉｏｏｉａ
ＣＤ，ｎｒｄｃｄ．ｈｎ，ｉｌｄｅｉｉｎｅｉｎｍｅｄｏｔｒｏａｉｎａｄｄｃｍａｏｌｒａｅｎａ￣ｆｄｆｔｒｉｐｅｅｔ，ＶＳｉｉｔｕｅＴｅａｓｓｏｍｐｅａｆｃｅｔｄｓｇ￣ｏｆｉｅｐｌｔｅｉｔｎｆｔｓｂｓｄｏｎｎｏｎｉｉｅｌｍｎｅｓｒｓｎｅｌｉｄ
音采用Ａ律６ｂｓＰＭ编码体制。为了解决这４ｋｐ的Ｃ２不同语音编码方式的通信网之间的互通，要种需对内插和抽取滤波器的设计以半带滤波器为基
础，特征介绍如下：其
① 半带ＦＲ滤波器的冲激响应除了零点不为Ｉ
ｈｌｎｅｔｏａｌｄｍｔｒｍｅｈｄｂａ
０引言
某专用网话音部分采用１ｂｓ的连续可变斜６ｋｐ
１内插和抽取滤波器的设计

抽取与内插滤波器

抽取与内插滤波器
2020年4月22日星期三
抽取滤波器和内插滤波器
抽取滤波器 • 2倍抽取滤波的矩阵表示
内插滤波器 2倍内插滤波的矩阵表示
抽取滤波器
X(ejW)
-p
可用理想低通滤波器滤除X(ejW)中的高频分量
W
p
但理想低通滤波器无法实现。
抽取滤波器
X(ejW)
W
-p
p
若Wm/M 为X(ejW)中需保留的最高频率分量，则有
第0列 h0[k] =h[2k] 第2列 h0[k-1] = h[2k-2]
第4列 h0[k-2] = h[2k-4] 第2n列 h0[k-n] = h[2k-2n]
第1列第3列
h-1[k]=h[2k -1] h-1[k-1] = h[2k-3]
第5列 h-1[k-2] = h[2k-5] 第2n+1列 h-1[k-n] = h[2k-(2n+1)]
2倍内插滤波器的时域表示
内插滤波器的时域表示
例：2倍抽取滤波器hR[k]=h[-k]的矩阵表示
例：2倍内插滤波器hR[k]=h[-k]的矩阵表示
利用MATLAB 计算抽样率变换
(1) 抽取 y = decimate(x,M)
用8阶Chebyshev I 型 IIR 低通滤波器进行滤波。为保证零相位，对序列进行正向和反向滤波。
第4行 h0[2-n]= h[4-2n]
第2k行 h0[k-n]= h[2k-2n]
第-1行第1行第3行第2k-1行
h-1[-n] = h[-1-2n] h-1[1-n]= h[1-2n ] h-1[2-n]= h[3-2n ] h-1[k-n] = h[2k-1-2n]

可视电话系统中用于图像抽取与内插的FIR滤波器设计

可视电话系统中用于图像抽取与内插的FIR 滤波器设计罗虹(宁夏工业设计院有限责任公司)【摘要】实现了可视电话系统中应用于视频格式CCIR601与QCIF 相互转换的抽取与内插的FIR 数字滤波器的设计,并在已实现的H.263编解码系统中使用了这些滤波器,实验结果获得了很好的图像质量。

【关键词】抽取;内插;FIR 滤波器【Ab str act 】An implementatio n of decimatio n and interpo latio n FIR digital filters is propo sed in this paper .T hese filters can be applied into the conv ersion between CCIR601PAL and QCIF.A test has been don e in a H.263codec sy stem and the result v erifies these filters are very suitable for this v ideo format co nversio n.【Key wo r ds 】Decim ation;Interpo lation;FIR filter1.前言在可视电话H.263编解码系统[5]中,需要完成视频格式CCIR601PAL(704*576)与QCIF(176*144)的相互转换,这一过程主要是要完成图像数据采样率的转换。

在数字采样率变换方法中,提高采样率称为内插;降低采样率称为抽取。

从CCIR601变换到QC IF,需要降低采样率,因而采用抽取,相反地,从QCIF 变换到CCIR601,采用内插。

对于抽取和内插,我们分别设计了FIR 数字滤波器,实践表明,这组滤波器非常适合CCIR601与QCIF 的相互转换,可获得很好的图像质量。

图1X(!)Y(!)抽取频谱图Fig.1Figu re s o f decima tion s pe ctrumX(!)Y(!)2.抽取与内插的工作原理2.1抽取的工作原理以4:1下取样为例,设输入信号(数字序列)为x (n),输出信号为y(n),则y (n)=x (4n)n=0,1,2,…根据Fourier 变换的性质,x(n)、y (n )的频谱Y(!)、X(!)满足如下关系[1]Y(!)=14[X(!4)+X(!4-"2)+X(!4-")+X(!4-32")]图1(a)、(b)分别为X(!)、Y(!)频谱图。

ch8_1信号的抽取与内插

4.
5.
掌握抽取滤波器和内插滤波器的多级实现的基本思想。
了解半带滤波器的基本特性以及设计方法。
6. 了解两通道滤波器组的基本概念，以及两通道FIR PR滤波器组的设计方法和主要步骤。
重点和难点
本章的重点是信号的抽取和内插的时域、频域及z域分析
本章的难点是两通道滤波器组的分析与设计
信号的内插与抽取

x[k / L]z k
x[n]z nL

n
X I ( z) X ( z L )
基本单元
XI(ej)= X(ejL)
XI(ej)= X(ejL)
L=5时内插序列的频谱
1 X(ej)
p
镜像 1
XI(ej)
镜像
p

p
p

p
p

基本单元的联接
M
N
L L
基本单元
xD[k]
0 1 2 3 k
基本单元
例： M倍抽取是时变系统。
x[k ]
1 0 2 3 4 5
x D [k ], M 2
k
0
1
2
k
y[k ] x[k 1]
2 0 1 3 4 5 6
y D [k ], M 2
k
0 1 2 3
k
利用Matlab实现对正弦序列抽取
N=40; w0=0.6*pi; M=2 k=0:N-1; x=sin(w0*k); y=x(1:M:end); subplot(2,1,1); stem(k,x); title('x[k]'); subplot(2,1,2); stem(0:length(y)-1,y); title('y[k]');

数字信号处理7-2抽取滤波器和内插滤波器

抽取滤波器的基本概念抽取滤波器的时域表示内插滤波器的基本概念内插滤波器的时域表示分数倍的抽样速率转换
M=2
抽取滤波器的基本概念
X(ej) 1
3 2/3 2/3
3
XD(ej)
1/2
3
序列抽取M倍不混叠的条件:
3
X(ej)=0，||>/M
x[k ]
H(z)
M
y[k ]
H(z) 2
D/A
fsam=32kHz
frec=16kHz
x(t)
x[k]
t
k
连续信号
抽样频率为32kHz的离散信号
问题解决：16 kHz 系统播放抽样频率 32 kHz信号
x[k]
w[k]
y(t)
x(t) A/D
H(z) 2
D/A
fsam=32kHz
frec=16kHz
w[k] k
频率转换后的离散信号
问题解决： 16 kHz 系统播放抽样频率 24 kHz信号
x(t)
x[k ]
A/D
2
fsam 24kHz
w[k ]
y(t)
H(z) 3
D/A
frec 16kHz
x(t)
连续信号号
问题解决： 16 kHz 系统播放抽样频率 24 kHz信号
...
/L /L
可用理想低通滤波器滤除内插后信号频谱XI(ej)中的镜像分量
H
(e
j
)
1, 0,
Ω π/L
π / L | | π
内插滤波器的基本概念
X(ej)
+m m
...
XI(ej)
m m
...

抽取及内插的频谱分析

抽取与内插的频谱分析工科试验班钟汇凯43我们知道，为了避免在抽样信号中出现混叠，抽样定理要求被抽样的信号是一个带限信号。

然而，在实际应用中，绝大多数信号都不能满足这个要求，为了减小混叠的影响以及放宽对滤波器性能指标的要求，在实际应用中往往采取一种提高抽样率的办法，使信号的抽样率远远大于限带滤波器通带频率的两倍。

例如，在下图中，当抽样频率略大于限带频率 ωm 的两倍时，混叠的影响还是很明显的，而当抽样频率远远大于两倍的 ωm 时，混叠的影响就非常之小了。

虽然提高抽样率可以减小混叠的影响，但是，在对连续时间信号进行处理的离散时间系统中，过高的抽样率将增加系统的成本，因为，过高的抽样率将要求离散时间系统以较高的速率工作，而高速率器件的成本一般都要贵于低速率的器件。

可以设想，如果能对信号的抽样率进行调整，使得在信号的抽样和恢复中使用较高的抽样率，在离散时间处理中使用较低的抽样率，那么，上述性能和成本的矛盾就可以得到适当的折中，而离散时间信号的抽取和内插就是一种调整信号抽样率的办法。

从技术性能层面来看。

这两种方法类似于连续时间信号的抽样和内插。

抽取离散时间信号的抽取包含信号抽样和尺度变换两个步骤：首先，以抽样间隔N 对离散时间信号进行抽样，然后再对抽样信号进行1/ N 的尺度压缩变换。

下图是离散时间信号的抽取过程，图中，x [ n ] 是离散时间信号，xs [ n ] 是抽样信号，抽样间隔N＝3，xd [ n ] 是抽取信号，它是xs [ n ] 进行1/N 尺度压缩变换后所得到的结果。

由图可见，在抽样信号xs [ n ] 和抽取信号xd [ n ] 之间存在以下关系:(1)由于抽样信号xs [ n ] 在N 的整数倍上和离散时间信号x [ n ] 相等，因此，式也可等效为：(2)虽然式(1)和式(2)在形式上完全相同，但两者的含义不同：式(1)的含义是，抽取信号xd [ n ] 是由抽样信号xs [ n ] 进行1/N 尺度压缩变换的结果；而式(2)的含义是，抽取信号xd [ n ] 是从离散时间信号x [ n ] 中每隔(N－1)个点取一个样本值所组成的一个新序列，这个过程就称为离散时间信号的抽取。

ch7_1信号的抽取与内插解读

x[k ]
L
y[ k ] x I [ k ]
x[ k / L ], k 0, L , 2 L , xI [ k ] 其他 0
x[k] 0 xI[k] 0 3
基本单元
1
2
3
k
6
9
k
利用MATLAB实现序列内插
N=20; w0=0.1*pi; L=3 k=0:N-1; x=sin(w0*k); y=zeros(1,L*length(x)); y(1:L:end)=x; subplot(2,1,1); stem(k,x); subplot(2,1,2); k=0:L*N-1; stem(k,y);
p
p
p
p
1

3XD(ej)
p
p
p
p

p p
序列抽取不混叠的条件
X(ej)=0，||>p/M
1
X(ej)
p
p
p
1

p X(ej)
p
p

p
p
p
1

p X(ej(p)
p
p

p
p
p

p
p
p

2XD(ej) 1
Y1 ( z ) H ( z )
1 M
M 1

l 0
l X ( z M WM )
Y2 ( z ) X ( z ) H ( z M ) M
1 M
M 1

l 0
1
1
l l X ( z M WM ) H (( z M WM )M )
H ( z ) M 1 l M X ( z W M) M l 0

抽取与内插浅说

内插和抽取浅说Freehardman， CD，2008-10-22（转载请注明出处）数字通信系统中由于要将模拟信号数字化，必须要进行A/D变换。

A/D变换所带来的问题不只是模拟信号数字化了。

更为关键的是如何最大限度的保持原有模拟信号中所包含的有用信息。

一般涉及到两个指标，一个是ADC的分辨率，简单的表示就是转换后的数字信号的位宽，另一个就是ADC的采样速率。

在ADC分辨率一定的情况下，提高采样速率是一种有效的提高采样质量的方法。

但是采样速率的提高会给后续的处理带来压力。

在这种情况下就产生了多速率信号处理技术。

所谓多速率处理就是对已有的数字样值的处理不是按照一个速率进行，而是多个速率，其本质我认为就是对数字样值数量进行的一个变换。

样值数量的改变就是增加新的样值或者减少已有的样值。

增加样值：如果样值数量增加，那么必然导致要有一个频率相应增加的时钟信号。

如果样值数量增加一倍，那么就要有一个2倍于采样时钟的信号。

依此类推。

这个过程用术语说就是内插。

其实相当于增加采样率。

减少样值：如果样值数量减少，那么必然导致要有一个频率相应降低的时钟信号。

如果样值数量减少1半，那么就要有一个频率是采样时钟频率一半的时钟信号。

依此类推。

这个过程用术语说就是抽取。

其实相当于减少采样率。

数字下变频与数字上变频：在数字通信系统中，一般信号输入端进行减少样值（抽取），在信号输出端进行增加样值（内插）。

由前面可知，抽取后，相应的时钟频率可以降低，所以这个过程也叫数字下变频；内插后，相应的时钟频率就要升高，因此这个过程又叫数字上变频。

疑问：既然要抽取，为什么不直接在ADC采样的时候就降低频率呢？我的理解是这样的：因为在硬件电路处理上，ADC的采样时钟确定后，在工作过程中就会再有采样时钟频率的变化。

如果要求采样时钟变化，则必然要有一个控制电路以及相应和指标都很好的DDS或其他可程控的时钟源，这在硬件设计和处理上都增加复杂度。

抽取对信号本身频率的影响：在时域和频域有相应的公式，可以很容易找到。

抽取和内插

多速率信号处理及抽取和内插一：多速率信号处理1、在信号处理系统中有时需要不同的抽样率，这样做的目的有时是为了适应不同系统之间的级联，以利于信号的处理、编码、传输和存储，有时则是为了节省计算工作量。

数据速率的转换两种途径：1）数字信号数模转换模拟信号模数转换另一抽样率抽样2）数字信号处理数字信号处理基本方法抽样率转换目的：改变原有数字信号的频率方法：抽取和内插，低通滤波。

低通滤波：抽取和内插的前提条件是信号频带内没有频谱混叠，实现这一点需要用到低通滤波。

2、多速率滤波器-->具有线性相位的FIR滤波器。

常用的多速率滤波器：多速率FIR滤波器，积分梳状滤波器（CIC）和半带滤波器(HB);3、常用多速率信号处理结构第一级：CIC滤波器。

用于实现抽取和低通滤波第二级：fir实现的半带滤波器优点：工作在较低频率下，且滤波器参数得到优化，更容易以较低阶数实现，达到节省资源，降低功耗的目的。

二：抽取概念：使抽样率降低的转换。

1、整数倍抽取当信号的抽取数据量太大时，为了减少数据量以便于处理和计算，我们把抽样数据每隔（D-1）个取一个，这里D是一个整数。

这样的抽取称为整数抽取，D称为抽取因子。

2、抽取后结果：信号的频谱：信号的频谱周期降低1/D;信号的时域：信号的时域每D个少了（D-1）信号。

3、抗混叠滤波：在抽取前，对信号进行低通滤波，把信号的频带限制在抽样后频率的一半以下，这样，整数倍抽取的的问题就变成了一个低通滤波的问题。

信号时域图信号频域图程序运行后所得到的滤波前后信号的时域图，滤波器的频率响应图如上图。

从图中可以看出，经半带滤波器滤波后的信号，与原信号相比，波形没有改变，但抽样速率降低了一半；半带滤波器通阻带容限相同，具有严格线性相位。

三：内插概念：使抽样率升高的转换。

1、整数倍内插：在已知的相邻抽样点之间等间隔插入（I-1）个零值点。

然后进行低通滤波，即可求得I倍内插的结果。

2、内插后结果：信号的时域：已知抽样序列的两相邻抽样点之间等间隔多了I-1个值信号的频谱：信号的频谱周期增加了I倍。

数字信号处理讲义-信号的抽取与内插

j2πl
X(e M
)
12
M倍抽取后频谱的变换规律
XD(ej)M 1M l01
2πl
j
X(e M )
X (e j
)
扩 M 倍
X
j
(e M
)
周期化 2π为
1 M1
2πl
j
X(e M )
M l0
13
证明
~M[k]
M1
1 kl WM
M l0
XD(z)x[kM ]zk
n
x[n]z M
k
n是M的整数倍
1X (ej( )
13 X D (ej )
序列抽取不混叠的条件 X(ej)=0，||>/M15
1 X(ej)
X(ej) 1
1
X(ej()
2XD(ej) 1
2倍抽取产生的频谱混叠
16
抽取和内插的变换域描述
(b) L倍内插
XI(z) xI[k]zk
Ml0
H(z)M1
M l0
1
X(zMWM l )
20
内插等式
x[k ] L
H (z L ) y3[k]
x[k ] H (z)
y4[k] L
Y3(z)X(zL)H(zL) Y4(z)X(z)H(z)LX(zL)H(zL)
21
基本单元的连接
x[k ]
L v1[k] M y1[k] ?
x[k] M v2[k] L
0
3
6
9
k
xD[k]
k
0
1
2
3
5
例： M倍抽取是时变系统。
x[k ]
xD [k], M 2

抽取和内插

数据速率的转换两种途径：1）数字信号→数模转换→模拟信号→模数转换→另一抽样率抽样2）数字信号处理→数字信号处理基本方法→抽样率转换目的：改变原有数字信号的频率方法：抽取和内插，低通滤波。

低通滤波：抽取和内插的前提条件是信号频带内没有频谱混叠，实现这一点需要用到低通滤波。

2、多速率滤波器-->具有线性相位的FIR滤波器。

常用的多速率滤波器：多速率FIR滤波器，积分梳状滤波器（CIC）和半带滤波器(HB);3、常用多速率信号处理结构第一级：CIC滤波器。

二：抽取概念：使抽样率降低的转换。

1、整数倍抽取当信号的抽取数据量太大时，为了减少数据量以便于处理和计算，我们把抽样数据每隔（D-1）个取一个，这里D是一个整数。

这样的抽取称为整数抽取，D称为抽取因子。

2、抽取后结果：信号的频谱：信号的频谱周期降低1/D;信号的时域：信号的时域每D个少了（D-1）信号。

信号时域图信号频域图程序运行后所得到的滤波前后信号的时域图，滤波器的频率响应图如上图。

三：内插概念：使抽样率升高的转换。

1、整数倍内插：在已知的相邻抽样点之间等间隔插入（I-1）个零值点。

然后进行低通滤波，即可求得I倍内插的结果。

2、内插后结果：信号的时域：已知抽样序列的两相邻抽样点之间等间隔多了I-1个值信号的频谱：信号的频谱周期增加了I倍。

ch7_1信号的抽取与内插

Y2 ( z) X ( z)H ( z M ) M

1
M 1
1
1
X ( z M WMl ) H (( z M WMl ) M )
M l0

H (z) M 1 M l0
1
X (z M WMl )
基本单元
内插等式
x[k ] L
H (z L ) y3[k]
x[k ] H (z)
y4[k] L
n
X I (z) X (z L ) XI(ej)= X(ejL)
基本单元
XI(ej)= X(ejL)
L=5时内插序列的频谱
1 X(ej)
镜像
1 XI(ej)
镜像

基本单元的连接
M
N
y[k]
x1[k]
1
x2[k]
2
M
y[k]
基本单元
基本单元的连接
x[k]
L v1[k] M y1[k] ?
x[k]
M
v2 [k ] L
y2[k]
如M和L互素，即M和L无公因子，则上述两种级联等价。
V1(ej ) X(ejL )
V2(ej )
1 M
M 1 k0
2πk j
X(e M )
Y1(ej )
Y3 (z) X (z L )H (z L )
Y4(z) X (z)H(z) L X (z L )H (z L )
基本单元
基本单元的连接
x[k]
L v1[k] M y1[k] ?
x[k]
M
v2 [k ] L
y2[k]
例: L=M=2

抽取与内插滤波器资料

n n
M倍抽取滤波的时域表示
x[k ]
H(z)
v[ k ]
M
y[k ]
v[k] x[n]h[k - n]
n
y[k ] v[k M ]
y[ k ]

n
x [ n ]h[ M k - n ]
内插滤波器(interpolation filter)
x[k ]
L
xI [k ]
H ( z)
矩阵[Ih]的第k行第n 列
[ I h ] k ,n h[k - 2n]
2倍内插滤波器的矩阵表示
y0 h0 y1 h1 y h 2 2 y3 h3 y4 y 5 y6 y7
h0 h1 h2 h0 h3 h1 h2 h0 h3 h1
x 0 x1 x2 x 3
y [Dh ]x
y0 h0 y h 1 2 y2 y3 h1 h3 h0 h2 h1 h3 h0 h2
h1 h0
2倍抽取滤波的矩阵表示
抽取矩阵[Dh]的行
h0 h2 h1 h0 [ Dh ] h3 h2 h1 h0 h3 h2 h1 h0
第0列 h0[k] =h[2k]
第2列 h0[k-1] = h[2k-2] 第4列 h0[k-2] = h[2k-4] 第2n列 h0[k-n] = h[2k-2n] 第1列

h-1[k]=h[2k -1]
第3列
h-1[k-1] = h[2k-3]
第5列 h-1[k-2] = h[2k-5] 第2n+1列 h-1[k-n] = h[2k-(2n+1)]

抽取与内插滤波器资料

滤波器的仿真工具
MATLAB
提供丰富的滤波器设计工具箱，支持多种滤波器类型和设计方法。
Python
使用SciPy、NumPy等库进行滤波器设计和仿真，具有强大的数据处理能力。
SPICE
电路仿真软件，可用于模拟电路中的滤波器设计和仿真。
仿真结果分析
幅频响应分析
观察滤波器的通带、阻带特性以及过渡带的陡峭程度。
抽取与内插滤波器资料
https://
REPORTING
• 引言 • 抽取滤波器 • 内插滤波器 • 抽取与内插滤波器的比较 • 滤波器的实现与仿真 • 总结与展望
目录
PART 01
引言
REPORTING
WENKU DESIGN
目的和背景
数字信号处理的发展
随着数字技术的飞速发展，数字信号处理已成为现代信号处理的主要手段。抽取与内插滤波器作为数字信号处理的重要组成部分，对于提高信号处理效率和质量具有重要意义。
PART 05
滤波器的实现与仿真
REPORTING
WENKU DESIGN
滤波器的实现方法
IIR滤波器实现
采用递归型结构，利用模拟滤波器设计方法进行设计，如巴特沃斯、切比雪夫等。
FIR滤波器实现
采用非递归型结构，通过窗函数法、频率采样法等进行设计。
自适应滤波器实现
根据输入信号的特性自适应地调整滤波器参数，如LMS算法、 RLS算法等。
图像处理
在图像处理中，抽取滤波器可用于图像缩放和图像压缩等领域，实现图像数据的降维和压缩。
PART 03
内插滤波器
REPORTING
WENKU DESIGN
内插滤波器的原理
采样定理