[发明专利]一种用于直扩系统的MPSK窄带干扰抑制方法

摘要

本发明涉及一种用于直扩系统的MPSK窄带干扰抑制方法，属于通信抗干扰技术领域。本发明采用频域和时域相结合的方法分别抑制MPSK窄带干扰信号的主瓣和旁瓣。首先利用直扩通信信号在频域的噪声特性和MPSK窄带干扰信号能量聚焦特性之间的差异，采用分段重叠处理的办法，尽可能多抑制MPSK窄带干扰信号主瓣的同时保障尽可能减少有用信号能量的损失；其次，利用MPSK窄带旁瓣频谱在时域表现为脉冲形式的这一特点在时域抑制MPSK窄带干扰信号的旁瓣，从而有效分离干扰和有用信号，以达到较好的误码率性能，保障直扩系统的可靠通信。另外，本发明方法计算复杂度低，易于系统实时实现。

主权项

一种用于直扩系统的MPSK窄带干扰抑制方法，其特征在于是通过以下步骤实现的：步骤1，将时域中长度为N，含有MPSK窄带干扰的通信信号f(x)分为k段，每段长度为M，各段数据记为{f1(x)，f2(x)，...，fk(x)}，kM＝N；步骤2，超前M/4处理：在步骤1中的通信信号f(x)前补长度为M/4的数据0，去除通信信号f(x)最后M/4个数据，形成新的数据a(x)，a(x)长度仍为N，将其分为k段，每段长度为M，各段数据记为{a1(x)，a2(x)，...，ak(x)}；步骤3，滞后M/4处理：去除步骤1中的通信信号f(x)前面M/4个数据，在通信信号f(x)后补长度为M/4的数据0，形成新的数据b(x)，b(x)长度仍为N，将其分为k段，每段长度为M，各段数据记为{b1(x)，b2(x)，...，bk(x)}；步骤4，将步骤2、步骤3得到的2k个时域数据段{a1(x)，a2(x)，...，ak(x)}、{b1(x)，b2(x)，...，bk(x)}分别乘以窗函数w(n)，然后做快速傅里叶变换，得到频域的数据{A1(ω)，A2(ω)，...，Ak(ω)}、{B1(ω)，B2(ω)，...，Bk(ω)}；步骤5，将步骤4得到的2k个频域数据段{A1(ω)，A2(ω)，...，Ak(ω)}、{B1(ω)，B2(ω)，...，Bk(ω)}分别求模值，计算每个频域数据段的自适应门限T＝K·θ(其中T为门限值，θ为每个数据段在频域的均值，K为调节因子)，将频域数据段中模值超过门限T的频域数据置为0，从而抑制干扰信号的主瓣；步骤6，将步骤5得到的经过门限处理后的2k个频域数据段做傅立叶逆变换，得到时域数据段{a′1(ω)，a′2(ω)，...，a′k(ω)}、{b′1(ω)，b′2(ω)，...，b′k(ω)}；步骤7，将步骤6得到的2k个时域数据段{a′1(ω)，a′2(ω)，...，a′k(ω)}、{b′1(ω)，b′2(ω)，...，b′k(ω)}进行去重叠处理，首先将2k个时域数据段排列为{a′1(ω)，b′1(ω)，a′2(ω)，b′2(ω)，...，a′k(ω)，b′k(ω)}，每个数据段长度为M，2k个时域数据段总长度为2kM＝2N；再去除每个时域数据段的前M/4和后M/4个数据，得到新的2k个时域数据段{a″1(ω)，b″1ω)，a″2(ω)，b″2(ω)，...，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，则利用时域取门限置零的方法。