CN108446504B - 基于凸优化的近场阵列天线方向图综合方法 - Google Patents
基于凸优化的近场阵列天线方向图综合方法 Download PDFInfo
- Publication number
- CN108446504B CN108446504B CN201810241246.6A CN201810241246A CN108446504B CN 108446504 B CN108446504 B CN 108446504B CN 201810241246 A CN201810241246 A CN 201810241246A CN 108446504 B CN108446504 B CN 108446504B
- Authority
- CN
- China
- Prior art keywords
- field
- array
- side lobe
- focusing
- convex optimization
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
- Aerials With Secondary Devices (AREA)
Abstract
本发明提供一种基于凸优化的近场阵列天线方向图综合方法,包括以下步骤:步骤1.确定主瓣形状、需进行副瓣控制的位置范围和峰值副瓣电平;步骤2.基于点聚焦平面天线阵列建立辐射近场区的辐射模型;步骤3.建立凸优化模型并求解得到天线阵列的激励;本发明具有能生成辐射近场区成任意形状、副瓣可控的方向图和计算速度快的优点;可用于近场平面阵列天线的优化设计。
Description
技术领域
本发明属于天线技术领域,具体是一种基于凸优化的近场阵列天线方向图综合方法。
背景技术
不同于一般工作在辐射远场区域的天线,对于某些特定场合,要求天线辐射的能量聚集在其近场区域的某一特定点或区域上,这就需要近场聚焦天线。针对近场天线的研究已经持续了很长时间。在特定区域具有高能量密度的近场聚焦天线不仅可以用于提高医学与工业检测系统的探测和成像能力,还能用于门禁系统、微波无线输能系统、射频识别系统以及微波无线电系统等方面,具有广泛的应用前景和研究价值。
跟成熟的远场方向图综合方法相比,近场去方向图综合方法较少,且大部分都集中在对点聚焦的研究上,如最初的直接把远场的低副瓣综合方法应用在近场,但这种办法生成的副瓣电平不可控制;又如对整个阵列进行优化的最小均方差方法,但其计算量巨大,对于阵元数一般较大的近场阵列天线来讲,针对所有阵元激励的优化必然增大计算量,同时该方法也无法对方向图主瓣的形状进行控制。
发明内容
本发明针对上述现有的近场阵列天线方向图综合方法存在的难以控制波束形状、副瓣电平较高且难以控制、计算时间长的问题,提出一种基于凸优化的在辐射近场区可生成任意波束形状、具有可控的副瓣电平的高效算法。
为实现上述目的,本发明技术方案如下:
一种基于凸优化的近场阵列天线方向图综合方法,包括以下步骤:
步骤1.确定主瓣形状、需进行副瓣控制的位置范围和峰值副瓣电平;
步骤2.基于点聚焦平面天线阵列建立辐射近场区的辐射模型;
步骤3.建立凸优化模型并求解得到天线阵列的激励。
作为优选方式,步骤2进一步为:
步骤2.基于点聚焦平面天线阵列建立辐射近场区的辐射模型;
近场聚焦平面阵列采用均匀布阵形式,共有N个阵元,其焦点与平面天线的距离为z,则近场电场强度被表示为:
其中r是近场观测点位置,rf是聚焦点位置,rn为第n个阵元的位置,an是第n个单元的复激励,其相位值由焦点位置决定,而幅度以需求的副瓣电平对应的远场切比雪夫或者泰勒低副瓣锥削作为初值;λ是自由空间中的工作波长,Fn(r)是第n个单元的归一化辐射方向图。
作为优选方式,步骤3进一步为:
步骤3.建立凸优化模型并求解得到天线阵列的激励
近场区域任意形状方向图的形成看作由M个不同聚焦位置的近场聚焦波束叠加生成,其焦点位置为(rf1,rf2,…,rfM),权重为(b1,b2,…,bM),M通过待生成方向图的主瓣宽度与点聚焦的3dB波束宽度的关系得到,具体的有
其中2BSx0.5与2BSy0.5是待生成方向图x、y维度的3dB波束宽度,2By0.5与2By0.5则是点聚焦的3dB波束宽度,此时有M=Mx×My;值得提出的是,聚焦点位置不同的点聚焦方向图的波束宽度实际上不相等但差别不大,本发明将其当做相等来处理;
此时待生成的任意形状近场方向图被表示为:
对方向图所在平面进行离散,假设离散后共有K个观测点rok,k=1,…,K,此时Enf写为以下的矩阵形式:
Enf=uAP
其中
u=(b1,b2,...,bM)
将Enf归一化后,凸优化问题表示为:
其中||·||2表示取二范数,Enf_ml和Enf_sl分别为Enf的主瓣和副瓣部分,Mml是待生成方向图的主瓣约束,Msl是待生成方向图的副瓣约束,主瓣部分采用均匀采样,其观测点数为K1,针对副瓣部分,为简化说明,以一维为例,在x<0与x>0的副瓣部分采样方式如下
其中x_negt与x_posp分布表示x<0与x>0副瓣区域观测点的位置,t=1,…,T,p=1,…,P,sl_neg与sl_pos分别表示为近场x<0与x>0副瓣区域的开始位置,副瓣区域观测点数为K2=T+P,且有K=K1+K2。
采用凸优化求解器对该问题进行求解得到u,阵列的激励I则通过I=uA得到。
本发明的有益效果为:本发明具有能生成辐射近场区成任意形状、副瓣可控的方向图和计算速度快的优点;可用于近场平面阵列天线的优化设计。
附图说明
图1为本发明的流程图。
图2为用本发明设计得到的阵列天线实例的辐射近场区方向图。
具体实施方式
以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。
一种基于凸优化的近场阵列天线方向图综合方法,包括以下步骤:
步骤1.确定主瓣形状、需进行副瓣控制的位置范围和峰值副瓣电平;
步骤2.基于点聚焦平面天线阵列建立辐射近场区的辐射模型;
近场聚焦平面阵列采用均匀布阵形式,共有N个阵元,其焦点与平面天线的距离为z,则近场电场强度被表示为:
其中r是近场观测点位置,rf是聚焦点位置,rn为第n个阵元的位置,an是第n个单元的复激励,其相位值由焦点位置决定,而幅度以需求的副瓣电平对应的远场切比雪夫或者泰勒低副瓣锥削作为初值;λ是自由空间中的工作波长,Fn(r)是第n个单元的归一化辐射方向图。
步骤3.建立凸优化模型并求解得到天线阵列的激励。
近场区域任意形状方向图的形成看作由M个不同聚焦位置的近场聚焦波束叠加生成,其焦点位置为(rf1,rf2,…,rfM),权重为(b1,b2,…,bM),M通过待生成方向图的主瓣宽度与点聚焦的3dB波束宽度的关系得到,具体的有
其中2BSx0.5与2BSy0.5是待生成方向图x、y维度的3dB波束宽度,2By0.5与2By0.5则是点聚焦的3dB波束宽度,此时有M=Mx×My;值得提出的是,聚焦点位置不同的点聚焦方向图的波束宽度实际上不相等但差别不大,本发明将其当做相等来处理;
此时待生成的任意形状近场方向图被表示为:
对方向图所在平面进行离散,假设离散后共有K个观测点rok,k=1,…,K,此时Enf写为以下的矩阵形式:
Enf=uAP
其中
u=(b1,b2,...,bM)
将Enf归一化后,凸优化问题表示为:
其中||·||2表示取二范数,Enf_ml和Enf_sl分别为Enf的主瓣和副瓣部分,Mml是待生成方向图的主瓣约束,Msl是待生成方向图的副瓣约束。本实例的综合区域为二维平面,主瓣部分采用均匀采样,其观测点数为K1,副瓣部分采样方式由上述的一维采样方式推广过来,采用如下matlab代码采样
其中sl为近场副瓣区域的开始位置,副瓣区域观测点数为K2,x_prime为副瓣观测点的位置,且有K=K1+K2×K2。
采用凸优化求解器对该问题进行求解得到u,阵列的激励I则通过I=uA得到。
本发明的效果通过以下仿真实验进一步说明:
1.设置仿真参数:本实例主瓣形状为主瓣宽度4λ×4λ平头波束,主瓣约束为Mml=1,其范围为|x|,|y|<2λ,副瓣约束为Msl=0.1,其范围为|x|,|y|≥2.8λ,即sl=2.8λ,K1=2601,K2=102,K=13005。采用的天线口径为10λ×10λ,单元间本实例采用口径为10λ×10λ,单元间距为0.5λ的等间距布阵,阵元总数为N=441,聚焦平面与天线口径的距离为10λ,an的幅度采用20dB远场泰勒低副瓣锥削,M=25。
2.仿真内容
根据设置,按照步骤3建立近场聚焦模型,再根据设计要求和近场聚焦模型,建立包括评价函数和约束条件的凸优化模型,再采用凸优化求解器进行求解得到阵列的激励I。使用阵列激励I计算阵列天线在距离天线口径10λ,大小为10λ×10λ的平面上的电场强度,将其归一化和取对数后画成图2。
由图2可知,得到近场阵列天线的方向图的参数与设计要求一致,主瓣宽度为4λ×4λ,峰值副瓣电平为21.89dB。在CPU为i5-3210M,内存8G的计算机上的计算时间仅为17.4秒。
上述实施例仅例示性说明本发明的原理及其功效,而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下,对上述实施例进行修饰或改变。因此,凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变,仍应由本发明的权利要求所涵盖。
Claims (1)
1.一种基于凸优化的近场阵列天线方向图综合方法,其特征在于包括以下步骤:
步骤1.确定主瓣形状、需进行副瓣控制的位置范围和峰值副瓣电平;
步骤2.基于点聚焦平面天线阵列建立辐射近场区的辐射模型;
近场聚焦平面阵列采用均匀布阵形式,共有N个阵元,其焦点与平面天线的距离为z,则近场电场强度被表示为:
其中r是近场观测点位置,rf是聚焦点位置,rn为第n个阵元的位置,an是第n个单元的复激励,其相位值由焦点位置决定,而幅度以需求的副瓣电平对应的远场切比雪夫或者泰勒低副瓣锥削作为初值;λ是自由空间中的工作波长,Fn(r)是第n个单元的归一化辐射方向图;
步骤3.建立凸优化模型并求解得到天线阵列的激励
近场区域任意形状方向图的形成看作由M个不同聚焦位置的近场聚焦波束叠加生成,其焦点位置为(rf1,rf2,…,rfM),权重为(b1,b2,…,bM),M通过待生成方向图的主瓣宽度与点聚焦的3dB波束宽度的关系得到,具体的有
其中2BSx0.5与2BSy0.5是待生成方向图x、y维度的3dB波束宽度,2By0.5与2By0.5则是点聚焦的3dB波束宽度,此时有M=Mx×My;值得提出的是,聚焦点位置不同的点聚焦方向图的波束宽度实际上不相等但差别不大,本方法将其当做相等来处理;
此时待生成的任意形状近场方向图被表示为:
对方向图所在平面进行离散,假设离散后共有K个观测点rok,k=1,…,K,此时Enf写为以下的矩阵形式:
Enf=uAP
其中
u=(b1,b2,...,bM)
将Enf归一化后,凸优化问题表示为:
其中||·||2表示取二范数,Enf_ml和Enf_sl分别为Enf的主瓣和副瓣部分,Mml是待生成方向图的主瓣约束,Msl是待生成方向图的副瓣约束,主瓣部分采用均匀采样,其观测点数为K1,针对副瓣部分,为简化说明,以一维为例,在x<0与x>0的副瓣部分采样方式如下
其中x_negt与x_posp分布表示x<0与x>0副瓣区域观测点的位置,t=1,…,T,p=1,…,P,sl_neg与sl_pos分别表示为近场x<0与x>0副瓣区域的开始位置,副瓣区域观测点数为K2=T+P,且有K=K1+K2;
采用凸优化求解器对该问题进行求解得到u,阵列的激励I则通过I=uA得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810241246.6A CN108446504B (zh) | 2018-03-22 | 2018-03-22 | 基于凸优化的近场阵列天线方向图综合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810241246.6A CN108446504B (zh) | 2018-03-22 | 2018-03-22 | 基于凸优化的近场阵列天线方向图综合方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108446504A CN108446504A (zh) | 2018-08-24 |
CN108446504B true CN108446504B (zh) | 2021-09-07 |
Family
ID=63196123
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810241246.6A Active CN108446504B (zh) | 2018-03-22 | 2018-03-22 | 基于凸优化的近场阵列天线方向图综合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108446504B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109541643B (zh) * | 2018-11-09 | 2023-02-03 | 电子科技大学 | 一种阵列天线的副瓣和交叉极化抑制方法 |
CN110717243B (zh) * | 2019-08-28 | 2021-05-14 | 西安电子科技大学 | 一种基于线性约束的宽带方向图综合方法 |
CN112615158B (zh) * | 2020-12-01 | 2022-01-28 | 厦门大学 | 超宽带扫描稀疏阵列天线的综合方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103825089A (zh) * | 2014-02-28 | 2014-05-28 | 电子科技大学 | 近场聚焦平面阵列天线 |
CN106898884A (zh) * | 2017-02-23 | 2017-06-27 | 金陵科技学院 | 基于混合算法的平面阵列天线和差波束方向图合成方法 |
CN107490729A (zh) * | 2017-08-18 | 2017-12-19 | 北京航空航天大学 | 一种天线近场无相位测量方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2378444B1 (de) * | 2010-04-13 | 2015-07-29 | CST-Computer Simulation Technology AG | Verfahren, Vorrichtung und Computerprogrammprodukt zur Bestimmung eines elektromagnetischen Nahfeldes einer Feldanregungsquelle eines elektrischen Systems |
-
2018
- 2018-03-22 CN CN201810241246.6A patent/CN108446504B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103825089A (zh) * | 2014-02-28 | 2014-05-28 | 电子科技大学 | 近场聚焦平面阵列天线 |
CN106898884A (zh) * | 2017-02-23 | 2017-06-27 | 金陵科技学院 | 基于混合算法的平面阵列天线和差波束方向图合成方法 |
CN107490729A (zh) * | 2017-08-18 | 2017-12-19 | 北京航空航天大学 | 一种天线近场无相位测量方法 |
Non-Patent Citations (2)
Title |
---|
Near-field beamformer design of non-uniform circular antenna array based on convex optimization;Yanping Liao 等;《ICIC Express Letters, Part B: Applications》;20150831;全文 * |
共形阵的优化设计及故障诊断算法研究;龙政斌;《中国优秀硕士学位论文全文数据库信息科技辑(月刊 )》;20180215;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN108446504A (zh) | 2018-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106886656B (zh) | 一种三维阵列天线方向图旁瓣抑制方法 | |
Panduro et al. | A comparison of genetic algorithms, particle swarm optimization and the differential evolution method for the design of scannable circular antenna arrays | |
Gies et al. | Particle swarm optimization for reconfigurable phase‐differentiated array design | |
Dib et al. | Application of Taguchi's optimization method and self-adaptive differential evolution to the synthesis of linear antenna arrays | |
Bucci et al. | Deterministic synthesis of uniform amplitude sparse arrays via new density taper techniques | |
CN108446504B (zh) | 基于凸优化的近场阵列天线方向图综合方法 | |
CN107729640A (zh) | 一种使用最少阵元的稀疏天线阵列综合布阵方法 | |
CN109033647B (zh) | 基于l1范数约束的近场稀布天线阵列优化方法 | |
CN106446469B (zh) | 基于二阶近似计算公式的空间网状天线电性能动力响应分析方法 | |
CN104112901B (zh) | 全息人工阻抗表面共形天线 | |
Li et al. | Improved GA and PSO culled hybrid algorithm for antenna array pattern synthesis | |
Koziel et al. | Rapid multi-objective optimization of antennas using nested kriging surrogates and single-fidelity EM simulation models | |
CN108446503B (zh) | 基于压缩感知与凸优化的近场稀布天线阵列优化方法 | |
CN105160115B (zh) | 基于近似与灵敏度分析的反射面天线机电集成优化设计方法 | |
Alvarez et al. | Design of antenna arrays for near-field focusing requirements using optimisation | |
CN104485513A (zh) | 宽频带天线阵列实时综合方法 | |
CN106021766B (zh) | 基于机电集成优化的圆极化网状天线高指向精度设计方法 | |
CN109271735B (zh) | 基于量子启发式引力搜索算法的阵列方向图综合方法 | |
CN106842951B (zh) | 面向电性能与控制的空间网状天线状态空间建模方法 | |
Zhang et al. | Multiple beamforming with null steering based on improved invasive weed optimization | |
CN106934097B (zh) | 面向电性能的空间网状天线关键动力学模态选取方法 | |
Yang et al. | An accurate near-field focusing of array antenna based on near-field active element pattern and infinitesimal dipole modeling | |
Mallahzadeh | Cosecant squared pattern synthesis for reflector antenna using a stochastic method | |
Radivojević et al. | Radiation pattern optimisation of an antenna array on the spherical surface by using a varying number of optimisation parameters | |
Cong et al. | A pattern synthesis method for large planar antenna array |
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 |