1.一种利用去模糊算子的功率信号滤波方法,其特征在于,包括:步骤1,输入实测的功率信号序列S;
步骤2,求取循环延迟矩阵DC,具体为:其中:
sn:所述信号序列S的第n个元素[n=1,2,…,N];
N:所述信号序列S的长度;
步骤3,求取去模糊矩阵D,具体为:T
D=[I+DC][SS‑UΣ1V];
其中:
新特征值矩阵;
所述新特征值矩阵的第i个特征值[i=1,2,…,N];
T
τi:矩阵SS的第i个特征值[i=1,2,…,N];
T
U:矩阵SS的左特征矢量矩阵;
T
V:矩阵SS的右特征矢量矩阵;
σS:所述信号序列S的均方差;
I:单位矩阵;
步骤4,求取模糊因子γ,具体为:其中:
T
矩阵[SS]的所有小于 的特征值中的最大值;
矩阵 的所有小于 的特征值中的最小值;
步骤5,求取去模糊因子μ,具体为:步骤6,求取去模糊调整因子α,具体为:其中:
mS:所述信号序列S的均值;
步骤7,迭代求取模糊化矢量WOPT和去模糊化矢量λOPT,具体为:第一步:迭代初始化,具体为:λ1=[S‑mS]:所述去模糊化矢量的初始化值;
W1=VS:所述模糊化矢量的初始化值;
k=1:迭代控制参数;
其中
T
V:矩阵SS的右特征矢量矩阵;
第二步:迭代更新,具体为:T
λk+1=λk+γ[SS‑I]Wk;
其中:
T ‑1
UO:矩阵[SS‑I] 的左特征矢量矩阵;
T ‑1
ΣO:矩阵[SS‑I] 的特征值矩阵;
第三步:迭代终止,具体为:迭代控制参数k加1,重复执行第二步,直至相邻两次迭代结果的差值小于0.001为止,此时k=K,WOPT=WK+1和λOPT=λK+1;
步骤8,对所述功率信号序列S进行滤除噪声处理,滤除噪声后的功率信号序列为SNEW;
具体为: