本发明属于通信技术领域，特别涉及一种基于畸变功率去除的窄带干扰检测方法及其通信装置。
背景技术：
通信设备和系统生存于纷繁复杂的电磁环境中，随着电子设备间的密度与日俱增，导致所占用的电磁频谱越来越宽，所传输的信息量越来越大，面临的电磁干扰也日益严重。
从电磁干扰信号频谱宽度，可以将电磁干扰分为宽带干扰源和窄带干扰源。在无线和有线通信系统中，远小于系统带宽的频带就称之为窄带干扰源。在这些干扰中，高功率的窄带干扰已经成为破坏通信系统顽存性最主要的因素之一，比如在无线局域网中，在相同频点上工作的蓝牙设备会对干扰范围内的wifi设备产生窄带干扰。在类如网线或电力载波通信的有线通信系统中，短波无线电会带来射频上的窄带干扰。
窄带干扰通常表征为在相对孤立的窄频段上具有的能量远高于传输信号的平均能量。正交频分复用(ofdm)由于具有较高的频谱效率和良好的抗多径能力，被广泛应用在实际通信系统中，对于使用正交频分复用(ofdm)技术的通信系统而言，有效信号带宽内的窄带干扰会导致部分子载波上传输的有效信息遭受不同程度的损失。此外，在使用ofdm技术的接收机中，使用有限点数的快速傅里叶变换(fft)将信号从时域变换到频域的过程会导致窄带干扰泄漏到相邻的子载波上。因此，在接收机中需要完成窄带干扰检测并将频域被干扰的子载波屏蔽。
在实际应用中，由于电磁环境的不可预知性，窄带干扰通常表征为干扰的数目不同、干扰的频率不同、干扰的带宽和强度不相同，较高功率的窄带干扰可能会掩盖那些低功率的窄带干扰信号。为了尽可能多的窄带干扰被检出，通常的检测是基于迭代的思想完成的,其大致过程如下：首先计算所有子载波上的平均功率；然后使用门限值nbithr去屏蔽受污染的子载波，若第i个子载波上的功率pmin>pavg*nbithr，那么我们就认为[i–nbiwidth,i nbiwidth]区间内的子载波上检测到了窄带干扰，其中nbiwidth为预设的单边带窄带干扰屏蔽宽度；最后，基于更新屏蔽信息后的子载波计算新的平均功率，并仿照上一步进行新一轮的窄带干扰检测。
上述常规方法的弊端在于每一次检测的可靠性低，通过多次迭代检测的处理延时不固定且不可预测，会导致硬件的实时处理难以控制。
技术实现要素：
本发明的目的在于提供一种基于畸变功率去除的窄带干扰检测方法及其通信装置，以解决上述问题。
为实现上述目的，本发明采用以下技术方案：
一种基于畸变功率去除的窄带干扰检测方法，包括以下步骤：
步骤1，获取噪声功率最小的子载波pmin，以及噪声功率最小的子载波号ind_min；
步骤2，计算剔除畸变功率的平均噪声功率；
步骤3，第一次窄带干扰检测；
1)计算子载波的功率pj与平均功率集合中元素个数n平均的乘积tmppowunmask；
2)计算平均噪声功率集合中总功率v平均与去除窄带干扰的可编程阈值d_nbithd的乘积tmppowthd；
3)若tmppowunmask大于等于tmppowthd，该子载波存在干扰，将数组nbi_null_vec[j-nbiwidth，j nbiwidth]写“0”；数组nbi_null_vec[0：n-1]用来表征子载波是否存在干扰，其中nbi_null_vec[i]为“0”代表第i个子载波存在干扰，“1”代表第i个子载波不存在干扰；
4)遍历子载波pnbiwidth、pnbiwidth 1……pn-1-nbiwidth之后，完成第一次检测；
步骤4，第二次窄带干扰检测；
1)对于第j个子载波，若步骤3中nbi_null_vec[j]值为“1”，则nsum＝nsum 1，psum＝psum pj，遍历所有子载波p0、p1……pn-1，得到第一次检测无干扰子载波总数nsum，第一次检测无干扰子载波总功率psum；
2)对子载波pnbiwidth、pnbiwidth 1……pn-1-nbiwidth，分别进行以下步骤；
3)将pj作为检测子载波，计算子载波功率pj与1)中第一次检测无干扰子载波总数nsum的乘积com2；
4)计算1)中第一次检测无干扰子载波功率psum与预设值th_glb的乘积com3；
5)若com2大于com3，该子载波存在干扰，更新数组nbi_null_vec[j-nbiwidth，j nbiwidth]为“0”；
6)遍历完子载波pnbiwidth、pnbiwidth 1……pn-1-nbiwidth后，输出所有存在窄带干扰的子载波。
进一步的，步骤2具体包括：
步骤s21，建立计算平均噪声功率集合；并且将pmin作为计算平均功率集合的第一个元素；
步骤s22，对子载波pind_min 1、pind_min 2……pn-1，分别进行步骤s23-s25；
步骤s23，将pj作为检测子载波，得到子载波pj前m个子载波的功率之和vj-m，以及子载波pj之后m个子载波的功率之和vj m，取pj vj-m和pj vj m中的最大值；
步骤s24，计算平均功率集合中小于等于m 1个子载波之和与去除功率的可编程阈值d_powthd的乘积com1；
步骤s25，若com1大于步骤s23中的最大值，则将子载波pj加入计算平均功率集合中；
步骤s26，遍历完子载波pind_min 1、pind_min 2……pn-1之后，得到最终的计算平均功率集合，并且计算该集合所有子载波的总功率v平均，计算该集合元素的个数n平均。
进一步的，步骤1中最小功率值为从噪声功率的起点开始取若干个子载波，得到噪声功率值最小子载波。
进一步的，步骤2中计算平均功率集合中的元素数量小于m 1时，计算集合中所有元素之和与去除功率的可编程阈值d_powthd的乘积com1；步骤2中计算平均功率集合中的元素数量大于m 1时，计算集合中最后m 1个元素之和与去除功率的可编程阈值d_powthd的乘积com1。
进一步的，步骤2中去除功率的可编程阈值d_powthd的值为5-15。
进一步的，步骤3中去除窄带干扰的可编程阈值d_nbithd的值为6-10。
进一步的，所述步骤4中遍历完pnbiwidth、pnbiwidth 1……pn-1-nbiwidth中每一个子载波，得到所有存在窄带干扰的子载波。
进一步的，在该方法运行的过程中，检测到数据包传输，则停止该方法的运行。
进一步的，一种通信装置，该通信装置具有数字信号处理芯片dsp或者专用集成电路芯片asic。
与现有技术相比，本发明有以下技术效果：
本发明提出的方法，其窄带干扰的检测过程所需要的运算复杂度固定，算法收敛快，检测性能卓越，适于实现asic芯片或dsp芯片。该方法在数据包传输的空闲时间段完成干扰的检测，避免了由于窄带信道造成的窄带干扰虚警，并且该方法与数据包检测结合紧密，提高了系统检测性能；该方法中第一次检测基于剔除畸变功率点的平均噪声功率集完成检测，当实际信道中包含窄带干扰强度比较低的情况，通过第二次检测搜索，确保信道中所有窄带干扰被正确检出。通过两次检测，在保证系统性能最优的情况下，实现最低的运算复杂度。
更进一步的，该方法在计算平均功率集合时剔除了畸变的功率点，提高了平均功率与实际噪声功率的相似度。
更进一步的，该方法中计算平均功率集合中元素之和依据集合中元素的数量进行相应的设置，使计算结果更可靠。
本发明的有益效果还在于：该方法能够适用于不同的通信芯片，包括专用电路芯片和数字信号处理芯片；该通信芯片能够检测多路窄带干扰，延时低且检测性能优越；该芯片的核心参数能够灵活配置，通过两次检测，在确保系统检测性能最优的情况下，实现最低的运算复杂度。
附图说明
图1为本发明方法的流程示意图；
图2为本发明方法中计算平均噪声功率的流程示意图；
图3为本发明中具有高噪声功率子载波的功率示意图。
具体实施方式
以下结合附图对本发明进一步说明：
请参阅图1至图3，一种基于畸变功率去除的窄带干扰检测方法，包括以下步骤：
步骤1，获取噪声功率最小的子载波pmin，以及噪声功率最小的子载波号ind_min；
步骤2，计算剔除畸变功率的平均噪声功率；
步骤3，第一次窄带干扰检测；
1)计算子载波的功率pj与平均功率集合中元素个数n平均的乘积tmppowunmask；
2)计算平均噪声功率集合中总功率v平均与去除窄带干扰的可编程阈值d_nbithd的乘积tmppowthd；
3)若tmppowunmask大于等于tmppowthd，该子载波存在干扰，将数组nbi_null_vec[j-nbiwidth，j nbiwidth]写“0”；数组nbi_null_vec[0：n-1]用来表征子载波是否存在干扰，其中nbi_null_vec[i]为“0”代表第i个子载波存在干扰，“1”代表第i个子载波不存在干扰；
4)遍历子载波pnbiwidth、pnbiwidth 1……pn-1-nbiwidth之后，完成第一次检测；
步骤4，第二次窄带干扰检测；
1)对于第j个子载波，若步骤3中nbi_null_vec[j]值为“1”，则nsum＝nsum 1，psum＝psum pj，遍历所有子载波p0、p1……pn-1，得到第一次检测无干扰子载波总数nsum，第一次检测无干扰子载波总功率psum；
2)对子载波pnbiwidth、pnbiwidth 1……pn-1-nbiwidth，分别进行以下步骤；
3)将pj作为检测子载波，计算子载波功率pj与1)中第一次检测无干扰子载波总数nsum的乘积com2；
4)计算1)中第一次检测无干扰子载波功率psum与预设值th_glb的乘积com3；
5)若com2大于com3，该子载波存在干扰，更新数组nbi_null_vec[j-nbiwidth，j nbiwidth]为“0”；
6)遍历完子载波pnbiwidth、pnbiwidth 1……pn-1-nbiwidth后，输出所有存在窄带干扰的子载波。
步骤2具体包括：
步骤s21，建立计算平均噪声功率集合；并且将pmin作为计算平均功率集合的第一个元素；
步骤s22，对子载波pind_min 1、pind_min 2……pn-1，分别进行步骤s23-s25；
步骤s23，将pj作为检测子载波，得到子载波pj前m个子载波的功率之和vj-m，以及子载波pj之后m个子载波的功率之和vj m，取pj vj-m和pj vj m中的最大值；
步骤s24，计算平均功率集合中小于等于m 1个子载波之和与去除功率的可编程阈值d_powthd的乘积com1；
步骤s25，若com1大于步骤s23中的最大值，则将子载波pj加入计算平均功率集合中；
步骤s26，遍历完子载波pind_min 1、pind_min 2……pn-1之后，得到最终的计算平均功率集合，并且计算该集合所有子载波的总功率v平均，计算该集合元素的个数n平均。
步骤1中最小功率值为从噪声功率的起点开始取若干个子载波，得到噪声功率值最小子载波。
步骤2中计算平均功率集合中的元素数量小于m 1时，计算集合中所有元素之和与去除功率的可编程阈值d_powthd的乘积com1；步骤2中计算平均功率集合中的元素数量大于m 1时，计算集合中最后m 1个元素之和与去除功率的可编程阈值d_powthd的乘积com1。
步骤2中去除功率的可编程阈值d_powthd的值为5-15。
步骤3中去除窄带干扰的可编程阈值d_nbithd的值为6-10。
所述步骤4中遍历完pnbiwidth、pnbiwidth 1……pn-1-nbiwidth中每一个子载波，得到所有存在窄带干扰的子载波。
在该方法运行的过程中，检测到数据包传输，则停止该方法的运行。
一种通信装置，该通信装置具有数字信号处理芯片dsp或者专用集成电路芯片asic。
实施例为：
本实施例中检测窄带干扰的方法在传输数据包传输的空闲时间段进行，若检测到数据包传输，则停止该检测方法。
传输带宽上有128个子载波，分别为：p0、p1…pi…p126、p127；选取前10个子载波，并且得到噪声功率最小的子载波p4。
然后构建平均功率集合，并且将p4作为该集合的第一个子载波。
从子载波p5开始遍历所有的子载波，检测后续子载波是否能够加入平均功率集合，具体过程如下：
计算子载波p5之前的2个子载波的功率之和得到v5-2，计算子载波p5之后的2个子载波的功率之和得到v5 2；其中p5之前的2个子载波包括p3、p4，p5之后的2个子载波包括p6和p7。在本发明中如果检测到最后两个子载波p126、p127，后两个子载波之和分别为v126 2＝p126 p127、v127 2＝p127 p127。
比较p5 v5 2和p5 v5-2的大小，得到其中的最大值。
计算平均功率集合中小于等于3个子载波和与去除功率的可编程阈值的乘积com1。
比较com1与上述p5 v5 2和p5 v5-2的最大值，若com1大于该最大值，则将子载波p5加入到计算平均功率集合中，否则，不更新平均功率集合。
在检测子载波p15是否能够加入计算平均功率集合时，此时计算平均功率集合有4个元素，按照先后顺序分别为p4、p6、p9和p12。
计算子载波p15之前2个子载波功率之和得到v15-2，计算子载波p15之后2个子载波功率之和得到v15 2；其中p15之前的2个子载波包括p13、p14，p15之后的2个子载波包括p16、p17。
比较p15 v15 2和p15 v15-2的大小，得到其中的最大值。
计算平均功率集合中3个子载波和与去除功率的可编程阈值的乘积com1。在检测子载波p15时，计算平均功率集合中有4个子载波，选择后更新的三个子载波p6、p9和p12，则com1为三个子载波和与去除功率的可编程阈值的乘积。
比较com1与上述p15 v15 2和p15 v15-2的最大值，若com1大于该最大值，p15更新到计算平均功率集合。
按照顺序对所有子载波完成检测，得到的计算平均功率集合中共有65个子载波，计算这65个子载波的功率之和v平均。
计算平均噪声功率集合的总功率v平均与去除干扰的可编程阈值的乘积，得到tmppowthd；将每一个子载波的功率乘以65与tmppowthd进行比较，若大于tmppowthd，则该子载波存在窄带干扰，更新表征子载波是否存在干扰的数组nbi_null_vec[i-5，i 5]中元素为“0”，得到第一次检测后所有子载波是否存在干扰的标志信息。
统计第一次检测后无干扰子载波的个数nsum和总功率psum。
重新开始检测，以第5个子载波为例，计算p5与nsum的乘积com2，计算psum与预设值的乘积com3，比较com2与com3，若com2大于com3，更新nbi_null_vec[0，10]为“0”，依次对所有子载波进行处理，得到最终的窄带干扰子载波。
本发明还提供了一种通信装置，该通信装置具有数字信号处理芯片(dsp)或者专用集成电路芯片(asic)。
1.一种基于畸变功率去除的窄带干扰检测方法，其特征在于，包括以下步骤：
步骤1，获取噪声功率最小的子载波pmin，以及噪声功率最小的子载波号ind_min；
步骤2，计算剔除畸变功率的平均噪声功率；
步骤3，第一次窄带干扰检测；
1)计算子载波的功率pj与平均功率集合中元素个数n平均的乘积tmppowunmask；
2)计算平均噪声功率集合中总功率v平均与去除窄带干扰的可编程阈值d_nbithd的乘积tmppowthd；
3)若tmppowunmask大于等于tmppowthd，该子载波存在干扰，将数组nbi_null_vec[j-nbiwidth，j nbiwidth]写“0”；数组nbi_null_vec[0：n-1]用来表征子载波是否存在干扰，其中nbi_null_vec[i]为“0”代表第i个子载波存在干扰，“1”代表第i个子载波不存在干扰；
4)遍历子载波pnbiwidth、pnbiwidth 1……pn-1-nbiwidth之后，完成第一次检测；
步骤4，第二次窄带干扰检测；
1)对于第j个子载波，若步骤3中nbi_null_vec[j]值为“1”，则nsum＝nsum 1，psum＝psum pj，遍历所有子载波p0、p1……pn-1，得到第一次检测无干扰子载波总数nsum，第一次检测无干扰子载波总功率psum；
2)对子载波pnbiwidth、pnbiwidth 1……pn-1-nbiwidth，分别进行以下步骤；
3)将pj作为检测子载波，计算子载波功率pj与1)中第一次检测无干扰子载波总数nsum的乘积com2；
4)计算1)中第一次检测无干扰子载波功率psum与预设值th_glb的乘积com3；
5)若com2大于com3，该子载波存在干扰，更新数组nbi_null_vec[j-nbiwidth，j nbiwidth]为“0”；
6)遍历完子载波pnbiwidth、pnbiwidth 1……pn-1-nbiwidth后，输出所有存在窄带干扰的子载波。
2.根据权利要求1所述的一种基于畸变功率去除的窄带干扰检测方法，其特征在于，步骤2具体包括：
步骤s21，建立计算平均噪声功率集合；并且将pmin作为计算平均功率集合的第一个元素；
步骤s22，对子载波pind_min 1、pind_min 2……pn-1，分别进行步骤s23-s25；
步骤s23，将pj作为检测子载波，得到子载波pj前m个子载波的功率之和vj-m，以及子载波pj之后m个子载波的功率之和vj m，取pj vj-m和pj vj m中的最大值；
步骤s24，计算平均功率集合中小于等于m 1个子载波之和与去除功率的可编程阈值d_powthd的乘积com1；
步骤s25，若com1大于步骤s23中的最大值，则将子载波pj加入计算平均功率集合中；
步骤s26，遍历完子载波pind_min 1、pind_min 2……pn-1之后，得到最终的计算平均功率集合，并且计算该集合所有子载波的总功率v平均，计算该集合元素的个数n平均。
3.根据权利要求1所述的一种基于畸变功率去除的窄带干扰检测方法，其特征在于，步骤1中最小功率值为从噪声功率的起点开始取若干个子载波，得到噪声功率值最小子载波。
4.根据权利要求2所述的一种基于畸变功率去除的窄带干扰检测方法，其特征在于，步骤2中计算平均功率集合中的元素数量小于m 1时，计算集合中所有元素之和与去除功率的可编程阈值d_powthd的乘积com1；步骤2中计算平均功率集合中的元素数量大于m 1时，计算集合中最后m 1个元素之和与去除功率的可编程阈值d_powthd的乘积com1。
5.根据权利要求2所述的一种基于畸变功率去除的窄带干扰检测方法，其特征在于，步骤2中去除功率的可编程阈值d_powthd的值为5-15。
6.根据权利要求1所述的一种基于畸变功率去除的窄带干扰检测方法，其特征在于，步骤3中去除窄带干扰的可编程阈值d_nbithd的值为6-10。
7.根据权利要求1所述的一种基于畸变功率去除的窄带干扰检测方法，其特征在于，所述步骤4中遍历完pnbiwidth、pnbiwidth 1……pn-1-nbiwidth中每一个子载波，得到所有存在窄带干扰的子载波。
8.根据权利要求1所述的一种基于畸变功率去除的窄带干扰检测方法，其特征在于，在该方法运行的过程中，检测到数据包传输，则停止该方法的运行。
9.一种通信装置，其特征在于，该通信装置具有数字信号处理芯片dsp或者专用集成电路芯片asic。
技术总结