20世纪70年代初,滤波器组首先在语音编码和多路复用中得到了应用。所谓滤波器组就是将信号按不同的频带分解成一组信号,并将分解得到的子带信号重新合成原信号的处理系统。80年代后期,M.J.Smith等人提出了运用谱分解法设计无失真地重建原始信号的滤波器方法,但存在相位失真。......
2023-06-23
为降低或解决截断误差,国内外有些专家和学者提出了一些改善措施。对截断信号进行“加窗”是一种常用的方法。以FIR滤波器设计为例,有限长度低通滤波器等于理想低通滤波器与矩形窗的乘积,频域则为差值函数与矩形窗频谱的卷积。
由图1-6可以看出,调整窗口长度N可以有效地控制过渡带宽度,但是波动幅度没有得到改善。而选择不同形状的窗函数使其谱函数的主瓣包含更多的能量,相应旁瓣的幅度就会减少,而旁瓣的减少则可使通带、阻带波动减少,从而加大阻带衰减,但这是以增大过渡带为代价的。
图1-6 理想低通滤波器加矩形窗
重叠法是由Malvar等人在DCT基础上提出的解决“方块效应”的方法。普通的块变换在把输入信号分块后,块与块之间相互独立地进行变换运算,而且其变换矩阵是方阵,即输出数与输入数相等。而在重叠式变换中,信号在被分块之后,每块信号需要和相邻块的部分信号一起做变换运算,因此其变换矩阵的输出数M和输入数L一般不等,通常的取法是使L=2M,即进行50%重叠处理。对于大多数图像,双正交重叠变换的压缩效果甚至好于Daubechies9/7小波。用DFT计算线性卷积中的重叠保留法是重叠算法的另一个应用。它将相临数据段的数据重叠。设输入数据段长度为N1,系统响应长度为N2。每个相继的输入段都由N1-N2+1个新点与前一段保留下来的N2-1个点组成。在每次运算得到的N1+N2-1个点中只有后面N2个点的结果正确,而舍去结果的前N1-1点。数据块的任何重叠率都是可行的,而其中以50%的重叠效率最高。重叠算法另一个重要应用是图像的二维重叠滤波。用常规的滤波器对图像处理时,在块的边缘会出现方块锯齿效应,虽然可以通过加大计算矩阵来减小锯齿,但这样就加大了计算量,同时锯齿效应仍然存在。这是由于进行二维离散变换时所取的方块尺寸总是有限的,用许多有限方块分别做处理后拼成一个整幅图像总会有边缘效应。二维重叠滤波将传统方法重叠进行,在原来处理矩阵的基础上向左右或上下移动一位再分别滤波,最后的输出是这些重叠滤波的叠加。由于向不同方向移动后块的高频误差是相反的,因此重叠后斜线处的方块锯齿效应消失。
全相位滤波法则是使用最大重叠的新方法,它考虑了对某个输入样本的所有长度为N的分段情况,如图1-7所示。
图1-7 全相位滤波的输入对应的所有输出
假定分段信号的长度为N,则其中任何一个输入点x(n)都可以出现在N个数据段中,x(n)在第k个数据段中的位置是k,因此,x(n)将遍历数据段中的所有时刻。第k段输入所对应x(n)的输出是x(n)在输入段中第k个位置时的结果,将每段N个输出结果中对应x(n)的输出算术平均作为系统输入为x(n)时的输出,这就是全相位滤波的基本思想。
显然,N阶全相位滤波中的输入段是采用(N-1)/N重叠的,即最大程度重叠。这样,对应数据的连续输入将产生连续的输出,且消除了常规方法中由于分段所产生的截断效应。它的另一个突出特点是便于对信号频率成分进行控制,对于满足一定条件的全相位滤波器,可实现完全线性相位,幅频特性则表现为通带和阻带波纹消失且衰减增大。这里所说的“全相”的概念是对时刻而言,因为滤波器的输入在遍历所有时刻后才得到相应的输出。
有关全相位数字信号处理方法及MATLAB实现的文章
20世纪70年代初,滤波器组首先在语音编码和多路复用中得到了应用。所谓滤波器组就是将信号按不同的频带分解成一组信号,并将分解得到的子带信号重新合成原信号的处理系统。80年代后期,M.J.Smith等人提出了运用谱分解法设计无失真地重建原始信号的滤波器方法,但存在相位失真。......
2023-06-23
图4-12 变换域信号处理图对于N阶系统,时刻n输入xn=[x,x(n-1),…其中X1是1维N阶全相位输入信号,即式就是全相位数字信号处理得到的输出与输入的数学关系表达。把式称为全相位变换,矩阵A称为全相位变换基,矩阵T则为全相位变换核。综合以上,变换域全相位数字信号处理一种可能实现如图4-16所示。......
2023-06-23
系统原理同图2-26所示,对应的输出为式中式中Tk的表达式与apDCT时相同。此外,apDWT具有与apDCT相同的传输特性函数形式,因此也具有严格零相位的特性。图4-11 基于apDWT子带分解信号图实际上,不论是何种正交变换,按照全相位信号处理方法得到的ap都可以对信号进行去除分段效应的滤波和处理,它是建立在新概念上的一种滤波系统。......
2023-06-23
国际标准建议CCITT H.261、JPEG、MPEG中都采用8×8的块做DCT,但基于DCT的数字滤波尚不普遍。为使DCT所具有的优良特性能在数字滤波中得到充分发挥,结合全相滤波的概念提出了DCT域ap系统设计。TN-2TN-1],其传输特性为显然,apDCT具有严格零相位。实现无窗、单窗和双窗apDCT与传统滤波器特性的MATLAB代码如下:显然,apDCT具有严格零相位。利用apDCT的各子带滤波器对图像的一行进行子带滤波的结果如图4-8所示。......
2023-06-23
由于不同的分段起点,可能存在N个不同的数据段都包含待插值的样点I,即将产生N个不同的插值结果,全相位插值就是取这N个插值的均值作为最终结果。正反变换定义如下:式中Walsh变换基只有±1离散值,傅里叶变换基有复数,因此选择DCT基作为插值的正交变换。图8-23 八点全相位DCT内插核时频图(N=5)对与上面实验中相同的Chirp信号进行插值重建,实验结果如图8-24所示。......
2023-06-23
在一维信号处理中,用来抑制噪声干扰并保留原有有用信号的滤波方法以及获得各种信号变换的数字滤波器,有着极其重要和广泛的应用。所以,统计结果表明二维图像信号谱能量90%以上集中在二维低频区域,如果进行2∶1或更高亚取样,谱也未发生重叠。为了克服边缘效应,一种有效的方法是对图像进行重叠滤波,在重叠滤波中的子图像在行列方向上逐点进行重叠。......
2023-06-23
为进一步改善系统性能,充分发挥窗函数在滤波器设计中的作用,综合上面两种情形得到在DFT前和IDFT后均加窗的全相位数字滤波器结构如图2-26所示。图2-28 加双窗对比单窗滤波结果对比上面介绍的3种带窗全相位数字滤波器具有完全相同的FIR滤波器等价结构,区别仅仅在于卷积窗的生成方法。为了便于讨论,将第1种类型称无窗全相位数字滤波器,第2种和第3种类型称单窗全相位数字滤波器,第4种类型称双窗全相位数字滤波器。......
2023-06-23
图2-22 加前单窗的全相位数字滤波器记前窗序列为F=[f(N-1)…因此,对比无窗全相位数字滤波器的唯一区别就是卷积窗的大小和形状。从图中可看出,单窗全相位数字滤波器极大地改善了通带及阻带内的波纹数量且降低了波纹幅度,最小旁瓣衰减也降低了6.8dB。图2-23 加前单窗对比无窗全相位数字滤波器特性在MATLAB代码中可以设置不同的窗函数来调整滤波器特性,设置不同的H来对信号进行滤波。......
2023-06-23
相关推荐