基于P范数的核最小对数绝对差自适应滤波算法

,  ,  , 

上一篇

下一篇

火元莲, 脱丽华, 齐永锋, 丁瑞博. 基于P范数的核最小对数绝对差自适应滤波算法[J]. 物理学报, 2022, 71(4): 048401-1. doi: 10.7498/aps.71.20211124
引用本文: 火元莲, 脱丽华, 齐永锋, 丁瑞博. 基于P范数的核最小对数绝对差自适应滤波算法[J]. 物理学报, 2022, 71(4): 048401-1. doi: 10.7498/aps.71.20211124
Yuan-Lian Huo, Li-Hua Tuo, Yong-Feng Qi, Rui-Bo Ding. Kernel least logarithm absolute difference algorithm based on P-norm[J]. Acta Physica Sinica, 2022, 71(4): 048401-1. doi: 10.7498/aps.71.20211124
Citation: Yuan-Lian Huo, Li-Hua Tuo, Yong-Feng Qi, Rui-Bo Ding. Kernel least logarithm absolute difference algorithm based on P-norm[J]. Acta Physica Sinica, 2022, 71(4): 048401-1. doi: 10.7498/aps.71.20211124

基于P范数的核最小对数绝对差自适应滤波算法

    通讯作者: E-mail: huoyuanlian@163.com; 
  • 中图分类号: 84.30.Vn, 05.40.Ca

Kernel least logarithm absolute difference algorithm based on P-norm

    Corresponding author: E-mail: huoyuanlian@163.com; 
  • MSC: 84.30.Vn, 05.40.Ca

  • 摘要: 为了进一步提高在α稳定分布噪声背景下非线性自适应滤波算法的收敛速度, 本文提出了一种新的基于p范数的核最小对数绝对差自适应滤波算法(kernel least logarithm absolute difference algorithm based on p-norm, P-KLLAD). 该算法结合核最小对数绝对差算法和p范数, 一方面利用最小对数绝对差准则保证了算法在α稳定分布噪声环境下良好的鲁棒性, 另一方面在误差的绝对值上添加p范数, 通过p范数和一个正常数a来控制算法的陡峭程度, 从而提高该算法的收敛速度. 在非线性系统辨识和Mackey-Glass混沌时间序列预测的仿真结果表明, 本文算法在保证鲁棒性能的同时提高了收敛速度, 并且在收敛速度和鲁棒性方面优于核最小均方误差算法、核分式低次幂算法、核最小对数绝对差算法和核最小平均p范数算法.
  • 加载中
  • 图 1  P-KLLAD在不同p下的MSE曲线

    Figure 1.  The MSE curves of P-KLLAD under different p.

    图 2  P-KLLAD在不同a下的MSE曲线

    Figure 2.  The MSE curves of P-KLLAD under different a.

    图 3  P-KLLAD在不同α下的MSE曲线

    Figure 3.  The MSE curves of P-KLLAD under different α.

    图 4  P-KLLAD在不同γ下的MSE曲线

    Figure 4.  The MSE curves of P-KLLAD under different γ.

    图 5  非线性系统辨识高斯噪声下的MSE曲线

    Figure 5.  The MSE curve of nonlinear system identification under Gaussian noise.

    图 6  非线性系统辨识α稳定分布噪声下的MSE曲线

    Figure 6.  The MSE curve of nonlinear system identification under α stably distributed noise.

    图 7  MG混沌时间序列预测高斯噪声下的MSE曲线

    Figure 7.  The MSE curve of MG chaotic time series prediction under Gaussian noise.

    图 8  P-KLLAD算法的预测值和期望值

    Figure 8.  Predicted values of the P-KLLAD algorithm and target values.

    图 9  MG混沌时间序列预测α稳定分布噪声下的MSE曲线

    Figure 9.  The MSE curve of MG chaotic time series prediction under α stably distributed noise.

    表 1  P-KLLAD算法

    Table 1.  P-KLLAD algorithm.

    P-KLLAD算法
    初始化:
      选择合适的参数a, p和核宽h,
       ${a_1}(1) = \eta \dfrac{ {ap{ {\left| {d\left( 1 \right)} \right|}^{2 p - 1} } } }{ {1{ { + } }{ {\left| {d\left( 1 \right)} \right|}^p} } }{ {\rm{sgn} } } \left( {d\left( 1 \right)} \right)$
    更新:
      每获得一对$ \left\{ {\boldsymbol u({ {i} }), d(i)} \right\} $, 计算输出值:
      $ {f_{i - 1} }(\boldsymbol u({\rm{i} })) = \eta \displaystyle \sum\limits_{j = 1}^{i - 1} { {a_j}({\rm{i} } - 1)\kappa (\boldsymbol u({\rm{i} }), \boldsymbol u(j)} ) $
    计算误差: $ e(i) = d(i) - {f_{i - 1} }(\boldsymbol u({\rm{i} })) $
    更新系数: ${a_i}(i) = \eta \dfrac{ {ap{ {\left| {e\left( i \right)} \right|}^{2 p - 1} } } }{ {1{ { + } }{ {\left| {e\left( i \right)} \right|}^p} } }{ {\rm{sgn} } } \left( {e\left( i \right)} \right)$
    下载: 导出CSV

    表 2  算法的复杂度对比

    Table 2.  Complexity of the algorithm.

    算法名称乘法或
    除法
    加法或
    减法
    指数
    运算
    符号
    函数
    分数幂
    运算
    KLMS$ 3 i - 1 $$ 2 i - 1 $$ i - 1 $——
    KFLP$ 3 i - 2 $$ 2 i - 1 $$ i - 1 $——
    KLMP$ 3 i - 2 $$ 2 i - 1 $$ i - 1 $11
    KLLAD$ 3 i + 1 $$ 2 i - 1 $$ i - 1 $——
    P-KLLAD$ 3 i + 1 $$ 2 i - 1 $$ i - 1 $12
    下载: 导出CSV

    表 3  各种算法的参数设定

    Table 3.  Parameter setting of various algorithms.

    KLMSKLMPKFLPKLLADP-KLLAD
    η0.10.10.10.10.1
    p1.20.91.1
    a53
    下载: 导出CSV

    表 4  均值 ± 标准偏差

    Table 4.  Mean ± standard deviation.

    算法MSE
    KLMS$ {{0}}{{.1952}} \pm {{0}}{{.1909}} $
    KLMP$ {{0}}{{.1921}} \pm {{0}}{{.01339}} $
    KFLP$ {{0}}{{.2311}} \pm {{0}}{{.02468}} $
    KLLAD$ {{0}}{{.1794}} \pm {{0}}{{.0209}} $
    P-KLLAD$ {{0}}{{.1792}} \pm {{0}}{{.0130}} $
    下载: 导出CSV

    表 5  均值 ± 标准偏差

    Table 5.  Mean ± standard deviation.

    算法MSE
    KLMS$ { {0} }{ {.0467} } \pm { {0} }{ {.0058} } $
    KLMP$ { {0} }{ {.0485} } \pm { {0} }{ {.0086} } $
    KFLP$ { {0} }{ {.0507} } \pm { {0} }{ {.0118} } $
    KLLAD$ { {0} }{ {.0492} } \pm { {0} }{ {.0354} } $
    P-KLLAD$ { {0} }{ {.0423} } \pm { {0} }{ {.0050} } $
    下载: 导出CSV
  • [1] Liu W F,  Príncipe J C,  Haykin S 2010 Kernel Adaptive Filtering: A Comprehensive Introduction (Hoboken, NJ, USA: John Wiley & Sons) pp16–32
    [2] 田中大, 高宪文, 石彤 2014 物理学报 63 70 doi: 10.7498/aps.63.160508 Tian Z D, Gao X W, Shi T 2014 Acta Phys. Sin. 63 70 doi: 10.7498/aps.63.160508
    [3] Zhang W, Zhu J 2020 Electronics. 9 940 doi: 10.3390/electronics9060940
    [4] Pauline S H, Samiappan D, Kumar R, Anand, A, Kar, A 2020 Applied Acoustics. 159 107074 doi: 10.1016/j.apacoust.2019.107074
    [5] 沈力华, 陈吉红, 曾志刚, 金健 2018 物理学报 67 030501 doi: 10.7498/aps.67.20171887 Shen L H, Chen J H, Ceng Z G, Jin J 2018 Acta Phys. Sin. 67 030501 doi: 10.7498/aps.67.20171887
    [6] 王世元, 史春芬, 钱国兵, 王万里 2018 物理学报 67 018401 doi: 10.7498/aps.67.20171803 Wang S Y, Shi C F, Qian G B, Wang W L 2018 Acta Phys. Sin. 67 018401 doi: 10.7498/aps.67.20171803
    [7] Wu Q, Li Y, Xue W 2019 Symmetry 11 1067 doi: 10.3390/sym11091067
    [8] Wen, F 2013 Electron. Lett. 49 1355 doi: 10.1049/el.2013.2331
    [9] 焦尚彬, 任超, 黄伟超, 梁炎明 2013 物理学报 62 210501 doi: 10.7498/aps.62.210501 Jiao S B, Ren C, Huang W C, Liang Y M 2013 Acta Phys. Sin. 62 210501 doi: 10.7498/aps.62.210501
    [10] Pelekanakis K, Chitre M 2014 IEEE Wirel. Commun. 13 3183 doi: 10.1109/TWC.2014.042314.131432
    [11] Aalo V A, Ackie A, Mukasa C 2019 Signal Process. 154 363 doi: 10.1016/j.sigpro.2018.09.025
    [12] Liu W, Pokharel P P, Principe J C 2008 IEEE Trans. Signal Process. 56 543 doi: 10.1109/TSP.2007.907881
    [13] Seo B 2011 Signal Process. 91 2623 doi: 10.1016/j.sigpro.2011.05.019
    [14] Ma W, Duan J, Man W, Zhao H, Chen B 2017 Eng. Appl. Artif. Intel. 58 101 doi: 10.1016/j.engappai.2016.11.010
    [15] Wu Z, Shi J, Xie Z, Ma W 2015 Signal Process. 117 11 doi: 10.1016/j.sigpro.2015.04.024
    [16] 赵知劲, 金明明 2017 计算机应用研究 34 3308 doi: 10.3969/j.issn.1001-3695.2017.11.024 Zhao J Z, Jin M M 2017 Appl. Res. Comp. 34 3308 doi: 10.3969/j.issn.1001-3695.2017.11.024
    [17] 董庆, 林云 2019 计算机科学 046 80 doi: 10.11896/jsjkx.181001925 Dong Q, Lin Y 2019 Comp. Sci. 046 80 doi: 10.11896/jsjkx.181001925
    [18] 火元莲, 王丹凤, 龙小强, 连培君, 齐永锋 2021 物理学报 70 415 doi: 10.7498/aps.70.20200954 Huo Y L, Wang D F, Long X Q, Lian P J, Qi Y F 2021 Acta Phys. Sin. 70 415 doi: 10.7498/aps.70.20200954
    [19] 林云, 雷洋, 曾俊俊 2016 电子技术应用 42 78 doi: 10.16157/j.issn.0258-7998.2016.02.021 Lin Y, Lei Y, Zeng J J 2016 Appl Electron. Tech. 42 78 doi: 10.16157/j.issn.0258-7998.2016.02.021
    [20] Sayin M O, Vanli N D, Kozat S S 2014 IEEE Trans. Signal Process. 62 4411 doi: 10.1109/TSP.2014.2333559
  • 加载中
图( 9) 表( 5)
计量
  • 文章访问数:  280
  • HTML全文浏览数:  280
  • PDF下载数:  0
  • 施引文献:  0
出版历程
  • 收稿日期:  2021-06-14
  • 刊出日期:  2022-02-20

基于P范数的核最小对数绝对差自适应滤波算法

    通讯作者: E-mail: huoyuanlian@163.com; 
  • 1. 西北师范大学, 物理与电子工程学院, 兰州 730000
  • 2. 西北师范大学, 计算机科学与工程学院, 兰州 730000

摘要: 为了进一步提高在α稳定分布噪声背景下非线性自适应滤波算法的收敛速度, 本文提出了一种新的基于p范数的核最小对数绝对差自适应滤波算法(kernel least logarithm absolute difference algorithm based on p-norm, P-KLLAD). 该算法结合核最小对数绝对差算法和p范数, 一方面利用最小对数绝对差准则保证了算法在α稳定分布噪声环境下良好的鲁棒性, 另一方面在误差的绝对值上添加p范数, 通过p范数和一个正常数a来控制算法的陡峭程度, 从而提高该算法的收敛速度. 在非线性系统辨识和Mackey-Glass混沌时间序列预测的仿真结果表明, 本文算法在保证鲁棒性能的同时提高了收敛速度, 并且在收敛速度和鲁棒性方面优于核最小均方误差算法、核分式低次幂算法、核最小对数绝对差算法和核最小平均p范数算法.

English Abstract

    • 核自适应滤波器[1] (kernel adaptive filtering, KAF)作为一种在重构核希尔伯特空间 (reproducing kernel Hilbert space, RKHS)中开发的高效非线性近似方法, 利用核函数[2]把原始空间的输入数据映射到RKHS空间, 从而高效解决非线性问题. 在非线性系统辨识[3]、回声消除[4]、时间序列预测[5-6]和信道均衡[7-8]等方面有很广泛的应用. 在真实的应用环境中都存在脉冲噪声和非高斯噪声, 这些噪声通常由一些自然现象和电子设备产生, 例如水下声学信号、低频大气噪声、各种类型的人为噪声, 在电力通信、控制系统、雷达探测、语音处理等实际领域的应用中较为常见. 这些噪声特性在统计意义下其概率密度分布呈现较为厚重的拖尾现象, 而α稳定分布可以对这类非高斯噪声进行很好的建模[9-11]. 普通的核最小均方(kernel least mean square, KLMS)类算法通常在高斯噪声下表现得很好[12], 但由于最小均方误差(mean square error, MSE)准则[13]只捕获误差信号的二阶统计量, 致使该类算法对异常值很敏感, 即缺乏鲁棒性. 于是一系列抗脉冲噪声干扰的核自适应滤波算法被提出. 文献[14]通过结合最小平均p范数(mean p-power error, MPE)准则[15-16]和核技巧, 提出了核最小平均p范数算法(kernel least mean p-power, KLMP)、核递归最小平均p范数算法(kernel recursive least mean p-power, KRLP), 在α稳定分布噪声环境下保证收敛速度的同时鲁棒性得到了提升. 文献[17]针对非高斯脉冲噪声环境提出了基于分数低阶统计误差准则的核分式低次幂算法(kernel kernel fractional lower power, KFLP). 文献[17]在针对文献[18]收敛速度慢的问题, 在此基础上结合Softplus函数的非线性饱和特性提出了非高斯冲激干扰下基于Softplus函数的核自适应滤波算法. 文献[19]结合线性最小对数绝对差算法[20](least logarithmic absolute difference, LLAD)和核技巧提出了核最小对数绝对差算法(kernel least logarithmic absolute difference, KLLAD), 使最小对数绝对差应用到非线性系统并且保证了其良好的抗脉冲干扰能力.

      KLLAD算法可以很好地处理异常值, 但存在收敛速度慢的问题, 为了改进KLLAD算法的收敛速度, 本文结合核最小对数绝对差算法和p范数, 提出了一种新的鲁棒核自适应滤波算法, 基于p范数的核最小对数绝对差算法(kernel least logarithm absolute difference algorithm based on p-norm, P-KLLAD). 该算法在对数的基础上添加了p范数, 通过p范数和一个正常数a来控制代价函数的陡峭程度, 从而控制算法的收敛速度, 使该算法在保证鲁棒性的基础上提升其收敛速度. 当背景噪声为α稳定分布噪声时, 在非线性系统辨识和Mackey-Glass混沌时间序列预测的仿真表明, 该算法的在收敛速度和鲁棒性方面优于核最小均方算法、核分式低次幂算法、核最小绝对差算法和核最小平均p范数算法.

    • α稳定分布和其他分布函数不同, 一般情况下没有闭式描述的概率密度函数表达式, 通常采用以下形式的特征函数来进行表示:

      其中j表示复数单位, sgn表示符号函数.

      α稳定分布主要有α, δ, β和γ四个参数共同决定.

      1) $ \alpha \in (0, 2] $: 特征指数, 它可以控制概率密度函数的拖尾的厚度. 当α取值越接近于0时, 信号或噪声出现脉冲的概率就越大. 当α为2时, 该分布退化为高斯分布.

      2) $ \delta \in R $: 位置参数, 它描述稳定分布的中心位置. 当$ \delta \in (1, 2] $时, δ表示稳定分布的均值. 当$ \delta \in $$ (0, 1] $时, δ表示稳定分布的中值.

      3) $ \beta \in \left[ { - 1, 1} \right] $: 对称系数, 它决定着稳定分布相对于中心δ的倾斜程度. 当$ \beta = 0 $时, 稳定分布是对称的. 当$ \beta \ne 0 $时, 稳定分布是非对称的.

      4) $ \gamma \in {R^ + } $: 分散系数, 它描述稳定分布相对于中心δ的扩散程度, γ越大预示着脉冲噪声强度越大.

      α稳定分布具有三种特殊的情形: 当$ \alpha =2 $, $ \beta =0 $时, α稳定分布对应于高斯分布; 当$ \alpha =1 $, $ \beta =0 $时, α稳定分布对应于柯西分布; 当$\alpha = {1}/{2}$, $ \beta ={{ - 1}} $时, α稳定分布对应于皮尔逊分布.

    • 核最小对数绝对差算法的代价函数:

      式中$ e\left( n \right) = d\left( n \right) - \boldsymbol\varOmega \left( {n - 1} \right) \boldsymbol \varphi \left( n \right) $是迭代n时刻的预测误差, $ d\left( n \right) $是观测到的期望信号, $ \boldsymbol \varOmega \left( {n - 1} \right) $是核自适应滤波器权重, $ \boldsymbol \varphi \left( n \right) $是核自适应滤波器输入. a是常数, 它决定了代价函数J的陡峭程度. 利用负随机梯度可以得到算法的权重更新公式:

      其中η表示步长, 则$ n + 1 $时刻的滤波器输出为

      根据Mercer定理, 任意核从输入空间U映射到特征空间F都可以表示为

      其中u表示输入信号在n时刻的值, $ u' $表示输入信号在下一时刻的值, (5)式被称为核技巧, 常用的核是高斯核.

      其中h为高斯核的核宽.

    • 本文提出的基于p范数的核最小对数绝对差自适应滤波算法是将核最小对数绝对差算法和p范数相结合, 构造一个新的代价函数. 在误差绝对值上添加了p范数, 并且添加了一个正常数来控制算法的陡峭程度, 使算法的收敛速度相比于核最小对数绝对差算法得到了提升, 同时利用最小对数绝对差准则保证了算法具有良好的鲁棒性. 本文算法的代价函数为

      p范数和正常数a一起决定了代价函数J的陡峭程度, 来控制收敛速度和稳态误差. 对(7)式用最速下降法求梯度得:

      利用负随机梯度可以得到算法的权重更新公式:

      其中η表示步长, 逐项递推得到如下形式:

      式中$\boldsymbol \varOmega \left( 0 \right)=0 $ , 则权重更新公式为

      本文p的取值为1.1. 从权重更新公式可以看出, 当出现脉冲异常值的时候, $ e(n) $会增大, 当$ e(n) $是小异常值时, 由于p的取值很接近于1, 那么权重更新公式会停止更新. 但是如果$ e(n) $是大异常值时, 那么分母项就会远大于分子项, 该算法就不具有很好的抗脉冲噪声效果. 综上, 滤波器$ n + 1 $时刻的输出为

      因此得到本文提出的P-KLLAD算法如表1所示.

    • 对本文提出的算法P-KLLAD进行均方收敛性能分析, 首先推导该算法的能量守恒关系, 再利用该能量守恒关系推导出最后的收敛公式. 针对非线性系统辨识, 可以得到未知系统输出为

      其中$ v\left( n \right) $是背景噪声, 则输出误差为

      将(13)式代入(14)式得到:

      其中${\tilde {\boldsymbol{\varOmega }}} \left( {n - 1} \right) = {\boldsymbol \varOmega _0} - \boldsymbol\varOmega \left( {n - 1} \right)$表示特征空间中权矢量的偏差, ${e_a}\left( n \right) = \tilde{\boldsymbol{\varOmega}} {\left( {n - 1} \right)^{\text{T}}}\boldsymbol\varphi \left( n \right)$表示第n次迭代时的先验误差. (9)式两边同时减去$ {\varOmega _0} $可得到:

      令${e_p}\left( n \right) = \tilde{\boldsymbol{\varOmega}}{\left( n \right)^{\text{T}}}\boldsymbol \varphi \left( n \right)$表示后验误差, 将(16)式代入可以表示为

      联合式(16)和式(17)消除$ e\left( n \right) $, 可得:

      对(18)式两边取内积:

      (19)式表示P-KLLAD算法的能量守恒关系. 由于${\big\| {\tilde{\boldsymbol{\varOmega}} (n)} \big\|^2}= \tilde{\boldsymbol{\varOmega}} {(n)^{\text{T}}}\boldsymbol{\tilde{\boldsymbol{\varOmega}}} (n)$, 则:

      将(17)式和(20)式代入(19)式, 可得:

      对(21)式两边取期望:

      如果算法收敛必须满足:

      所以保证该算法收敛的充分必要条件为

    • 对本文所提出的P-KLLAD算法的复杂度进行分析, 将复杂度分为乘除、加减、指数、符号函数和分数幂几种运算类型, 并与KLMS, KFLP, KLMP和KLLAD算法作比较. 所有算法在迭代次数为i时分别计算其不同运算类型的次数, 结果如表2所示. 由表2可以看出, 所有算法的运算次数都随着迭代次数i的增加而增加, 这是因为核自适应滤波算法的网络尺寸随着输入数据的增加而增大. 在第i次迭代时, 本文算法P-KLLAD相比于KLLAD算法增加了符号函数和分数幂的运算, 相比于KLMP算法增加了少量乘除和分数幂运算. 虽然P-KLLAD算法在一定程度上增加了计算复杂度, 但它在保持良好鲁棒性的同时又提高了收敛速度, 因此这并不妨碍该算法的实际应用.

    • 令$ a={{ 3}} $, 范数p分别取0.7, 0.9, 1.1, 1.3和1.5时P-KLLAD算法的学习曲线如图1所示.

      从图1可以看出, 随着p的增大, 本文算法的稳态性能和收敛速度都会变得更优, 但是当$ p \geqslant 1.3 $时该算法的抗脉冲噪声的能力会下降. 因此, 本文取p值为1.1, 这在一定程度上平衡了收敛速度和鲁棒性.

    • 令$ p=1.1 $, 正常数a分别取1, 2, 3, 4和5时P-KLLAD算法的学习曲线如图2所示.

      从图2可以看出本文算法的收敛速度随着a的增大而增大, 但稳态误差趋于不稳定. 为了平衡收敛速度和稳态误差, 本文中的a取值为3.

    • 首先讨论α稳定分布噪声的特征指数α取不同值时对算法性能的影响. 文中令$ \gamma =0.1 $, α分别取1.2, 1.5, 1.7和1.9, P-KLLAD算法的学习曲线如图3所示. 从图3可以看出, 本文的算法在不同α下都具有良好的抗脉冲干扰能力, 随着α的增大, 脉冲出现的概率也逐渐较小.

      图4给出了α稳定分布噪声的特征指数$ \alpha = $$ 1.3 $, 分散系数γ取不同值时对算法性能的影响, 令γ分别取0.1, 0.5, 1.0, 1.5和2.0. 从图4可以看出, 本文算法在不同γ下抗脉冲噪声干扰的效果也不同, 随着γ的增大噪声的脉冲会增大, 在小脉冲的时候算法具有良好的鲁棒性, 而在大脉冲的时候本文算法出现性能下降的问题.

      综上, 在后续实验中选择α稳定分布噪声的参数$ \alpha =1.3 $, $ \gamma =0.1 $来表示有脉冲干扰的噪声.

    • 为了进一步说明本文所提算法在实际应用中的有效性, 将本文所提出的算法与KLMS, KLMP, KLLAD, KLFP算法在非线性系统辨识和Mackey-Glass混沌时间序列中就收敛性和抗脉冲干扰性能进行对比. 所有实验高斯核核宽$ h = 1.0 $, 在两个仿真系统中, 每种算法均采用1000个样本作为训练数据, 100个样本用于测试, 得到的MSE性能曲线被用来评价算法性能.

    • 本实验用的非线性系统是由一个线性模型和一个非线性模型组合而成, 其中线性模型为: $ H\left( z \right) = 1 + 0.2{z^{ - 1}} $, 非线性模型为: $f\left( i \right) = x\left( i \right) - $$ 0.9 x{( i )^2}$, 则系统的期望输出: $d( i ) = x( i ) - 0.9 x{( i )^2} + $$ v( i)$. 其中$ v\left( i \right) $是环境噪声, 用α稳定分布噪声来表示. 为了使不同算法在各自最理想的情况下进行比较, 给出了不同算法的最优参数, 如表3所示.

      1) 5种算法在高斯分布噪声下的性能对比. 实验中令α稳定分布噪声参数为$ \alpha ={{ 2}} $, $ \delta =0 $, $ \beta = $$ 0 $和$ \gamma =0.1 $, 即此时该噪声是高斯分布噪声. 5种算法的MSE曲线如图5所示.

      从图5可以看出, 在高斯噪声环境下, 本文的算法和其他算法一样都有良好的稳态性能, 但是本文算法收敛速度比其他算法快, 说明本文算法具有更优的收敛性能. 其原因在于本文算法中的p范数和正常数a的主要作用就是用来提升该代价函数的陡峭程度, 从而提高算法的收敛性能.

      2)在脉冲噪声干扰下的性能比较. 实验中令α稳定分布噪声参数为$ \alpha {{ = 1}}{{.3}} $, $ \delta =0 $, $ \beta =0 $和$ \gamma =0.1 $, 此时该噪声存在脉冲干扰. 实验结果如图6所示. 表4给出了各比较算法达到稳态时候的稳态误差的均值和标准偏差.

      从图6和表4可以看出在α稳定分布噪声下本文算法和KLMP, KFLP, KLLAD都能达到比较好的稳态误差, 但是本文算法的收敛速度都比其他几种算法都快, 并且对脉冲噪声有良好的鲁棒性. 这是由于最小对数绝对差准则可以抑制脉冲噪声的无限增长.

    • Mackey-Glass混沌时间方程由以下延迟微分方程生成:

      其中, 参数设置为: $ a = 0.2 $, $ b = 0.1 $, $ {\text{c}} = 10 $. 时间序列通过一个6 s的采样周期进行离散. 根据Takens嵌入定理, 选择用最近的7个输入来预测当前的输入, 这个当前的输入就是期望信号. 这里将所有算法的参数设置为最优值, 其中KLMS和KFLP的步长分布为0.2和0.07, KLMP, KLLAD和P-KLLAD的步长均为0.12, 其他参数设置同上述实验一致.

      1)设置α稳定分布噪声参数为$ \alpha {{ = 2}} $, $ \delta =0 $, $ \beta =0 $和$ \gamma =0.1 $, 即此时该噪声是高斯分布噪声. 5种算法性能曲线比较结果如图7所示. 表5是图7中各算法达到稳态时稳态误差的均值和标准偏差. 图8给出了基于本文算法的时间序列预测中的预测值与期望值的比较结果.

      从图7可以看出, 在高斯噪声环境下, 本文的算法和KLMS算法都具有较快的收敛速度, 收敛性能优于KLMP, KFLP和KLLAD算法. 而从表5可以看出, 稳态时本文算法和KLMS算法的均值和标准偏差是最小的, 说明本文算法和KLMS都具有更优的收敛速度和稳态性能. 从图8可以看出本文算法对于MG混沌时间序列的预测值和期望值总体上基本保持一致.

      2)噪声存在脉冲异常值(α稳定分布噪声参数$ \alpha ={{ 1}}{{.3}} $, $ \delta =0 $, $ \beta =0 $和$ \gamma =0.1 $)时5种算法性能的比较. 结果如图9所示. 从图9可以看出, 存在脉冲噪声时, 本文算法比KLMS算法具有更好的稳态性能和收敛速度, 并且对脉冲噪声有良好的鲁棒性.

      综上, 本文的P-KLLAD算法虽牺牲了计算复杂度, 但在非线性系统辨识和Mackey-Glass混沌时间序列预测中不管是高斯噪声环境还是存在有小的脉冲噪声的情况下均具有良好的脉冲噪声抑制能力和较快的收敛速度.

    • 针对α稳定分布噪声模型下的非线性自适应滤波问题, 本文中提出了一种基于p范数的核最小对数绝对差自适应滤波算法P-KLLAD. 一方面利用最小对数绝对差准则来抑制α稳定分布噪声中的冲激干扰, 另一方面通过误差的p范数和一个正常数a来提高算法的收敛速度. 仿真结果表明, P-KLLAD算法在α稳定分布噪声环境下的非线性系统辨识和Mackey-Glass混沌时间序列预测中在保证鲁棒性能的同时提高了收敛速度. 但是算法在α稳定分布噪声为大异常值的时候出现性能下降问题, 有待进一步研究完善. 针对于计算复杂度问题, 后续将进一步研究有效的优化方案来降低算法复杂度.

    参考文献 (20)

目录

/

返回文章
返回