CN106501777A - Scattering center source diagnostic method based on ray tracing - Google Patents

Scattering center source diagnostic method based on ray tracing Download PDF

Info

Publication number
CN106501777A
CN106501777A CN201611131829.0A CN201611131829A CN106501777A CN 106501777 A CN106501777 A CN 106501777A CN 201611131829 A CN201611131829 A CN 201611131829A CN 106501777 A CN106501777 A CN 106501777A
Authority
CN
China
Prior art keywords
ray
scattering
target
point
ejection
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.)
Granted
Application number
CN201611131829.0A
Other languages
Chinese (zh)
Other versions
CN106501777B (en
Inventor
闫华
董天齐
李胜
崔闪
李楠
侯兆国
王超
周伟江
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Environmental Features
Original Assignee
Beijing Institute of Environmental Features
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Environmental Features filed Critical Beijing Institute of Environmental Features
Priority to CN201611131829.0A priority Critical patent/CN106501777B/en
Publication of CN106501777A publication Critical patent/CN106501777A/en
Application granted granted Critical
Publication of CN106501777B publication Critical patent/CN106501777B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/021—Auxiliary means for detecting or identifying radar signals or the like, e.g. radar jamming signals

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Image Generation (AREA)

Abstract

The scattering center source diagnostic method based on ray tracing is disclosed, belongs to electromagnetic property field.Methods described comprises the steps：According to position of the scattering center in target gridding model, its coordinate range of place pixel in target three-dimensional image is determined；Infinite ray is launched to target gridding model, is obtained, three-dimensional path is postponed the ejection ray for falling into the coordinate range as effective rays；Based on the projectile path of every effective rays, the target surface area corresponding to bin that all effective rayses are passed through is originated as the scattering of scattering center.The present invention scatters the projectile path that modeling technique follows the trail of the ejection ray per root using high-frequency electromagnetic, the target surface area corresponding to bin that all effective rayses are passed through is originated as the scattering of scattering center, the scattering source diagnostic message for obtaining can be directly related to target gridding model,, more completely with intuitively, computational efficiency is high for diagnostic result.

Description

Scattering center source diagnostic method based on ray tracing
Technical field
The present invention relates to electromagnetic property field, the scattering center source diagnostic method more particularly to based on ray tracing.
Background technology
Hereinafter the background of related of the present invention is illustrated, but these explanations might not constitute the existing of the present invention Technology.
Theoretical calculation and experiment measurement show that, in high frequency region, the total electromagnetic scattering of target is may be considered by some offices Synthesized by electromagnetic scattering on portion position, the scattering source of these localities is commonly known as equivalent multi-scattering centers, or abbreviation Scattering center.Target scattering center is one of essential characteristic that target is scattered in high frequency region, and it is in target scattering characteristics analysis, mesh The fields such as mark Feature extraction and recognition, echo signal high-speed simulation are respectively provided with important using value.
In target scattering characteristics analysis, the scattering source of target can be diagnosed by scattering center.According in scattering The position of heart distribution, it can be determined that the position of main scattering source on objective body, and then the local for judging to cause target mainly to scatter Structure and scattering mechanism.Target scattering center is diagnosed, and can be provided for applications such as the reduction of target RCS, target physical reproductions important Technical support.
The distribution usually in a series of isolated point of the multi-scattering centers of realistic objective.But, according to Electromagnetic Scattering Theory, The scattering when formation of each scattering center of target has not only with target at the scattering center location point is relevant, it is also possible to and In target, the scattering of certain area is all related, for example, the trihedral main scattering center of three orthogonal flat board compositions be by Multiple reflections between the certain surface area of each flat board and formed.Therefore, target scattering center diagnosis can not provide scattering The complete information in source, and complete scattering source-information needs accurately to provide the target for producing each scattering center contribution The set at all scattering positions in surface.At present, the correlation for proposing to realize this complete scattering source analysis is not yet disclosed Technology.
Content of the invention
It is an object of the invention to proposing the scattering center source diagnostic method based on ray tracing, can be diagnosed to be complete Scattering source information, diagnostic result is more directly perceived, and computational efficiency is higher.
Scattering center source diagnostic method of the present invention based on ray tracing, comprises the steps：
S1, the position according to scattering center in target gridding model, determine its place pixel in target three-dimensional image Coordinate range；
S2, to target gridding model launch infinite ray, obtain per root launch ejection of the ray in target gridding model Path, three-dimensional path is postponed the ejection ray for falling into the coordinate range as effective rays；The projectile path includes bullet Penetrate the position coordinates of number of times, the bin number that passes through and catapult point；
S3, the projectile path based on every effective rays, the object table corresponding to bin that all effective rayses are passed through Originate as the scattering of the scattering center in face region.
Preferably, obtaining ejection projectile path of the ray in target gridding model per root in step S2 includes：For every Root launches ray：
The light reflection law of initial incidence angle and geometrical optics approach GO based on ejection ray, determines that ejection ray exists Incidence angle on each bin and angle of reflection in target gridding model；
According to ejection incidence angle and angle of reflection of the ray on each bin, determine per root launch ray ejection number of times, Through bin number and catapult point position coordinates.
Preferably, further include before step S1：The scattering center of target is obtained as follows：
S01, to target gridding model launch infinite ray, obtain per root launch ray of the ray in target gridding model Pipe data, including：Ejection number of times, the electric field value of the bin number, the position coordinates of catapult point and emergent ray that passes through；
S02, the true data according to the ejection ray per root, determine the scattering of each location of pixels in target three-dimensional image Intensity and three-dimensional path postpone the ray tube numbering for falling into the location of pixels；
S03, by scattering strength in target three-dimensional image exceed preset strength threshold value extreme point place pixel centre bit Put as scattering center, the scattering strength of extreme point is the amplitude of scattering center.
Preferably, step S01 includes：
The light reflection law of initial incidence angle and geometrical optics approach GO based on ejection ray, true according to following relation Ejection number of times of the fixed ejection ray per root in target gridding model, the bin number that passes through, the position coordinates of catapult point：θi= θr, θiFor launching incidence angle of the ray on bin, θrFor launching angle of reflection of the ray on bin；
Based on the electric field value that perfact conductor boundary condition determines indirect ray：EiRepresent and enter on bin The electric-field intensity of ray, ErThe electric-field intensity of indirect ray on bin is represented,Direction arrow for reflecting surface normal at pip Amount；
Far field integral formula based on physical optical method determine last time launch after emergent ray to recipient to electricity Field value：
In formula, r is the distance that target-based coordinate system origin scatters site position to far field, and unit is m；EsR () is dissipated for far field Radio field intensity, unit are V/m；J is imaginary unit；K is in-field wave number；For scattering direction unit vector；E0For last Incident ray electric-field intensity at secondary catapult point, unit is V/m；S0For the surface area that the ray of target surface illuminates；R ' is S0 In optional position vector, dS' be r ' place Line Integral infinitesimal；The Z-direction of target-based coordinate system is radar line of sight direction.
Preferably, step S02 includes：
According to the true data that ray is launched per root, determine described ejection ray each pixel in target three-dimensional image Scattering strength；
For each pixel in target three-dimensional image, postpone all ejections for falling into the location of pixels to three-dimensional path Ray is sued for peace in the scattering strength of the location of pixels, obtains the scattering strength of each pixel in target three-dimensional image.
Preferably, the true data to launching ray per root carries out frequency domain and angular domain Fourier integral, obtains ejection and penetrates The scattering strength of line each position in target three-dimensional image.
Preferably, in target three-dimensional image, the scattering strength of each position meets following relation：
In formula, Image3D (x, y, z) is the scattering strength at (x, y, z) place in target three-dimensional image, and unit is V/m；Z is Radial distance, x for orientation lateral separation, y be pitching to lateral separation, unit is m；ziIt is to click and enter relative to phase zero Ray postpones in the total distance of radial direction, xiLateral separation for orientation postpones, yiFor pitching to lateral separation postpone, single Position is m.
Preferably, the lateral separation of orientation postpones xiIt is that first time catapult point in orientation and last time catapult point are horizontal To the mean value of distance, pitching to lateral separation postpone yiFor first time catapult point and the last time catapult point horizontal stroke of pitching up Mean value to distance.
Preferably, step S03 includes：
For each peak point in target three-dimensional image, according to order from big to small respectively according to equation below to working as Scattering point near previous peaks point is iterated：
(Residual Image3D)n+1=(Residual Image3D)n-[Anh(x-xn,y-yn,z-zn)]
In formula, (Residual Image3D)nFor peak point (xn,yn,zn) place scattering strength, (Residual Image3D)n+1For the scattering strength at (x, y, z) place near present peak value point, AnFor the scattering strength of present peak value point, unit is V/m；N is iterations；
As (Residual Image3D)n+1During less than default iteration threshold, stop iteration, reject scattering strength and be less than The scattering point of the iteration threshold, obtains the scattering center after rejecting clutter.
The present invention scatters modeling technique using high-frequency electromagnetic, launches infinite ray to target gridding model and obtains per root bullet The projectile path of ray, three-dimensional path is postponed to fall into the ejection ray of scattering center place pixel region as effectively penetrating Line, the target surface area corresponding to bin that all effective rayses are passed through are originated as the scattering of scattering center, obtain Scattering source diagnostic message can be directly related to target gridding model, and, more completely with intuitively, computational efficiency is high for diagnostic result.
Description of the drawings
By the specific embodiment part of offer referring to the drawings, the features and advantages of the present invention will become more Easy to understand, in the accompanying drawings：
Fig. 1 is the schematic flow sheet of scattering center source diagnostic method of the present invention based on ray tracing；
Fig. 2 is the generating principle schematic diagram for illustrating true data of the present invention.
Specific embodiment
With reference to the accompanying drawings the illustrative embodiments of the present invention are described in detail.Illustrative embodiments are retouched State merely for the sake of demonstration purpose, and be definitely not to the present invention and its application or the restriction of usage.
The flow chart that Fig. 1 shows scattering center source diagnostic method of the present invention based on ray tracing.Scattered in order to determine The position of the heart is hit, present invention position in step sl first according to scattering center in target gridding model determines its institute Position in pixel in target three-dimensional image.As pixel is not a point, but an image-region, therefore in the present invention Coordinate range with pixel in target three-dimensional image characterizes position of the scattering center place pixel in target three-dimensional image.
Modeling technique is scattered using high-frequency electromagnetic in step S2, infinite ray is launched to target gridding model, obtained per root Ejection projectile path of the ray in target gridding model, including launching the position of number of times, the bin number that passes through and catapult point Coordinate.Bin refers to that a facet or curved surface in target gridding model, multiple bins constitute target gridding model, bin Size can be designed according to radar test precision and other diagnosis requirements.
When infinite ray is incided on the bin of target gridding model, launch on the bin, a part of light is sent out Life reflects to form indirect ray, and each indirect ray can be incided on another bin, as shown in Figure 2 as incident ray. In the projectile path of record ejection ray, can along ray projectile path by the bin that sequentially passes through of ejection ray per root by It is numbered according to ascending order, ejection ray often through a bin, launches the number of times increase of ray once.Parallel per root by following the trail of Ray incides the ejection number of times after target gridding model and the bin number of process can obtain the institute for passing through each bin There is the position coordinates of ejection ray and catapult point.For the ease of description, three-dimensional path is postponed to fall into scattering center institute by the present invention Pixel coordinate range ejection ray as effective rays.The formation of scattering center is that have through scattering center place bin Multiple reflections of all effective rayses on multiple bins which passes through and formed, therefore, can be based in step S3 per The target surface area corresponding to bin that all effective rayses are passed through by the projectile path of root effective rays is used as scattering center Scattering source.
The present invention determines position of the scattering center place pixel in target three-dimensional image in step sl, in step S2 really Determine the effective rays of the coordinate range that three-dimensional path postpones to fall into scattering center place pixel.Those skilled in the art should manage Solution, during actual diagnosis, the order of step S1 and step S2 can also be overturned.
The present invention scatters modeling technique using high-frequency electromagnetic, according to scattering center place pixel in target three-dimensional image What position and following the trail of obtained launches the projectile path of ray diagnosing the scattering source information of scattering center per root, and obtained dissipates Penetrating source diagnostic message can be directly related to target gridding model, and, more completely with intuitively, computational efficiency is high for diagnostic result.
In certain embodiments, ejection projectile path bag of the ray in target gridding model per root is obtained in step S2 Include：For the ejection ray per root：The light reflection law of initial incidence angle and geometrical optics approach GO based on ejection ray, really Incidence angle and angle of reflection of the ray on each bin are launched surely；According to ejection incidence angle and reflection of the ray on each bin Angle, determines the position coordinates of the ejection number of times, the bin number that passes through and catapult point that launch ray per root, as shown in Figure 2.
In order to accurately obtain each scattering center of target, the integrality of scattering center and its scattering source, step S1 is improved May further include before：The scattering center of target is obtained as follows：
S01, to target gridding model launch infinite ray, obtain per root launch ray of the ray in target gridding model Pipe data, including：Ejection number of times, the electric field value of the bin number, the position coordinates of catapult point and emergent ray that passes through；
S02, the true data according to the ejection ray per root, determine the scattering strength of each pixel in target three-dimensional image The ray tube for postponing to fall into the location of pixels with three-dimensional path is numbered；
S03, by scattering strength in target three-dimensional image exceed preset strength threshold value extreme point place pixel centre bit Put as scattering center, the scattering strength of extreme point is the amplitude of scattering center.
Fig. 2 shows the generating principle schematic diagram of true data of the present invention.Target is obtained using the method for the present invention to dissipate The heart is hit, and needs to obtain the true data that ray is launched per root.When the scattering source to scattering center is diagnosed, by anti- Can obtain to the true data for following the trail of the ejection ray per root and fall into all effective of scattering center place location of pixels scope All bins that ray and every effective rays are passed through, simplify data processing step, greatly reduce data processing amount.It is preferred that Ground, step S01 include：
The light reflection law of initial incidence angle and geometrical optics approach GO based on ejection ray, true according to following relation Ejection number of times of the fixed ejection ray per root in target gridding model, the bin number that passes through, the position coordinates of catapult point：θi= θr, θiFor launching incidence angle of the ray on bin, θrFor launching angle of reflection of the ray on bin；
Based on the electric field value that perfact conductor boundary condition determines indirect ray：EiRepresent and enter on bin The electric-field intensity of ray, ErThe electric-field intensity of indirect ray on bin is represented,Direction arrow for reflecting surface normal at pip Amount；
Far field integral formula based on physical optical method PO determine last time launch after emergent ray to recipient to Electric field value：
In formula, r is the distance that target-based coordinate system origin scatters site position to far field, and unit is m；EsR () is dissipated for far field Radio field intensity, unit are V/m；J is imaginary unit；K is in-field wave number；For scattering direction unit vector；E0For last Incident ray electric-field intensity at secondary catapult point, unit is V/m；S0For the surface area that the ray of target surface illuminates；R ' is S0 In optional position vector, dS' be r ' place Line Integral infinitesimal；The Z-direction of target-based coordinate system is radar line of sight direction.
In certain embodiments, step S02 includes：According to the true data that ray is launched per root, determine the ejection The scattering strength of ray each location of pixels in target three-dimensional image；For each location of pixels in target three-dimensional image, The all ejection rays for postponing to fall into three-dimensional path the location of pixels are sued for peace in the scattering strength of the location of pixels, obtain The scattering strength of each location of pixels in target three-dimensional image.Preferably, the true data for launching ray per root can be entered Line frequency domain and angular domain Fourier integral, obtain the scattering strength for launching ray each location of pixels in target three-dimensional image.
In the present invention, the scattering strength of in target three-dimensional image each location of pixels can be determined according to following relation：
In formula, Image3D (x, y, z) is the scattering strength at (x, y, z) place in target three-dimensional image, and unit is V/m；Z is Radial distance, x for orientation lateral separation, y be pitching to lateral separation, unit is m；ziIt is to click and enter relative to phase zero Ray postpones in the total distance of radial direction, xiLateral separation for orientation postpones, yiFor pitching to lateral separation postpone, single Position is m.
Preferably, the lateral separation of orientation postpones xiIt is that first time catapult point in orientation and last time catapult point are horizontal To the mean value of distance, pitching to lateral separation postpone yiFor first time catapult point and the last time catapult point horizontal stroke of pitching up Mean value to distance.
Obtain in target three-dimensional image after the scattering strength of each location of pixels, directly scattering strength can be less than etc. Reject in the coordinate position of preset strength threshold value, so as to quickly and easily obtain each scattering center of target.Due to three-dimensional image The neighbouring standard compliant ideal point scattering source of each peak point point spread function, therefore for target three-dimensional image in every Individual peak point, can be carrying out to the scattering point near present peak value point according to equation below according to order from big to small respectively Iteration：
(Residual Image3D)n+1=(Residual Image3D)n-[Anh(x-xn,y-yn,z-zn)]
In formula, (Residual Image3D)nFor peak point (xn,yn,zn) place scattering strength, (Residual Image3D)n+1For the scattering strength at (x, y, z) place near present peak value point, AnFor the scattering strength of present peak value point, unit is V/m；N is iterations；
As (Residual Image3D)n+1During less than default iteration threshold, stop iteration, reject scattering strength and be less than The scattering point of the iteration threshold, obtains the scattering center after rejecting clutter.Adopt Scattering center, it is nearby scattering center due to being imaged the extension for causing by misdescription to prevent peak point or peak point.
Different from diagnostic method of the tradition based on RCS data imagings or scattering centers extraction, the present invention is as a result of height Frequency electromagnetic scattering modeling technique, the scattering source diagnostic message for obtaining can be directly related to target geometrical model, thus result is more For complete with directly perceived, and computational efficiency is also higher.This method has important engineering application value, can be applicable to realize target Character control, using Shape design target RCS reduce and target physical reproduction etc. aspect.
Although with reference to illustrative embodiments, invention has been described, but it is to be understood that the present invention does not limit to The specific embodiment that Yu Wenzhong is described in detail and illustrated, in the case of without departing from claims limited range, this Art personnel can make various changes to the illustrative embodiments.

Claims (9)

1. based on ray tracing scattering center source diagnostic method, it is characterised in that comprise the steps：
S1, the position according to scattering center in target gridding model, determine its seat of place pixel in target three-dimensional image Mark scope；
S2, to target gridding model launch infinite ray, obtain per root launch projectile path of the ray in target gridding model, Three-dimensional path is postponed the ejection ray for falling into the coordinate range as effective rays；The projectile path includes ejection time The position coordinates of number, the bin number that passes through and catapult point；
S3, the projectile path based on every effective rays, the target surface area corresponding to bin that all effective rayses are passed through Originate as the scattering of the scattering center in domain.
2. scattering center as claimed in claim 1 is originated diagnostic method, it is characterised in that obtains in step S2 and penetrates per root Projectile path of the line in target gridding model includes：For the ejection ray per root：
The light reflection law of initial incidence angle and geometrical optics approach GO based on ejection ray, determines ejection ray in target Incidence angle on each bin and angle of reflection in grid model；
According to ejection incidence angle and angle of reflection of the ray on each bin, determine ejection number of times, the process that ray is launched per root Bin number and catapult point position coordinates.
3. scattering center as claimed in claim 1 is originated diagnostic method, it is characterised in that further include before step S1： The scattering center of target is obtained as follows：
S01, to target gridding model launch infinite ray, obtain per root launch ray tube number of the ray in target gridding model According to, including：Ejection number of times, the electric field value of the bin number, the position coordinates of catapult point and emergent ray that passes through；
S02, the true data according to the ejection ray per root, determine the scattering strength of each location of pixels in target three-dimensional image The ray tube for postponing to fall into the location of pixels with three-dimensional path is numbered；
Make in S03, the center of the extreme point place pixel that scattering strength in target three-dimensional image is exceeded preset strength threshold value For scattering center, the scattering strength of extreme point is the amplitude of scattering center.
4. scattering center as claimed in claim 3 is originated diagnostic method, it is characterised in that step S01 includes：
The light reflection law of initial incidence angle and geometrical optics approach GO based on ejection ray, determines according to following relation every Root launches ejection number of times of the ray in target gridding model, the bin number that passes through, the position coordinates of catapult point：θi=θr, θi For launching incidence angle of the ray on bin, θrFor launching angle of reflection of the ray on bin；
Based on the electric field value that perfact conductor boundary condition determines indirect ray：EiRepresent that incidence is penetrated on bin The electric-field intensity of line, ErThe electric-field intensity of indirect ray on bin is represented,Direction vector for reflecting surface normal at pip；
Far field integral formula based on physical optical method determine last time launch after emergent ray to recipient to electric field value：
E s ( r ) = - j ke - j k r 2 π r ∫ S 0 ( e ^ r · n ^ ) E 0 e - jkr ′ · e ^ r dS ′
In formula, r is the distance that target-based coordinate system origin scatters site position to far field, and unit is m；EsR () is far field scattering electric field Intensity, unit are V/m；J is imaginary unit；K is in-field wave number；For scattering direction unit vector；E0Launch for last time Incident ray electric-field intensity at point, unit is V/m；S0For the surface area that the ray of target surface illuminates；R ' is S0In appoint Meaning position vector, dS' are the Line Integral infinitesimal at r ' places；The Z-direction of target-based coordinate system is radar line of sight direction.
5. scattering center as claimed in claim 3 is originated diagnostic method, it is characterised in that step S02 includes：
According to the true data that ray is launched per root, determine described ejection ray each location of pixels in target three-dimensional image Scattering strength；
For each location of pixels in target three-dimensional image, postpone all ejections for falling into the location of pixels to three-dimensional path Ray is sued for peace in the scattering strength of the location of pixels, obtains the scattering strength of each location of pixels in target three-dimensional image.
6. scattering center as claimed in claim 5 is originated diagnostic method, it is characterised in that to launching the ray tube of ray per root Data carry out frequency domain and angular domain Fourier integral, obtain ejection ray scattering of each location of pixels in target three-dimensional image strong Degree.
7. scattering center as claimed in claim 6 is originated diagnostic method, it is characterised in that each position in target three-dimensional image Scattering strength meet following relation：
Im a g e 3 D ( x , y , z ) = Σ i , r a y s E s ( r ) h ( x - x i , y - y i , z - z i )
In formula, Image3D (x, y, z) is the scattering strength at (x, y, z) place in target three-dimensional image, and unit is V/m；Z is for radially Distance, x for orientation lateral separation, y be pitching to lateral separation, unit is m；ziIt is to penetrate relative to phase zero point incidence Line postpones in the total distance of radial direction, xiLateral separation for orientation postpones, yiFor pitching to lateral separation postpone, unit is m.
8. scattering center as claimed in claim 7 is originated diagnostic method, it is characterised in that the lateral separation of orientation postpones xi For first time catapult point in orientation and the mean value of last time catapult point lateral separation, pitching to lateral separation postpone yi Mean value for pitch up first time catapult point and last time catapult point lateral separation.
9. scattering center as claimed in claim 8 is originated diagnostic method, it is characterised in that step S03 includes：
For each peak point in target three-dimensional image, according to order from big to small respectively according to equation below to working as leading peak Scattering point near value point is iterated：
(Residual Image3D)n+1=(Residual Image3D)n-[Anh(x-xn,y-yn,z-zn)]
In formula, (Residual Image3D)nFor peak point (xn,yn,zn) place scattering strength, (Residual Image3D)n+1 For the scattering strength at (x, y, z) place near present peak value point, AnFor the scattering strength of present peak value point, unit is V/m；N is for repeatedly Generation number；
As (Residual Image3D)n+1During less than default iteration threshold, stop iteration, scattering strength is rejected less than described The scattering point of iteration threshold, obtains the scattering center after rejecting clutter.
CN201611131829.0A 2016-12-09 2016-12-09 Scattering center source diagnostic method based on ray tracing Active CN106501777B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611131829.0A CN106501777B (en) 2016-12-09 2016-12-09 Scattering center source diagnostic method based on ray tracing

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611131829.0A CN106501777B (en) 2016-12-09 2016-12-09 Scattering center source diagnostic method based on ray tracing

Publications (2)

Publication Number Publication Date
CN106501777A true CN106501777A (en) 2017-03-15
CN106501777B CN106501777B (en) 2019-03-26

Family

ID=58329593

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611131829.0A Active CN106501777B (en) 2016-12-09 2016-12-09 Scattering center source diagnostic method based on ray tracing

Country Status (1)

Country Link
CN (1) CN106501777B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108594196A (en) * 2018-04-28 2018-09-28 北京环境特性研究所 A kind of target scattering center parameter extracting method and device
CN110083915A (en) * 2019-04-22 2019-08-02 武汉大学 A kind of high frequency region Radar Target Scatter central forward automation modeling method
CN117706490A (en) * 2024-02-06 2024-03-15 南京理工大学 Method for modeling coupling scattering center between metal targets based on single-station radar

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103713284A (en) * 2012-09-28 2014-04-09 中国航天科工集团第二研究院二O七所 SBR and PO technology-based strong scattering center calculation method
CN104656078A (en) * 2013-11-25 2015-05-27 北京环境特性研究所 Sea surface ship method based on ray tracing
CN105974386A (en) * 2016-05-05 2016-09-28 乐山师范学院 Multistatic radar multi-target imaging positioning method
CN106093932A (en) * 2016-02-29 2016-11-09 中国科学院国家空间科学中心 A kind of high-resolution radar scatterometer of scanning beam

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103713284A (en) * 2012-09-28 2014-04-09 中国航天科工集团第二研究院二O七所 SBR and PO technology-based strong scattering center calculation method
CN104656078A (en) * 2013-11-25 2015-05-27 北京环境特性研究所 Sea surface ship method based on ray tracing
CN106093932A (en) * 2016-02-29 2016-11-09 中国科学院国家空间科学中心 A kind of high-resolution radar scatterometer of scanning beam
CN105974386A (en) * 2016-05-05 2016-09-28 乐山师范学院 Multistatic radar multi-target imaging positioning method

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108594196A (en) * 2018-04-28 2018-09-28 北京环境特性研究所 A kind of target scattering center parameter extracting method and device
CN108594196B (en) * 2018-04-28 2020-07-03 北京环境特性研究所 Method and device for extracting parameters of target scattering center
CN110083915A (en) * 2019-04-22 2019-08-02 武汉大学 A kind of high frequency region Radar Target Scatter central forward automation modeling method
CN110083915B (en) * 2019-04-22 2022-08-19 武汉大学 Forward automatic modeling method for radar target scattering center in high-frequency region
CN117706490A (en) * 2024-02-06 2024-03-15 南京理工大学 Method for modeling coupling scattering center between metal targets based on single-station radar
CN117706490B (en) * 2024-02-06 2024-05-17 南京理工大学 Method for modeling coupling scattering center between metal targets based on single-station radar

Also Published As

Publication number Publication date
CN106501777B (en) 2019-03-26

Similar Documents

Publication Publication Date Title
CN105388465B (en) Sea clutter simulation method based on wave spectrum model
CN104914415B (en) Single-pulse radar coherent jamming method based on target range profile template matching
US7652620B2 (en) RCS signature generation for closely spaced multiple objects using N-point models
CN104215951B (en) System and method for detecting low-speed small target under sea cluster background
CN105486180B (en) Laser-fuze near-field return wave power calculation method based on wave beam decomposition and partial irradiation
CN106501777A (en) Scattering center source diagnostic method based on ray tracing
CN104614713A (en) Radar echo signal simulator suitable for onboard radar system
CN105259541B (en) A kind of method of multistation radar anti-active cheating formula interference
EP2887092A1 (en) Computing radar cross section
CN101464511B (en) Working wave position decision method for satellite-loaded synthetic aperture radar
CN106707254B (en) A kind of analogy method and device of radar echo signal
CN103728598A (en) Method for restraining flight path deception jamming through remotely configured initiative and passive radar net
CN105607045A (en) Optimized station distribution method for netted radar unit under deceptive jamming
CN108594228A (en) Spatial target posture method of estimation based on ISAR image refocusings
CN203930058U (en) A kind of synthetic-aperture radar Area Objects Echo Signal Simulator
CN105738887A (en) Airborne radar clutter power spectrum optimization method based on Doppler channel division
CN103439698B (en) Method for obtaining radar scattering area
CN103760544B (en) The scattering centers extraction method and system of radar target
CN105467369A (en) Target echo simulation method and apparatus
CN105372637B (en) The arrowband cheating interference method to radar based on moving vehicle target fine motion characteristic
CN104182768A (en) Quality classification method for inverse synthetic aperture radar images
CN105116377A (en) An FDOA positioning method based on an HHT instantaneous energy spectrum
CN103632036B (en) The electromagnetic hotspot distribution diagram construction method of target and system
CN110489714A (en) A kind of calculation method of two-dimensional random scattering from rough surface statistical moment
CN103675777A (en) Airborne radar clutter analogy method and device based on fitting method

Legal Events

Date Code Title Description
C06 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