专利名称：一种用于直扩系统的mpsk窄带干扰抑制方法
技术领域：
本发明涉及一种用于直扩系统的MPSK窄带干扰抑制方法，属于通信抗干扰技术领域。
背景技术：
直接序列扩频系统(简称直扩系统)具有抗干扰能力强、信息信号隐蔽、便于加 密、任意选址以及易于组网等独特优点，因此直扩技术目前在通信设备中被广泛应用，如目 前所采用的CDMA通信系统、GPS导航系统等都是直扩系统。扩频系统由于其自身的编码增 益和处理增益，具有一定的抗干扰能力，但是发射机/接收机的复杂性以及可供利用的带 宽限定了处理增益的上限，因而在遭受超出干扰容限的强干扰时，扩频系统性能会严重恶 化，使得通信质量得不到保证。MPSK(M-ary Phase Shift Keying)即M进制相移键控调制，是一种典型的数字调 制方式，它是利用载波的多种不同相位状态来表征数字信息的调制方式。例如BPSK(二相 相移键控)、QPSK(即4PSK，四相相移键控调制)等都是MPSK。MPSK窄带干扰很容易产生， 是一种常见的干扰信号，当干扰信号功率较大而且其主瓣或旁瓣进入通信信号带宽时，对 直扩系统的性能恶化明显。对MPSK窄带干扰抑制常用的处理方法一般有两大类一是利用模拟处理，对干扰 信号频率进行估计，在于扰信号频率处设置一个陷波器，对干扰信号进行抑制，实际中常用 声表面类的技术实现，需要额外增加专门的抗干扰设备，实现起来成本高。二是利用数字信 号处理的办法在频域滤波实现干扰抑制，然后将信号还原到时域。利用数字信号处理的办 法不需要加额外的装置，并且有更灵活的处理办法。专利“扩频系统中通过加窗处理消除窄带干扰的方法和装置”(专利号为 ZL200580047600. 2)所提出的扩频系统中通过加窗处理消除窄带干扰的方法和装置是利 用频域去干扰的办法去掉窄带干扰，虽有一定的效果，但是在干扰信号旁瓣泄露严重的情 况下抑制效果并不理想，例如对MPSK这种特殊的窄带干扰，只通过频域去干扰的方法所能 得到的效果不佳，因为MPSK窄带干扰信号频谱的旁瓣较高，并且旁瓣频谱分布的频段比较 宽，一般采用频域干扰抑制只能去掉干扰信号的主瓣频谱，而对旁瓣频谱无能为力，故抑制 效果不理想；若利用模拟处理的方法，同样不能滤掉MPSK窄带干扰信号频谱的大量旁瓣。 故只采用一般的滤波的方法难以减轻或消除MPSK窄带干扰信号旁瓣对直扩系统通信性能 的影响。
发明内容
在直接序列扩频系统中，为了解决传统方法难以完全抑制MPSK等旁瓣较大的窄 带干扰信号，使得误码率高，导致系统性能下降的问题，本发明提出了一种用于直扩系统的 MPSK窄带干扰抑制方法，采用频域抑制和时域抑制相结合的方法对MPSK窄带干扰进行有 效抑制，即首先在频域抑制MPSK窄带干扰信号主瓣能量，然后在时域抑制其旁瓣能量，以降低误码率，保障直接序列扩频系统的可靠通信。本发明首先利用直接序列扩频系统通信信号在傅立叶域的噪声特性和MPSK窄带 干扰信号在傅立叶域的能量聚集特性之间的差异，有效分离MPSK窄带干扰信号和有用通 信信号；其次，利用分段重叠加窗的办法，减小通信信号经快速傅里叶变换后的频谱泄漏， 尽量抑制MPSK窄带干扰信号的主瓣；最后，利用MPSK窄带干扰旁瓣频谱在时域表现为脉 冲形式的这一特点，通过在时域选择合适的门限，抑制MPSK窄带干扰信号的旁瓣，从而将 MPSK窄带干扰抑制得更加干净，以达到较好的误码率性能，保障直扩系统的可靠通信。本发明一种用于直扩系统的MPSK窄带干扰抑制方法，具体实现步骤为步骤1，将时域中长度为N、含有MPSK窄带干扰的通信信号f(x)分为k段，每段长 度为 M，各段数据记为 If1(X)^2(X), ...，fk(X)}，kM = N;步骤2，超前Μ/4处理在步骤1中的通信信号f(x)前补长度为Μ/4的数据0，去 除通信信号f(X)最后Μ/4个数据，形成新的数据a (X)，a (χ)长度仍为N，将其分为k段，每 段长度为M，各段数据记为{ai(X)，a2(X)，…，ak(x)}；步骤3，滞后M/4处理去除步骤1中的通信信号f (χ)前面Μ/4个数据，在通信信 号f(x)后补长度为Μ/4的数据0，形成新的数据b(x)，b(x)长度仍为N，将其分为k段，每 段长度为M，各段数据记为Ib1Uhb2(X),…，bk(x)}；步骤4，将步骤2、步骤3得到的2k个时域数据段Ia1 (χ)，a2 (χ)，…，ak (χ)}、Ib1 (χ)， b2(x)，…，bk(x)}分别乘以窗函数w(n)，然后做快速傅里叶变换，得到频域的数据{Α“ω)， A2(ω), ...，Α“ω)}、^1(GJ)jB2(CO), ...，Bk(co)};步骤5，将步骤4得到的2k个频域数据段(A1 (ω)，A2 (ω )，…，Ak (ω)}、(B1 (ω)， Β2(ω)，…，Bk( )}分别求模值，计算每个频域数据段的自适应门限T = Κ· θ (其中T为 门限值，θ为每个数据段在频域的均值，K为调节因子)，将频域数据段中模值超过门限T 的频域数据置为0，从而抑制干扰信号的主瓣；步骤6，将步骤5得到的经过门限处理后的2k个频域数据段做傅立叶逆变换，得到 时域数据段{a' !(ω)^' 2(ω)，-,a' k( )}、{b' ^ω)^' 2(ω)，…丄‘k( )}；步骤7，将步骤6得到的2k个时域数据段Ia' “ω)』' 2(ω), ...，a' “ω)}、 {b' !(ω), b' 2(ω)，…，b' k( )}进行去重叠处理，首先将2k个时域数据段排列为 Ia' ,(ω), b' !(ω), a' 2(ω), b' 2(ω)，…，a' k(co)，b' k ( ω )}，每个数据段长度 为M, 2k个时域数据段总长度为2kM = 2N ；再去除每个时域数据段的前M/4和后M/4个数 据，得到新的 2k 个时域数据段{a〃 ^ω)^" a" 2( ),b" 2(ω)，a" k( ), b" k( )}，记为信号f' (χ)，每个新的时域数据段长度变为M/2，2k个新的时域数据段总 长度为kM = N，信号f ‘ (χ)与处理前f (χ)信号的长度一致；步骤8，时域处理取步骤7所得到信号f' (x)的幅度均值的P倍作为时域脉冲 干扰抑制的门限G，将大于门限G的时域数据置零；其中P可以根据脉冲干扰的强弱进行调 整，典型取值为4 ；如果为了简化处理，并且在接收信号Eb/ΝΟ彡10db(l比特的信号能量比 噪声功率谱密度大于等于IOdb)的情况下，不采用时域置零的方法，而采用Ibit量化的方 法，即保证数据符号不变，正数用1表示，负数用-1表示，若数据是零，那仍然用零表示，如 果接收信号的信噪比低，即Eb/ΝΟ彡10db，则还是利用时域取门限置零的方法比较合适。本发明的理论依据和导出过程如下
(1)通信信号在频域的加窗特性对通信信号进行采样，取有限个数据进行傅立叶变换，信号截断的过程等价于将 信号乘以矩形窗函数。信号加窗是不可避免的，信号与窗函数相乘意味着总的变换是所期 望的变换与窗函数变换的卷积，如果信号的真正功率集中在一个频带内，此卷积运算将把 这个窄带功率扩展到临近的范围，造成频谱泄露。泄露效应是离散傅立叶变换所固有的，可 以采用窗函数加权技术，通过选择合适的窗函数使被加权信号在边缘比矩形窗函数圆滑而 减小了陡峭边缘所引起的旁瓣分量，达到抑制频谱泄露的目的。
以矩形窗为例，其第一旁瓣只比主瓣低13. 46dB,即旁瓣抑制度只有13. 46dB，对 于比有用信号大几十分贝的干扰来说，它的旁瓣也比信号大很多，在进行干扰抑制时，就造 成干扰消除的不彻底，或者是增大了消除的带宽范围，因而加重了对有用信号的损伤。为了 减小干扰的频谱泄露，必须采用旁瓣比较低的窗函数。但在旁瓣降低的同时，主瓣迅速变 宽，因此需要选择合适的窗函数。表1是各种窗函数的参数比较，根据直扩系统和MPSK干 扰的特点，本发明选用布莱克曼窗作为信号加权窗函数。表1各种窗函数 (2)分段重叠处理原理对接收信号进行加窗减小了边缘数据的不连续影响，但是加窗截断后会在一定程 度上扭曲信号，造成了信噪比损失。为了减小这种扭曲，可以采用重叠处理技术。采用两路信号处理通道，第一路首先 在原始信号前添加Μ/4个零，第二路的输入为原始信号Μ/4点的延迟信号，两路进行变换处 理之后，把前后各Μ/4的样点抛弃，保留中间Μ/2样点，然后将两路信号合成。这样做的目 的在于去掉由于加窗而扭曲较大的信号分量，保留损失很小的信号，当两路合成一路时，对 整个信号的扭曲就较小了。在频域的干扰抑制中，采用阈值检测方法，因为扩频信号加上高斯噪声的频谱较 为平坦，而MPSK干扰具有聚集特性，所以根据具体情况设定一个合理的阈值，对超过阈值 的数据点进行置零，该方法简单易行，适用于要求快速消除干扰影响的的场合。显然，阈值 选取的合适与否是该算法的关键。对此，设置门限为T = Κ· Θ，其中θ是接收信号在频 域的均值，K为随θ变化自适应选择的量化因子，根据性能需求选好相应的K就可以算出 干扰消除门限，完成对干扰的消除。(3) MPSK窄带干扰时域抑制原理首先以BPSK干扰为例说明时域抑制原理，由于频域干扰抑制可以认为是将超过门限的信号频谱置零，而超过门限部分频谱基本集中在干扰信号频谱主瓣内，可以认为是 将一个频域矩形窗内的干扰信号频谱置零，这个矩形窗所对应的时域冲击响应为h (η)，为s inc函数，图1中b图画出了 h(η)的主瓣部分。Pjs (n) = S(n)- S(n) Θ h(n)S(η)是窄带MPSK干扰的基带信号，如图1中a图所示。风… /^)表示S(n)与 h(n)卷积，如图1中c图所示。如图1中d图所示经过频域抑制之后，变换到时域所剩下的 没有抑制完的干扰(η)，也就是图1中a图减c图出来的结果。由上述推到过程可以看出经过频域干扰抑制，信号变回时域之后在干扰信号的 符号跳变点出现了类似于脉冲形状的干扰，如图1中d图所示，这些脉冲干扰是由于干扰信 号大量旁瓣没有被抑制掉所造成的。从物理概念上来讲，只要出现相位不连续，在频域都会 出现大量的副瓣，而频域无法将这些副瓣抑制干净，所以信号回到时域会在相位跳变点出 现类似与脉冲的干扰，干扰越强，相位跳变越大，这些脉冲干扰就会越强，而MPSK信号一般 是相位不连续的信号，必定在时域会产生不少脉冲干扰。所以如果将这些脉冲干扰置零，即 消除脉冲干扰，就可以抑制掉干扰信号的旁瓣，那么会提高直扩通信系统的抗干扰性能。(4)时域脉冲干扰抑制门限的选取方法在没有干扰的情况下，时域信号服从高斯分布，那么信号的幅度α服从瑞利分 布 其中α >0，瑞利分布的均值积累分布函数为l-exp(-x2/2a2)故当设置门限为X，那么α超过门限χ的概率为exp(-x2/2a2)所以如果用幅度的均值乘以系数P作为门限G，那么χ超过门限G的概率为门限G 1倍的均值2倍的均值3倍的均值4倍的均值 超过门限 Wexp(I) exp (- π ) exp(-—) exp (_4 π )概率=0. 4559 = 0. 0432 = 8. 5e_004 = 3. 5e_006由于时域的门限设置和信号的SNR (信噪比)密切相关，对应不同信噪比所选取的 门限就不一样，经过大量的仿真，当门限设置为信号幅度的均值的P倍时，基本不会将有用 信号置零，实际中可以结合当时的实际其情况适当的选取P的大小。其中P可以调整，根 据脉冲干扰的强弱，选择合适P，P取4是一个典型值。若存在较强的MPSK干扰时，由于频 域抑制之后会出现类似脉冲的干扰，这样会抬高门限，因此取4倍信号幅度的均值时这个 门限基本会比没有受到脉冲干扰影响位置的有用信号幅度大得多，因此有用信号不会有损 失，同时这个门限比大多数的脉冲干扰信号低得多，即因此可以抑制大部分的脉冲干扰。
针对干扰信号中未抑制彻底的副瓣能量在时域造成的影响，本文还提出了一种简化的处理办法。采用Ibit量化技术，即信号幅度大于0时用1表示，小于0时用-1表示。 该方法的目的是尽量消除时域幅度较大的脉冲干扰对扩频信号自相关特性的影响，采用了 一种类似平均化处理的办法。Ibit量化处理的方法可以明显提高信噪比，并且大大降低了 后续处理的运算量。同时，需要指出的是，Ibit量化处理会带来大概2dB左右的信噪比损 失，在实际系统中对扩频信号的捕获处理影响不大，这也是很多商用GPS采用Ibit量化的 原因。可参考文献《量化对GPS接收机捕获性能的影响》(作者是孙进芳等)可知对于GPS 接收系统，SNR(信号功率比噪声功率)大于_20dB时，基本上不会对系统产生大的影响， 换算成Eb/N0 = 10db(l比特的信号能量比噪声功率谱密度)，所以只要直扩系统接收Eb/ NO ^ 10db，本系统就可以采用Ibit量化处理时域脉冲干扰。有益效果本发明提出的一种用于直扩系统的MPSK窄带干扰抑制方法技术，相比于现有技 术其优点在于(1)提出联合频域和时域干扰抑制相结合的方法抑制MPSK窄带干扰信号，可以较 好地抑制MPSK窄带干扰主瓣和旁瓣，在不增加信号发射功率前提下有效提高直扩系统抗 MPSK干扰性能。(2)利用了直接序列扩频系统通信信号在频域的噪声特性和干扰干扰信号的能量 聚焦特性之间的差异，并使用重叠加窗的办法减轻频谱泄露，在频域充分抑制MPSK窄带干
扰信号的主瓣。(3)提出了在时域抑制MPSK窄带干扰信号旁瓣的方法，MPSK窄带干扰信号大量旁 瓣在频域无法抑制，所以在于扰抑制之后在时域会出现大量类似于脉冲的干扰，本发明提 出了抑制这些脉冲干扰的有效方法，从而抑制MPSK窄带干扰信号的旁瓣。(4)提出的方法能并行操作，可将通信信号分两路同时操作，有利于流水实现，适 合硬件平台的快速实现。因为分段重叠分成的两路信号，其处理方法一致，所以适合硬件流 水作业。(5)本发明方法计算复杂度低，易于系统实时实现。
图l-(a)是窄带BPSK干扰基带信号，(b)是矩形窗所对应的时域冲击响函数的主 瓣，即sine函数的主瓣，(c)是窄带BPSK干扰基带信号与sine函数的卷积，(d)是图(a) 减图(c)的结果；图2-本发明联合频域和时域干扰抑制原理框图；图3-本发明联合频域和时域干扰抑制具体实现框图；图4-具体实施例情况一中直扩信号叠加BPSK干扰的频谱；图5-具体实施例情况一中所述频域抑制之后信号频谱；图6-具体实施例情况一中所述频域抑制之后的信号时域波形；图7-具体实施例情况一中所述干扰信号时域波形；图8-具体实施例情况一中所述去掉脉冲干扰后信号的频谱；图9-具体实施例情况二中所述频域抑制之后信号频谱；
图10-具体实施例情况二中所述频域抑制之后的信号时域波形；图11-具体实施例情况二中所述干扰信号时域波形；图12-具体实施例情况二中所述去掉脉冲干扰后信号的频谱； 图13-具体实施例情况一中所述联合频域时域干扰抑制得的到的误码率曲线比较。
具体实施例方式根据前面“发明内容”部分中的论述，下面结合附图及具体实施例对本发明方法做 进一步详细说明，以更好地说明本发明的目的和优点。本发明提出一种用于直扩系统的MPSK窄带干扰抑制方法，其原理框图见图2和图 3，具体实施例如下假设主瓣带宽为16KHZ的BPSK干扰叠加在带宽为1. 25MHZ的直扩通信信号f (χ) 上，JSR(干扰信功率比有用信号功率)=45db，Eb/N0 = 10db，扩频码长2048，扩频倍数为 2048。Sf(X)为混合有1个BPSK干扰的直扩通信信号，长度N为2048点。按照本发明内容所述，通过如下步骤实现直扩系统的BPSK窄带干扰抑制步骤1，将长度为2048的通信信号f (χ)分为4段，每段长度M为512，各段数据记 为If1 (x),f2(x),…，f4(x)}，存放在FPGA (现场可编程门阵列)的RAM(随机存取存储器) 内；步骤2，在步骤1中的通信信号f (χ)前补长度为M/4 = 128的数据0，去除f (χ) 信号最后128个数据，形成新的数据a (χ)，a (χ)信号长度仍为2048，分为4段，每段长度为 512，各段数据记为(x),a2(x),…，a4(x)}； 步骤3，去除步骤1中的通信信号f (χ)前面128个数据，在通信信号f (χ)后补长度 为128个数据0，形成新的数据b (χ)，b (χ)信号长度仍为2048，分为4段，每段长度为512， 各段数据记为Ib1(X)^2(X),…，b4(x)};步骤4，将步骤2、步骤3得到的8个时域数据段Ia1 (χ)，a2 (χ)，…，a4 (χ)}、Ib1 (χ)， b2(x),…，b4(x)}分别乘布莱克曼窗函数w (η)，调用FPGA内的IP核做快速傅立叶变换得 到频域数据队(《)，4( )，…，Α4(ω)}、{BjcohBjco)，…，Β4(ω)};步骤5，对步骤4得到的8个频域数据段{、(《)，A2 (ω)，…，Α4(ω)}、{Β^ω), Β2(ω)，…，Β4(ω)}，分别求模值，计算每个频域数据段的自适应门限T = Κ· θ (其中T为 门限值，θ为均值，K为调节因子)，将频域数据段中超过门限的数据置为零，抑制BPSK的 主瓣干扰信号；步骤6，将步骤5中过门限处理后的8个数据段进行傅立叶逆变换得到时域数据
Ia'2(ω)，…，a' 4(co)}、{b' ^ω)^' 2(ω)，...，b' 4(ω)}；步骤7，对步骤6得到的8个时域数据段Ia' ^ω), a' 2(ω)，…，a' 4(ω)}、 {b' i( )，b' 2( ),...,b' 4(ω)}进行去重叠处理，首先将8个数据段排列为Ia' !(ω), b' !(ω), a' 2(ω), b' 2(ω)，…，a' 4(ω), b' 4 (ω)}，每个数据段长度为 512，8 个 数据段总长度为4096 ；再去除每个数据段的前128和后128个数据，得到新的8个数据段 {a" “cohb" “ω)，^ 2(co)，b" 2(ω)，...，a" 4(co)，b" 4(ω)}，记为 f' (χ)，每个 数据段长度变为256，8个数据段总长度为2048，f ‘ (χ)与处理前f (χ)信号的长度一致；
步骤8，将步骤7得到的信号取幅度均值的4倍时(这个系数可以适当调整)作为 时域脉冲干扰抑制的门限G，将大于门限G的数据置零。BPSK干扰信号加在直扩通信信号上一般有两种情况情况一 BPSK干扰信号中心频率刚好与直扩通信信号中心频率对齐。直扩信号叠加BPSK干扰的频谱如图4所示。经频域抑制之后的频谱图如图5所 示，可见大部分旁瓣没有被滤掉。经频域抑制之后的信号时域波形如图6所示。对应的干 扰信号时域波形如图7所示，可见图6中脉冲干扰出现在图7中相位跳变点的位置。通过 本发明的方法在时域中去掉脉冲干扰后，信号的频谱如图8所示可见抑制掉了干扰信号的 大部分旁瓣频谱。情况二 若BPSK干扰信号中心频率和直扩通信信号中心频率有一定频偏，设频偏 为 16kHZ。经频域抑制之后的信号频谱图如图9所示，可见大部分旁瓣没有被滤掉。信号在 频域抑制后变换到时域的波形如图10所示，对应的干扰信号时域波形如图11所示，可见图 10中的脉冲出现在图11中相位跳变点处出现。通过本发明的方法在时域中去掉脉冲干扰 后，信号的频谱图如图12所示，可见抑制掉了干扰信号的大部分旁瓣频谱。从图5和图9可见频域干扰抑制去掉了干扰的主瓣内的大部分频谱，而干扰信号 大量的旁瓣(即干扰信号的高频部分)没有抑制掉，将频域抑制后的信号变换到时域时，在 长度为2048的扩频码的一个周期内出现了多个类似于脉冲的干扰，而这些脉冲刚好对应 于干扰信号的相位跳变点，所以时域取门限将这些脉冲干扰致零，会大大减小干扰信号的 功率，提高直扩系统的抗干扰性能。图13为情况一中BPSK干扰信号中心频率刚好与直扩通信信号中心频率对齐，干 扰主瓣带宽为16kHZ，按照上述过程所得的到的误码率曲线，可见在干扰比较弱的情况下， 采用频域和时域的联合干扰抑制方法与单一的频域干扰抑制性能相差不大；但在干扰越 强的状态下，采用频域和时域的联合干扰抑制方法比单一的频域干扰抑制性能有显著的提高。以上的具体实施方式
描述，对发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保 护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本 发明的保护范围之内。
权利要求
一种用于直扩系统的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，则利用时域取门限置零的方法。
2.根据权利要求1所述的一种用于直扩系统的MPSK窄带干扰抑制方法，，其特征在于 步骤8所述信号f' (x)的幅度均值的P倍根据脉冲干扰的强弱进行调整，P的典型取值为 4。
全文摘要
本发明涉及一种用于直扩系统的MPSK窄带干扰抑制方法，属于通信抗干扰技术领域。本发明采用频域和时域相结合的方法分别抑制MPSK窄带干扰信号的主瓣和旁瓣。首先利用直扩通信信号在频域的噪声特性和MPSK窄带干扰信号能量聚焦特性之间的差异，采用分段重叠处理的办法，尽可能多抑制MPSK窄带干扰信号主瓣的同时保障尽可能减少有用信号能量的损失；其次，利用MPSK窄带旁瓣频谱在时域表现为脉冲形式的这一特点在时域抑制MPSK窄带干扰信号的旁瓣，从而有效分离干扰和有用信号，以达到较好的误码率性能，保障直扩系统的可靠通信。另外，本发明方法计算复杂度低，易于系统实时实现。
文档编号H04L27/18GK101841349SQ20101014379
公开日2010年9月22日 申请日期2010年4月12日 优先权日2010年4月12日
发明者吴葵, 郇浩, 陶然, 黄克武 申请人:北京理工大学