CN113806685A - 阵列天线方向图形成的快速方法 - Google Patents
阵列天线方向图形成的快速方法 Download PDFInfo
- Publication number
- CN113806685A CN113806685A CN202111142299.0A CN202111142299A CN113806685A CN 113806685 A CN113806685 A CN 113806685A CN 202111142299 A CN202111142299 A CN 202111142299A CN 113806685 A CN113806685 A CN 113806685A
- Authority
- CN
- China
- Prior art keywords
- directional diagram
- array
- field distribution
- complex field
- array antenna
- 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
Links
- 238000010586 diagram Methods 0.000 title claims abstract description 76
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000001514 detection method Methods 0.000 claims abstract description 13
- 238000003384 imaging method Methods 0.000 claims abstract description 10
- 238000004891 communication Methods 0.000 claims abstract description 6
- 238000012634 optical imaging Methods 0.000 claims abstract description 6
- 239000000126 substance Substances 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000005286 illumination Methods 0.000 claims description 3
- 230000003287 optical effect Effects 0.000 claims 1
- 238000012285 ultrasound imaging Methods 0.000 claims 1
- 238000004458 analytical method Methods 0.000 abstract description 3
- 230000005611 electricity Effects 0.000 abstract description 3
- 230000001678 irradiating effect Effects 0.000 abstract 1
- 230000015572 biosynthetic process Effects 0.000 description 7
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 230000035939 shock Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000007261 regionalization Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
- G06F17/142—Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Discrete Mathematics (AREA)
- Computing Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明涉及光学成像、微波成像、雷达探测、声呐、超声成像以及基于声、光、电等媒介的目标探测和无线通信技术领域,具体涉及一种阵列天线方向图形成的快速方法及其在上述各领域中的应用。本发明方法通过将阵列天线方向图分解为一系列不同方向、相位、强度的照射平面波的集合,计算出该平面波集合照射在阵列口径上的复场分布,该结果即为所需方向图对应的阵列天线单元幅度、相位加权系数。本发明采用非迭代的解析方法,可以快速模拟所需要的方向图,满足实时应用的要求,另外,本方法运算量小,对硬件环境和计算资源的配置要求低,硬件价格和运行成本得到了有效的控制。
Description
技术领域
本发明涉及光学成像、微波成像、雷达探测、声呐、超声成像以及基于声、光、电等媒介的目标探测和无线通信技术领域,具体涉及一种阵列天线方向图形成的快速方法及其在上述各领域中的应用。
背景技术
现有的阵列天线复杂方向图形成方法大部分采用优化方法,即利用各种优化算法对阵列单元的幅度、相位分布进行优化,进而合成所需的复杂形状方向图。
然而,上述优化方法具有明显的缺陷和不足,一方面,目前采用的优化算法大部分都是迭代类算法,需要经过多次迭代后才能获得符合要求的方向图,因此算法速度慢、实时性差,很难满足特定应用场景对实时复杂方向图的形成要求;另一方面,迭代类算法运算量大,对硬件环境和计算资源的配置要求高,故而造成硬件价格和运行成本均较高。
发明内容
为了克服现有阵列天线复杂方向图形成方法的上述缺陷和不足,本发明提供了一套解决方案:本发明采用非迭代的解析方法,可以快速模拟所需要的方向图,本方法运算量小,对硬件环境和计算资源的配置要求低,可满足实时应用的要求。
阵列天线理论表明,对阵列分布的复信号进行二维IFFT变换即可获得其远场方向图。而成像理论则证明,在远距离情况下,对阵列分布的复信号进行二维IFFT变换,其结果即为目标的像,综合二者可知,远场方向图与目标的像是一致的。
在上述理论的基础上,本发明人发现,对于复杂形状的阵列方向图,用一系列不同强度、不同入射方向、不同初始相位的平面波进行照射和叠加,即可获得所需的方向图。
方向图合成公式为:
=∑IFFT2D(W),
其中,G表示方向图,W表示单个平面波照射形成的阵列复场分布,符号IFFT2D表示二维快速傅里叶逆变换,符号∑表示求和运算。
根据傅里叶变换的线性性质,上式可改写为:
=IFFT2D{∑W},
具体实施时,先合成一个基本形状方向图G0,然后再进行空间扫描改变波束位置,并用基本形状方向图与冲击函数卷积的方式实现波束方向改变:
从而有:
其中,W0为基本形状方向图对应的阵列复场分布,Ws为扫描方向图对应的复场分布。
上式表明,基本形状方向图的阵列复场分布与针状扫描波束方向图的复场分布进行相乘,即可获得所需复杂形状的波束扫描方向图。
在上述认识的基础上,本发明提供了一种阵列天线方向图形成的快速方法,该方法通过将阵列天线方向图分解为一系列不同方向、相位、强度的照射平面波的集合,计算出该平面波集合照射在阵列口径上的复场分布,该结果即为所需方向图对应的阵列天线单元幅度、相位加权系数。
具体而言,本发明阵列天线方向图形成的快速方法其算法流程如下:
步骤一:对阵列天线方向图进行分解,获得照射波束中心指向以及基本形状方向图;
步骤二:对基本形状方向图进行分解,获得照射平面波集合;
步骤三:根据照射平面波集合,计算阵列口径处的基本形状方向图复场分布;
步骤四:根据照射波束中心指向,计算阵列口径处的扫描方向图复场分布;
步骤五:计算基本形状方向图复场分布与扫描方向图复场分布的乘积;
步骤六:根据基本形状方向图复场分布与扫描方向图复场分布的乘积结果,输出阵列天线单元幅度、相位加权系数。
进一步地,本发明方法步骤一中通过分解阵列天线方向图,获得照射波束中心指向以及基本形状方向图,其中:
所述基本形状方向图G0,是指不考虑照射波束中心指向时的阵列天线方向图形状,即阵列天线方向图中心方向为法向方向时的形状。
进一步地,本发明方法步骤二中对基本形状方向图进行分解,获得照射平面波集合,其具体步骤如下:
根据阵列天线同相馈电时默认针状波束的宽度,将基本形状方向图的图像离散化为间隔半功率波束宽度的点目标集合,每个点目标代表了一个照射平面波,从而获得了N个照射平面波形成的集合其中为代表入射平面波方向的球坐标角度参数,σ代表入射平面波强度,φ代表入射平面波相位;另外,点目标的角度坐标即为照射平面波的入射方向,点目标的强度代表了照射平面波的强度,当两个相邻的点目标之间无零陷时,对应的两个相邻照射平面波的相位差为零,当两个相邻的点目标之间有零陷时,对应的两个相邻照射平面波的相位差为180°。
进一步地,本发明方法步骤三中根据照射平面波集合,计算阵列口径处的基本形状方向图复场分布,其算法流程如下:
平面波到达阵列口径处的复场为:
E=σej(ψ+φ),
其中,j为虚数单位,e为欧拉常数;
基本形状方向图对应的复场分布为:
进一步地,本发明方法步骤四中根据照射波束中心指向,计算阵列口径处的扫描方向图复场分布,其算法流程如下:
扫描方向图复场分布为:
当无中心波束指向或存在多个针状波束的情况下,则直接选取:
Es=1。
进一步地,本发明方法步骤五中计算基本形状方向图复场分布与扫描方向图复场分布的乘积,算法如下:
E=E0Es;
步骤六中根据基本形状方向图复场分布与扫描方向图复场分布的乘积结果,输出阵列天线单元幅度、相位加权系数,具体算法如下:
WA=|E|,
此外,本发明还涉及上述方法在光学成像、微波成像、雷达探测、声呐、超声成像以及声、光、电目标探测和无线通信领域中的应用。
综上,本发明阵列天线方向图形成的快速方法具有以下优点:
(1)采用非迭代的解析方法,可以快速模拟所需要的方向图,可满足实时应用的要求。
(2)本方法运算量小,对硬件环境和计算资源的配置要求低,硬件价格和运行成本得到了有效的控制。
(3)具有良好的应用前景,可广泛应用于以声、光、电等为媒介的目标探测和无线通信技术领域,当探测媒介为电磁波时,本技术适用于微波成像、雷达探测、无线通信、合成孔径雷达、逆合成孔径雷达;当探测媒介为声波、超声波时,本技术适用于声呐、超声成像、合成孔径声呐;当探测媒介为光时,本技术适用于光学成像、合成孔径光学成像。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对本发明实施例描述中所需要使用的附图作简要介绍,显而易见地,以下附图仅仅是本发明中记载的一些实施例,对于本领域技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明方法的算法流程图。
图2为利用本发明方法获得的环形波束方向图。
图3为利用本发明方法获得的扇形波束方向图。
图4为利用本发明方法获得的差波束方向图。
图5为利用本发明方法获得的多波束方向图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合具体实施例及相应的附图对本发明的技术方案进行清楚、完整地描述。显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例,本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。
同时,应理解,本发明的保护范围并不局限于下述特定的具体实施方案;还应当理解,本发明实施例中使用的术语是为了描述特定的具体实施方案,而不是为了限制本发明的保护范围。
实施例1:一种阵列天线方向图形成的快速方法(参见附图1),包括:
步骤一:对阵列天线方向图进行分解,获得照射波束中心指向以及基本形状方向图,其中:
所述基本形状方向图G0,是指不考虑照射波束中心指向时的阵列天线方向图形状,即阵列天线方向图中心方向为法向方向时的形状。
步骤二:对基本形状方向图进行分解,获得照射平面波集合,其具体步骤如下:
根据阵列天线同相馈电时默认针状波束的宽度,将基本形状方向图的图像离散化为间隔半功率波束宽度的点目标集合,每个点目标代表了一个照射平面波,从而获得了N个照射平面波形成的集合其中为代表入射平面波方向的球坐标角度参数,σ代表入射平面波强度,φ代表入射平面波相位;另外,点目标的角度坐标即为照射平面波的入射方向,点目标的强度代表了照射平面波的强度,当两个相邻的点目标之间无零陷时,对应的两个相邻照射平面波的相位差为零,当两个相邻的点目标之间有零陷时,对应的两个相邻照射平面波的相位差为180°。
步骤三:根据照射平面波集合,计算阵列口径处的基本形状方向图复场分布,其算法流程如下:
平面波到达阵列口径处的复场为:
E=σej(ψ+φ),
其中,j为虚数单位,e为欧拉常数;
基本形状方向图对应的复场分布为:
步骤四:根据照射波束中心指向,计算阵列口径处的扫描方向图复场分布,其算法流程如下:
扫描方向图复场分布为:
当无中心波束指向或存在多个针状波束的情况下,则直接选取:
Es=1。
步骤五:计算基本形状方向图复场分布与扫描方向图复场分布的乘积,算法如下:
E=E0Es;
步骤六:根据基本形状方向图复场分布与扫描方向图复场分布的乘积结果,输出阵列天线单元幅度、相位加权系数,具体算法如下:
WA=|E|,
实施例2:环形波束方向图的快速形成(参见附图2)
利用实施例1的方法,快速形成环形波束方向图
工作频率30GHz,天线单元间距为半波长,阵列规模50*50,合成环形波束,波束中心指向为俯仰30°,方向图结果见图2。
实施例3:扇形波束方向图的快速形成(参见附图3)
利用实施例1的方法,快速形成扇形波束方向图
工作频率30GHz,天线单元间距为半波长,阵列规模50*50,合成扇形波束,波束中心指向为方位20°,方向图结果见图3。
实施例4:差波束方向图的快速形成(参见附图4)
利用实施例1的方法,快速形成差波束方向图
工作频率30GHz,天线单元间距为半波长,阵列规模50*50,合成差波束,波束中心指向为方位20°,方向图结果见图4。
实施例5:多波束方向图的快速形成(参见附图5)
利用实施例1的方法,快速形成多波束方向图
工作频率30GHz,天线单元间距为半波长,阵列规模50*50,合成多个波束,方向图结果见图5。
本发明中的各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似的部分互相参见即可。
以上所述仅为本发明的实施例而已,并不用于限制本发明。对于本领域技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原理之内所作的任何修改、替换等,均应包含在本发明的权利要求保护范围之内。
Claims (8)
1.一种阵列天线方向图形成的快速方法,其特征在于,所述方法通过将阵列天线方向图分解为一系列不同方向、相位、强度的照射平面波的集合,计算出该平面波集合照射在阵列口径上的复场分布,该结果即为所需方向图对应的阵列天线单元幅度、相位加权系数。
2.根据权利要求1所述的方法,其特征在于,所述方法的算法流程如下:
步骤一:对阵列天线方向图进行分解,获得照射波束中心指向以及基本形状方向图;
步骤二:对基本形状方向图进行分解,获得照射平面波集合;
步骤三:根据照射平面波集合,计算阵列口径处的基本形状方向图复场分布;
步骤四:根据照射波束中心指向,计算阵列口径处的扫描方向图复场分布;
步骤五:计算基本形状方向图复场分布与扫描方向图复场分布的乘积;
步骤六:根据基本形状方向图复场分布与扫描方向图复场分布的乘积结果,输出阵列天线单元幅度、相位加权系数。
4.根据权利要求2所述的方法,其特征在于,步骤二中对基本形状方向图进行分解,获得照射平面波集合,其具体步骤如下:
8.权利要求1-7任一项所述的方法在光学成像、微波成像、雷达探测、声呐、超声成像以及声、光、电目标探测和无线通信领域中的应用。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111142299.0A CN113806685B (zh) | 2021-09-28 | 2021-09-28 | 阵列天线方向图形成的快速方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111142299.0A CN113806685B (zh) | 2021-09-28 | 2021-09-28 | 阵列天线方向图形成的快速方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113806685A true CN113806685A (zh) | 2021-12-17 |
CN113806685B CN113806685B (zh) | 2024-03-26 |
Family
ID=78938710
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111142299.0A Active CN113806685B (zh) | 2021-09-28 | 2021-09-28 | 阵列天线方向图形成的快速方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113806685B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5081463A (en) * | 1989-04-13 | 1992-01-14 | Mitsubishi Denki Kabushiki Kaisha | Method and system for forming desired radiation pattern with array antenna |
US20090201214A1 (en) * | 2008-02-07 | 2009-08-13 | Saab Abse-581 | Wideband antenna pattern |
US20110025560A1 (en) * | 2006-09-22 | 2011-02-03 | Maurizio Crozzoli | Method and system for syntesizing array antennas |
CN108197362A (zh) * | 2017-12-23 | 2018-06-22 | 中国人民解放军战略支援部队信息工程大学 | Victs天线方向图和波束指向快速计算方法 |
CN112612024A (zh) * | 2020-12-13 | 2021-04-06 | 张艺恒 | 微波阵列快速成像方法 |
CN113361069A (zh) * | 2021-04-28 | 2021-09-07 | 西安电子科技大学 | 阵列天线方向图合成方法、系统、天线设计及测量方法 |
-
2021
- 2021-09-28 CN CN202111142299.0A patent/CN113806685B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5081463A (en) * | 1989-04-13 | 1992-01-14 | Mitsubishi Denki Kabushiki Kaisha | Method and system for forming desired radiation pattern with array antenna |
US20110025560A1 (en) * | 2006-09-22 | 2011-02-03 | Maurizio Crozzoli | Method and system for syntesizing array antennas |
US20090201214A1 (en) * | 2008-02-07 | 2009-08-13 | Saab Abse-581 | Wideband antenna pattern |
CN108197362A (zh) * | 2017-12-23 | 2018-06-22 | 中国人民解放军战略支援部队信息工程大学 | Victs天线方向图和波束指向快速计算方法 |
CN112612024A (zh) * | 2020-12-13 | 2021-04-06 | 张艺恒 | 微波阵列快速成像方法 |
CN113361069A (zh) * | 2021-04-28 | 2021-09-07 | 西安电子科技大学 | 阵列天线方向图合成方法、系统、天线设计及测量方法 |
Non-Patent Citations (4)
Title |
---|
FERNANDO RODRÍGUEZ VARELA ET AL: "Application of Nonuniform FFT to Spherical Near-Field Antenna Measurements", IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, vol. 68, no. 11, pages 7571 - 7579, XP011819317, DOI: 10.1109/TAP.2020.2998898 * |
徐秀成;刘强;: "基于FFT变换的阵列天线低副瓣方向图综合", 无线电通信技术, no. 06, pages 38 - 40 * |
李建新等: "基于FFT的阵列方向图快速计算", 微波学报, vol. 23, no. 1, pages 10 - 15 * |
陈宏超: "基于矩阵束算法的阵列方向图综合研究", 中国优秀硕士学位论文全文数据库 信息科技辑, no. 4 * |
Also Published As
Publication number | Publication date |
---|---|
CN113806685B (zh) | 2024-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112180329B (zh) | 一种基于阵元随机均匀分布球阵反卷积波束形成的汽车噪声源声成像方法 | |
AU2015292238B2 (en) | Planar sensor array | |
CN107884741B (zh) | 一种多球阵列多宽带声源快速定向方法 | |
CN113820711B (zh) | 阵列快速成像方法及其应用 | |
CN107515956B (zh) | 一种基于hfss单元法的大型有限平面阵列分析方法 | |
Xenaki et al. | Improving the resolution of three-dimensional acoustic imaging with planar phased arrays | |
CN109507643B (zh) | 旁瓣相消中展宽零陷的方法 | |
CN112162266B (zh) | 一种基于凸优化理论的共形阵二维波束优化方法 | |
CN114527427B (zh) | 一种基于球形麦克风阵列的低频波束形成声源定位方法 | |
CN108761394A (zh) | 一种基于时空处理的高分辨低旁瓣去卷积谱估计方法 | |
CN115038012B (zh) | 基于admm的麦克风阵列鲁棒频率不变波束形成方法 | |
CN109375227A (zh) | 一种解卷积波束形成三维声成像方法 | |
Chi et al. | High-resolution real-time underwater 3-D acoustical imaging through designing ultralarge ultrasparse ultra-wideband 2-D arrays | |
CN109541526B (zh) | 一种利用矩阵变换的圆环阵方位估计方法 | |
CN110554383B (zh) | 用于微波频段的mimo环形阵列方位向成像方法及装置 | |
CN113806685B (zh) | 阵列天线方向图形成的快速方法 | |
CN109669172A (zh) | 基于主瓣内强干扰抑制的弱目标方位估计方法 | |
CN1773307A (zh) | 小型天线阵口径扩展与空间信号处理方法 | |
Liu et al. | Beamforming of joint polarization-space matched filtering for conformal array | |
CN111551943A (zh) | 一种基于damas2的稀疏布阵高分辨三维声成像方法及系统 | |
CN117572435B (zh) | 一种基于反卷积的多波束合成孔径声呐高分辨成像方法 | |
Jiang et al. | Beam-time delay domain deconvolved scheme for high-resolution active localization of underwater targets | |
Huang | Theoretical and numerical simulation to control the wideband interfering noise by the broadband constant beam pattern method for a cylindrical array | |
Xie et al. | A deconvolved method for suppressing the sidelobes in the DOA estimation | |
Wang et al. | A Survey of Target Orientation Detection Algorithms Based on GPU Parallel Computing |
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 |