发明名称 一种子带自适应滤波器的权值向量更新方法
摘要 本发明公开了一种子带自适应滤波器的权值向量更新方法,属于数字滤波器设计领域。该方法首先采用子带输入信号和子带误差信号来估计子带误差信号功率和子带系统噪声功率,然后使用估计得到的子带误差信号功率和子带系统噪声功率来估计全带误差信号功率和全带系统噪声功率,最后通过比较全带误差信号功率和全带系统噪声功率的大小来确定重复利用权值向量的阶数。使用该方法的子带自适应滤波器,既可获得较快的收敛速度,还能获得较低的稳态失调。
申请公布号 CN103929150B 申请公布日期 2017.02.01
申请号 CN201410117428.4 申请日期 2014.03.27
申请人 苏州大学 发明人 倪锦根
分类号 H03H21/00(2006.01)I 主分类号 H03H21/00(2006.01)I
代理机构 苏州创元专利商标事务所有限公司 32103 代理人 范晴
主权项 一种子带自适应滤波器权值向量更新方法,其特征在于所述方法包括以下步骤:步骤1.利用子带自适应滤波器N个子带误差信号e<sub>i,D</sub>(k),i=0,1,…,N‑1,来估计N个子带误差信号功率<img file="FDA0001055325020000011.GIF" wi="158" he="78" />i=0,1,…,N‑1;所述子带误差信号功率<img file="FDA0001055325020000012.GIF" wi="157" he="78" />i=0,1,…,N‑1,由公式<img file="FDA0001055325020000013.GIF" wi="730" he="79" />来进行估计,其中λ为平滑因子,该平滑因子在[0.9,0.999]范围内取值;步骤2.利用子带自适应滤波器N个子带输入信号u<sub>i</sub>(n),i=0,1,…,N‑1,和N个子带误差信号e<sub>i,D</sub>(k),i=0,1,…,N‑1,来估计N个子带系统噪声的功率<img file="FDA0001055325020000014.GIF" wi="164" he="78" />i=0,1,…,N‑1;所述子带系统噪声功率<img file="FDA0001055325020000015.GIF" wi="163" he="78" />i=0,1,…,N‑1,由公式<img file="FDA0001055325020000016.GIF" wi="599" he="142" />来进行估计,其中r<sub>i</sub>(k)=λr<sub>i</sub>(k‑1)+(1‑λ)u<sub>i</sub>(k)e<sub>i,D</sub>(k),<img file="FDA0001055325020000017.GIF" wi="733" he="78" />且u<sub>i</sub>(k)=[u<sub>i</sub>(kN),u<sub>i</sub>(kN‑1),…,u<sub>i</sub>(kN‑M+1)]<sup>T</sup>,ξ为一个小的正常数,用来防止分母为0;步骤3.利用步骤1中得到的子带误差信号的功率<img file="FDA0001055325020000018.GIF" wi="162" he="78" />i=0,1,…,N‑1,来估计全带误差信号的功率<img file="FDA0001055325020000019.GIF" wi="147" he="63" />所述全带误差信号功率<img file="FDA00010553250200000110.GIF" wi="125" he="62" />由滤波器的所有子带误差信号功率<img file="FDA00010553250200000111.GIF" wi="162" he="78" />i=0,1,…,N‑1,合成,即<img file="FDA00010553250200000112.GIF" wi="404" he="127" />步骤4.利用步骤2中得到的子带系统噪声功率<img file="FDA00010553250200000113.GIF" wi="163" he="76" />i=0,1,…,N‑1,来估计全带系统噪声的功率<img file="FDA00010553250200000114.GIF" wi="147" he="71" />所述全带系统噪声功率<img file="FDA00010553250200000115.GIF" wi="123" he="77" />由滤波器的所有子带系统噪声功率<img file="FDA00010553250200000116.GIF" wi="163" he="79" />i=0,1,…,N‑1,合成,即<img file="FDA00010553250200000117.GIF" wi="403" he="127" />步骤5.根据步骤3中得到的全带误差信号功率<img file="FDA00010553250200000118.GIF" wi="123" he="63" />和步骤4中得到的全带系统噪声功率<img file="FDA00010553250200000119.GIF" wi="125" he="71" />来判定自适应滤波器权值向量重用的阶数R(k);所述自适应滤波器权值向量重用的阶数R(k)通过比较全带误差信号功率<img file="FDA00010553250200000120.GIF" wi="122" he="64" />与由全带系统噪声功率<img file="FDA00010553250200000121.GIF" wi="123" he="71" />表示的稳态均方误差<img file="FDA00010553250200000122.GIF" wi="395" he="70" />的大小得到,即<img file="FDA00010553250200000123.GIF" wi="653" he="151" />其中,<img file="FDA00010553250200000124.GIF" wi="230" he="61" />R<sub>max</sub>为预先选定的最大重用权值向量阶数;步骤6.用步骤5中得到的滤波器权值向量重用阶数R(k)来更新自适应滤波器的权值向量。
地址 215123 江苏省苏州市工业园区仁爱路199号