CN106990386A - 一种单偶极子极化敏感旋转阵列doa与极化参数联合估计方法 - Google Patents
一种单偶极子极化敏感旋转阵列doa与极化参数联合估计方法 Download PDFInfo
- Publication number
- CN106990386A CN106990386A CN201710239780.9A CN201710239780A CN106990386A CN 106990386 A CN106990386 A CN 106990386A CN 201710239780 A CN201710239780 A CN 201710239780A CN 106990386 A CN106990386 A CN 106990386A
- Authority
- CN
- China
- Prior art keywords
- array
- single dipole
- bay
- polarization
- vector
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/14—Systems for determining direction or deviation from predetermined direction
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明的目的在于提供一种单偶极子极化敏感旋转阵列DOA与极化参数联合估计方法,首先对所述极化敏感旋转阵列的阵列接收信号进行数学建模,之后对MUSIC算法的谱函数中的导向矢量的表达式进行修正,从而利用MUSIC算法实现入射信号DOA与极化参数的联合估计。本发明所述的旋转阵列不仅可以利用单偶极子构造极化敏感旋转阵列进行构造,还可以利用任意极化敏感天线单元或组合对其进行构造,具有很强的可移植性。该阵列可以有效降低***的通道数,极大的降低了***的成本,并且较少的阵元数有效的避免了阵元数较多带来的通道不一致性的问题。
Description
技术领域
本发明涉及的是一种阵列信号处理方法,具体地说是极化敏感阵列信号处理方法。
背景技术
基于极化敏感阵列的DOA与极化参数联合估计是阵列信号处理的一个重要研究方向,其主要任务是对极化敏感阵列接收到的信号的入射方向与极化状态进行估计。相对于基于标量阵列而言,极化敏感阵列可以同时敏感到入射信号的DOA与极化信息,增加了可利用的信息量,可以有效提高***性能。因此,如何利用极化敏感阵列对入射信号的DOA与极化参数进行估计,是将极化敏感阵列应用于实际测向***的基础。
阵列结构的选择直接影响着MUSIC算法的估计性能。均匀圆阵可以同时对入射信号的二维DOA信息进行估计,并且具有与方位角无关的波束宽度,是实际工程中应用常见的典型阵列之一。基于均匀圆阵的MUSIC算法的测向性能与阵列半径及阵元数相关,理论上,均匀圆阵的半径越大、阵元数越多,算法的估计精度越高。然而,在实际工程应用中,阵元数目的增多同时会带来通道不一致性增大的难题。并且阵元数的增多,意味着***的通道数的增多,以及接收***与信号处理***规模的增大,这极大的增加了***的成本。
发明内容
本发明的目的在于提供利用较少的阵元数估计入射信号的DOA与极化参数的一种单偶极子极化敏感旋转阵列DOA与极化参数联合估计方法。
本发明的目的是这样实现的:
本发明一种单偶极子极化敏感旋转阵列DOA与极化参数联合估计方法,其特征是:包括天线盘以及与天线盘边缘共形的单偶极子天线阵元1和单偶极子天线阵元2,天线盘以转速R匀速旋转;天线盘位于xoy平面内,将单偶极子天线阵元1与单偶极子天线阵元2分别放置于x轴的正半轴与负半轴,距离坐标原点的距离均为d;天线盘在旋转过程中,天线盘所在平面与旋转轴时刻保持垂直;
(1)在每个单偶极子天线旋转半个周期的过程中,按等时间间隔Δτ选取2M个阵列输出数据,构成虚拟单偶极子极化敏感均匀圆阵的接收数据矢量X=[X1,X2]T=ΦAS+N,S为K×1维的入射信号矢量,K为入射信源数,N=[N1,N2]T为2M×1维阵列接收噪声矢量,N1与N2分别为M×1维的阵元1与阵元2接收噪声矢量;
(2)根据接收数据矢量计算协方差矩阵Rx=E{XXH},E{·}表示求期望;
(3)对协方差矩阵Rx进行特征值分解,求得2M个特征值及对应特征向量;
(4)将Rx的2M个特征值进行降序排列,利用协方差矩阵Rx的后2M-K个小特征值对应特征向量张成噪声子空间UN;
(5)利用式构造二维DOA与极化参数联合估计的谱函数;
(6)改变谱函数中θ,φ,γ与η的取值,对构造的谱函数进行谱峰搜索,找出K个极大值点对应的DOA与极化参数的取值。
本发明还可以包括:
1、在阵列旋转过程中,选取的两个连续阵元之间的时延Δτ远大于接收机的数据采样周期,阵列旋转的转速满足:其中fs为接收机采样频率,L为快拍数。
2、X1与X2分别为由单偶极子天线阵元1与单偶极子天线阵元2的接收数据构成的M×1维阵列接收数据矢量;其中X1=Φ1A1S+N1,X2=Φ2A2S+N2,Φ=diag{Φ1,Φ2},Φ1=Φ2=diag{1,exp(-j2πfΔτ),…,exp(-j2πf(M-1)Δτ)},Φ1和Φ2分别表示单偶极子天线阵元1和单偶极子天线阵元2旋转后得到的附加相位矩阵;为2M×K维阵列流形矩,和分别表示单偶极子天线阵元1和单偶极子天线阵元2旋转后得到的M×K维阵列流形矩,为对应信源k(k=1,2,…K)的导向矢量,Uk为信源k入射到2M个虚拟阵元时的空域导向矢量矩阵,其计算式为B为虚拟阵列对应的极化敏感矩阵,其具体表达式为bm=[sin(2π(m-1)/2M),cos(2π(m-1)/2M)]为虚拟阵元m的极化敏感矢量;为对应信源k的极化空间域导向矢量。
本发明的优势在于:本发明通过对旋转中的阵元的输出数据进行采样,由于阵列的旋转,每个采样点所对应的阵元位置不同,从而构造出虚拟极化敏感均匀圆阵的接收数据矢量,对该虚拟极化敏感均匀圆阵进行数学建模,并利用MUSIC算法实现了入射信号DOA与极化参数的联合估计。极化敏感旋转阵列仅包含两个单偶极子天线,有效降低了接收***与信号处理***的规模,并且可以在一定程度上避免了通道不一致性对测向结果的影响。本发明的有益效果是:可以实现低成本、高精度的二维DOA与极化参数的联合估计。
本发明的优点如下:
首次提出利用单偶极子极化敏感天线构造极化敏感旋转阵列,对阵列旋转的转速及阵列的要求进行了说明;
首次提出利用所述极化敏感旋转阵列的二维DOA与极化参数联合估计方法,利用阵列的旋转构造虚拟极化敏感均匀圆阵,对该均匀圆阵的阵列模型进行数学建模,并利用MUSIC算法估计出入射信号的DOA与计划参数。
附图说明
图1为本发明的流程图;
图2为两天线阵列旋转示意图;
图3a为固定极化参数真值,当γ=40,η=50时,对应的空间谱图;图3b为固定极化参数真值,当γ=60,η=70时,对应的空间谱图;图3c为固定DOA真值,θ=125,φ=15时,对应的极化谱图;图3d为固定DOA真值,θ=55,φ=25时,对应的极化谱图;
图4为本发明方法与基于极化敏感阵列的联合参数估计算法通道不一致下的RMSE。
具体实施方式
下面结合附图举例对本发明做更详细地描述:
结合图1-4,本发明的一种由两个单偶极子天线构成的极化敏感旋转阵列,包括与天线盘边缘共形的两个单偶极子天线阵元,天线盘以一定的转速R(r/s)匀速旋转,其特征在于,天线盘位于xoy平面内,将单偶极子天线阵元1与阵元2分别放置于x轴的正半轴与负半轴,距离坐标原点的距离均为d;天线盘在旋转过程中天线盘所在平面与旋转轴时刻保持垂直;为了保证采样数据的稳定性,在阵列旋转过程中,选取的两个连续阵元之间的时延Δτ要远大于接收机的数据采样周期,且有进而可得阵列旋转的转速需要满足:其中fs为接收机采样频率,L为快拍数,2M为虚拟极化敏感均匀圆阵的阵元数。通常将转速R设置为10~20r/s,根据可得到时延Δτ的取值范围。
同时,本发明一种利用所述极化敏感旋转阵列进行二维DOA与极化参数联合估计的方法,具体步骤如下:
步骤(1):将两个单偶极子极化敏感天线分别放置于XOY平面内的(d,0)与(0,d)处,设置天线盘的转速为R,使天线盘以该转速逆时针匀速旋转,在旋转的同时保证天线盘与转轴垂直;
步骤(2):在每个单偶极子天线逆时针旋转半个周期的过程中,按等时间间隔Δτ选取2M(M为正整数)个阵列输出数据,构成虚拟单偶极子极化敏感均匀圆阵的接收数据矢量X=[X1,X2]T=ΦAS+N;
S为K×1维的入射信号矢量,K为入射信源数;
N=[N1,N2]T为2M×1维阵列接收噪声矢量,其中N1与N2分别为M×1维的阵元1与阵元2接收噪声矢量;
利用式Φ=diag{Φ1,Φ2}与Φ1=Φ2=diag{1,exp(-j2πfΔτ),…,exp(-j2πf(M-1)Δτ)}求得由于两单偶极子天线旋转而得到旋转阵列中的各阵元之间的采样延迟产生的附加相位矩阵Φ,其中,f为入射信号频率。
利用式与bm=[sin(2π(m-1)/2M),cos(2π(m-1)/2M)]求得虚拟极化敏感均匀圆阵的极化敏感矩阵B;利用式求得对应入射信号DOA与极化参数所有可能取值对应的极化空间域导向矢量;利用式求得对应入射信号DOA与极化参数所有可能取值对应的空域导向矢量矩阵;进而利用求得入射信号DOA与极化参数所有可能取值对应的导向矢量,用于后续的谱函数的构造。
步骤(2):根据接收数据矢量计算协方差矩阵Rx=E{XXH},E{·}表示求期望;
实际接收数据矩阵是有限长的,则可以利用最大似然估计计算协方差矩阵其中,L为快拍数。
步骤(3):对协方差矩阵Rx进行特征值分解,求得2M个特征值及对应特征向量;
对Rx进行特征值分解,可得
步骤(4):将Rx的2M个特征值进行降序排列,利用协方差矩阵Rx的后2M-K个小特征值对应特征向量张成噪声子空间UN;
将其特征值按降序排列为λ1≥λ2…≥λM,它们所对应的特征向量为ν1,ν2,…,νM,且各特征向量是相互正交的,则R的噪声子空间UN=[νK+1,νK+2,…,ν2M]。
步骤(5):利用式构造二维DOA与极化参数联合估计的谱函数,θ称为方位角,φ称为仰角,γ称为极化辅助角,η称为极化相位差;
步骤(6):改变谱函数中θ,φ,γ与η的取值,对构造的谱函数进行谱峰搜索,找出K个极大值点对应的DOA与极化参数的取值。
参照图2,是两天线阵列旋转示意图,空间水平位置上有单偶极子阵元1与阵元2两个阵元,两个阵元分别与圆周边缘共性,D个相互独立的信号源入射到半径为r的天线盘上,现假设基线1-2在XOY平面以Z轴为旋转轴按逆时针方向进行匀速旋转,阵元旋转周期为T,经过Δt(Δt<T)时间,阵元1、2分别旋转到阵元1'、2'的位置。
参照图3,是所述方法的估计谱图,其中图3a、图3b为对应极化参数真值的空间谱图,谱峰位置对应的方位角与仰角的取值即为入射信号的DOA估计值;图3c、图3d为对应DOA真值的极化谱图,谱峰位置对应的极化辅助角与极化相位差即为入射信号的极化参数估计值。
Claims (3)
1.一种单偶极子极化敏感旋转阵列DOA与极化参数联合估计方法,其特征是:包括天线盘以及与天线盘边缘共形的单偶极子天线阵元1和单偶极子天线阵元2,天线盘以转速R匀速旋转;天线盘位于xoy平面内,将单偶极子天线阵元1与单偶极子天线阵元2分别放置于x轴的正半轴与负半轴,距离坐标原点的距离均为d;天线盘在旋转过程中,天线盘所在平面与旋转轴时刻保持垂直;
(1)在每个单偶极子天线旋转半个周期的过程中,按等时间间隔Δτ选取2M个阵列输出数据,构成虚拟单偶极子极化敏感均匀圆阵的接收数据矢量X=[X1,X2]T=ΦAS+N,S为K×1维的入射信号矢量,K为入射信源数,N=[N1,N2]T为2M×1维阵列接收噪声矢量,N1与N2分别为M×1维的阵元1与阵元2接收噪声矢量;
(2)根据接收数据矢量计算协方差矩阵Rx=E{XXH},E{·}表示求期望;
(3)对协方差矩阵Rx进行特征值分解,求得2M个特征值及对应特征向量;
(4)将Rx的2M个特征值进行降序排列,利用协方差矩阵Rx的后2M-K个小特征值对应特征向量张成噪声子空间UN;
(5)利用式构造二维DOA与极化参数联合估计的谱函数;
(6)改变谱函数中θ,φ,γ与η的取值,对构造的谱函数进行谱峰搜索,找出K个极大值点对应的DOA与极化参数的取值。
2.根据权利要求1所述的一种单偶极子极化敏感旋转阵列DOA与极化参数联合估计方法,其特征是:
在阵列旋转过程中,选取的两个连续阵元之间的时延Δτ远大于接收机的数据采样周期,阵列旋转的转速满足:其中fs为接收机采样频率,L为快拍数。
3.根据权利要求2所述的一种单偶极子极化敏感旋转阵列DOA与极化参数联合估计方法,其特征是:
X1与X2分别为由单偶极子天线阵元1与单偶极子天线阵元2的接收数据构成的M×1维阵列接收数据矢量;其中X1=Φ1A1S+N1,X2=Φ2A2S+N2,Φ=diag{Φ1,Φ2},Φ1=Φ2=diag{1,exp(-j2πfΔτ),…,exp(-j2πf(M-1)Δτ)},Φ1和Φ2分别表示单偶极子天线阵元1和单偶极子天线阵元2旋转后得到的附加相位矩阵;为2M×K维阵列流形矩,和分别表示单偶极子天线阵元1和单偶极子天线阵元2旋转后得到的M×K维阵列流形矩,为对应信源k(k=1,2,…K)的导向矢量,Uk为信源k入射到2M个虚拟阵元时的空域导向矢量矩阵,其计算式为B为虚拟阵列对应的极化敏感矩阵,其具体表达式为bm=[sin(2π(m-1)/2M),cos(2π(m-1)/2M)]为虚拟阵元m的极化敏感矢量;为对应信源k的极化空间域导向矢量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710239780.9A CN106990386A (zh) | 2017-04-13 | 2017-04-13 | 一种单偶极子极化敏感旋转阵列doa与极化参数联合估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710239780.9A CN106990386A (zh) | 2017-04-13 | 2017-04-13 | 一种单偶极子极化敏感旋转阵列doa与极化参数联合估计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106990386A true CN106990386A (zh) | 2017-07-28 |
Family
ID=59415498
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710239780.9A Pending CN106990386A (zh) | 2017-04-13 | 2017-04-13 | 一种单偶极子极化敏感旋转阵列doa与极化参数联合估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106990386A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109738857A (zh) * | 2019-02-26 | 2019-05-10 | 中电科技扬州宝军电子有限公司 | 基于交叉偶极子阵列的非圆信号定位快速估计方法及装置 |
CN110007266A (zh) * | 2019-04-22 | 2019-07-12 | 哈尔滨工程大学 | 一种冲击噪声下的任意阵列相干源测向方法 |
CN110488097A (zh) * | 2019-07-02 | 2019-11-22 | 中国航空工业集团公司雷华电子技术研究所 | 基于线极化平台阵列的极化参数求解方法 |
CN110716171A (zh) * | 2019-08-28 | 2020-01-21 | 上海无线电设备研究所 | 一种基于遗传算法的极化doa联合估计方法 |
CN111065902A (zh) * | 2017-09-07 | 2020-04-24 | 德克萨斯仪器股份有限公司 | 使用量子分子旋转状态跃迁的压力感测 |
CN111698014A (zh) * | 2020-06-11 | 2020-09-22 | 上海磐启微电子有限公司 | 一种天线阵列*** |
US11258154B2 (en) | 2017-08-21 | 2022-02-22 | Texas Instruments Incorporated | Launch structures for a hermetically sealed cavity |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103323811A (zh) * | 2013-05-21 | 2013-09-25 | 西安电子科技大学 | 基于虚拟同心圆环阵列的参数估计方法 |
CN104375115A (zh) * | 2014-11-04 | 2015-02-25 | 中国航天科工集团第三研究院第八三五七研究所 | 极化敏感阵列下非圆信号doa与极化参数联合估计方法 |
CN104898085A (zh) * | 2015-05-14 | 2015-09-09 | 电子科技大学 | 一种极化敏感阵列参数估计的降维music算法 |
-
2017
- 2017-04-13 CN CN201710239780.9A patent/CN106990386A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103323811A (zh) * | 2013-05-21 | 2013-09-25 | 西安电子科技大学 | 基于虚拟同心圆环阵列的参数估计方法 |
CN104375115A (zh) * | 2014-11-04 | 2015-02-25 | 中国航天科工集团第三研究院第八三五七研究所 | 极化敏感阵列下非圆信号doa与极化参数联合估计方法 |
CN104898085A (zh) * | 2015-05-14 | 2015-09-09 | 电子科技大学 | 一种极化敏感阵列参数估计的降维music算法 |
Non-Patent Citations (2)
Title |
---|
司伟建 等: "基于阵列基线旋转的多目标 DOA 估计算法", 《***工程与电子技术》 * |
周欣: "基于电磁矢量传感器阵列的波达方向和极化参数同时估计的方法研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11258154B2 (en) | 2017-08-21 | 2022-02-22 | Texas Instruments Incorporated | Launch structures for a hermetically sealed cavity |
CN111065902A (zh) * | 2017-09-07 | 2020-04-24 | 德克萨斯仪器股份有限公司 | 使用量子分子旋转状态跃迁的压力感测 |
CN111065902B (zh) * | 2017-09-07 | 2022-02-01 | 德克萨斯仪器股份有限公司 | 使用量子分子旋转状态跃迁的压力感测 |
CN109738857A (zh) * | 2019-02-26 | 2019-05-10 | 中电科技扬州宝军电子有限公司 | 基于交叉偶极子阵列的非圆信号定位快速估计方法及装置 |
CN109738857B (zh) * | 2019-02-26 | 2020-09-01 | 中电科技扬州宝军电子有限公司 | 基于交叉偶极子阵列的非圆信号定位快速估计方法及装置 |
CN110007266A (zh) * | 2019-04-22 | 2019-07-12 | 哈尔滨工程大学 | 一种冲击噪声下的任意阵列相干源测向方法 |
CN110007266B (zh) * | 2019-04-22 | 2021-05-28 | 哈尔滨工程大学 | 一种冲击噪声下的任意阵列相干源测向方法 |
CN110488097A (zh) * | 2019-07-02 | 2019-11-22 | 中国航空工业集团公司雷华电子技术研究所 | 基于线极化平台阵列的极化参数求解方法 |
CN110488097B (zh) * | 2019-07-02 | 2021-10-12 | 中国航空工业集团公司雷华电子技术研究所 | 基于线极化平面阵列的极化参数求解方法 |
CN110716171A (zh) * | 2019-08-28 | 2020-01-21 | 上海无线电设备研究所 | 一种基于遗传算法的极化doa联合估计方法 |
CN111698014A (zh) * | 2020-06-11 | 2020-09-22 | 上海磐启微电子有限公司 | 一种天线阵列*** |
CN111698014B (zh) * | 2020-06-11 | 2023-09-26 | 上海磐启微电子有限公司 | 一种天线阵列*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106990386A (zh) | 一种单偶极子极化敏感旋转阵列doa与极化参数联合估计方法 | |
CN107015191B (zh) | 一种在多径干扰环境下单偶极子极化敏感阵列降维doa估计方法 | |
EP3254133B1 (en) | Direction finding using signal power | |
CN103091671B (zh) | 基于非同心电磁矢量阵列雷达的两维波达方向估计方法 | |
Nannuru et al. | Sparse Bayesian learning for beamforming using sparse linear arrays | |
CN102540138B (zh) | 一种多基线相位搜索式二维空间谱测向方法 | |
CN103018713B (zh) | 基于导航数字多波束接收阵列天线的卫星跟踪测角方法 | |
CN106483493B (zh) | 一种稀疏双平行线阵及二维波达方向估计方法 | |
CN105510874B (zh) | 一种近场信源多参数联合估计的降维music方法 | |
CN104699965B (zh) | 基于旋转干涉仪的近场源参数估计方法 | |
CN107919535B (zh) | 一种基于定向双圆阵的立体阵列天线及其构建方法 | |
CN103383452A (zh) | 分布式阵列目标到达角估计方法 | |
CN104360310A (zh) | 一种多目标近场源定位方法和装置 | |
CN112051540B (zh) | 一种快速高精度测向方法 | |
CN106526531A (zh) | 基于三维天线阵列的改进传播算子二维doa估计算法 | |
CN102830386A (zh) | 一种强干扰背景下的任意阵型微弱信号源角度估计方法 | |
CN106249196A (zh) | 三分量声矢量传感器稀疏阵列四元数解模糊方法 | |
CN109738853A (zh) | 一种电磁矢量互质阵基于旋转不变性的角度与极化估计方法 | |
CN106874642A (zh) | 基于旋转均匀圆阵的近场源参数估计解模糊方法 | |
CN103278796A (zh) | 锥面共形阵列多参数联合估计方法 | |
CN101252382B (zh) | 一种宽频段信号极化与doa估计方法及装置 | |
CN109696657A (zh) | 一种基于矢量水听器的相干声源定位方法 | |
CN110376547B (zh) | 基于二阶统计量的近场源定位方法 | |
CN215579031U (zh) | 一种旋转多基线相位干涉仪测向天线阵 | |
CN108318855A (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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170728 |