CN101741350A - 对无限冲激响应滤波器的非线性进行补偿的方法和装置 - Google Patents
对无限冲激响应滤波器的非线性进行补偿的方法和装置 Download PDFInfo
- Publication number
- CN101741350A CN101741350A CN200910241825A CN200910241825A CN101741350A CN 101741350 A CN101741350 A CN 101741350A CN 200910241825 A CN200910241825 A CN 200910241825A CN 200910241825 A CN200910241825 A CN 200910241825A CN 101741350 A CN101741350 A CN 101741350A
- Authority
- CN
- China
- Prior art keywords
- sequence
- linear
- thr
- iir filter
- phase deviation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Amplifiers (AREA)
Abstract
本发明提供一种对无限冲激响应滤波器的非线性进行补偿的方法和装置。所述方法包括:在无限冲激响应(IIR)滤波器的输入序列中增加线性预畸变分量,实现对所述输入序列的线性预畸变;从所述IIR滤波器的输出序列中去除所述线性预畸变分量,得到去线性畸变后的序列;对所述去线性畸变后的序列进行相位偏移,以消除所述IIR滤波器的小数倍相位偏移。本发明通过对无限冲激响应滤波器的暂态振荡和相位偏移进行补偿,能够在基本上不增加系统复杂度的前提下,实现稳定且无相位偏移的滤波过程。
Description
技术领域
本发明属于数字信号处理领域,特别涉及一种对无限冲激响应滤波器的非线性进行补偿的方法和装置。
背景技术
在信号处理中经常碰到滤波器的应用,随着互补金属氧化物半导体(CMOS)器件和数字信号处理(DSP)技术的飞速发展,数字滤波器得到了广泛的应用。
数字滤波器从实现方法上可分为无限冲激响应(IIR)滤波器和有限冲激响应(FIR)滤波器两种,其中,FIR滤波器没有极点,只有零点,因此系统比较稳定,而且能保证很好的线性相位,但是很难做到较好的通带特性和阻带衰减特性,为了达到好的特性通常阶数非常高,这使得计算量和系统复杂度特别高。
IIR滤波器很容易实现通带特性和带外的衰减特性,并且通常阶数很低,节省系统资源,但是,IIR滤波器的稳定性不好,而且带内相位是非线性的。
发明内容
本发明所要解决的技术问题是提供一种对无限冲激响应滤波器的非线性进行补偿的方法和装置,使得在基本上不增加系统复杂度的前提下,实现稳定且无相位偏移的滤波过程。
为解决上述技术问题,本发明提供技术方案如下:
一种对无限冲激响应滤波器的非线性进行补偿的方法,包括:
在IIR滤波器的输入序列中增加线性预畸变分量,实现对所述输入序列的线性预畸变;
从所述IIR滤波器的输出序列中去除所述线性预畸变分量,得到去线性畸变后的序列;
对所述去线性畸变后的序列进行相位偏移,以消除所述IIR滤波器的小数倍相位偏移。
上述的方法,其中,所述在IIR滤波器的输入序列中增加线性预畸变分量为:
xd(n)=x(n)-d(n)
其中,x(n)为所述输入序列,d(n)为所述线性预畸变分量,xd(n)为线性预畸变后的输入序列,n=1,...,N,N为输入序列的长度;
并且,所述线性预畸变分量d(n)为如下的斜坡函数:
其中,Athr为门限电平。
上述的方法,其中,所述门限电平Athr根据对所述斜坡函数的谱泄漏对滤波性能的影响,以及所述IIR滤波器的暂态振荡对滤波性能的影响进行折中得到。
上述的方法,其中,所述从所述IIR滤波器的输出序列中去除所述线性预畸变分量为:
xo(n)=xf(n)+d(n)
其中,xf(n)为所述IIR滤波器的输出序列,xo(n)为所述去线性畸变后的序列。
上述的方法,其中,所述对所述去线性畸变后的序列进行相位偏移包括:
估计所述IIR滤波器的相位偏移值;
根据所述相位偏移值对所述去线性畸变后的序列进行相位偏移。
上述的方法,其中,所述估计所述IIR滤波器的相位偏移值包括:
A1:选取长度为Np的数据作为训练序列,由所述IIR滤波器对所述训练序列进行滤波后得到第一滤波序列xl(n),对所述训练序列进行理想滤波后得到第二滤波序列xref(n);
A2:按照如下公式对所述第一滤波序列xl(n)进行相位偏移:
其中,xt(n)为相位偏移后序列,P为相位偏移因子,Pref为初始相位偏差,(-M,M)为相位偏移扫描区间;
A3:计算相位偏移后序列xt(n)与第二滤波序列xref(n)的均方误差,并获取其中的最小值;
A4:判断所述最小值是否小于误差门限,若是,将对应的相位偏移因子作为所述相位偏移值,否则,增加M的值后返回步骤A2。
上述的方法,其中,所述根据所述相位偏移值对所述去线性畸变后的序列进行相位偏移为:
其中,Popt为相位偏移值,xp(n)为对所述去线性畸变后的序列进行相位偏移后的序列。
一种对无限冲激响应IIR滤波器的非线性进行补偿的装置,包括:
线性预畸变单元,用于在IIR滤波器的输入序列中增加线性预畸变分量,实现对所述输入序列的线性预畸变;
去线性畸变单元,用于从所述IIR滤波器的输出序列中去除所述线性预畸变分量,得到去线性畸变后的序列;
相位偏移单元,用于对所述去线性畸变后的序列进行相位偏移,以消除所述IIR滤波器的小数倍相位偏移。
上述的装置,其中,所述在IIR滤波器的输入序列中增加线性预畸变分量为:
xd(n)=x(n)-d(n)
其中,x(n)为所述输入序列,d(n)为所述线性预畸变分量,xd(n)为线性预畸变后的输入序列,n=1,...,N,N为输入序列的长度;
并且,所述线性预畸变分量d(n)为如下的斜坡函数:
其中,Athr为门限电平。
上述的装置,其中,所述从所述IIR滤波器的输出序列中去除所述线性预畸变分量为:
xo(n)=xf(n)+d(n)
其中,xf(n)为所述IIR滤波器的输出序列,xo(n)为所述去线性畸变后的序列。
上述的装置,其中,所述相位偏移单元进一步用于:
估计所述IIR滤波器的相位偏移值;
根据所述相位偏移值对所述去线性畸变后的序列进行相位偏移。
本发明实施例通过对无限冲激响应滤波器的暂态振荡和相位偏移进行补偿,能够在基本上不增加系统复杂度的前提下,实现稳定且无相位偏移的滤波过程。
附图说明
图1为本发明实施例的对无限冲激响应滤波器的非线性进行补偿的方法流程图;
图2为本发明实施例中对输入数据进行线性畸变前后的序列示意图;
图3为本发明实施例中相位补偿的原理示意图;
图4为本发明实施例中估计IIR滤波器的相位偏移值的方法流程图;
图5为本发明实施例中对训练序列的均方误差计算结果示意图;
图6为本发明实施例的对无限冲激响应滤波器的非线性进行补偿的装置结构图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图及具体实施例对本发明进行详细描述。
参照图1,本发明实施例的对无限冲激响应(IIR)滤波器的非线性进行补偿的方法,主要包括如下步骤:
步骤101:在IIR滤波器的输入序列中增加线性预畸变分量,实现对所述输入序列的线性预畸变;
具体地,采用如下公式在IIR滤波器的输入序列中增加线性预畸变分量:
xd(n)=x(n)-d(n)
其中,x(n)为所述输入序列,d(n)为所述线性预畸变分量,xd(n)为线性预畸变后的输入序列,n=1,...,N,N为输入序列的长度;
并且,所述线性预畸变分量d(n)采用如下的斜坡函数:
其中,Athr为门限电平。
所述门限电平Athr根据对所述斜坡函数的谱泄漏对滤波性能的影响,以及所述IIR滤波器的暂态振荡对滤波性能的影响进行折中得到:
若所述斜坡函数的谱泄漏对滤波性能影响较为明显,则降低所述门限电平,直到其影响接近IIR滤波器的暂态振荡对滤波性能的影响;
若IIR滤波器的暂态振荡对滤波性能的影响较为明显,则提高所述门限电平,直到其影响接近斜坡函数的谱泄漏对滤波性能的影响。
以长期演进(LTE)系统中信道估计(CHE)频域插值(6X)过程中的抗镜像滤波为例,线性预畸变结果如图2所示。
所述线性预畸变过程完成后,将xd(n)输入到所述IIR滤波器中进行滤波,并且,滤波过程包括去除整数倍相位偏移。也就是说,IIR滤波器的输出序列为已经消除了整数倍相位偏移的序列。
步骤102:从所述IIR滤波器的输出序列中去除所述线性预畸变分量,得到去线性畸变后的序列;
具体地,采用如下公式从所述IIR滤波器的输出序列中去除所述线性预畸变分量:
xo(n)=xf(n)+d(n)
其中,xf(n)为所述IIR滤波器的输出序列,xo(n)为所述去线性畸变后的序列。
通过上述步骤,消除了IIR滤波器的暂态响应带来的非理想振荡。经实验验证,暂态响应补偿后的CHE性能跟补偿前的性能相比,均方误差从0.36%降低到了0.0046%。
步骤103:对所述去线性畸变后的序列进行相位偏移,以消除所述IIR滤波器的小数倍相位偏移。
本步骤中,可以采用公知的方法进行所述相位偏移。
另外,本发明实施例还提供如下的优选方法来进行所述相位偏移:
首先,估计所述IIR滤波器的相位偏移值;
然后,根据所述相位偏移值对所述去线性畸变后的序列进行相位偏移。
其中,该优选方法进行相位补偿的原理如图3所示:利用序列相邻采样点的相关性来计算无偏值。
参照图4,所述估计所述IIR滤波器的相位偏移值包括如下步骤:
步骤401:选取长度为Np的数据作为训练序列,由所述IIR滤波器对所述训练序列进行滤波后得到第一滤波序列xl(n),对所述训练序列进行理想滤波后得到第二滤波序列xref(n);
步骤402:按照如下公式对所述第一滤波序列xl(n)进行相位偏移:
其中,xt(n)为相位偏移后序列,P为相位偏移因子,Pref为初始相位偏差,(-M,M)为相位偏移扫描区间;
步骤403:计算相位偏移后序列xt(n)与第二滤波序列xref(n)的均方误差,并获取其中的最小值;
具体地,均方误差的计算公式为:
参照图5,不同的m对应不同的均方误差值,扫描所有的均方误差值,从中找出最小的均方误差值,并记录此时m的取值Ma。
步骤404~406:判断所述最小值是否小于误差门限,若是,将对应的相位偏移因子作为所述相位偏移值,否则,增加M的值后返回步骤402。
当所述最小值小于误差门限时,将对应的相位偏移因子作为所述相位偏移值:
在得到所述相位偏移值后,就可以采用如下的公式对所述去线性畸变后的序列进行相位偏移:
其中,xp(n)为对所述去线性畸变后的序列进行相位偏移后的序列。
经实验验证,进行相位补偿后的CHE性能跟补偿前的性能相比,均方误差进一步降低到了0.00099%
参照图6,本发明实施例的对无限冲激响应滤波器的非线性进行补偿的装置,包括,线性预畸变单元、去线性畸变单元和相位偏移单元。其中:
线性预畸变单元,用于在IIR滤波器的输入序列中增加线性预畸变分量,实现对所述输入序列的线性预畸变。
具体地,采用如下公式在IIR滤波器的输入序列中增加线性预畸变分量:
xd(n)=x(n)-d(n)
其中,x(n)为所述输入序列,d(n)为所述线性预畸变分量,xd(n)为线性预畸变后的输入序列,n=1,...,N,N为输入序列的长度;
并且,所述线性预畸变分量d(n)采用如下的斜坡函数:
其中,Athr为门限电平。
所述门限电平Athr根据对所述斜坡函数的谱泄漏对滤波性能的影响,以及所述IIR滤波器的暂态振荡对滤波性能的影响进行折中得到:
若所述斜坡函数的谱泄漏对滤波性能影响较为明显,则降低所述门限电平,直到其影响接近IIR滤波器的暂态振荡对滤波性能的影响;
若IIR滤波器的暂态振荡对滤波性能的影响较为明显,则提高所述门限电平,直到其影响接近斜坡函数的谱泄漏对滤波性能的影响。
去线性畸变单元,用于从所述IIR滤波器的输出序列中去除所述线性预畸变分量,得到去线性畸变后的序列。
具体地,采用如下公式从所述IIR滤波器的输出序列中去除所述线性预畸变分量:
xo(n)=xf(n)+d(n)
其中,xf(n)为所述IIR滤波器的输出序列,xo(n)为所述去线性畸变后的序列。
相位偏移单元,用于对所述去线性畸变后的序列进行相位偏移,以消除所述IIR滤波器的小数倍相位偏移。
所述相位偏移单元可以采用公知的方法进行所述相位偏移。
另外,本发明实施例还提供如下的优选方法来进行所述相位偏移:
首先,估计所述IIR滤波器的相位偏移值;
然后,根据所述相位偏移值对所述去线性畸变后的序列进行相位偏移。
其中,相位偏移值的估计方法可以参照上文方法实施例中的描述以及图4所示。在得到所述相位偏移值后,就可以采用如下的公式对所述去线性畸变后的序列进行相位偏移:
其中,xp(n)为对所述去线性畸变后的序列进行相位偏移后的序列。
综上所述,本发明实施例通过对无限冲激响应滤波器的暂态振荡和相位偏移进行补偿,能够在基本上不增加系统复杂度的前提下,实现稳定且无相位偏移的滤波过程。
最后应当说明的是,以上实施例仅用以说明本发明的技术方案而非限制,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神范围,其均应涵盖在本发明的权利要求范围当中。
Claims (11)
1.一种对无限冲激响应IIR滤波器的非线性进行补偿的方法,其特征在于,包括:
在IIR滤波器的输入序列中增加线性预畸变分量,实现对所述输入序列的线性预畸变;
从所述IIR滤波器的输出序列中去除所述线性预畸变分量,得到去线性畸变后的序列;
对所述去线性畸变后的序列进行相位偏移,以消除所述IIR滤波器的小数倍相位偏移。
2.如权利要求1所述的方法,其特征在于,所述在IIR滤波器的输入序列中增加线性预畸变分量为:
xd(n)=x(n)-d(n)
其中,x(n)为所述输入序列,d(n)为所述线性预畸变分量,xd(n)为线性预畸变后的输入序列,n=1,...,N,N为输入序列的长度;
并且,所述线性预畸变分量d(n)为如下的斜坡函数:
其中,Athr为门限电平。
3.如权利要求2所述的方法,其特征在于:
所述门限电平Athr根据对所述斜坡函数的谱泄漏对滤波性能的影响,以及所述IIR滤波器的暂态振荡对滤波性能的影响进行折中得到。
4.如权利要求2所述的方法,其特征在于,所述从所述IIR滤波器的输出序列中去除所述线性预畸变分量为:
xo(n)=xf(n)+d(n)
其中,xf(n)为所述IIR滤波器的输出序列,xo(n)为所述去线性畸变后的序列。
5.如权利要求4所述的方法,其特征在于,所述对所述去线性畸变后的序列进行相位偏移包括:
估计所述IIR滤波器的相位偏移值;
根据所述相位偏移值对所述去线性畸变后的序列进行相位偏移。
6.如权利要求5所述的方法,所述估计所述IIR滤波器的相位偏移值包括:
A1:选取长度为Np的数据作为训练序列,由所述IIR滤波器对所述训练序列进行滤波后得到第一滤波序列x1(n),对所述训练序列进行理想滤波后得到第二滤波序列xref(n);
A2:按照如下公式对所述第一滤波序列x1(n)进行相位偏移:
其中,xt(n)为相位偏移后序列,P为相位偏移因子,Pref为初始相位偏差,(-M,M)为相位偏移扫描区间;
A3:计算相位偏移后序列xt(n)与第二滤波序列xref(n)的均方误差,并获取其中的最小值;
A4:判断所述最小值是否小于误差门限,若是,将对应的相位偏移因子作为所述相位偏移值,否则,增加M的值后返回步骤A2。
7.如权利要求6所述的方法,其特征在于,所述根据所述相位偏移值对所述去线性畸变后的序列进行相位偏移为:
其中,Popt为相位偏移值,xp(n)为对所述去线性畸变后的序列进行相位偏移后的序列。
8.一种对无限冲激响应IIR滤波器的非线性进行补偿的装置,其特征在于,包括:
线性预畸变单元,用于在IIR滤波器的输入序列中增加线性预畸变分量,实现对所述输入序列的线性预畸变;
去线性畸变单元,用于从所述IIR滤波器的输出序列中去除所述线性预畸变分量,得到去线性畸变后的序列;
相位偏移单元,用于对所述去线性畸变后的序列进行相位偏移,以消除所述IIR滤波器的小数倍相位偏移。
9.如权利要求8所述的装置,其特征在于,所述在IIR滤波器的输入序列中增加线性预畸变分量为:
xd(n)=x(n)-d(n)
其中,x(n)为所述输入序列,d(n)为所述线性预畸变分量,xd(n)为线性预畸变后的输入序列,n=1,...,N,N为输入序列的长度;
并且,所述线性预畸变分量d(n)为如下的斜坡函数:
其中,Athr为门限电平。
10.如权利要求9所述的装置,其特征在于,所述从所述IIR滤波器的输出序列中去除所述线性预畸变分量为:
xo(n)=xf(n)+d(n)
其中,xf(n)为所述IIR滤波器的输出序列,xo(n)为所述去线性畸变后的序列。
11.如权利要求10所述的装置,其特征在于,所述相位偏移单元进一步用于:
估计所述IIR滤波器的相位偏移值;
根据所述相位偏移值对所述去线性畸变后的序列进行相位偏移。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009102418251A CN101741350B (zh) | 2009-12-09 | 2009-12-09 | 对无限冲激响应滤波器的非线性进行补偿的方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009102418251A CN101741350B (zh) | 2009-12-09 | 2009-12-09 | 对无限冲激响应滤波器的非线性进行补偿的方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101741350A true CN101741350A (zh) | 2010-06-16 |
CN101741350B CN101741350B (zh) | 2012-03-28 |
Family
ID=42464338
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009102418251A Active CN101741350B (zh) | 2009-12-09 | 2009-12-09 | 对无限冲激响应滤波器的非线性进行补偿的方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101741350B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102082559A (zh) * | 2010-12-02 | 2011-06-01 | 广东宝莱特医用科技股份有限公司 | 线性相位iir滤波器的一种实现方法 |
CN102857453A (zh) * | 2012-09-18 | 2013-01-02 | 苏州东奇信息科技有限公司 | 用于消除absk信号冲击滤波响应起始振荡的方法 |
CN103152013A (zh) * | 2013-03-27 | 2013-06-12 | 北京众谱达科技有限公司 | 车载在用无线电设备检测系统滤波器单元 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3611472B2 (ja) * | 1999-02-02 | 2005-01-19 | 松下電器産業株式会社 | 適応等化回路 |
CN100593900C (zh) * | 2007-03-23 | 2010-03-10 | 展讯通信(上海)有限公司 | 直流偏差去除装置和方法 |
CN101505139B (zh) * | 2008-11-17 | 2011-08-17 | 深圳市云海通讯股份有限公司 | 一种数字线性预畸变方法和系统 |
-
2009
- 2009-12-09 CN CN2009102418251A patent/CN101741350B/zh active Active
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102082559A (zh) * | 2010-12-02 | 2011-06-01 | 广东宝莱特医用科技股份有限公司 | 线性相位iir滤波器的一种实现方法 |
CN102082559B (zh) * | 2010-12-02 | 2013-08-21 | 广东宝莱特医用科技股份有限公司 | 线性相位iir滤波器的一种实现方法 |
CN102857453A (zh) * | 2012-09-18 | 2013-01-02 | 苏州东奇信息科技有限公司 | 用于消除absk信号冲击滤波响应起始振荡的方法 |
CN102857453B (zh) * | 2012-09-18 | 2015-09-30 | 苏州东奇信息科技股份有限公司 | 用于消除absk信号冲击滤波响应起始振荡的方法 |
CN103152013A (zh) * | 2013-03-27 | 2013-06-12 | 北京众谱达科技有限公司 | 车载在用无线电设备检测系统滤波器单元 |
CN103152013B (zh) * | 2013-03-27 | 2016-03-30 | 北京众谱达科技有限公司 | 车载在用无线电设备检测系统滤波器单元 |
Also Published As
Publication number | Publication date |
---|---|
CN101741350B (zh) | 2012-03-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107346974B (zh) | 减少模数转换器中的非线性的方法和装置 | |
US9521024B2 (en) | Device of crest factor reduction | |
CN101378251B (zh) | 一种提高低采样率情况下信号显示质量的方法 | |
KR20140085556A (ko) | 블록 기반 파고율 저감 | |
CN101741350B (zh) | 对无限冲激响应滤波器的非线性进行补偿的方法和装置 | |
US9496965B2 (en) | Chromatic dispersion compensator with integrated anti-aliasing filter and resampler | |
CN104143967A (zh) | 用于级联积分梳状抽取器的补偿滤波器的设计方法和系统 | |
WO2012145955A1 (en) | Channel impulse response (cir)/dc offset (dco) joint estimation block and method | |
CN105471800A (zh) | 一种基于叠接相加的f-ofdm多子带频域滤波器 | |
US8284877B2 (en) | Noise filtering inverse transformation | |
WO2019120562A1 (en) | Altering filtering by changing mixing frequency when interferer present | |
CN103002197B (zh) | 一种信号数据处理方法、装置及智能终端 | |
CN107342750B (zh) | 适用于多奈奎斯特区的分数延迟优化方法及其实现结构 | |
CN103401831A (zh) | 一种基于fpga的lte数字中频削峰方法及装置 | |
CN103001907A (zh) | 一种数字削峰方法和装置 | |
CN102891662A (zh) | 一种通用的速率下变换、上变换装置及方法 | |
CN110971197B (zh) | 数字预失真装置及方法 | |
US9595925B2 (en) | Distortion-compensating power amplifier and method for compensating for distortion to amplify power | |
CN102347920A (zh) | 一种增强的信道估计方法及装置 | |
EP3235147B1 (en) | Chromatic dispersion compensation filter | |
CN101771638B (zh) | 信道估计的方法和装置 | |
CN106936757B (zh) | 前导符号信号长度的选取方法 | |
Liang et al. | A Wiener model based Post-calibration of ADC nonlinear distortion | |
JP2011510525A (ja) | フィルタ装置 | |
Liu et al. | A novel oversampling scheme for design of hybrid filter bank based ADCs |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |