CN117706490B - Method for modeling coupling scattering center between metal targets based on single-station radar - Google Patents
Method for modeling coupling scattering center between metal targets based on single-station radar Download PDFInfo
- Publication number
- CN117706490B CN117706490B CN202410167220.7A CN202410167220A CN117706490B CN 117706490 B CN117706490 B CN 117706490B CN 202410167220 A CN202410167220 A CN 202410167220A CN 117706490 B CN117706490 B CN 117706490B
- Authority
- CN
- China
- Prior art keywords
- ray
- representing
- scattering
- scattering center
- radar
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 68
- 230000008878 coupling Effects 0.000 title claims abstract description 25
- 238000010168 coupling process Methods 0.000 title claims abstract description 25
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 25
- 239000002184 metal Substances 0.000 title claims abstract description 14
- 238000004364 calculation method Methods 0.000 claims abstract description 13
- 230000009466 transformation Effects 0.000 claims abstract description 5
- 239000013598 vector Substances 0.000 claims description 43
- 230000001419 dependent effect Effects 0.000 claims description 12
- 230000007246 mechanism Effects 0.000 claims description 9
- 238000000605 extraction Methods 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 7
- 238000003384 imaging method Methods 0.000 claims description 7
- 230000010354 integration Effects 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 3
- 238000012512 characterization method Methods 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 230000005684 electric field Effects 0.000 claims description 3
- 230000009191 jumping Effects 0.000 claims description 3
- 230000035699 permeability Effects 0.000 claims description 3
- 238000010408 sweeping Methods 0.000 claims description 3
- 238000002592 echocardiography Methods 0.000 claims description 2
- 230000008569 process Effects 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 6
- 240000008791 Antiaris toxicaria Species 0.000 description 1
- 241000321453 Paranthias colonus Species 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
Classifications
-
- 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
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a method for modeling a coupling scattering center between metal targets based on a single-station radar, which comprises the following steps: firstly, geometrically modeling the whole object by using computer simulation software, decomposing and splitting components of the object, numbering different components and surface elements in sequence, and carrying out shielding judgment and ray tracing on all the surface elements by using a bouncing ray method; classifying the rays according to the bouncing times and paths of the rays, solving the equivalent points of each ray according to the Fourier transformation relation between the target echo and the shape function by using the principle of phase equivalence, and obtaining the distribution of the equivalent points; then calculating the projection sizes of equivalent points in each set in the distance direction and the azimuth direction, and establishing the type of a scattering center; and finally, respectively establishing relevant parameters of the position, the amplitude, the length and the frequency factor according to the type of the scattering center to obtain an attribute scattering center model. The method has the advantages of high calculation process efficiency, concise model description and high accuracy of the established coupling scattering center model.
Description
Technical Field
The invention belongs to the technical field of electromagnetic calculation, and particularly relates to a method for modeling a coupling scattering center between metal targets based on a single-station radar.
Background
Since the fifth sixty of the last century, radar target scattering center models have attracted considerable attention as radar identification systems have evolved. The methods of constructing scattering center models can be generally divided into two categories: the method is a reverse method based on parameter estimation, is not limited by the target and the radar resolution, but the parameter estimation method is low in precision and consumes huge computing resources; the other type is a forward method based on a target geometric model, reflects fine physical information such as the shape and the size of the target, can approach the real target more accurately, and is more efficient.
The attribute scattering center model contains the dependence of the scattering field of the target on frequency and azimuth angle, has definite physical significance, and can explain the shape, the gesture and the position information of the target. In recent years, various methods of modeling the attribute scattering center have been proposed. Document 1(H. Kang, L. Xu, H. Zhao, J. Li, Q. Guo and J. Zhang, "A Robust Scattering Center Extraction Method and Its Application in In ISAR Imaging," in IEEE Sensors Journal, vol. 23, no. 21, pp. 26302-26310, 1 Nov.1, 2023.) proposes a scattering center extraction method based on distance-azimuth profile reconstruction, but the algorithm is a reverse method, which requires higher calculation time cost and limits the application in the extraction of the target coupling scattering center. Document 2(K. Huang, S. -Y. He, Y. -H. Zhang, H. -C. Yin, Z. -D. Bian and G. -Q. Zhu, "Composite Scattering Analysis of the Ship on a Rough Surface Based on the Forward Parametric Scattering Center Modeling Method," in IEEE Antennas and Wireless Propagation Letters, vol. 18, no. 12, pp. 2493-2497, Dec. 2019.) proposes a forward modeling method for extracting a scattering center from a rough surface combination model of a target, which acquires a model parameter of the scattering center of the target from a CAD model, but adopts the principle of equivalent optical path difference when the coupling scattering center is processed, and requires equivalent propagation paths of rays for a plurality of times when the coupling scattering center is processed, so that the calculation efficiency is reduced. Patent CN110083915A discloses a forward automatic modeling method for a radar target scattering center in a high-frequency region, and the forward automatic calculation of the scattering center attribute parameter has clear corresponding relation with a target structure, thereby being more beneficial to target identification. However, this method is not suitable for this type of target because of its low accuracy in processing targets with strongly coupled scattering centers.
Disclosure of Invention
The invention aims to provide a method for modeling the coupling scattering center between metal targets based on single-station radar, which has the advantages of high model description accuracy, high calculation process efficiency, accurate position and simple model description.
The technical solution for realizing the purpose of the invention is as follows: a method for modeling a coupling scattering center between metal targets based on a single-station radar comprises the following steps:
Step 1, geometrically modeling the whole object by using computer simulation software, then decomposing and splitting the component of the object by taking a normal vector of the surface element and discontinuity of the surface element as the basis, numbering different components and surface elements in sequence, and carrying out shielding judgment and ray tracing on all the surface elements by using a bouncing ray method;
step 2, classifying rays according to the bouncing times and paths of the rays, solving equivalent points of each ray according to Fourier transformation relation between the target echo and the shape function by using a phase equivalent principle, and obtaining the distribution of the equivalent points;
step 3, calculating the projection sizes of equivalent points in each set in the distance direction and the azimuth direction, and establishing the type of the scattering center;
And 4, respectively establishing relevant parameters of the position, the amplitude, the length and the frequency factor according to the type of the scattering center to obtain an attribute scattering center model, and realizing target attribute feature extraction.
In step 1, the whole object is geometrically modeled by using computer simulation software, then the object is decomposed and split by taking the normal vector of the surface element and the discontinuity of the surface element as the basis, different parts and surface elements are numbered in sequence, and all the surface elements are subjected to shielding judgment and ray tracing by using a bouncing ray method, which is specifically as follows:
Step 1.1, setting a plurality of surfaces with different scattering mechanisms of all targets, and geometrically modeling the whole target by using computer simulation software;
Step 1.2, decomposing a part of the target, and numbering each surface in sequence;
Step 1.3, dividing all the faces into triangular surface elements, and numbering all the triangular surface elements and points;
step 1.4, shielding judgment is carried out on all the surface elements by using a bouncing ray method;
And 1.5, carrying out ray tracing on each triangular surface element in sequence.
Further, in step 1.2, the target is subjected to component decomposition, following the following principles:
(1) Decomposing the target to obtain a set of all entity parts;
(2) Taking abrupt change of normal vector of the target surface as a basis to decompose the component;
(3) The scattering mechanism of each part of the split target is clear, and one part does not have multiple scattering mechanisms.
Further, in step 1.3, all the faces are split into triangle surface elements, and all the triangle surface elements and points are numbered as follows:
dividing all the faces into triangular surface elements, wherein the side length of the triangular surface elements is not less than ，/>To solve for the wavelength used in the scattering center amplitude using physical optics, all triangle bins and points are numbered and the part number of each triangle bin, the three vertex numbers of each triangle bin and the three-dimensional coordinates of each vertex are recorded.
Further, in step 1.5, ray tracing is performed on each triangle surface element in sequence, which is specifically as follows:
step 1.5.1 determining an incident ray according to the radar incident direction Wherein/>And/>Respectively a pitch angle and an azimuth angle of the radar incidence direction under a spherical coordinate system;
step 1.5.2 incident rays for each triangle bin The first reflection point of the triangle is the geometric center of the triangle surface element contacted by the ray for the first time, and the reflected ray in the first reflection direction is calculated according to the geometric optics principle;
step 1.5.3, judging whether the reflected ray and other triangle surface elements have intersection points: if the intersection point exists, recording the part number of the triangle surface element and the number of the triangle surface element, calculating the direction of the next reflection, and entering a step 1.5.4; if the intersection point does not exist, directly entering the step 1.5.5;
step 1.5.4, judging whether the ejection times of the same incident ray are more than 6 times: if not more than 6 times, returning to the step 1.5.3; if more than 6 times, discard the ray traced this time, and jump to analysis of the next ray;
And 1.5.5, judging whether the included angle between the last reflection direction and the radar incidence direction is larger than 90 degrees, if not smaller than 90 degrees, discarding the ray traced at the time, jumping to analysis of the next ray, and if smaller than 90 degrees, considering the ray to be effective.
Further, in step 2, the rays are classified according to the bouncing times and paths of the rays, and the equivalent points of each ray are solved by using the principle of phase equivalence according to the fourier transform relationship between the target echo and the shape function, so as to obtain the distribution of the equivalent points, which is specifically as follows:
Step 2.1, respectively solving equivalent points which can generate the same electric field contribution with the ray bounced for multiple times according to the bounced times of the ray, wherein if the bounced times are 1, the geometric center of the triangle surface element is the equivalent point; if the bounce times is greater than 1, solving the equivalent point of each ray according to the principle of phase equivalence, wherein the method comprises the following steps:
According to the high frequency approximation method, the far field fringe field generated by the induced current is expressed as:
(1)
wherein, Representing ray-generated far-field scatter fields; /(I)Representing imaginary units; /(I)Representing wave numbers; /(I)Representing a field point location vector in far field computation; /(I)Representing a source point position vector,/>Representing the angular frequency of the incident wave; /(I)Representing the amplitude of the incident wave; /(I)And/>The normal unit vector and the tangential unit vector of the surface element contacted by the ray for the last time are respectively represented; represents permeability in vacuum; /(I) Representing phase terms due to multiple reflections of rays,/>Represents the first/>, of the bar raySecondary reflection,/>Representing the total number of reflections,/>Representing the incident plane wave vector,/>Represents the/>Reflected wave vector after the secondary contact surface element,/>Reflected wave vector representing last contact surface source,/>Is ray and/>A position vector of the intersection of the secondary contact bins; Representing a three-dimensional green's function;
step 2.2, in the inverse aperture radar imaging, the radar echo of multiple catapulting is regarded as a Fourier transform form of a shape function:
(2)
Wherein the method comprises the steps of Is a shape function of the ray; wave vector/>, is carried out on echo obtained by sweeping sweep angleIs a shape function obtained by two-dimensional Fourier transform of (a) and has/>, for a single-station signalThe radar echo is rewritten as:
(3)
Wherein the method comprises the steps of =/>；/>A position vector representing an equivalent point;
Step 2.3, respectively making radar incident direction And/>Rotating an angle to set 0.001 DEG, the last reflected wave vectors are/>, respectivelyAnd/>；/>、/>And/>A group of bases constituting a three-dimensional space, so/>With and without only a unique solution, the contact point and the direction ejected in the middle of the target are not changed, the equivalent point position is unchanged, and the following equation set is solved to obtain/>：
(4)
Further, in step3, the projection sizes of the equivalent points in each set in the distance direction and the azimuth direction are calculated, and the scattering center type is established, specifically as follows:
Step 3.1, according to attribute scattering center theory, total scattering field of the target Expressed as:
(5)
wherein, Represents the/>The magnitude of the individual scattering centers; /(I)Representing radar frequency,/>Representing the radar center frequency; Representing a frequency dependent factor,/> Representing the dependence of azimuth angle,/>；/>Represents the/>The length of the individual scattering centers; /(I)Representing the propagation velocity of the representative wave in vacuum; /(I)Representing the center azimuth of the radar; /(I)Represents the/>The location of the individual scattering centers;
step 3.2, traversing all rays in each set, respectively calculating projections of equivalent points in the set in azimuth directions and obtaining maximum and minimum values of the projections And/>The projection length of the set in the azimuth direction is=/>；
Step 3.3, setting a decision threshold in the azimuth directionWhen/></>When the scattering center model corresponds to a local scattering center,/>Representing the dependence of scattering amplitude and azimuth angle of the local scattering center;
When (when) When the scattering center model corresponds to the distributed scattering center, the dependence of the scattering center model on the azimuth angle is represented by the scattering center length/>And tilt angle/>Characterization.
Further, in step 4, according to the type of the scattering center, relevant parameters of position, amplitude, length and frequency factor are respectively established to obtain an attribute scattering center model, and extraction of target attribute features is achieved, specifically as follows:
Step 4.1, length of local scattering center Frequency dependent factor/>All 0, position/>Sum amplitudeIs calculated from the following formula:
(6)
wherein, Representing the position of equivalent points in a collection,/>Representing the magnitude of the equivalent points in the set,；
Step 4.2, calculating projections of equivalent points in the set in the distance direction and obtaining maximum and minimum values of the projections respectively, wherein the length of the distributed scattering center in the azimuth direction is not 0 and the projections depend on the azimuth angle and rays in the setAnd/>The projection of the set in the distance direction is/>=/>；
Step 4.3, setting a threshold value for distance upWhen/></>Time,/>=/>Frequency dependent factor/>The values of (2) are as follows: when the equivalent point distribution is a curve,/>; When the equivalent points are distributed as straight lines or hyperbolas,; When the equivalent point distribution is cylindrical,/>; When the equivalent point distribution is planar,/>; Position/>Sum amplitude/>The calculation formula of (2) is as follows:
(7)
Wherein the method comprises the steps of And/>The positions of the two points of minimum and maximum projection in azimuth;
When (when) >/>When the scattered field is:
(8)
Wherein the method comprises the steps of Representing the integration region,/>Representing the magnitude function at the integrated,/>Representing a dot/>The included angle between the scattering direction and the radar incidence direction, the scattered field is written as a fractional integral:
(9)
Wherein the method comprises the steps of Representation/>At/>N-th derivative of the position,/>Two endpoints of the integral area in azimuth;
Due to Only get/< < 1)In one aspect, the fringe field is reduced to:
(10)
The scattering center of the collection is equivalent to two points And/>Length/>、/>Frequency dependent factor/>、All 0, position/>、/>Sum amplitude is/>、/>The calculation formula of (2) is as follows:
(11)
wherein, And/>Coordinates of two endpoints of the integration region respectively; /(I)And/>Is the fringe field amplitude of the corresponding endpoint.
Compared with the prior art, the invention has the remarkable advantages that: (1) Considering the coupling scattering centers among the targets, the description of the model is more accurate; (2) Solving the equivalent points by a phase equivalent method, so that the calculation process is more efficient and the position is more accurate; (3) And the scattering centers are classified according to the point cloud of the equivalent points, so that the description of the model is simpler.
Drawings
FIG. 1 is a flow chart of a method for modeling coupled scattering centers between metallic targets based on a single-station radar of the present invention.
FIG. 2 is a schematic diagram of the geometry of SLICY models in an embodiment of the present invention.
FIG. 3 is a schematic diagram of an exploded geometry of a SLICY model in an embodiment of the present invention.
FIG. 4 is a schematic diagram of a ray tracing process in an embodiment of the invention.
Fig. 5 is a schematic view of a local scattering center in an embodiment of the invention.
FIG. 6 is a schematic diagram of a distributed scattering center when the distance projection is less than a threshold in an embodiment of the present invention.
FIG. 7 is a schematic diagram of a distributed scattering center when the distance projection is greater than a threshold in an embodiment of the present invention.
FIG. 8 is a schematic view of a scattering center of SLICY model in an embodiment of the present invention.
FIG. 9 is a diagram of SLICY model-coupled scattering center ISAR in an embodiment of the invention.
Detailed Description
The invention is described in further detail below with reference to the drawings and the specific examples.
With reference to a flowchart shown in fig. 1, the modeling method of the coupling scattering center between metal targets based on the single-station radar comprises the following steps:
Step 1, geometrically modeling the whole object by using computer simulation software, then decomposing and splitting the component of the object by taking a normal vector of the surface element and discontinuity of the surface element as the basis, numbering different components and surface elements in sequence, and carrying out shielding judgment and ray tracing on all the surface elements by using a bouncing ray method;
step 2, classifying rays according to the bouncing times and paths of the rays, solving equivalent points of each ray according to Fourier transformation relation between the target echo and the shape function by using a phase equivalent principle, and obtaining the distribution of the equivalent points;
step 3, calculating the projection sizes of equivalent points in each set in the distance direction and the azimuth direction, and establishing the type of the scattering center;
And 4, respectively establishing relevant parameters of the position, the amplitude, the length and the frequency factor according to the type of the scattering center to obtain an attribute scattering center model, and realizing target attribute feature extraction.
As a specific example, in step 1, the computer simulation software is used to perform geometric modeling on the whole object, then the component decomposition subdivision is performed on the object based on the normal vector of the surface element and the discontinuity of the surface element, different components and surface elements are numbered in sequence, the shielding judgment and ray tracing are performed on all the surface elements by using the bouncing ray method, the size of SLICY model is shown as figure 2, the lengths and the angles of the opposite angles of each side of the model are marked in the figure, the different components of SLICY model are numbered as 1-15 in the figure in combination with the figure 3, only the bright area component under the incidence of plane waves is marked, and the invention is described by using Sandia Laboratories Implementation of Cylinders (SLICY) model, which is specifically as follows:
Step 1.1, setting a plurality of surfaces with different scattering mechanisms of all targets, and geometrically modeling the whole target by using computer simulation software;
Step 1.2, decomposing a part of the target, and numbering each surface in sequence;
Step 1.3, dividing all the faces into triangular surface elements, and numbering all the triangular surface elements and points;
step 1.4, shielding judgment is carried out on all the surface elements by using a bouncing ray method;
And 1.5, carrying out ray tracing on each triangular surface element in sequence.
As a specific example, in step 1.2, when the target is decomposed in parts, the following principle needs to be followed:
(1) Decomposing the target to obtain a set of all entity parts;
(2) Taking abrupt change of normal vector of the target surface as a basis to decompose the component;
(3) Ensuring that the scattering mechanism of the individual components of the split target is sufficiently simple that there are no multiple scattering mechanisms present in one component.
As a specific example, in step 1.3, all the faces are subdivided into triangular bins, and all the triangular bins and points are numbered as follows:
dividing all the faces into triangular surface elements, wherein the side length of the triangular surface elements is not less than ，/>To solve for the wavelength used in the scattering center amplitude using physical optics, all triangle bins and points are numbered and the part number of each triangle bin, the three vertex numbers of each triangle bin and the three-dimensional coordinates of each vertex are recorded.
As a specific example, in step 1.5, each triangle is followed in turn by ray tracing, in conjunction with FIG. 4, whereRepresenting the incident plane wave vector,/>Reflected wave vector representing four reflections,/>The intersection of four bins is represented as follows:
step 1.5.1 determining an incident ray according to the radar incident direction Wherein/>And/>Respectively a pitch angle and an azimuth angle of the radar incidence direction under a spherical coordinate system;
step 1.5.2 incident rays for each triangle bin The first reflection point of the triangle is the geometric center of the triangle surface element contacted by the ray for the first time, and the reflected ray in the first reflection direction is calculated according to the geometric optics principle;
step 1.5.3, judging whether the reflected ray and other triangle surface elements have intersection points: if the intersection point exists, recording the part number of the triangle surface element and the number of the triangle surface element, calculating the direction of the next reflection, and entering a step 1.5.4; if the intersection point does not exist, directly entering the step 1.5.5;
step 1.5.4, judging whether the ejection times of the same incident ray are more than 6 times: if not more than 6 times, returning to the step 1.5.3; if more than 6 times, discard the ray traced this time, and jump to analysis of the next ray;
And 1.5.5, judging whether the included angle between the last reflection direction and the radar incidence direction is larger than 90 degrees, if not smaller than 90 degrees, discarding the ray traced at the time, jumping to analysis of the next ray, and if smaller than 90 degrees, considering the ray to be effective.
As a specific example, in step 2, the rays are classified according to the number of bounces and paths of the rays, and according to the fourier transform relationship between the target echo and the shape function, the equivalent points of each ray are solved by using the principle of phase equivalence, so as to obtain the distribution of the equivalent points, which is specifically as follows:
Step 2.1, respectively solving equivalent points which can generate the same electric field contribution with the ray bounced for multiple times according to the bounced times of the ray, wherein if the bounced times are 1, the geometric center of the triangle surface element is the equivalent point; if the bounce times is greater than 1, solving the equivalent point of each ray according to the principle of phase equivalence, wherein the method comprises the following steps:
According to the high frequency approximation method, the far field fringe field generated by the induced current can be expressed as:
(1)
wherein, Representing ray-generated far-field scatter fields; /(I)Representing imaginary units; /(I)Representing wave numbers; /(I)Representing a field point location vector in far field computation; /(I)Representing a source point position vector,/>Representing the angular frequency of the incident wave; /(I)Representing the amplitude of the incident wave; /(I)And/>The normal unit vector and the tangential unit vector of the surface element contacted by the ray for the last time are respectively represented; represents permeability in vacuum; /(I) Representing phase terms due to multiple reflections of rays,/>Represents the first/>, of the bar raySecondary reflection,/>Representing the total number of reflections,/>Representing the incident plane wave vector,/>Represents the/>Reflected wave vector after the secondary contact surface element,/>Reflected wave vector representing last contact surface source,/>Is ray and/>A position vector of the intersection of the secondary contact bins; Representing a three-dimensional green's function;
step 2.2, in inverse aperture radar imaging, the radar echoes of multiple shots can be regarded as fourier transform versions of their shape functions:
(2)
Wherein the method comprises the steps of Is a shape function of the ray; wave vector/>, is carried out on echo obtained by sweeping sweep angleIs a shape function obtained by two-dimensional Fourier transform of (a) and has/>, for a single-station signalThe radar echo is rewritten as:
(3)
Wherein the method comprises the steps of =/>；/>A position vector representing an equivalent point;
Step 2.3, respectively making radar incident direction And/>Rotating a very small angle, setting to 0.001 DEG, the last reflected wave vectors are/>, respectivelyAnd/>；/>、/>And/>A group of bases constituting a three-dimensional space, so/>There is only one solution, because the deflection angle is extremely small, it can be considered that the contact point and the direction ejected in the middle of the target are not changed, the equivalent point position is unchanged, and the following equation system can be solved to obtain/>：
(4)
As a specific example, in step 3, the projection sizes of the equivalent points in each set in the distance direction and the azimuth direction are calculated, and the scattering center type is established, specifically as follows:
Step 3.1, according to attribute scattering center theory, total scattering field of the target Can be expressed as:
(5)
wherein, Represents the/>The magnitude of the individual scattering centers; /(I)Representing radar frequency,/>Representing the radar center frequency; Representing a frequency dependent factor,/> Representing the dependence of azimuth angle,/>；/>Represents the/>The length of the individual scattering centers; /(I)Representing the propagation velocity of the representative wave in vacuum; /(I)Representing the center azimuth of the radar; /(I)Represents the/>The location of the individual scattering centers;
Step 3.2, traversing all rays in each set, calculating projections of equivalent points in the set in azimuth directions and obtaining maximum and minimum values of the projections And/>The projected length of the set in azimuth is/>=/>；
Step 3.3, setting a decision threshold in the azimuth directionWhen/></>When the scattering center model corresponds to a local scattering center,/>Representing the dependence of scattering amplitude and azimuth angle of the local scattering center;
When (when) When the scattering center model corresponds to the distributed scattering center, the dependence of the scattering center model on the azimuth angle is represented by the scattering center length/>And tilt angle/>Characterization.
As a specific example, in step 4, according to the type of the scattering center, the relevant parameters of the position, the amplitude, the length and the frequency factor are respectively established, so as to obtain an attribute scattering center model, which is specifically as follows:
step 4.1, FIG. 5 shows the equivalent dot matrix of the corner reflector portion in SLICY model, the length of the local scattering center Frequency dependent factor/>All 0, position/>Sum amplitude/>Is calculated from the following formula:
(6)
wherein, Representing the position of equivalent points in a collection,/>Representing the magnitude of the equivalent points in the set,；
Step 4.2, calculating projections of equivalent points in the set in the distance direction and obtaining maximum and minimum values of the projections respectively, wherein the length of the distributed scattering center in the azimuth direction is not 0 and the projections depend on the azimuth angle and rays in the setAnd/>The projection of the set in the distance direction is/>=/>；
Step 4.3, setting a threshold value for distance upAs shown in FIG. 6, when/></>Time,/>=/>Frequency dependent factor/>The values of (2) are as follows: when the equivalent point distribution is a curve,/>; When the equivalent point distribution is straight or hyperboloid,/>; When the equivalent point distribution is cylindrical,/>; When the equivalent point distribution is planar,/>; Position ofSum amplitude/>The calculation formula of (2) is as follows:
(7)
Wherein the method comprises the steps of And/>The positions of the two points of minimum and maximum projection in azimuth;
as shown in fig. 7, when >/>When the scattered field is:
(8)/>
Wherein the method comprises the steps of Representing the integration region,/>Representing the magnitude function at the integrated,/>Representing a dot/>The included angle between the scattering direction and the radar incidence direction, the scattered field is written as a fractional integral:
(9)
Wherein the method comprises the steps of Representation/>At/>N-th derivative of the position,/>Two endpoints of the integral area in azimuth;
Due to Only/< <1 >, can be takenOne of the time points, the scattered field can be simplified as:
(10)
The scattering center of the collection is equivalent to two points And/>Length/>、/>Frequency dependent factor/>、All 0, position/>、/>Sum amplitude is/>、/>The calculation formula of (2) is as follows:
(11)
wherein, And/>Coordinates of two endpoints of the integration region respectively; /(I)And/>Is the fringe field amplitude of the corresponding endpoint.
Examples
Aiming at the problems of difficult modeling and low imaging precision of the coupling scattering centers among metal targets, the traditional scattering center model modeling method has long time and no wide applicability.
In this embodiment, SLICY models are taken as an example, and radar parameters are as follows: radar center frequencyBandwidth ofThe number of frequency points is 81, and the pitch angle/>, of the radarRadar sweep angle width/>The number of scan angles is 51, and the polarization mode VV is polarized. Threshold in distance upward/>And determination threshold in azimuth/>Taking 0.1 and 0.3 respectively, the scattering center modeling results and ISAR imaging results are shown in fig. 8 and 9 respectively, there are 5 coupling scattering centers and one specular scattering center: SC1 is the coupling scattering center generated between part 3 and part 4,/>，，/>,/>；/>SC2 is the coupling scattering center generated between part 7 and part 8,/>，/>,/>; SC3 is the coupling scattering center generated between part 9 and part 12,/>，/>，,/>，/>，/>，，/>，/>，/>，,/>; SC4 is the coupling scattering center generated between the component 9 and the component 14,，/>，/>,/>，，/>，/>，/>，，/>，/>,/>; SC5 is the coupling scattering center generated between the components 9, 10, 11,/>，/>，，/>. The scattering center of the single reflection can acquire parameters according to the traditional attribute scattering center method, and SC6 is the scattering center generated by the specular reflection of the component 2,/>，/>，，/>. The ISAR imaging result is shown in fig. 9, and from the result, it can be seen that the scattering center model constructed by the method is identical to the actual model structure, and compared with the traditional method, the modeling accuracy of the target scattering center is higher.
The foregoing is merely a preferred embodiment of the present invention and it should be noted that modifications and adaptations to those skilled in the art may be made without departing from the principles of the present invention, which are intended to be comprehended within the scope of the present invention.
Claims (7)
1. A method for modeling a coupling scattering center between metal targets based on a single-station radar is characterized by comprising the following steps:
Step 1, geometrically modeling the whole object by using computer simulation software, then decomposing and splitting the component of the object by taking a normal vector of the surface element and discontinuity of the surface element as the basis, numbering different components and surface elements in sequence, and carrying out shielding judgment and ray tracing on all the surface elements by using a bouncing ray method;
step 2, classifying rays according to the bouncing times and paths of the rays, solving equivalent points of each ray according to Fourier transformation relation between the target echo and the shape function by using a phase equivalent principle, and obtaining the distribution of the equivalent points;
step 3, calculating the projection sizes of equivalent points in each set in the distance direction and the azimuth direction, and establishing the type of the scattering center;
step 4, respectively establishing relevant parameters of position, amplitude, length and frequency factor according to the type of the scattering center to obtain an attribute scattering center model, and realizing target attribute feature extraction;
In the step 2, the rays are classified according to the bouncing times and paths of the rays, and the equivalent points of each ray are solved by using the principle of phase equivalence according to the Fourier transformation relation between the target echo and the shape function, so as to obtain the distribution of the equivalent points, and the method specifically comprises the following steps:
Step 2.1, respectively solving equivalent points which can generate the same electric field contribution with the ray bounced for multiple times according to the bounced times of the ray, wherein if the bounced times are 1, the geometric center of the triangle surface element is the equivalent point; if the bounce times is greater than 1, solving the equivalent point of each ray according to the principle of phase equivalence, wherein the method comprises the following steps:
According to the high frequency approximation method, the far field fringe field generated by the induced current is expressed as:
(1)
wherein, Representing ray-generated far-field scatter fields; /(I)Representing imaginary units; /(I)Representing wave numbers; /(I)Representing a field point location vector in far field computation; /(I)Representing the angular frequency of the incident wave; /(I)Representing the amplitude of the incident wave; /(I)And/>The normal unit vector and the tangential unit vector of the surface element contacted by the ray for the last time are respectively represented; /(I)Represents permeability in vacuum; representing phase terms due to multiple reflections of rays,/> Representing ray numberSecondary reflection,/>Representing the total number of reflections,/>Representing the incident plane wave vector,/>Represents the/>Reflected wave vector after the secondary contact surface element,/>Reflected wave vector representing last contact surface source,/>Is ray and/>A position vector of the intersection of the secondary contact bins; /(I)Representing a three-dimensional green's function;
step 2.2, in the inverse aperture radar imaging, the radar echo of multiple catapulting is regarded as a Fourier transform form of a shape function:
(2)
Wherein the method comprises the steps of Is a shape function of the ray; concerning wave vectors of echoes obtained by sweeping anglesIs a shape function obtained by two-dimensional Fourier transform of (a) and has/>, for a single-station signalThe radar echo is rewritten as:
(3)
Wherein the method comprises the steps of =/>；/>A position vector representing an equivalent point;
Step 2.3, respectively making radar incident direction And/>Rotating an angle to set 0.001 DEG, the last reflected wave vectors are/>, respectivelyAnd/>；/>、/>And/>A group of bases constituting a three-dimensional space, so/>With only one solution, the contact point and the direction ejected in the middle of the target are not changed, the equivalent point position is unchanged, and the following equation set is solved to obtain：
(4)
2. The modeling method of coupling scattering centers between metal targets based on single-station radar according to claim 1, wherein in step 1, the computer simulation software is used to perform geometric modeling on the whole target, then the target is subjected to component decomposition and subdivision based on the normal vector of the surface element and the discontinuity of the surface element, different components and surface elements are numbered in sequence, and all the surface elements are subjected to shielding judgment and ray tracing by using a bouncing ray method, which is specifically as follows:
Step 1.1, setting a plurality of surfaces with different scattering mechanisms of all targets, and geometrically modeling the whole target by using computer simulation software;
Step 1.2, decomposing a part of the target, and numbering each surface in sequence;
Step 1.3, dividing all the faces into triangular surface elements, and numbering all the triangular surface elements and points;
step 1.4, shielding judgment is carried out on all the surface elements by using a bouncing ray method;
And 1.5, carrying out ray tracing on each triangular surface element in sequence.
3. The method for modeling coupling scattering centers between metal targets based on single-station radar according to claim 2, wherein the targets are decomposed in step 1.2 by following the following principles:
(1) Decomposing the target to obtain a set of all entity parts;
(2) Taking abrupt change of normal vector of the target surface as a basis to decompose the component;
(3) The scattering mechanism of each part of the split target is clear, and one part does not have multiple scattering mechanisms.
4. A method for modeling a coupled scattering center between metallic targets based on a single-station radar as claimed in claim 3, wherein in step 1.3, all faces are subdivided into triangular face elements, and all triangular face elements and points are numbered as follows:
dividing all the faces into triangular surface elements, wherein the side length of the triangular surface elements is not less than ，/>To solve for the wavelength used in the scattering center amplitude using physical optics, all triangle bins and points are numbered and the part number of each triangle bin, the three vertex numbers of each triangle bin and the three-dimensional coordinates of each vertex are recorded.
5. The modeling method of coupling scattering centers between metal targets based on single-station radar according to claim 4, wherein in step 1.5, ray tracing is performed on each triangular surface element in sequence, specifically as follows:
step 1.5.1 determining an incident ray according to the radar incident direction Wherein/>And/>Respectively a pitch angle and an azimuth angle of the radar incidence direction under a spherical coordinate system;
step 1.5.2 incident rays for each triangle bin The first reflection point of the triangle is the geometric center of the triangle surface element contacted by the ray for the first time, and the reflected ray in the first reflection direction is calculated according to the geometric optics principle;
step 1.5.3, judging whether the reflected ray and other triangle surface elements have intersection points: if the intersection point exists, recording the part number of the triangle surface element and the number of the triangle surface element, calculating the direction of the next reflection, and entering a step 1.5.4; if the intersection point does not exist, directly entering the step 1.5.5;
step 1.5.4, judging whether the ejection times of the same incident ray are more than 6 times: if not more than 6 times, returning to the step 1.5.3; if more than 6 times, discard the ray traced this time, and jump to analysis of the next ray;
And 1.5.5, judging whether the included angle between the last reflection direction and the radar incidence direction is larger than 90 degrees, if not smaller than 90 degrees, discarding the ray traced at the time, jumping to analysis of the next ray, and if smaller than 90 degrees, considering the ray to be effective.
6. The method for modeling coupling scattering centers between metal targets based on single-station radar according to claim 5, wherein the projection size of equivalent points in each set in the distance direction and the azimuth direction is calculated in step 3, and the type of the scattering centers is established, specifically as follows:
Step 3.1, according to attribute scattering center theory, total scattering field of the target Expressed as:
(5)
wherein, Represents the/>The magnitude of the individual scattering centers; /(I)Representing radar frequency,/>Representing the radar center frequency; /(I)Representing a frequency dependent factor,/>The dependence of the scattering amplitude on the azimuth angle of the local scattering center is shown,；/>Represents the/>The length of the individual scattering centers; /(I)Representing the propagation velocity of the representative wave in vacuum; Representing the center azimuth of the radar; /(I) Represents the/>The location of the individual scattering centers;
step 3.2, traversing all rays in each set, respectively calculating projections of equivalent points in the set in azimuth directions and obtaining maximum and minimum values of the projections And/>The projection length of the set in the azimuth direction is/>=；
Step 3.3, setting a decision threshold in the azimuth directionWhen/></>When the scattering center model corresponds to the local scattering center;
When (when) When the scattering center model corresponds to the distributed scattering center, the dependence of the scattering center model on the azimuth angle is represented by the scattering center length/>And tilt angle/>Characterization.
7. The modeling method of coupling scattering centers between metal targets based on single-station radar according to claim 6, wherein in step 4, according to the type of the scattering centers, the relevant parameters of position, amplitude, length and frequency factor are respectively established to obtain an attribute scattering center model, and the target attribute feature extraction is realized, specifically as follows:
Step 4.1, length of local scattering center Frequency dependent factor/>All 0, position/>Sum amplitude/>Is calculated from the following formula:
(6)
wherein, Representing the position of equivalent points in a collection,/>Representing the magnitude of equivalent points in the set,/>；
Step 4.2, calculating projections of equivalent points in the set in the distance direction and obtaining maximum and minimum values of the projections respectively, wherein the length of the distributed scattering center in the azimuth direction is not 0 and the projections depend on the azimuth angle and rays in the setAndThe projection of the set in the distance direction is/>=/>；
Step 4.3, setting a threshold value for distance upWhen/></>Time,/>=/>Frequency dependent factor/>The values of (2) are as follows: when the equivalent point distribution is a curve,/>; When the equivalent point distribution is straight or hyperboloid,/>; When the equivalent point distribution is cylindrical,/>; When the equivalent point distribution is planar,/>; Position/>Sum amplitude/>The calculation formula of (2) is as follows:
(7)
Wherein the method comprises the steps of And/>The positions of the two points of minimum and maximum projection in azimuth;
When (when) >/>When the scattered field is:
(8)
Wherein the method comprises the steps of Representing the integration region,/>Representing the magnitude function at the integrated,/>Representing a dot/>The included angle between the scattering direction and the radar incidence direction, the scattered field is written as a fractional integral:
(9)
Wherein the method comprises the steps of Representation/>At/>N-th derivative of the position,/>Two endpoints of the integral area in azimuth;
Due to Only get/< < 1)In one aspect, the fringe field is reduced to:
(10)
The scattering center of the collection is equivalent to two points And/>Length/>、/>Frequency dependent factor/>、/>All 0, position/>、/>Sum amplitude is/>、/>The calculation formula of (2) is as follows:
(11)
wherein, And/>Coordinates of two endpoints of the integration region respectively; /(I)And/>Is the fringe field amplitude of the corresponding endpoint.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410167220.7A CN117706490B (en) | 2024-02-06 | 2024-02-06 | Method for modeling coupling scattering center between metal targets based on single-station radar |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410167220.7A CN117706490B (en) | 2024-02-06 | 2024-02-06 | Method for modeling coupling scattering center between metal targets based on single-station radar |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117706490A CN117706490A (en) | 2024-03-15 |
| CN117706490B true CN117706490B (en) | 2024-05-17 |
Family
ID=90153869
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410167220.7A Active CN117706490B (en) | 2024-02-06 | 2024-02-06 | Method for modeling coupling scattering center between metal targets based on single-station radar |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117706490B (en) |
Citations (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3099833A (en) * | 1958-03-17 | 1963-07-30 | Nat Res Dev | Method of radio scanning |
| US5608527A (en) * | 1995-03-08 | 1997-03-04 | Optical Dimensions, Llc | Apparatus and method for dynamic measurement of surface roughness |
| US6646593B1 (en) * | 2002-01-08 | 2003-11-11 | Science Applications International Corporation | Process for mapping multiple-bounce ghosting artifacts from radar imaging data |
| US7259715B1 (en) * | 2002-01-08 | 2007-08-21 | Science Applications International Corporation | Method for developing and using an image reconstruction algorithm for multipath scattering |
| EP2088449A1 (en) * | 2008-02-07 | 2009-08-12 | Saab Ab | Side lobe suppression |
| CN106324579A (en) * | 2016-10-31 | 2017-01-11 | 上海无线电设备研究所 | RCS extraction algorithm for strong scattering sources in target |
| CN106501777A (en) * | 2016-12-09 | 2017-03-15 | 北京环境特性研究所 | Scattering center source diagnostic method based on ray tracing |
| CN106680791A (en) * | 2017-01-04 | 2017-05-17 | 北京环境特性研究所 | Radar echo simulating method based on broadband sweep frequency data |
| RU2628997C1 (en) * | 2016-06-14 | 2017-08-24 | Российская Федерация, от имени которой выступает Министерство обороны Российской Федерации | Method of obtaining two-dimensional radar images of object at multi-frequency pulse sensing and inverse device synthesis with iterative distance reconciliation from equivalent antenna phase center to synthesization point |
| CN110083915A (en) * | 2019-04-22 | 2019-08-02 | 武汉大学 | A kind of high frequency region Radar Target Scatter central forward automation modeling method |
| CN110988868A (en) * | 2019-11-19 | 2020-04-10 | 南京理工大学 | Equivalent fractional order mode vortex electromagnetic wave generation and imaging method |
| CN111830500A (en) * | 2020-04-30 | 2020-10-27 | 南京理工大学 | Radar image simulation method of sea surface ship target based on improved SBR (sequencing batch reactor) rapid imaging technology |
| CN112711888A (en) * | 2021-01-08 | 2021-04-27 | 北京理工大学 | Combined scattering joint calculation method for bidirectional reflection distribution function and scattering center |
| CN112799063A (en) * | 2020-12-25 | 2021-05-14 | 北京环境特性研究所 | Double-station scattering center modeling method based on bounce ray method |
| CN114924245A (en) * | 2022-04-13 | 2022-08-19 | 北京理工大学 | Method for directly establishing scattering center model based on geometric structure |
| CN115128555A (en) * | 2022-06-06 | 2022-09-30 | 中国人民解放军空军工程大学 | Broadband echo signal simulation generation method based on scattering |
| CN115657017A (en) * | 2022-08-27 | 2023-01-31 | 南京理工大学 | Increment length diffraction theory ray imaging method for quick imaging of electrically large-size target |
| CN115754961A (en) * | 2022-11-16 | 2023-03-07 | 上海无线电设备研究所 | Echo generation method based on attribute scattering center near field correction model |
| CN116401844A (en) * | 2023-03-24 | 2023-07-07 | 武汉大学 | Radar target characteristic planning method based on forward model of scattering center |
| CN116430350A (en) * | 2023-06-14 | 2023-07-14 | 南京理工大学 | Modeling method and system for space dynamic group target attribute scattering center |
| CN116522589A (en) * | 2023-03-24 | 2023-08-01 | 武汉大学 | High-frequency modeling method for anisotropic medium target electromagnetic scattering |
| CN116577748A (en) * | 2023-03-13 | 2023-08-11 | 西安电子科技大学 | Scattering center parameter extraction method based on microwave photon signals |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6959869B2 (en) * | 1999-06-07 | 2005-11-01 | Metrologic Instruments, Inc. | Automatic vehicle identification (AVI) system employing planar laser illumination and imaging (PLIIM) based subsystems |
| EP2142921B1 (en) * | 2007-04-13 | 2020-01-01 | Centre Hospitalier De L'Universite de Montreal | Method of ultrasound scatterer characterization |
| US7652620B2 (en) * | 2007-08-09 | 2010-01-26 | Raytheon Company | RCS signature generation for closely spaced multiple objects using N-point models |
| US20230258770A1 (en) * | 2022-02-11 | 2023-08-17 | Aptiv Technologies Limited | Single-point Radar Cross Section Approaches for Radar Simulation |
-
2024
- 2024-02-06 CN CN202410167220.7A patent/CN117706490B/en active Active
Patent Citations (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3099833A (en) * | 1958-03-17 | 1963-07-30 | Nat Res Dev | Method of radio scanning |
| US5608527A (en) * | 1995-03-08 | 1997-03-04 | Optical Dimensions, Llc | Apparatus and method for dynamic measurement of surface roughness |
| US6646593B1 (en) * | 2002-01-08 | 2003-11-11 | Science Applications International Corporation | Process for mapping multiple-bounce ghosting artifacts from radar imaging data |
| US7259715B1 (en) * | 2002-01-08 | 2007-08-21 | Science Applications International Corporation | Method for developing and using an image reconstruction algorithm for multipath scattering |
| EP2088449A1 (en) * | 2008-02-07 | 2009-08-12 | Saab Ab | Side lobe suppression |
| RU2628997C1 (en) * | 2016-06-14 | 2017-08-24 | Российская Федерация, от имени которой выступает Министерство обороны Российской Федерации | Method of obtaining two-dimensional radar images of object at multi-frequency pulse sensing and inverse device synthesis with iterative distance reconciliation from equivalent antenna phase center to synthesization point |
| CN106324579A (en) * | 2016-10-31 | 2017-01-11 | 上海无线电设备研究所 | RCS extraction algorithm for strong scattering sources in target |
| CN106501777A (en) * | 2016-12-09 | 2017-03-15 | 北京环境特性研究所 | Scattering center source diagnostic method based on ray tracing |
| CN106680791A (en) * | 2017-01-04 | 2017-05-17 | 北京环境特性研究所 | Radar echo simulating method based on broadband sweep frequency data |
| CN110083915A (en) * | 2019-04-22 | 2019-08-02 | 武汉大学 | A kind of high frequency region Radar Target Scatter central forward automation modeling method |
| CN110988868A (en) * | 2019-11-19 | 2020-04-10 | 南京理工大学 | Equivalent fractional order mode vortex electromagnetic wave generation and imaging method |
| CN111830500A (en) * | 2020-04-30 | 2020-10-27 | 南京理工大学 | Radar image simulation method of sea surface ship target based on improved SBR (sequencing batch reactor) rapid imaging technology |
| CN112799063A (en) * | 2020-12-25 | 2021-05-14 | 北京环境特性研究所 | Double-station scattering center modeling method based on bounce ray method |
| CN112711888A (en) * | 2021-01-08 | 2021-04-27 | 北京理工大学 | Combined scattering joint calculation method for bidirectional reflection distribution function and scattering center |
| CN114924245A (en) * | 2022-04-13 | 2022-08-19 | 北京理工大学 | Method for directly establishing scattering center model based on geometric structure |
| CN115128555A (en) * | 2022-06-06 | 2022-09-30 | 中国人民解放军空军工程大学 | Broadband echo signal simulation generation method based on scattering |
| CN115657017A (en) * | 2022-08-27 | 2023-01-31 | 南京理工大学 | Increment length diffraction theory ray imaging method for quick imaging of electrically large-size target |
| CN115754961A (en) * | 2022-11-16 | 2023-03-07 | 上海无线电设备研究所 | Echo generation method based on attribute scattering center near field correction model |
| CN116577748A (en) * | 2023-03-13 | 2023-08-11 | 西安电子科技大学 | Scattering center parameter extraction method based on microwave photon signals |
| CN116401844A (en) * | 2023-03-24 | 2023-07-07 | 武汉大学 | Radar target characteristic planning method based on forward model of scattering center |
| CN116522589A (en) * | 2023-03-24 | 2023-08-01 | 武汉大学 | High-frequency modeling method for anisotropic medium target electromagnetic scattering |
| CN116430350A (en) * | 2023-06-14 | 2023-07-14 | 南京理工大学 | Modeling method and system for space dynamic group target attribute scattering center |
Non-Patent Citations (7)
| Title |
|---|
| Composite Scattering Analysis of the Ship on a Rough Surface Based on the Forward Parametric Scattering Center Modeling Method;Huang, Kai等;IEEE ANTENNAS AND WIRELESS PROPAGATION LETTERS;20200116;第18卷(第12期);全文 * |
| Off-grid One-dimensional Scattering Center Extraction Method Based on Compressed Sensing;Huang, Xiao等;2021 PHOTONICS & ELECTROMAGNETICS RESEARCH SYMPOSIUM (PIERS 2021);20220601;全文 * |
| Theoretical evaluation of phase center of artificial metal structure in NASA(JPL) TOPSAR data;Kim, JS等;IGARSS 2005: IEEE International Geoscience and Remote Sensing Symposium, Vols 1-8, Proceedings;20050101;全文 * |
| 基于散射中心模型的目标电磁特性智能生成网络研究;孙圣凯等;电波科学学报;20231231;38(5);全文 * |
| 复杂目标电磁散射中心建模及其在SAR目标解译识别中的应用;张磊;中国博士学位论文全文数据库基础科学辑;20230115;全文 * |
| 复杂目标近场电磁散射的建模方法;李静, 李向军;西安工业学院学报;20031230(第04期);全文 * |
| 雷达目标散射中心正向自动化建模方法研究与实现;何思远等;系统工程与电子技术;20210731;第43卷(第10期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117706490A (en) | 2024-03-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| He et al. | A forward approach to establish parametric scattering center models for known complex radar targets applied to SAR ATR | |
| Schuler et al. | Extraction of virtual scattering centers of vehicles by ray-tracing simulations | |
| Chiang et al. | Model-based classification of radar images | |
| CN112612024A (en) | Microwave array rapid imaging method | |
| Jianxiong et al. | Global scattering center model extraction of radar targets based on wideband measurements | |
| CN113376597A (en) | Complex terrain electromagnetic scattering rapid simulation method based on digital elevation map and GPU | |
| CN117706490B (en) | Method for modeling coupling scattering center between metal targets based on single-station radar | |
| CN116430350B (en) | Modeling method and system for space dynamic group target attribute scattering center | |
| Vatsavayi et al. | Efficient ISAR image classification using MECSM representation | |
| CN115205354B (en) | Phased array laser radar imaging method based on RANSAC and ICP point cloud registration | |
| Zhang et al. | Forward derivation and analysis for 3-D scattering center position of radar target | |
| Perotoni et al. | Numerical evaluation of an air-to-air missile radar cross section signature at X-band | |
| Perälä et al. | M-SPURT–Compressing the Target Characterization for a Fast Monostatic RCS Simulation | |
| Douchin et al. | SE-Workbench-RF: Performant and high-fidelity raw data generation for various radar applications | |
| Larsson | Near to far field transformation of RCS using a compressive sensing method | |
| CN113075659B (en) | Self-adaptive partitioning preprocessing method and system for grid model of flat scene | |
| CN111414801B (en) | Classification and identification method for electrically large non-cooperative target with uncertain shape | |
| Yoo et al. | A radar target DB construction method using 3D scattering centers | |
| Kobayashi | RCS Evaluation by Image-based Near-field to Far-field Transformation | |
| Zhang et al. | Three-dimensional scattering centers extraction of radar targets using high resolution techniques | |
| Wang et al. | Three‐dimensional point cloud reconstruction of inverse synthetic aperture radar image sequences based on back projection and iterative closest point fusion | |
| Guo et al. | Research in the Effect of Multi-Path Coupling Scattering in Radar Imaging | |
| Chen et al. | Scattering center modeling for low-detectable targets | |
| CN113567949A (en) | Method and device for associating scattering centers | |
| Liu et al. | Scattering centers diagnosis and parameters modification of the complex targets' geometry model based on the limited observed data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |