[發明專利]一種天文時序信號的噪聲估計方法

有效

• 申請號: 202010540771.5
• 申請日: 2020-06-15
• 公開（公告）號: CN111863016B
• 公開（公告）日: 2022-09-02
• 發明（設計）人: 于蘭;馮松;孟瑤
• 申請（專利權）人: 云南國土資源職業學院
• 主分類號: G10L21/0216
• 分類號: G10L21/0216;G10L25/21;G10L25/27;G10L25/45
• 代理公司: 昆明人從眾知識產權代理有限公司 53204
• 代理人: 陳波

1.一種天文時序信號的噪聲估計方法，其特征在于：首先獲取天文時序信號，對該時序信號進行歸一化處理；接著，將歸一化后的信號和窗函數相乘并運用周期圖法得到信號功率譜D；然后根據功率譜D定義該信號功率譜的噪聲模型M和模型中的參數θ；定義該功率譜D在模型M和參數θ下的似然函數，以及參數θ的先驗分布函數，根據貝葉斯公式定義參數θ在功率譜D和模型M下的后驗概率分布函數；運用漢密爾頓蒙特卡洛算法來對后驗概率分布函數進行采樣，來獲取模型M中參數θ的值，獲得對該功率譜D的噪聲模型M的最佳擬合曲線方程；對最佳擬合曲線方程進行卡方檢驗，計算置信水平。

2.根據權利要求1所述的天文時序信號的噪聲估計方法，其特征在于：所述方法具體步驟如下：

步驟1：獲取天文時序信號，并利用該信號的均值和標準差對該時序信號進行歸一化處理；

步驟2：將歸一化后的信號和海明窗函數相乘并運用周期圖法得到該信號功率譜D；

步驟3：根據功率譜D設計該功率譜的噪聲模型M和模型參數θ；定義該功率譜D在模型M和參數θ下的似然函數p(D|θ，M)，以及參數θ的先驗分布函數p(θ|M)，根據貝葉斯公式定義參數θ在功率譜D和模型M下的后驗概率分布函數p(θ|D，M)＝p(D|θ，M)p(θ|M)；

步驟4：根據后驗概率分布函數p(θ|D，M)，利用漢密爾頓蒙特卡洛算法對后驗概率分布進行采樣，來獲取模型M中參數θ的值，獲得對該功率譜D的噪聲模型M的最佳擬合曲線方程；

步驟5：對最佳擬合曲線方程進行卡方檢驗，計算99％的置信水平：將位于99％置信水平之上的頻率成分認為是真正的振蕩信號，在99％置信水平之下的認為是噪聲。