一种基于尺度迭代和频谱补偿的宽带信号检测方法

1.本发明属于无线通信技术领域，尤其涉及一种基于尺度迭代和频谱补偿的宽带信号检测方法。


背景技术：

2.宽带信号检测的目的是在宽带频谱上快速找到信号，并估计出各个信号中心频率和带宽，为后续信号分析与处理提供先验信息。在各类民用和军用通信系统中均有大量应用，例如卫星频谱检测、雷达信号扫描以及战场信息的感知等均需要对宽带频谱进行检测。
3.然而，在非合作通信中，通常要在先验信息缺失、电磁环境复杂以及存在不平坦有色噪声基底等条件下完成频谱检测，这就要求宽带信号检测算法具备在低信噪比下无需先验信息实现盲检测并具备抗不平坦噪底的功能。
4.针对宽带信号的检测算法通常有如下几种：一是信号能量检测法，如urkowitz(1967)提出未知确定性信号的能量检测方法、黎严(2017)提出改进的多天线信号能量检测方案等，这类方法具有抗频偏的优点，但易受噪声影响，在低信噪比下检测效果较差；二是循环平稳检测法，如a.tani(2016)提出的基于二阶循环平稳特性算法，优点是检测性能较好，但复杂度较高，检测时间长；三是匹配滤波器检测法，如z.zhang(2010)提出的提出了平行匹配滤波器与分段滤波器相结合的检测方法，这类方法在理论上是一种最优的信号检测方法，但是需要载频、带宽、调制方式等大量先验信息，在实际中很难实现，不适用于非协作通信。另外，近年来还提出了一些新型算法，信号检测性能较经典算法有了很大的提升。张洋(2016)提出的梯度双阈值算法对宽带功率谱分段计算梯度,再进行自适应双阈值检测。齐佩汉(2014)提出的功率谱对消法利用分段频带内部分谱线强度和与全部谱线强度和的比值作为检验统计量进行信号存在性的判断。


技术实现要素：

5.针对现有的宽带信号的检测算法易受噪声影响、在低信噪比下检测效果较差、复杂度较高，检测时间长、需要大量先验信息等问题，而提出了一种基于尺度迭代和频谱补偿的宽带信号检测方法。
6.为了实现上述目的，本发明采用了如下技术方案：
7.一种基于尺度迭代和频谱补偿的宽带信号检测方法，包括如下步骤：
8.步骤a:接收信号经过ad采样得到宽带采样信号，经过数字信道化将宽带信号划分为合适带宽的多段信号；
9.步骤b:取数段宽带信号，据各段信号功率谱密度计算最大功率谱；
10.步骤c:利用尺度迭代形态学滤波算法估计并去除有色噪声基底，并计算出白噪声功率门限；
11.步骤d:检测功率谱中幅度最高的信号，估计出其中心频率、带宽以及信号频带起止位置，并根据补偿算法将信号剔除并完成补偿每次检出一个信号后，根据两次信号检测
前后残余信号能量差值高于预设能量门限，则以补偿后的功率谱再次进行步骤d，否则信号检测结束。
12.优选的，利用形态学滤波带宽筛选的特性，对功率谱进行形态学滤波迭代，使噪声基底估计更准确，利用平滑迭代克服信号中强脉冲和白噪声对信号参数估计造成的误差，通过频谱补偿改善频带混叠信号的检测性能。
13.优选的，所述步骤a具体包括如下步骤：
14.a1：计算实信号信道化子信道的中心频率d为数据抽取倍数，整个频带被分为d个实信号的对称子带；
15.a2：抗混叠滤波器h
lp
(n)为低通fir滤波器,k为信道个数,d为抽取因子,且n与d具有整数倍关系,即k＝fd,则传统信道化结构中的第k个子信道的输出为：
[0016][0017]
a3：将上式改写为多相滤波结构，得到第k个子信道的输出表达式为：
[0018][0019]
a4：将中心频率ωk代入式中,此时第k个子信道的输出表达式为：
[0020][0021]
优选的，所述步骤b中，将接收数据进行m次分段，每段长度为l，段间间隔为k，提取数据可以矩阵形式表示：
[0022][0023]
对每段数据以周期图法估计出各段数据的功率谱pi(f)，即
[0024][0025]
其中，w(n)是所选窗函数,求出各分段信号功率谱后在每个频率点处求出功率谱的最大值，即
[0026][0027]
优选的，所述步骤c具体包括如下步骤：
[0028]
b1:尺度b
bi
取值最小检测带宽b
min
为进行一次形态学滤波pn1(f)；
[0029]
b2:尺度b
bi
取值上一次滤波尺度b
bi-1
的二倍进行一次形态学滤波，若尺度大于最大检测带宽b
max
，则以b
max
为尺度进行一次形态学滤波，第i次滤波结果为pni(f)；
[0030]
b3:将第i次和第i-1次相邻两次滤波结果差分，差分结果为δpni，逐频点比较两次差分结果，若δpni(f)＞δpn
i-1
(f)，则将频点f处尺度更新为b
bi
，否则尺度保持不变；
[0031]
b4:若b
bi
＝b
max
，则停止迭代；
[0032]
b5:各个频点f处以迭代确定的尺度b
bi
(f)进行一次形态学滤波估计出有色噪声基底；
[0033]
b6:以原功率谱减b5估计出噪声基底得到噪底平坦的功率谱；
[0034]
b7:将功率谱转换为对数功率谱，在功率谱上由0至最大值划定适当个数的功率带并统计落在各带内的频点数，功率带个数越多，白噪声门限估计越准确，相邻功率带频点数分别差分，以差分结果最大的功率带的上限作为白噪声功率门限。
[0035]
优选的，所述步骤d中具体包括如下步骤：
[0036]
c1:在功率谱中寻找最大值位置f
max
作为中心频率f
ctr
，在中心频率幅度p(f
ctr
)-3db处计算出3db带宽b
3db
，带宽左边界为f
3dbl
，右边界为f
3dbr
；
[0037]
c2:将功率谱在[f
3dbl
,f
3dbr
]以尺度λ
sm
＝0.1*b
3db
进行平滑滤波处理后返回c1，循环五次后进入c3；
[0038]
c3:经过五次循环平滑后，所得带宽即为3db带宽，在最终的[f
3dbl
,f
3dbr
]范围内根据下式计算中心频率；
[0039][0040]
c4:在f
3dbl
以左，f
3dbr
以右寻找功率谱首个极小值点或零点作为信号主瓣边界，左边界为f
ls
，右边界为f
rs
，在[f
ls
,f
rs
]范围内将信号归零剔除，根据下式计算带外斜率k
ls
、k
rs
：
[0041][0042][0043]
c5:在f
ls
，f
rs
处分别以斜率k
ls
、k
rs
向已估计信号带内[f
ls
,f
rs
]延伸，至两线相交为止，以两条延长线代替已估计信号的带内谱线，减小对为未估计信号的损伤；
[0044]
c6:若c4已估计能量低于最小判定能量或带宽小于最小检测带宽，则将其舍弃，若残余功率谱最大幅度不高于白噪声门限，则信号检测完成。
[0045]
本发明所具有的有益效果为：针对宽带信号，利用形态学滤波带宽筛选的特性，对功率谱进行形态学滤波迭代，使噪声基底估计更准确，避免不平坦噪底造成误检；利用平滑迭代克服信号中强脉冲和白噪声对信号参数估计造成的误差，通过频谱补偿改善频带混叠信号的检测性能，实现宽带信号的盲检测，提高检测准确度。
附图说明
[0046]
图1为一种基于尺度迭代和频谱补偿的宽带信号检测方法的步骤流程图；
[0047]
图2为数字信道化结构；
[0048]
图3为宽带信号频谱示意图；
[0049]
图4为信号检测环节流程；
[0050]
图5为有色噪声基底估计结果示意图；
[0051]
图6为信号检测结果示意图；
[0052]
图7为在不同信噪比下的信号正确概率。
具体实施方式
[0053]
下面将对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。
[0054]
实施例一
[0055]
如图1-7所示，种基于尺度迭代和频谱补偿的宽带信号检测方法，包括如下步骤：
[0056]
步骤a:接收信号经过ad采样得到宽带采样信号，经过数字信道化将宽带信号划分为合适带宽的多段信号；
[0057]
步骤b:取数段宽带信号，据各段信号功率谱密度计算最大功率谱；
[0058]
步骤c:利用尺度迭代形态学滤波算法估计并去除有色噪声基底，并计算出白噪声功率门限；
[0059]
步骤d:检测功率谱中幅度最高的信号，估计出其中心频率、带宽以及信号频带起止位置，并根据补偿算法将信号剔除并完成补偿每次检出一个信号后，根据两次信号检测前后残余信号能量差值高于预设能量门限，则以补偿后的功率谱再次进行步骤d，否则信号检测结束。
[0060]
本实施例中，所述步骤a具体包括如下步骤：
[0061]
a1：计算实信号信道化子信道的中心频率d为数据抽取倍数，整个频带被分为d个实信号的对称子带；
[0062]
a2：抗混叠滤波器h
lp
(n)为低通fir滤波器,k为信道个数,d为抽取因子,且n与d具有整数倍关系,即k＝fd,则传统信道化结构中的第k个子信道的输出为：
[0063][0064]
a3：将上式改写为多相滤波结构，得到第k个子信道的输出表达式为：
[0065][0066]
a4：将中心频率ωk代入式中,此时第k个子信道的输出表达式为：
[0067][0068]
本实施例中，所述步骤b中，将接收数据进行m次分段，每段长度为l，段间间隔为k，提取数据可以矩阵形式表示：
[0069]
[0070]
对每段数据以周期图法估计出各段数据的功率谱pi(f)，即
[0071][0072]
其中，w(n)是所选窗函数，求出各分段信号功率谱后在每个频率点处求出功率谱的最大值，即
[0073][0074]
本实施例中，所述步骤c具体包括如下步骤：
[0075]
b1:尺度b
bi
取值最小检测带宽b
min
为进行一次形态学滤波pn1(f)；
[0076]
b2:尺度b
bi
取值上一次滤波尺度b
bi-1
的二倍进行一次形态学滤波，若尺度大于最大检测带宽b
max
，则以b
max
为尺度进行一次形态学滤波，第i次滤波结果为pni(f)；
[0077]
b3:将第i次和第i-1次相邻两次滤波结果差分，差分结果为δpni，逐频点比较两次差分结果，若δpni(f)＞δpn
i-1
(f)，则将频点f处尺度更新为b
bi
，否则尺度保持不变；
[0078]
b4:若b
bi
＝b
max
，则停止迭代；
[0079]
b5:各个频点f处以迭代确定的尺度b
bi
(f)进行一次形态学滤波估计出有色噪声基底；
[0080]
b6:以原功率谱减b5估计出噪声基底得到噪底平坦的功率谱；
[0081]
b7:将功率谱转换为对数功率谱，在功率谱上由0至最大值划定适当个数的功率带并统计落在各带内的频点数，功率带个数越多，白噪声门限估计越准确，相邻功率带频点数分别差分，以差分结果最大的功率带的上限作为白噪声功率门限。
[0082]
本实施例中，所述步骤d中具体包括如下步骤：
[0083]
c1:在功率谱中寻找最大值位置f
max
作为中心频率f
ctr
，在中心频率幅度p(f
ctr
)-3db处计算出3db带宽b
3db
，带宽左边界为f
3dbl
，右边界为f
3dbr
；
[0084]
c2:将功率谱在[f
3dbl
,f
3dbr
]以尺度λ
sm
＝0.1*b
3db
进行平滑滤波处理后返回c1，循环五次后进入c3；
[0085]
c3:经过五次循环平滑后，所得带宽即为3db带宽，在最终的[f
3dbl
,f
3dbr
]范围内根据下式计算中心频率；
[0086][0087]
c4:在f
3dbl
以左，f
3dbr
以右寻找功率谱首个极小值点或零点作为信号主瓣边界，左边界为f
ls
，右边界为f
rs
，在[f
ls
,f
rs
]范围内将信号归零剔除，根据下式计算带外斜率k
ls
、k
rs
：
[0088][0089][0090]
c5:在f
ls
，f
rs
处分别以斜率k
ls
、k
rs
向已估计信号带内[f
ls
,f
rs
]延伸，至两线相交为
止，以两条延长线代替已估计信号的带内谱线，减小对为未估计信号的损伤；
[0091]
c6:若c4已估计能量低于最小判定能量或带宽小于最小检测带宽，则将其舍弃，若残余功率谱最大幅度不高于白噪声门限，则信号检测完成。
[0092]
实施例二
[0093]
如图1-7所示，一种基于尺度迭代和频谱补偿的宽带信号检测方法，包括如下步骤：
[0094]
步骤a:接收信号经过ad采样得到宽带采样信号，经过数字信道化将宽带信号划分为合适带宽的多段信号；
[0095]
步骤b:取数段宽带信号，据各段信号功率谱密度计算最大功率谱；
[0096]
步骤c:利用尺度迭代形态学滤波算法估计并去除有色噪声基底，并计算出白噪声功率门限；
[0097]
步骤d:检测功率谱中幅度最高的信号，并根据补偿算法将信号剔除并完成补偿每次检出一个信号后，根据两次信号检测前后残余信号能量差值高于预设能量门限，则以补偿后的功率谱再次进行步骤d，否则信号检测结束。
[0098]
本实施例中，所述步骤a具体包括如下步骤：
[0099]
a1：计算实信号信道化子信道的中心频率d为数据抽取倍数，整个频带被分为d个实信号的对称子带；
[0100]
a2：抗混叠滤波器h
lp
(n)为低通fir滤波器,k为信道个数,d为抽取因子,且n与d具有整数倍关系,即k＝fd,则传统信道化结构中的第k个子信道的输出为：
[0101][0102]
a3：将上式改写为多相滤波结构，得到第k个子信道的输出表达式为：
[0103][0104]
a4：将中心频率ωk代入式中,此时第k个子信道的输出表达式为：
[0105][0106]
本实施例中，所述步骤b中，将接收数据进行m次分段，每段长度为l，段间间隔为k，提取数据可以矩阵形式表示：
[0107][0108]
对每段数据以周期图法估计出各段数据的功率谱pi(f)，即
[0109]
[0110]
其中，w(n)是所选窗函数，求出各分段信号功率谱后在每个频率点处求出功率谱的最大值，即
[0111][0112]
本实施例中，所述步骤c具体包括如下步骤：
[0113]
b1:尺度b
bi
取值最小检测带宽b
min
为进行一次形态学滤波pn1(f)；
[0114]
b2:尺度b
bi
取值上一次滤波尺度b
bi-1
的二倍进行一次形态学滤波，第i次滤波结果为pni(f)；
[0115]
b3:将第i次和第i-1次相邻两次滤波结果差分，差分结果为δpni，若δpni(f)＞δpn
i-1
(f)，则将频点f处尺度更新为b
bi
，否则尺度保持不变；
[0116]
b4:若b
bi
＝b
max
，则停止迭代；
[0117]
b5:各个频点f处以迭代确定的尺度b
bi
(f)进行一次形态学滤波估计出有色噪声基底；
[0118]
b6:以原功率谱减b5估计出噪声基底得到噪底平坦的功率谱；
[0119]
b7:将功率谱转换为对数功率谱，在功率谱上由0至最大值划定适当个数的功率带并统计落在各带内的频点数，功率带个数越多，白噪声门限估计越准确，相邻功率带频点数分别差分，以差分结果最大的功率带的上限作为白噪声功率门限。
[0120]
本实施例中，所述步骤d中具体包括如下步骤：
[0121]
c1:在功率谱中寻找最大值位置f
max
作为中心频率f
ctr
，在中心频率幅度p(f
ctr
)-3db处计算出3db带宽b
3db
，带宽左边界为f
3dbl
，右边界为f
3dbr
；
[0122]
c2:将功率谱在[f
3dbl
,f
3dbr
]以尺度λ
sm
＝0.1*b
3db
进行平滑滤波处理后返回c1，循环五次后进入c3；
[0123]
c3:经过五次循环平滑后，所得带宽即为3db带宽，在最终的[f
3dbl
,f
3dbr
]范围内根据下式计算中心频率；
[0124][0125]
c4:在f
3dbl
以左，f
3dbr
以右寻找功率谱首个极小值点或零点作为信号主瓣边界，在[f
ls
,f
rs
]范围内将信号归零剔除，根据下式计算带外斜率k
ls
、k
rs
：
[0126][0127][0128]
c5:在f
ls
，f
rs
处分别以斜率k
ls
、k
rs
向已估计信号带内[f
ls
,f
rs
]延伸，至两线相交为止，减小对为未估计信号的损伤；
[0129]
c6:若c4已估计能量低于最小判定能量或带宽小于最小检测带宽，则将其舍弃，若残余功率谱最大幅度不高于白噪声门限，则信号检测完成。
[0130]
实施例三
[0131]
如图1-7所示，一种基于尺度迭代和频谱补偿的宽带信号检测方法，包括如下步
骤：
[0132]
步骤a:接收信号经过ad采样得到宽带采样信号，经过数字信道化将宽带信号划分为合适带宽的多段信号；
[0133]
步骤b:取数段宽带信号，据各段信号功率谱密度计算最大功率谱；
[0134]
步骤c:利用尺度迭代形态学滤波算法估计并去除有色噪声基底，并计算出白噪声功率门限；
[0135]
步骤d:检测功率谱中幅度最高的信号，根据两次信号检测前后残余信号能量差值高于预设能量门限，则以补偿后的功率谱再次进行步骤d，否则信号检测结束。
[0136]
本实施例中，所述步骤a具体包括如下步骤：
[0137]
a1：计算实信号信道化子信道的中心频率d为数据抽取倍数，整个频带被分为d个实信号的对称子带；
[0138]
a2：抗混叠滤波器h
lp
(n)为低通fir滤波器,k为信道个数,d为抽取因子,且n与d具有整数倍关系,即k＝fd,则传统信道化结构中的第k个子信道的输出为：
[0139][0140]
a4：将中心频率ωk代入式中,此时第k个子信道的输出表达式为：
[0141][0142]
本实施例中，所述步骤b中，将接收数据进行m次分段，每段长度为l，段间间隔为k，提取数据可以矩阵形式表示：
[0143][0144]
对每段数据以周期图法估计出各段数据的功率谱pi(f)，即
[0145][0146]
其中，w(n)是所选窗函数，求出各分段信号功率谱后在每个频率点处求出功率谱的最大值，即
[0147][0148]
本实施例中，所述步骤c具体包括如下步骤：
[0149]
b1:尺度b
bi
取值最小检测带宽b
min
为进行一次形态学滤波pn1(f)；
[0150]
b2:尺度b
bi
取值上一次滤波尺度b
bi-1
的二倍进行一次形态学滤波，若尺度大于最大检测带宽b
max
，则以b
max
为尺度进行一次形态学滤波，第i次滤波结果为pni(f)；
[0151]
b3:将第i次和第i-1次相邻两次滤波结果差分，差分结果为δpni，逐频点比较两次差分结果，若δpni(f)＞δpn
i-1
(f)，则将频点f处尺度更新为b
bi
，否则尺度保持不变；
[0152]
b4:各个频点f处以迭代确定的尺度b
bi
(f)进行一次形态学滤波估计出有色噪声基底；
[0153]
b5:以原功率谱减b5估计出噪声基底得到噪底平坦的功率谱；
[0154]
b6:将功率谱转换为对数功率谱，在功率谱上由0至最大值划定适当个数的功率带并统计落在各带内的频点数，以差分结果最大的功率带的上限作为白噪声功率门限。
[0155]
本实施例中，所述步骤d中具体包括如下步骤：
[0156]
c1:在功率谱中寻找最大值位置f
max
作为中心频率f
ctr
，在中心频率幅度p(f
ctr
)-3db处计算出3db带宽b
3db
，带宽左边界为f
3dbl
，右边界为f
3dbr
；
[0157]
c2:将功率谱在[f
3dbl
,f
3dbr
]以尺度λ
sm
＝0.1*b
3db
进行平滑滤波处理后返回c1，循环五次后进入c3；
[0158]
c3:所得带宽即为3db带宽，在最终的[f
3dbl
,f
3dbr
]范围内根据下式计算中心频率；
[0159][0160]
c4:在f
3dbl
以左，f
3dbr
以右寻找功率谱首个极小值点或零点作为信号主瓣边界，左边界为f
ls
，右边界为f
rs
，在[f
ls
,f
rs
]范围内将信号归零剔除，根据下式计算带外斜率k
ls
、k
rs
：
[0161][0162][0163]
c5:在f
ls
，f
rs
处分别以斜率k
ls
、k
rs
向已估计信号带内[f
ls
,f
rs
]延伸，至两线相交为止，以两条延长线代替已估计信号的带内谱线，减小对为未估计信号的损伤；
[0164]
c6:若c4已估计能量低于最小判定能量或带宽小于最小检测带宽，则将其舍弃，若残余功率谱最大幅度不高于白噪声门限，则信号检测完成。
[0165]
实施例四
[0166]
如图1-7所示，一种基于尺度迭代和频谱补偿的宽带信号检测方法，包括如下步骤：
[0167]
步骤a:接收信号经过ad采样得到宽带采样信号；
[0168]
步骤b:取数段宽带信号，据各段信号功率谱密度计算最大功率谱；
[0169]
步骤c:利用尺度迭代形态学滤波算法估计并去除有色噪声基底，并计算出白噪声功率门限；
[0170]
步骤d:检测功率谱中幅度最高的信号，估计出其中心频率、带宽以及信号频带起止位置，并根据补偿算法将信号剔除并完成补偿每次检出一个信号后，根据两次信号检测前后残余信号能量差值高于预设能量门限，则以补偿后的功率谱再次进行步骤d，否则信号检测结束。
[0171]
本实施例中，所述步骤a具体包括如下步骤：
[0172]
a1：计算实信号信道化子信道的中心频率d为数据抽取倍数，整个频带被分为d个实信号的对称子带；
[0173]
a2：抗混叠滤波器h
lp
(n)为低通fir滤波器,k为信道个数,d为抽取因子,且n与d具有整数倍关系,即k＝fd,则传统信道化结构中的第k个子信道的输出为：
[0174][0175]
a3：将上式改写为多相滤波结构，得到第k个子信道的输出表达式为：
[0176][0177]
a4：将中心频率ωk代入式中,此时第k个子信道的输出表达式为：
[0178][0179]
本实施例中，所述步骤b中，将接收数据进行m次分段，每段长度为l，段间间隔为k，提取数据可以矩阵形式表示：
[0180][0181]
对每段数据以周期图法估计出各段数据的功率谱pi(f)，即
[0182][0183]
其中，w(n)是所选窗函数，求出各分段信号功率谱后在每个频率点处求出功率谱的最大值，即
[0184][0185]
本实施例中，所述步骤c具体包括如下步骤：
[0186]
b1:尺度b
bi
取值最小检测带宽b
min
为进行一次形态学滤波pn1(f)；
[0187]
b2:尺度b
bi
取值上一次滤波尺度b
bi-1
的二倍进行一次形态学滤波，第i次滤波结果为pni(f)；
[0188]
b3:将第i次和第i-1次相邻两次滤波结果差分，若δpni(f)＞δpn
i-1
(f)，则将频点f处尺度更新为b
bi
，否则尺度保持不变；
[0189]
b4:若b
bi
＝b
max
，则停止迭代；
[0190]
b5:各个频点f处以迭代确定的尺度b
bi
(f)进行一次形态学滤波估计出有色噪声基底；
[0191]
b6:以原功率谱减b5估计出噪声基底得到噪底平坦的功率谱；
[0192]
b7:将功率谱转换为对数功率谱，功率带个数越多，白噪声门限估计越准确，相邻功率带频点数分别差分，以差分结果最大的功率带的上限作为白噪声功率门限。
[0193]
本实施例中，所述步骤d中具体包括如下步骤：
[0194]
c1:在功率谱中寻找最大值位置f
max
作为中心频率f
ctr
，在中心频率幅度p(f
ctr
)-3db处计算出3db带宽b
3db
，带宽左边界为f
3dbl
，右边界为f
3dbr
；
[0195]
c2:将功率谱在[f
3dbl
,f
3dbr
]以尺度λ
sm
＝0.1*b
3db
进行平滑滤波处理后返回c1，循环五次后进入c3；
[0196]
c3:经过五次循环平滑后，所得带宽即为3db带宽，在最终的[f
3dbl
,f
3dbr
]范围内根据下式计算中心频率；
[0197][0198]
c4:在f
3dbl
以左，f
3dbr
以右寻找功率谱首个极小值点或零点作为信号主瓣边界，左边界为f
ls
，右边界为f
rs
，在[f
ls
,f
rs
]范围内将信号归零剔除，根据下式计算带外斜率k
ls
、k
rs
：
[0199][0200][0201]
c5:在f
ls
，f
rs
处分别以斜率k
ls
、k
rs
向已估计信号带内[f
ls
,f
rs
]延伸，至两线相交为止，以两条延长线代替已估计信号的带内谱线，减小对为未估计信号的损伤；
[0202]
c6:若c4已估计能量低于最小判定能量或带宽小于最小检测带宽，则将其舍弃，若残余功率谱最大幅度不高于白噪声门限，则信号检测完成。
[0203]
实施例五
[0204]
如图1-7所示，一种基于尺度迭代和频谱补偿的宽带信号检测方法，包括如下步骤：
[0205]
步骤a:接收信号经过ad采样得到宽带采样信号，经过数字信道化将宽带信号划分为合适带宽的多段信号；
[0206]
步骤b:取数段宽带信号，据各段信号功率谱密度计算最大功率谱；
[0207]
步骤c:利用尺度迭代形态学滤波算法估计并去除有色噪声基底，并计算出白噪声功率门限；
[0208]
步骤d:检测功率谱中幅度最高的信号，估计出其中心频率、带宽以及信号频带起止位置，根据两次信号检测前后残余信号能量差值高于预设能量门限，则以补偿后的功率谱再次进行步骤d，否则信号检测结束。
[0209]
本实施例中，所述步骤a具体包括如下步骤：
[0210]
a1：计算实信号信道化子信道的中心频率整个频带被分为d个实信号的对称子带；
[0211]
a2：抗混叠滤波器h
lp
(n)为低通fir滤波器,n与d具有整数倍关系,即k＝fd,则传统信道化结构中的第k个子信道的输出为：
[0212][0213]
a3：将上式改写为多相滤波结构，得到第k个子信道的输出表达式为：
[0214][0215]
a4：将中心频率ωk代入式中,此时第k个子信道的输出表达式为：
[0216][0217]
本实施例中，所述步骤b中，将接收数据进行m次分段，每段长度为l，段间间隔为k，提取数据可以矩阵形式表示：
[0218][0219]
对每段数据以周期图法估计出各段数据的功率谱pi(f)，即
[0220][0221]
其中，w(n)是所选窗函数，求出各分段信号功率谱后在每个频率点处求出功率谱的最大值，即
[0222][0223]
本实施例中，所述步骤c具体包括如下步骤：
[0224]
b1:尺度b
bi
取值最小检测带宽b
min
为进行一次形态学滤波pn1(f)；
[0225]
b2:尺度b
bi
取值上一次滤波尺度b
bi-1
的二倍进行一次形态学滤波，以b
max
为尺度进行一次形态学滤波，第i次滤波结果为pni(f)；
[0226]
b3:将第i次和第i-1次相邻两次滤波结果差分，若δpni(f)＞δpn
i-1
(f)，则将频点f处尺度更新为b
bi
，否则尺度保持不变；
[0227]
b4:若b
bi
＝b
max
，则停止迭代；
[0228]
b5:各个频点f处以迭代确定的尺度b
bi
(f)进行一次形态学滤波估计出有色噪声基底；
[0229]
b6:以原功率谱减b5估计出噪声基底得到噪底平坦的功率谱；
[0230]
b7:在功率谱上由0至最大值划定适当个数的功率带并统计落在各带内的频点数，白噪声门限估计越准确，以差分结果最大的功率带的上限作为白噪声功率门限。
[0231]
本实施例中，所述步骤d中具体包括如下步骤：
[0232]
c1:在功率谱中寻找最大值位置f
max
作为中心频率f
ctr
，在中心频率幅度p(f
ctr
)-3db处计算出3db带宽b
3db
，带宽左边界为f
3dbl
，右边界为f
3dbr
；
[0233]
c2:将功率谱在[f
3dbl
,f
3dbr
]以尺度λ
sm
＝0.1*b
3db
进行平滑滤波处理后返回c1，循环五次后进入c3；
[0234]
c3:经过五次循环平滑后，所得带宽即为3db带宽，在最终的[f
3dbl
,f
3dbr
]范围内根据下式计算中心频率；
[0235][0236]
c4:在f
3dbl
以左，f
3dbr
以右寻找功率谱首个极小值点或零点作为信号主瓣边界，左边界为f
ls
，右边界为f
rs
，在[f
ls
,f
rs
]范围内将信号归零剔除，根据下式计算带外斜率k
ls
、k
rs
：
[0237][0238][0239]
c5:在f
ls
，f
rs
处分别以斜率k
ls
、k
rs
向已估计信号带内[f
ls
,f
rs
]延伸，至两线相交为止，减小对为未估计信号的损伤；
[0240]
c6:若c4已估计能量低于最小判定能量或带宽小于最小检测带宽，则将其舍弃，若残余功率谱最大幅度不高于白噪声门限，则信号检测完成。
[0241]
由此，信号检测环节流程如图4所示。
[0242]
算法性能仿真结果如图5-图7所示，仿真参数设定：信号采样率为8mhz，符号速率为2mbaud，调制方式为qpsk。采样时间设为25ms，共有3个tdma信号，每个信号突发时长为0.36ms。信道环境为加性高斯白噪声。
[0243]
为验证本专利的有色噪声基底估计准确度，图5给出了某宽带功率谱有色噪声基底估计结果及剔除噪底后效果。信号带宽37.5mhz，中心频率为963mhz，采样率为75mhz。可以看出，噪声基底同信号起伏贴合紧密，去除有色噪声后功率谱噪底平坦，仅残余白噪声，且信号个数无变化，说明该算法具有良好的噪底估计性能，适合在起伏噪底条件下使用。
[0244]
为验证本专利的有多信号检测性能，图6给出了某宽带功率谱信号检测结果。可以看出，强弱信号均可检出，且当信号重叠时，强信号并未对弱信号的检出造成影响。
[0245]
为进一步验证本文算法性能的稳定性，设置实验条件为：宽带信号内有25个子信号，调制方式包括bpsk、qpsk、8psk、8qam、16qam，信号带宽范围为2k～5mhz，高斯白噪声信道，信噪比为-14～6db。由于能量检测算法、梯度双阈值算法以及功率谱对消法均受不平坦噪底影响较大，功率谱噪底均为平坦噪底。如图7可以看出，本文算法正确检出概率始终远高于经典的能量检测算法，在信噪比极低时同梯度双阈值算法以及功率谱对消法相差不大，在-10db开始高于三种算法，在snr≥0db时检测概率高于90％。
[0246]
具体的基于尺度迭代和频谱补偿的宽带信号检测与测量方法系统计算过程可参见上述实施例，本发明实施例在此不再赘述。
[0247]
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换，但这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和
范围。