一种鲁棒性自适应滤波系统

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


背景技术：

2.传统的滤波器是由电容、电感和电阻组成的滤波电路。滤波器可以对电源线中特定频率的频点或该频点以外的频率进行有效滤除，得到一个特定频率的电源信号，或消除一个特定频率后的电源信号。
3.但其不能实现自适应滤波，本发明提供一种鲁棒性自适应滤波系统能够实现自适应滤波。


技术实现要素：

4.本发明所要解决的技术问题是现有技术中存在的不能实现自适应滤波的技术问题。提供一种新的鲁棒性自适应滤波系统，该鲁棒性自适应滤波系统具有自适应滤波、鲁棒性强的特点。
5.为解决上述技术问题，采用的技术方案如下：
6.一种鲁棒性自适应滤波系统，所述鲁棒性自适应滤波系统为时分复用滤波，包括数字滤波器和广义滤波权值预测控制单元；
7.数字滤波器为fir滤波器，包含路数为m的多路复用器，多路复用器的输入连接ram存储器、采样数据输入口，输出连接单路fir滤波器单元，单路fir滤波器单元的输入为第m路的输入x
m
(n)，输出为第m
‑
1路采样数据x
m
‑1(n)和第m路滤波结果y
m
(n)；
8.m为小于等于m的正整数，n为fir横向滤波器阶数，i为小于n的整数；
9.所述广义滤波权值预测控制单元执行下列步骤，完成对数字滤波器的输出估计：
10.步骤1，定义第m路的滤波器单元滤波权值为w
m
(n)，定义第m单路fir滤波器单元的输出y
m
(n)与目标信号d
m
(n)之差e
m
(n)的均方值为均方差误差
[0011][0012]
其中，r＝e[x
m
(n)xt
m
(n)]，p＝e[d
m
(n)x
m
(n)]；
[0013]
步骤2，采用最小均方差方法计算出均方差误差最小时的w
m
(n)作为最优滤波权值；
[0014]
步骤3，定义最优滤波权值调整估计时长j，在时长j内预估并控制调整滤波权值至最优滤波权值。
[0015]
上述方案中，未优化，进一步地，所述步骤3还包括：
[0016]
步骤3.1，定义干扰估计时长p，定义干扰估计系数λ1，最优滤波权值调整系数为λ2，第m单路fir滤波器单元的滤波权值输出估计模型为：
[0017][0018][0019]
1＝e
1，j
(q
‑1)(1+m1q
‑1)δ+q
‑
j
f
1，j
(q
‑1)；
[0020]
1＝e
2，j
(q
‑1)(1+m2q
‑1)δ+q
‑
j
f
2，p
(q
‑1)；
[0021]
其中，w(k)为第k时刻滤波权值，u(k)第k时刻调整滤波权值的输入值，z(k)干扰值，m1(q
‑1)为预定义的q
‑1多项式，n1(q
‑1)为预定义的q
‑1多项式，m2(q
‑1)为预定义的q
‑1多项式，n2(q
‑1)为预定义的q
‑1多项式；c1＝1，c2＝1，k为时刻值，q
‑1为预定义的后移系数，δ＝1
‑
q
‑1为差分系数，ε1(k)和ε2(k)为随机干扰项因子；
[0022]
步骤3.2，采集数据w(k)，z(k)，w(k
‑
1)，z(k
‑
1)，δu(k
‑
1)...δu(k
‑
20)；δu(k一1)...δu(k
‑
c2)为过去调整滤波权值的输入值的差分量；
[0023]
步骤3.3，定义调整滤波权值的输入值的差分量
[0024][0025]
其中，b1为最优滤波权值调整系数，b2干扰模型系数；
[0026]
定义定义计算出
[0027]
步骤3.4，根据u(k)＝u(k
‑
1)+δu(k)计算出控制量；
[0028]
步骤3.5，定义k＝k+1，返回步骤3.2。
[0029]
进一步地，定义干扰估计时长p为两倍最优滤波权值调整估计时长j。
[0030]
进一步地，最优滤波权值调整估计时长j为一个滤波周期。
[0031]
进一步地，所述fir滤波器为fir横向滤波器。
[0032]
本发明的有益效果：为了实现高鲁棒性自适应滤波，本发明使用fpga或者matlab实现fir横向滤波器，通过复用单路滤波器单元实现系统的简化和高鲁棒性。在此基础上，对于权值的调整存在滞后性，进而采用了预估计算法在设定的调整时间内进行权值调整过程和结果的预估计。另外，考虑干扰对于权值系数调整的影响，也进行了估计。
附图说明
[0033]
下面结合附图和实施例对本发明进一步说明。
[0034]
图1，鲁棒性自适应滤波系统示意图。
具体实施方式
[0035]
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
[0036]
实施例1
[0037]
本实施例提供一种鲁棒性自适应滤波系统，如图1，所述鲁棒性自适应滤波系统为时分复用滤波，包括数字滤波器和广义滤波权值预测控制单元；
[0038]
数字滤波器为fir滤波器，包含路数为m的多路复用器，多路复用器的输入连接ram存储器、采样数据输入口，输出连接单路fir滤波器单元，单路fir滤波器单元的输入为第m路的输入xm(n)，输出为第m
‑
1路采样数据x
m
‑1(n)和第m路滤波结果y
m
(n)；
[0039]
m为小于等于m的正整数，n为fir横向滤波器阶数，i为小于n的整数；
[0040]
所述广义滤波权值预测控制单元执行下列步骤，完成对数字滤波器的输出估计：
[0041]
步骤1，定义第m路的滤波器单元滤波权值为w
m
(n)，定义第m单路fir滤波器单元的输出y
m
(n)与目标信号d
m
(n)之差e
m
(n)的均方值为均方差误差
[0042][0043][0044]
其中，r＝e[x
m
(n)xt
m
(n)]，p＝e[d
m
(n)x
m
(n)]；
[0045]
步骤2，采用最小均方差方法计算出均方差误差最小时的w
m
(n)作为最优滤波权值；
[0046]
步骤3，定义最优滤波权值调整估计时长j，在时长j内预估并控制调整滤波权值至最优滤波权值。
[0047]
优选地，所述步骤3还包括：
[0048]
步骤3.1，定义干扰估计时长p，定义干扰估计系数λ1，最优滤波权值调整系数为λ2，第m单路fir滤波器单元的滤波权值输出估计模型为：
[0049][0050][0051]
1＝e
1，j
(q
‑1)(1+m1q
‑1)δ+q
‑
j
f
1，j
(q
‑1)；
[0052]
1＝e
2，j
(q
‑1)(1+m2q
‑1)δ+q
‑
j
f
2，p
(q
‑1)；
[0053]
其中，w(k)为第k时刻滤波权值，u(k)第k时刻调整滤波权值的输入值，z(k)干扰值，m1(q
‑1)为预定义的q
‑1多项式，n1(q
‑1)为预定义的q
‑1多项式，m2(q
‑1)为预定义的q
‑1多项式，n2(q
‑1)为预定义的q
‑1多项式；c1＝1，c2＝1，k为时刻值，q
‑1为预定义的后移系数，δ＝1
‑
q
‑1为差分系数，ε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]056][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]
尽管上面对本发明说明性的具体实施方式进行了描述，以便于本技术领域的技术人员能够理解本发明，但是本发明不仅限于具体实施方式的范围，对本技术领域的普通技术人员而言，只要各种变化只要在所附的权利要求限定和确定的本发明精神和范围内，一切利用本发明构思的发明创造均在保护之列。