CN109861770B - 一种基于波束形成输出功率组合的宽带信号检测方法 - Google Patents
一种基于波束形成输出功率组合的宽带信号检测方法 Download PDFInfo
- Publication number
- CN109861770B CN109861770B CN201910204868.6A CN201910204868A CN109861770B CN 109861770 B CN109861770 B CN 109861770B CN 201910204868 A CN201910204868 A CN 201910204868A CN 109861770 B CN109861770 B CN 109861770B
- Authority
- CN
- China
- Prior art keywords
- signal
- array
- output power
- beam forming
- broadband
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 claims abstract description 22
- 230000010287 polarization Effects 0.000 claims abstract description 20
- 238000005070 sampling Methods 0.000 claims abstract description 13
- 238000012360 testing method Methods 0.000 claims abstract description 9
- 238000004364 calculation method Methods 0.000 abstract description 2
- 238000004422 calculation algorithm Methods 0.000 description 8
- 238000004088 simulation Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 238000000342 Monte Carlo simulation Methods 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000003657 Likelihood-ratio test Methods 0.000 description 1
- 230000001149 cognitive effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003071 parasitic effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明提供一种基于波束形成输出功率组合的宽带信号检测方法，具体过程为：对阵列接收到的远场宽带完全极化信号进行采样，将当前感兴趣的扫描角所对应的空域时延消去，得到仅包含载频对应的时延部分并将其作为波束形成的权矢量，根据所述权矢量估算波束形成输出功率，再将不同扫描角对应的波束形成输出功率进行组合构建信号检验统计量；将所述统计量与波束扫描次数的倒数作对比，根据对比结果判定是否存在信号。该方法计算复杂度低，且能够准确的实现对信号进行检测。
Description
技术领域
本发明属于阵列信号处理技术领域，尤其涉及一种基于波束形成输出功率组合的宽带信号检测方法。
背景技术
目前，阵列信号检测被广泛应用于许多领域，如雷达、声纳、电子监视和认知无线电等。一般地，对于信号的检测均是基于恒虚警概率检测(CFAR)或广义似然比检验(GLRT)方法，常用的方法是基于波束扫描或空域滤波方法，在每个空域间隔进行信号检测。除此之外，许多学者也对基于能量的信号检测方法进行研究。
传统算法对阵列信号检测多以能量或波束输出功率为基础，提出基于上述输出能量的检验统计量进行信号检测。文献“Blind Energy-based Detection for SpatialSpectrum Sensing”提出基于电子导向寄生阵列(ESPAR)对接收信号进行检测，ESPAR相当于利用硬件实现波束扫描，然后在每个波束扫描区域求取波束形成的输出功率，将所有区域的输出功率进行整理得最大功率和最小功率作(MMEN)为检验统计量。该方法由于仅使用了最大和最小功率，遗漏了其余功率信息，导致信号检测性能有所损失。
发明内容
为解决上述问题，本发明提供了一种基于波束形成输出功率组合的宽带信号检测方法，该方法能够准确对信号进行检测。
实现本发明的技术方案如下：
一种基于波束形成输出功率组合的宽带信号检测方法，具体过程为：
对阵列接收到的远场宽带完全极化信号进行采样，将当前感兴趣的扫描角所对应的空域时延消去，得到仅包含载频对应的时延部分并将其作为波束形成的权矢量，根据所述权矢量估算波束形成输出功率，再将不同扫描角对应的波束形成输出功率进行组合构建信号检验统计量；
将所述统计量与波束扫描次数的倒数作对比，根据对比结果判定是否存在信号。
进一步地，本发明估算出的波束形成输出功率
为：
其中，
表示指向角度Ψl的波束扫描权矢量，符号||·||表示
范数，T表示采样得到的数据总数，“H”表示共轭转置，x(t)表示各个快拍数据构造的观测矢量，所述快拍数据由阵列接收到的远场宽带完全极化信号按时间进行采样得到；
所述统计量表示为
当统计量与1/L之差小于设定的阈值时，判定不存在信号，否则，判定存在信号，由此完成了信号的检测，L为波束扫描次数。
进一步地，本发明
当所述阵列为标量阵列时，
当所述阵列为极化阵列时，
其中，θk,φk,γk,ηk分别为信号的方位角、俯仰角、极化辅助角和极化相位差；
sk(t)为入射的第k个远场宽带完全极化信号在第t时刻采样数据，
为第k个宽带信号、第m个阵元相对于参考阵元的传播时延，f0为载波频率。
有益效果：
第一，本发明是针对宽带信号进行检测，较之传统的宽带信号检测方法，NBPD方法基于导向延迟技术，将当前感兴趣的扫描角所对应的空域时延人为消去，得到仅包含载频对应的时延部分并将其作为波束形成的权矢量，然后利用柯西-施瓦茨不等式，将不同扫描角对应的波束形成输出功率进行组合以构建信号检验统计量，且理论阈值仅和扫描次数相关，该方法计算复杂度低，且能够准确的实现对信号进行检测。
第二，本发明通过仿真计算，证明本发明的宽带信号检测方法优于最大最小功率方法。
附图说明
图1为本发明采用的传感器(标量/极化)组成的阵列示意图，标量条件下仅考虑单个偶极子(实线表示)，xoy平面角度θ；极化接收模型考虑交叉偶极子(实线加虚线表示)，三维空间内角度(θ,φ)；
图2为本发明基于波束形成输出功率组合的宽带信号检测方法的流程图；
图3为本发明在标量阵列下宽带信号检测方法随信噪比变化的性能仿真图；
图4为本发明在标量阵列下宽带信号检测的接收特性曲线图；
图5为本发明在极化阵列下宽带信号检测方法随信噪比变化的性能仿真图；
图6为本发明在极化阵列下宽带信号检测的接收特性曲线图。
具体实施方式
下面结合附图并举实施例，对本发明进行详细描述。
波束形成输出功率近似表达推导：
假设接收阵列为M元，可以任意且无模糊的模型放置，考虑K个远场宽带完全极化信号入射到如图1所示的M元阵列，阵列对接收到的远场宽带完全极化信号按时间进行采样得到T个快拍数据，并根据各个快拍数据构造观测矢量x(t)，具体的：
其中，n(t)为零均值复高斯白噪声；sk(t)为入射的第k个远场宽带完全极化信号在第t时刻采样。
一般情况下，若采用标量阵列，sk(t)可表示如下
若采用极化阵列(交叉偶极子阵列)，sk(t)可表示如下
其中，
为第k个宽带信号(θk,φk,γk,ηk分别为信号的方位角、俯仰角、极化辅助角和极化相位差)、第m(m＝0,1...M-1)个阵元相对于参考阵元的传播时延，f0为载波频率。
令
为指向角度Ψl(Ψl为第l次扫描角度的统称，l＝1,…,L)的波束扫描权矢量，那么对当前的角度Ψl进行预导向延迟，(2)式可写为
其中，
为第l次扫描所对应的入射角度、第m(m＝0,1...M-1)个阵元相对于参考阵元的延迟。
当扫描角度刚好是信号的入射角时，即
此时sk(t)不为0，(4)式可写为
由此可得，此时标量接收模型下的波束扫描权矢量为
由于导向延迟与极化部分无关，因此极化接收模型下的波束扫描权矢量可直接写为
获取
之后，归一化波束形成输出功率可表示为
其中“H”表示共轭转置，“E”表示数学期望，
Rxx＝E{x(t)xH(t)} (8)
为阵列输出协方差矩阵。
在二元假设检验情况下，
表示不存在感兴趣信号；
表示存在感兴趣信号，它们对应的波束形成输出功率为，
其中σ2为噪声功率，且
Rvv＝E{v(t)vH(t)}＝σ2I (11)
I表示单位矩阵。
实际上，Rxx通过估计得到
因此，上述波束形成输出功率近似可表示为，
“||·||”表示
范数。
存在目标信号与不存在目标信号的情况推导：
对于任意的扫描角度Ψl，均有
成立。
根据柯西-施瓦茨不等式可得
其中，βl为一实数，不等式(15)中等号成立的充要条件为，
κ为定常数。
令β1＝β2＝…＝βL＝1可得，
通过式(5)可知，当存在感兴趣信号时，条件(17)一般不会成立，因此
相反，不存在感兴趣信号时，由于
成立，因此
基于上述两个理论推导所得到的结论，本发明提供了一种基于波束形成输出功率组合的宽带信号检测方法，具体过程为：
估算波束形成输出功率，并根据所述输出功率计算信号检验的统计量；
将所述统计量与1/L做对比，根据对比结果判定是否存在信号。
本发明实施例估算出的波束形成输出功率
为：
其中，
表示指向角度Ψl的波束扫描权矢量，符号||·||表示
范数，T表示采样得到的数据总数，“H”表示共轭转置，x(t)表示各个快拍数据构造的观测矢量，所述快拍数据由阵列接收到的远场宽带完全极化信号按时间进行采样得到；
所述统计量表示为
由公式(19)可知，检测阈值可近似选择为1/L，即当统计量与1/L之差小于设定的阈值时，判定不存在信号，否则，判定存在信号，由此完成了信号的检测。
当所述阵列为标量阵列时，
当所述阵列为极化阵列时，
其中，θk,φk,γk,ηk分别为信号的方位角、俯仰角、极化辅助角和极化相位差；
sk(t)为入射的第k个远场宽带完全极化信号在第t时刻采样，
为第k个宽带信号、第m个阵元相对于参考阵元的传播时延，f0为载波频率。
本发明的效果可以通过以下的仿真进一步说明：
标量阵列接收模型考虑由8个单偶极子组成均匀线阵，设由两个远场宽带BPSK信号入射到此阵列，信号来向分别为20°和40°，相位分别为π/6和π/3，中心频率为2e9，采样频率为4.4e9，下面的仿真中蒙特卡洛实验重复次数为1000，虚警概率为0.01，谱峰搜索的步长选为5°。
图3展示了本发明在标量阵列接收模型下，对宽带信号的检测概率随信噪比变化曲线，由图可知，本发明所提的NBPD算法与MMEN算法相比，其检测概率在同一信噪比条件下优于MMEN方法，且随着信噪比的升高NBPD方法始终优于MMEN方法；图4展示了本发明在标量阵列接收模型下，对宽带信号检测的接收特性曲线，此时信噪比为-18dB，由图可知，本发明所提的NBPD算法与MMEN算法相比，其检测概率在同一虚警概率条件下优于MMEN方法。
多极化阵列接收模型考虑由4个三极子组成均匀线阵，设由两个远场宽带BPSK信号入射到此阵列，信号来向分别为20°和40°，相位分别为π/6和π/3，中心频率为2e9，采样频率为4.4e9，下面的仿真中蒙特卡洛实验重复次数为1000，虚警概率为0.01，谱峰搜索的步长选为5°。
图5展示了本发明在极化阵列接收模型下，对宽带信号的检测概率随信噪比变化曲线，由图可知，本发明所提的NBPD算法与MMEN算法相比，其检测概率在同一信噪比条件下优于MMEN方法，且随着信噪比的升高NBPD方法始终优于MMEN方法；图6展示了本发明在极化阵列接收模型下，对宽带信号检测的接收特性曲线，此时信噪比为-16dB，由图可知，本发明所提的NBPD算法所提的MMEN算法相比，其检测概率在同一虚警概率条件下优于MMEN方法。
在本方法中，标量阵列或多极化阵列均可用来接收信号，基于上述阵列对待检测区域进行波束扫描，组合各个扫描角度的波束形成功率，获取信号检测的检验统计量。在此基础上，发明了一种基于波束形成输出功率组合的宽带信号检测方法(NBPD)。本方法相比于MMEN方法，具有相同的计算复杂度，且利用了所有波束扫描的输出功率，针对宽带信号在标量阵列和多极化阵列模型下进行检测。同时仿真表面，本发明所提算法对宽带信号检测性能优于现有方法。
当然，本发明还可有其他多种实施例，在不背离本发明精神及其实质的情况下，熟悉本领域的技术人员当可根据本发明作出各种相应的改变和变形，但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (3)
1.一种基于波束形成输出功率组合的宽带信号检测方法，其特征在于，具体过程为：
对阵列接收到的远场宽带完全极化信号进行采样，将当前感兴趣的扫描角所对应的空域时延消去，得到仅包含载频对应的时延部分并将其作为波束形成的权矢量，根据所述权矢量估算波束形成输出功率，再将不同扫描角对应的波束形成输出功率进行组合构建信号检验统计量；
将所述统计量与波束扫描次数的倒数作对比，根据对比结果判定是否存在信号。
2.根据权利要求1所述基于波束形成输出功率组合的宽带信号检测方法，其特征在于，估算出的波束形成输出功率
为：
其中，
表示指向角度Ψl的波束扫描权矢量，符号||·||表示l2范数，T表示采样得到的数据总数，“H”表示共轭转置，x(t)表示各个快拍数据构造的观测矢量，所述快拍数据由阵列接收到的远场宽带完全极化信号按时间进行采样得到；
所述统计量表示为
当统计量与1/L之差小于设定的阈值时，判定不存在信号，否则，判定存在信号，由此完成了信号的检测，L为波束扫描次数。
3.根据权利要求1所述基于波束形成输出功率组合的宽带信号检测方法，其特征在于，当所述阵列为标量阵列时，
当所述阵列为极化阵列时，
其中，θk,φk,γk,ηk分别为信号的方位角、俯仰角、极化辅助角和极化相位差；
其中，sk(t)为入射的第k个远场宽带完全极化信号在第t时刻采样数据，
为第k个宽带信号、第m个阵元相对于参考阵元的传播时延，f0为载波频率。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910204868.6A CN109861770B (zh) | 2019-03-18 | 2019-03-18 | 一种基于波束形成输出功率组合的宽带信号检测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910204868.6A CN109861770B (zh) | 2019-03-18 | 2019-03-18 | 一种基于波束形成输出功率组合的宽带信号检测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109861770A CN109861770A (zh) | 2019-06-07 |
| CN109861770B true CN109861770B (zh) | 2020-07-24 |
Family
ID=66901121
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910204868.6A Expired - Fee Related CN109861770B (zh) | 2019-03-18 | 2019-03-18 | 一种基于波束形成输出功率组合的宽带信号检测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109861770B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112994766B (zh) * | 2019-12-16 | 2022-12-27 | 华为技术有限公司 | 一种波束传输方法及相关设备 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1674462A (zh) * | 2004-03-25 | 2005-09-28 | 电子科技大学 | 利用干涉原理估计相干信号doa的方法 |
| CN103616661B (zh) * | 2013-10-31 | 2016-01-20 | 电子科技大学 | 一种稳健的远场窄带信号源个数估计方法 |
| US9979457B2 (en) * | 2015-05-08 | 2018-05-22 | Lg Electronics Inc. | Multi-beamforming training |
| CN106100775B (zh) * | 2016-08-23 | 2018-05-04 | 桂林电子科技大学 | 基于邻接矩阵的ofdm频谱感知方法 |
| CN108445486A (zh) * | 2018-03-13 | 2018-08-24 | 南京理工大学 | 基于协方差矩阵重建和导向矢量修正的波束形成方法 |
-
2019
- 2019-03-18 CN CN201910204868.6A patent/CN109861770B/zh not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN109861770A (zh) | 2019-06-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ye et al. | On the resiliency of MUSIC direction finding against antenna sensor coupling | |
| Liao et al. | New approaches to direction-of-arrival estimation with sensor arrays in unknown nonuniform noise | |
| US8669901B2 (en) | Method for determining azimuth and elevation angles of arrival of coherent sources | |
| CN109407055A (zh) | 基于多径利用的波束形成方法 | |
| CN103616661B (zh) | 一种稳健的远场窄带信号源个数估计方法 | |
| Xie et al. | A recursive angle-Doppler channel selection method for reduced-dimension space-time adaptive processing | |
| CN106291449A (zh) | 对称稳定分布噪声下波达方向角估计新方法 | |
| Qi et al. | Time-frequency DOA estimation of chirp signals based on multi-subarray | |
| Chen et al. | Enhanced target detection for HFSWR by 2-D MUSIC based on sparse recovery | |
| CN110196417B (zh) | 基于发射能量集中的双基地mimo雷达角度估计方法 | |
| CN101252382A (zh) | 一种宽频段信号极化与doa估计方法及装置 | |
| CN109861770B (zh) | 一种基于波束形成输出功率组合的宽带信号检测方法 | |
| CN103901421B (zh) | 基于对角减载的水声阵列smi-mvdr空间谱估计方法 | |
| Ping et al. | Study of 2D DOA estimation for uniform circular array in wireless location system | |
| He et al. | DOA estimation of wideband signals based on iterative spectral reconstruction | |
| CN114325565B (zh) | 一种基于子空间关系的阵列超分辨测向方法 | |
| CN115980721A (zh) | 一种无误差协方差矩阵分离的阵列自校正方法 | |
| Hameed et al. | Low complexity single snapshot DoA method | |
| CN108692718A (zh) | 基于盲波束的稳健导航抗干扰方法及其系统 | |
| CN104375116B (zh) | 一种基于无线传感器阵列的到达方向检测方法 | |
| Pesavento et al. | A theoretical and experimental performance study of a root-MUSIC algorithm based on a real-valued eigendecomposition | |
| Feng et al. | Wideband direction of arrival estimation based on the principal angle between subspaces | |
| Liu et al. | Compressive sensing for very high frequency radar with application to low-angle target tracking under multipath interference | |
| CN114609581B (zh) | 一种比较扫描方向矢量投影极大值的阵列超分辨测向方法 | |
| Tayem et al. | Capon root-MUSIC-like direction of arrival estimation based on real data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200724 |