CN101807902A - 复系数线性相位无限脉冲响应数字滤波器 - Google Patents
复系数线性相位无限脉冲响应数字滤波器 Download PDFInfo
- Publication number
- CN101807902A CN101807902A CN 201010121871 CN201010121871A CN101807902A CN 101807902 A CN101807902 A CN 101807902A CN 201010121871 CN201010121871 CN 201010121871 CN 201010121871 A CN201010121871 A CN 201010121871A CN 101807902 A CN101807902 A CN 101807902A
- Authority
- CN
- China
- Prior art keywords
- impulse response
- filter
- infinite impulse
- digital filter
- phase
- 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
- Networks Using Active Elements (AREA)
- Filters That Use Time-Delay Elements (AREA)
Abstract
本发明属于数字信号处理技术领域,具体为一种复系数线性相位无限脉冲响应数字滤波器。在数字滤波器中,无限脉冲响应数字滤波器很难得到线性相位特性,然而无限脉冲响应数字滤波器由于其较低的阶数受到广泛的应用。本发明用时间反转的方法设计复系数线性相位无限脉冲响应数字滤波器。本发明在推导得出复系数线性相位无限脉冲响应数字滤波器的条件的基础上,用时间反转的方法实现满足这种条件的滤波器。
Description
技术领域
本发明属于数字信号处理技术领域,具体涉及一种复系数线性相位无限脉冲响应数字滤波器。
背景技术
线性相位是滤波器的一个很重要的特性,具有线性相位的滤波器在滤波时不会对信号造成畸变。系数对称的有限脉冲响应滤波器具有严格的线性相位特性,而无限脉冲响应滤波器却很难得到严格的线性相位特性,但是无限脉冲响应滤波器由于其较低的阶数而被广泛应用,其中复系数无限脉冲响应滤波器也有着广泛的应用。
在已有的文献中,已经有一些设计复系数无限脉冲响应滤波器的方法,比如文献【1】【2】中基于全通滤波器的近似线性相位滤波器,这种方法只能得到近似的线性相位,阶数通常也比较高,而且很难对系数进行实时的重配置;再比如文献【3】-【6】中用限制优化的技术获得近似的线性相位,这也只能得到近似的线性相位;文献【7】中用局部时间反转的方法实现设计线性相位无限脉冲响应滤波器,但是不能用于复系数滤波器。
发明内容
本发明的目的在于提供一种能够用于复系数的线性相位无限脉冲响应滤波器。
考虑两个N阶复系数滤波器的传递函数:
和
其中,M和N是正整数,且M≤N,H1(z)和H2(z)的零点数目和极点数目分别相等,即分别为M和N。如果他们的零点和极点满足
这里,z1m *是复数z1m的共轭,p1n *为复数p1n的共轭。
则将这两个滤波器级联得到的滤波器具有严格的线性相位特性。上述结论的证明如下:
级联得到的滤波器的传递函数为
为了分析该滤波器的频域特性,用exp(jω)代替z,其中ω是归一化的角频率,并用rzm exp(jφzm)代替z1m,rpnexp(jφpn)替代p1n,这样就可以得到该合成滤波器的频响传输函数:
为了便于分析,将该传输函数重新分解成下面的M+N个函数:
Hzm[exp(jω)]和Hpn[exp(jω)]的相频响应为:
综上所述可得到:
可以看出该滤波器具有严格的线性相位特性。但是我们会发现上面得到的滤波器是不稳定的,线性相位要求该滤波器的极点在z平面上关于单位圆对称,那么就必定会有极点在单位圆之外。尽管如此,我们可以根据文献【8】提到的时间反转的方法来解决稳定性问题。
假设H1(z)的所有极点位于单位圆内,那么相应地,H2(z)的所有极点位于单位圆外。用文献【8】中提到的方法,我们构造极点都在单位圆之内的传输函数H3(z):
H3(z)用来代替H2(z),但是H3(z)的输入输出都需要做一次时间反转.这样,用输入输出都做一次时间反转的H3(z)和H1(z)级联得到的滤波器不仅具有线性相位特性而且稳定。
附图说明
图1为复系数线性相位滤波器的架构图。
图2为复系数线性相位滤波器的幅频和相频曲线。
具体实施方式
1.根据滤波器指标设计一个基本的稳定的无限脉冲响应滤波器(滤波器I),其传输函数为
2.将滤波器I的传输函数转化成零极点的形式:
3.根据发明内容,应该构造出滤波器II与滤波器I级联,滤波器II的传输函数为H3(z)=H2(z-1),H2(z)和H1(z)的零极点满足关系式
则
滤波器II的传输函数为
4.滤波过程如图1,输入数据先由滤波器I滤波,然后经过时间反转缓冲器I后由滤波器II滤波,最后再做一次时间反转后输出。
参考文献:
【1】A.Fernandez-Vazquez and G.Jovanovic-Dolecek,“Design of real and complex linearphase IIR filter banks,”IEEE Int.Symp.On Communications and Information Technology,vol.1,Oct.2005,pp.305-308.
【2】F.Argenti,V.Cappellini,A.Sciorpes and A.N.Venetsanopoulos,“Design of IIRlinear-phase QMF banks based on complex allpass sections,”IEEE Trans.Acoust.,Speech,Signal Processing,vol.44,May 1996,pp.1262-1267.
【3】J.P.Thiran,“Recursive digital filters with maximally flat group delay,”IEEE Trans.Circuit Theory,vol.18,Nov.1971,pp.659-664.
【4】G.Cortelazzo and M.Lightner,“Simultaneous design in both magnitude andgroup-delay of IIR and FIR filters based on multiple criterion optimization,”IEEE Trans.Acoust.,Speech,Signal Processing,vol.32,Oct.1984,pp.949-967.
【5】R.Mikhael,P.Agathoklis and C.Xiao,“Design of linear-phase recursive filters byoptimization of model reduced non-recursive filters,”Proc.PACRIM.2003,vol.1,Aug.2003,pp.94-97.
【6】M.A.Al-Alaoui,“Linear Phase Low-Pass IIR Digital Differentiators,”IEEE Trans.Signal Processing,vol.55,Feb.2007,pp.697-706
【7】Scott R.Powell and Paul M.Chau,“A Technique for Realizing Linear PhaseIIR Filters,”IEEE Trans.Signal Processing,vol.39,Nov.1991,pp.2425-2435.J.J.Kormylo and V.K.Jain,“Two-pass recursive digital flter with zero phase shift,”IEEETrans.Acoust.,Speech,Signal Processing,vol.22,Oct.1974,pp.384-387.
Claims (1)
1.一种复系数线性相位无限脉冲响应数字滤波器,其特征在于其传递函数如下:
其中:
M和N均为正整数,且M≤N,p1n *为复数p1n的共轭,z1m *为复数z1m的共轭。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010121871 CN101807902B (zh) | 2010-03-11 | 2010-03-11 | 复系数线性相位无限脉冲响应数字滤波器 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010121871 CN101807902B (zh) | 2010-03-11 | 2010-03-11 | 复系数线性相位无限脉冲响应数字滤波器 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101807902A true CN101807902A (zh) | 2010-08-18 |
CN101807902B CN101807902B (zh) | 2013-06-12 |
Family
ID=42609535
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010121871 Expired - Fee Related CN101807902B (zh) | 2010-03-11 | 2010-03-11 | 复系数线性相位无限脉冲响应数字滤波器 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101807902B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102082559A (zh) * | 2010-12-02 | 2011-06-01 | 广东宝莱特医用科技股份有限公司 | 线性相位iir滤波器的一种实现方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1206958A (zh) * | 1997-07-29 | 1999-02-03 | 德国汤姆逊-布朗特公司 | 有限脉冲响应滤波器的自适应方法 |
CN1452323A (zh) * | 2002-04-19 | 2003-10-29 | 松下电器产业株式会社 | 有限脉冲响应滤波器、通信发射设备和通信接收设备 |
CN1765051A (zh) * | 2004-01-30 | 2006-04-26 | 索尼株式会社 | 取样速率转换器及其方法和音频装置 |
-
2010
- 2010-03-11 CN CN 201010121871 patent/CN101807902B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1206958A (zh) * | 1997-07-29 | 1999-02-03 | 德国汤姆逊-布朗特公司 | 有限脉冲响应滤波器的自适应方法 |
CN1452323A (zh) * | 2002-04-19 | 2003-10-29 | 松下电器产业株式会社 | 有限脉冲响应滤波器、通信发射设备和通信接收设备 |
CN1765051A (zh) * | 2004-01-30 | 2006-04-26 | 索尼株式会社 | 取样速率转换器及其方法和音频装置 |
Cited By (2)
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滤波器的一种实现方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101807902B (zh) | 2013-06-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Vaidyanathan et al. | Design of doubly-complementary IIR digital filters using a single complex allpass filter, with multirate applications | |
Kormylo et al. | Two-pass recursive digital filter with zero phase shift | |
Vaidyanathan et al. | Efficient reconstruction of band-limited sequences from nonuniformly decimated versions by use of polyphase filter banks | |
WO2013060138A1 (zh) | 基于逻辑电路的法罗滤波器及其实现方法 | |
Elamaran et al. | CIC for decimation and interpolation using Xilinx system generator | |
CN103647523B (zh) | 一种降低窄过渡带fir低通滤波器实现复杂度的方法 | |
CN101807902B (zh) | 复系数线性相位无限脉冲响应数字滤波器 | |
Kumar et al. | Near perfect reconstruction quadrature mirror filter | |
Shah et al. | Efficient implementations of discrete wavelet transforms using fpgas | |
Zhou et al. | Fast design of 2D fully oversampled DFT modulated filter bank using Toeplitz-block Toeplitz matrix inversion | |
Löwenborg et al. | A class of two-channel IIR/FIR filter banks | |
CN101072018B (zh) | 数字信号的分频滤波方法及其系统 | |
EP2651033B1 (en) | Filter system | |
CN102571031B (zh) | 一种用于实现小波变换的宽带低插损声表面波滤波器组 | |
Misra et al. | Quadrature mirror filter bank with canonical signed digit representation using linear optimization algorithm | |
Elias et al. | Tree-structured IIR/FIR uniform-band and octave-band filter banks with very low-complexity analysis or synthesis filters | |
Moazzen et al. | Design of filterbanks using a fast optimization approach | |
Vaishnavi et al. | Implementation of CIC filter for DUC/DDC | |
Vandenbussche et al. | Multiplicative finite impulse response filters: implementations and applications using field programmable gate arrays | |
Chakrabarti et al. | An IFIR approach for designing M-band NPR Cosine Modulated Filter Bank with CSD | |
Živaljević et al. | Transitional Butter–Elliptic filter suitable for complementary decomposition | |
Lowenborg et al. | First-order sensitivity of complementary diplexers | |
ElGarewi et al. | Analysis of optimization algorithms for non-uniform filter bank design | |
Yang et al. | Minimax design of orthogonal filter banks with sparse coefficients | |
Paruzel et al. | Efficient fractional delay Hilbert transform filter in the Farrow structure |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130612 Termination date: 20160311 |
|
CF01 | Termination of patent right due to non-payment of annual fee |