[發(fā)明專利]基于菱形積分雙譜的通信電臺個體識別方法及裝置

無效

• 申請?zhí)枺?/td>: 201110354259.2
• 申請日: 2011-11-10
• 公開（公告）號: CN102510314A
• 公開（公告）日: 2012-06-20
• 發(fā)明（設(shè)計）人: 徐書華;徐麗娜;徐爭光;莫益軍;鄒寧
• 申請（專利權(quán)）人: 華中科技大學(xué)
• 主分類號: H04B17/00
• 分類號: H04B17/00
• 代理公司: 北京市德權(quán)律師事務(wù)所 11302
• 代理人: 劉麗君

1.一種基于菱形積分雙譜的通信電臺個體識別方法，其特征在于，包括：

對通信電臺個體信號進(jìn)行雙譜特征分析得到菱形積分雙譜RHIB；

對所述菱形積分雙譜RHIB的特征按照Fisher可分離度準(zhǔn)則進(jìn)行優(yōu)化，并進(jìn)行通信電臺個體特征融合；

采用一對一的支撐矢量機(jī)作為模式分類器對通信電臺個體進(jìn)行識別。

2.根據(jù)權(quán)利要求1所述的識別方法，其特征在于，對通信電臺個體信號進(jìn)行雙譜特征分析得到菱形積分雙譜RHIB包括：

計算通信電臺個體信號{x⁽ⁱ⁾(n)}的傅立葉變換矢量X⁽ⁱ⁾(ω)，其中，所述頻率ω∈Ω，所述T為信號長度，所述Z表示整數(shù)域；

計算通信電臺個體信號{x⁽ⁱ⁾(n)}的直接雙譜B⁽ⁱ⁾(ω)，所述B⁽ⁱ⁾(ω)＝B⁽ⁱ⁾(ω₁，ω₂)＝X⁽ⁱ⁾(ω₁)X⁽ⁱ⁾(ω₂)X⁽ⁱ⁾(ω₁+ω₂)，所述ω＝(ω₁，ω₂)；

計算通信電臺個體信號{x⁽ⁱ⁾(n)}的菱形積分雙譜RHIB⁽ⁱ⁾(ω)，所述，所述B(ω₁，ω₂)表示菱形路徑P_r上通信信號的雙譜數(shù)值，所述RHIB⁽ⁱ⁾(ω)是第i類通信電臺信號樣本在頻率ω處的菱形積分雙譜。

3.根據(jù)權(quán)利要求2所述的識別方法，其特征在于，所述對所述菱形積分雙譜RHIB的特征按照Fisher可分離度準(zhǔn)則進(jìn)行優(yōu)化，并進(jìn)行通信電臺個體特征融合包括：

假定是第i類通信電臺信號樣本在頻率ω處的菱形積分雙譜集合，是第j類通信電臺信號樣本在頻率ω處的菱形積分雙譜集合；其中，所述i和j表示通信電臺信號的類別，所述中的k表示所監(jiān)測到的第i類通信電臺的第k個監(jiān)測通信信號樣本，所述中的k表示所監(jiān)測到的第j類通信電臺的第k個監(jiān)測通信信號樣本，所述N_i是第i類通信電臺監(jiān)測信號樣本的個數(shù)，所述N_j是第j類通信電臺監(jiān)測信號樣本的個數(shù)；

將第i類和第j類個體信號特征間的Fisher可分離度按照如下公式進(jìn)行計算：

m(i,j)(ω)=Σl=i,jp(l)[meank(RHIBk(l)(ω))-meanl[meank(RHIBk(l)(ω))]]2Σl=i,jp(l)vark(RHIBk(l)(ω)),]]>i≠j，

其中，所述m^(i，j)(ω)表示第i類和第j類個體信號特征間在頻率ω處的可分離度；所述p^(l)為第l類電臺通信信號發(fā)生的先驗概率；所述mean_k(RHIB_k^(l)(ω))代表第l類電臺通信信號在頻率ω處所有樣本菱形積分雙譜的均值；所述var_k(RHIB_k^(l)(ω))代表第l類電臺通信信號在頻率ω處所有樣本菱形積分雙譜的b方差；所述mean_l[mean_k(RHIB_k^(l)(ω))]表示所有類型信號全體樣本在頻率ω處的菱形積分雙譜均值；

對于所監(jiān)測到的第i類通信電臺的所有N_i個信號樣本，根據(jù)第i類和其他類別個體信號特征間的Fisher可分離度篩選出第i類電臺的菱形積分雙譜信號特征，并標(biāo)記為RHIB_i(l)；

采用快速傅里葉變換方法對第i類通信信號個體進(jìn)行分析，得到通信電臺個體的發(fā)射載頻特征，標(biāo)記為f_c(i)；

對所述菱形積分雙譜信號特征RHIB_i(l)和所得到的通信電臺個體的發(fā)射載頻特征f_c(i)按照OSF_i＝{f_c(i)，RHIB_i(l)}進(jìn)行特征融合，然后對所述OSF_i進(jìn)行歸一化處理，從而得到第i類通信電臺的個體信號特征SF_i。