[发明专利]一种利用海斯矩阵的变压器振声信号重构方法和系统

权利要求书

1.一种利用海斯矩阵的变压器振声信号重构方法，其特征在于，包括：

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

步骤2，根据海斯矩阵对所述变压器振声信号序列S进行重构，重构后的信号序列为S_NEW；具体为，

其中，D为海斯预测矩阵；W_OPT为海斯最佳权重矩阵；λ_OPT为海斯最佳调整矩阵；α为海斯调整因子；μ为海斯因子。

2.根据权利要求1所述的方法，其特征在于，所述步骤2之前，所述方法还包括：

步骤3，求取所述海斯预测矩阵D、海斯最佳权重矩阵W_OPT、海斯最佳调整矩阵λ_OPT、海斯调整因子α和海斯因子μ。

3.根据权利要求2所述的方法，其特征在于，所述步骤3包括：

步骤301，求取循环延迟矩阵D_C，具体为：

其中：

s_n：所述信号序列S的第n个元素[n＝1,2,…,N]

N：所述信号序列S的长度

步骤302，求取所述海斯预测矩阵D，具体为：

D＝[S^TS-I][ΣV+I]D_C

其中：

Σ：矩阵[D_C+I]^-1的特征值矩阵

V：矩阵[D_C+I]^-1的右特征矢量矩阵

I：单位矩阵

步骤303，求取所述海斯调整因子α，具体为：

其中：

矩阵[S^TS][D+I]^-1的最大特征值

矩阵[S^TS][D+I]^-1的最小特征值

步骤304，求取所述海斯因子μ，具体为：

其中：

L_MAX：矩阵[S^TS][D+I]^-1中所有元素绝对值的最大值

L_MIN：矩阵[S^TS][D+I]^-1中所有元素绝对值的最小值

步骤305，迭代求取所述海斯最佳权重矩阵W_OPT、海斯最佳调整矩阵λ_OPT，具体为：

第一步：迭代初始化，具体为：

λ₁＝[S-m_S]^T[S-m_S]：所述海斯最佳调整矩阵的初始化值

W₁＝ΨV:所述海斯最佳权重矩阵的初始化值

k＝1：迭代控制参数

其中

所述W₁的特征值矩阵[i＝1,2,…,N_Ψ]

所述矩阵W₁的特征值矩阵中的第i个特征值

ζ_iL：所述循环延迟矩阵D_C的第i_L个特征值

N_Ψ：所述矩阵W₁的特征值矩阵Ψ中非零特征值的个数

m_S：所述信号序列S的均值

第二步：迭代更新，具体为：

其中：

d：用于求取最小值的第一中间参量

c：用于求取最小值的第二中间参量

v：用于求取最小值的第三中间参量

U_O：矩阵[S^TS-I]^-1的左特征矢量矩阵

Σ_O：矩阵[S^TS-I]^-1的特征值矩阵

第三步：迭代终止，具体为迭代控制参数k加1，重复执行第二步，直至相邻两次迭代结果的差值小于0.001为止，此时k＝K，W_OPT＝W_K+1和λ_OPT＝λ_K+1。