一种基于子空间特征值分解的实测频响函数降噪方法

专利主权项内容

1.一种基于子空间特征值分解的实测频响函数降噪方法，其特征在于：包括如下步骤：步骤1、频响函数的离散傅里叶反变换，子空间算法在时域内进行，根据傅里叶反变换将频响函数变为时域内的脉冲响应函数对于带噪的频响函数，h为脉冲响应函数，H与h均为离散观测序列，n为脉冲响应序列的长度, 据离散傅里叶反变换将离散傅里叶反变换后的带噪脉冲响应函数写成由纯净的脉冲响应函数h'加上噪声n构成kkh＝h'+n, 0≤k≤n-1 (2)kkk假定频域噪声信号N的实部与虚部独立，服从均值为零的高斯分布，且方差相等即：sReal(N)～G(0, σ), Imag(N)～G(0, σ) (3)s2s2其中N是n的复数域表达，根据离散傅里叶变换对(DFT)：sk将n展开成实部与虚部k利用欧拉公式展开，经推导可证明n也服从高斯分布，期望与方差分别为 : kE[n]＝0 (6)k即对于频域带噪信号，若其噪声信号的实部与虚部独立，并服从均值为零的高斯分布，则经过傅里叶反变换后的脉冲响应函数同样服从高斯分布，期望与方差如上式；步骤2、构造线性滤波器L，将由纯净的脉冲响应函数h'加上噪声n构成的带噪脉冲响应函数成向量的形式：kkh＝h′+n (8)上式中向量的长度均为n×1，假定h′的估计可通过h的线性滤波器L构成：将残余误差写成由频响函数失真误差ε和残余噪声误差ε两个部分构成h′n其中，频响函数失真误差定义为：残余噪声误差定义为：由上述式(6)、(7)、以及aa＝trace(aa)推导出频响函数失真能量为：TT残余噪声能量为：其中，R为协方差矩阵，trace为矩阵的迹，由上述推导可得到以下时域约束最优问题用于求解最优线性估计器：*其中, δ为正常数，即满足残余噪声能量小于δ时，频响函数失真最小，引入拉格朗日乘子μ，μ≥0，构造拉格朗日乘式：根据KKT条件 : 以及/>最终可得到最优估计器：步骤3、特征值分解，假定嵌入维数为m，延迟时间为τ，利用延迟坐标法重构脉冲响应函数的相空间矩阵，重构的相空间矩阵Y写成：对降噪分析而言，直接取延迟时间τ＝1，重构的相空间矩阵Y写成：其中，Y为Toeplitz矩阵，假设频响函数中存在加性噪声，则相空间矩阵Y可写成：Y＝X+N (20)其中，Y、X和N分别表示原始带噪、纯净信号和噪声信号的脉冲响应函数的重构矩阵，带噪脉冲响应函数的协方差可写成：R＝R+R (21)hh′n采用下式对R进行估计：h步骤4、用R对R进行估计，R为m×m矩阵, 将纯净脉冲响应函数的协方差矩阵R进行EVD分解：hh′hh′将式(23)带入式(21)：假定R的秩为r，其中r＜m，则Λ中的前r个对角元素不为零，即h′h′将式(25)带入(24)，则原始脉冲响应的协方差矩阵R可写成分块形式：h直接对通过估计获得的R进行EVD分解，并按(26)进行分块：h对比式(26)、(27)，可知：将式(28)带入式(25)，可知：步骤5、得出最优估计器，根据式(17)、(23)以及(29)，最终得最优估计器LTDC当μ＝1，式(30)即为所谓线性最小二乘误差(LMMSE)估计器：当μ＝0，即为所谓最小二乘(LS)估计器：步骤6、求解纯净脉冲响应函数，根据式(30)可知纯净脉冲响应函数估计表示成步骤7、求解降噪后频响函数，在获得纯净脉冲响应函数估计后，再利用傅里叶变换即可得到降噪后频响函数；步骤8、求解降噪结果，依据特征值分析结果，对比保留不同特征值个数对应的降噪效果，确定协方差特征值个数，估算其后验信噪比SNR，并取拉格朗日乘子μ，依据最优估计器L求解得到降噪结果。dBTDC