CN109670228A - 一种阵列天线单元方向图快速插值计算方法 - Google Patents
一种阵列天线单元方向图快速插值计算方法 Download PDFInfo
- Publication number
- CN109670228A CN109670228A CN201811515015.6A CN201811515015A CN109670228A CN 109670228 A CN109670228 A CN 109670228A CN 201811515015 A CN201811515015 A CN 201811515015A CN 109670228 A CN109670228 A CN 109670228A
- Authority
- CN
- China
- Prior art keywords
- data
- directional diagram
- coordinate
- interpolation
- 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 30
- 238000004364 calculation method Methods 0.000 title claims abstract description 8
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000004422 calculation algorithm Methods 0.000 abstract description 4
- 230000000694 effects Effects 0.000 abstract description 3
- 238000006243 chemical reaction Methods 0.000 abstract 1
- 238000000034 method Methods 0.000 description 11
- 238000013507 mapping Methods 0.000 description 2
- 230000007935 neutral effect Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Discrete Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Computing Systems (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公布了一种阵列天线单元方向图快速插值计算方法,针对传统单元方向图插值计算时间长和算法复杂的缺点,通过获取单元方向图投影数据再进行离散插值以达到同样效果。本发明不仅能够简便、快速转化天线单元方向图用于方向图乘积定理计算,而且通用性好,不限于单元方向图还适用于任意方向图转化为与坐标无关的数据。
Description
技术领域
本发明属于雷达天线技术领域,涉及平面相控阵天线单元方向图插值计算方法。
背景技术
现代雷达技术的发展对天线各项指标的要求越来越高,相控阵天线体制对于满足这些要求提供了巨大的潜力,它已经广泛应用于各种星载、机载、舰载和地面等雷达领域。天线方向图可以描述为二重级数求和的形式,对于电尺寸较大(几十、几百个波长甚至更大)的相控阵天线(如天基雷达大口径相控阵、大型稀疏阵天线等),直接根据公式按级数求和计算相控阵天线方向图的方法计算速度慢、效率低,不能适应工程设计过程中大量仿真计算需求。自从傅立叶在1822年提出傅立叶级数以来,傅立叶变换得到了很大的发展。1965年Cooley-Tukey提出了快速傅立叶变换(FFT)算法,很大程度上提高了傅立叶变换的计算速度,在许多科学领域获得了广泛和成功的应用。例如在信号处理中的应用,在平面近场测量中为了加速天线参数评估,采用FFT技术提高了效率。
平面阵列可由(M+1)×(N+1)维矩阵表征相控阵天线在U-V空间的二维方向图,取值点落在正方形区域内,与圆形可见区域并不重合。由于阵因子方向图具有周期性特征,从而可将数据拓展到整个U-V平面内。针对不同的单元间距,FFT变换后的数据点需要进行补点或删点。取整函数定义为g(x)=sgn(x)·int(|x|),由此确定阵因子数据的两个维度取值为[2g(Mdx/λ),2g(Ndy/λ)]。随着取值总点数的成倍增加,级数求和公式计算方向图的时间成倍增加。用FFT计算的时间以点数增加倍数的乘方增加。随着单元数的增加,在相同取值总点数的情况下,用级数求和公式计算阵列天线方向图的时间大致按相同倍率增加。但采用FFT计算方向图时,计算时间只随着变换点数的增加而增加,与单元数无关。由此可以看出,利用FFT计算的速度要比求和算法快得多,尤其是单元数越多,差别越明显。
现有商用仿真软件中得到天线阵中单元方向图数据分布在U-V坐标系下是不均匀的,无法直接使用。因此需要利用离散数据进行插值,才能与经FFT逆变换得到的阵因子相乘。工程应用中离散数据插值通常使用反距离加权法,具体为:在计算插值点取值时按距离越近权重越大的原则,用若干临近离散点数据的线性加权来拟合点的值。表达式为式中,Z即是插值点的值,Zi(i=1,2,…,t)是选取的离散点的值,di是插值点与离散点的距离,p为距离的幂。根据反距离加权法,选取p=2,即采用距离平方反比法,同时选择t=4,即采用四点插值法。由此可见,该插值算法的计算效率并不高,特别不利于进行多频点批量计算。本发明通过获取单元方向图投影数据后再进行离散插值以达到同样效果,计算效率改善显著,且不限于单元方向图还适用于任意方向图转化为与坐标无关的数据。
发明内容
针对现有的阵列天线单元方向图插值技术缺点,本发明为提高方向图插值效率,提出一种快速的平面阵列天线单元方向图插值方法,具体包含以下步骤:
第一步:已知U-V坐标下单元方向图数据A0(U,V)以及对应的U、V坐标数据;
第二步:将U、V坐标数据转化为整数,以满足最小整数为1,最大整数为方向图数据规模的最小值;
第三步:将A0(U,V)对应第二步中取整的U、V坐标顺序重新排列,转化为U-V平面上的投影数据A1,该矩阵与U、V坐标无关;
第四步:对投影数据A1先去除零元素,再进行离散插值,再将空元素补零得到A2;
第五步:对A2按照快速傅里叶变换计算的阵因子数据规模进行插值,得到可用于方向图乘积计算的单元方向图数据A3。
本发明与现有技术方法相比,其有益效果是:
①针对阵列天线单元方向图插值耗时的缺点,本发明通过获取单元方向图投影数据后再进行离散插值达到同样效果,提高了插值计算的效率。
②本发明通用性好,本发明不限于单元方向图还适用于任意方向图转化为与坐标无关的数据。
下面结合附图对本发明作进一步详细描述。
附图说明
图1为阵列天线单元方向图快速插值方法流程图。
图2为Vivaldi天线阵中单元方向图在U-V坐标系下三维图。
图3为U、V坐标数据变换示意图。
图4为Vivaldi天线阵中单元方向图的投影图。
图5为Vivaldi天线阵中单元方向图的投影插值图。
图6为Vivaldi天线阵中单元方向图的投影加密插值图。
具体实施方式
本实施例以CST仿真导出的8×8规模Vivaldi天线阵中单元方向图数据为例进行快速插值计算方法的说明,以下是对方向图数据快速插值方法的步骤。
第一步:一般仿真软件是在球坐标系下建模,因此导出的方向图数据按照公式变换成U-V坐标系下的方向图A0(U,V)。本例中的已知数据规模均为180×180的矩阵,由于A0(U,V)数据和U、V坐标有关,所以图2是采用Matlab中作图函数meshc(u,v,A0)画出的单元方向图形状,可以看出方向图在U-V平面的投影也包含了所需方向图数据信息;
第二步:将U、V坐标数据转化为整数,以满足最小整数为1,最大整数为方向图数据规模的最小值,本例中为180,如图3所示;
第三步:利用第二步中整数化的U、V坐标数据,将A0(U,V)转化为U-V平面上的投影数据A1,该数据与U、V坐标无关,所以图4是采用Matlab中作图函数mesh(A1)画出的单元方向图形状,可以看出该方向图中有很多零元素;
第四步:将投影数据A1先去除零元素,再进行离散插值,再将空元素补零得到A2,如图5所示;
第五步:对A2按照快速傅里叶变换计算的阵因子数据规模进行插值,得到可用于方向图乘积定理计算的单元方向图数据A3。本例中阵因子数据规模为485×1025,经过该方法计算后的单元方向图如图6所示。
Claims (1)
1.一种阵列天线单元方向图快速插值计算方法,其特征在于包含以下步骤:
第一步:已知U-V坐标下单元方向图数据A0(U,V)以及对应的U、V坐标数据;
第二步:将U、V坐标数据转化为整数,以满足最小整数为1,最大整数为方向图数据规模的最小值;
第三步:将A0(U,V)对应第二步中取整的U、V坐标顺序重新排列,转化为U-V平面上的投影数据A1,该矩阵与U、V坐标无关;
第四步:对投影数据A1先去除零元素,再进行离散插值,再将空元素补零得到A2;
第五步:对A2按照快速傅里叶变换计算的阵因子数据规模进行插值,得到可用于方向图乘积计算的单元方向图数据A3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811515015.6A CN109670228B (zh) | 2018-12-12 | 2018-12-12 | 一种阵列天线单元方向图快速插值计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811515015.6A CN109670228B (zh) | 2018-12-12 | 2018-12-12 | 一种阵列天线单元方向图快速插值计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109670228A true CN109670228A (zh) | 2019-04-23 |
CN109670228B CN109670228B (zh) | 2023-12-15 |
Family
ID=66143737
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811515015.6A Active CN109670228B (zh) | 2018-12-12 | 2018-12-12 | 一种阵列天线单元方向图快速插值计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109670228B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112348793A (zh) * | 2020-11-05 | 2021-02-09 | 中国电子科技集团公司第二十九研究所 | 一种自动识别并计算天线方向图凹坑的方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104391183A (zh) * | 2014-10-11 | 2015-03-04 | 中国船舶重工集团公司第七二四研究所 | 近场测量天线远区场特性快速计算方法 |
-
2018
- 2018-12-12 CN CN201811515015.6A patent/CN109670228B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104391183A (zh) * | 2014-10-11 | 2015-03-04 | 中国船舶重工集团公司第七二四研究所 | 近场测量天线远区场特性快速计算方法 |
Non-Patent Citations (5)
Title |
---|
刘恒,赵宏伟,李维梅,刘波: "平面六边形稀疏矩阵天线的优化", vol. 55, no. 12, pages 1365 - 1370 * |
刘恒;刘波;谢广钱;赵宏伟;: "迭代傅里叶算法用于六边形稀疏阵列天线", vol. 30, no. 06 * |
刘明罡: "平面天线阵方向图计算中的仿射坐标 FFT 方法", 《微波学报》, vol. 27, no. 2, pages 8 - 12 * |
刘明罡;: "平面天线阵方向图计算中的仿射坐标FFT方法", no. 02 * |
郜参观;周希朗;蒋凡杰;高为;: "相控微带阵馈电的抛物柱天线方向图快速计算", no. 09 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112348793A (zh) * | 2020-11-05 | 2021-02-09 | 中国电子科技集团公司第二十九研究所 | 一种自动识别并计算天线方向图凹坑的方法 |
CN112348793B (zh) * | 2020-11-05 | 2023-04-25 | 中国电子科技集团公司第二十九研究所 | 一种自动识别并计算天线方向图凹坑的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109670228B (zh) | 2023-12-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108872929B (zh) | 基于内插虚拟阵列协方差矩阵子空间旋转不变性的互质阵列波达方向估计方法 | |
CN112989680B (zh) | 减少网格使用量的fvfd远场积分边界条件计算方法 | |
CN108896954B (zh) | 互质阵中一种基于联合实值子空间的波达角估计方法 | |
CN102156764B (zh) | 一种分析天线辐射和电磁散射的多分辨预条件方法 | |
CN111610486B (zh) | 基于平面互质阵列虚拟域张量空间谱搜索的高分辨精确二维波达方向估计方法 | |
CN102081690B (zh) | 复杂电路的矩阵分解结合新奇异值分解方法 | |
JP7327840B2 (ja) | 相互相関テンソルに基づく三次元の互いに素のキュービックアレイの到来方向推定方法 | |
CN104391183A (zh) | 近场测量天线远区场特性快速计算方法 | |
CN106199524B (zh) | 基于基追踪去噪的远场宽带rcs数据采集与压缩方法 | |
CN110146842B (zh) | 基于欠采样的信号载频与二维doa参数估计方法 | |
CN109670228A (zh) | 一种阵列天线单元方向图快速插值计算方法 | |
CN108614234B (zh) | 基于多采样快拍互质阵列接收信号快速傅里叶逆变换的波达方向估计方法 | |
Chen et al. | Analysis of electromagnetic scattering from PEC targets using improved fast dipole method | |
Bruno et al. | A fast high-order solver for EM scattering from complex penetrable bodies: TE case | |
Caputo et al. | Neural network characterization of microstrip patches for reflectarray optimization | |
Zhang et al. | Combination of asymptotic phase basis functions and matrix interpolation method for fast analysis of monostatic RCS | |
Blankrot et al. | FMM-accelerated source-model technique for many-scatterer problems | |
CN106503349A (zh) | 一种类周期结构目标电磁散射特性快速计算方法 | |
Meng et al. | Target localization for monostatic MIMO radar based on a double polynomial rooting technique | |
Song et al. | On the investigation of triangular interpolation-based multilevel fast multipole algorithm | |
Ergul et al. | Fast and accurate solutions of large-scale scattering problems with parallel multilevel fast multipole algorithm | |
CN107622034B (zh) | 基于3d-fft的用于快速计算大型机会阵方向图的方法 | |
Laviada et al. | Efficient calculation of the reduced matrix in the Characteristic Basis Functions Method | |
Sun et al. | A novel method for fast analysis of electromagnetic scattering characteristics of pec objects | |
RU213021U9 (ru) | Устройство определения направления прихода радиосигнала |
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 | ||
CB02 | Change of applicant information |
Address after: 210003 No. 346, Zhongshan North Road, Jiangsu, Nanjing Applicant after: 724 Research Institute of China Shipbuilding Corp. Address before: 210003 No. 346, Zhongshan North Road, Jiangsu, Nanjing Applicant before: 724TH RESEARCH INSTITUTE OF CHINA SHIPBUILDING INDUSTRY Corp. |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |