一种宽带bpsk信号的单音和多音干扰去噪方法
【专利摘要】本发明公开的一种宽带BPSK信号的单音和多音干扰去噪方法，包含以下步骤：接收机接收在某一跳时域上的数字信号；对时域接收信号进行傅立叶变换，获得频域信号；通过数字滤波器对频域信号的绝对值进行滤波；对滤波后的信号进行分帧，并分别计算每个子帧信号的能量大小；根据滤波后其子帧信号的能量分布情况，评估出接收信号中受干扰影响最小的频域区间；基于该段信号计算出一个幅度阈值，然后在频域上对接收信号进行去噪处理。本发明的方法，能够有效降低接收信号的解调错误率。
【专利说明】
一种宽带BPSK信号的单音和多音干扰去噪方法
技术领域
[0001] 本发明涉及数字信号处理领域，特别涉及一种宽带BPSK信号的单音和多音干扰去 噪方法。
【背景技术】
[0002] 通信对抗作为现代战争的重要组成部分，在战争进程中的作用越来越大。随着数 据传输量的日益扩大，通信信号的带宽正在不断拓展，因此跳频宽带信号正日益成为战术 通信采用的一种关键波形技术，因此研究相应的信号干扰去噪算法显得十分的重要。
[0003] 单音干扰也称为点频干扰，干扰信号在一个频率上发射，因此干扰信号是一个单 频连续波音调。由于信号频带非常窄，因此当干扰频率落在接收机频段内时，干扰的能量可 几乎全部进入接收机，可以形成对有用信号的有效干扰。
[0004] 多音干扰也可称为多频干扰，这种信号本质上是单音干扰的简单叠加。由于单音 干扰对于信道中信号的影响只能破坏其某一个频率点上的频谱，所以造成的误码率提升效 果有限，因此可以在同一时间发送多个单音干扰破坏有用信号多个频域点上的波形，既增 加了干扰概率，又使其在接收端不能正确接收发送的信号，进一步提高了误码率。多音干扰 可以说在单音干扰的基础上同时发送多个单音，尽管发送功率要求比较高，但也达到了有 效提高有用信号误码率的目的。
[0005] 高速跳频机制是对抗跟踪式干扰的一种重要手段，但如果通信信号的带宽较大， 干扰方在干扰功率有限的前提下，可以在一个非常宽的通信频段上采用间隔较大的梳状多 音阻塞式干扰方式来对通信信号进行干扰，一旦有一个及以上的干扰频点信号落在通信信 号的带宽内，就可以有效的对该信号进行干扰，提高接收方的信号误码率。
【发明内容】
[0006] 本发明的目的在于克服现有技术的缺点与不足，提供一种宽带BPSK信号的单音和 多音干扰去噪方法。
[0007] 本发明的目的通过以下的技术方案实现：
[0008] 一种宽带BPSK信号的单音和多音干扰去噪方法，包含以下步骤：
[0009] 接收机接收在某一跳时域上的数字信号；
[0010]对时域接收信号进行傅立叶变换，获得频域信号；
[0011] 通过数字滤波器对频域信号的绝对值进行滤波；
[0012] 对滤波后的信号进行分帧，并分别计算每个子帧信号的能量大小；
[0013] 根据滤波后其子帧信号的能量分布情况，评估出接收信号中受干扰影响最小的频 域区间；
[0014] 基于该段信号计算出一个幅度阈值，然后在频域上对接收信号进行去噪处理。
[0015] 所述宽带BPSK信号的单音和多音干扰去噪方法，具体包含以下步骤：
[0016]设接收机在某一跳时域上的数字接收信号为
[0017] X[n] = S[n]+W[n] + j[n] ,n=l ,2, . . . ,N (I)
[0018] 其中S[n]为含N个符号的BPSK信号，w[n]为高斯白噪声信号，j[n]为单音或多音干扰 信号；
[0019] 对时域接收信号进行傅立叶变换，获得频域信号：
[0020] X[n]=FFT(x[n]) ,n = l,2,...,N (2)
[0021 ]通过3个数字滤波器对频域信号X[n]的绝对值I X[n] I进行滤波，其滤波器脉冲响应 数列依次为
[0022]
[0023]
[0024]
[0025] 设滤波后的3个子信号依次为Xl[m]、X2[m]和X3[m]，m=l，2,. . .，N+2,对这3个子信号 同时进行分帧，每个子帧含64个符号样点，相邻子帧之间重叠3/4个帧，即48个符号样点； [0026]分别计算每个子帧的能量大小，记为？1 [15]、？2[1{]和?3[1{]氺=1，2，...，1(;其中1(为子 帧的数量i
良示上取整；
[0027] 由于3个滤波器的脉冲响应序列长度都只有3,根据数字滤波器的工作原理知，对 一个长度为N的数字序列滤波后的子信号长度将为N+2;
[0028]
[0029]
[0030] 设D= [dhcb，...，dK]，则对里面的元素重复16次，然后在队列尾部增加48个dK，得 到如下长度为N的新队列
[0031]
[0032] 再对队列DD进行平滑处理得到一个新队列DDD，具体算法如下：[0nWI
[0034]基于队列DDD，以64个为单位，计算这些元素的和，得到 (8)
[0035] G[r]=DDD[64Xr-63]+DDD[64Xr-62]+. . .+DDD[64Xr] , r = I , . . . ,N/64 (9)
[0036]
[0037] 设向量G[r]中的第rr个的元素值最小，则认为在频域信号X[n]中，区间[64 Xrr-63， 64Xrr-62,. . .，64Xrr]上的信号受到的频域干扰影响最小；
[0038] 对频域信号X[n]在该区间上的信号[X[64Xrr-63]，X[64Xrr-62]，· · ·，X[64Xrr]]以幅度从大 到小进行排序，挑选幅度值排名第4的元素的幅度值作为一个阈值，记为C;令
[0039] Δ [n] = sign( |X[n] | ~C) (10)
[0040]
[0041]其中函数sign定义为
[0042]
[0043] 然后对频域信号X[n]进行如下去噪处理： 「00441
[0047]本发明与现有技术相比，具有如下优点和有益效果：
[0048] 理论分析和仿真可以证明，对于任何一段数字信号Y，分别通过上述的3个滤波器 后，3个滤波信号的能量的和将等于原信号的能量，SP
[0049]
[0050] 其中符号?代表卷积运算；
[00511同时对于BPSK信号s[n]而言，如果它仅仅受到加性高斯白噪声信号w[n]的影响，即 接收信号为X[n] = s[n]+w[n]，则在任意信噪比值下，对接收信号进行频域变换后X = FFT(X)取 绝对值|X|，对该信号进行滤波，然后以64个符号为长度对三个滤波信号进行分帧，经过大 量仿真统计后发现，相同位置上的每个子帧的能量分布比例情况非常稳定，特征为：
[0052]
[0053] 其中符号E(.)代表随机变量的均值。
[0054]针对宽带BPSK调制信号的接收机，本发明根据长期的仿真和试验测试结果，针对 单音干扰和多音干扰提供了一种新型的去噪算法，它把接收的频域信号取绝对值后，利用 一个含3个都仅含3个非零元素的响应系数的数字滤波器组成的，可以保持信号能量和保持 不变的滤波器组进行滤波，再根据滤波后其子信号能量分布情况，评估出接收信号中受干 扰影响最小的频域区间。并基于该段信号计算出一个幅度阈值，然后在频域上对接收信号 进行去噪处理，以此大幅度减少单音信号的干扰能力，从而提升信号的接收性能。
【附图说明】
[0055]图1是本发明所述的一种宽带BPSK信号的单音和多音干扰去噪方法的流程图。 [0056]图2是长度为2048个符号的BPSK信号、高斯白噪声、1个单音信号和混合信号的频 域仿真结果图。
[0057]图3是长度为2048个符号的BPSK信号、高斯白噪声、2个单音信号和混合信号的频 域仿真结果图。
[0058]图4是受1个单音干扰的未去噪信号、去噪后信号和OdB无干扰信号的BER曲线图。 [0059]图5是受2个单音干扰的未去噪信号、去噪后信号和OdB无干扰信号的BER曲线图。
【具体实施方式】
[0060]下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限 于此。
[0061 ] 一种宽带BPSK信号的单音和多音干扰去噪方法，包含以下步骤：
[0062]接收机接收在某一跳时域上的数字信号；
[0063] 对时域接收信号进行傅立叶变换，获得频域信号；
[0064] 通过数字滤波器对频域信号的绝对值进行滤波；
[0065] 对滤波后的信号进行分帧，并分别计算每个子帧信号的能量大小；
[0066] 根据滤波后其子帧信号的能量分布情况，评估出接收信号中受干扰影响最小的频 域区间；
[0067] 基于该段信号计算出一个幅度阈值，然后在频域上对接收信号进行去噪处理。 [0068]所述宽带BPSK信号的单音和多音干扰去噪方法，具体包含以下步骤：
[0069] 设接收机在某一跳时域上的数字接收信号为
[0070] X[n] = S[n]+W[n] + j[n] ,n=l ,2, . . . ,N (I)
[0071] 其中s[n]为含N个符号的BPSK信号，w[n]为高斯白噪声信号，j[n]为单音或多音干扰 信号；
[0072] 对时域接收信号进行傅立叶变换，获得频域信号：
[0073] X[n]=FFT(x[n]) ,n = l,2,...,N (2)
[0074] 通过3个数字滤波器对频域信号X[n]的绝对值|X[n] I进行滤波，其滤波器脉冲响应 数列依次为
[0075]
[0076]
[0077]
[0078]设滤波后的3个子信号依次为Xl[m]、X2[m]和X3[m]，m=l，2,. . .，N+2,对这3个子信号 同时进行分帧，每个子帧含64个符号样点，相邻子帧之间重叠3/4个帧，即48个符号样点； [0079]分别计算每个子帧的能量大小，记为？1 [15]、？2[1{]和?3[1{]氺=1，2，...，1(;其中1(为子 帧的数量且
表示上取整；
[0080] 由于3个滤波器的脉冲响应序列长度都只有3,根据数字滤波器的工作原理可知， 对一个长度为N的数字序列滤波后的子信号长度将为N+2;
[0081]
[0082]
[0083] 设D= [dhcb，...，dK]，则对里面的元素重复16次，然后在队列尾部增加48个dK，得 到如下长度为N的新队列 [0084]
[0085] 再对队列DD进行平滑处理得到一个新队列DDD，具体算法如下：
[0086]
[0087]基于队列DDD，以64个为单位，计算这些元素的和，得到
[0088] G[r]=DDD[64Xr-63]+DDD[64Xr-62]+. . .+DDD[64Xr] , r = I , . . . ,N/64 (9)
[0089]
[0090] 设向量G[r]中的第rr个的元素值最小，则认为在频域信号X[n]中，区间[64 Xrr-63， 64Xrr-62,. . .，64Xrr]上的信号受到的频域干扰影响最小；
[0091] 对频域信号X[n]在该区间上的信号[X[64Xrr-63]，X[64Xrr-62]，…，X[64Xrr]]以幅度从大 到小进行排序，挑选幅度值排名第4的元素的幅度值作为一个阈值，记为C;令
[0092] Δ [n] = sign( |X[n] | ~C) (10)
[0093]
[0094]其中函数sign定义为 [0095]
[0096] 然后对频域信号X[n]进行如下去噪处理：
[0097]
[0100] 具体过程如图1所示，理论分析和仿真可以证明，对于任何一段数字信号Υ，分别通 过上述的3个滤波器后，3个滤波信号的能量的和将等于原信号的能量，即
[0101]
[0102] 其中符号發代表卷积运算；
[0103]同时对于BPSK信号s[n]而言，如果它仅仅受到加性高斯白噪声信号w[n]的影响，即 接收信号为X[n] = S[n]+W[n]，则在任意信噪比值下，对接收信号进行频域变换后X = FFT(X)取 绝对值Ixl，对该信号进行滤波，然后以64个符号为长度对三个滤波信号进行分帧，经过大 量仿真统计后发现，相同位置上的每个子帧的能量分布比例情况非常稳定，特征为
[0104]
[0105] 其中符号E(.)代表随机变量的均值。
[0106] 图2分别给出了一段符号长度为2048个，功率归一化的纯BPSK信号的频域信号、功 率为1的高斯白噪声信号的频域信号、功率为100的随机相位单音干扰的频域信号以及三种 信号相加后的频域信号；图3分别给出了一段符号长度为2048个，功率归一化的纯BPSK信号 的频域信号、功率为1的高斯白噪声信号的频域信号、功率为100的随机相位单音干扰的频 域信号以及三种信号相加后的频域信号；图4则分别给出了受1个随机相位单音干扰且包含 和BPSK信号功率相同的白噪声的未去噪信号，去噪处理后的信号在不同的信干比值下的 BER性能曲线，同时给出了仅包含和BPSK信号功率相同的白噪声的信号的BER性能曲线进行 比较；图5则分别给出了受2个随机相位单音干扰且包含和BPSK信号功率相同的白噪声的未 去噪信号，去噪处理后的信号在不同的信干比值下的BER性能曲线，同时给出了仅包含和 BPSK信号功率相同的白噪声的信号的BER性能曲线进行比较。
[0107] 从图4和图5的结果可以看出，由于单音信号的带宽非常窄，一旦落入通信信号的 带宽内，几乎可以把所有的干扰功率都注入到BPSK接收信号中，导致信号BER性能急剧恶 化;而上述提供的去噪算法利用了滤波器组对单音干扰信号频率附件的能量比例情况发现 明显变化的特性，可以非常精准的评估出一个合适的频域信号的幅度去噪阈值，然后进行 接收信号的频域信号连续大于阈值的频点区域或幅度值过大的单频点进行消幅处理，从而 在频域上删除单音或多音信号的干扰能量，从而有效降低接收信号的解调错误率。
[0108] 以上对本发明所提供的一种针对BPSK信号接收机的信号单音或多音去噪算法进 行了详细介绍，本文中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实 施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术 人员，依据本发明的思想，在【具体实施方式】及应用范围上均会有改变之处。综上所述，本说 明书内容不应理解为对本发明的限制。
【主权项】
1. 一种宽带BPSK信号的单音和多音干扰去噪方法，其特征在于，包含W下步骤： 接收机接收在某一跳时域上的数字信号； 对时域接收信号进行傅立叶变换，获得频域信号； 通过数字滤波器对频域信号的绝对值进行滤波； 对滤波后的信号进行分帖，并分别计算每个子帖信号的能量大小； 根据滤波后其子帖信号的能量分布情况，评估出接收信号中受干扰影响最小的频域区 间； 基于该段信号计算出一个幅度阔值，然后在频域上对接收信号进行去噪处理。2. 根据权利要求1所述宽带BPSK信号的单音和多音干扰去噪方法，其特征在于，具体包 含W下步骤： 设接收机在某一跳时域上的数字接收信号为 X[n] = S[n]+W[n] + j[n] ,H=I ,2, ... ,N (1) 其中SW为含N个符号的BPSK信号，WW为高斯白噪声信号，j[n]为单音或多音干扰信号； 对时域接收信号进行傅立叶变换，获得频域信号： X[n]=FFT(x[n]) ,n=l ,2,. . . ,N (2) 通过3个数字滤波器对频域信号X[n]的绝对值|X[n] I进行滤波，其滤波器脉冲响应数列 依次为设滤波后的3个子信号依次为Xl[m]、X2[m]和X3[m] ,m=l ,2, . . . ,N+2,对运3个子信号同时 进行分帖，每个子帖含64个符号样点，相邻子帖之间重叠3/4个帖，即48个符号样点； 分别计算每个子帖的能量大小，记为Pl[k]、P2[k]和?3山，4=1，2，...，1(;其中1(为子帖的 数量_符号表示上取整； 由于3个滤波器的脉冲响应序列长度都只有3,根据数字滤波器的工作原理知，对一个 长度为N的数字序列滤波后的子信号长度将为化2;化U=LCii,02,...,加」，训化I生凹的凡余里巧io化，'说；化化I)y\yu;毛巧W旨加，1#至Ij女日 下长度为N的新队列 a) \ t t' (8) 基于队列D孤，W64个为单位，计算运些元素的和，得到 G[r]=孤D[64Xr-63]+D孤[64Xr-62]+. . .+D孤[64Xr] ,r=l , . . . ,N/64 (9) 设向量G[r]中的第rr个的元素值最小，则认为在频域信号X[n]中，区间[64 X rr-63，64 X rr-62,. . .，64Xrr]上的信号受到的频域干扰影响最小； 对频域f目号Xw在该区间上的f目号[X[64Xrr-63] ,X[64Xrr-62] , . . . ,X[64Xrr]] W幅度从大到小 进行排序，挑选幅度值排名第4的元素的幅度值作为一个阔值，记为C;令 A [n] = si即（|x[n] |-C) (10) 其中函数sign定义为 然后对频域信号x[n]进;fjyu r去^采^巧：
【文档编号】H04L25/03GK105915484SQ201610236747
【公开日】2016年8月31日
【申请日】2016年4月15日
【发明人】谢映海, 郑晨熹
【申请人】广州海格通信集团股份有限公司