[0001] 本发明涉及滤波器领域，具体涉及一种鲁棒性自适应滤波系统。

[0002] 传统的滤波器是由电容、电感和电阻组成的滤波电路。滤波器可以对电源线中特定频率的频点或该频点以外的频率进行有效滤除，得到一个特定频率的电源信号，或消除一个特定频率后的电源信号。

[0003] 但其不能实现自适应滤波，本发明提供一种鲁棒性自适应滤波系统能够实现自适应滤波。

[0035] 为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。

[0036] 实施例1

[0037] 本实施例提供一种鲁棒性自适应滤波系统，如图1，所述鲁棒性自适应滤波系统为时分复用滤波，包括数字滤波器和广义滤波权值预测控制单元；

[0038] 数字滤波器为FIR滤波器，包含路数为M的多路复用器，多路复用器的输入连接RAM存储器、采样数据输入口，输出连接单路FIR滤波器单元，单路FIR滤波器单元的输入为第m路的输入xm(n)，输出为第m‑1路采样数据xm‑1(n)和第m路滤波结果ym(n)；

[0039] m为小于等于M的正整数，N为FIR横向滤波器阶数，i为小于N的整数；

[0040] 所述广义滤波权值预测控制单元执行下列步骤，完成对数字滤波器的输出估计：

[0041] 步骤1，定义第m路的滤波器单元滤波权值为Wm(n)，定义第m单路FIR滤波器单元的输出ym(n)与目标信号dm(n)之差em(n)的均方值为均方差误差

[0042]

[0043]

[0044] 其中，R＝E[Xm(n)XTm(n)]，P＝E[dm(n)Xm(n)]；

[0045] 步骤2，采用最小均方差方法计算出均方差误差最小时的Wm(n)作为最优滤波权值；

[0046] 步骤3，定义最优滤波权值调整估计时长j，在时长j内预估并控制调整滤波权值至最优滤波权值。

[0047] 优选地，所述步骤3还包括：

[0048] 步骤3.1，定义干扰估计时长p，定义干扰估计系数λ1，最优滤波权值调整系数为λ2，第m单路FIR滤波器单元的滤波权值输出估计模型为：

[0049]

[0050]

[0051] 1＝E1，j(q‑1)(1+M1q‑1)Δ+q‑jF1，j(q‑1)；

[0052] 1＝E2，j(q‑1)(1+M2q‑1)Δ+q‑jF2，p(q‑1)；

[0053] 其中，w(k)为第k时刻滤波权值，u(k)第k时刻调整滤波权值的输入值，z(k)干扰‑1 ‑1 ‑1 ‑1 ‑1 ‑1值，M1(q )为预定义的q 多项式，N1(q )为预定义的q 多项式，M2(q )为预定义的q 多项‑1 ‑1 ‑1
式，N2(q )为预定义的q 多项式；C1＝1，C2＝1，k为时刻值，q 为预定义的后移系数，Δ＝1‑‑1
q 为差分系数，ε1(k)和ε2(k)为随机干扰项因子；

[0054] 步骤3.2，采集数据w(k)，z(k)，w(k‑1)，z(k‑1)，Δu(k‑1)...Δu(k‑20)；Δu(k‑1)...Δu(k‑c2)为过去调整滤波权值的输入值的差分量；

[0055] 步骤3.3，定义调整滤波权值的输入值的差分量

[0056]

[0057] 其中，b1为最优滤波权值调整系数，b2干扰模型系数；

[0058] 定义计算出

[0059] 步骤3.4，根据u(k)＝u(k‑1)+Δu(k)计算出控制量；

[0060] 步骤3.5，定义k＝k+1，返回步骤3.2。

[0061] 具体地，定义干扰估计时长p为两倍最优滤波权值调整估计时长j。干扰一般不会过于频繁的变化，这能够减小对于干扰估计带来的开销增加。

[0062] 具体地，最优滤波权值调整估计时长j为一个滤波周期。通过时序的衔接，能够兼顾鲁棒性、开销。

[0063] 具体地，所述FIR滤波器为FIR横向滤波器。

[0064] 尽管上面对本发明说明性的具体实施方式进行了描述，以便于本技术领域的技术人员能够理解本发明，但是本发明不仅限于具体实施方式的范围，对本技术领域的普通技术人员而言，只要各种变化只要在所附的权利要求限定和确定的本发明精神和范围内，一切利用本发明构思的发明创造均在保护之列。