余弦调制滤波器组的稀疏FIR原型滤波器的设计方法技术

本发明专利技术公开了一种余弦调制滤波器组的稀疏线性相位FIR原型滤波器的设计方法。具体步骤如下：第1、余弦调制滤波器组的稀疏线性相位FIR原型滤波器设计参数的初始化；第2、迭代计算满足完全重建条件的余弦调制滤波器组的稀疏线性相位FIR原型滤波器，包括单位脉冲响应的非零抽头系数数目、位置以及具系数数值的确定。本发明专利技术可设计低非零抽头数的原型滤波器，滤波器的稀疏性可使其实现所用的加法器乘法器数目减少，从而能提高其运算速度、减小运算误差和降低能耗，进而降低生产成本。

Design method of sparse FIR prototype filter for cosine modulated filter banks

The invention discloses a design method of a sparse linear phase FIR prototype filter of a cosine modulated filter bank. The specific steps are as follows: first, cosine modulated filter banks with sparse linear phase FIR filter prototype parameters initialization; second, to meet the conditions of the iterative calculation of perfect reconstruction cosine modulated filter banks with sparse linear phase FIR filter prototype, including the determination of the unit impulse response of the number of non zero coefficients, and the coefficients of the position. The prototype filter of the invention can design the low non-zero number of taps, sparse filter can make the number of adders and multipliers used to achieve reduction, which can improve the operation speed and reduce the computation errors and reduce energy consumption, thereby reducing the cost of production.

【技术实现步骤摘要】
余弦调制滤波器组的稀疏FIR原型滤波器的设计方法
本专利技术属于数字信号处理
，提供了一种稀疏、高效的余弦调制滤波器组的线性相位FIR(有限脉冲响应)原型滤波器的设计方法。
技术介绍
多速率滤波器组的理论和设计因为其在通信、语音和图像编码/压缩、系统辨识、快速计算等领域有广泛的应用而备受关注。而滤波器组的构成形式一般可分为DFT滤波器组和余弦调制滤波器组。余弦调制滤波器组通过对低通原型滤波器进行优化设计，并通过快速离散余弦变换(DCT)得到分析和综合滤波器组，因其具有计算复杂度低和设计过程简单等优点，在信号处理、通信和生物医学工程等领域有广泛应用。具有稀疏的线性相位FIR(有限脉冲响应)原型滤波器的余弦调制滤波器组是每一通道的滤波器系数具有稀疏特性(非零抽头系数的数目小于滤波器阶数)的滤波器组。稀疏的滤波器其实现所用的加法和乘法器数目远少于与其滤波效果相当的同类滤波器，因此，稀疏的滤波器具有运算速度高、运算误差小和能耗低等优点。目前对于余弦调制滤波器组的设计方法主要分为分别设计滤波器组中的分析滤波器组与综合滤波器组和单独设计低通原型滤波器后经过余弦调制得到滤波器组，而对于已经提出的余弦滤波组的设计方法中，P.P.Vaidyanathan和R.D.Koilpillai的设计法可算是经典，通过分析方法，估算出满足幅度失真和混叠失真标准最小的滤波器组。
技术实现思路
本专利技术目的是设计实现小纹波，低抽头数，低幅度失真和混叠失真的余弦调制滤波器组的线性相位FIR原型滤波器，并提供一种全新的设计方法——可设计稀疏，高效的余弦调制滤波器组的线性相位FIR原型滤波器的方法。本专利技术提供的余弦调制滤波器组的稀疏线性相位FIR原型滤波器的设计方法具体步骤如下：第1、余弦调制滤波器组的稀疏线性相位FIR原型滤波器设计参数的初始化；第2、迭代计算满足完全重建条件的余弦调制滤波器组的稀疏线性相位FIR原型滤波器，包括单位脉冲响应的非零抽头系数数目、位置以及具系数数值的确定。(下面以II型线性相位FIR滤波器为例)：(一)根据设计要求构造初始参数：本专利技术根据余弦调制滤波器组的通道数M，选择通带、过渡带和阻带分别对应的采样数Lp，Lt，Ls和纹波值δp，δt，δs，确定线性相位FIR原型滤波器的初始阶数N，线性相位FIR原型滤波器的抽头系数用向量h表示为：h＝2[h1,h2,…,hm…,hN/2]T(1)其中hm(1≤m≤N/2)代表FIR原型滤波器的第m个抽头系数；将余弦调制滤波器组的稀疏线性相位FIR原型滤波器设计问题转化为如下的数学优化问题：s.t.|Bh-d|≤e(2b)其中||h||0代表0-范数运算，即表示抽头系数向量中非零抽头的个数；“min”与“s.t.”组合的公式(2a)-(2c)表示求解满足(2b)和(2c)要求的||h||0的最小值；采样矩阵B表示为B＝[Bp；Bt；Bs]，其中Bp、Bt和Bs分别代表通带、过渡带和阻带的采样矩阵，表示为其中表示维行向量；(Lp+Lt+Ls)×1维的向量d是离散化的理想频率响应，表示为：其中f(ω)是过渡带要逼近的频率响应函数，表示为：其中ω0＝απ/2M(0≤α≤1)，公式(7)满足完全重建条件，表示过渡带的频率采样点；误差向量e为(Lp+Lt+Ls)×1维的列向量，表示为：e＝[δp…δpδt…δtδs…δs]T,(8)(二)设定(Lp+Lt+Ls)×1维权重向量初始值为w(1)＝[1,1,…,1]T,本专利技术在第k(1≤k≤N/2)次迭代中，对B矩阵的列向量做归一化处理：其中利用OMP算法求解如下问题：s.t.||h(k)||0≤k(10b)并且计算(Lp+Lt+Ls)×1维的残差向量r(k)，表示为：r(k)＝Φ(k)s(k)-d(11)其中k×1维的s(k)为公式(10)的运行结果，表示OMP算法从B(k)中选出的列向量集合，集合Λ(k)＝{n1,n2,…,nk}表示非零抽头系数的指标集。(三)本专利技术利用得到的非零抽头系数的指标集Λ(k)，求解如下的线性规划问题中：s.t.|B(k)h(k)-d|≤e+μ·1L×1(11b)判断μ是否小于零，如果大于零，则更新权重向量w，其更新公式表示为：其中w(k+1)(l)表示新权重向量w(k+1)中的值，rl(k)表示残差向量r(k)中的值；将新权重向量w(k+1)带入到第2步中循环计算；如果μ小于等于零，则停止迭代运算，计算得到的即为最终的稀疏的线性相位FIR原型滤波器。本专利技术具有如下有益效果：1、本专利技术首次提供了一种稀疏、高效的余弦调制滤波器组的线性相位FIR原型滤波器设计方法。2、本专利技术可设计低非零抽头数的原型滤波器，滤波器的稀疏性可使其实现所用的加法器乘法器数目减少，从而能提高其运算速度、减小运算误差和降低能耗，进而降低生产成本。3、仿真结果表明，在相同设计指标的要求下，本专利技术的非零抽头系数的数目比国内外最佳的同类滤波器的数目少35％以上。附图说明图1是实现本专利技术的余弦调制滤波器组的稀疏线性相位FIR原型滤波器设计方法流程图；图2根据函数eam(ω)计算得出近似完全重建的余弦调制滤波器组的幅度失真图；图3是根据函数ea(ω)计算得出近似完全重建的余弦调制滤波器组的混叠失真图；图4是画出表-2中余弦调制滤波器组的稀疏线性相位FIR原型滤波器的频域响应图。具体实施方式实施例1：本专利技术提供的余弦调制滤波器组的稀疏线性相位FIR原型滤波器设计方法具体步骤如下：第1、余弦调制滤波器组的稀疏线性相位FIR原型滤波器设计参数的初始化；第2、迭代计算满足完全重建条件的余弦调制滤波器组的稀疏线性相位FIR原型滤波器，包括单位脉冲响应的非零抽头系数数目、位置以及具系数数值的确定。为了验证该滤波器组设计方法的有效性，对该方法进行了计算机模拟仿真。设计要求：利用文献：(F.Tan,etal.:“Optimaldesignofcosinemodulatedfilterbanksusingquantum-behavedparticleswarmoptimizationalgorithm,”4thInternationalCongressonImageandSignalProcessing,vol.5,pp.2280-2284,2011.)(F.Tan,etal.:“基于量子粒子群优化算法的余弦调制滤波器组优化设计,”第四界国际图像与信号处理会议,vol.5,pp.2280-2284,2011.)中所给出的设计指标，滤波器组通道数M＝16，原型滤波器初始系数个数N＝256，并设计通带、过渡带和阻带分别对应的采样数Lp＝4，Lt＝7，Ls＝94，纹波值δp＝δt＝δs＝1·10-8带入计算。本专利技术用IROMP算法设计余弦滤波器组的原型滤波器，根据权重值的计算，在得到的变量集B中，选择非零的系数位置对应的列向量Bj，带入到IROMP算法中迭代计算得到的滤波器系数。步骤一：根据余弦调制滤波器组的稀疏线性相位FIR原型滤波器的设计参数要求将各设计参数带入初始化条件，得到待求解问题：s.t.|Bh-d|≤e(2b)将设计的通带、过渡带和阻带分别对应的采样数Lp＝4，Lt＝7，Ls＝94带入公式(6)得到(4+7+94)×1维的离散化的理想频率响应d，本文档来自技高网...

【技术保护点】
一种余弦调制滤波器组的稀疏线性相位FIR原型滤波器的设计方法，其特征在于,该方法按照下述步骤进行：第1、根据余弦调制滤波器组的通道数M，选择通带、过渡带和阻带分别对应的采样数Lp，Lt，Ls和纹波值δp，δt，δs，确定线性相位FIR原型滤波器的初始阶数N，线性相位FIR原型滤波器的抽头系数用向量h表示为：h＝2[h1,h2,…,hm…,hN/2]

【技术特征摘要】
1.一种余弦调制滤波器组的稀疏线性相位FIR原型滤波器的设计方法，其特征在于,该方法按照下述步骤进行：第1、根据余弦调制滤波器组的通道数M，选择通带、过渡带和阻带分别对应的采样数Lp，Lt，Ls和纹波值δp，δt，δs，确定线性相位FIR原型滤波器的初始阶数N，线性相位FIR原型滤波器的抽头系数用向量h表示为：h＝2[h1,h2,…,hm…,hN/2]T(1)其中hm(1≤m≤N/2)代表FIR原型滤波器的第m个抽头系数；将余弦调制滤波器组的稀疏线性相位FIR原型滤波器设计问题转化为如下的数学优化问题：s.t.|Bh-d|≤e(2b)其中||h||0代表0-范数运算，即表示抽头系数向量中非零抽头的个数；“min”与“s.t.”组合的公式(2a)-(2c)表示求解满足(2b)和(2c)要求的||h||0的最小值；采样矩阵B表示为B＝[Bp；Bt；Bs]，其中Bp、Bt和Bs分别代表通带、过渡带和阻带的采样矩阵，表示为其中表示维行向量；(Lp+Lt+Ls)×1维的向量d是离散化的理想频率响应，表示为：其中f(ω)是过渡带要逼近的频率响应函数，表示过渡带的频率采样点；误差向量e为(Lp+Lt+Ls)×...

【专利技术属性】
技术研发人员：徐微，李怡，缪竟鸿，李安宇，张瑞华，
申请(专利权)人：天津工业大学，
类型：发明
国别省市：天津,12