CN114741646A - Method for rapidly solving wide-angle electromagnetic scattering characteristics of conductor target - Google Patents
Method for rapidly solving wide-angle electromagnetic scattering characteristics of conductor target Download PDFInfo
- Publication number
- CN114741646A CN114741646A CN202210213525.8A CN202210213525A CN114741646A CN 114741646 A CN114741646 A CN 114741646A CN 202210213525 A CN202210213525 A CN 202210213525A CN 114741646 A CN114741646 A CN 114741646A
- Authority
- CN
- China
- Prior art keywords
- characteristic
- target
- basis function
- solving
- mode
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 239000004020 conductor Substances 0.000 title claims abstract description 18
- 230000006870 function Effects 0.000 claims abstract description 39
- 239000011159 matrix material Substances 0.000 claims abstract description 22
- 230000009467 reduction Effects 0.000 claims description 3
- 108020001568 subdomains Proteins 0.000 claims description 3
- 238000010276 construction Methods 0.000 abstract description 10
- 238000004364 calculation method Methods 0.000 abstract description 9
- 230000005284 excitation Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000000903 blocking effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Operations Research (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a method for quickly solving wide-angle electromagnetic scattering characteristics of a conductor target, relates to the field of electromagnetic numerical calculation, and can effectively improve the speed of solving the wide-angle electromagnetic scattering characteristics of an electrically large-size conductor target by a characteristic basis function method. Firstly, a generalized eigenvalue equation is constructed according to the self-impedance matrix of the expanded subdomains, and the eigenmode of each subdomain is solved. Then, the effective mode is intercepted according to the mode significance and is used as a characteristic basis function for linearly expressing the target surface current information. The invention provides a new method for constructing the characteristic basis function by the characteristic basis function method, greatly improves the construction efficiency of the characteristic basis function, generates fewer characteristic basis functions, and improves the construction and solving speed of the reduced matrix equation.
Description
Technical Field
The invention relates to the technical field of electromagnetic numerical calculation, in particular to a method for quickly solving wide-angle electromagnetic scattering characteristics of an electrically large-size conductor target.
Background
The analysis of the wide-angle electromagnetic scattering characteristics of the target plays a crucial role in the technical development of radar detection, target recognition and the like. The Radar Cross Section (RCS) of a conductor target is calculated, and a surface integral equation is usually established and solved by a moment method (MoM). However, when the target electrical size is large, the computation complexity of the moment method and the memory requirement are both increased dramatically.
A fast algorithm based on a moment method, namely a Characteristic Basis Function Method (CBFM), can effectively reduce the number of unknowns to be solved, reduce the calculation time and the memory consumption, and is an important method for analyzing the electromagnetic scattering characteristics of an electrically large-size target. The characteristic basis function method is based on the region blocking principle, divides a target into a plurality of smaller subdomains, constructs a group of characteristic basis functions for each subdomain, and finally constructs a reduced matrix equation which can be solved by a direct method. However, for the analysis of the single-station RCS, the construction of the characteristic basis functions depends on a large number of incident plane excitations, the singular value decomposition time is long, the number of generated characteristic basis functions is large, and the solution efficiency is still low.
Disclosure of Invention
The invention aims to provide a method for quickly solving wide-angle electromagnetic scattering characteristics of a conductor target, so as to solve the problems in the background technology. The method is based on a characteristic basis function method, and takes a characteristic mode of each sub-domain as a characteristic basis function, so that the construction of the characteristic basis function is accelerated, the number of the characteristic basis functions is reduced, and the solving efficiency of the wide-angle electromagnetic scattering characteristic of the target is improved.
In order to achieve the purpose, the technical scheme of the invention comprises the following basic steps:
step 1: according to the basic principle of a characteristic basis function method, the whole conductor target surface is divided into a plurality of smaller sub-domains, each sub-domain is expanded, and then a matrix equation is constructed.
Step 2: and according to the characteristic mode theory of the conductor target, constructing a generalized characteristic value equation of each sub-domain, and solving the characteristic mode of each sub-domain.
And step 3: and setting a threshold condition according to the mode significance, and intercepting the effective mode meeting the threshold condition.
And 4, step 4: and removing the extension part of the intercepted characteristic module, and using the characteristic module as a characteristic basis function to linearly represent the target surface current information.
And 5: and constructing a reduction matrix equation and solving the target surface current.
And 6: and solving the value of the RCS of the target single station according to the result obtained in the step 5.
Compared with the prior art, the invention has the advantages that:
1. because the characteristic mode is irrelevant to excitation, the characteristic basis function and the reduced matrix constructed by the method can be repeatedly used for excitation at different incidence angles, and the efficiency of solving the wide-angle electromagnetic scattering characteristic of the target is guaranteed.
2. The invention adopts the characteristic mode theory of the conductor to directly solve the characteristic basis function from the impedance matrix. The construction procedure of the characteristic base function is simplified, and the construction time of the characteristic base function can be greatly reduced.
3. The number of the generated characteristic basis functions is obviously reduced, and the dimension of the reduced matrix can be effectively reduced, so that the construction and solving efficiency of the reduced matrix equation is improved.
Drawings
FIG. 1 is a basic flow diagram of the process of the present invention.
FIG. 2 is a diagram illustrating the calculation results of the method of the present invention.
Detailed Description
The basic flow diagram of the method of the invention is shown in fig. 1, and the implementation of the technical scheme is further described in detail in the following with the attached drawings:
step 1: the target is first divided into M adjacent sub-fields, and each sub-field is then divided into NiA unit (i ═ 1, 2.., M).Then, each sub-field is expanded, and the number of units of each sub-field after expansion is
Constructing a matrix equation according to the principle of a characteristic basis function method:
wherein ZiiAnd ZijRespectively a self-impedance matrix and a transimpedance matrix, ViFor incident plane wave excitation, JiRepresented as a corresponding coefficient matrix.
Step 2: according to the eigenmode theory of the conductor target, a self-impedance matrix Z of the expanded subdomains is usediiConstructing a generalized eigenvalue equation:
wherein, XiAnd RiAre each ZiiImaginary and real parts of, lambdaiIn order to be a value of the characteristic,
as a characteristic value λiThe corresponding feature vector is the feature mode of the expanded subfield i. Solving this equation to obtain λiAnd
and step 3: according to pattern significance
And setting a threshold value tau to meet the condition that the MS is more than or equal to the characteristic value of tau and the corresponding characteristic mode thereof are effective modes.
And 4, step 4: and removing an expanded part from the intercepted characteristic mode, and using the characteristic mode as a characteristic basis function to linearly represent target surface current information, wherein the surface current information on the sub-field i can be represented as:
wherein,
for the coefficients corresponding to the characteristic basis functions,
representing the p characteristic basis function of the ith sub-field.
And 5: bringing equation (3) into equation (1) in accordance with
Constructing a reduced matrix equation:
wherein
Is a matrix of coefficients of the characteristic basis functions,
in order to reduce the impedance matrix,
is that
The transposing of (1). And solving the reduction matrix equation to obtain the coefficient of the characteristic basis function, and solving the current of the target surface.
Step 6: and (5) solving the value of the RCS of the target single station according to the target surface current obtained in the step 5.
The method of the present invention is further described below in the context of a specific example, it being understood that the examples described are only a portion of the examples of the invention, and not all examples. All other embodiments, which can be obtained by a person skilled in the art without making any creative effort based on the embodiments in the present invention, belong to the protection scope of the present invention.
The invention is exemplified by a single station RCS that calculates a conductor cube of 1m side length. The target surface is subdivided at a pitch of 0.1 λ, where λ is the wavelength of the incident plane wave. The incident angle of the plane wave is
The frequency is 800 MHz. The target is divided into 8 subfields each extending for a width of 0.15 lambda, yielding 30322 unknowns in total. And respectively adopting the method and the traditional CBFM to solve the RCS of the target single station, wherein when the method is adopted, the threshold tau is 0.0005, and when the traditional CBFM is adopted, the SVD threshold is 0.001.
As can be seen from the attached figure 2, the calculation result of the method of the invention is well matched with MoM and CBFM, and has higher calculation precision. Table 1 shows the computation time, the number of characteristic basis functions, and the relative error of the method of the present invention and the method using the conventional CBFM method. It can be seen that, under the condition of obtaining almost the same calculation accuracy, the number of the characteristic basis functions generated by the method of the invention is reduced by 29%, the construction time of the characteristic basis functions is reduced by 96%, and the total calculation time is reduced by 75%. The example verifies the accuracy and the high efficiency of the invention.
TABLE 1
In summary, the invention solves the problem of low efficiency of the conventional CBFM characteristic basis function construction by solving the characteristic model of each sub-domain and using the characteristic model as the characteristic basis function of the sub-domain, and the number of generated characteristic basis functions is significantly reduced, thereby improving the construction and solving speed of the reduced matrix. The method has higher calculation efficiency and precision when solving the wide-angle electromagnetic scattering characteristic of the large-size conductor target.
The foregoing shows and describes the general principles, essential features, and advantages of the invention. It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, which are described in the specification and illustrated only to illustrate the principle of the present invention, but that various changes and modifications may be made therein without departing from the spirit and scope of the present invention, which fall within the scope of the invention as claimed.
Claims (4)
1. A method for rapidly solving wide-angle electromagnetic scattering characteristics of a conductor target is characterized by comprising the following steps:
step 1: according to the basic principle of a characteristic basis function method, the whole conductor target surface is divided into a plurality of smaller sub-domains, each sub-domain is expanded, and then a matrix equation is constructed.
Step 2: and according to the characteristic mode theory of the conductor target, constructing a generalized characteristic value equation of each sub-domain, and solving the characteristic mode of each sub-domain.
And step 3: and setting a threshold condition according to the mode significance, and intercepting the effective mode meeting the threshold condition.
And 4, step 4: and removing the extension part of the intercepted characteristic module, and using the characteristic module as a characteristic basis function to linearly represent the target surface current information.
And 5: and constructing a reduction matrix equation and solving the target surface current.
Step 6: and solving the value of the target single-station RCS according to the result obtained in the step 5.
2. The method as claimed in claim 1, wherein in step 2, the self-impedance matrix Z of the expanded sub-domains is used according to the eigenmode theory of the conductor targetiiConstructing a generalized eigenvalue equation:
in the formula, XiAnd RiAre each ZiiImaginary and real parts of, lambdaiIn order to be a characteristic value of the image,
as a characteristic value λiThe corresponding feature vector is the feature mode of the expanded subfield i.
3. The method as claimed in claim 1, wherein in step 3, the method is based on the significance of the pattern
And setting a threshold tau to meet the condition that the MS is more than or equal to the characteristic value of the tau and the corresponding characteristic mode is an effective mode.
4. The method for rapidly solving the wide-angle electromagnetic scattering property of the conductor target according to claim 1, wherein in step 4, the truncated characteristic mode is removed from the extension portion and is used as a characteristic basis function for linearly representing the surface current information of the target, and the surface current information on the sub-field i can be represented as:
wherein,
for the coefficients corresponding to the characteristic basis functions,
representing the p characteristic basis function of the ith sub-field.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210213525.8A CN114741646A (en) | 2022-03-04 | 2022-03-04 | Method for rapidly solving wide-angle electromagnetic scattering characteristics of conductor target |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210213525.8A CN114741646A (en) | 2022-03-04 | 2022-03-04 | Method for rapidly solving wide-angle electromagnetic scattering characteristics of conductor target |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114741646A true CN114741646A (en) | 2022-07-12 |
Family
ID=82275375
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210213525.8A Withdrawn CN114741646A (en) | 2022-03-04 | 2022-03-04 | Method for rapidly solving wide-angle electromagnetic scattering characteristics of conductor target |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114741646A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116502524A (en) * | 2023-04-18 | 2023-07-28 | 北京航空航天大学 | RCS reduction method for metal structure under broadband scanning |
-
2022
- 2022-03-04 CN CN202210213525.8A patent/CN114741646A/en not_active Withdrawn
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116502524A (en) * | 2023-04-18 | 2023-07-28 | 北京航空航天大学 | RCS reduction method for metal structure under broadband scanning |
| CN116502524B (en) * | 2023-04-18 | 2024-01-30 | 北京航空航天大学 | RCS reduction method for metal structure under broadband scanning |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110162847B (en) | Machine learning auxiliary antenna design method based on feature adding strategy | |
| CN111276822B (en) | Antenna pattern main lobe controllable antenna array design method | |
| CN105786765B (en) | A kind of method generating excitation extraneous features basic function quick self-adaptedly | |
| CN114741646A (en) | Method for rapidly solving wide-angle electromagnetic scattering characteristics of conductor target | |
| CN107145732B (en) | Method for rapidly solving target single-station electromagnetic scattering characteristics based on improved CBFM | |
| CN110276086A (en) | Stealthy super surface Radar cross-section redu-ction method based on space reflection | |
| CN107729627B (en) | Quick selection method for unit of opportunistic array antenna | |
| CN111931353A (en) | Scattered field solving method applied to simulated FSS structure | |
| CN114755652A (en) | Method for acquiring electrically large-size target broadband RCS (radar cross section) based on ACA (advanced communication architecture) and CAT (CAT) | |
| CN112784459A (en) | Electromagnetic simulation method based on compression type finite element tearing and butt joint method | |
| Guney et al. | Input resistance calculation for circular microstrip antennas using adaptive neuro-fuzzy inference system | |
| CN107305536B (en) | Mixed order time domain discontinuous Galois field gold method | |
| CN116882149A (en) | Antenna array synthesis method based on hybrid differential drosophila optimization algorithm | |
| Liu et al. | Application of Chaos Mind Evolutionary Algorithm in Antenna Arrays Synthesis. | |
| CN113567941B (en) | Rapid simulation method for near-field Doppler echo of target | |
| CN107346985B (en) | Interference alignment method combined with transmitting antenna selection technology | |
| CN114491992B (en) | Efficient electromagnetic scattering method based on equivalent dipole moment and physical optical method | |
| CN115688475A (en) | Efficient prediction method for low-pressure discharge threshold of microwave device of payload | |
| CN113051836B (en) | Rapid modeling method for antenna under machine learning-assisted array environment | |
| CN118278227B (en) | Efficient grid refinement self-adaptive moment method suitable for universal error indicator | |
| CN113240065A (en) | Passive radar station distribution method based on improved particle swarm optimization algorithm | |
| CN112836375B (en) | High-efficiency target electromagnetic scattering simulation method | |
| CN114722589B (en) | Method for rapidly solving three-dimensional target double-station RCS | |
| CN118114481A (en) | Method for obtaining wide-angle RCS of coating target based on greedy algorithm | |
| CN116879962B (en) | Efficient electromagnetic analysis method for dynamic group targets in sea surface environment |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WW01 | Invention patent application withdrawn after publication | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20220712 |