CN110708039B - 一种farrow滤波器的系数拟合方法 - Google Patents

一种farrow滤波器的系数拟合方法 Download PDF

Info

Publication number
CN110708039B
CN110708039B CN201910874068.5A CN201910874068A CN110708039B CN 110708039 B CN110708039 B CN 110708039B CN 201910874068 A CN201910874068 A CN 201910874068A CN 110708039 B CN110708039 B CN 110708039B
Authority
CN
China
Prior art keywords
fitting
coefficients
filter
wavelet
coefficient
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.)
Active
Application number
CN201910874068.5A
Other languages
English (en)
Other versions
CN110708039A (zh
Inventor
李非
肖琨
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangxi Normal University
Original Assignee
Guangxi Normal University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Guangxi Normal University filed Critical Guangxi Normal University
Priority to CN201910874068.5A priority Critical patent/CN110708039B/zh
Publication of CN110708039A publication Critical patent/CN110708039A/zh
Application granted granted Critical
Publication of CN110708039B publication Critical patent/CN110708039B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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/0248—Filters characterised by a particular frequency response or filtering method
    • H03H17/028—Polynomial filters
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/40—Arrangements for reducing harmonics

Landscapes

  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Algebra (AREA)
  • General Physics & Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Filters That Use Time-Delay Elements (AREA)

Abstract

本发明公开了一种基于小波函数拟合的FARROW滤波器系数拟合方法，首先，对低通滤波器进行多相分解得到多相滤波器系数矩阵；接着，对该系数矩阵的每一列系数进行小波拟合，通过最小二乘法计算得到小波拟合系数；最后，将拟合得到的系数加载到FARROW滤波器结构中，输入信号通过该结构即可输出期望位置处的重采样输出值。本发明能够降低滤波器系数拟合时出现的震荡现象，与传统多项式拟合相比在相同拟合阶数情况下，改变尺度因子可以有效的提高拟合精度，提高通用解调器的计算效率。

Description

一种FARROW滤波器的系数拟合方法
技术领域
本发明涉及数字信号处理技术领域，具体涉及一种FARROW滤波器的系数拟合方法。
背景技术
适应不同调制和码率的通用解调器符合未来通信技术发展的需求，而采样率转换是实现通用解调器的关键技术之一。对于不同调制和码率的信号，其符号率也是变化的，为了提高定时误差估计和载波频偏估计的精度，通常要求输入信号的等效采样率为符号率的整数倍，因此在对信号进行此类估计之前必须对信号的等效采样率进行处理。
为了提高系统的计算效率，降低资源的消耗，通用解调器的采样率转换主要利用多相滤波结构来实现，包括利用多项式分段拟合滤波器系数。这种多项式分段拟合的缺陷在于当拟合点数多时，多项式次数高，容易产生振荡现象，从而影响拟合精度。
发明内容
本发明提供一种基于小波函数拟合的FARROW滤波器系数拟合方法，其目的在于降低滤波器系数拟合时出现的震荡现象，与传统多项式拟合相比在相同拟合阶数情况下，本发明能够有效的在通用解调器的采样率转换中克服拟合精度产生的数据误差的影响，有助于节约计算资源和提高通用解调器的计算效率。
本发明的技术方案是：首先，对低通滤波器进行多相分解得到多相滤波器系数矩阵；接着，对该系数矩阵的每一列系数进行小波拟合，通过最小二乘法计算得到小波拟合系数；最后，将拟合得到的系数加载到FARROW滤波器结构中，输入信号通过该结构即可输出期望位置处的重采样输出值。
详细技术方案如下：
(1)对低通滤波器进行多相分解得到多相滤波器系数矩阵：
对于具有KI个抽头系数的FIR滤波器，抽头系数表示成I行K列的形式：
h(0) h(I) ... h((K-1)I)
h(1) h(I+1) ... h((K-1)I+1)
... ... ... ...
h(i) h(I+i) ... h((K-1)I+i)；
h(i+1) h(I+i+1) ... h((K-1)I+i+1)
... ... ...
h(I-1) h(2I-1) ... h(KI-1)
其中,i＝0,1,…,I-1。
(2)对所述FIR滤波器系数矩阵的每一列系数进行小波拟合得到小波拟合系数：
设拟合函数为满足容许性条件的一维小波，尺度系数d∈Z，s为位移系数，U为拟合阶数，as,k为对第k列的拟合系数，对任意的k，k＝0,1,…,K-1,定义表达式：其中权系数δi＞0，U≤I；
取得极值的必要条件为其中s＝0,1,…,U；通过以上公式求解得到小波拟合系数a0,k、a1,k、…、aU,k。
(3)将拟合得到的系数加载到FARROW滤波器结构中，输入信号通过该结构即可输出期望位置处的重采样输出值：
当t＝nTx+ΔTx，0≤Δ≤1时，通过得到多相子滤波器的系数；在nTx+ΔTx时刻的输出为：/>
具体实施方式
采样率转换的问题可归纳为一个重采样的过程，先将采样信号x(nTx)重建为模拟信号：其中Tx是采样周期，在t＝mTy处对y(t)重采样得到：对信号插值的过程就是采样率转换的过程，常见的插值方法为整数倍内插，即在两个原始采样点之间插入I-1个零值,I为内插因子；然而，实际中的符号定时与采样时钟是独立的，因此采样因子不成整数倍的关系。理想内插滤波器的冲击响应由于内插滤波器是非因果的，只能无限地逼近其理想特性，为了提高计算效率和节省硬件资源，本发明采用多相滤波结构实现一种FARROW滤波器的系数拟合。
实施例的具体步骤如下：
(1)对于具有KI个抽头系数的FIR滤波器，抽头系数可以表示成I行K列的形式：
h(0) h(I) ... h((K-1)I)
h(1) h(I+1) ... h((K-1)I+1)
... ... ... ...
h(i) h(I+i) ... h((K-1)I+i)；
h(i+1) h(I+i+1) ... h((K-1)I+i+1)
... ... ...
h(I-1) h(2I-1) ... h(KI-1)
其中,i＝0,1,…,I-1。
(2)设拟合函数为满足容许性条件的一维小波，尺度系数d∈Z，s为位移系数，U为拟合阶数，as,k为对第k列的拟合系数。对任意的k，k＝0,1,…,K-1,定义表达式：其中权系数δi＞0，U≤I。于是，取得极值的必要条件为/>其中s＝0,1,…,U。由以上公式可求解得到a0,k、a1,k、…、aU,k。
(3)当t＝nTx+ΔTx，0≤Δ≤1时，通过得到多相子滤波器的系数，于是，在nTx+ΔTx时刻的输出为：/>

Claims (1)

1.一种FARROW滤波器的系数拟合方法，所述方法包括：
(1)对低通滤波器进行多相分解得到多相滤波器系数矩阵：
对于具有KI个抽头系数的FIR滤波器，抽头系数表示成I行K列的形式：
其中,i＝0,1,…,I-1；
(2)对所述FIR滤波器系数矩阵的每一列系数进行小波拟合得到小波拟合系数：
设拟合函数 为满足容许性条件的一维小波，尺度系数d∈Z，s为位移系数，U为拟合阶数，as,k为对第k列的拟合系数，对任意的k，k＝0,1,…,K-1,定义表达式：
其中权系数δi＞0，U≤I；
取得极值的必要条件为其中s＝0,1,…,U；通过以上公式求解得到小波拟合系数a0,k、a1,k、…、aU,k；
(3)将拟合得到的系数加载到FARROW滤波器结构中，输入信号通过该结构即可输出期望位置处的重采样输出值：
当t＝nTx+ΔTx，0≤Δ≤1时，通过得到多相子滤波器的系数；在nTx+ΔTx时刻的输出为：/>
CN201910874068.5A 2019-09-17 2019-09-17 一种farrow滤波器的系数拟合方法 Active CN110708039B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910874068.5A CN110708039B (zh) 2019-09-17 2019-09-17 一种farrow滤波器的系数拟合方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910874068.5A CN110708039B (zh) 2019-09-17 2019-09-17 一种farrow滤波器的系数拟合方法

Publications (2)

Publication Number Publication Date
CN110708039A CN110708039A (zh) 2020-01-17
CN110708039B true CN110708039B (zh) 2023-08-08

Family

ID=69196102

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910874068.5A Active CN110708039B (zh) 2019-09-17 2019-09-17 一种farrow滤波器的系数拟合方法

Country Status (1)

Country Link
CN (1) CN110708039B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111327297A (zh) * 2020-03-09 2020-06-23 华北电力大学 一种基于窗函数设计的自适应重采样方法
CN118316420B (zh) * 2024-06-03 2024-08-27 中国人民解放军国防科技大学 一种高精度可变时延控制方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000039952A1 (en) * 1998-12-28 2000-07-06 Motorola Inc., Method and apparatus for implementing wavelet filters in a digital system
EP2618269A1 (en) * 2012-01-17 2013-07-24 Siemens Aktiengesellschaft A method for processing of measurements from several sensors
CN108777569A (zh) * 2018-05-23 2018-11-09 成都玖锦科技有限公司 基于多相滤波器的任意延时方法
CN110146922A (zh) * 2018-08-16 2019-08-20 中铁二院工程集团有限责任公司 高速铁路地震预警系统单双地震计干扰识别方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7206359B2 (en) * 2002-03-29 2007-04-17 Scientific Research Corporation System and method for orthogonally multiplexed signal transmission and reception

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000039952A1 (en) * 1998-12-28 2000-07-06 Motorola Inc., Method and apparatus for implementing wavelet filters in a digital system
EP2618269A1 (en) * 2012-01-17 2013-07-24 Siemens Aktiengesellschaft A method for processing of measurements from several sensors
CN108777569A (zh) * 2018-05-23 2018-11-09 成都玖锦科技有限公司 基于多相滤波器的任意延时方法
CN110146922A (zh) * 2018-08-16 2019-08-20 中铁二院工程集团有限责任公司 高速铁路地震预警系统单双地震计干扰识别方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
一种巴特沃斯型FIR数字滤波器的新设计;林战平;;广西物理(第03期);全文 *

Also Published As

Publication number Publication date
CN110708039A (zh) 2020-01-17

Similar Documents

Publication Publication Date Title
CN110708039B (zh) 一种farrow滤波器的系数拟合方法
CN102414988B (zh) 用于从不均匀样本重构均匀样本的方法或结构
CN1161921C (zh) 数字解调器的符号定时恢复电路
US8564462B2 (en) Digital correction techniques for data converters
JP2014023164A (ja) 時間インタリーブａｄｃの不整合補正方法
JP2011223570A (ja) Ａｄ変換回路
CN110266311A (zh) 一种tiadc系统失配误差校准方法、装置、设备及介质
CN111194077B (zh) 一种低采样率下的定时同步方法
TWI699092B (zh) 餘數產生系統、類比數位轉換器以及產生殘餘訊號之方法
CN104393872A (zh) 一种多通道并行adc系统的采样时间误差校正方法
CN109088617B (zh) 比率可变数字重采样滤波器
CN107294534A (zh) 用于窄带信号采样的双通道tiadc频响失配实时校正方法
CN101895354B (zh) 数字域测量信号误差矢量幅度的方法
JP2011234154A (ja) アナログデジタル変換器およびその動作方法
Pei et al. Fractional bilinear transform for analog-to-digital conversion
US9954514B2 (en) Output range for interpolation architectures employing a cascaded integrator-comb (CIC) filter with a multiplier
CN107342750B (zh) 适用于多奈奎斯特区的分数延迟优化方法及其实现结构
CN116594593B (zh) 一种基于理查森外推法的数字微分器及其优化结构
US10498312B2 (en) Glitch immune cascaded integrator comb architecture for higher order signal interpolation
CN114928349A (zh) 连续时间流水线模数转换器及其数字重建滤波器
Liu et al. A novel oversampling scheme for design of hybrid filter bank based ADCs
CN114189229B (zh) 一种基于自适应分段算法的小数延时滤波方法
US20130002462A1 (en) Analog to Digital Converter with Generalized Beamformer
CN110460554B (zh) 用于全数字接收机的并行符号同步系统及方法
CN110135004B (zh) 一种宽带射频功放的动态非线性建模系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant