A radar target monostatic RCS fast squint method

Abstract

The application discloses a radar target single-station RCS fast scanning angle method, which comprises the following steps: firstly, a triangle grid is used to divide a created radar target geometric model; secondly, an ACA algorithm is used to perform low-rank compression processing on an excitation matrix; thirdly, a double-layer mixed compression method combining an ACA algorithm and QR decomposition is used; fourthly, QR decomposition is performed on the compressed matrix, and ACA algorithm is used for re-compression; fifthly, the number of excitations required for iterative solving in the solving process is further reduced; and finally, a multilevel fast multipole method is used to combine a near-zone block diagonal preconditioning technology to construct a matrix equation, so that the number of iterative solving times of the matrix equation is greatly reduced, the solving of the matrix equation is accelerated, and the solving efficiency of the radar target single-station RCS under the excitation of a plane wave at multiple angles is improved.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of radar data processing, and particularly relates to a radar target single-station RCS fast scanning angle method. BACKGROUND

[0002] With the rapid development of computer technology and electronic science and technology, numerical calculation plays a crucial role in the research of many practical engineering problems. In the field of electromagnetism, for example, large antenna design and simulation, radar target stealth performance analysis, signal integrity and electromagnetic compatibility analysis and many other practical engineering applications require efficient establishment of three-dimensional geometric models and electromagnetic field numerical simulation analysis.

[0003] Since the middle of the twentieth century, efficient and accurate numerical calculation methods have emerged, and many scholars have applied them to electromagnetic field simulation analysis. Electromagnetic field numerical simulation software has also developed rapidly. For example, High Frequency Structure Simulator (HFSS), Feldberechnung bei Korpern mit beliebiger Oberflache (FEKO), and Computer Simulation Technology Studio Suite (CST) have become important tools in electromagnetic aided design and analysis. The development of such simulation software has also promoted the rapid development of many fields such as electronics, communications, radar, and remote sensing. It also plays a crucial role in cross-disciplines such as optics, quantum, and biology. Although existing electromagnetic simulation software is powerful and can solve many common scientific and practical engineering electromagnetic simulation problems, and with the increasing complexity of large-scale electromagnetic scattering models, it is still extremely challenging to quickly calculate the wide-angle domain electromagnetic scattering characteristics of radar targets under the premise of ensuring solution accuracy. For the problem of wide-angle domain, multi-angle radar target electromagnetic scattering calculation, due to the large electrical size of the radar scattering target, the scattering cross section (RCS) of the single-station radar scattering target at a single angle of incidence is often solved using an iterative method based on the traditional Method of Moments (MOM).

[0004] However, in practical engineering applications, the conductor target is usually of large electrical size or complex structure, and it is often necessary to solve the radar scattering target at multiple angles of incidence. For this type of electromagnetic scattering problem, the number of unknowns in the solution process using the Method of Moments is large, and the memory occupied and the time required for solving the matrix equation obtained by the Method of Moments are long, which poses a great challenge to computer resources. SUMMARY

[0005] In order to solve the above problems existing in the prior art, the present application provides a radar target single station RCS fast scan angle method. The technical problem to be solved by the present application is realized through the following technical scheme:

[0006] The radar target single station RCS fast scan angle method provided by the present application comprises:

[0007] creating a radar target geometric model; wherein, materials are set for the radar target geometric model, and plane wave excitations of multiple angles are set; an excitation matrix is constructed according to the plane wave excitations of multiple angles;

[0008] triangular mesh division is performed on the radar target geometric model to generate a plurality of triangular meshes, and an impedance matrix is constructed according to the plurality of triangular meshes;

[0009] ACA algorithm is used to compressively decompose the excitation matrix to obtain a first full-rank matrix and a second full-rank matrix;

[0010] QR decomposition is performed on the first full-rank matrix to obtain a first triangular matrix and a first orthogonal matrix, and QR decomposition and Hermite matrix operation are performed on the second full-rank matrix to obtain a second triangular matrix and a second orthogonal matrix;

[0011] matrix multiplication operation is performed on the first triangular matrix and the second triangular matrix to obtain a first matrix;

[0012] ACA algorithm is used to compressively decompose the first matrix again to obtain a third full-rank matrix and a fourth full-rank matrix;

[0013] matrix multiplication operation is performed on the third full-rank matrix and the first orthogonal matrix to obtain a second matrix, and matrix multiplication operation is performed on the fourth full-rank matrix and the second orthogonal matrix to obtain a third matrix;

[0014] based on the multilevel fast multipole method, the second matrix is taken as a new excitation matrix, a matrix equation is constructed according to the impedance matrix and the new excitation matrix in combination with a near-zone block diagonal preconditioning technique, and a generalized minimum residual method is used to solve the matrix equation to obtain a first current matrix;

[0015] matrix multiplication operation is performed on the first current matrix and the third matrix to obtain a second current matrix, and radar target single station RCS is calculated according to the second current matrix.

[0016] In an embodiment of the present application, ACA algorithm is used to compressively decompose the excitation matrix, and the formula is represented as:

[0017] V N×M≈A N×r ·B r×M ;

[0018] wherein, V N×M represents an excitation matrix constructed under a plurality of angle plane wave excitations, A N×r represents a first full rank matrix, B r×M represents a second full rank matrix, N represents the number of unknowns, M represents the number of excitation vectors, r represents the effective rank number of the excitation matrix, and r < min(N, M).

[0019] In an embodiment of the present application, QR decomposition is performed on the first full rank matrix to obtain a first triangular matrix and a first orthogonal matrix, and the formula is represented as:

[0020] A N×r = AQ N×r · AR r×r ;

[0021] wherein, A N×r represents the first full rank matrix, AQ N×r represents the first orthogonal matrix, and AR r×r represents the first triangular matrix.

[0022] In an embodiment of the present application, QR decomposition and Hermitian matrix operation are performed on the second full rank matrix to obtain a second triangular matrix and a second orthogonal matrix, and the formula is represented as:

[0023]

[0024] wherein, represents the Hermitian matrix of the second full rank matrix, BR r×r represents the second triangular matrix, BQ r×M represents the second orthogonal matrix, and (·) H represents Hermitian matrix operation.

[0025] In an embodiment of the present application, matrix multiplication operation is performed on the first triangular matrix and the second triangular matrix to obtain a first matrix, and the formula is represented as:

[0026] ABR r×r = AR r×r · BR r×r ;

[0027] wherein, ABR r×r represents the first matrix, AR r×r represents the first triangular matrix, and BR r×r represents the second triangular matrix.

[0028] In an embodiment of the present application, ACA algorithm is used to perform re-compression decomposition on the first matrix, and the formula is represented as:

[0029] ABR r×r =A1 r×k ·B1 k×r ;

[0030] wherein, ABR r×r represents the first matrix, A1 r×k represents the third full-rank matrix, B1 k×r represents the fourth full-rank matrix, k represents the effective rank number of the first matrix, and k < r.

[0031] In an embodiment of the present application, the third full-rank matrix and the first orthogonal matrix are subjected to matrix multiplication operation to obtain a second matrix, which is represented by the formula:

[0032]

[0033] wherein, represents the second matrix, A1 r×k represents the third full-rank matrix, AQ N×r represents the first orthogonal matrix.

[0034] In an embodiment of the present application, the fourth full-rank matrix and the second orthogonal matrix are subjected to matrix multiplication operation to obtain a third matrix, which is represented by the formula:

[0035]

[0036] wherein, represents the third matrix, B1 k×r represents the fourth full-rank matrix, BQ r×M represents the second orthogonal matrix.

[0037] In an embodiment of the present application, based on the multilayer fast multipole method, the second matrix is taken as a new excitation matrix, and a matrix equation is constructed according to the impedance matrix and the new excitation matrix in combination with the near-zone block diagonal preconditioning technique, which is represented by the formula:

[0038]

[0039] wherein, Z N×N represents the impedance matrix, M N×N represents a preconditioning matrix obtained by the near-zone block diagonal preconditioning technique, (·) -1 represents an inverse operation, I N×k represents a first current matrix, represents the second matrix.

[0040] In an embodiment of the present application, the first current matrix and the third matrix are subjected to matrix multiplication operation to obtain a second current matrix, which is represented by the formula:

[0041]

[0042] wherein I N×M denotes a second current matrix, I N×k denotes a first current matrix, denotes a third matrix.

[0043] The present application has the following beneficial effects:

[0044] The radar target single-station RCS fast scanning angle method provided by the present application solves the problem of too long time required by the existing multi-layer fast multipole method for solving multiple angles of plane wave excitation, is a new idea for realizing radar target single-station RCS fast scanning, and the specific implementation process comprises: creating a radar target geometric model; wherein, materials are set for the radar target geometric model, and multiple angles of plane wave excitation are set; an excitation matrix is constructed according to the multiple angles of plane wave excitation; the radar target geometric model is triangularly meshed to generate a plurality of triangular meshes, and an impedance matrix is constructed according to the plurality of triangular meshes; ACA algorithm is used to compress and decompose the excitation matrix to obtain a first full-rank matrix and a second full-rank matrix; QR decomposition is performed on the first full-rank matrix to obtain a first triangular matrix and a first orthogonal matrix, and QR decomposition and Hermite matrix operation are performed on the second full-rank matrix to obtain a second triangular matrix and a second orthogonal matrix; matrix multiplication is performed on the first triangular matrix and the second triangular matrix to obtain a first matrix; ACA algorithm is used to compress and decompose the first matrix again to obtain a third full-rank matrix and a fourth full-rank matrix; matrix multiplication is performed on the third full-rank matrix and the first orthogonal matrix to obtain a second matrix, and matrix multiplication is performed on the fourth full-rank matrix and the second orthogonal matrix to obtain a third matrix; based on the multi-layer fast multipole method, the second matrix is taken as a new excitation matrix, a matrix equation is constructed according to the impedance matrix and the new excitation matrix in combination with the near-zone block diagonal preconditioning technique, and the generalized minimum residual method is used to solve the matrix equation to obtain a first current matrix; matrix multiplication is performed on the first current matrix and the third matrix to obtain a second current matrix, and radar target single-station RCS is calculated according to the second current matrix. It can be seen that, firstly, the created radar target geometric model is triangularly meshed, then ACA algorithm is used to compress and process the excitation matrix at a low rank, and then a double-layer mixed compression method combining ACA algorithm and QR decomposition is used to compress the compressed matrix by QR decomposition and ACA algorithm, further reducing the number of excitations required for iterative solution in the solving process, and finally the multi-layer fast multipole method is used to construct a matrix equation in combination with the near-zone block diagonal preconditioning technique, which greatly reduces the number of iterations of the matrix equation, speeds up the solution of the matrix equation, and improves the solving efficiency of radar target single-station RCS under multiple angles of plane wave excitation.

[0045] The application will be further described in detail below with reference to the accompanying drawings and embodiments. BRIEF DESCRIPTION OF DRAWINGS

[0046] Figure 1 is a flowchart of a radar target monostatic RCS fast scan angle method provided by an embodiment of the application;

[0047] Figure 2 is a radar target geometric model schematic diagram provided by an embodiment of the application;

[0048] Figure 3 is a schematic diagram of setting materials, multiple angle plane wave excitations on a radar target geometric model provided by an embodiment of the application;

[0049] Figure 4 is a triangular row grid schematic diagram after triangular meshing of a radar target geometric model provided by an embodiment of the application;

[0050] Figure 5 is a network topology structure schematic diagram corresponding to a multilayer fast multipole method based on the method of moments provided by an embodiment of the application;

[0051] Figure 6 is an implementation flowchart of a QR decomposition combined with an ACA hybrid compression acceleration algorithm provided by an embodiment of the application;

[0052] Figure 7 is a radar target monostatic RCS comparison result schematic diagram of an existing multilayer fast multipole method and the method provided by an embodiment of the application;

[0053] Figure 8 is a structure schematic diagram of an electronic device provided by an embodiment of the application. DETAILED DESCRIPTION

[0054] The application will be further described in detail below with reference to the accompanying drawings and embodiments, but the embodiments of the application are not limited thereto.

[0055] In order to improve the ability to process large size and complex target scattering characteristics, please refer to Figure 1 An embodiment of the application provides a radar target monostatic RCS fast scan angle method, specifically including the following steps:

[0056] S10, creating a radar target geometric model; wherein, materials are set on the radar target geometric model, and multiple angle plane wave excitations are set; and an excitation matrix is constructed according to the multiple angle plane wave excitations.

[0057] Select the radar target scattering model with large size and complex structure, perform geometric modeling on the selected structure according to the design size of the selected structure and the continuity of the boundary surface between different structures, and create a radar target geometric model corresponding to the design requirements; Figure 2 A schematic diagram of the geometric model of an aircraft is given.

[0058] The created radar target geometric model is marked with material properties corresponding to actual design requirements, and the geometric bodies that require special design in the created radar target geometric model are marked with material properties corresponding to the actual design of the radar target scattering model whose electromagnetic characteristics are to be solved, and plane wave excitations at multiple angles are set; Figure 3 In the figure, number 1 is set as an ideal conductor material, and number 2 is a schematic diagram of plane wave excitation at multiple angles. Using existing methods, an excitation matrix is ​​constructed based on plane wave excitation at multiple angles, which is denoted as V N×M , N represents the number of unknowns, and M represents the number of excitation vectors.

[0059] S20. Perform triangular mesh subdivision on the radar target geometric model to generate a plurality of triangular meshes, and construct an impedance matrix based on the plurality of triangular meshes.

[0060] Perform triangular meshing on the created radar target geometric model to generate several triangular meshes; Figure 2 The aircraft model shown is divided into triangle meshes, and the resulting triangle mesh structures are as follows: Figure 4 shown.

[0061] The multi-layer fast multipole method is based on the method of moments, such as Figure 5 As shown in the figure, by grouping some adjacent nodes, and then using different switches to establish point-to-point communication links, and performing three processes: aggregation, transfer, and configuration; where aggregation refers to aggregating the current action in each group to the group center for expression; transfer refers to transferring the amount aggregated at the group center to the far-field group center; and configuration refers to configuring the value of the far-field group center to each current item. Based on the existing multi-layer fast multipole method of moment method, an impedance matrix is ​​constructed based on several triangular grids, denoted as Z N×N The specific process will not be described here.

[0062] Based on the design of S10 and S20, and also based on the existing multi-layer fast multipole method of moment method, the matrix equation can be constructed:

[0063] Z N×N I N×M =V N×M (1)

[0064] Among them, Z N×N represents the impedance matrix, I N×Mdenotes a current matrix to be solved, i.e., a second current matrix mentioned later, V N×M denotes an excitation matrix constructed under plane wave excitation of multiple angles. A conventional method directly solves formula (1), and there are problems such as that there are many unknowns in the solving process, and that the memory occupied and the time required for calculation of a matrix equation obtained through the moment method are long. Through research and analysis on the existing method, the present application proposes a process for quickly calculating formula (1), as shown in the following formula (2): Figure 6 Compared with the conventional moment method, the ACA mixed compression acceleration algorithm based on the multi-layer fast multipole method improves the ability to process large-size and complex target scattering characteristics, and specifically:

[0065] S30, the ACA algorithm is used to compress and decompose the excitation matrix to obtain a first full-rank matrix and a second full-rank matrix.

[0066] An optional solution provided by the present application is to compress and decompose the excitation matrix by using an adaptive cross approximation (ACA) algorithm, and the formula is represented as:

[0067] V N×M ≈A N×r ·B r×M (2)

[0068] Wherein, V N×M denotes an excitation matrix constructed under plane wave excitation of multiple angles, A N×r denotes a first full-rank matrix, B r×M denotes a second full-rank matrix, N denotes the number of unknowns, M denotes the number of excitation vectors, and r denotes the effective rank of the excitation matrix, and r < min(N, M).

[0069] S40, QR decomposition is performed on the first full-rank matrix to obtain a first triangular matrix and a first orthogonal matrix, and QR decomposition and Hermite matrix operation are performed on the second full-rank matrix to obtain a second triangular matrix and a second orthogonal matrix.

[0070] An optional solution provided by the present application is to perform QR decomposition on the first full-rank matrix A N×r obtained in S30 to obtain a first triangular matrix and a first orthogonal matrix, and the formula is represented as:

[0071] A N×r =AQ N×r ·AR r×r (3)

[0072] Wherein, A N×r denotes a first full-rank matrix, AQ N×r denotes a first orthogonal matrix, and AR r×rdenotes a first triangular matrix.

[0073] The application provides an alternative scheme, wherein the second full-rank matrix B r×M is subjected to QR decomposition and Hermite matrix operation to obtain a second triangular matrix and a second orthogonal matrix, and the formula is as follows:

[0074]

[0075] wherein, denotes a Hermite matrix of the second full-rank matrix, and BR r×r denotes the second triangular matrix, and BQ r×M denotes the second orthogonal matrix, and (·) H denotes obtaining a Hermite matrix.

[0076] S50, matrix multiplication is performed on the first triangular matrix and the second triangular matrix to obtain a first matrix.

[0077] The application provides an alternative scheme, wherein the first triangular matrix AR r×r obtained in S40 and the second triangular matrix BR r×r are subjected to matrix multiplication to obtain a first matrix, and the formula is as follows:

[0078] ABR r×r = AR r×r · BR r×r (5)

[0079] wherein, ABR r×r denotes the first matrix, AR r×r denotes the first triangular matrix, and BR r×r denotes the second triangular matrix.

[0080] S60, ACA algorithm is adopted to perform re-compression decomposition on the first matrix to obtain a third full-rank matrix and a fourth full-rank matrix.

[0081] The application provides an alternative scheme, wherein ACA algorithm is adopted to perform re-compression decomposition on the first matrix ABR r×r obtained in S50, and the formula is as follows:

[0082] ABR r×r = A1 r×k · B1 k×r (6)

[0083] wherein, ABR r×r denotes the first matrix, A1 r×k denotes the third full-rank matrix, B1 k×r denotes the fourth full-rank matrix, k denotes the effective rank number of the first matrix, and k < r.

[0084] S70, matrix multiplication operation is performed on the third full rank matrix and the first orthogonal matrix to obtain a second matrix, and matrix multiplication operation is performed on the fourth full rank matrix and the second orthogonal matrix to obtain a third matrix.

[0085] The application provides an alternative scheme, the third full rank matrix A1 r×k obtained in S60 is subjected to matrix multiplication operation with the first orthogonal matrix AQ N×r obtained in S40 to obtain a second matrix, and the formula is as follows:

[0086]

[0087] wherein, the second matrix is represented as A1 r×k the third full rank matrix is represented as AQ N×r the first orthogonal matrix is represented as.

[0088] The application provides an alternative scheme, the fourth full rank matrix B1 k×r obtained in S60 is subjected to matrix multiplication operation with the second orthogonal matrix BQ r×M obtained in S40 to obtain a third matrix, and the formula is as follows:

[0089]

[0090] wherein, the third matrix is represented as B1 k×r the fourth full rank matrix is represented as BQ r×M the second orthogonal matrix is represented as.

[0091] Through the above S30-S70, the excitation matrix V N×M under the excitation of the plane wave of multiple angles can be compressed and decomposed as follows:

[0092]

[0093] wherein, k < r < min (N, M).

[0094] S80, based on the multilayer fast multipole method, the second matrix is taken as a new excitation matrix, a matrix equation is constructed by combining a near-zone block diagonal preconditioning technique according to the impedance matrix and the new excitation matrix, and a generalized minimum residual method is used to solve the matrix equation to obtain a first current matrix.

[0095] The application provides an alternative scheme, based on the multilayer fast multipole method, the second matrix obtained in S70 is taken as a new excitation matrix, a matrix equation is constructed by combining a near-zone block diagonal preconditioning technique according to the impedance matrix obtained in S20 and the new excitation matrix, and the formula is as follows:

[0096]

[0097] wherein, Z N×N represents an impedance matrix, M N×N represents a preconditioning matrix obtained by a near-zone block-diagonal preconditioning technique, the preconditioning matrix M N×N , (·) -1 represents an inverse operation, I N×k represents a first current matrix, represents a second matrix.

[0098] The matrix equation constructed for formula (10) is solved by using a generalized minimal residual (GMRES) method to obtain a first current matrix I N×k , and the first current matrix is a partial current matrix.

[0099] S90, the first current matrix and the third matrix are subjected to matrix multiplication to obtain a second current matrix, and a radar target single-station RCS is calculated according to the second current matrix.

[0100] An alternative provided by the present application is that the first current matrix I N×k obtained by S80 and the third matrix obtained by S70 are subjected to matrix multiplication to obtain a second current matrix, and the formula is represented as:

[0101]

[0102] wherein, I N×M represents the second current matrix, I N×k represents the first current matrix, represents the third matrix. The second current matrix I N×M is the final complete current matrix, which is used to calculate the radar target single-station RCS.

[0103] According to the second current matrix I N×M , a radar target single-station RCS is calculated by using an existing method, and the detailed calculation process is not described herein.

[0104] In order to verify the effectiveness of the radar target single-station RCS fast angle scanning method provided by the embodiment of the present application, the following experiments are performed for verification.

[0105] 1. Simulation experiment conditions

[0106] The hardware platform of the simulation experiment of the present application is a workstation with a model of Intel Xeon Gold 5118R dual processor, a main frequency of 2.3 GHz, a single processor core number of 24, and a memory of 512 GB.

[0107] The software platforms for the simulation experiment of the present invention are: Win10 operating system and Visual Studio 2017.

[0108] The radar target scattering model selected for the simulation experiment of the present invention is Figure 2 The aircraft model shown has 3D dimensions of 15941.7 mm × 10625.45 mm × 2138.5 mm, a computation frequency of 500 MHz, and eight CPU cores. The geometric material is set to a perfect conductor. A plane wave is incident on the model at 721 angles in the polar coordinate system, with θ = 90° and φ = 0° to 360°, with φ intervals of 0.5. Horizontal polarization (in the +Y direction) is used for excitation. The monostatic RCS of the model is calculated for each of these 721 plane wave excitation angles.

[0109] 2. Simulation content and results analysis

[0110] Using the method proposed by the present invention, calculate Figure 2 Radar target scattering characteristics of the aircraft model shown under plane wave excitation at multiple angles. Figure 7 A schematic diagram is given comparing the results of the single-station RCS of the aircraft model under plane wave excitation at multiple angles using the existing multi-layer fast multipole method for discrete angle calculation and the method proposed in this invention for rapid calculation. Figure 7 The horizontal axis represents φ (Phi) in polar coordinates, in degrees (deg), and the vertical axis represents the radar target single-station RCS (Monostatic RCS Total), in decibel square meter (dBsm). Figure 7 It can be seen that the calculation results of the method proposed in the present invention are in good agreement with the angle-by-angle discrete calculation results of the existing multi-layer fast multipole method, which verifies the calculation accuracy of the technology of the present invention.

[0111] Table 1 shows the comparison results of the angle-by-angle discrete calculation of relevant parameters using the method proposed in the present invention and the existing multi-layer fast multipole method.

[0112] Table 1 Comparison results of two methods for calculating relevant parameters

[0113]

[0114] As can be seen from Table 1, the accelerated solution method of the present invention reduces the number of matrix equation solutions to 20.11% of that using the existing multi-layer fast multipole method, and the total calculation time is also reduced to 21.56% of the existing multi-layer fast multipole method. This greatly reduces the number of matrix equation solutions and the total calculation time, enabling rapid single-station RCS scanning of radar targets.

[0115] In summary, the radar target single station RCS fast scanning angle method provided by the embodiment of the application solves the problem of long time required by the existing multilayer fast multipole method for solving multiple angle plane wave excitation, is a new idea for realizing radar target single station RCS fast scanning, and the specific implementation process comprises: creating a radar target geometric model; wherein, materials are set for the radar target geometric model, and multiple angle plane wave excitation is set; an excitation matrix is constructed according to the multiple angle plane wave excitation; the radar target geometric model is triangularly meshed to generate a plurality of triangular meshes, and an impedance matrix is constructed according to the plurality of triangular meshes; ACA algorithm is used to compress and decompose the excitation matrix to obtain a first full rank matrix and a second full rank matrix; QR decomposition is performed on the first full rank matrix to obtain a first triangular matrix and a first orthogonal matrix, and QR decomposition and Hermite matrix operation are performed on the second full rank matrix to obtain a second triangular matrix and a second orthogonal matrix; the first triangular matrix and the second triangular matrix are multiplied to obtain a first matrix; ACA algorithm is used to compress and decompose the first matrix again to obtain a third full rank matrix and a fourth full rank matrix; the third full rank matrix and the first orthogonal matrix are multiplied to obtain a second matrix, and the fourth full rank matrix and the second orthogonal matrix are multiplied to obtain a third matrix; based on the multilayer fast multipole method, the second matrix is taken as a new excitation matrix, a matrix equation is constructed according to the impedance matrix and the new excitation matrix in combination with the near-zone block diagonal preconditioning technique, and the generalized minimum residual method is used to solve the matrix equation to obtain a first current matrix; the first current matrix and the third matrix are multiplied to obtain a second current matrix, and the radar target single station RCS is calculated according to the second current matrix. It can be seen that, firstly, the created radar target geometric model is triangularly meshed, then ACA algorithm is used to compress and process the excitation matrix, then a double-layer hybrid compression method combining ACA algorithm and QR decomposition is used, QR decomposition is performed on the compressed matrix, and ACA algorithm is used for re-compression, the number of excitations required for iterative solution in the solving process is further reduced, finally, the multilayer fast multipole method is used, and the matrix equation is constructed in combination with the near-zone block diagonal preconditioning technique, which greatly reduces the number of iterative solutions of the matrix equation, speeds up the solution of the matrix equation, and improves the solving efficiency of the radar target single station RCS under multiple angle plane wave excitation.

[0116] The embodiment of the present application realizes the radar target single-station RCS fast simulation solution, overcomes the problem of the prior art that only the radar target double-station RCS fast calculation is realized, so that the embodiment of the present application has the advantage of meeting the engineering actual demand; since the QR decomposition and the ACA algorithm are combined in the embodiment of the present application, only the current vector of the radar target under a few excitation vectors is solved, the problem of the prior art that only one angle excitation is accelerated is overcome, so that the method proposed in the present application has the advantage of stronger universality; since the ACA algorithm is used to compress the excitation matrix generated by the multi-layer fast multipole method in the embodiment of the present application, the problems of the prior art that the moment method is used to solve the long time and large memory occupation are overcome, so that the present application has the advantage of processing the scattering characteristics of large-size and complex targets; the embodiment of the present application can be used in colleges, research institutes, enterprises and other institutions to quickly calculate the radar target single-station RCS of the scattering target under the excitation of multiple angle plane waves, and has higher engineering application value.

[0117] Please refer to Figure 8 The embodiment of the present application provides an electronic device, which comprises a processor 801, a communication interface 802, a memory 803 and a communication bus 804, wherein the processor 801, the communication interface 802 and the memory 803 complete mutual communication through the communication bus 804.

[0118] The memory 803 is used for storing a computer program.

[0119] The processor 801 is used for executing the program stored in the memory 803, so as to realize the steps of the radar target single-station RCS fast scanning angle method.

[0120] The embodiment of the present application provides a computer readable storage medium, which stores a computer program, and the computer program is executed by a processor to realize the steps of the radar target single-station RCS fast scanning angle method.

[0121] For the device / electronic device / storage medium embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and the related part can be referred to the part of the method embodiment.

[0122] In the description of the present application, it should be understood that the terms "first", "second" are only used for the purpose of description, and cannot be understood as indicating or implying relative importance or implicitly indicating the number of the indicated technical features. Therefore, the features limited by "first", "second" can explicitly or implicitly include one or more features. In the description of the present application, the meaning of "multiple" is two or more, unless otherwise specifically limited.

[0123] Although the present application has been described in connection with certain embodiments, persons skilled in the art will understand and appreciate that many modifications can be made to the described embodiments and these modifications can be made without departing from the scope of the application in practical matter. Throughout the specification, the word "comprising" does not exclude the presence of other elements or additional steps. Throughout the specification, the word "a" or "an" shall not exclude a plurality of items. Several of the stated embodiments have been described above, but many modifications, alterations, and permutations can be made to these embodiments without departing from the application. Accordingly, the above description does not exhaust the scope of the application.

[0124] The above description is further detailed in connection with specific preferred embodiments of the application, and it is not to be construed that the specific implementation of the application is limited to these descriptions. For those skilled in the art, without departing from the concept of the application, a number of simple deductions or substitutions can be made, and all of these should be considered as falling within the scope of the protection of the application.

Claims

1. A radar target single-station RCS fast scanning method, characterized by: include: Creating a radar target geometric model; wherein, setting materials for the radar target geometric model and setting plane wave excitations at multiple angles; constructing an excitation matrix based on the plane wave excitations at multiple angles; Performing triangular mesh subdivision on the radar target geometric model to generate a plurality of triangular meshes, and constructing an impedance matrix based on the plurality of triangular meshes; Using an ACA algorithm to compress and decompose the excitation matrix to obtain a first full rank matrix and a second full rank matrix; Performing QR decomposition on the first full rank matrix to obtain a first triangular matrix and a first orthogonal matrix, and performing QR decomposition and Hermitian matrix operation on the second full rank matrix to obtain a second triangular matrix and a second orthogonal matrix; Performing a matrix multiplication operation on the first triangular matrix and the second triangular matrix to obtain a first matrix; Using the ACA algorithm to perform compression decomposition on the first matrix again to obtain a third full rank matrix and a fourth full rank matrix; Performing a matrix multiplication operation on the third full rank matrix and the first orthogonal matrix to obtain a second matrix, and performing a matrix multiplication operation on the fourth full rank matrix and the second orthogonal matrix to obtain a third matrix; Based on the multi-layer fast multipole method, the second matrix is ​​used as a new excitation matrix, a matrix equation is constructed according to the impedance matrix and the new excitation matrix in combination with a near-block diagonal preconditioning technique, and the matrix equation is solved using a generalized minimum residual method to obtain a first current matrix; A second current matrix is ​​obtained by performing a matrix multiplication operation on the first current matrix and the third matrix, and a single-station RCS of the radar target is obtained by calculating according to the second current matrix.

2. The radar target single-station RCS fast scanning angle method according to claim 1, characterized in that: The ACA algorithm is used to compress and decompose the excitation matrix, and the formula is expressed as follows: V N×M ≈A N×r ·B r×M ; Among them, V N×M represents the excitation matrix constructed under plane wave excitation at multiple angles, A N×r represents the first full rank matrix, B r×M represents the second full rank matrix, N represents the number of unknowns, M represents the number of excitation vectors, r represents the effective rank of the excitation matrix, and r <min(N,M)。 3. The radar target single-station RCS fast scanning angle method according to claim 2, characterized in that: Perform QR decomposition on the first full rank matrix to obtain a first triangular matrix and a first orthogonal matrix, which can be expressed as follows: A N×r =AQ N×r ·AR r×r ; Among them, A N×r represents the first full rank matrix, AQ N×r Represents the first orthogonal matrix, AR r×r represents the first triangular matrix.

4. The radar target single-station RCS fast scanning angle method according to claim 3, characterized in that: Perform QR decomposition and Hermitian matrix operations on the second full rank matrix to obtain a second triangular matrix and a second orthogonal matrix, which are expressed as follows: in, represents the Hermitian matrix of the second full rank matrix, BR r×r Represents the second triangular matrix, BQ r×M represents the second orthogonal matrix, (·) H Indicates the Hermitian matrix.

5. The radar target single-station RCS fast scanning angle method according to claim 4, characterized in that: The first triangular matrix and the second triangular matrix are multiplied to obtain a first matrix, which is expressed as follows: APR r×r =AR r×r ·BR r×r ; Among them, ABR r×r Represents the first matrix, AR r×r Represents the first triangular matrix, BR r×r represents the second triangular matrix.

6. The radar target single-station RCS fast scanning angle method according to claim 5, characterized in that: The ACA algorithm is used to compress and decompose the first matrix again, and the formula is expressed as: APR r×r =A1 r×k ·B1 k×r ; Among them, ABR r×r Represents the first matrix, A1 r×k represents the third full rank matrix, B1 k×r represents the fourth full rank matrix, k represents the effective rank of the first matrix, and k <r。 7. The radar target single-station RCS fast scanning angle method according to claim 6, characterized in that: The third full rank matrix and the first orthogonal matrix are multiplied to obtain a second matrix, which is expressed as follows: in, Represents the second matrix, A1 r×k represents the third full rank matrix, AQ N×r Represents the first orthogonal matrix.

8. The radar target single-station RCS fast scanning angle method according to claim 7, characterized in that: The fourth full rank matrix and the second orthogonal matrix are multiplied to obtain a third matrix, which is expressed as follows: in, Represents the third matrix, B1 k×r represents the fourth full rank matrix, BQ r×M Represents the second orthogonal matrix.

9. The radar target single-station RCS fast scanning angle method according to claim 8, characterized in that: Based on the multi-layer fast multipole method, the second matrix is ​​used as a new excitation matrix. According to the impedance matrix and the new excitation matrix, the matrix equation is constructed in combination with the near-block diagonal preconditioning technology. The formula is expressed as follows: Among them, Z N×N represents the impedance matrix, M N×N represents the preconditioned matrix obtained by the near-block diagonal preconditioning technique, (·) -1 Indicates the inverse operation, I N×k represents the first current matrix, Represents the second matrix.

10. The radar target single-station RCS fast scanning angle method according to claim 9, characterized in that: The first current matrix and the third matrix are multiplied to obtain a second current matrix, which is expressed as follows: Among them, I N×M Represents the second current matrix, I N×k represents the first current matrix, Represents the third matrix.

Images