CN106130505A - 一种可变分数时延数字滤波器的优化方法 - Google Patents
一种可变分数时延数字滤波器的优化方法 Download PDFInfo
- Publication number
- CN106130505A CN106130505A CN201610459173.9A CN201610459173A CN106130505A CN 106130505 A CN106130505 A CN 106130505A CN 201610459173 A CN201610459173 A CN 201610459173A CN 106130505 A CN106130505 A CN 106130505A
- Authority
- CN
- China
- Prior art keywords
- digital filter
- frequency characteristic
- optimization method
- variable fractional
- fractional delay
- 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.)
- Pending
Links
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/0009—Time-delay networks
- H03H17/0018—Realizing a fractional delay
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/02—Frequency selective networks
- H03H17/0223—Computation saving measures; Accelerating measures
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H2017/0072—Theoretical filter design
Landscapes
- Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Mathematical Physics (AREA)
- Filters That Use Time-Delay Elements (AREA)
Abstract
本发明公开了一种可变分数时延数字滤波器的优化方法,属于信号处理与通信技术的领域。本发明采用二阶锥规划方法拟合数字延时滤波器的频率特性,根据拟合的频率特性确定群时延函数,采用梯度方法优化群时延函数,优化得到的数字滤波器同时具有最优均方误差和群时延。
Description
技术领域
本发明公开了一种可变分数时延数字滤波器的优化方法,属于信号处理与通信技术的领域。
背景技术
可变分数时延数字滤波器在各种数字信号处理领域具有广泛的应用。例如,数字接收机的时延补偿,COMB数字滤波器设计,采样率转换,语音编码,时延估计,数字信号插值等等。事实上,在那些频率特性需要可调的数字信号处理领域中都有广泛的应用。
在数字通信***中,***的整数时延可以通过简单的移位寄存器完成,而对非整数时延,可以先通过前述方法完成整数时延,而对分数部分需要设计专门的滤波器,这就提出了可变分数时延(VFD)滤波器。可变分数时延滤波器旨在通过拟合设计出具有如下频率特性的滤波器:
Hd(ω,p)=e-jωp,ω∈[0,2π],p∈[0,1]
拟合的数字滤波器结构为:
要求拟合出的滤波器频率特性H(ω,p)与目标滤波器频率特性Hd(ω,p)尽可能相等。国内外的研究人员一般通过系数关系优化设计
minimaxW(ω,p)|H(ω,p)-Hd(ω,p)|
但这类方法得到的滤波器具有较大的峰值时延偏移。也有部分方法试图优化
但直接优化以上目标函数存在非线性问题,因而计算效率低下、效果差。还有部分方法通过融合以上两个目标函数,这样可以同时兼顾峰值偏移和均方误差,但同样存在非线性问题。
发明内容
本发明所要解决的技术问题是针对上述背景技术的不足,提供了一种可变分数时延数字滤波器的优化方法,优化得到的数字滤波器具有最优均方误差和群时延,解决了现有技术不能同时优化均方误差和群时延的技术问题。
本发明为实现上述发明目的采用如下技术方案:
一种可变分数时延数字滤波器的优化方法,包括如下步骤:
采用二阶锥规划方法拟合数字延时滤波器的频率特性;
根据拟合的频率特性确定群时延函数,采用梯度方法优化群时延函数。
进一步的,所述可变分数时延数字滤波器的优化方法中,拟合数字延时滤波器频率特性的目标函数为:minimax W(ω,p)|H(ω,p)-Hd(ω,p)|,其中,H(ω,p)为拟合的频率特性,Hd(ω,p)为目标频率特性,W(ω,p)是权重因子,ω为频率,ω∈[0,2π],p为随机数,p∈[0,1]。
作为所述可变分数时延数字滤波器优化方法的进一步优化方案,Hd(ω,p)=e-jωp,amn为系数,M、N均为正整数。
再进一步的,所述可变分数时延数字滤波器的优化方法中,根据拟合的频率特性确定群时延函数τ(ω,p),具体为:
更进一步的,所述可变分数时延数字滤波器的优化方法中,优化群时延函数的目标函数为:
本发明采用上述技术方案,具有以下有益效果:
(1)采用本发明技术方案优化得到的数字滤波器同时具有最优均方误差和群时延;
(2)采用二阶锥规划的作用是将频率特性优化这一非线性问题转化为线性问题,再以二阶锥规划得到的线性解做为群时延优化的初始值,自然而然地把群时延优化问题转化为线性问题,提高了运算效率。
具体实施方式
下面对发明的技术方案进行详细说明。
依序进行两个步骤的优化设计,步骤二以步骤一的结果作为初始值,步骤二的结果为最终的优化结果。步骤一采用SOCP方法进行,可以采用广泛使用的共享软件包SeDuMi;步骤二采用梯度计算方法进行。
步骤一:采用二阶锥规划方法计算:minimax|H(ω,p)-Hd(ω,p)|,以拟合数字延时滤波器的频率特性,Hd(ω,p)=e-jωp,ω∈[0,2π],p∈[0,1];H(ω,p)为拟合的频率特性,Hd(ω,p)为目标频率特性,W(ω,p)是权重因子,ω为频率,ω∈[0,2π],p为随机数,p∈[0,1],amn为系数,M、N均为正整数。
步骤二:根据拟合的频率特性确定群时延函数τ(ω,p),
采用梯度方法计算:以优化群时延函数。
二阶锥规划方法将频率特性优化这一非线性问题转化为线性问题,再以二阶锥规划得到的线性解做为群时延优化的初始值,自然而然地把群时延优化问题转化为线性问题,使得优化后的数字滤波器同时具有最优均方误差和群时延,提高了运算效率。
Claims (5)
1.一种可变分数时延数字滤波器的优化方法,其特征在于,包括如下步骤:
采用二阶锥规划方法拟合数字延时滤波器的频率特性;
根据拟合的频率特性确定群时延函数,采用梯度方法优化群时延函数。
2.根据权利要求1所述一种可变分数时延数字滤波器的优化方法,其特征在于,拟合数字延时滤波器频率特性的目标函数为:minimax W(ω,p)|H(ω,p)-Hd(ω,p)|,其中,H(ω,p)为拟合的频率特性,Hd(ω,p)为目标频率特性,W(ω,p)是权重因子,ω为频率,ω∈[0,2π],p为随机数,p∈[0,1]。
3.根据权利要求2所述一种可变分数时延数字滤波器的优化方法,其特征在于,Hd(ω,p)=e-jωp,amn为系数,M、N均为正整数。
4.根据权利要求2或3所述一种可变分数时延数字滤波器的优化方法,其特征在于,根据拟合的频率特性确定群时延函数τ(ω,p),具体为:
5.根据权利要求4所述一种可变分数时延数字滤波器的优化方法,其特征在于,优化群时延函数的目标函数为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610459173.9A CN106130505A (zh) | 2016-06-22 | 2016-06-22 | 一种可变分数时延数字滤波器的优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610459173.9A CN106130505A (zh) | 2016-06-22 | 2016-06-22 | 一种可变分数时延数字滤波器的优化方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106130505A true CN106130505A (zh) | 2016-11-16 |
Family
ID=57269184
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610459173.9A Pending CN106130505A (zh) | 2016-06-22 | 2016-06-22 | 一种可变分数时延数字滤波器的优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106130505A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108051785A (zh) * | 2017-11-24 | 2018-05-18 | 电子科技大学 | 宽带数字阵列雷达接收通道的优化设计方法 |
CN113037248A (zh) * | 2020-12-04 | 2021-06-25 | 北京理工大学 | 一种基于分段频域优化的小数延时滤波器设计方法 |
CN113607446A (zh) * | 2021-05-20 | 2021-11-05 | 西安交通大学 | 一种机械设备早期故障诊断方法、***、设备及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101425183A (zh) * | 2008-11-13 | 2009-05-06 | 上海交通大学 | 基于二阶锥规划的变形体三维跟踪方法 |
CN102426562A (zh) * | 2011-08-15 | 2012-04-25 | 天津大学 | 一种基于支持向量机的核矩阵近似方法 |
US20150333732A1 (en) * | 2014-05-13 | 2015-11-19 | Hong Kong Applied Science And Technology Research Institute Co., Ltd. | Compensation filter for cascaded-integrator-comb decimator |
-
2016
- 2016-06-22 CN CN201610459173.9A patent/CN106130505A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101425183A (zh) * | 2008-11-13 | 2009-05-06 | 上海交通大学 | 基于二阶锥规划的变形体三维跟踪方法 |
CN102426562A (zh) * | 2011-08-15 | 2012-04-25 | 天津大学 | 一种基于支持向量机的核矩阵近似方法 |
US20150333732A1 (en) * | 2014-05-13 | 2015-11-19 | Hong Kong Applied Science And Technology Research Institute Co., Ltd. | Compensation filter for cascaded-integrator-comb decimator |
Non-Patent Citations (1)
Title |
---|
WU-SHENG LU: "A Second-Order Cone Programming Approach for Minimax Design of 2-D FIR Filters with Low Group Delay", 《IEEE INTERNATIONAL SYMPOSIUM ON CIRCUITS AND SYSTEMS》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108051785A (zh) * | 2017-11-24 | 2018-05-18 | 电子科技大学 | 宽带数字阵列雷达接收通道的优化设计方法 |
CN113037248A (zh) * | 2020-12-04 | 2021-06-25 | 北京理工大学 | 一种基于分段频域优化的小数延时滤波器设计方法 |
CN113037248B (zh) * | 2020-12-04 | 2022-11-11 | 北京理工大学 | 一种基于分段频域优化的小数延时滤波器设计方法 |
CN113607446A (zh) * | 2021-05-20 | 2021-11-05 | 西安交通大学 | 一种机械设备早期故障诊断方法、***、设备及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101295969B (zh) | 高阶有限冲击响应数字滤波器设计方法 | |
CN106130505A (zh) | 一种可变分数时延数字滤波器的优化方法 | |
CN103152012B (zh) | 一种降低通带群延迟误差的fir滤波器设计方法 | |
CN105137757A (zh) | 一种具有频率自适应能力的重复控制器及控制方法 | |
CN107294511B (zh) | 一种低复杂度的可变分数时延滤波方法及滤波器 | |
US20100185450A1 (en) | Method and system for efficient optimization of audio sampling rate conversion | |
CN108092644A (zh) | 一种陷波频率精准可调的稀疏二维fir陷波滤波器的设计方法 | |
CN105159062A (zh) | 一种基于***式快速重复控制器的复合控制方法及*** | |
CN103888104A (zh) | Fir数字滤波器设计方法和*** | |
CN103346756A (zh) | 一种空时自适应滤波方法及装置 | |
CN102435849A (zh) | 使用数字滤波器提高频率测量精度的方法 | |
CN103956987B (zh) | 一种数字滤波器设计方法 | |
KR101943529B1 (ko) | 오디오 신호를 처리하기 위한 방법 및 장치 | |
CN106982045A (zh) | 一种基于socp优化的eir‑cmfb结构的设计方法 | |
CN107342750A (zh) | 适用于多奈奎斯特区的分数延迟优化方法及其实现结构 | |
Ferdi | Impulse invariance-based method for the computation of fractional integral of order 0< α< 1 | |
CN102411301A (zh) | 控制对象频域相位特性辩识方法和装置 | |
CN105337588A (zh) | 一种组合积分滤波*** | |
Devate et al. | Wideband IIR digital integrator and differentiator using curve fitting technique | |
CN106026974A (zh) | 阻带响应加权通带总体响应误差约束空域矩阵滤波器设计方法 | |
Li et al. | Compensation method for the CIC filter in digital down converter | |
CN107453348A (zh) | 多电飞机电网频率估计方法及装置、航空电气***、飞机 | |
CN103634027B (zh) | 超宽带信号的数字正交调制实时处理方法 | |
CN107786184A (zh) | 音频频谱转换处理方法 | |
Chan et al. | A new method for least-squares and minimax group-delay error design of allpass variable fractional-delay digital filters |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20161116 |
|
WD01 | Invention patent application deemed withdrawn after publication |