CN109343004B - 提高平面相控阵天线波束指向精度的迭代馈相计算方法 - Google Patents
提高平面相控阵天线波束指向精度的迭代馈相计算方法 Download PDFInfo
- Publication number
- CN109343004B CN109343004B CN201811084409.0A CN201811084409A CN109343004B CN 109343004 B CN109343004 B CN 109343004B CN 201811084409 A CN201811084409 A CN 201811084409A CN 109343004 B CN109343004 B CN 109343004B
- Authority
- CN
- China
- Prior art keywords
- phase
- array
- feed
- beam pointing
- actual
- 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.)
- Active
Links
Images
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/66—Radar-tracking systems; Analogous systems
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明公开了一种提高平面相控阵天线波束指向精度的迭代馈相计算方法。该方法首先以平面相控阵天线的中心对称点为原点确定各阵元坐标,将所有阵元根据中心对称特性依次编号;取前一半阵元,依次计算各阵元的理论馈相值和实际馈相值,实际馈相值由理论馈相值经进位量化或舍位量化得到,量化方式的选取遵循使已计算馈相阵元累积波束指向误差最小的原则;最后根据中心对称特性完成后一半阵元的馈相。本发明有效提高了平面相控阵天线的波束指向精度,同时相位计算的运算量和复杂度适中,满足相控阵天线***中快速波束扫描和跟踪的要求。
Description
技术领域
本发明属于相控阵天线领域,是一种提高平面相控阵天线波束指向精度的阵元迭代馈相计算方法。
背景技术
相控阵天线广泛应用于雷达***中,以实现空间快速扫描以及目标检测。相控阵天线的波束指向通过调整各阵元的激励电流来控制,包括激励电流的幅度和相位。移相器是调整阵元相位的关键器件,包括模拟移相器和数字移相器,其中,数字移相器因其集成度高、功耗低、移相值稳定等优点而被广泛采用。然而,数字移相器的馈相精度受移相器位数的限制,仅能提供最小相移步进整数倍的移相量,由此产生的实际馈相与理论馈相之间的误差称为馈相误差。馈相误差引起波束指向误差,即实际波束指向偏离预定波束指向。
目前,国内外有不少专家学者针对提高平面相控阵天线波束指向精度问题进行了深入研究。传统的馈相计算方法包括进位法、舍尾法和四舍五入法,然而这种方法会带来较大的周期性波束指向误差。郭燕昌等提出了随机馈相法和适当随机馈相法,包括预加相位法、二可能值法、适当预加相位法、适当二可能值法等,破坏了传统馈相方法的周期性误差,从统计角度消除了波束指向误差。但随机馈相法的馈相结果具有随机性和不稳定性,往往需要多次计算馈相值并从中选取性能最优的馈相方案。为了获得更加稳定的馈相结果,董亮等提出了基于智能优化算法的随机馈相计算方法,通过随机搜索和多次迭代获得全局较优解,然而该方法存在计算量大、计算复杂度高的问题,难以满足实际***中的快速馈相要求。Chris Hemmi等人提出了二维量化方法,该方法利用了与波束扫描方向正交的阵元的额外自由度来提高指向精度,但对一般平面阵列缺乏通用性。
发明内容
本发明的目的在于提供一种利用平面相控阵天线的中心对称性,实现高精度波束指向的迭代馈相计算方法。
实现本发明的技术解决方案为:一种提高平面相控阵天线波束指向精度的迭代馈相计算方法,包括以下步骤:
步骤1、以平面相控阵天线的中心对称点为原点,以阵面为xoy平面,垂直阵面向上为z轴建立空间直角坐标系,平面相控阵阵元总数为2N并依次编号为No.0~No.(2N-1),其中No.i阵元与No.(2N-1-i)阵元的坐标关于原点奇对称,0≤i≤(N-1);
步骤2、对于预定的波束指向和工作波长,取前一半阵元,即No.0~No.(N-1)阵元,从i=0至i=(N-1)依次计算No.i阵元的理论馈相值,并根据理论馈相值计算实际馈相值,具体为由理论馈相值经进位量化或舍位量化得到,量化方式的选取遵循使已计算馈相阵元累积波束指向误差最小的原则;
步骤3、根据中心对称特性完成后一半阵元的馈相计算,即No.N~No.(2N-1)阵元,具体为:No.i阵元的实际馈相值与其关于原点成中心对称的No.(2N-1-i)阵元的实际馈相值互为相反数,N≤i≤(2N-1)。
本发明与现有技术相比,其显著优点为:(1)本发明为确定性的馈相计算方法,馈相结果稳定;(2)本发明通过对阵元的依次迭代馈相,提高了平面相控阵的波束指向精度;(3)本发明相位计算方法的复杂度和计算量适中,满足相控阵***中快速波束扫描和跟踪的要求。
下面结合附图对本发明作进一步详细描述。
附图说明
图1是本发明提高平面相控阵波束指向精度的阵元迭代馈相计算流程图。
图2是平面相控阵阵元排布图。
图3是空间坐标系定义及离轴角、旋转角、波束指向误差定义图。
图4是阵元编号示意图。
图5是本方法馈相后的天线方向图,其中图5(a)中配置预定波束指向为(0°,0°),图5(b)中配置预定波束指向为(40°,40°)。
具体实施方式
一种提高平面相控阵天线波束指向精度的迭代馈相计算方法,包括以下步骤:
步骤1、以平面相控阵天线的中心对称点为原点,以阵面为xoy平面,垂直阵面向上为z轴建立如图3所示空间直角坐标系,平面相控阵阵元总数为2N并依次编号为No.0~No.(2N-1),其中No.i阵元与No.(2N-1-i)阵元的坐标关于原点奇对称,0≤i≤(N-1);
步骤2、对于预定的波束指向和工作波长,取前一半阵元,即No.0~No.(N-1)阵元,从i=0至i=(N-1)依次计算No.i阵元的理论馈相值,并根据理论馈相值计算实际馈相值,具体为由理论馈相值经进位量化或舍位量化得到,量化方式的选取遵循使已计算馈相阵元累积波束指向误差最小的原则;
前一半阵元,即No.0~No.(N-1)阵元,从i=0至i=(N-1)依次计算No.i阵元的理论馈相值,并根据理论馈相值计算实际馈相值的具体步骤为:
步骤2-1、计算No.i阵元的理论馈相值,具体为
步骤2-2、No.i阵元的理论馈相值按进位法或舍位法量化为数字移相器最小移相步进的整数倍,从而得到实际馈相值φi,具体为:
No.i阵元的馈相误差表示为:
δφi=φi-φtheo,i
理论馈相值按进位法量化得到的实际馈相值和馈相误差分别记为φcarry,i,δφcarry,i,按舍位法量化得到的实际馈相值和馈相误差分别记为φtruncate,i,δφtruncate,i;
步骤2-3、计算由No.0~No.i阵元及其关于原点成中心对称的No.(2N-1-i)~No.(2N-1)阵元导致的累积波束指向误差γi,具体为:
根据中心对称特性,No.i阵元的实际馈相值与其关于原点成中心对称的No.(2N-1-i)阵元的实际馈相值互为相反数;由No.0~No.i阵元及其关于原点成中心对称的No.(2N-1-i)~No.(2N-1)阵元导致的累积波束指向误差γi表示为:
按进位法量化φtheo,i导致的累积波束指向误差记为γcarry,i,按舍位法量化φtheo,i导致的累积波束指向误差记为γtruncate,i。
步骤2-4、选择使已计算馈相阵元导致的累积波束指向误差最小的量化方法作为φtheo,i的实际量化方法,具体为:
比较γcarry,i和γtruncate,i,如果γcarry,i<γtruncate,i,则No.i阵元的实际馈相值由φtheo,i进位量化得到,即否则No.i阵元的实际馈相值由φtheo,i舍位量化得到,即
步骤3、根据中心对称特性完成后一半阵元的馈相计算,即No.N~No.(2N-1)阵元,具体为:No.i阵元的实际馈相值与其关于原点成中心对称的No.(2N-1-i)阵元的实际馈相值互为相反数,N≤i≤(2N-1)。
实施例1
如图1所示,本发明提高平面相控阵天线波束指向精度的迭代馈相计算方法,步骤如下:
仿真条件:具有中心对称特性的平面相控阵天线阵元排布如图2所示,由182(14×13)个阵元组成的三角形栅格平面相控阵,平面相控阵的中心对称点为点A,任意相邻阵元间的距离d=0.56λ,λ为工作波长,每个阵元后接一个6位数字移相器,各阵元激励电流的幅度为1。
步骤1、首先以平面相控阵的中心对称点A为原点,以阵面为xoy平面,垂直阵面向上为z轴建立如图3所示空间直角坐标系。平面相控阵总阵元数为182并将其依次如图4所示编号为No.0~No.181,其中No.i阵元与No.(181-i)阵元的坐标关于原点奇对称,0≤i≤90。
(1)计算No.i阵元的理论馈相值φtheo,i,通过公式:
(2)将No.i阵元的理论馈相值按进位法或舍位法量化为数字移相器最小移相步进的整数倍,从而得到实际馈相值φi,具体为:
按进位法量化时,No.i阵元的实际馈相值φcarry,i和馈相误差δφcarry,i为
δφcarry,i=φcarry,i-φtheo,i
按舍位法量化时,No.i阵元的实际馈相值φtruncate,i和馈相误差δφtruncate,i为
δφtruncate,i=φtruncate,i-φtheo,i
其中数字移相器的最小移相步进Δ=2π/26。
(3)计算由No.0~No.i阵元及其关于原点成中心对称的No.(181-i)~No.181阵元导致的累积波束指向误差γi,其中编号No.0~No.(i-1)阵元及其关于原点成中心对称的阵元的实际馈相值和馈相误差已确定。
分别计算按进位法量化φtheo,i导致的累积波束指向误差γcarry,i和按舍位法量化φtheo,i导致的累积波束指向误差γtruncate,i,通过公式
其中
(4)比较γcarry,i和γtruncate,i,如果γcarry,i<γtruncate,i,则No.i阵元的实际馈相值由φtheo,i进位量化得到,即否则No.i阵元的实际馈相值由φtheo,i舍位量化得到,即
(5)i=i+1,重复步骤2-1~2-4,计算No.i阵元的实际馈相值,直到i=N-1。
步骤3、根据中心对称特性完成后一半阵元的馈相计算,即No.91~No.181阵元,具体为:No.i阵元与其关于原点成中心对称的No.(181-i)阵元的实际馈相值互为相反数,φi=-φ181-i,i=91,92,...,181。
通过计算,当预定波束指向分别配置为(0°,0°)和(40°,40°)时,经本方法馈相后的天线方向图如图5所示,其实际波束指向分别为(0°,0°)和(39.99°,40°),预定波束指向与实际波束指向的离轴角偏差在0.01°以内,预定波束指向与实际波束指向的旋转角偏差也在0.01°以内。
由上可知,本发明能有效提高平面相控阵天线的波束指向精度,同时相位计算方法的复杂度适中,满足相控阵***中快速波束扫描和跟踪的要求。
Claims (3)
1.一种提高平面相控阵天线波束指向精度的迭代馈相计算方法,其特征在于,包括以下步骤:
步骤1、以平面相控阵天线的中心对称点为原点,以阵面为xoy平面,垂直阵面向上为z轴建立空间直角坐标系,平面相控阵阵元总数为2N并依次编号为No.0~No.(2N-1),其中No.i阵元与No.(2N-1-i)阵元的坐标关于原点奇对称,0≤i≤(N-1);
步骤2、对于预定的波束指向和工作波长,取前一半阵元,即No.0~No.(N-1)阵元,从i=0至i=(N-1)依次计算No.i阵元的理论馈相值,并根据理论馈相值计算实际馈相值,具体为由理论馈相值经进位量化或舍位量化得到,量化方式的选取遵循使已计算馈相阵元累积波束指向误差最小的原则,步骤2中前一半阵元,即No.0~No.(N-1)阵元,从i=0至i=(N-1)依次计算No.i阵元的理论馈相值,并根据理论馈相值计算实际馈相值,具体步骤为:
步骤2-1、计算No.i阵元的理论馈相值,具体为:
步骤2-2、No.i阵元的理论馈相值按进位法或舍位法量化为数字移相器最小移相步进的整数倍,从而得到实际馈相值φi,具体为:
No.i阵元的馈相误差表示为:
δφi=φi-φtheo,i
理论馈相值按进位法量化得到的实际馈相值和馈相误差分别记为φcarry,i,δφcarry,i,按舍位法量化得到的实际馈相值和馈相误差分别记为φtruncate,i,δφtruncate,i;
步骤2-3、计算由No.0~No.i阵元及其关于原点成中心对称的No.(2N-1-i)~No.(2N-1)阵元导致的累积波束指向误差γi,具体为:
根据中心对称特性,No.i阵元的实际馈相值与其关于原点成中心对称的No.(2N-1-i)阵元的实际馈相值互为相反数;由No.0~No.i阵元及其关于原点成中心对称的No.(2N-1-i)~No.(2N-1)阵元导致的累积波束指向误差γi表示为:
按进位法量化φtheo,i导致的累积波束指向误差记为γcarry,i,按舍位法量化φtheo,i导致的累积波束指向误差记为γtruncate,i;
步骤2-4、选择使已计算馈相阵元导致的累积波束指向误差最小的量化方法作为φtheo,i的实际量化方法,具体为:
比较γcarry,i和γtruncate,i,如果γcarry,i<γtruncate,i,则No.i阵元的实际馈相值由φtheo,i进位量化得到,即否则No.i阵元的实际馈相值由φtheo,i舍位量化得到,即
步骤3、根据中心对称特性完成后一半阵元的馈相计算,即No.N~No.(2N-1)阵元,具体为:No.i阵元的实际馈相值与其关于原点成中心对称的No.(2N-1-i)阵元的实际馈相值互为相反数,N≤i≤(2N-1)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811084409.0A CN109343004B (zh) | 2018-09-18 | 2018-09-18 | 提高平面相控阵天线波束指向精度的迭代馈相计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811084409.0A CN109343004B (zh) | 2018-09-18 | 2018-09-18 | 提高平面相控阵天线波束指向精度的迭代馈相计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109343004A CN109343004A (zh) | 2019-02-15 |
CN109343004B true CN109343004B (zh) | 2023-01-20 |
Family
ID=65305842
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811084409.0A Active CN109343004B (zh) | 2018-09-18 | 2018-09-18 | 提高平面相控阵天线波束指向精度的迭代馈相计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109343004B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110940310B (zh) * | 2019-11-15 | 2021-12-07 | 北京遥测技术研究所 | 弹箭载中继测控终端相控阵天线波束指向角度计算方法 |
CN112928495B (zh) * | 2021-01-21 | 2022-09-20 | 南京理工大学 | 基于fpga的k位数字移相器高精度顺序馈相方法 |
CN114491876B (zh) * | 2022-04-13 | 2022-07-19 | 杭州兆华电子股份有限公司 | 一种漏损声波检测阵列的自适应重构组阵方法及*** |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102405603A (zh) * | 2009-04-21 | 2012-04-04 | 马维尔国际贸易有限公司 | 具有选择性波束衰减的多点机会波束成形 |
CN102810742A (zh) * | 2012-08-01 | 2012-12-05 | 四川九洲电器集团有限责任公司 | 一种提高相控阵波束指向精度的方法 |
CN106096246A (zh) * | 2016-06-06 | 2016-11-09 | 华北电力大学(保定) | 基于pm2.5和pm10的气溶胶光学厚度估计方法 |
CN107066786A (zh) * | 2016-11-22 | 2017-08-18 | 深圳职业技术学院 | 基于神经网络的气溶胶光学厚度反演算法 |
CN107329134A (zh) * | 2017-06-29 | 2017-11-07 | 电子科技大学 | 一种基于阵元馈电波形控制的波控阵超宽带雷达天线阵列 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9964482B2 (en) * | 2015-01-26 | 2018-05-08 | The United States Of America As Represented By The Secretary Of The Army | Detecting clouds using polarized sunlight |
CN105674901B (zh) * | 2016-04-08 | 2017-12-26 | 电子科技大学 | 基于统计学分段的大气气溶胶反演方法 |
-
2018
- 2018-09-18 CN CN201811084409.0A patent/CN109343004B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102405603A (zh) * | 2009-04-21 | 2012-04-04 | 马维尔国际贸易有限公司 | 具有选择性波束衰减的多点机会波束成形 |
CN102810742A (zh) * | 2012-08-01 | 2012-12-05 | 四川九洲电器集团有限责任公司 | 一种提高相控阵波束指向精度的方法 |
CN106096246A (zh) * | 2016-06-06 | 2016-11-09 | 华北电力大学(保定) | 基于pm2.5和pm10的气溶胶光学厚度估计方法 |
CN107066786A (zh) * | 2016-11-22 | 2017-08-18 | 深圳职业技术学院 | 基于神经网络的气溶胶光学厚度反演算法 |
CN107329134A (zh) * | 2017-06-29 | 2017-11-07 | 电子科技大学 | 一种基于阵元馈电波形控制的波控阵超宽带雷达天线阵列 |
Also Published As
Publication number | Publication date |
---|---|
CN109343004A (zh) | 2019-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109343004B (zh) | 提高平面相控阵天线波束指向精度的迭代馈相计算方法 | |
Keizer | Linear array thinning using iterative FFT techniques | |
CN106407723B (zh) | 面向低副瓣的稀疏排布阵列天线激励电流幅度的确定方法 | |
CN109633526B (zh) | 基于方向函数的非均匀圆阵相位干涉仪测向解模糊方法 | |
CN108896954B (zh) | 互质阵中一种基于联合实值子空间的波达角估计方法 | |
CN107944133B (zh) | 基于多目标量子蜘蛛群演化机制的环形天线阵列稀疏方法 | |
CN106483493B (zh) | 一种稀疏双平行线阵及二维波达方向估计方法 | |
CN107658573B (zh) | 一种基于量子搜寻者搜索机制的圆环阵方向图综合方法 | |
CN107918108A (zh) | 一种均匀圆阵二维波达方向快速估计方法 | |
CN111985145B (zh) | 一种大间距相控阵天线栅瓣抑制方法及抑制*** | |
CN110632559B (zh) | 一种基于应用特征的mimo脉冲雷达波形设计与优化方法 | |
CN104931919A (zh) | 一种二维平面圆环阵列的测向方法 | |
CN105301564A (zh) | 一种基于双线性拟合的相控阵天线波束指向误差补偿方法 | |
CN106096218A (zh) | 一种移动卫星通信稀布平面天线阵列优化方法 | |
CN110232228B (zh) | 一种多联骨牌非规则子阵阵面优选设计方法 | |
CN104950290A (zh) | 基于加权k均值聚类的大型相控阵天线子阵划分方法 | |
CN103777179A (zh) | 一种用于雷达立体共形阵列的子阵降维方法 | |
CN107238829A (zh) | 一种双层圆环阵稳健超指向性波束形成方法 | |
CN115037352A (zh) | 一种面向卫星波束对准增强的降噪方法 | |
CN111914427B (zh) | 一种基于面积归一化策略的多约束矩形阵列稀布优化方法 | |
Li et al. | Novel subarray partition algorithm for solving the problem of too low beam collection efficiency caused by dividing a few subarrays | |
CN110557180B (zh) | 一种多分查找波束成形方法、装置、***及存储介质 | |
CN109670228B (zh) | 一种阵列天线单元方向图快速插值计算方法 | |
CN109190193B (zh) | 基于果蝇算法的子阵级相控阵天线方向图综合方法 | |
CN111537948A (zh) | 一种基于esprit算法的分离式长电偶极子的角度估计方法 |
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 |