CN101931124A - 一种对数螺旋阵列天线布阵方法 - Google Patents
一种对数螺旋阵列天线布阵方法 Download PDFInfo
- Publication number
- CN101931124A CN101931124A CN2009102633900A CN200910263390A CN101931124A CN 101931124 A CN101931124 A CN 101931124A CN 2009102633900 A CN2009102633900 A CN 2009102633900A CN 200910263390 A CN200910263390 A CN 200910263390A CN 101931124 A CN101931124 A CN 101931124A
- Authority
- CN
- China
- Prior art keywords
- array
- antenna
- logarithmic spiral
- logarithmic
- spiral
- 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
Images
Landscapes
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
一种对数螺旋阵列天线布阵方法,天线阵采用N条对数螺旋布阵,单条对数螺旋线方程为,天线单元处于对数螺旋方程上,为对数螺旋线上任意一点在极坐标下的位置坐标,r为点到坐标原点的距离(单位:波长)即相对波长数,为相对于极轴(0°)的角度,a是对数螺旋线方程的参变量,为阵列设计时可优化的常数,每条对数螺旋上的天线单元n沿着螺旋线方程位于处,相邻单元间相位差为;N条螺旋线依次绕旋转轴旋转360°/N,各条对数螺旋线初始起点的r相同,每条螺旋线上首个天线单元初始起点的相位为,每条螺旋线上从初始起点依次相隔各分别放置m个阵元,阵列中心即原点处可放置一个阵元(亦可不放),当参数a、确定时,阵列流形就完全确定,N条对数螺旋阵列总的天线单元数为Nm+1。
Description
技术领域
本发明涉及超分辨空间谱(波达方向(Direction of Arrival,DOA))估计技术,尤其涉及一种对数螺旋阵列天线布阵方法,该方法可以实现高分辨二维波达方向(DOA)估计和目标测量、成像,可应用于雷达系统、通信系统和天文观测等领域。
背景技术
超分辨空间谱估计技术又称为波达方向(Direction of Arrival,DOA)估计技术,能够准确估计空间信号源的方向及空间分布,突破传统测向体制瑞利限的限制,在无线通信、雷达、声纳、导航、地震探测和医学等领域有着广泛的应用。空间谱的获取是通过天线阵列信号的迭加和抽取实现的,因此,阵列流形即阵列的几何结构、阵元的幅相分布、阵元间的耦合、以及周围散射体的影响等对空间谱估计至关重要。目前天线的分析和综合方法主要还是集中在天线的电磁辐射特性方面,较少关注对空间谱获取的影响;而在阵列信号处理理论的研究中,虽然对一维线阵的阵列流形问题有较深入的研究,如等距阵、非等距阵、虚拟阵列等方法,但对更加符合实际情况和更有针对性的二维阵列研究较少。常规的采用增大天线单元间距提高阵列分辨率的方法不可避免的形成栅瓣,从而出现虚假谱峰造成分辨率下降。
发明内容
为克服现有技术的缺陷,本发明提出了一种对数螺旋阵列天线布阵方法,通过一种新颖的二维天线阵列流形-对数螺旋阵,突破常规测向体制瑞利限的限制,实现高分辨率的二维DOA估计、目标定向和成像。
为实现上述目的,本发明采用以下技术方案:一种对数螺旋阵列天线布阵方法,其特征在于:天线阵采用N条对数螺旋布阵,(N=2,3,4...),单条对数螺旋线方程为 天线单元处于对数螺旋方程 上, 为对数螺旋线上任意一点在极坐标下的位置坐标,r为点到坐标原点的距离(单位:波长)即相对波长数, 为相对于极轴(0°)的角度,a是对数螺旋线方程的参变量,阵列设计时可根据天线方向图的波束宽度、副瓣电平等技术指标的需要进行优化,每条对数螺旋上的天线单元n沿着螺旋线方程位于 处,相邻单元间相位差为 N条螺旋线依次绕旋转轴旋转360°/N,各条对数螺旋线初始起点的r相同,每条螺旋线上首个天线单元初始起点的相位为 每条螺旋线上从初始起点相隔 各分别放置m个阵元,阵列中心即原点处放置一个阵元(亦可不放),当参数a、 确定时,阵列流形就完全确定,N条对数螺旋阵列总的天线单元数为Nm+1(阵列中心即原点处不放置阵元时为Nm)。
为了获得最佳的阵列布局,避免出现虚假谱峰,采用优化方法对天线单元的配置,即采用包括遗传算法(Genetic Algorithm,GA)、蚁群算法(Ant Colony Optimization,ACO)、粒子群优化算法(Particle Swarm Optimization,PSO)在内的各种算法,对影响天线方向图变化的参量a、 进行优化。
本发明的优点及显著效果:
1)天线阵采用多条(2、3、4...)对数螺旋布阵,增大了天线单元间距,获得了大的天线孔径,锐化了天线波束,提高了整个天线的分辨率。
2)由于采用对数螺旋阵,在x,y平面单元排列呈现非周期特性,因而,消除了栅瓣产生的条件,突破了半波长的瑞丽极限,在未出现虚假谱峰的基础上,提高了天线的分辨率,结合现代阵列信号处理技术,可以获得二维超分辨信号估计。
3)天线波束宽度窄、天线单元数少。天线波束宽度窄提高了天线的定向辐射和接收能力,提高了天线的角分辨能力,而与其他形式的阵列(如线阵、圆阵、二维平面阵)相比,其天线单元数显著下降,这样就大大减少了设备量,降低了建设和维护成本。
4)对数螺旋天线阵形式,其天线单元间距远大于半波长(瑞丽限),但经过优化布阵确能消除天线方向图中的栅瓣。天线的栅瓣(幅度接近主瓣电平的副瓣)是产生虚假谱峰,造成目标定向和成像模糊的根本原因。
5)对数螺旋天线阵可以方便的和现代阵列信号处理技术相结合,产生高分辨的功率谱,实现目标的精确定位和成像。在加大天线单元间距,增加天线综合孔径的同时,又削除了天线阵列的栅瓣,突破了传统测向体制瑞利限的限制,提高了天线阵的分辨率。仿真结果表明这种对数螺旋天线阵实现二维DOA估计的方法准确有效,并可将其推广应用于MIMO(multiple input multiple output)多天线移动通信系统、MIMO雷达系统和射电天文研究等领域。
附图说明
图1是本发明对数螺旋阵列示意图;
图2是螺旋阵列二维波达方向示意图;
图3是螺旋阵列天线方向图;
图4是螺旋阵列二维空间谱估计图;
图5是螺旋阵列相邻来波二维空间谱估计图;
图6是螺旋阵列二维空间谱估计(SNR=0dB)图。
具体实施方式
下面结合附图及具体实施例对本发明作进一步说明
天线阵根据天线单元的排列情况有很多种形式,如线阵、圆阵、二维平面阵等,本发明提出一种对数螺旋阵列并通过它实现信号的二维DOA估计。以三条螺旋阵列为例(如图1所示),图中,A1,A2,...A10是同一平面上天线单元所处的位置,除原点处A1外,各点满足如下方程:
其中r为螺旋线上任一点到坐标原点的距离(相对波长数), 为相对于极轴(0°)的角度, 为对数螺旋线上任意一点在极坐标下的位置坐标,a是对数螺旋线方程的参变量,阵列设计时可根据天线方向图的波束宽度、副瓣电平等技术指标的需要进行优化。N条螺旋线依次绕旋转轴旋转360°/N,三条螺旋线依次绕旋转轴旋转120°选择的控制常数,每条螺旋线上的天线单元n沿着螺旋线方程位于 处,相邻单元间相位差为 各条螺旋线初始起点的r相同,螺旋线上首个天线单元初始起点的相位为 每条螺旋线上从初始起点相隔 各分别放置m个阵元,阵列中心放置一个阵元,当参数a、 确定时,阵列流形就完全确定。三条螺旋阵列总的天线单元数为3m+1。由于平面螺旋阵单元数相对少,单元间距按指数展开大于半波长,使得天线孔径增大,分辨力提高。又因为单元排列的随机性破坏了出现栅瓣的周期性,所以经过优化阵列流形可以抑制栅瓣。相对于线阵,对数螺旋阵能提供俯仰、方位二维信息;相对于圆阵,对数螺旋阵能提供更低的副瓣电平和更高的分辨率。
为了获得最佳的阵列布局,避免出现虚假谱峰,可以采用优化方法对天线单元的配置,图2、3采用优化算法,如遗传算法(Genetic Algorithm,GA)、蚁群算法(Ant ColonyOptimization,ACO),粒子群优化算法(Particle Swarm Optimization,PSO)对对数螺旋面阵天线单元的位置参数进行综合优化影响天线方向图变化的参量a、 作为优化变量,优化的目标是在上半平面,天线方向图除主瓣区外,所有副瓣均低于-12dB,经PSO优化得到10单元对数螺旋阵(Log spiral Array,LSA)的优化参数:a=0.1783, 计算出10单元LSA阵的天线方向图如图3所示,天线波束宽度θ3dB=12.8°,可视范围内副瓣电平SLL<-12dB。
以下是对数螺旋阵天线的天线方向图仿真结果:
图4、5采用阵列信号处理MUSIC方法得到规一化的二维空间谱,与“谱峰”对应的所有θ和 即给出波达方向的估计。当波达方向为θ1=30°,φ1=250°,θ2=35°,φ2=120°, 快拍数为512,信噪比为10dB时,相应参数的对数螺旋阵列二维空间谱估计的仿真结果见图4。信号谱高于噪声信号27dB左右。当两信号源比较接近,波达方向为θ1=30°,φ1=145°;θ2=35°,φ2=130°,快拍数为512,信噪比为10dB时,相应参数的对数螺旋阵列二维空间谱估计仿真结果见图5。
图6当信噪比为0dB,即信号完全淹没于噪声中间时,依然可以获得清晰的二维空间谱。在下面的算例中,选取快拍数为256,波达方向为θ1=30°,φ1=145°;θ2=35°,φ2=130°,获得的二维空间谱如图6所示。由图可见,信号谱峰仍高于噪声信号10dB左右。
本发明对数螺旋阵的工作原理及工作过程:
设阵元数为N,其中N>P,(P为来波的数量,即目标数)原点指向各阵元的向量模值为r(n)(n=1,2,...,N),方位角为 各阵元的直角坐标为An(xn,yn),则
设波达方向俯仰角为θ1,θ2,...,θP,方位角为φ1,φ2,......φP。以位于坐标原点的阵元为基准,各信号源在基准点的复包络分别为s1(t),s2(t)...,sP(t)。如图1所示,第n个阵元在xy平面上的位置用矢量表示为 波达方向的单位矢量为(ui,vi,cosθi),ui=sinθicosφi,vi=sinθisinφi,θi为俯仰角,φi为方位角(i=1,2,...,P)。Dn在波达方向矢量上的投影为
对数螺旋阵的阵因子为:
阵因子方程中的幅度加权函数In(n=1,2,...,N)可以是相等的均匀阵,也可以是按某种分布变化的(如圆泰勒分布)。对于均匀阵,其激励幅度In=I0=1,激励相位αn=0。由于该阵列在x,y方向呈现非周期特性,因而,消除了栅瓣产生的条件。
Claims (3)
1.一种对数螺旋阵列天线布阵方法,其特征在于:天线阵采用N条对数螺旋布阵,N=2,3,4...,单条对数螺旋线方程为天线单元处于对数螺旋方程上,r,为对数螺旋线上任意一点在极坐标下的位置坐标,r为点到坐标原点的距离即相对波长数,为相对于极轴0°的角度,a是对数螺旋线方程的参变量,阵列设计时根据天线方向图的波束宽度、副瓣电平技术指标的需要进行优化,a>0.1;每条对数螺旋上的天线单元n沿着螺旋线方程位于处,相邻单元间相位差为N条螺旋线依次绕旋转轴旋转360°/N,各条对数螺旋线初始起点的r相同,每条螺旋线上首个天线单元初始起点的相位为每条螺旋线上从初始起点依次相隔各分别放置m个阵元,当参数a、确定时,阵列流形就完全确定,N条对数螺旋阵列总的天线单元数为Nm。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2009102633900A CN101931124A (zh) | 2009-12-18 | 2009-12-18 | 一种对数螺旋阵列天线布阵方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2009102633900A CN101931124A (zh) | 2009-12-18 | 2009-12-18 | 一种对数螺旋阵列天线布阵方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN101931124A true CN101931124A (zh) | 2010-12-29 |
Family
ID=43370153
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2009102633900A Pending CN101931124A (zh) | 2009-12-18 | 2009-12-18 | 一种对数螺旋阵列天线布阵方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101931124A (zh) |
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103513225A (zh) * | 2013-09-22 | 2014-01-15 | 西安电子科技大学 | 一种基于空间增益的稀疏平面阵形优化方法 |
| CN104267400A (zh) * | 2014-10-21 | 2015-01-07 | 内蒙古工业大学 | 用于mimo-sar成像的微波信号收发系统、方法及成像系统 |
| CN104409853A (zh) * | 2014-11-27 | 2015-03-11 | 中国船舶重工集团公司第七二四研究所 | 平面阵列天线电扫波束形状控制方法 |
| CN104900988A (zh) * | 2015-05-21 | 2015-09-09 | 电子科技大学 | 一种改变辅助栅格圆半径来设计稀布圆形天线阵列的方法 |
| CN106443571A (zh) * | 2016-09-05 | 2017-02-22 | 东南大学 | 一种保角变换旋转不变子空间波达方向估计方法 |
| CN107636896A (zh) * | 2015-03-05 | 2018-01-26 | 集美塔公司 | 用于圆柱形馈送天线的天线元件布置 |
| CN108256696A (zh) * | 2018-03-16 | 2018-07-06 | 电子科技大学 | 一种结合状态预测和粒子群优化的组网雷达天线配置方法 |
| CN105515629B (zh) * | 2015-12-23 | 2018-09-28 | 东南大学 | 一种基于改进粒子群优化的lte网络覆盖优化方法 |
| CN109612572A (zh) * | 2018-11-14 | 2019-04-12 | 国网上海市电力公司 | 用于快速识别高压电抗器异响声源位置的装置及方法 |
| CN110276127A (zh) * | 2019-06-20 | 2019-09-24 | 清华大学深圳研究生院 | 一种螺线型气动软体致动器静力学控制方法 |
| CN111969332A (zh) * | 2020-07-30 | 2020-11-20 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 大间距宽角扫描毫米波相控阵天线稀布阵方法 |
| CN113848555A (zh) * | 2021-09-30 | 2021-12-28 | 沈阳工业大学 | 一种用于对数螺旋阵的相干信号doa估计的预处理方法 |
| CN114491876A (zh) * | 2022-04-13 | 2022-05-13 | 杭州兆华电子股份有限公司 | 一种漏损声波检测阵列的自适应重构组阵方法及系统 |
| CN115377702A (zh) * | 2022-07-27 | 2022-11-22 | 香港理工大学深圳研究院 | 一种天线阵列优化方法、装置、终端及存储介质 |
| WO2023279737A1 (zh) * | 2021-07-06 | 2023-01-12 | 东南大学 | 一种多约束毫米波车载mimo雷达天线阵列设计方法与装置 |
| US20230194644A1 (en) * | 2021-12-17 | 2023-06-22 | The Boeing Company | Antennas for producing a variable phase response, angle-of-arrival sensors and methods for determining angle of arrival |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1453901A (zh) * | 2002-04-25 | 2003-11-05 | 哈里公司 | 串馈式螺旋卷阵列天线 |
| CN201364955Y (zh) * | 2009-02-12 | 2009-12-16 | 任国滨 | 螺旋阵列天线 |
-
2009
- 2009-12-18 CN CN2009102633900A patent/CN101931124A/zh active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1453901A (zh) * | 2002-04-25 | 2003-11-05 | 哈里公司 | 串馈式螺旋卷阵列天线 |
| CN201364955Y (zh) * | 2009-02-12 | 2009-12-16 | 任国滨 | 螺旋阵列天线 |
Non-Patent Citations (1)
| Title |
|---|
| 华光 等: "基于对数螺旋阵的波达方向估计", 《电波科学学报》 * |
Cited By (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103513225B (zh) * | 2013-09-22 | 2015-06-17 | 西安电子科技大学 | 一种基于空间增益的稀疏平面阵形优化方法 |
| CN103513225A (zh) * | 2013-09-22 | 2014-01-15 | 西安电子科技大学 | 一种基于空间增益的稀疏平面阵形优化方法 |
| CN104267400A (zh) * | 2014-10-21 | 2015-01-07 | 内蒙古工业大学 | 用于mimo-sar成像的微波信号收发系统、方法及成像系统 |
| CN104409853A (zh) * | 2014-11-27 | 2015-03-11 | 中国船舶重工集团公司第七二四研究所 | 平面阵列天线电扫波束形状控制方法 |
| CN107636896A (zh) * | 2015-03-05 | 2018-01-26 | 集美塔公司 | 用于圆柱形馈送天线的天线元件布置 |
| CN104900988A (zh) * | 2015-05-21 | 2015-09-09 | 电子科技大学 | 一种改变辅助栅格圆半径来设计稀布圆形天线阵列的方法 |
| CN105515629B (zh) * | 2015-12-23 | 2018-09-28 | 东南大学 | 一种基于改进粒子群优化的lte网络覆盖优化方法 |
| CN106443571A (zh) * | 2016-09-05 | 2017-02-22 | 东南大学 | 一种保角变换旋转不变子空间波达方向估计方法 |
| CN108256696B (zh) * | 2018-03-16 | 2021-10-26 | 电子科技大学 | 一种结合状态预测和粒子群优化的组网雷达天线配置方法 |
| CN108256696A (zh) * | 2018-03-16 | 2018-07-06 | 电子科技大学 | 一种结合状态预测和粒子群优化的组网雷达天线配置方法 |
| CN109612572A (zh) * | 2018-11-14 | 2019-04-12 | 国网上海市电力公司 | 用于快速识别高压电抗器异响声源位置的装置及方法 |
| CN110276127A (zh) * | 2019-06-20 | 2019-09-24 | 清华大学深圳研究生院 | 一种螺线型气动软体致动器静力学控制方法 |
| CN110276127B (zh) * | 2019-06-20 | 2023-04-07 | 清华大学深圳研究生院 | 一种螺线型气动软体致动器静力学控制方法 |
| CN111969332A (zh) * | 2020-07-30 | 2020-11-20 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 大间距宽角扫描毫米波相控阵天线稀布阵方法 |
| WO2023279737A1 (zh) * | 2021-07-06 | 2023-01-12 | 东南大学 | 一种多约束毫米波车载mimo雷达天线阵列设计方法与装置 |
| CN113848555A (zh) * | 2021-09-30 | 2021-12-28 | 沈阳工业大学 | 一种用于对数螺旋阵的相干信号doa估计的预处理方法 |
| US20230194644A1 (en) * | 2021-12-17 | 2023-06-22 | The Boeing Company | Antennas for producing a variable phase response, angle-of-arrival sensors and methods for determining angle of arrival |
| US12072427B2 (en) * | 2021-12-17 | 2024-08-27 | The Boeing Company | Antennas for producing a variable phase response, angle-of-arrival sensors and methods for determining angle of arrival |
| CN114491876A (zh) * | 2022-04-13 | 2022-05-13 | 杭州兆华电子股份有限公司 | 一种漏损声波检测阵列的自适应重构组阵方法及系统 |
| CN115377702A (zh) * | 2022-07-27 | 2022-11-22 | 香港理工大学深圳研究院 | 一种天线阵列优化方法、装置、终端及存储介质 |
| CN115377702B (zh) * | 2022-07-27 | 2025-04-22 | 香港理工大学深圳研究院 | 一种天线阵列优化方法、装置、终端及存储介质 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101931124A (zh) | 一种对数螺旋阵列天线布阵方法 | |
| CN106848546B (zh) | 一种宽带双极化天线阵列装置及高分辨测向方法 | |
| CN102904069B (zh) | 基于Sinc函数的圆形口径场分布的阵列天线的设计方法 | |
| CN112858994B (zh) | 一种基于均匀圆阵的比幅测向方法 | |
| CN106772260A (zh) | 基于凸优化算法的雷达阵列和差波束方向图优化方法 | |
| CN102608580B (zh) | 数字阵列超低副瓣自适应数字波束形成方法 | |
| CN110988868B (zh) | 一种等效分数阶模式涡旋电磁波产生及成像方法 | |
| Chen et al. | Two-dimensional direction of arrival estimation for improved archimedean spiral array with MUSIC algorithm | |
| CN111710992B (zh) | 一种基于等效线阵的弧形阵列天线的波束扫描方法及装置 | |
| Lu et al. | Optimum spatial arrangement of array elements for suppression of grating-lobes of radar cross section | |
| CN106329152A (zh) | 一种针对半球覆盖波束成形的阵列设计方法 | |
| CN106443571A (zh) | 一种保角变换旋转不变子空间波达方向估计方法 | |
| CN112505689B (zh) | 一种电磁涡旋雷达同心圆环阵列设计方法 | |
| Nechaev et al. | Evaluation of the influence of directivity factor of directive elements of conformal antenna arrays on the performances of azimuth-elevation DOA estimation | |
| CN112014835A (zh) | 分布式稀疏阵列雷达在栅瓣模糊下的目标跟踪方法和装置 | |
| CN109271735B (zh) | 基于量子启发式引力搜索算法的阵列方向图综合方法 | |
| CN102820527B (zh) | 一种雷达天线以及雷达系统 | |
| CN106990385B (zh) | 基于最小冗余线阵的非均匀面阵设计方法 | |
| Zou et al. | Beamforming with distortionless co-polarisation for conformal arrays based on geometric algebra | |
| CN102544732B (zh) | 一种近场相控阵波束聚焦方法 | |
| Riddolls | A Canadian perspective on Arctic and Polar over-the-horizon radar | |
| Ji et al. | Large-Scale Sparse Antenna Array Optimization for RCS Reduction with an AM-FCSN | |
| CN100380134C (zh) | 小型天线阵口径扩展与空间信号处理方法 | |
| Hansen | Conformal array realizations of the near and far fields of a complex point source | |
| Deng et al. | Frequency diverse arc array beampattern synthesis analysis with nonlinear frequency offset |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20101229 |