CN112068101B - Target scattering separation method based on mode filtering - Google Patents

Target scattering separation method based on mode filtering Download PDF

Info

Publication number
CN112068101B
CN112068101B CN202010943150.1A CN202010943150A CN112068101B CN 112068101 B CN112068101 B CN 112068101B CN 202010943150 A CN202010943150 A CN 202010943150A CN 112068101 B CN112068101 B CN 112068101B
Authority
CN
China
Prior art keywords
wave
separated
mode
target
directional diagram
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
Application number
CN202010943150.1A
Other languages
Chinese (zh)
Other versions
CN112068101A (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.)
Xidian University
Original Assignee
Xidian University
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 Xidian University filed Critical Xidian University
Priority to CN202010943150.1A priority Critical patent/CN112068101B/en
Publication of CN112068101A publication Critical patent/CN112068101A/en
Application granted granted Critical
Publication of CN112068101B publication Critical patent/CN112068101B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/418Theoretical aspects
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/72Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables
    • G01N27/82Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables for investigating the presence of flaws
    • G01N27/83Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables for investigating the presence of flaws by investigating stray magnetic fields
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/411Identification of targets based on measurements of radar reflectivity
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/414Discriminating targets with respect to background clutter

Landscapes

  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Analytical Chemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Electrochemistry (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention provides a target scattering separation method based on mode filtering, which mainly solves the problem that the prior art can not accurately separate the scattering of each area of a target through single-frequency point measurement. The implementation scheme is as follows: acquiring field data of the whole space of the whole scatterer and preprocessing the field data; unfolding a far-field directional diagram on the whole space based on spherical waves to obtain spherical wave mode coefficients; and filtering high-order mode items in the spherical wave mode coefficient through a mode filter, and finally obtaining the filtered scattering field of each region of the whole scatterer. The invention only needs single-frequency point scattering data of the whole scatterer, so that the scattering field of each region of the whole scatterer can be separated from measurement scenes with extremely narrow bandwidth such as UHF/VHF, and meanwhile, the scattering measurement efficiency of each region is improved. Can be used for target subregion scattering measurement.

Description

Target scattering separation method based on mode filtering
Technical Field
The invention belongs to the technical field of electromagnetism, and particularly relates to a target scattering separation method which can be used for target subarea scattering measurement.
Background
The scattering measurement of the airplane is one of the main methods for judging the coating failure of the airplane body, the airplane may damage a coating at a certain part after executing a flight task once, so that the scattering performance changes, the airplane body can be divided into regions conventionally, the regions are detected for multiple times, and the scattering performance of each region before and after executing the task is compared respectively to judge the damaged region of the coating. However, this method may cause secondary damage during disassembly and transportation inspection. When scattering is measured outdoors, a large measurement error may be formed due to the harsh electromagnetic environment.
For example, a patent application with application publication number CN201910767035.0 entitled "a time domain gate transformation method and apparatus" discloses a method for separating and removing stray influence outside a target region, which is mainly based on a time domain gate technology, and is characterized by acquiring frequency domain RCS data of a target to be detected, converting the frequency domain RCS data into time domain RCS data according to a frequency domain time domain transformation method, acquiring a time domain gate function for time domain gating, applying the time domain gate function to the time domain RCS data to obtain time domain gating data, converting the time domain gating data into frequency domain data according to a time domain frequency domain transformation method, and performing renormalization processing on the frequency domain data to remove data aliasing.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provides a target scattering separation method based on mode filtering, which aims to solve the problem of accurately separating the scattering of each region of a target through single-frequency point measurement and improve the scattering measurement efficiency of each region.
In order to achieve the purpose, the technical scheme adopted by the invention comprises the following steps:
(1) Sampling scattered field data of the whole scatterer to obtain a far field directional diagram covering scattered field data of a spherical area of a part
Figure BDA0002674356370000011
To far field directional diagram
Figure BDA0002674356370000012
Coordinate translation operation and directional diagram zero filling operation are carried out to obtain a far field directional diagram on the whole space
Figure BDA0002674356370000021
Wherein theta is a pitch angle in a rectangular coordinate system of a measuring space,
Figure BDA0002674356370000022
measuring the azimuth angle in the space rectangular coordinate system;
(2) For far field directional diagram in the whole space
Figure BDA0002674356370000023
Based on the spherical wave expansion, obtaining a spherical wave mode expansion coefficient:
(2a) For far field directional diagram in the whole space
Figure BDA0002674356370000024
Based on spherical wave expansion, the expansion expression is:
Figure BDA0002674356370000025
wherein s represents an emergent wave, a represents an incident wave, m is a mode term mark having a value of 0 to infinity, n is a mode term mark having a value of 1 to infinity,
Figure BDA0002674356370000026
the m-n spherical wave expansion coefficient of the outgoing wave of the transverse electric wave TE,
Figure BDA0002674356370000027
is the m.n term spherical wave expansion coefficient of TE incident wave,
Figure BDA0002674356370000028
is an m-n term spherical wave expansion coefficient of the outgoing wave of the transverse magnetic wave TM,
Figure BDA0002674356370000029
is the m & n spherical wave expansion coefficient of TM incident wave,
Figure BDA00026743563700000210
m.n terms of TE emergent wave mode,
Figure BDA00026743563700000211
is a TE incident wave mode of m-n terms,
Figure BDA00026743563700000212
is an m.n term TM emergent wave mode,
Figure BDA00026743563700000213
is an m.n term TM incident wave mode, and j is an imaginary number symbol;
(2b) Cutting out the above
Figure BDA00026743563700000214
The A term and the B term in the spherical wave expansion expression of the method obtain a complete scattering information far-field directional diagram capable of containing a target scatterer to be separated
Figure BDA00026743563700000215
Figure BDA00026743563700000216
Wherein B =2 nextpow2(MN) ,MN=[kR 0 ]+ u, A =2B, u is a self-defined integer, R 0 Is the radius of the smallest sphere which takes the origin of the measurement coordinate as the center and surrounds the object to be separated;
(2c) For the complete scattering information far-field directional diagram containing the target scatterer to be separated
Figure BDA00026743563700000217
Performing expansion to obtain spherical wave expansion coefficient of TE wave
Figure BDA00026743563700000218
And the spherical wave expansion coefficient of TM wave
Figure BDA00026743563700000219
Figure BDA00026743563700000220
Figure BDA00026743563700000221
Wherein i = s or a, k is the wave number;
(3) Filtering out high-order mode items of spherical waves to obtain a far-field directional diagram of the target to be separated after filtering
Figure BDA00026743563700000222
(3a) Obtaining the number of truncation modes N of the spherical wave function 0 :N 0 =[kR 1 ]+v,
Wherein [ kR ] 1 ]Is kR 1 Integer rounded up, R 1 Taking the aperture center of the target to be separated as an origin and the radius of a minimum ball surrounding the target to be separated, wherein v is a self-setting integer;
(3b) Selecting a mode Filter (n) to Filter out the spherical wave expansion coefficient of the TE wave
Figure BDA0002674356370000031
And spherical wave expansion coefficient of TM wave
Figure BDA0002674356370000032
Obtaining the TE wave mode coefficient representing the target to be separated from the medium-high order mode coefficient
Figure BDA0002674356370000033
And TM wave mode coefficient
Figure BDA0002674356370000034
(3c) Intercepting the complete scattering information far field directional diagram containing the target scatterer to be separated
Figure BDA0002674356370000035
2N in spherical wave expansion expression 0 2 Term, obtaining far-field directional diagram of target to be separated after filtering
Figure BDA0002674356370000036
Figure BDA0002674356370000037
(4) Obtaining far field directional diagram of target to be separated in array
Figure BDA0002674356370000038
And the far-field directional diagram of the target to be separated after filtering is obtained
Figure BDA0002674356370000039
And far field pattern of object to be separated in array
Figure BDA00026743563700000310
And comparing, namely successfully separating the scattering field of the target to be separated from the scattering field of the whole scatterer when the horizontal section pattern curves or the vertical section pattern curves of the two can be matched.
Compared with the prior art, the invention has the following advantages:
according to the invention, single-frequency point scattered field data of the whole scatterer are sampled, and a scattered field directional diagram is expanded and filtered based on spherical waves, so that scattered fields of all parts in the whole scatterer can be obtained through a mathematical post-processing method, the requirement on a broadband measurement system is reduced, and the scattering measurement efficiency and precision of all parts are improved; meanwhile, the requirement on a broadband measurement system is reduced, so that the scattering fields of all parts of the whole scatterer can be separated for UHF/VHF frequency band measurement scenes.
Drawings
FIG. 1 is a flow chart of an implementation of an embodiment of the present invention;
FIG. 2 is a diagram of a three-unit metal patch array model in a scatterer to be measured according to the present invention;
fig. 3 is a two-station RCS cross-sectional pattern of a three-element metal patch array in a scatterer under test of the present invention;
fig. 4 is an in-array two-station RCS cross-sectional pattern of one cell in a three-cell metal patch array in a scatterer to be measured according to the present invention;
fig. 5 is a cross-sectional view of a dual station RCS of a cell isolated from a three-cell metal patch array in a scatterer under test using the present invention;
FIG. 6 is a comparison diagram of the total field pattern of the three-unit metal patch array in the scatterer to be measured, the in-array pattern of one unit, and the unit pattern separated from the array in the vertical azimuth plane;
fig. 7 is a comparison graph of the total field pattern of the three-element metal patch array coated with the wave-absorbing material in the scatterer to be measured of the present invention, the in-array pattern of one element, and the element pattern separated from the array in the vertical azimuth plane.
Detailed Description
Embodiments and effects of the present invention are described in further detail below with reference to the accompanying drawings.
Referring to FIG. 2, the scatterer used in this example is three identical metal patch units M, M placed in the XOY plane from the three-dimensional coordinate system XYZ 2 ,M 3 Formed patch array M = [ M = 1 ,M 2 ,M 3 ]Directly above it is a unit amplitude plane wave, M 2 The aperture center of the unit is located at the origin of coordinates of XYZ of a three-dimensional coordinate system, M 1 And M 3 The distance from the aperture center of the unit to the coordinate origin of the three-dimensional coordinate system XYZ is p, and the size of each patch unit is lambda multiplied by lambda, wherein lambda is the wavelength.
Referring to fig. 1, the implementation steps of the invention are as follows:
step 1, scattering field measurement data of a scatterer are obtained and preprocessed to obtain a far field pattern on the whole space
Figure BDA0002674356370000041
(1a) Sampling scattered field data of the whole scatterer, namely sampling scattered field data of a three-unit metal patch array M, wherein the sampling of the scattered field data needs to be performed by a scanning mode to perform area sampling:
the existing scanning modes include plane scanning, cylindrical scanning and spherical scanningThe existing scanning area has a near field, a quasi far field and a far field, the scanning mode of the embodiment adopts but is not limited to plane scanning, and the sampling area adopts but is not limited to the far field; obtaining the far field directional diagram of the three-unit metal patch array M
Figure BDA0002674356370000042
Wherein theta is a pitch angle in a rectangular coordinate system of a measuring space,
Figure BDA0002674356370000043
measuring the azimuth angle in the space rectangular coordinate system;
(1b) Will direct far field pattern
Figure BDA0002674356370000044
Translating to the aperture center of the patch unit to be separated to obtain a translated far-field directional diagram
Figure BDA0002674356370000045
Figure BDA0002674356370000046
Wherein, (x, y, z) is the caliber center coordinate of the patch unit to be separated,
Figure BDA0002674356370000047
representing the phase of coordinate translation, j is an imaginary number symbol, and k is a wave number;
(1c) For far field directional diagram obtained after translation
Figure BDA0002674356370000051
Carrying out directional diagram zero filling operation to obtain a far field directional diagram on the whole space
Figure BDA0002674356370000052
Step 2, aiming far field directional diagram on the whole space
Figure BDA0002674356370000053
Based on the spherical wave expansion, obtaining a spherical wave mode expansion coefficient:
(2a) For far field directional diagram in the whole space
Figure BDA0002674356370000054
Based on the spherical wave expansion, namely expanding into the superposition of a series of spherical wave mode terms, the expansion expression is as follows:
Figure BDA0002674356370000055
wherein s represents an emergent wave, a represents an incident wave, m is a mode term mark having a value of 0 to infinity, n is a mode term mark having a value of 1 to infinity,
Figure BDA0002674356370000056
the m-n spherical wave expansion coefficient of the outgoing wave of the transverse electric wave TE,
Figure BDA0002674356370000057
is the m.n term spherical wave expansion coefficient of TE incident wave,
Figure BDA0002674356370000058
is an m-n term spherical wave expansion coefficient of the outgoing wave of the transverse magnetic wave TM,
Figure BDA0002674356370000059
is the m & n spherical wave expansion coefficient of TM incident wave,
Figure BDA00026743563700000510
is a TE emergent wave mode of m-n terms,
Figure BDA00026743563700000511
is a TE incident wave mode of m.n terms,
Figure BDA00026743563700000512
is an m.n term TM emergent wave mode,
Figure BDA00026743563700000513
is m.n term TM incident wave mode;
the front A term and the front B term of a series of spherical wave mode terms in the expansion expression contain complete scattering information of the three-unit metal patch array M;
(2b) Cutting out the above
Figure BDA00026743563700000514
The first A term and the B term in the spherical wave expansion expression are used for obtaining a complete scattering information far field directional diagram capable of containing a to-be-three-unit metal patch array M
Figure BDA00026743563700000515
Figure BDA00026743563700000516
Wherein B =2 nextpow2(MN) ,MN=[kR 0 ]+ u, A =2B, u is a self-defined integer, R 0 The radius of the minimum sphere which takes the origin of the measurement coordinate as the center and surrounds the paster unit to be separated; when separating M 1 ,M 3 When the scattered field of the patch unit is present,
Figure BDA00026743563700000517
when separating M 2 When the scattered field of the patch unit is present,
Figure BDA00026743563700000518
(2c) According to the complete scattering information far-field directional diagram containing the target scatterer to be separated
Figure BDA00026743563700000519
Obtaining a spherical wave expansion coefficient of the TE wave
Figure BDA00026743563700000520
And the spherical wave expansion coefficient of TM wave
Figure BDA00026743563700000521
Figure BDA0002674356370000061
Figure BDA0002674356370000062
Wherein i = s or a.
Step 3, filtering high-order mode items of spherical waves to obtain a far field directional diagram of the target to be separated after filtering
Figure BDA0002674356370000063
(3a) Obtaining the number of truncated modes N of the spherical wave function 0
(3a1) Setting an integer v, wherein the integer v depends on the coordinate position relation of the source point and the field point and the required precision;
(3a2) Acquiring the radius of a minimum ball which takes the aperture center of the target to be separated as an origin and surrounds the target to be separated:
Figure BDA0002674356370000064
(3a3) Calculating the number of truncated modes N of the spherical wave function from the results of (3 a 1) and (3 a 2) 0
N 0 =[kR 1 ]+ v, wherein, [ kR 1 ]Is kR 1 An integer rounded up;
(3b) Selecting a mode Filter (n) to Filter out the spherical wave expansion coefficient of the TE wave
Figure BDA0002674356370000065
And the spherical wave expansion coefficient of TM wave
Figure BDA0002674356370000066
Higher order mode coefficient of medium:
existing modeThe Filter has a rectangular window, a cosine square window, a Hanning window, a Hamming window, a Blackman window, etc., and the mode Filter (n) of the present embodiment adopts but is not limited to a rectangular window, i.e., a rectangular window
Figure BDA0002674356370000067
TE wave mode coefficient representing target to be separated is obtained through rectangular window
Figure BDA0002674356370000068
And TM wave mode coefficient
Figure BDA0002674356370000069
Figure BDA00026743563700000610
Figure BDA00026743563700000611
(3c) Intercepting the spherical wave expansion expression in the step (2 b)
Figure BDA00026743563700000612
Spherical wave mode front 2N in 0 2 Item, obtaining far field directional diagram of the filtered patch unit to be separated
Figure BDA00026743563700000613
Figure BDA00026743563700000614
And 4, carrying out visual analysis on the far-field directional diagram of the to-be-separated patch unit obtained after filtering, and determining whether to finally separate the scattered field of the to-be-separated patch unit from the whole scattered field of the three-unit metal patch M.
(4.1) for the Patch Unit to be separated in the three-Unit Metal Patch array MThe scattered field data is sampled, in this example, the first element M in a three-element metal patch array M 1 Sampling the scattered field of the unit to obtain the far field directional diagram of the unit in the array
Figure BDA0002674356370000071
And the scattered field of the to-be-separated patch unit is used as a reference for judging whether the scattered field is separated from the whole scattered field of the three-unit metal patch array M or not;
(4.2) obtaining the far field directional diagram of the to-be-separated patch unit after filtering
Figure BDA0002674356370000072
With reference far field pattern
Figure BDA0002674356370000073
And comparing, namely successfully separating the scattered field of the to-be-separated patch unit from the integral scattered field of the three-unit metal patch array M when the vertical azimuth plane directional diagram curves of the three-unit metal patch array M can be matched.
The technical effects of the invention are further explained by combining simulation experiments as follows:
1. the experimental environment is as follows:
experimental software: FEKO + WinProp2018, MATLAB R2017a,
configuration of an experimental computer: intel (R) Core (TM) i7-8700K CPU 3.70GHZ, windows 10 (Pro)
The experimental conditions are as follows: a plane wave of unit amplitude with a frequency of 2GHZ is irradiated to a three-unit metal patch array M along a Z-axis direction of a three-dimensional coordinate system XYZ.
2. The experimental contents are as follows:
experiment 1, measuring a scattering field of a three-unit metal patch array M by using FEKO + WinProp2018 software to obtain a double-station RCS (radar cross section) directional diagram of the three-unit metal patch array M, as shown in FIG. 3;
experiment 2, the second unit M in the three-unit metal patch array M is processed by FEKO + WinProp2018 software 2 Measuring the scattered field of the unit to obtain M 2 The dual station RCS cross-sectional pattern of the cell, as shown in fig. 4;
experiment 3, using MATLAB R2017a software, the scattered field of the three-unit metal patch array M is expanded and mode-filtered based on spherical waves, and M separated from fig. 3 is obtained 2 The dual station RCS cross-sectional pattern of the cell, as shown in fig. 5;
the far-field patterns in the vertical azimuth plane obtained in experiment 1, experiment 2 and experiment 3 were compared, as shown in fig. 6.
Comparing fig. 4 and 5, the shadow distribution of the two is substantially the same, M in the three-unit metal patch array M in fig. 6 2 Far field directional diagram of unit in vertical azimuth plane and M separated from three-unit metal patch array M 2 The unit is jointed with the far field directional diagram curve on the vertical azimuth surface, the goodness of fit is high, and M is successfully matched 2 The scattered field of the unit is separated from the whole scattered field of the three-unit metal patch array M.
Experiment 4, M in three-unit metal patch array 1 ,M 3 The patch unit is coated with the wave-absorbing material, and the experiment 1, the experiment 2 and the experiment 3 are repeated, and the result is shown in fig. 7.
As can be seen from FIG. 7, M in the three-unit metal patch array coated with the wave-absorbing material 2 Far field directional pattern of unit in vertical azimuth plane and M separated from three-unit metal patch array coated with wave absorbing material 2 The unit is jointed with the far field directional diagram curve on the vertical azimuth surface, the goodness of fit is high, and M is successfully matched 2 The scattered field of the unit is separated from the integral scattered field of the three-unit metal patch array coated with the wave-absorbing material.

Claims (6)

1. A target scattering separation method based on mode filtering is characterized by comprising the following steps:
(1) Sampling scattered field data of the whole scatterer to obtain a far field directional diagram covering scattered field data of a spherical area of a part
Figure FDA0002674356360000011
To far field directional diagram
Figure FDA0002674356360000012
Coordinate translation operation and directional diagram zero filling operation are carried out to obtain a far field directional diagram on the whole space
Figure FDA0002674356360000013
Wherein theta is a pitch angle in a rectangular coordinate system of a measuring space,
Figure FDA0002674356360000014
measuring the azimuth angle in the space rectangular coordinate system;
(2) For far field directional diagram in the whole space
Figure FDA0002674356360000015
Based on the spherical wave expansion, obtaining a spherical wave mode expansion coefficient:
(2a) For far field directional diagram in the whole space
Figure FDA0002674356360000016
Based on spherical wave expansion, the expansion expression is:
Figure FDA0002674356360000017
wherein s represents an emergent wave, a represents an incident wave, m is a mode term mark having a value of 0 to infinity, n is a mode term mark having a value of 1 to infinity,
Figure FDA0002674356360000018
the m-n spherical wave expansion coefficient of the outgoing wave of the transverse electric wave TE,
Figure FDA0002674356360000019
is the m.n term spherical wave expansion coefficient of TE incident wave,
Figure FDA00026743563600000110
is an m-n term spherical wave expansion coefficient of the outgoing wave of the transverse magnetic wave TM,
Figure FDA00026743563600000111
is the m.n spherical wave expansion coefficient of the TM incident wave,
Figure FDA00026743563600000112
m.n terms of TE emergent wave mode,
Figure FDA00026743563600000113
is a TE incident wave mode of m.n terms,
Figure FDA00026743563600000114
is an m.n term TM emergent wave mode,
Figure FDA00026743563600000115
is an m.n term TM incident wave mode, and j is an imaginary number symbol;
(2b) Cutting out the above
Figure FDA00026743563600000116
The A term and the B term in the spherical wave expansion expression of the method are used for obtaining a complete scattering information far-field directional diagram capable of containing a target scatterer to be separated
Figure FDA00026743563600000117
Figure FDA00026743563600000118
Wherein B =2 nextpow2(MN) ,MN=[kR 0 ]+ u, A =2B, u is a self-defined integer, R 0 Is the radius of the smallest sphere which takes the origin of the measurement coordinate as the center and surrounds the object to be separated;
(2c) For the complete scattering information far-field directional diagram containing the target scatterer to be separated
Figure FDA00026743563600000119
Processing to obtain spherical wave expansion coefficient of TE wave
Figure FDA00026743563600000120
And spherical wave expansion coefficient of TM wave
Figure FDA00026743563600000121
Figure FDA0002674356360000021
Figure FDA0002674356360000022
Wherein i = s or a, k is the wave number;
(3) Filtering out spherical wave high-order mode items to obtain a far field directional diagram of the target to be separated after filtering
Figure FDA0002674356360000023
(3a) Obtaining the number of truncation modes N of the spherical wave function 0 :N 0 =[kR 1 ]+v,
Wherein, [ kR ] 1 ]Is kR 1 Integer rounded up, R 1 Taking the aperture center of the target to be separated as an origin and the radius of a minimum ball surrounding the target to be separated, wherein v is a self-setting integer;
(3b) Selecting a mode Filter (n) to Filter out the spherical wave expansion coefficient of the TE wave
Figure FDA0002674356360000024
And the spherical wave expansion coefficient of TM wave
Figure FDA0002674356360000025
Obtaining the TE wave mode coefficient representing the target to be separated from the medium-high order mode coefficient
Figure FDA0002674356360000026
And TM wave mode coefficient
Figure FDA0002674356360000027
(3c) Intercepting the complete scattering information far field directional diagram containing the target scatterer to be separated
Figure FDA0002674356360000028
2N in spherical wave expansion expression 0 2 Term, obtaining far field pattern of the filtered object to be separated
Figure FDA0002674356360000029
Figure FDA00026743563600000210
(4) Obtaining far field patterns of objects to be separated in the whole scatterer
Figure FDA00026743563600000211
And the far field directional diagram of the target to be separated after filtering
Figure FDA00026743563600000212
Far field pattern in array with object to be separated
Figure FDA00026743563600000213
And comparing, namely successfully separating the scattering field of the target to be separated from the scattering field of the whole scatterer when the horizontal section pattern curves or the vertical section pattern curves of the two can be matched.
2. The method of claim 1, wherein the scattered field data of the entire scatterer is sampled in (1) in a scanning manner selected from one of a plane, a cylinder, and a sphere, and the sampled region is selected from one of a near field, a quasi-far field, and a far field.
3. The method of claim 1, wherein (1) the far-field pattern is aligned
Figure FDA0002674356360000031
Coordinate translation operation is carried out by converting far-field directional diagram
Figure FDA0002674356360000032
Translating to the aperture center of the target to be separated to obtain a translated scattered field directional diagram
Figure FDA0002674356360000033
Is represented as follows:
Figure FDA0002674356360000034
wherein (x, y, z) is the aperture center coordinate of the object to be separated,
Figure FDA0002674356360000035
denotes the phase of the coordinate shift, j is the sign of the imaginary number, and k is the wave number.
4. The method of claim 1, wherein the mode Filter (n) in (3 b) is selected from any one of a rectangular window, a cosine square window, a hanning window, a hamming window, and a blackman window.
5. The method according to claim 1, wherein the TE wave mode coefficients characterizing the objects to be separated are obtained in (3 b)
Figure FDA0002674356360000036
Is represented as follows:
Figure FDA0002674356360000037
wherein Filter (n) is the selected mode Filter,
Figure FDA0002674356360000038
is the spherical wave expansion coefficient of the TE wave.
6. The method according to claim 1, wherein TM wave mode coefficients characterizing the objects to be separated are obtained in (3 b)
Figure FDA0002674356360000039
Is represented as follows:
Figure FDA00026743563600000310
wherein Filter (n) is the selected mode Filter,
Figure FDA00026743563600000311
is the spherical wave expansion coefficient of the TM wave.
CN202010943150.1A 2020-09-09 2020-09-09 Target scattering separation method based on mode filtering Active CN112068101B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010943150.1A CN112068101B (en) 2020-09-09 2020-09-09 Target scattering separation method based on mode filtering

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010943150.1A CN112068101B (en) 2020-09-09 2020-09-09 Target scattering separation method based on mode filtering

Publications (2)

Publication Number Publication Date
CN112068101A CN112068101A (en) 2020-12-11
CN112068101B true CN112068101B (en) 2023-03-14

Family

ID=73663243

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010943150.1A Active CN112068101B (en) 2020-09-09 2020-09-09 Target scattering separation method based on mode filtering

Country Status (1)

Country Link
CN (1) CN112068101B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113589051A (en) * 2021-08-04 2021-11-02 西安电子科技大学 Clutter suppression field electromagnetic measurement device and method
CN114488040B (en) * 2021-12-31 2024-09-06 西安电子科技大学 Phased array amplitude inversion method based on mode filtering

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1621857A (en) * 2003-11-24 2005-06-01 中国科学院声学研究所 Method and system for measuring carrier to bottom velocity by correlated speed measuring sound radar
CN106526562A (en) * 2016-10-17 2017-03-22 西北工业大学 Antenna directional pattern-based method for correcting RCS of calculation target
CN106646395A (en) * 2016-09-30 2017-05-10 西安电子科技大学 Radar echo deduction method for flight target
CN110632399A (en) * 2019-09-19 2019-12-31 电子科技大学 Correction method for spherical near-field measurement data and antenna directional pattern measurement method

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2526783A (en) * 2014-05-23 2015-12-09 Toshiba Res Europ Ltd Method for localising scattering elements in a 3D environment

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1621857A (en) * 2003-11-24 2005-06-01 中国科学院声学研究所 Method and system for measuring carrier to bottom velocity by correlated speed measuring sound radar
CN106646395A (en) * 2016-09-30 2017-05-10 西安电子科技大学 Radar echo deduction method for flight target
CN106526562A (en) * 2016-10-17 2017-03-22 西北工业大学 Antenna directional pattern-based method for correcting RCS of calculation target
CN110632399A (en) * 2019-09-19 2019-12-31 电子科技大学 Correction method for spherical near-field measurement data and antenna directional pattern measurement method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Efficient RCS Computation Over a Broad Frequency Band Using Subdomain MoM;XING WANG等;《IEEE Access》;20200214;全文 *
阵列天线辐射和散射分析与控制方法研究;张帅;《中国博士学位论文全文数据库》;20130415;全文 *

Also Published As

Publication number Publication date
CN112068101A (en) 2020-12-11

Similar Documents

Publication Publication Date Title
CN112068101B (en) Target scattering separation method based on mode filtering
CN103777186B (en) Decompose and moving target near field radar return characteristic computing method under local irradiation based on wave beam
CN105911533B (en) A kind of down-sampled fast scanning method of three-dimensional imaging based on flat scanning structure
Satterwhite Diffraction by a quarter plane, the exact solution, and some numerical results
CN108959806B (en) Equivalent radiation modeling method based on spherical surface near-field measurement and spherical mode source
CN113376597A (en) Complex terrain electromagnetic scattering rapid simulation method based on digital elevation map and GPU
CN108646228A (en) A kind of Radar Cross Section Calculating based on PO methods
CN110274920B (en) Method for analyzing transient characteristics of metal medium target by body surface subdivision bounce ray method
Culotta-López et al. Radiation center estimation from near-field data using a direct and an iterative approach
CN111460900B (en) Complex electromagnetic environment construction equivalence quantitative evaluation method
CN111735996A (en) Multipath interference suppression method and device for mathematical wave-absorbing darkroom construction
Tian et al. Examination of the effectiveness of mode orthogonalisation and filtering for scattering suppression in antenna measurements through computational electromagnetic simulation
Fang et al. Visualization and quantitative evaluation of delamination defects in GFRPs via sparse millimeter-wave imaging and image processing
CN116401844A (en) Radar target characteristic planning method based on forward model of scattering center
CN115524541A (en) Arbitrary polarization line antenna array directional diagram measuring method based on quasi-far-field mode filtering
CN115169170A (en) Composite target scattering semi-analytic rapid calculation method based on non-uniform grid model
Meschino et al. A practical guide on using SPOT-GPR, a freeware tool implementing a SAP-DoA technique
Gregson et al. Examination of the effectiveness of FF MARS in a CATR through electromagnetic simulation
Le Goff et al. A novel and innovative near field system for testing radomes of commercial aircrafts
Wang et al. Edge diffraction in NURBS-UTD method
KR102561279B1 (en) Sub-dimensional three dimension near-field to far-field radar cross section transformation method
Tian et al. Examination of spherical antenna far-field scattering suppression through electromagnetic simulation
Hu et al. Efficient wideband MRCS simulation for radar HRRP target recognition based on MSIB and PCA
CN117147977B (en) Test quality estimation method for unmanned aerial vehicle external field composite scattering measurement
Craeye et al. Modulated oscillations appearing in the scan impedance of a finite phased array

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