Full-channel opening array antenna far field calibration method, electronic device and storage medium

Abstract

The application discloses a full-channel opening array antenna far field calibration method, electronic equipment and a storage medium, and comprises the following steps: establishing a two-dimensional plane mathematical model of an array antenna to be calibrated; acquiring coordinate data of each array unit in the array antenna model; dividing the overall directional diagram by the active unit directional diagram to obtain a first array factor; converting the first array factor to a () space to obtain a second array factor; determining sampling points according to the second array factor; calculating a weighted coefficient matrix according to the coordinate data of the sampling points and the coordinate data of the array units; and calculating array excitation information according to the weighted coefficient matrix and the second array factor to complete array backstepping calibration. The application can perform full-channel opening backstepping calibration on array antennas with non-half-wavelength, irregularly arranged array spacings, and almost has no requirement for the array factor, so that the application range is obviously expanded.

Description

Technical Field

[0001] This invention relates to the field of array antenna reverse calibration, and more particularly to a far-field calibration method, electronic device, and storage medium for a fully open array antenna. Background Technology

[0002] In active phased arrays, it is common for each antenna element to have a corresponding T / R module, meaning one channel corresponds to one antenna element. A typical planar array antenna is as follows: Figure 1 As shown. This method enables independent and precise control of each antenna element, allowing for the calibration of the array antenna using various methods such as the rotating electric vector method and the commutation measurement method. However, for an array antenna with all channels open, the switching of individual elements cannot be controlled, thus traditional calibration methods (rotating electric vector method) cannot be used for calibration. In this case, an FFT-based synthesis method can be used for back-calculation, but it has some drawbacks: Existing full-channel open array antenna calibration methods have the following shortcomings: (1) The FFT-based inverse synthesis method has strict requirements on the array spacing. For a linear array, the array factor expression is as follows: (1) In the formula, Let x represent the coordinates in the u-direction in the u, v space, which can be understood as x in f(x). express The array factor, where n represents the nth element in the array; N represents the array element count; This represents the current amplitude and phase corresponding to the nth unit; k represents the coordinate index in the u direction, i.e., the kth point in the u direction; K represents the total number of grid points in the u direction; d represents the array spacing. Indicates wavelength; The expression for the inverse Fourier transform is as follows: (2) In the formula, the inverse Fourier transform converts the data from the frequency domain to the time domain. This represents a discrete sequence in the time domain, where n is the index of the time-domain sampling point; represents a discrete sequence in the frequency domain, represents the complex amplitude of different frequency components, and k is the index of the frequency domain sampling point.

[0003] Based on equations (1) and (2), if the FFT transformation can be applied to the calculation of array antenna patterns, the following correspondence must be satisfied: (3) because and ,but This condition must be met, namely, it can only handle reverse calibration of array antennas with an array spacing of half a wavelength.

[0004] (2) The FFT-based synthesis method has strict requirements on the array layout. As mentioned in point (1) above, if you want to use the FFT synthesis method to calibrate the array antenna, the array spacing must be half a wavelength. This means that the array antenna can only be arranged in a regular pattern, and it cannot handle array back-reasoning under sparse array layout.

[0005] (3) The FFT-based comprehensive inverse reasoning method has certain requirements for the array factor. This method obtains the array excitation by performing an FFT transformation on the input array factor. The ideal range of the array factor is -1≤u≤1, -1≤v≤1 (for a one-dimensional array, it generally only refers to u in the u,v space; for a two-dimensional array, it needs to include the entire u,v two-dimensional space; it can be understood as x relative to xy, i.e., the relationship between a one-dimensional line vector and a two-dimensional plane). However, the array factor in the (u,v) space obtained in actual testing is limited to the range u^2+v^2≤1, as detailed in the appendix. Figure 2 As shown in the figure, compared with the ideal matrix factor (u,v) space range, the (u,v) space given by the actual test is reduced. This means that the calibration method has certain requirements for the matrix factor, that is, the effective information of the matrix factor should be distributed in the range of u^2+v^2≤1 as much as possible. This also limits the application scenarios of the calibration method to some extent. Summary of the Invention

[0006] The purpose of this invention is to overcome the shortcomings of the prior art and provide a far-field calibration method, electronic device and storage medium for a full-channel open array antenna.

[0007] The objective of this invention is achieved through the following technical solution: A first aspect of the present invention provides a far-field calibration method for a fully open array antenna, comprising the following steps: Establish a two-dimensional planar mathematical model of the array antenna to be calibrated, and determine the global rectangular coordinate system of the array antenna model; Obtain the coordinates of each array element in the global Cartesian coordinate system in the array antenna model. Coordinate data; Obtain the overall radiation pattern of the entire array antenna to be calibrated by actual measurement. The active element radiation pattern of the center array element of the array antenna model was obtained through simulation. ; Overall orientation map Divide by the radiation pattern of the active element Obtain the first matrix factor ,Right now ;in Indicates pitch angle, Indicates azimuth; First matrix factor Depend on( Space transformation to ( The second matrix factor is obtained in the space. ,in ; Based on the second matrix factor exist The spatial distribution is determined by selecting P sampling points and defining the location of each sampling point. spatial Coordinate data, and the second matrix factor corresponding to the sampling points. ; Based on sampling points Coordinate data and array units From the coordinate data, the weighting coefficient matrix A is calculated; Based on the weighting coefficient matrix A and the second matrix factor =b calculates the array excitation information x, and the calculation formula is Ax=b; thus completing the array back-tracking calibration.

[0008] Furthermore, the simulation method includes HFSS or CST.

[0009] Furthermore, the step of using the weighted coefficient matrix A and the second matrix factor... =b calculates the array excitation information x, including: Perform singular value decomposition on the weighted coefficient matrix A, and solve for the singular values ​​of the matrix based on the SVD method. ,in It is a left singular component; It is a singular value; for ; Singularity Singular values ​​are handled using different methods depending on the different properties of the system. ; Project the top right corner, and calculate the orthogonal projection onto the left singular vector space. , ; The scaling transformation is performed by applying the corresponding singular value processing method to the singular values ​​of each non-zero singular component according to the different properties of the system. Reconstruct the right singular vector to obtain the final solution.

[0010] Furthermore, the method of handling singular values ​​according to different properties of the system includes: For stable systems, the direct pseudo-inverse method is used; For ill-conditioned systems, regularization is used. method.

[0011] Furthermore, the ill-conditioned system refers to a system of equations that is extremely sensitive to minute changes in the input data, including the coefficient matrix and the right-hand side, leading to a significant amplification of the error in the solution; the stable system is the opposite.

[0012] A second aspect of the present invention provides an electronic device including a storage unit and a processing unit, wherein the storage unit stores computer instructions executable on the processing unit, and the processing unit executes the steps of the full-channel open array antenna far-field calibration method as described in the first aspect when executing the computer instructions.

[0013] A third aspect of the present invention provides a storage medium having computer instructions stored thereon, which, when executed, perform the steps of the full-channel open array antenna far-field calibration method as described in the first aspect.

[0014] The beneficial effects of this invention are: The calibration method in this exemplary embodiment clarifies the relationship between array factors and array excitation through mathematical derivation, and establishes a mathematical model for solving the array excitation algorithm by inverse calculation of array factors. This method breaks away from the traditional FFT synthesis concept, thus avoiding the problems of array spacing requiring half a wavelength and array arrangement requiring regular layout in traditional FFT synthesis inverse calculation algorithms.

[0015] Based on this calibration method, full-channel reverse calibration can be performed on array antennas with non-half-wavelength and irregular array spacing, and there are almost no requirements for the pairing factor, which significantly expands the application range of this method (compared to traditional calibration methods). Attached Figure Description

[0016] Figure 1 This is a planar array antenna of a general form in the prior art; Figure 2 In the existing technology Spatial diagram; Figure 3 This is a flowchart of a far-field calibration method for a full-channel open array antenna provided in an exemplary embodiment of the present invention. Detailed Implementation

[0017] The technical solution of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0018] In the description of this invention, it should be noted that the directions or positional relationships indicated by terms such as "center," "upper," "lower," "left," "right," "vertical," "horizontal," "inner," and "outer" are based on the directions or positional relationships shown in the accompanying drawings and are only for the convenience of describing this invention and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, or be constructed and operated in a specific orientation, and therefore should not be construed as a limitation of this invention. Furthermore, the terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance.

[0019] In the description of this invention, it should be noted that, unless otherwise explicitly specified and limited, the terms "installation," "connection," and "joining" should be interpreted broadly. For example, they can refer to fixed connections, detachable connections, or integral connections; they can refer to mechanical connections or electrical connections; they can refer to direct connections or indirect connections through an intermediate medium; and they can refer to the internal communication between two components. Those skilled in the art can understand the specific meaning of the above terms in this invention based on the specific circumstances.

[0020] Furthermore, the technical features involved in the different embodiments of the present invention described below can be combined with each other as long as they do not conflict with each other.

[0021] First, the implementation logic of this method will be explained, taking a two-dimensional planar array as an example: For a two-dimensional planar array with M rows and N columns, assume that the array has Q elements and P sampling points in the u, v space. Its array factor... With the position of the array element The relationships between the incentive settings are as follows: (4) In equation (4), The matrix factor is a set of matrix factors. direction and A two-dimensional function of direction; This represents the excitation (including amplitude and phase) of the m-row n-column unit. This indicates the position of the m-row, n-column cell in the x-direction; This indicates the position of the m-row, n-column cell in the x-direction; , k represents the wave number ( (), wavelength-dependent. For a fixed... and Expanding equation (4) gives: (5) In equation (5), , All of these are known quantities related to m and n. The physical meaning of is the first in the array Unit pair In space the first The weighted effect of each sampling point , To set up the array excitation, further expanding equation (4) yields: (6) The system of equations (6) can be further simplified to: (7) It can be seen from equation (7) that if If the matrix is ​​known, it can be solved to obtain the solution. (Array excitation) The specific solution formula is as follows: (8) It can be seen from equation (8) that for Any one in the matrix If the corresponding and , and Given the information, the corresponding weighting coefficients can be calculated.

[0022] See Figure 3 , Figure 3 A flowchart of a far-field calibration method for a full-channel open array antenna provided in an exemplary embodiment of the present invention is shown, including the following steps: S1: Establish a two-dimensional planar mathematical model of the array antenna to be calibrated, and determine the global rectangular coordinate system of the array antenna model.

[0023] In one specific exemplary embodiment, the array antenna element size is 64 (Q=64), the array spacing is 15mm, and the operating frequency is 8GHz-12GHz (i.e., the ratio of array spacing to wavelength is 0.4-0.6). A global Cartesian coordinate system is established with the geometric center of the array antenna in the two-dimensional plane as the origin of the global coordinate system. , For the direction of dissemination, The plane is the plane on which the array is located.

[0024] S2: Obtain the coordinates of each array element in the global Cartesian coordinate system of the array antenna model. Coordinate data.

[0025] In a specific exemplary embodiment, each array element is abstracted as a point mass on a plane, and the position of each point mass is determined by the position of the physical geometric center of each element in the array in the global coordinate system. For each element q, its position in the global coordinate system can be obtained. Representation of position information on a plane .

[0026] S3: Obtain the overall radiation pattern of the entire array antenna to be calibrated from actual measurements. The active element radiation pattern of the center array element of the array antenna model was obtained through simulation. ; Overall orientation map Divide by the radiation pattern of the active element Obtain the first matrix factor ,Right now in Indicates pitch angle, Indicates the azimuth angle.

[0027] In one specific exemplary embodiment, the overall radiation pattern of the array antenna is obtained by actual measurement in an anechoic chamber. The radiation pattern of the active element of a certain central element of the array is obtained through HFSS or CST simulation. Finally, it can be determined by... Calculated First matrix factor in space .

[0028] S4: Factor the first matrix Depend on( Space transformation to ( The second matrix factor is obtained in the space. ,in .

[0029] S5: Based on the second matrix factor exist The spatial distribution is selected using a number of sampling points, P. And determine the sampling points at spatial Coordinate information of coordinate data and corresponding , Q is the number of array elements in the antenna array to be calibrated; S6: Based on sampling points Coordinate data and array units From the coordinate data, the weighting coefficient matrix A is calculated; The expression is as follows: in, The solution formula is as follows: in, Indicates wave number, The antenna operating frequency, The speed of light; and They represent the nth in the array Each unit in direction and Coordinates of direction; and They represent In space the first Each sampling point is direction and The coordinates of the direction, where k represents the wave number. The physical meaning of is the first in the array Unit pair In space the first The weighting effect of each sampling point.

[0030] S7: Based on the weighted coefficient matrix A and the second matrix factor =b calculates the array excitation information x, with the formula Ax=b; thus completing the array back-reasoning calibration. A system of relevant equations is established based on the principle of array factor solving. The system of equations is expressed as follows: Wherein, the coefficient matrix Given the matrix factor matrix Given the array excitation information matrix The solution to x can then be obtained by solving the system of equations.

[0031] In summary, the calibration method in this exemplary embodiment clarifies the relationship between array factors and array excitation through mathematical derivation, and establishes a mathematical model for solving the array excitation algorithm by inversely deducing the array factors. This method breaks away from the traditional FFT synthesis approach, thus avoiding the problems encountered in traditional FFT synthesis inverse algorithms, such as the array spacing having to be half a wavelength and the array arrangement having to be regular.

[0032] Based on this calibration method, full-channel reverse calibration can be performed on array antennas with non-half-wavelength and irregular array spacing, and there are almost no requirements for the pairing factor, which significantly expands the application range of this method (compared to traditional calibration methods).

[0033] More preferably, in an exemplary embodiment, the step of basing the weighted coefficient matrix A on the second matrix factor... =b calculates the array excitation information x, including: S8: Perform singular value decomposition on the weighted coefficient matrix A, and solve for the singular values ​​of the matrix based on the SVD method. ,;in for An orthogonal matrix representing the left singular component; for Diagonal array, elements , represents a singular value; for An orthogonal matrix represents the right singular component; The values ​​of each element can be obtained by solving the corresponding functions in MATLAB.

[0034] S9: Singular Values Singular values ​​are handled using different methods depending on the different properties of the system. In a specific exemplary embodiment, the method of processing singular values ​​according to different properties of the system includes: For stable systems, the direct pseudo-inverse method is used: For ill-conditioned systems, regularization is used. method: In the formula, Indicates wavelength.

[0035] It should be noted that the ill-conditioned system refers to a system of equations that is extremely sensitive to small changes in the input data, including the coefficient matrix and the right-hand side, leading to a significant amplification of the error in the solution; the stable system is the opposite. Generally, in mathematics, it can be directly judged by the norm condition number of the coefficient matrix A; the larger the value, the more severe the ill-conditioned system. However, in practical engineering applications, whether a system is ill-conditioned is determined by a preset standard threshold (i.e., if the condition number is greater than the preset standard threshold, it is considered an ill-conditioned system).

[0036] S10: Project the top right end and calculate the orthogonal projection of the left singular vector space; In step S10, the vector c is obtained, laying the foundation for solving for vector d in step 11. Here, b corresponds to the b in Ax=b. It is obtained from the singular value decomposition of the matrix in step S8; c is an intermediate variable obtained after projection, which is convenient for subsequent use.

[0037] It should be noted that for this step: First, performing singular value decomposition on matrix A yields: ; Wherein, the left singular matrix U is an orthogonal matrix, which is an orthogonal basis of the output space; the right singular matrix V is an orthogonal basis of the input space.

[0038] Based on this, the top right end of the projection... .

[0039] (Therefore, the left singular vector is introduced.) Its mathematical essence is to transform the problem into the optimal coordinate system, thereby decoupling the original system of equations, making it easier to solve and regularize. Here, b is the same as b in Ax=b.

[0040] S11: Descaling transformation: Based on the different properties of the system and the corresponding singular value processing method, scale transformation is performed on the singular values ​​of each non-zero singular component. The essence of this step is to suppress the negative impact of small singular values ​​on the solution by modifying the inverse of the singular values, thereby achieving the separation of the physical dominant mode from the noise mode.

[0041] S12: Reconstruct the right singular vector to obtain the final solution.

[0042] In mathematical logic, step 12, x is obtained by multiplying the vector V and the vector d. The vector V is obtained through the singular value decomposition of the matrix in step S8; the vector d is obtained by solving the scaling transformation in step S11.

[0043] It should be noted that: First, the coordinates c of the right-hand top b in the left singular space U are obtained through S10; second, d is obtained by solving for the right projected signal c through S11, the physical essence of which is to realize the transformation from signal space to source space, and the numerical function can control the propagation of the ill-conditioned system through scaling; finally, the final solution is obtained by solving for d and V through S12. Here, V is the right singular matrix (obtained from S8); d is the scaled coefficient vector (obtained from S11). Having the same inventive concept as the exemplary embodiments described above, another exemplary embodiment of the present invention provides an electronic device, including a storage unit and a processing unit, wherein the storage unit stores computer instructions executable on the processing unit, and the processing unit executes the computer instructions as follows: Figure 1 The steps of the described method for far-field calibration of a full-channel open array antenna.

[0044] Electronic devices are manifested in the form of general-purpose computing devices. Components of electronic devices may include, but are not limited to: at least one processing unit, at least one storage unit, and a bus connecting different system components (including storage units and processing units).

[0045] The storage unit stores program code that can be executed by the processing unit, causing the processing unit to perform the steps described in the "Exemplary Methods" section above, based on various exemplary embodiments of the present invention. For example, the processing unit can perform actions such as... Figure 1 The method shown in the figure.

[0046] The storage unit may include readable media in the form of volatile storage units, such as random access memory (RAM) and / or cache storage units, and may further include read-only memory (ROM).

[0047] The storage unit may also include a program / utility having a set (at least one) of program modules, including but not limited to: an operating system, one or more application programs, other program modules, and program data, each or some combination of these examples may include an implementation of a network environment.

[0048] A bus can represent one or more of several types of bus structures, including a memory cell bus or memory cell controller, a peripheral bus, a graphics acceleration port, a processing unit, or a local bus that uses any of the various bus structures.

[0049] The electronic device can also communicate with one or more external devices (e.g., keyboards, pointing devices, Bluetooth devices, etc.), one or more devices that enable a user to interact with the electronic device, and / or any device that enables the electronic device to communicate with one or more other computing devices (e.g., routers, modems, etc.). This communication can be achieved through input / output (I / O) interfaces. Furthermore, the electronic device can communicate with one or more networks (e.g., local area networks (LANs), wide area networks (WANs), and / or public networks, such as the Internet) via a network adapter. The network adapter communicates with other modules of the electronic device via a bus. It should be understood that other hardware and / or software modules can be used in conjunction with the electronic device, including but not limited to: microcode, device drivers, redundant processing units, external disk drive arrays, RAID systems, tape drives, and data backup storage systems.

[0050] Through the above description, those skilled in the art will readily understand that the exemplary embodiments described herein can be implemented by software or by combining software with necessary hardware. Therefore, the technical solution according to this exemplary embodiment can be embodied in the form of a software product, which can be stored in a non-volatile storage medium (such as a CD-ROM, USB flash drive, external hard drive, etc.) or on a network, including several instructions to cause a computing device (such as a personal computer, server, terminal device, or network device, etc.) to execute the method according to this exemplary embodiment.

[0051] Having the same inventive concept as the exemplary embodiments described above, another exemplary embodiment of the present invention provides a storage medium storing computer instructions that, when executed, perform the steps of the full-channel open array antenna far-field calibration method as described in the first aspect.

[0052] Based on this understanding, the technical solution of this embodiment, or the part that contributes to the prior art, or part of the technical solution, can be embodied in the form of a software product (program product). The computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present invention.

[0053] The program product may employ any combination of one or more readable media. A readable medium may be a readable signal medium or a readable storage medium. A readable storage medium may be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination thereof. More specific examples of readable storage media (a non-exhaustive list) include: an electrical connection having one or more wires, a portable disk, a hard disk, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fiber, portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination thereof.

[0054] Computer-readable signal media may include data signals propagated in baseband or as part of a carrier wave, carrying readable program code. Such propagated data signals may take various forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination thereof. A readable signal medium may also be any readable medium other than a readable storage medium, capable of sending, propagating, or transmitting programs for use by or in conjunction with an instruction execution system, apparatus, or device.

[0055] The program code contained on the readable medium may be transmitted using any suitable medium, including but not limited to wireless, wired, optical fiber, RF, etc., or any suitable combination thereof.

[0056] Program code for performing the operations of this invention can be written in any combination of one or more programming languages, including object-oriented programming languages ​​such as Java and C++, and conventional procedural programming languages ​​such as C or similar languages. The program code can execute entirely on the user's computing device, partially on the user's device, as a standalone software package, partially on the user's computing device and partially on a remote computing device, or entirely on a remote computing device or server. In cases involving remote computing devices, the remote computing device can be connected to the user's computing device via any type of network, including a local area network (LAN) or a wide area network (WAN), or it can be connected to an external computing device (e.g., via the Internet using an Internet service provider).

[0057] Obviously, the above embodiments are merely illustrative examples for clear explanation and are not intended to limit the implementation. Those skilled in the art can make other variations or modifications based on the above description. It is neither necessary nor possible to exhaustively list all possible implementations. However, obvious variations or modifications derived therefrom are still within the scope of protection of this invention.

Claims

1. A far-field calibration method for a fully open array antenna, characterized in that: Includes the following steps: Establish a two-dimensional planar mathematical model of the array antenna to be calibrated, and determine the global rectangular coordinate system of the array antenna model; Obtain the coordinates of each array element in the global Cartesian coordinate system in the array antenna model. Coordinate data; Obtain the overall radiation pattern of the entire array antenna to be calibrated by actual measurement. The active element radiation pattern of the central array element of the array antenna model was obtained through simulation. ; Overall orientation map Divide by the radiation pattern of the active element Obtain the first matrix factor ,Right now ;in Indicates pitch angle, Indicates azimuth; First matrix factor Depend on( Space transformation to ( The second matrix factor is obtained in the space. ,in ; Based on the second matrix factor exist The spatial distribution is determined by selecting P sampling points and defining the location of each sampling point. spatial Coordinate data, and the second matrix factor corresponding to the sampling points. ; Based on sampling points Coordinate data and array units From the coordinate data, the weighting coefficient matrix A is calculated; Based on the weighting coefficient matrix A and the second matrix factor =b calculates the array excitation information x, and the calculation formula is Ax=b; thus completing the array back-tracking calibration. Among them, the basis of sampling points Coordinate data and array units The coordinate data is used to calculate the weighted coefficient matrix A. The expression is as follows: ； in, The solution formula is as follows: ； in, Indicates wave number, The antenna operating frequency, The speed of light; and They represent the nth in the array Each unit in direction and Coordinates of direction; and They represent In space the first Each sampling point is direction and The coordinates of direction, where k represents the wave number; The physical meaning of is the first in the array Unit pair In space the first The weighting effect of each sampling point; Q is the number of array elements in the array antenna to be calibrated.

2. The far-field calibration method for a full-channel open array antenna according to claim 1, characterized in that: The simulation methods include HFSS or CST.

3. The far-field calibration method for a full-channel open array antenna according to claim 1, characterized in that: The weighting coefficient matrix A and the second matrix factor are used as follows: =b calculates the array excitation information x, including: Perform singular value decomposition on the weighted coefficient matrix A, and solve for the singular values ​​of the matrix based on the SVD method. ,in It is a left singular component; It is a singular value; for ; Singularity Singular values ​​are handled using different methods depending on the different properties of the system. ; Project the right-hand side and calculate the orthogonal projection onto the left singular vector space. , ; The scaling transformation is performed by applying the corresponding singular value processing method to the singular values ​​of each non-zero singular component according to the different properties of the system. Reconstruct the right singular vector to obtain the final solution.

4. The far-field calibration method for a full-channel open array antenna according to claim 3, characterized in that: The method of handling singular values ​​according to different properties of the system includes: For stable systems, the direct pseudo-inverse method is used; For ill-conditioned systems, regularization is used. method.

5. The far-field calibration method for a full-channel open array antenna according to claim 4, characterized in that: The ill-conditioned system refers to a system of equations that is extremely sensitive to minute changes in the input data, including the coefficient matrix and right-hand side terms, which leads to a significant amplification of the error in the solution; the stable system is the opposite.

6. An electronic device comprising a storage unit and a processing unit, wherein the storage unit stores computer instructions executable on the processing unit, characterized in that: When the processing unit executes the computer instructions, it performs the steps of the far-field calibration method for a full-channel open array antenna as described in any one of claims 1 to 5.

7. A storage medium storing computer instructions thereon, characterized in that: When the computer instructions are executed, they perform the steps of the far-field calibration method for a full-channel open array antenna as described in any one of claims 1 to 5.

Images