CN111551908A - Method for reducing complexity of array element activation algorithm of phased array system - Google Patents

Method for reducing complexity of array element activation algorithm of phased array system Download PDF

Info

Publication number
CN111551908A
CN111551908A CN202010233508.1A CN202010233508A CN111551908A CN 111551908 A CN111551908 A CN 111551908A CN 202010233508 A CN202010233508 A CN 202010233508A CN 111551908 A CN111551908 A CN 111551908A
Authority
CN
China
Prior art keywords
array
activation
subarray
target
included angle
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
Application number
CN202010233508.1A
Other languages
Chinese (zh)
Other versions
CN111551908B (en
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.)
Southwest Electronic Technology Institute No 10 Institute of Cetc
Original Assignee
Southwest Electronic Technology Institute No 10 Institute of Cetc
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 Southwest Electronic Technology Institute No 10 Institute of Cetc filed Critical Southwest Electronic Technology Institute No 10 Institute of Cetc
Priority to CN202010233508.1A priority Critical patent/CN111551908B/en
Publication of CN111551908A publication Critical patent/CN111551908A/en
Application granted granted Critical
Publication of CN111551908B publication Critical patent/CN111551908B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • 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
    • G06—COMPUTING OR CALCULATING; 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)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Mathematical Optimization (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Remote Sensing (AREA)
  • Data Mining & Analysis (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention provides a method for reducing the complexity of an array element activation algorithm of a phased array system, and aims to provide a method capable of reducing the calculated amount and remarkably reducing the calculation complexity. The invention is realized by the following technical scheme: firstly, dividing every M adjacent array elements into a sub-array, dividing the array elements into N/M sub-arrays according to the total number N of the array elements on the basis of sub-array division, and constructing an activation calculation module; then, the maximum included angle between the direction vector pointed by each subarray and the pointing direction of all array elements in the subarray is calculated off-line and stored in a storage unit of an activation calculation module; when the activation area is calculated, aiming at each target wave beam, the activation calculation module reads the direction vector pointed by each subarray from the storage unit, then calculates the included angle between each subarray and the antenna wave beam pointing target in sequence, and when some subarrays cannot judge whether the activation is carried out, the included angles between the array elements in the subarrays and the antenna wave beam pointing target are calculated.

Description

降低相控阵系统阵元激活算法复杂度的方法A Method for Reducing the Complexity of Element Activation Algorithm in Phased Array System

技术领域technical field

本发明涉及阵列信号处理技术领域，降低相控阵系统阵元激活算法复杂度的方法，尤其是适用于大型共形相控阵系统中降低计算复杂度的天线阵元激活方法。The invention relates to the technical field of array signal processing, and a method for reducing the complexity of an array element activation algorithm of a phased array system, in particular to an antenna array element activation method suitable for reducing the computational complexity in a large conformal phased array system.

背景技术Background technique

阵列信号处理是信号处理领域的一个重要分支，具有广泛的应用。随着技术的发展，共形相控阵利用阵列信号处理技术，可以实现多波束、全空域的目标搜索、跟踪。相比较于传统的相控阵技术，共形相控阵由于阵元数目多、硬件复杂，因此极大的提高了计算复杂度，对阵列信号处理算法提出了更高的要求。阵列具有最大响应的方向被称为波束指向方向，即在这一方向阵列具有最大增益。对一个线性阵列，当信号在增益与相位上不做改变就进行合并时的波束指向面为阵列的宽面，并垂直于阵元的连线。阵列方向图在波束指向方向上的任一侧都将衰减到零值，即在该位置阵列响应为零，通常称为零陷。在波束指向方向上的两侧零陷间的方向图被称为主瓣。主瓣两个功率点间的宽度称为半功率波束宽度。所谓“相控阵”即通过相位控制实现波束扫描的阵列，其相位值可以通过计算机灵活改变。正是由于它的这种灵活性，使得相控阵雷达能够利用同一天线口径形成多个独立的发射波束和接收波束，且能根据实际情况的需求确定雷达的最佳工作方式，以获得满足各种要求的复杂波束。对于大型的相控阵雷达，通常包含数以百计甚至上千的天线阵元。为解决大型相控阵雷达阵元级波束形成技术面临的硬件系统复杂、实时性低等问题，一般采用子阵级波束形成技术来处理。但子阵级处理往往会破坏静态方向图的性能。在信号处理过程中，如果采用阵元级数字波束形成方式，则其运算量是非常庞大的，且其对应的硬件系统很复杂，成本也相对较高。如果采用传统稀疏阵列算法，如遗传算法、模拟退火算法等，具有随机性，需要很长的时间才能得到结果，对于大型阵列这些算法很难实现。Array signal processing is an important branch in the field of signal processing and has a wide range of applications. With the development of technology, conformal phased array can achieve multi-beam, full-space target search and tracking by using array signal processing technology. Compared with the traditional phased array technology, the conformal phased array greatly increases the computational complexity due to the large number of array elements and complex hardware, and puts forward higher requirements for the array signal processing algorithm. The direction in which the array has the greatest response is called the beam pointing direction, that is, the direction in which the array has the greatest gain. For a linear array, when the signals are combined without changing the gain and phase, the beam pointing surface is the broad surface of the array and is perpendicular to the line connecting the array elements. The array pattern will attenuate to a value of zero on either side of the beam pointing direction, ie the array response is zero at that location, commonly referred to as nulling. The pattern between the nulls on both sides in the beam pointing direction is called the main lobe. The width between the two power points of the main lobe is called the half-power beamwidth. The so-called "phased array" is an array that realizes beam scanning through phase control, and its phase value can be flexibly changed by a computer. It is precisely because of its flexibility that the phased array radar can use the same antenna aperture to form multiple independent transmit beams and receive beams, and can determine the best working mode of the radar according to the needs of the actual situation, so as to meet the needs of each complex beams required. For large phased array radars, it usually contains hundreds or even thousands of antenna elements. In order to solve the problems of complex hardware system and low real-time performance faced by the large-scale phased array radar array element-level beamforming technology, the sub-array-level beamforming technology is generally used to deal with it. But subarray-level processing tends to destroy the performance of static patterns. In the process of signal processing, if the array element-level digital beamforming method is adopted, the amount of computation is very large, and the corresponding hardware system is very complicated and the cost is relatively high. If traditional sparse array algorithms, such as genetic algorithm, simulated annealing algorithm, etc., are random, it takes a long time to get results, which is difficult to implement for large arrays.

相控阵系统阵面天线的波束指向由波束控制系统来执行，它主要通过对阵面各单元相位和增益的控制实现波束空间指向的变化。其中各单元相位变化对确定的阵列天线而言主要取决于天线波束指向角的变化。波束控制计算机根据波束指向要求对阵面各单元点的相位、幅度进行统一运算后将相位、幅度等数据分别传输至阵面各点，当相控阵面单元较多时，其计算量大，运算时间影响了波束扫描的速度。同时，随着相控阵面规模的增大，波束控制系统也越来越复杂，此时常规集中计算方法对数字信号处理器产生巨大压力，严重影响了波束控制的响应时间。由于现代相控阵电子系统对波束控制的速度要求越来越高，因此对系统的波束运算、数据传输等要求也相应提高。考虑到共形相控阵具有遮挡效应，即对于某个指向的目标，有些阵元对于该目标是可见的，有些阵元是不可见的。因此，对于某个指向的目标，并不是所有阵元都能接收到该目标的信号。在形成波束时，不可见的阵元或边缘的阵元贡献很小，传统的方法是依次计算每个阵元与目标的夹角。对于大型共形相控阵系统，由于阵元数目庞大，目标波束数目多，一般采用阵元激活的方式，来确定哪些阵元参与到波束合成，不激活的阵元不参与波束合成。传统的阵元激活的策略是依次判断每个阵元与目标的夹角，若夹角小于激活角度，则该阵元为激活阵元，否则为不激活阵元。采用传统的激活判断方法，存在计算量大的问题。The beam pointing of the front antenna of the phased array system is carried out by the beam control system, which mainly realizes the change of the beam spatial direction by controlling the phase and gain of each element of the plane. The phase change of each element mainly depends on the change of the antenna beam pointing angle for the determined array antenna. The beam control computer performs a unified operation on the phase and amplitude of each element point on the surface according to the beam pointing requirements, and then transmits the phase, amplitude and other data to each point of the array respectively. When there are many phased array elements, the calculation amount is large and the operation time Affects the speed of beam scanning. At the same time, as the size of the phased array increases, the beam steering system becomes more and more complex. At this time, the conventional centralized computing method puts a huge pressure on the digital signal processor, which seriously affects the response time of the beam steering. Since modern phased array electronic systems have higher and higher requirements for beam control speed, the requirements for beam operation and data transmission of the system are also increased accordingly. Considering that the conformal phased array has an occlusion effect, that is, for a certain pointing target, some array elements are visible to the target, and some array elements are invisible. Therefore, for a pointed target, not all array elements can receive the signal of the target. When forming a beam, the invisible elements or edge elements contribute very little. The traditional method is to calculate the angle between each element and the target in turn. For a large-scale conformal phased array system, due to the large number of array elements and the large number of target beams, the array element activation method is generally used to determine which array elements participate in beamforming, and the inactive array elements do not participate in beamforming. The traditional array element activation strategy is to judge the included angle between each array element and the target in turn. If the included angle is smaller than the activation angle, the array element is an active array element, otherwise it is an inactive array element. Using the traditional activation judgment method has the problem of a large amount of calculation.

发明内容SUMMARY OF THE INVENTION

本发明的目的是针对现有大型共形相控阵系统阵元激活存在的不足之处，提供一种能够减少计算量，计算速度快，收敛快，性能良好，并能显著降低相控阵系统阵元激活算法复杂度的方法。The purpose of the present invention is to provide a method that can reduce the amount of calculation, has fast calculation speed, fast convergence, good performance, and can significantly reduce the number of phased array system array elements. A method for meta-activation algorithmic complexity.

本发明的上述目的可以通过以下措施来达到。一种降低相控阵系统阵元激活算法复杂度的方法，具有如下技术特征：首先将相控阵系统进行子阵划分，即将每M个相邻阵元划分为一个子阵，在子阵划分基础上，根据阵元的总数N划分为N/M个子阵，并划分激活区域，构建激活计算模块；然后离线计算出每个子阵指向的方向矢量与子阵内所有阵元指向的最大夹角，并存储在激活计算模块的存储单元中；计算激活区域时，针对每个目标波束，激活计算模块从存储单元中读取出每个子阵指向的方向矢量，然后依次计算每个子阵与天线波束指向目标的夹角，当某些子阵无法判断是否激活时，再计算这些子阵内的阵元与天线波束指向目标的夹角。The above objects of the present invention can be achieved by the following measures. A method for reducing the complexity of an array element activation algorithm of a phased array system has the following technical features: first, the phased array system is divided into sub-arrays, that is, every M adjacent array elements are divided into a sub-array, and the sub-array is divided into sub-arrays. On the basis, according to the total number N of array elements, it is divided into N/M sub-arrays, and the activation area is divided to build an activation calculation module; then the direction vector pointed by each sub-array and the maximum angle pointed by all the array elements in the sub-array is calculated offline. , and stored in the storage unit of the activation calculation module; when calculating the activation area, for each target beam, the activation calculation module reads the direction vector pointed to by each subarray from the storage unit, and then calculates each subarray and the antenna beam in turn. The included angle pointing to the target, when some subarrays cannot be determined to be activated, then calculate the included angle between the array elements in these subarrays and the antenna beam pointing to the target.

本发明相比于现有技术的有益效果是：Compared with the prior art, the present invention has the following beneficial effects:

显著降低计算复杂度。本发明将相控阵系统进行子阵划分，经过子阵划分后，大型相控阵雷达的接收通道数明显较少，子阵数目显著减少，从而大大降低了系统的复杂度、数据处理量和硬件成本。划分的子阵可以克服栅瓣和栅零点的产生，提高了波束性能；在子阵划分基础上，可以结合子阵加权对子阵级自适应和波束进行优化，得到了副瓣较低、干扰抑制较好的自适应波束。阵列结构在子阵级自适应波束形成，使干扰抑制效果较好。由于子阵数目显著减少，可以显著减少与目标夹角的计算次数，降低了计算复杂度。Significantly reduces computational complexity. The invention divides the phased array system into sub-arrays. After the sub-array division, the number of receiving channels of the large phased array radar is significantly reduced, and the number of sub-arrays is significantly reduced, thereby greatly reducing the complexity of the system, the amount of data processing and the hardware cost. The divided sub-arrays can overcome the generation of grating lobes and grating zeros, and improve the beam performance; on the basis of sub-array division, sub-array level adaptation and beam optimization can be combined with sub-array weighting, resulting in lower side lobes, less interference Better suppression of adaptive beams. The array structure is adaptive to beamforming at the sub-array level, so that the interference suppression effect is better. Since the number of sub-arrays is significantly reduced, the number of calculations of the included angle with the target can be significantly reduced, and the calculation complexity is reduced.

减少计算量，收敛变快。本发明采用离线计算出每个子阵阵列信号指向的方向矢量，并将每个子阵指向的方向矢量存储在激活计算模块的存储单元中；计算激活区域时，针对每个目标波束，激活计算模块从存储单元中读取出每个子阵指向的方向矢量，然后依次计算每个子阵与天线波束指向目标的夹角，当某些子阵无法判断是否激活时，再计算这些子阵内的阵元与目标的夹角。因此可以显著减少计算量，收敛变快。Reduce the amount of calculation, the convergence is faster. The invention adopts offline calculation to calculate the direction vector pointed by each sub-array array signal, and stores the direction vector pointed to by each sub-array in the storage unit of the activation calculation module; when calculating the activation area, for each target beam, the activation calculation module starts from The direction vector pointed to by each sub-array is read out from the storage unit, and then the angle between each sub-array and the antenna beam pointing to the target is calculated in turn. target angle. Therefore, the amount of calculation can be significantly reduced, and the convergence is faster.

本发明的适用于大型共形相控阵系统的阵元激活判断。The invention is suitable for the activation judgment of the array element of the large-scale conformal phased array system.

附图说明Description of drawings

图1是本发明降低相控阵系统阵元激活算法复杂度的流程图。FIG. 1 is a flow chart of the present invention for reducing the complexity of an array element activation algorithm of a phased array system.

图2是采用本发明前后阵元激活判断计算量的对比示意图。FIG. 2 is a schematic diagram showing the comparison of the calculation amount of the activation judgment of array elements before and after using the present invention.

下面结合附图和实施实例对本发明进一步说明。The present invention will be further described below in conjunction with the accompanying drawings and implementation examples.

具体实施方式Detailed ways

参阅图1。根据本发明，首先将相控阵系统进行子阵划分，即将每M个相邻阵元划分为一个子阵，在子阵划分基础上，根据阵元的总数N划分为N/M个子阵，并划分激活区域，构建激活计算模块；然后离线计算出每个子阵指向的方向矢量与子阵内所有阵元指向的最大夹角，并存储在激活计算模块的存储单元中；计算激活区域时，针对每个目标波束，激活计算模块从存储单元中读取出每个子阵指向的方向矢量，然后依次计算每个子阵与天线波束指向目标的夹角，当某些子阵无法判断是否激活时，再计算这些子阵内的阵元与天线波束指向目标的夹角。See Figure 1. According to the present invention, the phased array system is firstly divided into sub-arrays, that is, every M adjacent array elements are divided into a sub-array, and on the basis of sub-array division, N/M sub-arrays are divided according to the total number N of array elements, And divide the activation area to build the activation calculation module; then calculate offline the direction vector pointed by each sub-array and the maximum angle pointed to by all the array elements in the sub-array, and store it in the storage unit of the activation calculation module; when calculating the activation area, For each target beam, the activation calculation module reads the direction vector pointed by each sub-array from the storage unit, and then calculates the angle between each sub-array and the antenna beam pointing to the target in turn. When some sub-arrays cannot be determined to be activated, Then calculate the angle between the array elements in these sub-arrays and the antenna beam pointing to the target.

在可选的实施例中，根据相控阵系统子阵划分获得的总共N/M个子阵，将相控阵系统阵元激活分为子阵激活判断和阵元激活判断两个阶段。在子阵激活阶段，激活计算模块计算每个子阵与目标的夹角，若子阵与目标的夹角小于某个门限时，则激活计算模块判断该子阵里的M个阵元全部激活；若子阵与目标的夹角大于某个门限时，则激活计算模块判断该子阵里的M个阵元全部不激活；在不满足以上子阵与目标的夹角＜或＞某个门限的两个条件时，则激活计算模块进入阵元激活判断阶段。在阵元激活判断阶段，激活计算模块对该子阵内的M个阵元依次判断与目标的夹角，当某个阵元与目标夹角小于激活角度时，则判断该阵元激活，否则不激活。In an optional embodiment, according to a total of N/M sub-arrays obtained by dividing the phased array system into sub-arrays, the phased array system array element activation is divided into two stages: sub-array activation judgment and array element activation judgment. In the subarray activation stage, the activation calculation module calculates the angle between each subarray and the target. If the angle between the subarray and the target is less than a certain threshold, the activation calculation module determines that all M array elements in the subarray are activated; When the angle between the array and the target is greater than a certain threshold, the activation calculation module determines that all M array elements in the subarray are not activated; When the conditions are met, the activation calculation module enters the array element activation judgment stage. In the phase of judging the activation of the array element, the activation calculation module judges the angle between the M array elements in the sub-array and the target in turn. When the angle between an array element and the target is smaller than the activation angle, the array element is judged to be activated, otherwise Not activated.

为了减少计算量，子阵的划分采用约定好的固定划分方式，每M个相邻阵元划分为一个子阵，然后离线计算出每个子阵指向的方向矢量和每个子阵指向与所有与子阵内所有阵元指向的最大夹角，并进行初始化，存储在激活计算模块的存储单元中。激活计算模块从存储单元中读取出每个子阵指向的方向矢量，设置子阵编号n＝1，计算第n个子阵目标的夹角，判断目标夹角小于激活角与最大夹角θmax之差，是则全部激活子阵内所有阵元，子阵编号n＝n+1，否则，判断子阵指向与目标的夹角大于激活角度与最大夹角θmax之和，是则全部激活子阵内所有阵元，子阵编号n＝n+1，否则，计算个子阵内每个阵元与目标的夹角，根据计算结果，判断夹角是否小于激活角，是则该子阵被激活，完成该子阵内所有阵元的激活判断，否则，该子阵不被激活，再继续计算个子阵内每个阵元与目标的夹角，直到完成该子阵内所有阵元的激活判断。子阵编号n＝n+1，判断子阵编号n是否大于子阵数目，是则结束程序，否则，返回计算第n个子阵目标的夹角，直到子阵编号n大于子阵数目结束程序。In order to reduce the amount of calculation, the division of the sub-array adopts the agreed fixed division method. Every M adjacent array elements is divided into a sub-array, and then the direction vector pointed to by each sub-array and the direction vector pointed to by each sub-array are calculated offline. The maximum included angle pointed to by all array elements in the array is initialized and stored in the storage unit of the activation calculation module. The activation calculation module reads the direction vector pointed to by each sub-array from the storage unit, sets the sub-array number n=1, calculates the included angle of the nth sub-array target, and judges that the target included angle is less than the activation angle and the maximum included angle θ max . If the difference is true, all array elements in the sub-array are activated, and the sub-array number is n=n+1. Otherwise, it is judged that the angle between the sub-array and the target is greater than the sum of the activation angle and the maximum angle θ max . If yes, all the sub-arrays are activated. All array elements in the array, the subarray number n=n+1, otherwise, calculate the angle between each array element in the subarray and the target, according to the calculation result, judge whether the angle is smaller than the activation angle, if yes, the subarray is activated , complete the activation judgment of all the array elements in the sub-array, otherwise, the sub-array will not be activated, and then continue to calculate the angle between each array element in the sub-array and the target, until the activation judgment of all the array elements in the sub-array is completed . Subarray number n=n+1, judge whether the subarray number n is greater than the number of subarrays, if so, end the program, otherwise, return to calculate the angle of the nth subarray target, until the subarray number n is greater than the number of subarrays to end the program.

激活计算模块根据子阵指向的方向矢量由子阵内包含的M个阵元指向的方向矢量求和，并作归一化处理，得到子阵指向的方向矢量

Figure BDA0002430194920000041
即The activation calculation module sums the direction vectors pointed by the M array elements contained in the sub-array according to the direction vector pointed by the sub-array, and normalizes it to obtain the direction vector pointed by the sub-array
Figure BDA0002430194920000041
which is

Figure BDA0002430194920000042
Figure BDA0002430194920000042

式中，

Figure BDA0002430194920000043
为阵元i与相控阵原点连线的单位向量，即阵元指向的方向矢量，
Figure BDA0002430194920000044
为
Figure BDA0002430194920000045
的2范数。In the formula,
Figure BDA0002430194920000043
is the unit vector connecting the array element i and the origin of the phased array, that is, the direction vector pointed by the array element,
Figure BDA0002430194920000044
for
Figure BDA0002430194920000045
2 norm of .

激活计算模块根据最大夹角θmax计算公式Activate the calculation module to calculate the formula according to the maximum angle θ max

Figure BDA0002430194920000046
Figure BDA0002430194920000046

离线计算每个子阵指向与该子阵内M个阵元指向的最大夹角θmax，并存储在激活计算模块的存储单元中。The maximum angle θ max pointed to by each sub-array and the M array elements in the sub-array is calculated offline, and stored in the storage unit of the activation calculation module.

对于每个目标波束，计算激活区域时，激活计算模块首先进行子阵激活判断，对于每个子阵，激活计算模块从存储单元中读取该子阵指向的方向矢量，然后计算与目标的夹角，即

Figure BDA0002430194920000047
其中，上式中的
Figure BDA0002430194920000048
表示目标指向的方向矢量。For each target beam, when calculating the activation area, the activation calculation module firstly determines the activation of the sub-array. For each sub-array, the activation calculation module reads the direction vector pointed by the sub-array from the storage unit, and then calculates the angle between the target and the target. ,Right now
Figure BDA0002430194920000047
Among them, in the above formula
Figure BDA0002430194920000048
Represents the direction vector to which the target is pointing.

若子阵与目标的夹角小于激活角度与θmax之差，则激活计算模块判断该子阵内的M个阵元全部激活；若子阵与目标的夹角大于激活角度与θmax之和，则激活计算模块判断该子阵内的M个阵元全部不激活；否则，该子阵内的M个阵元部分激活，此时激活计算模块进行阵元激活判断，依次判断子阵内的M个阵元是否激活。If the angle between the sub-array and the target is less than the difference between the activation angle and θ max , the activation calculation module determines that all M array elements in the sub-array are activated; if the angle between the sub-array and the target is greater than the sum of the activation angle and θ max , then The activation calculation module judges that the M array elements in the sub-array are all inactive; otherwise, the M array elements in the sub-array are partially activated. At this time, the activation calculation module judges the activation of the array elements, and then judges the M array elements in the sub-array in turn. Whether the array element is activated.

当激活计算模块需要依次判断子阵内的M个阵元是否需要激活时，激活计算模块从存储单元中依次读取M个阵元指向的方向矢量，并计算阵元与目标的夹角，即When the activation calculation module needs to judge in turn whether the M array elements in the subarray need to be activated, the activation calculation module reads the direction vectors pointed to by the M array elements from the storage unit in turn, and calculates the angle between the array element and the target, that is

Figure BDA0002430194920000049
Figure BDA0002430194920000049

Figure BDA0002430194920000051
Figure BDA0002430194920000051

Figure BDA0002430194920000052
若阵元与目标的夹角小于激活角度，则激活计算模块判断该阵元激活，否则不激活。
Figure BDA0002430194920000052
If the angle between the array element and the target is smaller than the activation angle, the activation calculation module determines that the array element is activated, otherwise it is not activated.

参阅图2。图2给出了采用本发明前后对于不同目标激活计算复杂度的仿真结果，该仿真基于某项目的球形相控阵系统，包含105个子阵，每个子阵16个阵元，激活角度为60度。从仿真结果可以看出，采用本发明后，可以显著降低与目标夹角的计算次数，平均计算次数为传统方法的16.7％，从而显著降低了计算复杂度。See Figure 2. Fig. 2 shows the simulation results of different target activation computational complexity before and after adopting the present invention. The simulation is based on a spherical phased array system of a certain project, which includes 105 sub-arrays, each sub-array has 16 array elements, and the activation angle is 60 degrees. . It can be seen from the simulation results that the calculation times of the included angle with the target can be significantly reduced after the invention is adopted, and the average calculation times are 16.7% of the traditional method, thereby significantly reducing the calculation complexity.

Claims (10)

1. A method for reducing the computational complexity of an active array element of a phased array system has the following technical characteristics: firstly, carrying out subarray division on a phased array system, namely dividing every M adjacent array elements into a subarray, on the basis of the subarray division, dividing the phased array system into N/M subarrays according to the total number N of the array elements, dividing an activation area, and constructing an activation calculation module; then, the maximum included angle between the direction vector pointed by each subarray and the pointing direction of all array elements in the subarray is calculated off-line and stored in a storage unit of an activation calculation module; when the activation area is calculated, aiming at each target wave beam, the activation calculation module reads the direction vector pointed by each subarray from the storage unit, then calculates the included angle between each subarray and the antenna wave beam pointing target in sequence, and when some subarrays cannot judge whether the activation is carried out, the included angles between the array elements in the subarrays and the antenna wave beam pointing target are calculated.
2. The method for reducing the complexity of the phased array system array element activation algorithm according to claim 1, wherein: and according to the phased array system subarray division, obtaining N/M total subarrays, and dividing the phased array system array element activation into a subarray activation judgment stage and an array element activation judgment stage.
3. The method for reducing the complexity of the phased array system array element activation algorithm according to claim 2, wherein: in the sub-array activation stage, an activation calculation module calculates an included angle between each sub-array and a target, and if the included angle between each sub-array and the target is smaller than a certain threshold, the activation calculation module judges that all M array elements in the sub-array are activated; and if the included angle between the subarray and the target is larger than a certain threshold, the activation calculation module judges that all M array elements in the subarray are not activated.
4. The method for reducing the complexity of the phased array system array element activation algorithm according to claim 3, wherein: and when the two conditions that the included angle between the subarray and the target is less than or greater than a certain threshold are not met, the activation calculation module enters an array element activation judgment stage, in the array element activation judgment stage, the activation calculation module sequentially judges the included angle between M array elements in the subarray and the target, when the included angle between a certain array element and the target is less than the activation angle, the activation of the array element is judged, and otherwise, the activation is not carried out.
5. The method for reducing the complexity of the phased array system array element activation algorithm according to claim 1, wherein: the division of the subarrays adopts a well-agreed fixed division mode, every M adjacent array elements are divided into one subarray, then the direction vector pointed by each subarray and the maximum included angle between the pointing direction of each subarray and the pointing direction of all the array elements in the subarrays are calculated off line, initialization is carried out, and the vector is stored in a storage unit of an activation calculation module.
6. The method for reducing the complexity of the phased array system array element activation algorithm according to claim 1, wherein: the activation calculation module reads out the direction vector pointed by each subarray from the storage unit, sets the subarray number n as 1, calculates the included angle between the nth subarray and the target, judges that the included angle between the nth subarray and the target is smaller than the difference between the activation angle and the maximum included angle, if yes, all array elements in the subarrays are activated, the subarray number n as n +1, otherwise, judges that the included angle pointed by the subarray and the target is larger than the sum of the activation angle and the maximum included angle, if so, not activating all array elements in the sub-array, otherwise, calculating the included angle between each array element in the sub-array and the target, judging whether the included angle is smaller than the activation angle or not according to the calculation result, if so, activating the subarray to finish the activation judgment of all array elements in the subarray, otherwise, and the subarray is not activated, and the included angle between each array element in the subarray and the target is continuously calculated until the activation judgment of all the array elements in the subarray is completed.
7. The method for reducing the complexity of the phased array system array element activation algorithm according to claim 5, wherein: and (4) judging whether the sub-array number n is greater than the sub-array number or not, if so, ending the program, otherwise, returning to calculate the included angle of the nth sub-array target, and ending the program until the sub-array number n is greater than the sub-array number.
8. The method for reducing the complexity of the phased array system array element activation algorithm according to claim 1, wherein: the activation calculation module sums the direction vectors pointed by the M array elements in the subarray according to the direction vector pointed by the subarray, and performs normalization processing to obtain the direction vector pointed by the subarray
Figure FDA0002430194910000021
Namely, it is
Figure FDA0002430194910000022
In the formula (I), the compound is shown in the specification,
Figure FDA0002430194910000023
is a unit vector of a connecting line of an array element i and the phased array origin, namely a direction vector pointed by the array element,
Figure FDA0002430194910000024
is composed of
Figure FDA0002430194910000025
2 norm of (d).
9. The method for reducing the complexity of the phased array system array element activation algorithm according to claim 1, wherein: activating the calculation module according to the maximum included angle thetamaxFormula for calculation
Figure FDA0002430194910000026
Calculating the maximum included angle theta between the direction of each subarray and the directions of M array elements in the subarray off linemaxAnd stored in a memory unit of the active computing module.
10. The method for reducing the complexity of the phased array system array element activation algorithm according to claim 1, wherein: for each target beam, when an activation area is calculated, firstly, subarray activation judgment is carried out; for each subarray, activating a calculation module to read the direction vector pointed by the subarray from a storage unit, then calculating the included angle between the subarray and the target,
namely, it is
Figure FDA0002430194910000027
Wherein, in the above formula
Figure FDA0002430194910000028
A direction vector representing a target pointing direction; when the activation calculation module needs to sequentially judge whether M array elements in the subarray need to be activated, the activation calculation module sequentially reads the direction vectors pointed by the M array elements from the storage unit and calculates the included angle between the array elements and the target, namely
Figure FDA0002430194910000029
Figure FDA00024301949100000210
Figure FDA0002430194910000031
And if the included angle between the array element and the target is smaller than the activation angle, the activation calculation module judges that the array element is activated, otherwise, the activation calculation module does not activate the array element.
CN202010233508.1A 2020-03-30 2020-03-30 Method for reducing complexity of phased array system array element activation algorithm Active CN111551908B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010233508.1A CN111551908B (en) 2020-03-30 2020-03-30 Method for reducing complexity of phased array system array element activation algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010233508.1A CN111551908B (en) 2020-03-30 2020-03-30 Method for reducing complexity of phased array system array element activation algorithm

Publications (2)

Publication Number Publication Date
CN111551908A true CN111551908A (en) 2020-08-18
CN111551908B CN111551908B (en) 2022-07-29

Family

ID=72003771

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010233508.1A Active CN111551908B (en) 2020-03-30 2020-03-30 Method for reducing complexity of phased array system array element activation algorithm

Country Status (1)

Country Link
CN (1) CN111551908B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112327728A (en) * 2021-01-06 2021-02-05 成都天锐星通科技有限公司 An array element operation control method, device, electronic device and storage medium
CN114447635A (en) * 2022-04-11 2022-05-06 西安星通通信科技有限公司 Method and system for improving conformal phased array antenna EIRP
CN114492615A (en) * 2022-01-21 2022-05-13 北京信息科技大学 Phased array unit control method for rapid target identification based on deep learning
CN116231335A (en) * 2021-12-06 2023-06-06 广州海格通信集团股份有限公司 Phased array antenna beam precise adjustment method and device

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6460102A (en) * 1987-08-31 1989-03-07 Mitsubishi Electric Corp Antenna system
EP0480737A2 (en) * 1990-10-10 1992-04-15 Texas Instruments Incorporated Digital radar system and method
US5600326A (en) * 1991-12-16 1997-02-04 Martin Marietta Corp. Adaptive digital beamforming architecture and algorithm for nulling mainlobe and multiple sidelobe radar jammers while preserving monopulse ratio angle estimation accuracy
JP2010032497A (en) * 2008-07-02 2010-02-12 Toshiba Corp Radar apparatus and method for forming reception beam of the same
CN102110883A (en) * 2010-12-01 2011-06-29 西安空间无线电技术研究所 Beam-forming method for forming array antenna of variable beam
CN104251989A (en) * 2014-10-09 2014-12-31 哈尔滨工程大学 Compression spatial spectrum-based single base MIMO (Multiple Input Multiple Output) radar target DOA (Direction of Arrival) estimation method
CN105353356A (en) * 2015-11-16 2016-02-24 中国电子科技集团公司第十研究所 Conformal digital beam forming system based on array slide window
CN107064888A (en) * 2017-01-13 2017-08-18 北京理工大学 A kind of method of large-scale conformal phased array antenna active region selection
CN108037692A (en) * 2017-11-22 2018-05-15 西南电子技术研究所（中国电子科技集团公司第十研究所） The beam-steering methods of extensive conformal the Digital Phased Array
CN108141256A (en) * 2016-02-04 2018-06-08 华为技术有限公司 The phase adjusting method and device of aerial array
CN108808266A (en) * 2018-06-12 2018-11-13 电子科技大学 A kind of four-dimensional antenna array combined optimization method for irregular subarray arrangement
CN108987941A (en) * 2018-05-22 2018-12-11 中国科学院国家空间科学中心 A kind of compressed sensing based one-dimensional Antenna Subarray Division
CN109344456A (en) * 2018-09-11 2019-02-15 西安空间无线电技术研究所 A Phase Compensation and Verification Method of Highly Integrated Subarray

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6460102A (en) * 1987-08-31 1989-03-07 Mitsubishi Electric Corp Antenna system
EP0480737A2 (en) * 1990-10-10 1992-04-15 Texas Instruments Incorporated Digital radar system and method
US5600326A (en) * 1991-12-16 1997-02-04 Martin Marietta Corp. Adaptive digital beamforming architecture and algorithm for nulling mainlobe and multiple sidelobe radar jammers while preserving monopulse ratio angle estimation accuracy
JP2010032497A (en) * 2008-07-02 2010-02-12 Toshiba Corp Radar apparatus and method for forming reception beam of the same
CN102110883A (en) * 2010-12-01 2011-06-29 西安空间无线电技术研究所 Beam-forming method for forming array antenna of variable beam
CN104251989A (en) * 2014-10-09 2014-12-31 哈尔滨工程大学 Compression spatial spectrum-based single base MIMO (Multiple Input Multiple Output) radar target DOA (Direction of Arrival) estimation method
CN105353356A (en) * 2015-11-16 2016-02-24 中国电子科技集团公司第十研究所 Conformal digital beam forming system based on array slide window
CN108141256A (en) * 2016-02-04 2018-06-08 华为技术有限公司 The phase adjusting method and device of aerial array
CN107064888A (en) * 2017-01-13 2017-08-18 北京理工大学 A kind of method of large-scale conformal phased array antenna active region selection
CN108037692A (en) * 2017-11-22 2018-05-15 西南电子技术研究所（中国电子科技集团公司第十研究所） The beam-steering methods of extensive conformal the Digital Phased Array
CN108987941A (en) * 2018-05-22 2018-12-11 中国科学院国家空间科学中心 A kind of compressed sensing based one-dimensional Antenna Subarray Division
CN108808266A (en) * 2018-06-12 2018-11-13 电子科技大学 A kind of four-dimensional antenna array combined optimization method for irregular subarray arrangement
CN109344456A (en) * 2018-09-11 2019-02-15 西安空间无线电技术研究所 A Phase Compensation and Verification Method of Highly Integrated Subarray

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
CANGSANG ZHAO 等: ""An 18 Mb, 12.3 GB/s CMOS pipeline-burst cache SRAM with 1.54 Gb/s/pin"", 《1999 IEEE INTERNATIONAL SOLID-STATE CIRCUITS CONFERENCE. DIGEST OF TECHNICAL PAPERS. ISSCC. FIRST EDITION (CAT. NO.99CH36278)》 *
YUNQIANG PAN 等: ""A new sparse channel estimation for 2D MIMO-OFDM systems based on compressive sensing"", 《2014 SIXTH INTERNATIONAL CONFERENCE ON WIRELESS COMMUNICATIONS AND SIGNAL PROCESSING (WCSP)》 *
李鹏凯 等: ""一种基于子阵优化的相控阵天线设计方法"", 《2019年全国天线年会论文集（中册）》 *
汪伟 等: ""双极化缝隙波导滤波天线阵"", 《2019年全国天线年会论文集（中册）》 *
王成 等: ""低副瓣相控阵面天线波束扫描及互耦特性分析"", 《雷达与对抗》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112327728A (en) * 2021-01-06 2021-02-05 成都天锐星通科技有限公司 An array element operation control method, device, electronic device and storage medium
CN116231335A (en) * 2021-12-06 2023-06-06 广州海格通信集团股份有限公司 Phased array antenna beam precise adjustment method and device
CN114492615A (en) * 2022-01-21 2022-05-13 北京信息科技大学 Phased array unit control method for rapid target identification based on deep learning
CN114447635A (en) * 2022-04-11 2022-05-06 西安星通通信科技有限公司 Method and system for improving conformal phased array antenna EIRP

Also Published As

Publication number Publication date
CN111551908B (en) 2022-07-29

Similar Documents

Publication Publication Date Title
CN111551908A (en) Method for reducing complexity of array element activation algorithm of phased array system
CN108508423B (en) Subarray digital sum-difference monopulse angle measurement method based on special-shaped array
CN108037692B (en) Beam control method of large-scale conformal digital phased array
CN104991236B (en) A kind of single base MIMO radar not rounded signal coherence source Wave arrival direction estimating method
CN109946664A (en) A kind of array radar target seeker Monopulse estimation method under major lobe suppression
CN109635240A (en) Large scale digital array null broadens Adaptive beamformer method
CN114488027A (en) Wave beam zero setting and widening method of conformal array antenna
CN110045334B (en) Low sidelobe nulling beam forming method
CN113252998A (en) Flatness optimization method for phased array antenna and difference beam signal level
CN113311397A (en) Large array rapid self-adaptive anti-interference method based on convolutional neural network
CN109188386B (en) MIMO radar high-resolution parameter estimation method based on improved two-dimensional ESPRIT algorithm
CN118364696A (en) Thin-cloth rectangular planar array antenna synthesis based on improved gray wolf optimization algorithm and phase compensation
CN111487594A (en) A Circular Array Beamforming Method Based on Particle Swarm Optimization
CN107064888B (en) A Method of Active Area Selection for Large Conformal Phased Array Antenna
CN115133291B (en) Irregular antenna subarray, phased array antenna and design method of phased array antenna
CN114384502B (en) A coherent gain processing method based on sparse array
CN116184412A (en) A Near Field Passive Target Tracking Method Based on Feedback Mechanism
CN119519783B (en) A Virtual Array Forming Method for a Mathematical Multibeam Spherical Phased Array Antenna
CN111896929B (en) DOD/DOA estimation algorithm of non-uniform MIMO radar
CN107153175B (en) Phase weighting sidelobe suppression method based on alternative projection
CN120822321A (en) A signal source localization method based on a novel non-uniform array
CN111293438A (en) Grating lobe suppression method and system based on uniform linear array broadband beam forming
CN115494458B (en) Sparse array direction diagram sidelobe suppression method based on optimized vector weighting
CN111257863A (en) High-precision multi-point linear constraint self-adaptive monopulse direction finding method
CN112164898B (en) A dual-hole coprime array antenna structure

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