[發明專利]基于極值似然比的復雜調制信號處理結果可信性校驗方法

有效

• 申請號: 201810618697.7
• 申請日: 2018-06-14
• 公開（公告）號: CN108988987B
• 公開（公告）日: 2021-02-26
• 發明（設計）人: 胡國兵;吳珊珊;董振華;趙嬪嬌;楊莉
• 申請（專利權）人: 金陵科技學院
• 主分類號: H04L1/00
• 分類號: H04L1/00;H04L27/32
• 代理公司: 南京鐘山專利代理有限公司 32252
• 代理人: 戴朝榮

1.一種基于極值似然比的復雜調制信號處理結果可信性校驗方法，其特征在于，包括以下步驟：

步驟一、參考信號構建：對觀測信號進行調制方法識別，根據識別結果對應的模型，估計相應參數并構建參考信號；所述步驟一中，建立疊加了高斯白噪聲的LFM/BPSK混合調制信號模型為：

x(n)＝s(n)+w(n)

＝Aexp[j(2πf₀Δtn+πlΔt²n²+θ(n)+θ₀)]+w(n),0≤n≤N-1

其中，A為信號幅度，j為虛數單位，f₀為起始頻率，Δt為采樣間隔，l為調頻斜率，BPSK分量的相位函數θ(n)＝πd₂(n)，d₂(n)為二元編碼信號，其碼元寬度為T_c，碼元個數N_c，碼字為c_m,m＝1,...,N_c，θ₀為初相位，N為樣本點數，w(n)為零均值加性復高斯白噪聲過程，方差為2σ²；

構建參考信號y(n)，建立假設檢驗問題H₀和H₁，H₀表示調制方式識別結果正確且無解碼錯誤，H₁表示調制方式識別錯誤或存在至少一位錯誤解碼；

步驟二、噪聲方差估計：對觀測信號的噪聲方差進行估計；所述步驟二中，利用二階四階矩方法對接收信號的方差σ²進行計算，而后利用公式計算出表示噪聲方差；

步驟三、修正相關譜計算：將參考信號與接收到的觀測信號相關，去直流并取模后得到修正相關譜；所述步驟三具體包括：

步驟3.1、計算參考信號和觀測信號的相關序列z(n)＝x(n)y(n)；

步驟3.2、計算修正的相關譜：對相關序列作DFT變換并取模，得到相關譜Z(k)＝DFT[z(n)]，將Z(k)中直流分量濾除后取模值，得到修正的相關譜Z_m(k)；

步驟四、相關譜標準化；所述步驟四具體包括：

步驟4.1、定義隨機變量

在H₀假設下，R(k)為自由度為2的獨立同分布指數分布隨機序列，其概率密度及分布函數分別為：

其中，r表示隨機變量；

在H₁假設下，L為相關譜中峰值的個數，R(k)表示為：

由上式可見，R(k)分成兩組，其中一組是L個峰值，k₀為峰值的起點位置，即R_i,i∈(k₀,k₀+L-1)，服從自由度為2、參數的獨立非同分布的非中心卡方分布，記為R₁₁(k)，長度為N₁₁＝L，其概率密度及分布函數分別為：

其中，I₀(x)是零階修正貝塞爾函數，Q₁(a,b)是一階通用Q函數；

R(k)的另一組是除L個峰值之外的噪聲譜，服從自由度為2的獨立同分布指數分布，記為R₁₀(k)，長度為N₁₀＝N-L-1；

步驟4.2、相關譜標準化

在H₀假設下，將修正的相關譜作標準化處理，即：

其中，u表示隨機變量，N₀＝N-1，歸一化系數為：

在H₁假設下，將修正的相關譜進行分類，對于R₁₀(k)及R₁₁(k)分別通過及進行標準化處理；

其中，歸一化常數分別由下列公式計算得到：

步驟五、統計量計算；

步驟六、門限計算：在給定的虛警概率下確定門限；

步驟七、可信性判決：通過比較統計量與門限，校驗復雜調制信號處理結果。

2.如權利要求1所述的一種基于極值似然比的復雜調制信號處理結果可信性校驗方法，其特征在于：

所述步驟五中，計算統計量