[發明專利]一種用于直擴系統的MPSK窄帶干擾抑制方法

無效

• 申請號: 201010143796.8
• 申請日: 2010-04-12
• 公開（公告）號: CN101841349A
• 公開（公告）日: 2010-09-22
• 發明（設計）人: 陶然;吳葵;郇浩;黃克武
• 申請（專利權）人: 北京理工大學
• 主分類號: H04B1/707
• 分類號: H04B1/707;H04B1/10;H04L25/02;H04L27/18
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.一種用于直擴系統的MPSK窄帶干擾抑制方法，其特征在于是通過以下步驟實現的：

步驟1，將時域中長度為N，含有MPSK窄帶干擾的通信信號f(x)分為k段，每段長度為M，各段數據記為{f₁(x)，f₂(x)，...，f_k(x)}，kM＝N；

步驟2，超前M/4處理：在步驟1中的通信信號f(x)前補長度為M/4的數據0，去除通信信號f(x)最后M/4個數據，形成新的數據a(x)，a(x)長度仍為N，將其分為k段，每段長度為M，各段數據記為{a₁(x)，a₂(x)，...，a_k(x)}；

步驟3，滯后M/4處理：去除步驟1中的通信信號f(x)前面M/4個數據，在通信信號f(x)后補長度為M/4的數據0，形成新的數據b(x)，b(x)長度仍為N，將其分為k段，每段長度為M，各段數據記為{b₁(x)，b₂(x)，...，b_k(x)}；

步驟4，將步驟2、步驟3得到的2k個時域數據段{a₁(x)，a₂(x)，...，a_k(x)}、{b₁(x)，b₂(x)，...，b_k(x)}分別乘以窗函數w(n)，然后做快速傅里葉變換，得到頻域的數據{A₁(ω)，A₂(ω)，...，A_k(ω)}、{B₁(ω)，B₂(ω)，...，B_k(ω)}；

步驟5，將步驟4得到的2k個頻域數據段{A₁(ω)，A₂(ω)，...，A_k(ω)}、{B₁(ω)，B₂(ω)，...，B_k(ω)}分別求模值，計算每個頻域數據段的自適應門限T＝K·θ(其中T為門限值，θ為每個數據段在頻域的均值，K為調節因子)，將頻域數據段中模值超過門限T的頻域數據置為0，從而抑制干擾信號的主瓣；

步驟6，將步驟5得到的經過門限處理后的2k個頻域數據段做傅立葉逆變換，得到時域數據段{a′₁(ω)，a′₂(ω)，...，a′_k(ω)}、{b′₁(ω)，b′₂(ω)，...，b′_k(ω)}；

步驟7，將步驟6得到的2k個時域數據段{a′₁(ω)，a′₂(ω)，...，a′_k(ω)}、{b′₁(ω)，b′₂(ω)，...，b′_k(ω)}進行去重疊處理，首先將2k個時域數據段排列為{a′₁(ω)，b′₁(ω)，a′₂(ω)，b′₂(ω)，...，a′_k(ω)，b′_k(ω)}，每個數據段長度為M，2k個時域數據段總長度為2kM＝2N；再去除每個時域數據段的前M/4和后M/4個數據，得到新的2k個時域數據段{a″₁(ω)，b″₁ω)，a″₂(ω)，b″₂(ω)，...，a″_k(ω)，b″_k(ω)}，記為信號f′(x)，每個新的時域數據段長度變為M/2，2k個新的時域數據段總長度為kM＝N，信號f′(x)與處理前f(x)信號的長度一致；

步驟8，時域處理：取步驟7所得到信號f≥(x)的幅度均值的P倍作為時域脈沖干擾抑制的門限G，將大于門限G的時域數據置零；如果為了簡化處理，并且在接收信號Eb/N0≥10db即1比特的信號能量比噪聲功率譜密度大于等于10db的情況下，不采用時域置零的方法，而采用1bit量化的方法，即保證數據符號不變，正數用1表示，負數用-1表示，若數據是零，仍然用零表示，如果接收信號的信噪比低即Eb/N0≤10db，則利用時域取門限置零的方法。

2.根據權利要求1所述的一種用于直擴系統的MPSK窄帶干擾抑制方法，，其特征在于：步驟8所述信號f′(x)的幅度均值的P倍根據脈沖干擾的強弱進行調整，P的典型取值為4。