[发明专利]一种利用特征选择的变压器状态振声检测信号滤波方法和系统

权利要求书

1.一种利用特征选择的变压器状态振声检测信号滤波方法，其特征在于，包括：

步骤101获取按时间顺序采集的信号序列S；

步骤102求取PCA延迟常数，具体为：所述PCA延迟常数κ为平均矩阵B的非零特征值个数；其中，所述平均矩阵B的求取公式为B＝[S-m₀]^T[S-m₀]；m₀为所述信号序列S的均值；

步骤103求取PCA延迟序列，具体为：所述PCA延迟序列S^pca中的第n个元素计算公式为其中，n为元素序号，其取值范围为n＝1,2,…,N；N为所述信号序列S的长度；为所述信号序列S的第|n+κ|_N个元素；||_N表示以N为模下取整运算；

步骤104迭代控制参数初始化，所述迭代控制参数k的初始化值为k＝0；

步骤105PCA滤波矩阵初始化，具体为：所述PCA滤波矩阵的初始化值为P⁰，其第i行第j列元素的计算公式为其中，i为行序号，其取值范围为i＝1,2,…,N；j为列序号，其取值范围为j＝1,2,…,N；s_i为所述信号序列S的第i个元素；s_j为所述信号序列S的第j个元素；λ_min为所述平均矩阵B的非零最小特征值；λ_max为所述平均矩阵B的最大特征值；

步骤106迭代更新，具体为：所述PCA滤波矩阵的迭代更新值为P^k+1，其计算公式为其中，U为所述平均矩阵B的左特征矢量矩阵；Γ为所述平均矩阵B的特征值矩阵；

步骤107迭代结束，具体为：计算迭代误差ε，其计算公式为ε＝||P^k+1-P^k||；如果所述迭代误差ε大于或者等于迭代阈值则所述迭代控制参数k的值加1；并返回所述步骤106和所述步骤107，直至所述迭代误差ε小于所述迭代阈值并将所述PCA滤波矩阵的当前值P^k+1赋值给PCA最佳滤波矩阵P_OPT；其中，所述迭代阈值的取值为

步骤108求取滤除噪声后的信号序列，具体为：所述滤除噪声后的信号序列S_new的计算公式为S_new＝P_OPTS。

2.一种利用特征选择的变压器状态振声检测信号滤波系统，其特征在于，包括：

模块201获取按时间顺序采集的信号序列S；

模块202求取PCA延迟常数，具体为：所述PCA延迟常数κ为平均矩阵B的非零特征值个数；其中，所述平均矩阵B的求取公式为B＝[S-m₀]^T[S-m₀]；m₀为所述信号序列S的均值；

模块203求取PCA延迟序列，具体为：所述PCA延迟序列S^pca中的第n个元素计算公式为其中，n为元素序号，其取值范围为n＝1,2,…,N；N为所述信号序列S的长度；为所述信号序列S的第|n+κ|_N个元素；||_N表示以N为模下取整运算；

模块204迭代控制参数初始化，所述迭代控制参数k的初始化值为k＝0；

模块205PCA滤波矩阵初始化，具体为：所述PCA滤波矩阵的初始化值为P⁰，其第i行第j列元素的计算公式为其中，i为行序号，其取值范围为i＝1,2,…,N；j为列序号，其取值范围为j＝1,2,…,N；s_i为所述信号序列S的第i个元素；s_j为所述信号序列S的第j个元素；λ_min为所述平均矩阵B的非零最小特征值；λ_max为所述平均矩阵B的最大特征值；

模块206迭代更新，具体为：所述PCA滤波矩阵的迭代更新值为P^k+1，其计算公式为其中，U为所述平均矩阵B的左特征矢量矩阵；Γ为所述平均矩阵B的特征值矩阵；

模块207迭代结束，具体为：计算迭代误差ε，其计算公式为ε＝||P^k+1-P^k||；如果所述迭代误差ε大于或者等于迭代阈值则所述迭代控制参数k的值加1；并返回所述模块206和所述模块207，直至所述迭代误差ε小于所述迭代阈值并将所述PCA滤波矩阵的当前值P^k+1赋值给PCA最佳滤波矩阵P_OPT；其中，所述迭代阈值的取值为

模块208求取滤除噪声后的信号序列，具体为：所述滤除噪声后的信号序列S_new的计算公式为S_new＝P_OPTS。