专利名称：：一种基于加权次梯度投影的数字助听器回声路径估计方法
技术领域：
：本发明涉及一种数字助听器回声路径自适应估计方法。二.
背景技术：
：回声是困扰数字助听器使用的一个普遍问题。通常回声消除的主要思想是先估计回声路径，然后根据回声路径估计回声信号，并将其从输入信号中减去。自适应滤波器法是最常见的回声消除方法。这种方案的核心是设计自适应滤波器估计回声路径，进而估计出回声信号，并将其从输入信号中减去。因为回声路径随环境而变化，故自适应滤波器应能快速跟踪这种变化。数字助听器回声路径的自适应估计算法有很多种。在本发明之前，最常用的是是LMS(LeastMeanSquare,最小均方误差)算法、NLMS(NormalLeastSquare,归一化最小均方误差)算法和RLS(RecursiveLeastSquare,迭代最小二乘)算法。2000年J.M.Kates在U.S.patentUS6072884(Feedbackcancellationapparatusandmethods)中介绍了LMS算法和NLMS算法在数字助听器回声消除系统的应用，这类算法通过最小化均方误差目标函数的方法估计回声路径，主要优点是计算量较小，有利于实时实现，但是算法的收敛速度和收敛精度都不高，并且在低信噪比和高有色噪声的环境下性能明显下降。S.Haykin在AdaptiveFilterTheory.(UpperSaddleRiver,NJ:Prentice-Hall,1996)—书中介绍了RLS自适应估计算法，该算法虽然能够取得较高的收敛性能，但计算复杂性随滤波器阶数快速上升，难以在数字助听器中实时实现。近年来AP(AffineProjection,仿射投影)算法和POCS(ProjectionofOntoConvexSets,凸集投影)算法等投影算法也被应用于通信系统回声路径估计，在这类应用中，AP算法改善了经典NLMS算法在有色信号输入时的收敛性能，而POCS算法相比AP算法的优点在于POCS利用投影算子向特定的闭合凸集多次投影，在低信噪比情况下取得了比AP算法更快的收敛速度。但是，这类投影算法的缺点在于它们所使用的投影算子常常难以定义，而且也没有在数字助听器回声消除系统中应用的研究。三.
发明内容1.所要解决的技术问题本发明的目的就在于克服上述缺陷，针对数字助听器回声消除模型，设计实现了一种新的加权次梯度投影回声路径估计方法。2.技术方案一种数字助听器回声路径自适应估计方法，其主要技术步骤如下(1)初始化自适应滤波器h。，该滤波器系数在算法中将迭代更新，收敛于实际反馈路径h.;(2)构造数字助听器输出信号矩阵U,=[ut,ut—,，…,L+,r，这里k表示离散信号时间序号，r为算法的步长，114=["(""("1)，...,"("^+1)广为助听器输出语音信号向量，N为自适应滤波器长度；(3)获取数字助听器输入信号向量dr-Uh'+n广[rf(A),^fc-1),…^(""l)r，其中h'表示实际反馈路径，！^表示助听器中原始输入信号S^和叠加噪声的总和；(7)(4)计算凸函数g(10—llfh广dU-P，其中参数p使用公式p:-m广mH十算，这里CT为I^的方差；(5)计算次梯度t=Vg(h')=(U:h'-dt);(6)根据已有的数字助听器系统先验信息设定指数衰减权重矩阵A=A......00......aw=BB，这里A和B满秩，N为自适应滤波器阶数,、l,……y为指数衰减因子，满足0<7<1，y的值通过先验知识获得或经过测量得到；计算加权次梯度投影<formula>formulaseeoriginaldocumentpage4</formula>(式l)(8)更新滤波器系数h",ht+^(尸矿(w(h,)—ht),^e、根据图l所示的模型，该向量d,:-U〖h'+i^，其中h'表示实际反馈路径，1^表示助听器中原始输入信号S^和叠加噪声的总和。(4)为了在第k个时刻通过Updt，使自适应滤波器系数h向实际反馈路径h'逼近，采取加权次梯度投影的逼近方法，其公式推导如下。定义凸集Cmin^lu:h-d*|,(式3){;"1"..，附}满足如下约束Co):={"W:|u:h_dA『-ps0}f2、'(式4)={/z":|u〖h-d*||、/>}这样对实际反馈路径h'的估计问题就转变成了向凸集C的投影问题。定义凸函数g(h):=|urh-d『-AVheY(式5)以及梯度算子t:=Vg(h)=2U(Urh-d),Vhe^,(式6)那么可以定义基于该凸函数的半空间/T(h4):={he;W:(h-hjs+g(hj,，(式7)当VxeC;(p)有g(x)SO，并且可以得到(x—y)rs+g(y)5g(X)S0,(式8)(式8)意味着xe/T(y)，因此C^(A)c/r(hJ，即对CJp)的投影可以扩展为对半空间/T(^)的投影，此为次梯度投影，其投影公式为P/r(h4)(h*):(式9)参数p的取值直接影响了lTeC(P)的概率。较大的p的取值意味着约束集合包含真实反馈路径概率较大，但收敛速度较慢；相反，较小的p的取值意味着收敛速度较快，但约束集合包含真实反馈路径的概率较小，可能导致自适应系统收敛到错误的反馈路径。因此p的取值应与相应的噪声环境相吻合。假定随机变量《=||^||2满足,2p.d.f.,可以使用简单的公式yo:=附《=ro"2计算/0，这里CT表示i^的方差。由于在助听器回声路径估计时泄漏在耳道内的脉冲响应是遵循指数衰减规律的，因此充分考虑这样的衰减特性可以提高算法的收敛速度和收敛精度。下面介绍利用回声路径的先验信息设定权值矩阵，使自适应滤波器的系数按照权值矩阵的规律收敛到理想路径的加权次梯度投影方法。这里设定指数衰减权重矩阵A:0:BB，A和B满秩，N为自适应滤波器阶数，a,.、/",..….r为指数衰减因子，满足0<7<1,y的值通过先验知识获得或经过测量得到。当考虑投影的权重矩阵A后，定义ii-B—^禾卩tJ:BU，则构造凸函数g(fi)定义为g(ii):=|tjrii—d|『—p'(式io)并且注意到禍H卜,b-'h—4—P，(式u)=|urh-d『-P=g(h)因此凸集(^(^):=5;;￡//:8(￡)^0}实际上和(式4)定义的凸集完全相同。重新定义半空间g-(^)使其能够分离^=B-、和凸集&Go)，根据次梯度投影原理可以得到t):={":(fi-吣)r々g(fik)+g(、),，这里fg(￡k)表示在Kk点凸函数对￡的次梯度，并且々g(ii)=2帥q-d),Vfie:W。这里注意到禾口々g(fi)=2BU(lfBH-=BVg(h)A-(ht):={hW.(B-1h-B、)rBVg(hk)+g(h*)《0}={hKh-ht)r(B-1)rBVg(hk)+g(ht)S0}={lH(h-)rVg(hk)+g(ht),Kh4)(式12)(式13)(式14)(式15)7根据(式9)从iik向半空间5—(fi4)的投影可以定义为h4+(式16)根据定义fi-B-ih禾口6二BU,我们可以将(式16)B'h*UB—、)_利用(式ii)、(式M)和(式i6)，(式n)等价于写为g(B-、)S0g(B、)〉0'(式17)即、+—BBg(ht)(V&(ht))'BBV&(hJ▽&(h4)g(hj>0一她)g(、)>0(式18)(式19)(V&(ht))'AV&(hJ(式19)定义了加权次梯度投影公式，表示为WAPSM(WeightedAdaptiveProjectionSubgradientMethod,加权次梯度投影算法)。它意味着h,的更新方向由设定的权矩阵A修正，该投影的几何意义解释(简单起见，设M-2)由图2所示。本方法根据(式19)计算滤波器系数、的加权次梯度投影。首先计算凸函数g(ht):=|u、-d|f-P，其中参数/7使用公式p:-^:mH十算，这里tr为i^的方差；然后计算次梯度t-Vg(hJ-2(U:h'-dt);最后根据(式19)计算计算加权次梯投影/V,h、(ht)。(5)更新滤波器系数:h*+W-(h,>(h*)—h"4e.(式20)(6)*=/t+l,返回步骤(2)，直至信号序列处理完毕。3.性能评价分析以下的实验比较了本发明提出的WAPSM(WeightedAdaptiveProjectionSubgradientMethod,加权次梯度投影)算法与NLMS算法和APSM算法的性能差异。实验使用了两种回声路径，第一种回声路径hl(图3)为空间尺寸5mX4mX3m的真实房间作为回声路径，第二种回声路径gl(图4)由下式给出W)e掘"1""200，(式21)这里^w)为均值0方差0.5的高斯白噪声，之所以选用这种路径是因为它非常类似于助听器扬声器经耳道和耳罩泄露到麦克风的真实测量回声路径。实验考虑两种类型的助听器输入信号，第一种使用高斯白噪声作为输入信号，第二种使用TIMIT数据库中的"anl01-mtms-senn"语音段作为输入信号，该语音片段采样率为16KHz。实验系统的前向路径采用固定增益G。，并具有单一延迟单元，其传递函数为G(z)=G。z-'。实验针对G。=3,5,10多种取值进行，不考虑自动增益控制(AGC)单元。在实验中NLMS算法取归一化收敛因子//=0.5，在hi路径估计时^=0.99，在gl路径估计时y=0.95123。回声路径估计方法的性能由自适应滤波器失准系数(Misalignment)，其定义为附/》="L;;~这里《和d)为真实反馈路径的参数及估计路径的参数，为范数算子。比较WAPSM算法和NLMS算法以及APSM算法在固定增益G。情况下的失准性能。实际结果如图5所示。图中实验结果显示，不论是白噪声作为输入信号还是语音片段作为输入信号，WAPSM算法都取得了比NLMS算法和APSM算法更稳定的收敛性能和更小的失准系数。对hl反馈路径，由于路径的指数衰减并不明显，因此WAPSM算法和APSM算法性能相似，二者在白噪声输入的情况下，都取得了比NLMS算法低15dB的最终失准系数；对语音片段输入，二者比NLMS算法失准系数降低了6dB。对gl反馈路径，由于路径有明显的指数衰减特征，WAPSM算法比APSM和NLMS算法的性能更优越。在白噪声输入的情况下，WAPSM算法的失准系数比APSM算法下降3dB，比NLMS算法下降4dB;在语音片段作为输入的情况下，WAPSM算法的失准系数比APSM下降了ldB，比NLMS算法下降了5dB。另外，WAPSM算法的失准系数下降的收敛速度和收敛稳定性得到显著提升。特别从图6(c)和图6(d)可以看出，在实际语音片段作为助听器输入信号的情况下，实验结果显示WAPSM算法取得了卓越的失准性能，而绝大多数常用算法，如NLMS算法，在实际语音信号作为输入时性能都有显著下降。权利要求1.一种数字助听器回声路径自适应估计方法，其特征是使用加权次梯度投影算法估计回声路径，其主要步骤为(1)初始化自适应滤波器h0，该滤波器系数在算法中将迭代更新，收敛于实际反馈路径h*；(2)构造数字助听器输出信号矩阵Uk＝[uk，uk-1，...，uk-r+1]T∈RN×r，这里k表示离散信号时间序号，r为算法的步长，uk＝[u(k)，u(k-1)，...，u(k-N+1)]T为助听器输出语音信号向量，N为自适应滤波器长度；(3)获取数字助听器输入信号向量<mathsid="math0001"num="0001"><math><![CDATA[<mrow><msub><mi>d</mi><mi>k</mi></msub><mo>:</mo><mo>=</mo><msubsup><mi>U</mi><mi>k</mi><mi>T</mi></msubsup><msup><mi>h</mi><mo>*</mo></msup><mo>+</mo><msub><mi>n</mi><mi>k</mi></msub><mo>=</mo><msup><mrow><mo>[</mo><mi>d</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo><mi>d</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>d</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mi>r</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>]</mo></mrow><mi>T</mi></msup><mo>,</mo></mrow>]]></math>id="icf0001"file="A2008100240450002C1.tif"wi="72"he="4"top="91"left="88"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>其中h*表示实际反馈路径，nk表示助听器中原始输入信号sk和叠加噪声的总和；(4)计算凸函数g(hk)＝||UThk-d||2-ρ，其中参数ρ使用公式ρ＝mξ＝rσ2计算，这里σ为nk的方差；(5)计算次梯度<mathsid="math0002"num="0002"><math><![CDATA[<mrow><mi>t</mi><mo>=</mo><mo>▿</mo><mi>g</mi><mrow><mo>(</mo><msub><mi>h</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mrow><mn>2</mn><mi>U</mi></mrow><mi>k</mi></msub><mrow><mo>(</mo><msubsup><mi>U</mi><mi>k</mi><mi>T</mi></msubsup><msub><mi>h</mi><mi>k</mi></msub><mo>-</mo><msub><mi>d</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>;</mo></mrow>]]></math>id="icf0002"file="A2008100240450002C3.tif"wi="42"he="4"top="133"left="58"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>(6)根据已有的数字助听器系统先验信息设定指数衰减权重矩阵id="icf0003"file="A2008100240450002C4.tif"wi="46"he="20"top="151"left="39"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>这里A和B满秩，N为自适应滤波器阶数，αi＝α0γi-1，i＝1，......，N，γ为指数衰减因子，满足0＜γ＜1，γ的值通过先验知识获得或经过测量得到；(7)计算加权次梯度投影<mathsid="math0003"num="0003"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><msup><mi>H</mi><mo>-</mo></msup><mrow><mo>(</mo><msub><mi>h</mi><mi>k</mi></msub><mo>)</mo></mrow></mrow></msub><mrow><mo>(</mo><msub><mi>h</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>=</mo><mfencedopen='{'close=''><mtable><mtr><mtd><msub><mi>h</mi><mi>k</mi></msub><mo>,</mo></mtd><mtd><mi>g</mi><mrow><mo>(</mo><msub><mi>h</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>≤</mo><mn>0</mn></mtd></mtr><mtr><mtd><msub><mi>h</mi><mi>k</mi></msub><mo>-</mo><mfrac><mrow><mi>Ag</mi><mrow><mo>(</mo><msub><mi>h</mi><mi>k</mi></msub><mo>)</mo></mrow></mrow><mrow><msup><mrow><mo>(</mo><mo>▿</mo><mi>g</mi><mrow><mo>(</mo><msub><mi>h</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mi>T</mi></msup><mi>A</mi><mo>▿</mo><mi>g</mi><mrow><mo>(</mo><msub><mi>h</mi><mi>k</mi></msub><mo>)</mo></mrow></mrow></mfrac><mo>▿</mo><mi>g</mi><mrow><mo>(</mo><msub><mi>h</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>,</mo></mtd><mtd><mi>g</mi><mrow><mo>(</mo><msub><mi>h</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>></mo><mn>0</mn></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math>id="icf0004"file="A2008100240450002C5.tif"wi="82"he="17"top="206"left="51"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>(式1)(8)更新滤波器系数<mathsid="math0004"num="0004"><math><![CDATA[<mrow><msub><mi>h</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mi>h</mi><mi>k</mi></msub><mo>+</mo><msub><mi>λ</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mi>P</mi><mrow><msup><mi>H</mi><mo>-</mo></msup><mrow><mo>(</mo><msub><mi>h</mi><mi>k</mi></msub><mo>)</mo></mrow></mrow></msub><mrow><mo>(</mo><msub><mi>h</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>h</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>,</mo><msub><mi>λ</mi><mi>k</mi></msub><mo>∈</mo><mo>[</mo><mn>0,2</mn><mo>]</mo><mo>;</mo></mrow>]]></math>id="icf0005"file="A2008100240450002C6.tif"wi="62"he="5"top="236"left="62"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>(式2)(9)k＝k+1，返回步骤(2)，直至信号序列处理完毕。全文摘要本发明提出了一种用于数字助听器回声路径估计的加权次梯度投影自适应系统估计方法。本发明通过推导得出加权次梯度投影公式，并将其用于数字助听器自适应回声消除系统之中。其具体方法是利用数字助听器连续输出信号序列和输入信号序列，通过加权次梯度投影算法使自适应滤波器系数h逼近实际回声路径h<sup>*</sup>。当算法的权值矩阵设置合理的情况下，算法的收敛速度和精度相比较以前的算法有很大提高。并且，无论在白噪声环境还是有色噪声环境下，算法都取得了卓越的失准性能。文档编号H04R25/00GK101568058SQ20081002404公开日2009年10月28日申请日期2008年4月25日优先权日2008年4月25日发明者杰乔,华余,吉奚,王青云,力赵,邹采荣申请人:王青云