[发明专利]一种变压器运行状态振声检测信号的滤波方法和装置

权利要求书

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

步骤3，生成投影矩阵D_t；所述步骤3进一步包括：

步骤31，输入实测的训练用样本数据S＝[S₁,S₂,......,S_N]^T,N为振动信号序列的长度；[*]^T表示矢量*的转置；所述训练用样本数据的选择原则是：振声信号的信噪比大于预定阈值的数据；

步骤32，估计样本数据的分布概率；所述步骤32包括：

步骤321，确定样本数据的分布范围Ω：Ω＝[S_MIN S_MAX]

最大值S_MAX：

最小值S_MIN：

步骤322，确定区间个数；

区间个数M：M＝[50～80]N，即区间个数是数据个数的50到80倍；

步骤323，根据所述定区间个数，计算样本数据的分布概率；所述步骤323包括：将分布范围Ω划分为M个区间，每个区间的分布范围为

统计落入区间Ω_m的训练数据的个数C_m，并计算分布概率F_m：

步骤33，初始化迭代参数，该步骤具体为：

生成第一参数矩阵A₀：A₀＝[a_ij]_K×K＝rand[0 1]，表示第一参数矩阵中的每一个元素a_ij从[0 1]范围内均匀随机抽取；

生成第二参数矩阵B₀：B₀＝[b_ij]_K×K＝rand[0 1]，表示第二参数矩阵中的每一个元素b_ij从[0 1]范围内均匀随机抽取；

生成初始的投影矩阵D₀：D₀＝[d_ij]_N×K＝rand[0 1]，表示初始的投影矩阵中的每一个元素d_ij从[0 1]范围内均匀随机抽取；

表示上取整

步骤34，设置t＝1,开始迭代；

步骤35，根据样本数据的分布概率F_m，从分布范围Ω＝[S_MIN S_MAX]中随机抽取N次，所抽取的数据组成随机序列

步骤36，计算修正系数α_t:

其中

α_t为第t步迭代得到的修正系数矢量，为一K维矢量

S_t为第t步迭代得到的随机抽取的样本序列

为*的L₂模的平方

||*||₁为*的L₁模

λ为小于1的任意正数

步骤37，根据所述修正系数，更新迭代参数；

其中：表示α_t的转置

步骤38，根据所述迭代参数，更新投影矩阵；

其中

tr(*):矩阵*的迹

矩阵D_t-1的转置

步骤39，根据所述投影矩阵，判断迭代是否结束；具体为：

如果则迭代结束，进入滤波步骤：

如果则t＝t+1，返回重新迭代；

步骤1，输入一实测的振声信号序列S；

步骤2，对所述振声信号序列S进行滤除噪声处理，生成滤除噪声后的数据序列S_new；具体为：

其中，为矩阵D_t的转置。