CN107241131A - 一种利用信号非圆特性的波束形成方法 - Google Patents

一种利用信号非圆特性的波束形成方法 Download PDF

Info

Publication number
CN107241131A
CN107241131A CN201710408514.4A CN201710408514A CN107241131A CN 107241131 A CN107241131 A CN 107241131A CN 201710408514 A CN201710408514 A CN 201710408514A CN 107241131 A CN107241131 A CN 107241131A
Authority
CN
China
Prior art keywords
signal
rounded
vector
azimuth
phase
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.)
Granted
Application number
CN201710408514.4A
Other languages
English (en)
Other versions
CN107241131B (zh
Inventor
万群
刘江波
董芳圆
申艳妮
常煜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huizhongxingzhi Technology Beijing Co ltd
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201710408514.4A priority Critical patent/CN107241131B/zh
Publication of CN107241131A publication Critical patent/CN107241131A/zh
Application granted granted Critical
Publication of CN107241131B publication Critical patent/CN107241131B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/08—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
    • H04B7/0837—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using pre-detection combining
    • H04B7/0842—Weighted combining
    • H04B7/086—Weighted combining using weights depending on external parameters, e.g. direction of arrival [DOA], predetermined weights or beamforming

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本发明属于电子信息技术领域中的波束形成方法，特别是一种利用信号非圆特性的波束形成方法。本发明方法首先对接收的非圆信号进行时域采样得到观测矢量；然后，对观测矢量取实部和虚部获得了一个扩展的观测矢量的信号模型；再将搜索空域范围和非圆相位按角度划分网格获得权向量矩阵，并将其与扩展的观测矢量做内积，画出接收信号的三维方向图，该图最大值点对应的方位角和非圆相位即为非圆信号的波达方向和非圆相位的估计值；最后把来波方向估计值与信号的非圆系数结合，构造目标信号加权向量函数，然后求得期望信号方向图。基于信号非圆特性波束形成方法利用非圆特性做选择，滤除干扰信号，从而显著减少干扰对感兴趣信号DOA估计的影响。

Description

一种利用信号非圆特性的波束形成方法
技术领域
本发明属于电子信息技术领域中的波束形成方法，特别是一种利用信号非圆特性的波束形成方法。
背景技术
波束形成利用空间多传感阵列所构成的系统对空间信号进行发射或接收。通过对阵列接收信息的加权处理，使得阵列输出中期望信号的成分尽可能被有效的接收和增强，同时干扰和噪声得到有效抑制。在基于这种原理的现有方法中，典型的有最小方差无畸变(MVDR，Minimum Variance Distortionless Response)波束形成、线性多约束最小方差(LCMV，Linear Constraint Minimum Variance)波束形成算法等。
传统的线性波束形成器一般是基于期望信号、干扰和环境噪声都是平稳的随机过程的简单假设，其概率分布满足复圆对称高斯分布。当波束形成技术应用到一些需要人工调制信号(非圆信号)的领域时，如无线通信领域，传统的基于各个准则的线性波束形成器不再是最优的，因为这些人工调制信号是非平稳的随机过程，其统计特性与传统的线性波束形成器中关于信号的假设不再相同，仅利用传统的最优波束形成器对非平稳过程进行波束形成并不能完全利用观测矢量的全部信息。
对于均匀线阵而言，增加阵元数，能够提高常规波束形成方法的分辨率，但是不能降低旁瓣水平。在不增加阵元数的情况下，若进行非均匀加权处理，虽然可以提高常规波束形成方法的分辨率，但是波束的旁瓣水平会增加，同时会伴随信噪比增益损失。
发明内容
为解决现有技术的不足，本发明提出一种高分辨率的非圆信号波束形成方法。利用空域滤波器对接收的直达波进行时域采样，将信号的扩展序列与非圆相位和空间角度的联合矢量做内积，估计出信号来波方向和信号的非圆相位值，再由联合估计值更新空域滤波器的权向量，进而改变阵列方向图的指向和形状，使阵列的空域滤波器工作在最佳状态，最终滤出期望信号。
本发明的技术方案为：
首先，对接收的非圆信号进行时域采样得到观测矢量；然后，对观测矢量取实部和虚部获得了一个扩展的观测矢量的信号模型；再将搜索空域范围和非圆相位按角度划分网格获得权向量矩阵，并将其与扩展的观测矢量做内积，画出接收信号的三维方向图，该图最大值点对应的方位角和非圆相位即为非圆信号的波达方向和非圆相位的估计值；最后把来波方向估计值与信号的非圆系数结合，构造目标信号加权向量函数，然后求得期望信号方向图。
一种利用信号非圆特性的波束形成方法，具体步骤如下：
S1、接收机分别在每个采样时刻对所接收的信号进行采样，第n个采样时刻的采样序列表示为x(n)＝As(n)+n(n)＝[x1(n) x2(n) … xM(n)]T，其中，n为不为零的自然数，M为阵元数；
S2、对于每个采样时刻，取信号的实部Re(x(n))和虚部Im(x(n))，构成采样向量
S3、将扫描的空域按方位角划分成L个网格点，每个网格点对应的角度方向为θ1,θ2,...,θL，
同时，将非圆相位角划分成K个网格点，每个网格点对应的非圆相位角为φ1,φ2,...,φK，
对于方位角为θi，相位角为φj的导向矢量可表示为
同时，分别取ai,j的实部和虚部构成扩展向量
其中，
S4、计算并画出ci,j的二维幅度谱，分别求出幅度谱极大值点对应的方位角和相位角，即信号的来波方向和非圆相位：
S5、由非圆信号的非圆相位估计值和来波方向为的信号进行空域滤波，则加权向量为同时对加权向量也按实部和虚部进行扩展
S6、确定接收信号的方位角，即，将S2所述向量与S5所述加权向量为做内积，并画出相应的波形图。
本发明的有益效果是：
基于信号非圆特性波束形成方法利用非圆特性做选择，滤除干扰信号，从而显著减少干扰对感兴趣信号DOA估计的影响。并且，对具有不同非圆特性的多个信号同时入射的情况，可通过多个非圆相位的相关运算(均匀加权处理)，实现对入射信号的DOA分组估计，从而估计出感兴趣信号的DOA。
对于均匀线阵而言，增加阵元数，能够提高常规波束形成方法的分辨率，但是不能降低旁瓣水平。在不增加阵元数的情况下，若进行非均匀加权处理，虽然可以提高常规波束形成方法的分辨率，但是波束的旁瓣水平会增加，同时会伴随信噪比增益损失。本专利在不增加阵元数的情况下,通过利用信号的非圆特性进行均匀加权处理，不仅能提高分辨率，还能降低波束的旁瓣水平，并且不会引起信噪比增益损失。
具体实施方式
下面将结合实施例，对本发明方法进行进一步说明。
本实施方式以8个阵元的均匀线性天线阵列，窄带信号源场景为例。阵元序号分别为1、2、3、4、5、6、7、8，并默认为1阵元为参考阵元，各阵元相对于1阵元的距离分别为0.5米，1.0米，1.5米，2.0米，2.5米，3米，3.5米。设定空间中有两个互不相干的信号源，其入射方向分别为θ1＝7°,θ2＝12°的信是非圆相位分别为φs1＝72°,φs2＝240°的BPSK信号。理想空间白噪声为复高斯白噪声，信噪比为10dB。采样长度为256。
本发明的具体实施方式流程如下：
步骤1.接收机分别在每个采样时刻对所接收的信号进行采样，在n＝1时刻的采样序列可表示为
-0.2849 -1.5153i
-0.4061 -1.9216i
0.5365 -1.2301i
0.5782 -0.2108i
0.6328 -0.0278i
0.7663 -0.2512i
0.1910 +0.1125i
0.3622 -0.0236i
在n＝2,3,...,256时刻对所有接收的信号进行时域采样，得到信号时域序列x(n)；
步骤2.对于每个采样时刻，取信号的实部Re(x(n))和虚部Im(x(n))
构成采样向量如下给出第1，2列和最后两列前8行的值值；
-0.2849 -0.3418 … -0.7884 0.0442
-0.4061 0.2475 … 0.5831 -0.0895
0.5365 -0.5379 … -0.9193 -0.6744
0.5782 -0.0272 … -0.2618 -0.5453
0.6328 0.0718 … -0.1683 0.0480
0.7663 -0.2217 … 0.3917 0.1526
0.1910 0.5060 … -0.3524 0.1019
0.3622 -0.1618 … 0.3445 0.3186
步骤3.将扫描的空域按方位角划分成1801个网格点，每个网格点对应的角度方向为θ1,θ2,...,θ1801；同时，将非圆相位角划分成360个网格点，每个网格点对应的非圆相位角为φ1,φ2,...,φ360，对于方位角为θi，相位角为φj的导向矢量可表示为 其中，其中同时，分别取ai,j的实部和虚部构成扩展向量如下给出第一列a0,0的值
0
-3.1416
-6.2832
-9.4248
-12.5664
-15.7080
-18.8496
-21.9911
步骤4.计算并画出ci,j的二维幅度谱，分求出幅度谱两个极大值点分别对应的方位角和相位角也即信号的来波方向和非圆相位：
步骤5.由非圆信号的非圆相位估计值和来波方向为的信号进行空域滤波，则加权向量为同时对加权向量也按实部和虚部进行扩展下面为的值
步骤6.确定接收信号的方位角：将步骤2所述向量与步骤5所述加权向量为做内积，并画出相应的波形图。
仿真实验表明常规均匀波束形成法只有一个主瓣，不能分辨两个入射方向靠的很近的信号，而基于非圆特性的波束形成算法有两个主瓣，其估计值与理论值一致，能清楚的分辨两个信号，且该算法旁瓣水平比于常规均匀波束形成算法低。综上所述，相比于常规均匀波束形成算法，信号非圆特性波束形成算法通过利用信号的非圆特性进行均匀加权处理，不仅能提高分辨率，还能降低波束的旁瓣水平，并且不会引起信噪比增益损失。

Claims (1)

1.一种利用信号非圆特性的波束形成方法，其特征在于，具体步骤如下：
S1、接收机分别在每个采样时刻对所接收的信号进行采样，第n个采样时刻的采样序列表示为x(n)＝As(n)+n(n)＝[x1(n)x2(n)…xM(n)]T，其中，n为不为零的自然数，M为阵元数；
S2、对于每个采样时刻，取信号的实部Re(x(n))和虚部Im(x(n))，构成采样向量
S3、将扫描的空域按方位角划分成L个网格点，每个网格点对应的角度方向为θ1,θ2,...,θL，
同时，将非圆相位角划分成K个网格点，每个网格点对应的非圆相位角为φ1,φ2,...,φK，
对于方位角为θi，相位角为φj的导向矢量可表示为
同时，分别取ai,j的实部和虚部构成扩展向量
其中，
S4、计算并画出ci,j的二维幅度谱，分别求出幅度谱极大值点对应的方位角和相位角，即信号的来波方向和非圆相位：
S5、由非圆信号的非圆相位估计值和来波方向为的信号进行空域滤波，则加权向量为同时对加权向量也按实部和虚部进行扩展
S6、确定接收信号的方位角，即，将S2所述向量与S5所述加权向量为做内积，并画出相应的波形图。
CN201710408514.4A 2017-06-02 2017-06-02 一种利用信号非圆特性的波束形成方法 Active CN107241131B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710408514.4A CN107241131B (zh) 2017-06-02 2017-06-02 一种利用信号非圆特性的波束形成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710408514.4A CN107241131B (zh) 2017-06-02 2017-06-02 一种利用信号非圆特性的波束形成方法

Publications (2)

Publication Number Publication Date
CN107241131A true CN107241131A (zh) 2017-10-10
CN107241131B CN107241131B (zh) 2020-11-13

Family

ID=59985321

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710408514.4A Active CN107241131B (zh) 2017-06-02 2017-06-02 一种利用信号非圆特性的波束形成方法

Country Status (1)

Country Link
CN (1) CN107241131B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112671436A (zh) * 2019-10-15 2021-04-16 三星电子株式会社 通信装置及其数据接收方法
CN113221337A (zh) * 2021-04-23 2021-08-06 电子科技大学 一种基于模式搜索的近场稀疏阵设计方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101483280A (zh) * 2009-02-23 2009-07-15 重庆大学 一种稳健波束合成器权值求解方法
CN105302936A (zh) * 2015-08-31 2016-02-03 中国科学院声学研究所 基于相关计算和协方差矩阵重构的自适应波束形成方法
CN106788653A (zh) * 2016-11-25 2017-05-31 西安电子科技大学 一种基于协方差矩阵重构的自适应波束形成方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101483280A (zh) * 2009-02-23 2009-07-15 重庆大学 一种稳健波束合成器权值求解方法
CN105302936A (zh) * 2015-08-31 2016-02-03 中国科学院声学研究所 基于相关计算和协方差矩阵重构的自适应波束形成方法
CN106788653A (zh) * 2016-11-25 2017-05-31 西安电子科技大学 一种基于协方差矩阵重构的自适应波束形成方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JIANGBO LIU 等: "Robust adaptive beamforming for noncircular signal against array steering vector mismatch and interference nonstationarity", 《2015 IEEE CHINA SUMMIT AND INTERNATIONAL CONFERENCE ON SIGNAL AND INFORMATION PROCESSING (CHINASIP)》 *
孙心宇 等: "一种基于非圆相干信号的波束成形算法", 《宇航学报》 *
石运梅: "基于非圆信号的自适应波束形成算法研究", 《万方学位论文》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112671436A (zh) * 2019-10-15 2021-04-16 三星电子株式会社 通信装置及其数据接收方法
CN113221337A (zh) * 2021-04-23 2021-08-06 电子科技大学 一种基于模式搜索的近场稀疏阵设计方法
CN113221337B (zh) * 2021-04-23 2022-09-09 电子科技大学 一种基于模式搜索的近场稀疏阵设计方法

Also Published As

Publication number Publication date
CN107241131B (zh) 2020-11-13

Similar Documents

Publication Publication Date Title
US10571544B2 (en) Direction finding using signal power
CN109407055A (zh) 基于多径利用的波束形成方法
Liu et al. Low angle estimation in MIMO radar
CN107643532B (zh) 卫星导航信号处理方法及装置、计算机设备及存储介质
CN101644760B (zh) 一种适用于高分辨阵列的快速鲁棒的信源个数检测方法
WO2019200663A1 (zh) 一种基于频率分集阵列天线的波束控制方法
CN106646529A (zh) 一种基于多波束优选的gnss天线阵抗干扰方法
CN102664666A (zh) 一种高效的宽带稳健自适应波束形成方法
CN108880586B (zh) 一种宽带弱信号增强方法与装置
KR101807499B1 (ko) 상호간섭효과를 보상한 원형 배열 안테나 기반의 도래각 추정 장치 및 추정 방법
CN110708103A (zh) 一种无需预延迟的宽带波束形成方法
CN103293517A (zh) 基于脊参数估计的对角加载稳健自适应雷达波束形成方法
CN107241131A (zh) 一种利用信号非圆特性的波束形成方法
CN111880198A (zh) 基于交替极化敏感阵列的空时极化抗干扰方法
CN112834981B (zh) 一种脉冲噪声背景下的零陷阵波达方位估计方法
CN104346532B (zh) 一种mimo雷达降维自适应波束形成方法
CN106338742B (zh) 基于互谱准则的降维自适应多波束gps信号抗干扰方法
Cho et al. Efficient implementation of the Capon beamforming using the Levenberg-Marquardt scheme for two dimensional AOA estimation
CN108416105A (zh) 脉冲和高斯噪声下稳健的自适应波束形成算法
JP2004336390A (ja) アダプティブアレーおよび測位装置
Zhou et al. Research on interferometer direction finding technology based on digital beam forming
RU2431862C1 (ru) Способ поляризационно-независимого пеленгования многолучевых радиосигналов
Islam et al. A computationally efficient near field broadband beamformer
Adrian-Ionut et al. A speed convergence Least Squares Constant Modulus Algorithm for smart antenna beamforming
CN109699069B (zh) 一种基于三维波束形成的联合多普勒频偏和载波频偏的估计方法

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210913

Address after: Block B, jingxinyuan, No. 25, beiwucun Road, Haidian District, Beijing 100088

Patentee after: HUIZHONGXINGZHI TECHNOLOGY (BEIJING) Co.,Ltd.

Address before: 611731, No. 2006, West Avenue, Chengdu hi tech Zone (West District, Sichuan)

Patentee before: University of Electronic Science and Technology of China

CP02 Change in the address of a patent holder
CP02 Change in the address of a patent holder

Address after: Room 209, building 1, East District, yard 10, northwest Wangdong Road, Haidian District, Beijing 100081

Patentee after: HUIZHONG XINGZHI TECHNOLOGY (BEIJING) CO.,LTD.

Address before: Block B, jingxinyuan, No. 25, beiwucun Road, Haidian District, Beijing 100088

Patentee before: HUIZHONG XINGZHI TECHNOLOGY (BEIJING) CO.,LTD.

CP02 Change in the address of a patent holder
CP02 Change in the address of a patent holder

Address after: 09 / F, 13 / F, building 683, zone 2, No. 5, Zhongguancun South Street, Haidian District, Beijing 100081

Patentee after: HUIZHONG XINGZHI TECHNOLOGY (BEIJING) CO.,LTD.

Address before: Room 209, building 1, East District, yard 10, northwest Wangdong Road, Haidian District, Beijing 100081

Patentee before: HUIZHONG XINGZHI TECHNOLOGY (BEIJING) CO.,LTD.