CN114840052B - 一种通用化二维数字波束形成技术的ip核及生成方法 - Google Patents

一种通用化二维数字波束形成技术的ip核及生成方法 Download PDF

Info

Publication number
CN114840052B
CN114840052B CN202210782939.2A CN202210782939A CN114840052B CN 114840052 B CN114840052 B CN 114840052B CN 202210782939 A CN202210782939 A CN 202210782939A CN 114840052 B CN114840052 B CN 114840052B
Authority
CN
China
Prior art keywords
antenna
spacing
calculating
carrier frequency
vector
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
Application number
CN202210782939.2A
Other languages
English (en)
Other versions
CN114840052A (zh
Inventor
赵佳琪
冯浩轩
全英汇
肖国尧
孙宗正
张宇
刘成
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xidian University
Original Assignee
Xidian University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xidian University filed Critical Xidian University
Priority to CN202210782939.2A priority Critical patent/CN114840052B/zh
Publication of CN114840052A publication Critical patent/CN114840052A/zh
Application granted granted Critical
Publication of CN114840052B publication Critical patent/CN114840052B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F1/00—Details not covered by groups G06F3/00 - G06F13/00 and G06F21/00
    • G06F1/02—Digital function generators
    • G06F1/03—Digital function generators working, at least partly, by table look-up
    • G06F1/0321—Waveform generators, i.e. devices for generating periodical functions of time, e.g. direct digital synthesizers
    • 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
    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03L—AUTOMATIC CONTROL, STARTING, SYNCHRONISATION OR STABILISATION OF GENERATORS OF ELECTRONIC OSCILLATIONS OR PULSES
    • H03L7/00—Automatic control of frequency or phase; Synchronisation
    • H03L7/06—Automatic control of frequency or phase; Synchronisation using a reference signal applied to a frequency- or phase-locked loop
    • H03L7/16—Indirect frequency synthesis, i.e. generating a desired one of a number of predetermined frequencies using a frequency- or phase-locked loop
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

本发明涉及一种通用化二维数字波束形成技术的IP核及生成方法，生成方法包括：配置阵列天线的天线行数、天线列数、天线行间距、天线列间距和相位零点，并计算各阵元距离相位零点的水平间距和竖直间距；通过载频系数配置载频本振，并利用期望波束形成方向、载频本振、水平间距和竖直间距计算期望方向的空域导向矢量；利用天线行数、天线列数和高斯窗口宽度参数计算目标高斯窗；将空域导向矢量与目标高斯窗相乘，计算得到波束形成权系数；对上述步骤所形成的程序进行综合与时序优化，并生成可调用的通用化IP核。该IP核生成方法有效提高了该方法对于不同阵列天线排布的适用性，提高了该方法的通用性；该方法的可移植性高，操作简单，实现简单易懂。

Description

一种通用化二维数字波束形成技术的IP核及生成方法
技术领域
本发明属于信号处理技术领域，具体涉及一种通用化二维数字波束形成技术的IP核及生成方法。
背景技术
当多个不同位置的天线接收同一位置的信号时，信号到达天线的时候会出现相位差。数字波束形成(Digital Beam Forming，DBF)技术在阵列雷达领域应用十分广泛，是数字阵列雷达理论的重要组成部分。DBF技术的原理就是通过生成权系数将同一信号到达不同位置天线时将因为路程不同而导致的相位差进行补偿，让所有不同位置天线形成的波束都指向同一个方向，让其在期望方向的能量接收最大。DBF技术可以通过改变权值让阵列波束指向不同的方向，实现波束扫描，并且可以选用合适的窗函数来降低副瓣。
然而，DBF技术对于阵列天线的排布等因素的影响较大，所以在工程实现中一般都是通过FPGA的Verilog语言配合乘法器等IP核进行算法实现。但是该方法具有可移植性差、通用性低的缺点。
发明内容
为了解决现有技术中存在的上述问题，本发明提供了一种通用化二维数字波束形成技术的IP核及生成方法。本发明要解决的技术问题通过以下技术方案实现：
本发明实施例提供了一种通用化二维数字波束形成技术的IP核生成方法，包括步骤：
配置阵列天线的天线行数、天线列数、天线行间距、天线列间距和相位零点，并计算所述阵列天线中各阵元距离所述相位零点的水平间距和竖直间距；
通过载频系数配置载频本振，并利用期望波束形成方向、所述载频本振、所述水平间距和所述竖直间距计算期望方向的空域导向矢量；
利用所述天线行数、所述天线列数和高斯窗输出点数、高斯窗口宽度参数计算目标高斯窗；
将所述空域导向矢量与所述目标高斯窗相乘，计算得到波束形成权系数；
对上述步骤所形成的程序进行综合与时序优化，并将优化后的程序生成可调用的通用化IP核。
在本发明的一个实施例中，计算所述阵列天线中各阵元距离相位零点的水平间距和竖直间距，包括：
利用所述天线行数、所述天线列数、所述天线行间距、所述天线列间距计算各阵元距离所述相位零点的所述水平间距和所述竖直间距。
在本发明的一个实施例中，所述相位零点包括所述阵列天线的中心位置或者所述阵列天线的左下角天线。
在本发明的一个实施例中，当所述相位零点为所述阵列天线的中心位置时，所述水平间距和所述竖直间距为：
Figure 65215DEST_PATH_IMAGE001
Figure 429200DEST_PATH_IMAGE002
其中，
Figure 665140DEST_PATH_IMAGE003
表示各阵元距离相位零点的水平间距，
Figure 865309DEST_PATH_IMAGE004
表示各阵元距离相位零点的竖直间距，
Figure 297427DEST_PATH_IMAGE005
表示天线列数，
Figure 1072DEST_PATH_IMAGE006
表示天线行数，
Figure 860443DEST_PATH_IMAGE007
表示天线列间距，
Figure 344646DEST_PATH_IMAGE008
表示天线行间距。
在本发明的一个实施例中，当所述相位零点为所述阵列天线的左下角天线时，所述水平间距和所述竖直间距为：
Figure 331187DEST_PATH_IMAGE009
Figure 404185DEST_PATH_IMAGE010
其中，
Figure 185191DEST_PATH_IMAGE011
表示各阵元距离相位零点的水平间距，
Figure 343640DEST_PATH_IMAGE012
表示各阵元距离相位零点的竖直间距，
Figure 868293DEST_PATH_IMAGE005
表示天线列数，
Figure 795798DEST_PATH_IMAGE006
表示天线行数，
Figure 747704DEST_PATH_IMAGE007
表示天线列间距，
Figure 206499DEST_PATH_IMAGE008
表示天线行间距。
在本发明的一个实施例中，通过载频系数配置载频本振，并利用期望波束形成方向、所述载频本振、所述水平间距和所述竖直间距计算期望方向的空域导向矢量，包括：
通过所述载频系数配置所述载频本振：
Figure 518531DEST_PATH_IMAGE013
其中，
Figure 316854DEST_PATH_IMAGE014
表示射频本振，
Figure 688930DEST_PATH_IMAGE015
表示基底射频，
Figure 572703DEST_PATH_IMAGE016
表示载频系数，
Figure 501476DEST_PATH_IMAGE017
表示步进频率；
利用所述载频本振、所述水平间距、所述竖直间距和所述期望波束形成方向计算期望方向在方位向的空域导向矢量和期望方向在俯仰向的空域导向矢量：
Figure 403573DEST_PATH_IMAGE018
Figure 697282DEST_PATH_IMAGE019
其中，
Figure 52040DEST_PATH_IMAGE020
表示期望方向在方位向的空域导向矢量，
Figure 722187DEST_PATH_IMAGE021
表示期望方向在俯仰向的空域导向矢量，
Figure 280121DEST_PATH_IMAGE022
，
Figure 993999DEST_PATH_IMAGE023
，
Figure 852365DEST_PATH_IMAGE024
表示各阵元距离相位零点的水平间距，
Figure 309891DEST_PATH_IMAGE025
表示各阵元距离相位零点的竖直间距，
Figure 734050DEST_PATH_IMAGE026
表示利用载频本振计算出的发射波波长；
将所述方位向的空域导向矢量和所述俯仰向的空域导向矢量进行罗内克积，得到所述期望方向的空域导向矢量：
Figure 635141DEST_PATH_IMAGE027
其中，
Figure 964491DEST_PATH_IMAGE028
表示期望方向的空域导向矢量。
在本发明的一个实施例中，利用所述天线行数、所述天线列数和高斯窗口宽度参数计算目标高斯窗，包括：
将所述高斯窗输出点数分别取所述天线行数和所述天线列数，并利用所述高斯窗口宽度参数计算出第一高斯向量和第二高斯向量，计算公式为：
Figure 976441DEST_PATH_IMAGE029
其中，
Figure 442057DEST_PATH_IMAGE030
表示高斯窗点数，
Figure 514050DEST_PATH_IMAGE031
表示高斯窗口宽度参数，
Figure 143745DEST_PATH_IMAGE032
，
Figure 943074DEST_PATH_IMAGE033
取值
Figure 341825DEST_PATH_IMAGE030
次；
将所述第一高斯向量和所述第二高斯向量进行克罗内克积，计算出
Figure 584719DEST_PATH_IMAGE034
个高斯窗点数，得到所述目标高斯窗。
在本发明的一个实施例中，其特征在于，所述IP核生成方法利用Vivado-HLS工具采用C++语言实现。
本发明的另一个实施例提供了一种通用化二维数字波束形成的IP核，由上述任一实施例所述的通用化二维数字波束形成技术的IP核生成方法生成得到，包括：
天线阵列生成模块，用于配置阵列天线的天线行数、天线列数、天线行间距、天线列间距和相位零点，并计算所述阵列天线中各阵元距离所述相位零点的水平间距和竖直间距；
期望方向导向矢量生成模块，用于通过载频系数配置载频本振，并利用期望波束形成方向、所述载频本振、所述水平间距和所述竖直间距计算期望方向的空域导向矢量；
高斯窗生成模块，用于利用所述天线行数、所述天线列数和高斯窗口宽度参数计算目标高斯窗；
波束形成权系数生成模块，用于将所述空域导向矢量与所述目标高斯窗相乘，计算得到波束形成权系数；
IP核生成模块，对上述步骤所形成的程序进行综合与时序优化，并将优化后的程序生成可调用的通用化IP核。
与现有技术相比，本发明的有益效果：
1、本发明的IP核生成方法可以自定义相位零点，通过进行天线排布自定义并计算各阵元距离相位零点的距离可以实现阵列天线排布自定义，通过载频系数配置载频本振可以实现载频自定义，有效提高了该方法对于不同阵列天线排布的适用性，提高了该方法的通用性；结合天线行数、天线列数、高斯窗口宽度计算目标高斯窗可以实现窗函数自定义，有效提高了天线的综合性能，形成高增益的窄波束；另外，生成可调用的通用化IP核，可以提高该方法的可移植性，并且操作简单，方法实现简单易懂。
2、本发明的IP核生成方法利用Vivado-HLS工具采用C++语言实现，并对算法进行并行优化与时序优化，在保证计算精度的前提下，有效减小资源占用，提高了实时性要求。
附图说明
图1为本发明实施例提供的一种通用化二维数字波束形成技术的IP核生成方法的流程示意图；
图2为本发明实施例提供的一种通用化天线阵列排布图；
图3为本发明实施例提供的一种交叉天线阵列排布图。
具体实施方式
下面结合具体实施例对本发明做进一步详细的描述，但本发明的实施方式不限于此。
实施例一
请参见图1，图1为本发明实施例提供的一种通用化二维数字波束形成技术的IP核生成方法的流程示意图。该通用化二维数字波束形成技术的IP核生成方法利用Vivado HLS工具，采用C++语言编写将波束形成方法进行通用化设计，具体包括步骤：
S1、配置阵列天线的天线行数、天线列数、天线行间距、天线列间距和相位零点，并计算阵列天线中各阵元距离相位零点的水平间距和竖直间距。
请参见图2，图2为本发明实施例提供的一种通用化天线阵列排布图。
步骤S1具体包括步骤：
S11、获取阵列天线的配置参数，配置参数包括：天线行数、天线列数、天线行间距、天线列间距和相位零点。
具体的，相位零点可以为阵列天线中的任意位置，即相位零点可以为某个天线所在的位置上，也可以为无天线设置的位置上。优选的，相位零点包括阵列天线的中心位置或者阵列天线的左下角天线这两种模式。
S12、利用天线行数、天线列数、天线行间距、天线列间距计算各阵元距离相位零点的水平间距和竖直间距。
在一个具体实施例中，当相位零点为阵列天线的中心位置，即坐标为(M/2，N/2)的位置为相位零点时，水平间距和竖直间距为：
Figure 623082DEST_PATH_IMAGE035
(1)
其中，
Figure 39151DEST_PATH_IMAGE036
表示各阵元距离相位零点的水平间距，
Figure 495671DEST_PATH_IMAGE037
表示各阵元距离相位零点的竖直间距，
Figure 893155DEST_PATH_IMAGE038
表示天线列数，
Figure 497443DEST_PATH_IMAGE039
表示天线行数，
Figure 389306DEST_PATH_IMAGE040
表示天线列间距，
Figure 949601DEST_PATH_IMAGE041
表示天线行间距。
在一个具体实施例中，当相位零点为阵列天线的左下角天线，即坐标为(1,1)的天线为相位零点时，水平间距和竖直间距为：
Figure 268718DEST_PATH_IMAGE042
(2)
其中，
Figure 547252DEST_PATH_IMAGE036
表示各阵元距离相位零点的水平间距，
Figure 242807DEST_PATH_IMAGE043
表示各阵元距离相位零点的竖直间距，
Figure 657608DEST_PATH_IMAGE044
表示天线列数，
Figure 413205DEST_PATH_IMAGE045
表示天线行数，
Figure 913457DEST_PATH_IMAGE046
表示天线列间距，
Figure 158842DEST_PATH_IMAGE047
表示天线行间距。
S2、通过载频系数配置载频本振，并利用期望波束形成方向、载频本振、水平间距和竖直间距计算期望方向的空域导向矢量。具体包括步骤：
S21、通过载频系数配置载频本振。
具体的，获取载频系数
Figure 693728DEST_PATH_IMAGE048
，通过载频系数配置载频本振，计算公式如下：
Figure 354648DEST_PATH_IMAGE049
(3)
其中，
Figure 342196DEST_PATH_IMAGE050
表示射频本振，
Figure 379553DEST_PATH_IMAGE051
表示基底射频，
Figure 768946DEST_PATH_IMAGE052
表示载频系数，
Figure 866346DEST_PATH_IMAGE053
表示步进频率。
本实施例中，通过配置载频系数的大小可以实现射频本振的频率在基底射频的基础上步进达到使用者的射频要求，其中，基底射频与步进频率是在形成IP核前设置好的固定参数。通过该方法有效实现了载频可调，可调范围的最低值为
Figure 154239DEST_PATH_IMAGE054
，最高值可通过载频系数与步进进行设置。例如，本实例中设置基底射频为14.5GHz，步进频率为50MHz，则可以通过配置载频系数F为30使射频本振达到16GHz。
S22、利用载频本振、水平间距、竖直间距和期望波束形成方向计算期望方向在方位向的空域导向矢量和期望方向在俯仰向的空域导向矢量。
具体的，获取期望波束形成的方位角
Figure 978976DEST_PATH_IMAGE055
与俯仰角
Figure 973607DEST_PATH_IMAGE056
，利用载频本振、发射波波长、水平间距、竖直间距、期望波束形成的方位角
Figure 491176DEST_PATH_IMAGE057
、俯仰角
Figure 204049DEST_PATH_IMAGE058
计算期望方向在方位向的空域导向矢量和期望方向在俯仰向的空域导向矢量。
期望方向在方位向的空域导向矢量为：
Figure 645525DEST_PATH_IMAGE059
(4)
期望方向在俯仰向的空域导向矢量为：
Figure 743931DEST_PATH_IMAGE060
(5)
其中，
Figure 183134DEST_PATH_IMAGE061
表示期望方向在方位向的空域导向矢量，
Figure 632570DEST_PATH_IMAGE062
表示期望方向在俯仰向的空域导向矢量，
Figure 549842DEST_PATH_IMAGE063
，
Figure 768333DEST_PATH_IMAGE064
，
Figure 378437DEST_PATH_IMAGE065
表示步骤S12计算出的各阵元距离相位零点的水平间距，
Figure 49590DEST_PATH_IMAGE066
表示步骤S12计算出的各阵元距离相位零点的竖直间距，
Figure 36132DEST_PATH_IMAGE067
表示利用步骤S21配置出的载频本振计算出的发射波波长。
S23、将方位向的空域导向矢量和俯仰向的空域导向矢量进行罗内克积，得到期望方向的空域导向矢量。
具体的，将
Figure 109130DEST_PATH_IMAGE068
与
Figure 890135DEST_PATH_IMAGE069
进行克罗内克积（Kronecker积）便得到了期望方向的空域导向矢量，用
Figure 861634DEST_PATH_IMAGE070
表示期望方向的空域导向矢量：
Figure 369975DEST_PATH_IMAGE071
(6)
其中，
Figure 313792DEST_PATH_IMAGE072
表示期望方向的空域导向矢量。
具体的，期望方向的空域导向矢量的点数为
Figure 514966DEST_PATH_IMAGE073
。
S3、利用天线行数、天线列数和高斯窗口宽度参数计算目标高斯窗。具体包括步骤：
S31、将高斯窗输出点数分别取天线行数和天线列数，并利用高斯窗口宽度参数计算出第一高斯向量和第二高斯向量。
计算第一高斯向量：
具体的，首先获取高斯窗输出点数K：将高斯窗输出点数K取值为天线行数M，此时输出为M个值，计算M个数的离散程度得到M个数的标准差；然后获取高斯窗口宽度参数
Figure 911443DEST_PATH_IMAGE031
，高斯窗口宽度参数
Figure 223476DEST_PATH_IMAGE031
与M个数的标准差的倒数成正比，与窗宽成反比；最后，利用高斯窗计算第一高斯向量：
Figure 10080DEST_PATH_IMAGE074
(7)
其中，
Figure 382156DEST_PATH_IMAGE075
表示高斯窗点数，
Figure 265929DEST_PATH_IMAGE031
表示高斯窗口宽度。
上式中，
Figure 381653DEST_PATH_IMAGE076
的取值范围为：
Figure 34482DEST_PATH_IMAGE077
(8)
其中，
Figure 577459DEST_PATH_IMAGE076
取值
Figure 682949DEST_PATH_IMAGE075
次。
在一个具体实施例中，n取值为取值范围的最小值依次加1得到。例如，K=5，则n取值范围为-2到2，此时n取值为-2、-1、0、1、2；K=6，则n取值范围为-2.5到2.5，此时n取值为-2.5、-1.5、-0.5、0.5、1.5、2.5。
计算第二高斯向量：
具体的，首先获取高斯窗输出点数K：将高斯窗输出点数K取值为天线列数N，此时输出为N个值，计算N个数的离散程度得到N个数的标准差；然后获取计算第一高斯向量时所用的高斯窗口宽度参数
Figure 602364DEST_PATH_IMAGE031
，高斯窗口宽度
Figure 172016DEST_PATH_IMAGE031
与N个数的标准差的倒数成正比，与窗宽成反比；最后，利用式(7)和式(8)计算第二高斯向量。
S32、将第一高斯向量和第二高斯向量进行克罗内克积，计算出
Figure 636627DEST_PATH_IMAGE078
个高斯窗点数，得到目标高斯窗。
具体的，将两组高斯向量进行克罗内克积，计算出
Figure 744260DEST_PATH_IMAGE078
个的高斯窗点数，得到目标高斯窗。
需要说明的是，步骤S3中，窗函数的输出点数需要与阵列的行数列数相对应，例如阵列天线排布的行数为M，列数为N，则需要将K值分别取M与N计算出两组高斯点组成的向量，将这两组向量进行克罗内克积，得到
Figure 952519DEST_PATH_IMAGE078
个点。
本实施例中，如果高斯窗可以较好的与阵列排布相匹配，且当不需要高斯窗函数时可以设置宽度参数
Figure 563629DEST_PATH_IMAGE031
为0，则输出高斯窗值全为1，相当于对步骤S2计算的空域导向矢量没做加窗处理。
S4、将空域导向矢量与目标高斯窗相乘，计算得到波束形成权系数。
具体的，将步骤S2计算得到的
Figure 464720DEST_PATH_IMAGE079
个点的空域导向矢量与步骤S3计算得到的
Figure 794070DEST_PATH_IMAGE079
个点的目标高斯窗分别点乘得到点数为
Figure 524128DEST_PATH_IMAGE079
的波束形成权系数，该波束形成权系数为降低副瓣的波束形成权系数。
S5、对上述步骤所形成的程序进行综合与时序优化，并将优化后的程序生成可调用的通用化IP核。具体包括步骤：
S51、对各阵元距离相位零点的水平间距和竖直间距的计算进行优化。
具体的，在计算各阵元距离相位零点的水平间距和竖直间距时，多次使用for循环，可用PIPELINE指令对循环进行流水线设计，减少循环使用时间，因为该方法用到行数与列数内嵌for循环，所以可以通过UNROLL指令对内循环进行展开，降低for循环消耗时间，使得水平间距和竖直间距的计算可以分别进行，降低计算时间。
S52、对执行步骤的顺序进行优化。
具体的，对于步骤S3计算高斯窗点数组成的向量时，除了输入点数以外参数都相同，则可将高斯窗生成模块写为函数体，重读调用两次，实现并行运算，可降低时间，或者利用循环分时复用同一模块减小资源消耗。
本实施例中，通过上述步骤，可以使得步骤S2和步骤S3并行，从而降低计算时间，减小资源消耗。
S53、对输入输出接口进行接口优化。
具体的，对输入输出接口进行接口优化，在计算结果精度保证的前提下尽可能的减小寄存器的使用，极大了节省了FPGA的资源消耗。
S54、将优化后的程序生成可调用的通用化IP核。
具体的，将上述优化形成的程序生成可以被verilog直接调用的IP核，便于使用者直接使用，避免了开发者重复劳动，并且IP核使用方便简单易懂，可移植性高更有助于缩短算法开发周期。
进一步的，生成的通用化IP核可用于对不同排布的天线阵列计算波束形成权系数，从而进行相位补偿，便于后续对目标的测角测速。
请参见图3，图3为本发明实施例提供的一种交叉天线阵列排布图。当天线排布为一类非全部占满
Figure 271636DEST_PATH_IMAGE080
个天线位置的情况时，同样可以通过本实施例的通用化IP进行计算生成权系数。具体实现方法是同样设置天线行数为M，天线列数为N，计算出
Figure 343628DEST_PATH_IMAGE081
个点数的权矢量，因为输出结果为一组向量，该向量中每个元素的位置与天线阵列上每一个天线的位置为一一对应关系，所以对于特殊排布的阵列天线只需选取输出数组中天线对应位置的元素再重新组合成一个新的权矢量向量即可。
具体的，对每个天线接收到的数字信号乘以一个复数以进行相位补偿，该复数即为补偿值，多个补偿复数形成权矢量；每个天线所对应的补偿值不一样，计算M*N个点数的矩阵，得到补偿值所组成的矩阵，每个值在矩阵中的位置代表某个天线，从而可以通过把某个补偿值去掉的方式，组合成新的天线阵列。
本实施例的IP核生成方法具有以下优点：1）可以实现载频可调。对于不同的发射波形，其载频有可能存在一定差异，会导致一般波束形成方法对于不同载频波形的通用性较差，本实施例的载频可调极大地提高算法对于不同在的适用性。2）实现阵列天线排布自定义。对于不同排布的阵列天线，波束形成算法求权系数的过程存在差异，本实施例可以自定义天线的行数、列数、行间距、列间距与相位零点，能满足各式各样的不同排布方式的阵列天线，极大地提高了本实施例方法对于不同阵列天线排布的适用性。3）自定义配置窗函数。波束形成方法中添加窗函数可以有效降低副瓣，本实施例的自定义配置窗函数有效提高了天线的综合性能，形成高增益的窄波束。4）本实施例基于Xilinx的Vivado HLS工具采用C++语言编写，并对方法进行并行优化与时序优化，在保证计算精度的前提下，有效减小资源占用和提高实时性要求；基于Xilinx的Vivado HLS工具可以生成适用于Vivado工程的IP核，提高了算法的可移植性。5）本实施例技术成本低，算法耗时短，有效提高了实时性能。
综上，本实施例的本发明的IP核生成方法可以自定义相位零点，通过计算各阵元距离相位零点的距离可以实现阵列天线排布自定义，通过载频系数配置载频本振可以实现载频自定义，有效提高了该方法对于不同阵列天线排布的适用性，提高了该方法的通用性；结合天线行数、天线列数、高斯窗口宽度计算目标高斯窗可以实现窗函数自定义，有效提高了天线的综合性能，形成高增益的窄波束；另外，生成可调用的通用化IP核，可以提高该方法的可移植性，并且操作简单，方法实现简单易懂。本实施例的IP核生成方法利用Vivado-HLS工具采用C++语言实现，并对算法进行并行优化与时序优化，在保证计算精度的前提下，有效减小资源占用，提高了实时性要求。
实施例二
在实施例一和实施例二的基础上，本实施例提供了一种通用化二维数字波束形成的IP核，该IP核由实施例一或者实施例二的生成方法生成得到，包括天线阵列生成模块、期望方向导向矢量生成模块、高斯窗生成模块、波束形成权系数生成模块和IP核生成模块组成。
具体的，天线阵列生成模块用于配置阵列天线的天线行数、天线列数、天线行间距、天线列间距和相位零点，并计算阵列天线中各阵元距离相位零点的水平间距和竖直间距。期望方向导向矢量生成模块用于通过载频系数配置载频本振，并利用期望波束形成方向、载频本振、水平间距和竖直间距计算期望方向的空域导向矢量。高斯窗生成模块用于利用天线行数、天线列数和高斯窗口宽度参数计算目标高斯窗。波束形成权系数生成模块用于将空域导向矢量与高斯窗相乘，计算得到波束形成权系数。IP核生成模块用于对上述步骤所形成的程序进行综合与时序优化，并将优化后的程序生成可调用的通用化IP核。
上述各个模块中的执行步骤请参见实施例一，本实施例不再赘述。
本实施例的IP核由天线阵列生成模块、期望方向导向矢量生成模块、高斯窗生成模块、波束形成权系数生成模块和IP核生成模块组成，有效提高了该方法的通用性，其通用性主要包括，载频可自定义、天线阵列排布可自定义、相位零点可自定义、窗函数可自定义；并且在进行约束与优化之后，可以有效降低该核资源消耗与时间占用。本实施例生成可被调用的IP核既能提高可移植性，并且操作简易，方法实现简单易懂。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护范围。

Claims (7)

1.一种通用化二维数字波束形成技术的IP核生成方法，其特征在于，包括步骤：
配置阵列天线的天线行数、天线列数、天线行间距、天线列间距和相位零点，并计算所述阵列天线中各阵元距离所述相位零点的水平间距和竖直间距，所述相位零点包括所述阵列天线的中心位置或者所述阵列天线的左下角天线；
通过载频系数配置载频本振，并利用期望波束形成方向、所述载频本振、所述水平间距和所述竖直间距计算期望方向的空域导向矢量；包括：
通过所述载频系数配置所述载频本振：
Fc＝F0+F×F1
其中，Fc表示射频本振，F0表示基底射频，F表示载频系数，F1表示步进频率；
利用所述载频本振、所述水平间距、所述竖直间距和所述期望波束形成方向计算期望方向在方位向的空域导向矢量和期望方向在俯仰向的空域导向矢量：
Figure FDA0003812133330000011
Figure FDA0003812133330000012
其中，ay表示期望方向在方位向的空域导向矢量，az表示期望方向在俯仰向的空域导向矢量，
Figure FDA0003812133330000013
Dy表示各阵元距离相位零点的水平间距，Dz表示各阵元距离相位零点的竖直间距，λ表示利用载频本振计算出的发射波波长；
将所述方位向的空域导向矢量和所述俯仰向的空域导向矢量进行罗内克积，得到所述期望方向的空域导向矢量：
ω＝kron(ay，az)
其中，ω表示期望方向的空域导向矢量；
利用所述天线行数、所述天线列数和高斯窗口宽度参数计算目标高斯窗；
将所述期望方向的空域导向矢量与所述目标高斯窗相乘，计算得到波束形成权系数；
对上述步骤所形成的程序进行综合与时序优化，并将优化后的程序生成可调用的通用化IP核。
2.根据权利要求1所述的通用化二维数字波束形成技术的IP核生成方法，其特征在于，计算所述阵列天线中各阵元距离所述相位零点的水平间距和竖直间距，包括：
利用所述天线行数、所述天线列数、所述天线行间距、所述天线列间距计算各阵元距离所述相位零点的所述水平间距和所述竖直间距。
3.根据权利要求1所述的通用化二维数字波束形成技术的IP核生成方法，其特征在于，当所述相位零点为所述阵列天线的中心位置时，所述水平间距和所述竖直间距为：
Figure FDA0003812133330000021
Figure FDA0003812133330000022
其中，Dy表示各阵元距离相位零点的水平间距，Dz表示各阵元距离相位零点的竖直间距，N表示天线列数，M表示天线行数，dy表示天线列间距，dz表示天线行间距。
4.根据权利要求1所述的通用化二维数字波束形成技术的IP核生成方法，其特征在于，当所述相位零点为所述阵列天线的左下角天线时，所述水平间距和所述竖直间距为：
Dy＝[0，1，2，3…(N-1)]×dy
Dz＝[0，1，2，3…(M-1)]×dz
其中，Dy表示各阵元距离相位零点的水平间距，Dz表示各阵元距离相位零点的竖直间距，N表示天线列数，M表示天线行数，dy表示天线列间距，dz表示天线行间距。
5.根据权利要求1所述的通用化二维数字波束形成技术的IP核生成方法，其特征在于，利用所述天线行数、所述天线列数和高斯窗口宽度参数计算目标高斯窗，包括：
将高斯窗输出点数分别取所述天线行数和所述天线列数，并利用所述高斯窗口宽度参数计算出第一高斯向量和第二高斯向量，计算公式为：
Figure FDA0003812133330000023
其中，K表示高斯窗点数，α表示高斯窗口宽度参数，
Figure FDA0003812133330000024
n取值K次；
将所述第一高斯向量和所述第二高斯向量进行克罗内克积，计算出M×N个高斯窗点数，得到所述目标高斯窗。
6.根据权利要求1所述的通用化二维数字波束形成技术的IP核生成方法，其特征在于，所述IP核生成方法利用Vivado-HLS工具采用C++语言实现。
7.一种通用化二维数字波束形成的IP核，其特征在于，由如权利要求1-6任一项所述的通用化二维数字波束形成技术的IP核生成方法生成得到，包括：
天线阵列生成模块，用于配置阵列天线的天线行数、天线列数、天线行间距、天线列间距和相位零点，并计算所述阵列天线中各阵元距离所述相位零点的水平间距和竖直间距；
期望方向导向矢量生成模块，用于通过载频系数配置载频本振，并利用期望波束形成方向、所述载频本振、所述水平间距和所述竖直间距计算期望方向的空域导向矢量；
高斯窗生成模块，用于利用所述天线行数、所述天线列数和高斯窗口宽度参数计算目标高斯窗；
波束形成权系数生成模块，用于将所述空域导向矢量与所述目标高斯窗相乘，计算得到波束形成权系数；
IP核生成模块，用于对上述步骤所形成的程序进行综合与时序优化，并将优化后的程序生成可调用的通用化IP核。
CN202210782939.2A 2022-07-05 2022-07-05 一种通用化二维数字波束形成技术的ip核及生成方法 Active CN114840052B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210782939.2A CN114840052B (zh) 2022-07-05 2022-07-05 一种通用化二维数字波束形成技术的ip核及生成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210782939.2A CN114840052B (zh) 2022-07-05 2022-07-05 一种通用化二维数字波束形成技术的ip核及生成方法

Publications (2)

Publication Number Publication Date
CN114840052A CN114840052A (zh) 2022-08-02
CN114840052B true CN114840052B (zh) 2022-09-30

Family

ID=82574575

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210782939.2A Active CN114840052B (zh) 2022-07-05 2022-07-05 一种通用化二维数字波束形成技术的ip核及生成方法

Country Status (1)

Country Link
CN (1) CN114840052B (zh)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09261008A (ja) * 1996-01-17 1997-10-03 Atr Kodenpa Tsushin Kenkyusho:Kk ディジタルビーム形成装置のためのｆｉｒ型ディジタルフィルタのフィルタ係数演算装置、ディジタルビーム形成装置のためのｆｉｒ型ディジタルフィルタ及びディジタルビーム形成装置

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1075117A (ja) * 1996-08-29 1998-03-17 Mitsubishi Electric Corp アンテナ装置
CN101995566B (zh) * 2010-10-15 2013-01-23 西安电子科技大学 二维数字阵列雷达数字波束形成系统及方法
CN107356943B (zh) * 2017-06-01 2020-04-28 西南电子技术研究所（中国电子科技集团公司第十研究所） 数字波束形成和相位拟合方法
CN107979404B (zh) * 2017-10-27 2020-08-11 西安电子科技大学 基于虚拟阵列零陷加宽的自适应波束形成方法
WO2020028584A1 (en) * 2018-08-03 2020-02-06 Viasat, Inc. Antenna array system with disparate beam forming networks and non-linear filtering to mitigate interference
CN111398907B (zh) * 2019-01-02 2023-04-07 中国人民解放军空军预警学院 一种相控阵雷达空域二维数字差波束形成方法
US11217889B2 (en) * 2019-08-13 2022-01-04 The Boeing Company Method to optimize beams for phased array antennas
US11916305B2 (en) * 2020-07-01 2024-02-27 Linquest Corporation Systems and methods for massive phased arrays via beam-domain processing
CN111865385B (zh) * 2020-07-20 2021-07-02 电子科技大学 一种基于fpga的二维平面阵数字波束形成方法
CN111812627B (zh) * 2020-07-20 2022-05-17 电子科技大学 一种基于fpga的可控多波束形成方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09261008A (ja) * 1996-01-17 1997-10-03 Atr Kodenpa Tsushin Kenkyusho:Kk ディジタルビーム形成装置のためのｆｉｒ型ディジタルフィルタのフィルタ係数演算装置、ディジタルビーム形成装置のためのｆｉｒ型ディジタルフィルタ及びディジタルビーム形成装置

Also Published As

Publication number Publication date
CN114840052A (zh) 2022-08-02

Similar Documents

Publication Publication Date Title
CN108508423B (zh) 基于异型阵的子阵数字和差单脉冲测角方法
CN108051782B (zh) 基于子阵划分的大规模相控阵差波束形成系统
CN105606906B (zh) 一种毫米波相控阵测试标定方法
CN110161452B (zh) 基于互质式l型电磁矢量传感器阵列的波达方向估计方法
CN104391183A (zh) 近场测量天线远区场特性快速计算方法
CN113162670A (zh) 数字多波束校正与合成方法
Gu et al. Design of wide scanning sparse planar array using both matrix-pencil and space-mapping methods
CN113126047A (zh) 基于远场点源的自适应通道校正方法
CN110244273A (zh) 一种基于均匀分布式阵列的目标角度估计方法
CN110837075A (zh) 一种低复杂度的极化参数估计跟踪装置及方法
CN113419208B (zh) 干扰源实时测向方法、装置、存储介质及电子设备
CN108667489B (zh) 多波束波形发射方法及系统
CN114840052B (zh) 一种通用化二维数字波束形成技术的ip核及生成方法
CN115728722A (zh) 4d雷达的天线阵列、数据探测方法及4d雷达
CN115825944B (zh) 基于外辐射源雷达的单快拍多目标来波方向估计方法
CN109490846B (zh) 基于空时联合优化的多输入多输出雷达波形设计方法
CN111162878A (zh) 一种基于子阵降维带约束的多域联合抗干扰方法
CN102904015B (zh) 一种短波小型圆形接收天线阵
CN202871978U (zh) 一种短波小型圆形接收天线阵
Ghafi et al. Reconfigurable linear antenna arrays for beam-pattern matching in collocated MIMO radars
CN114239395A (zh) 机载柔性端射阵列天线的波束赋形与稀布重构优化方法
CN115037352A (zh) 一种面向卫星波束对准增强的降噪方法
CN114859311A (zh) 一种基于Vivado-HLS的LCMV旁瓣抑制方法
CN114563760A (zh) 一种基于sca阵型的二阶超波束形成方法、设备及介质
CN111859272A (zh) 一种面向大规模天线的快速music谱分解方法、装置及设备

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