FPGA-based rectangular planar phased array antenna rapid phase matching calculation system and method
Technical Field
The invention relates to the technical field of phased array antenna beam pointing control, in particular to an efficient phase matching solution scheme based on an FPGA, which is suitable for a large-scale rectangular planar phased array antenna in high dynamic application scenes such as low-orbit satellite communication, communication in motion and the like.
Background
With the development of low-orbit satellite communications and 5G/6G networks, terminals are placing demands on phased array antennas for miniaturization, low cost, fast beam switching (in the order of microseconds). The traditional phase matching resolving method based on Digital Signal Processing (DSP) or software is faced with the problems of large calculation amount, high delay, large resource occupation and the like in a large-scale array (more than 1024 array elements, for example), and is difficult to meet the requirement of dynamic satellite tracking.
In the prior art, phase matching calculation needs to calculate the phase difference for each array element independently and perform the residual operation, so that a large number of multipliers and storage units are consumed when the FPGA is realized. For example, a conventional scheme may require N multipliers (N is the number of elements) and trigonometric function look-up table resources increase linearly with array size.
Disclosure of Invention
Aiming at the technical problems, the invention aims to provide an efficient phase matching resolving system and algorithm based on an FPGA, which can obviously reduce hardware resource consumption while guaranteeing microsecond computing speed, and is suitable for a large-scale rectangular plane phased array antenna system with low cost and low power consumption.
In order to achieve the aim of the invention, the technical scheme adopted by the invention is that the rectangular plane phased array antenna rapid phase matching and resolving device based on the FPGA comprises:
input module for receiving beam pointing pitch angle And azimuth angle;
Trigonometric function table look-up module, which adopts symmetrical compression storage strategy to store onlySine and cosine values in the enclosure are expanded to a full period through address mapping, and the method is used for quickly acquiring trigonometric function parameters required by phase calculation;
the row-column decomposition calculation module is used for decomposing the two-dimensional rectangular plane phased array into one-dimensional linear arrays in the directions of an x axis and a y axis, and calculating the phase difference between adjacent array elements according to the pitch angle and the azimuth angle of beam pointing;
The phase difference between adjacent array elements in the x-axis direction is as follows:
(1)
adjacent array element in y-axis direction the phase difference between them is:
(2)
Wherein the method comprises the steps of ,Is the distance between the adjacent array elements,For transmitting signal carrier wavelengths.
The phase accumulation module is used for accumulating adjacent phase differences obtained by column-row decomposition to generate row/column phase vectors and storing the row/column phase vectors, wherein the accumulation operation is expressed as multiplication operation in a formula, and the phase difference between an m-th array element and a 0-coordinate array element in the row vectors is expressed as follows:
(3)
The phase difference between the nth element and the 0-coordinate element in the column vector is expressed as:
(4)
The array element phase calculation module is used for carrying out combination operation according to the row/column phase vector obtained by the phase accumulation module and the formula (5) to obtain a phase value of each array element in the array;
(5)
the residual calculating module converts the traditional residual calculating operation into normalized numerical residual calculating processing of the numerical value 1 so as to simplify the calculation time and the resource consumption of the FPGA.
The traditional calculation phase formula is:
(6)
in calculating the phase according to the nature of the remainder Modulo with subsequent remainder operationThe formula is rewritten, however:
(7)
and the quantization output module is used for carrying out quantization processing on the phase after the remainder according to the bit width of the digital phase shifter and generating an array element phase matching control code.
As a preferred mode of the invention, the row-column decomposition calculation module further comprises a parallel pipeline architecture, which respectively and independently processes phase difference accumulation operation in the x-axis direction and the y-axis direction, and stores intermediate results through a dual-port BRAM to support subsequent array element phase synthesis.
As a preferred mode of the invention, the specific implementation process of the residual module comprises the steps of converting absolute phase values in a formula into normalized coefficients and determining final phase values by calculating a remainder of 1.
As an optimal mode of the invention, the trigonometric function table look-up module adopts 16bit quantization precision, wherein the high order is a sign bit, and the parallel resolving of the phases of multiple array elements is realized through time sequence control.
The invention also provides a resolving method based on the rapid phase matching resolving device, which comprises the following steps:
step one, will 、The input trigonometric function table lookup module obtains the normalized trigonometric function value、And;
Step two, calculating a row direction normalization coefficient according to a formula (10) and a formula (11):
(10)
(11)
Wherein the method comprises the steps of 、The distance between array elements in the row-column direction, c is the speed of light,Is the carrier frequency;
Step three, calculating the phase difference between adjacent array elements in the row/column direction according to the results of the step one and the step two through a formula (12) and a formula (13):
(12)
(13)
Generating a basic phase sequence through a row/column accumulator according to the result of the step three and then according to a formula (14) and a formula (15), and storing the basic phase sequence;
(14)
(15)
Step five, according to the result of step four, carry on the phase synthesis of phase matching to each array element (m, n) according to the formula (16);
(16)
Step six, carrying out remainder of 1 on the synthesized phase to obtain 0-2 pi internal phase;
And step seven, quantizing the phase according to the quantization precision of the phase shifter to obtain a phase shifter control code.
The invention also provides a low-orbit satellite communication terminal antenna system which comprises the FPGA-based rectangular planar phased array antenna rapid phase matching and resolving device, and the configuration meets the technical requirement that the beam pointing switching period is less than or equal to 100ms and the gain is more than or equal to 25 dBi.
Compared with the prior art, the technical scheme provided by the invention has the following technical effects:
1. The phase matching resolving system provided by the invention supports the rectangular plane phased array with the number of the adaptive array elements being more than or equal to 1024 and the single phase matching resolving time being less than or equal to 20 mu s under the condition that the clock frequency is more than or equal to 100 MHz.
2. The invention realizes the hardware optimization of the phase matching resolving system, and utilizes the trigonometric function property to realize the followingIs mapped to the trigonometric function value of (2)And saving 87.5% of storage resources.
3. The invention decomposes phase calculation into the phases of table lookup, multiplication, accumulation, array element synthesis, remainder calculation, quantization and the like, and improves throughput rate.
4. The invention is particularly suitable for the fast beam pointing switching requirement of the low-orbit satellite communication terminal, and meets the requirement of high real-time performance while reducing the hardware cost.
Drawings
Fig. 1 is a block diagram of a rectangular planar phased array antenna phase matching calculation system based on an FPGA in this embodiment.
Fig. 2 is a flow chart of a rectangular planar phased array antenna phase matching calculation method based on an FPGA in this embodiment.
Detailed Description
For a better understanding of the invention, reference will now be made in detail to the drawings and to the accompanying examples.
As shown in fig. 1, a block diagram of a rectangular planar phased array antenna phase matching resolving system based on an FPGA specifically includes:
input module for receiving beam pointing pitch angle And azimuth angle;
Trigonometric function table look-up module, which adopts symmetrical compression storage strategy to store onlySine and cosine values in the range are expanded to a full period through address mapping, and the method is used for quickly acquiring trigonometric function parameters required by phase calculation;
the row-column decomposition calculation module is used for decomposing the two-dimensional rectangular plane phased array into one-dimensional linear arrays in the directions of an x axis and a y axis, and calculating the phase difference between adjacent array elements according to the pitch angle and the azimuth angle of beam pointing;
The phase difference between adjacent array elements in the x-axis direction is as follows:
(1)
adjacent array element in y-axis direction the phase difference between them is:
(2)
Wherein the method comprises the steps of ,Is the distance between the adjacent array elements,For transmitting signal carrier wavelengths.
The phase accumulation module is used for accumulating adjacent phase differences obtained by column-row decomposition to generate row/column phase vectors and storing the row/column phase vectors, wherein the accumulation operation is expressed as multiplication operation in a formula, and the phase difference between an m-th array element and a 0-coordinate array element in the row vectors is expressed as follows:
(3)
The phase difference between the nth element and the 0-coordinate element in the column vector is expressed as:
(4)
The array element phase calculation module is used for carrying out combination operation according to the row/column phase vector obtained by the phase accumulation module and the formula (5) to obtain a phase value of each array element in the array;
(5)
the residual calculating module converts the traditional residual calculating operation into normalized numerical residual calculating processing of the numerical value 1 so as to simplify the calculation time and the resource consumption of the FPGA.
The traditional calculation phase formula is:
(6)
in calculating the phase according to the nature of the remainder Modulo with subsequent remainder operationThe formula is rewritten, however:
(7)
and the quantization output module is used for carrying out quantization processing on the phase after the remainder according to the bit width of the digital phase shifter and generating an array element phase matching control code.
As shown in fig. 2, a flow chart of a rectangular planar phased array antenna phase matching resolving method based on an FPGA specifically includes the following steps:
step one, will 、The input trigonometric function table lookup module obtains the normalized trigonometric function value、And;
Step two, calculating a row direction normalization coefficient according to a formula (10) and a formula (11):
(10)
(11)
Wherein the method comprises the steps of 、The distance between array elements in the row-column direction, c is the speed of light,Is the carrier frequency;
Step three, calculating the phase difference between adjacent array elements in the row/column direction according to the results of the step one and the step two through a formula (12) and a formula (13):
(12)
(13)
Generating a basic phase sequence through a row/column accumulator according to the result of the step three and then according to a formula (14) and a formula (15), and storing the basic phase sequence;
(14)
(15)
Step five, according to the result of step four, carry on the phase synthesis of phase matching to each array element (m, n) according to the formula (16);
(16)
Step six, carrying out remainder of 1 on the synthesized phase to obtain 0-2 pi internal phase;
And step seven, quantizing the phase according to the quantization precision of the phase shifter to obtain a phase shifter control code.
The actual measurement data shows that when the Xilinx Zynq-7045 chip is adopted, only 4 DSP multipliers and 70.5BRAM resources are occupied, the resource amount is small, and the Xilinx Zynq-7045 chip can be transplanted to a low-end FPGA (such as Artix series).
The speed is about 12 for completing 1024 array element phase matching calculation by 100MHz clock。
The expansibility is that the expansion of the array scale is supported by adjusting the bit width of the accumulator, and the expansion can be extended to other polygonal arrays by discarding the calculated value of a specific array element.
The present invention is not limited to the above-mentioned embodiments, and any changes or substitutions that are easily contemplated by those skilled in the art within the scope of the present invention are intended to be included in the scope of the present invention. Therefore, the protection scope of the present invention should be subject to the protection scope of the claims.