发明名称 一种变步长LMS自适应滤波方法
摘要 本发明公开了一种变步长LMS自适应滤波方法,包括以下步骤,步骤一,采用二阶AR模型产生输入信号;步骤二,将各时刻的输入信号与其对应的二阶线性预测滤波器的二阶预测系数的乘机相加获得该时刻的输出信号;步骤三,将输入信号与输出信号做差获得误差值;步骤四,将误差值、步长值以及输入信号的乘积作为瞬时变化量更新二阶线性预测滤波器的二阶预测系数;步骤五,迭代步骤二至步骤四使得w<sub>1</sub>(n),w<sub>2</sub>(n)收敛,输出的输出信号即为所求。本发明同时考虑了指数为误差的一次和二次幂的两项,在保持收敛速度的同时,获得更优的稳态误差。
申请公布号 CN104283528A 申请公布日期 2015.01.14
申请号 CN201410478096.2 申请日期 2014.09.18
申请人 河海大学 发明人 李昌利;张师明
分类号 H03H21/00(2006.01)I 主分类号 H03H21/00(2006.01)I
代理机构 南京纵横知识产权代理有限公司 32224 代理人 董建林
主权项 一种变步长LMS自适应滤波方法,其特征在于:包括以下步骤,步骤一,采用二阶AR模型产生输入信号;输入信号公式为,x(n)+a<sub>1</sub>x(n‑1)+a<sub>2</sub>x(n‑2)=v(n)其中x(n)为n时刻的输入信号,x(n‑1)为n‑1时刻的输入信号,x(n‑2)为n‑2时刻的输入信号,a<sub>1</sub>=‑0.195,a<sub>2</sub>=0.95,v(n)为均值为零、方差为0.0965的正态分布噪声;步骤二,将各时刻的输入信号与其对应的二阶线性预测滤波器的二阶预测系数的乘积相加获得该时刻的输出信号;输出信号公式为,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>w</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mi>x</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><msub><mi>w</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mi>x</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000572437990000011.GIF" wi="1035" he="95" /></maths>其中,<img file="FDA0000572437990000012.GIF" wi="145" he="86" />为n时刻的输出信号,w<sub>1</sub>(n)和w<sub>2</sub>(n)是二阶预测系数;步骤三,将输入信号与输出信号做差获得误差值;误差值公式为,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>e</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mi>x</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000572437990000013.GIF" wi="488" he="84" /></maths>其中,e(n)表示n时刻的误差值;步骤四,将误差值、步长值以及输入信号的乘积作为瞬时变化量更新二阶线性预测滤波器的二阶预测系数;二阶线性预测滤波器的二阶预测系数公式为,<maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>w</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msub><mi>w</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>+</mo><mi>&mu;</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mi>x</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>e</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>w</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msub><mi>w</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>+</mo><mi>&mu;</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mi>x</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>2</mn><mo>)</mo></mrow><mi>e</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000572437990000014.GIF" wi="942" he="195" /></maths>其中,μ(n)为变步长,μ(n)是e(n)的Sigmoid函数;步骤五,迭代步骤二至步骤四使得w<sub>1</sub>(n)和w<sub>2</sub>(n)收敛,输出的输出信号即为所求。
地址 211100 江苏省南京市江宁区佛城西路8号
您可能感兴趣的专利