[發(fā)明專利]一種基于直接法求解瞬時頻率的毛刺定位與消除方法

有效

• 申請?zhí)枺?/td>: 201711344919.2
• 申請日: 2017-12-15
• 公開（公告）號: CN108134604B
• 公開（公告）日: 2021-01-15
• 發(fā)明（設(shè)計）人: 張慶;賈林山;何小高
• 申請（專利權(quán)）人: 西安交通大學
• 主分類號: H03L7/085
• 分類號: H03L7/085;H04L27/36
• 代理公司: 西安智大知識產(chǎn)權(quán)代理事務(wù)所 61215
• 代理人: 賀建斌

1.一種基于直接法求解瞬時頻率的毛刺定位與消除方法，其特征在于，包括以下步驟：

1)使用局部均值分解方法將原始信號x(t)分解為乘積函數(shù)分量，從而得到乘積函數(shù)的調(diào)頻信號s(t)，將調(diào)頻信號s(t)所有的局部極大值點MaxE_i置為1，所有的局部極小值點MinE_i置為-1，且滿足-1≤s(t)≤1；

2)將調(diào)頻信號s(t)取反余弦，得到非單調(diào)的瞬時相位

3)對非單調(diào)的瞬時相位進行展開，得到單調(diào)遞增的瞬時相位其展開方法為：將非單調(diào)的瞬時相位的所有局部極值點進行編號；如果第一個局部極值點是局部極大值點，那么，如果局部極值點的編號i為奇數(shù)，第i極值點與第i+1極值點之間的相位重新定義為如果i為偶數(shù)，第i與i+1個極值點之間的相位重新定義為如果第一個局部極值點是局部極小值點，那么，如果局部極值點的編號i為奇數(shù)，第i極值點與第i+1極值點之間的相位重新定義為如果i為偶數(shù)，第i與i+1個極值點之間的相位重新定義為最終得到單調(diào)遞增的瞬時相位

4)對單調(diào)遞增的瞬時相位進行前向差分，設(shè)調(diào)頻信號s(t)的采樣頻率為f_s，得到調(diào)頻信號s(t)的帶有毛刺的瞬時頻率f₀(t)：

上式中diff表示前向差分運算；

5)確定帶有毛刺的瞬時頻率f₀(t)的毛刺位置：確定非單調(diào)的瞬時相位的局部極值點E，針對第i個局部極值點為E_i，將E_i位置及其兩側(cè)各2個采樣點位置共計5個點按幅值的絕對值大小進行降序排序，分別標識為D₁～D₅，其對應(yīng)采樣時刻分別為

設(shè)如果Z＜0，那么，若則將與存儲至毛刺位置序列T中，若則將與存儲至毛刺位置序列T中；如果Z＞0，若則將與存儲至毛刺位置序列T中，若則將與存儲至毛刺位置序列T中；最終得到毛刺位置序列T；

6)將步驟4)得到的帶有毛刺的瞬時頻率f₀(t)在步驟5)得到的毛刺位置序列T處的數(shù)據(jù)刪除，并使用三次樣條插值補出毛刺位置序列T處的瞬時頻率值，最終得到無毛刺的瞬時頻率f(t)。