[發明專利]一種頻率估計方法、裝置、系統及計算機可讀存儲介質

有效

• 申請號: 201910749415.1
• 申請日: 2019-08-14
• 公開（公告）號: CN110442827B
• 公開（公告）日: 2023-05-26
• 發明（設計）人: 馬志峰;曲芳;李洪濤;沈黎明;孔凡軍;楊忠禮;蘇高峰;劉東;孫志強;王新銘;王斌;牛元泰;趙慧光;王智慧
• 申請（專利權）人: 國網河南省電力公司新鄉供電公司;國家電網有限公司
• 主分類號: G06F18/2131
• 分類號: G06F18/2131;G06F18/10;G01K11/00;G06F17/14;G06F17/16
• 代理公司: 北京集佳知識產權代理有限公司 11227
• 代理人: 侯珊

1.一種頻率估計方法，其特征在于，包括：

對聲表面波傳感器的回波信號進行采樣，獲得信息矩陣；所述回波信號為所述聲表面波傳感器內的諧振器返回的信號，返回的這個信號與被測物體的溫度有關；

對所述信息矩陣重構獲得重構矩陣；

對所述重構矩陣進行降噪處理，獲得降噪矩陣；

對所述降噪矩陣進行補零處理；

將所述降噪矩陣進行快速傅里葉變換FFT，并根據FFT變換的結果計算所述回波信號的頻率；

所述對所述信息矩陣重構獲得重構矩陣的過程具體為：

將所述信息矩陣利用時間序列重構吸引子軌跡矩陣；

根據所述吸引子軌跡矩陣，獲得重構矩陣；

所述對所述重構矩陣進行降噪處理，獲得降噪矩陣的過程具體為：

利用奇異值分解SVD對所述信息矩陣進行降噪處理，獲得降噪矩陣；

獲得降噪矩陣的過程包括：

重構矩陣B形式如下：

取一系列不同的行數m構造矩陣，根據下式：

η_i＝λ_i/(λ₁+λ₂+…+λ_p),i＝1,2,…p

利用重構矩陣的奇異值計算各分量信號的信息量，其中，λ_i為相應奇異值，η_i為該奇異值的信息量；

確定降噪矩陣的行數m＝i，則n＝int(N/m)，其中，N為時間序列點；

選取前兩個奇異值作為信號部分，其余奇異值作為噪聲部分濾掉，則得到降噪矩陣；

將所述降噪矩陣進行快速傅里葉變換FFT，并根據FFT變換的結果計算所述回波信號的頻率的過程包括：

將降噪矩陣還原成時間序列，并利用FFT變換公式對降噪處理后的時間序列信號進行N點的FFT，可用下式表示：

其中，k＝0,1,2,…,N-1；e為指數，X_k為回波信號的頻率；x_n為第n個信息分量，N為時間序列點；

求得令|X_k|取最大值時的k值，記為k_max，則f＝k_max*f_s/N；

其中，f_s為采樣頻率，f即為所求頻率。

2.根據權利要求1所述的頻率估計方法，其特征在于，在所述根據FFT變換的結果計算所述回波信號的頻率之后，該方法還包括：

對所述頻率進行校正。

3.根據權利要求2所述的頻率估計方法，其特征在于，所述對所述頻率進行校正的過程具體為：

利用高斯擬合算法對所述頻率進行校正。