申请/专利权人：开沃新能源汽车集团有限公司

申请日：2017-10-20

公开（公告）日：2024-03-22

公开（公告）号：CN107729845B

主分类号：G06F18/10

分类号：G06F18/10;G06F17/14

优先权：

专利状态码：有效-授权

法律状态：2024.03.22#授权;2020.05.22#实质审查的生效;2018.02.23#公开

摘要：本发明公开了一种基于子空间特征值分解的实测频响函数降噪方法，主要包括1初步推导最优估计器，并进行特征值分解；2求解最优估计器，并得到降噪结果。本发明方法不同于传统的降噪方法，是在完成频响函数测试以后进一步对频响函数进行降噪。为了验证方法的有效性，设计六阶系统，并人为引入了噪声，分别利用两种算法对带噪频响函数进行了降噪分析，证明了子空间降噪算法在频响函数降噪中的有效性，最后将降噪算法应用在实测频响函数的降噪中，取得了较好的降噪效果。

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

全文数据：一种基于子空间特征值分解的实测频响函数降噪方法技术领域：[0001]本发明涉及一种基于子空间特征值分解的实测频响函数降噪方法，其属于频响函数降噪处理领域。背景技术：[0002]模态识别、有限元模型修正、故障诊断、结构损伤定位、振动噪声传递路径分析，以及动态载荷识别等分析过程都直接需要对频响函数或者频响函数矩阵进行微分或求逆运算。在子结构频响函数综合中，也需要对界面测点的频响函数矩阵进行求逆运算。求逆或者微分过程对信号噪声十分敏感。因此如何降低实测频响函数中噪声水平成为上述分析能否成功的一个关键。针对实际测量中频响函数的激励与响应存在的噪声，工程上一般采用ft、H2、H3、H4以及Hv等经典的方法，或者Hc估计等其它一些方法对实测频响函数进行估计或者说降噪处理。一般来说，工程应用中使用Hv估计方法较多。以上这些经典的方法属于在测量频响函数时为了排除噪声的影响而采用的估计方法，其中大都采用平均技术。[0003]Sanliturk等利用奇异值分解SVD方法对测试中获得频响函数进行了进一步降噪处理。其基本思路是对带噪频响函数进行相空间重构，利用奇异值分解方法将重构的相空间矩阵分解为信号子空间和噪声子空间，认为前若干个较大的奇异值对应的就是信号子空间，而其它奇异值对应的是噪声子空间，并将其置为零。再利用奇异值分解的逆变换对纯净信号矩阵进行估计，最后由相空间重构的逆过程，通过平均的方法获得降噪后的频响函数。孙鑫晖等人利用同样的方法对频响函数进行了降噪。利用奇异值分解降噪是一种有效的降噪分析方法，在对信号增强与信号特征提取中有较多应用。基于奇异值分解方法实际上属于子空间降噪算法中的一种。子空间方法不同于其它时频降噪方法，如维纳滤波或者卡尔曼滤波等等。子空间降噪算法在语音增强、障诊断以及图像降噪方面均有应用。Ephraim等通过对协方差矩阵的特征值分解EVD，带噪语音信号进行子空间分解，并与谱减法进行了对比，证明其在降噪方面的优势。Hu等将纯净的语音信号与噪声信号的协方差同时对角化，将Ephraim方法推广到了有色噪声的降噪处理中。Hermus等人对目前基于子空间语音增强的方法进行了总结，并将各种子结构语音增强算法进行了对比，指出与谱减法相比，子空间方法减少了音乐噪声的残留，提高了语音识别的准确性。钱征文等利用子空间中的SVD分解对信号进行了降噪处理，利用FFT变换结果确定重构相空间矩阵的秩。吕勇等利用子空间降噪算法对齿轮故障信号进行了分析与降噪，提取出相应的故障特征信号。[0004]综上，H1-H4以及Hv等经典的方法，均是在完成激励与响应点响应测量之后，在获得频响函数之前通过互谱平均的方法对频响函数进行估计，以获得信噪比较高的频响函数。但通常利用这些经典的方法估计后获得的频响函数依然会有噪声残留，这些残留噪声依然会影响频响函数逆运算或微分运算的精度。发明内容：[0005]本发明是为了解决上述现有技术存在的问题而提供一种基于子空间特征值分解的实测频响函数降噪方法，能够有效地除去频响函数的背景噪声。[0006]本发明采用如下技术方案：一种基于子空间特征值分解的实测频响函数降噪方法，包括如下步骤：[0007]步骤1、频响函数的离散傅里叶反变换，子空间算法在时域内进行，根据傅里叶反变换将频响函数变为时域内的脉冲响应函数对于带噪的频响函数，h为脉冲响应函数，H与h均为离散观测序列，η为脉冲响应序列的长度，据离散傅里叶反变换[0009]将离散傅里叶反变换后的带噪脉冲响应函数写成由纯净的脉冲响应函数h’k加上噪声nk构成[0010]hk=h’k+nk，Xk^in-l2[0011]假定频域噪声信号Ns的实部与虚部独立，服从均值为零的高斯分布，且方差相等即：[0012]RealNs〜G0，σ2，ImagNs〜G0，σ2⑶[0013]其中Ns是nk的复数域表达，根据离散傅里叶变换对DFT:[0015]将nk展开成实部与虚部[0017]利用欧拉公式展开，经推导可证明nk也服从高斯分布，期望与方差分别为：[0020]即对于频域带噪信号，若其噪声信号的实部与虚部独立，并服从均值为零的高斯分布，则经过傅里叶反变换后的脉冲响应函数同样服从高斯分布，期望与方差如上式；[0021]步骤2、构造线性滤波器L，将由纯净的脉冲响应函数h’k加上噪声nk构成的带噪脉冲响应函数成向量的形式：[0022]h=!+]!⑻[0023]上式中向量的长度均为ηX1，假定l·的估计可通过h的线性滤波器L构成：[0024]h=Lh9[0025]将残余误差写成由频响函数失真误差W和残余噪声误差£„两个部分构成[0026]ε-h-h'^L-Eh'+Ln=Sb,+eh10[0027]其中，频响函数失真误差定义为：[0028]Sw=L-EIt']][0029]残余噪声误差定义为：[0030]εΑ=Ln12[0031]由上述式⑹、（7、=O^ise五以及aTa=traceaaT推导出频响函数失真能量为：[0033]残余噪声能量为：[0035]其中，R*为协方差矩阵，trace为矩阵的迹，由上述推导可得到以下时域约束最优问题用于求解最优线性估计器：[0038]其中，δ为正常数，即满足残余噪声能量小于δ时，频响函数失真最小，引入拉格朗日乘子μ，μ彡0,构造拉格朗日乘式：16[0040]根据KKT条件：最终可得到最优估计器：[0042]步骤3、特征值分解，假定嵌入维数为m，延迟时间为τ，利用延迟坐标法重构脉冲响应函数的相空间矩阵，重构的相空间矩阵Y写成：[0044]对降噪分析而言，直接取延迟时间τ=1，重构的相空间矩阵Y写成：[0046]其中，Y为Toeplitz矩阵，假设频响函数中存在加性噪声，则相空间矩阵Y可写成：[0047]Y=X+N20[0048]其中，Υ、Χ和N分别表示原始带噪、纯净信号和噪声信号的脉冲响应函数的重构矩阵，带噪脉冲响应函数的协方差可写成：[0049]Rh=Rh+Rn21[0050]采用下式对Rh进行估计：[0052]步骤4、用Rh对Rh'进行估计，Rh为mXm矩阵，将纯净脉冲响应函数的协方差矩阵IV进行EVD分解：[0054]将式23带入式21:[0056]假定Rh'的秩为r，其中rm，则Ah'中的前r个对角元素不为零，即[0058]将式25带入24，则原始脉冲响应的协方差矩阵Rh可写成分块形式：[0060]直接对通过估计获得的Rh进行EVD分解，并按（26进行分块：[0062]对比式26、（27，可知：[0064]将式28带入式25，可知：[0066]步骤5、得出最优估计器，根据式17、（23以及29，最终得最优估计器Ltdc[0068]当μ=1，式30即为所谓线性最小二乘误差LMMSE估计器：[0070]当μ=0,即为所谓最小二乘LS估计器：[0072]步骤6、求解纯净脉冲响应函数，根据式30可知纯净脉冲响应函数估计表示成[0074]步骤7、求解降噪后频响函数，在获得纯净脉冲响应函数估计I后，再利用傅里叶变换即可得到降噪后频响函数；[0075]步骤8、求解降噪结果，依据特征值分析结果，对比保留不同特征值个数对应的降噪效果，确定协方差特征值个数，估算其后验信噪比SNRdB，并取拉格朗日乘子μ，依据最优估计器L™求解得到降噪结果。[0076]本发明具有如下有益效果:考虑实测频响函数噪声问题，本发明提出了基于子空间特征值分解的实测频响函数降噪方法，该方法不同于传统的降噪方法，是在完成频响函数测试以后进一步对频响函数进行降噪。为了验证方法的有效性，设计六阶系统，并人为引入了噪声，分别利用两种算法对带噪频响函数进行了降噪分析，证明了子空间降噪算法在频响函数降噪中的有效性，最后将降噪算法应用在实测频响函数的降噪中，取得了较好的降噪效果。附图说明：[0077]图1为纯净与带噪脉冲响应函数。[0078]图2为纯净与带噪频响函数幅值与相位对比。[0079]图3为协方差矩阵特征值。[0080]图4为不同μ的降噪效果幅值与相位对比。[0081]图5为SVD与EVD降噪效果对比。[0082]图6为实测频响函数Hs2x,8ζ幅值与相位。[0083]图7为脉冲响应函数。[0084]图8为特征值分析结果。[0085]图9为降噪前后频响函数HS2X,8Z幅值与相位对比。具体实施方式：[0086]下面结合附图对本发明作进一步的说明。[0087]本发明基于子空间特征值分解的实测频响函数降噪方法，主要包括：[0088]1初步推导最优估计器，并进行特征值分解；[0089]2求解最优估计器，并得到降噪结果。[0090]进一步地，具体包括如下步骤：[0091]步骤1:频响函数的离散傅里叶反变换。子空间算法一般在时域内进行，根据傅里叶反变换将频响函数变为时域内的脉冲响应函数对于带噪的频响函数，h为脉冲响应函数，这里H与h均为离散观测序列，η为脉冲响应序列的长度，据离散傅里叶反变换[0093]将离散傅里叶反变换后的带噪脉冲响应函数写成由纯净的脉冲响应函数h’k加上噪声nk构成[0094]hk=h’k+nk，Xk^in-l2[0095]假定频域噪声信号Ns的实部与虚部独立，服从均值为零的高斯分布，且方差相等即：[0096]RealNs〜G0，σ2，ImagNs〜G0，σ2⑶[0097]其中Ns是nk的复数域表达，根据离散傅里叶变换对DFT:[0099]将nk展开成实部与虚部[0101]利用欧拉公式展开，经推导可证明nk也服从高斯分布，期望与方差分别为：[0104]即对于频域带噪信号，若其噪声信号的实部与虚部独立，并服从均值为零的高斯分布，则经过傅里叶反变换后的脉冲响应函数同样服从高斯分布，期望与方差如上式。[0105]步骤2:构造线性滤波器L。将由纯净的脉冲响应函数h’k加上噪声nk构成的带噪脉冲响应函数成向量的形式：[0106]h=!+]!⑻[0107]上式中向量的长度均为ηX1。假定h的估计可通过h的线性滤波器L构成：[0108]it二Ui^9[0109]将残余误差可写成由频响函数失真误差W和残余噪声误差^两个部分构成[0111]其目的是保证残余误差低于一定阈值的同时使频响函数失真度最小。显然该方式是频响函数失真和残余误差程度之间的一种折中。其中，[0112]频响函数失真误差定义为：[0114]残余噪声误差定义为：[0115]Sn^Ln12[0116]由上述式6、（7、βη'=O^ise疋以及aTa=traceaaT可推导出频响函数失真能量为：[0118]残余噪声能量为：[0120]其中，R*为协方差矩阵，trace为矩阵的迹。[0121]由上述推导可得到以下时域约束最优问题用于求解最优线性估计器：[0124]其中，δ为正常数，即满足残余噪声能量小于δ时，频响函数失真最小。引入拉格朗日乘子μ，μ彡0,构造拉格朗日乘式：16[0126]根据KKT条件：最终可得到最优估计器：[0128]可见，如果已知Rh'，即可得到最优估计器LTDC，但通常来说Rh'是未知的，需要利用Rh进行估计。[0129]步骤3:特征值分解。假定嵌入维数为m，延迟时间为τ，利用延迟坐标法重构脉冲响应函数的相空间矩阵。重构的相空间矩阵Y可以写成：[0131]对降噪分析而言，一般直接取延迟时间τ=1，重构的相空间矩阵Y可以写成：[0133]其中，Y为Toeplitz矩阵，当然也可以将脉冲响应序列重构成Hankel矩阵形式。假设频响函数中存在加性噪声，则相空间矩阵Y可写成：[0134]Y=X+N20[0135]其中，Υ、Χ和N分别表示原始带噪、纯净信号和噪声信号的脉冲响应函数的重构矩阵。带噪脉冲响应函数的协方差可写成：[0136]Rh=Rh+Rn21[0137]可采用下式对Rh进行估计：[0139]步骤4:用Rh对Rh'进行估计。Rh为mXm矩阵，将纯净脉冲响应函数的协方差矩阵IV进行EVD分解：[0140]i?h,=UiilAvUv-23[0141]将式23带入式21:[0143]假定Rh的秩为r，其中rm，则Ah中的前r个对角元素不为零，即[0145]将式25带入24，则原始脉冲响应的协方差矩阵Rh可写成分块形式：[0147]直接对通过估计获得的Rh进行EVD分解，并按26进行分块：[0149]对比式26、（27，可知：[0150]RhAhl=Ab,,+σ；ΰίχΕτ28；[0151]将式28带入式25，可知：[0153]步骤5:得出最优估计器。根据式17、（23以及29，最终可得最优估计器Ltdc[0155]事实上，当μ=1，式30即为所谓线性最小二乘误差LMMSE估计器：[0157]当μ=0,即为所谓最小二乘LS估计器：[0158]Lw—Lxqc^0=UmUu32[0159]步骤6:求解纯净脉冲响应函数。根据式30可知纯净脉冲响应函数估计可表示成[0160]h=I^mh33[0161]步骤7:求解降噪后频响函数。在获得纯净脉冲响应函数估计i后，再利用傅里叶变换即可得到降噪后频响函数。[0162]步骤8:求解降噪结果。依据特征值分析结果，对比保留不同特征值个数对应的降噪效果，确定协方差特征值个数，估算其后验信噪比SNRdB，并取拉格朗日乘子μ，依据最优估计器L™求解得到降噪结果。[0163]算例一：[0164]为验证基于子空间特征值分解的算法在频响函数降噪中的效果，设计一个六阶系统，其传递函数为-其中s为Laplace算子。系统参数设定如下:mi=l，m2=0.5，m3=0.8，ki=403i2，k2=4003i2，k3=40003T2，.，给定采样率为100Hz，数据样本长度为1000,令传递函数中S=jω，计算纯净频响函数的单位脉冲响应，并叠加“加性”高斯白噪声，Matlab的命令为awgnX，10,’measured’），信噪比SNR为10dB。[0165]加入噪声前、后的脉冲响应函数如图1。[0166]经过傅里叶变换后的频响函数幅值与相位如图2。[0167]步骤1:特征值分解。利用式（22估计脉冲响应函数重构相空间矩阵的协方差矩阵，并进行特征值分解。图3给出了前50个特征值的分析结果。[0168]步骤2:最优估计器与降噪结果。根据图3,取协方差特征值个数为6。注意式30，最优估计器Ltdc与拉格朗日乘子μ有关。图4给出μ=0〜10时降噪结果。[0169]由图4可知，显然降噪后频响函数无论是幅值与相位，其精度大幅度提高，特别是对频响函数的反共振点的估计。但降噪效果与μ的大小并非是单调的，观察发现当με〇,5时，降噪效果最佳。目前关于μ取值尚无统一的计算公式，这里给出语音增强算法中μ取值的一种估计方法供参考。[0171]其中，μ=4.2,c=6.25,SNRdB为后验信噪比估计。这里SNRdB约为10dB，计算可得μ=2.6,与前述分析的με〇,5的最优取值范围是一致的。[0172]步骤3:降噪效果对比。为了检验在低信噪比情况下基于EVD与SVD两种算法的降噪效果，将信噪比下调至5dB。[0173]图5为采用SVD_LS方法与EVDy=3两种方法降噪后的结果。其中，SVD_LS为基于子空间SVD最小二乘意义下的降噪算法。[0174]显然即使在信噪比达到5dB时两种方法依然具有较好的降噪效果，而从降噪后的相位对比来看，EVD方法降噪效果比SVD_LS方法更优。[0175]算例二：[0176]以某微车的实测跨点频响函数为例，利用EVD方法的时域约束估计方法对实测频响函数进行降噪。其中频响函数通过锤击实验法得到，采样率为1024Hz，频率分辨率为0.125Hz，频响函数的估计方法采用H1估计。[0177]图6给出了车身左前悬架+z激励下驾驶室内右前导轨+X测点响应的实测频响函数。[0178]步骤1:特征值分解。图7给出了经过傅里叶反变换得到的实测频响函数的脉冲响应函数。构造Topelitz的频响函数相空间矩阵Y，并对Y进行特征值分解.[0179]步骤2:最优估计器与降噪结果及对比。图8给出了脉冲响应函数重构相空间矩阵的协方差矩阵的前100个特征值的分析结果。根据图8特征值分析结果，并对比保留不同特征值个数对应的降噪效果，最终确定协方差特征值个数为30,经过估算其后验信噪比SNRdB约为16dB，最终取拉格朗日乘子μ=1.6，并由式30，确定最优估计器Ltdc。[0180]图9给出了TDC降噪结果，由图可知采用EVD方法的时域约束估计方法降噪结果。显然，结果降噪后能明显提高频响函数的精度，特别是对反共振点的改善。[0181]以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下还可以作出若干改进，这些改进也应视为本发明的保护范围。

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

百度查询： 开沃新能源汽车集团有限公司 一种基于子空间特征值分解的实测频响函数降噪方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD