1.一种利用多优化理论的变压器振声信号重构方法,其特征在于,包括:步骤1,输入实测的变压器振声信号序列S;
步骤2,根据多优化理论对所述变压器振声信号序列S进行重构,重构后的信号序列为SNEW;具体为,SNEW=LBOPT+[STS-I]-1(S-uOPT);其中,BOPT为最佳预测矢量;uOPT为最佳修正矢量;L为系统矩阵。
2.根据权利要求1所述的方法,其特征在于,所述步骤2之前,所述方法还包括:步骤3,求取所述最佳预测矢量BOPT、最佳修正矢量uOPT和系统矩阵L。
3.根据权利要求2所述的方法,其特征在于,所述步骤3包括:步骤301,求取循环延迟矩阵DC,具体为:其中:
sn:所述信号序列S的第n个元素[n=1,2,…,N]N:所述信号序列S的长度步骤302,求取所述循环延迟矩阵DC的奇异值分解,具体为:DC=UΣV
其中:
U:所述循环延迟矩阵DC的左特征矢量矩阵Σ:所述循环延迟矩阵DC的特征值矩阵V:所述循环延迟矩阵DC的右特征矢量矩阵步骤303,求取所述系统矩阵L,具体为:L=[STS][D+I]-1其中:
[D+I]-1:所述矩阵[D+I]的逆矩阵步骤304,求取所述最佳预测矢量BOPT和最佳修正矢量uOPT,具体为:第一步:迭代初始化,具体为:B1=S:所述最佳预测矢量的初始化值u1=UΨVS:所述最佳修正矢量的初始化值k=1:迭代控制参数
其中
系统特征值矩阵[i=1,2,…,NΨ]所述系统特征值矩阵中的第i个特征值所述循环延迟矩阵DC的第iL个特征值NΨ:所述系统特征值矩阵Ψ中非零特征值的个数σS:所述信号序列S的均方差第二步:迭代更新,具体为:其中:
d:用于求取最小值的第一中间参量c:用于求取最小值的第二中间参量v:用于求取最小值的第三中间参量修正因子
第三步:迭代终止,具体为迭代控制参数k加1,重复执行第二步,直至相邻两次迭代结果的差值小于0.001为止,此时k=K,BOPT=BK+1,uOPT=uK+1。
4.一种利用多优化理论的变压器振声信号重构系统,其特征在于,包括:获取模块,输入实测的变压器振声信号序列S;
滤波模块,根据多优化理论对所述变压器振声信号序列S进行重构,重构后的信号序列为SNEW;具体为,SNEW=LBOPT+[STS-I]-1(S-uOPT);其中,BOPT为最佳预测矢量;uOPT为最佳修正矢量;L为修正矩阵。
5.根据权利要求4所述的系统,其特征在于,还包括:计算模块,求取所述最佳预测矢量BOPT、最佳修正矢量uOPT和系统矩阵L。