[發(fā)明專利]頻譜泄漏校正的多頻實(shí)信號(hào)頻率估計(jì)算法

在審

• 申請(qǐng)?zhí)枺?/td>: 202011531792.7
• 申請(qǐng)日: 2020-12-16
• 公開（公告）號(hào): CN112881796A
• 公開（公告）日: 2021-06-01
• 發(fā)明（設(shè)計(jì)）人: 陳鵬;蘇欣;賈琦
• 申請(qǐng)（專利權(quán)）人: 中國空氣動(dòng)力研究與發(fā)展中心設(shè)備設(shè)計(jì)及測(cè)試技術(shù)研究所
• 主分類號(hào): G01R23/16
• 分類號(hào): G01R23/16
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.頻譜泄漏校正的多頻實(shí)信號(hào)參數(shù)估計(jì)算法，其特征在于：適用對(duì)象為含噪多頻實(shí)信號(hào)的參數(shù)估計(jì)。

該方法包括以下步驟：

第一步：利用式Y(jié)(k)＝FFT(x(n))對(duì)采樣信號(hào)x(n)進(jìn)行預(yù)處理，并求取頻譜中M個(gè)極大值點(diǎn)的索引k_m。

式中：n表示采樣信號(hào)x(n)時(shí)刻點(diǎn)，n＝0，1，…，N-1，N表示信號(hào)長(zhǎng)度，F(xiàn)FT(·)表示對(duì)·進(jìn)行快速傅里葉變換；k_m表示頻譜中能量極大值的索引，m＝1，2，…，M，M表示頻率分量個(gè)數(shù)。

第二步：首先設(shè)定頻譜偏移量初值復(fù)幅值初值然后利用對(duì)信號(hào)進(jìn)行頻譜插值，插值間隔0.5；并根據(jù)估計(jì)值構(gòu)造非i分量正頻率插值點(diǎn)頻譜和所有負(fù)頻率的插值點(diǎn)頻譜

其次利用進(jìn)行頻譜泄漏校正，利用計(jì)算頻譜偏移量，并利用求取復(fù)幅值。

最后通過循環(huán)計(jì)算，并利用式得到各頻率分量粗略的參數(shù)估計(jì)值。

式中：下標(biāo)i表示第i頻率分量，1≤i≤M，表示參數(shù)·的估計(jì)值，和分別表示采樣信號(hào)、第i分量正頻率、非i分量正頻率和所有負(fù)頻率的插值點(diǎn)頻譜，復(fù)幅值A(chǔ)與A^*互為共軛，Re[·]、Im[·]、|·|和∠·分別表示取復(fù)數(shù)·的實(shí)部、虛部、模和角度。第三步：重復(fù)第二步，得到各頻率分量精確的參數(shù)估計(jì)值。