CN103487791A - RCS conversion method based on scattering center matrix - Google Patents

RCS conversion method based on scattering center matrix Download PDF

Info

Publication number
CN103487791A
CN103487791A CN201310436825.3A CN201310436825A CN103487791A CN 103487791 A CN103487791 A CN 103487791A CN 201310436825 A CN201310436825 A CN 201310436825A CN 103487791 A CN103487791 A CN 103487791A
Authority
CN
China
Prior art keywords
target
scattering
rcs
formula
matrix
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
CN201310436825.3A
Other languages
Chinese (zh)
Other versions
CN103487791B (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.)
Shanghai Radio Equipment Research Institute
Original Assignee
Shanghai Radio Equipment Research Institute
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 Shanghai Radio Equipment Research Institute filed Critical Shanghai Radio Equipment Research Institute
Priority to CN201310436825.3A priority Critical patent/CN103487791B/en
Publication of CN103487791A publication Critical patent/CN103487791A/en
Application granted granted Critical
Publication of CN103487791B publication Critical patent/CN103487791B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G01S7/41—Details 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Image Generation (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention relates to an RCS conversion method based on a scattering center matrix. The method comprises the steps of adopting a moment method or an iterative physical optical method to obtain the approximation relation between an induction source and an exposure field according to a set target rough geometric model, and setting the scattering center conversion matrix between the exposure field and the RCS with a target to be measured comprising an undetermined coefficient; according to the observation direction required to conduct RCS conversion, using the direction of an observation center as a fixed incident angle, and obtaining double-station scattering performance data under the different scattering angles of the target to be measured; adopting the least square method to determine the undetermined coefficient in the conversion matrix; achieving the RCS conversion of the target. The method is not limited in the high-frequency area, the suitable angle range is larger, the required data size is small, and engineering realization is easy.

Description

A kind of RCS conversion method based on the scattering center matrix
Technical field
The present invention relates to the electromagnetic characteristic of scattering data processing field, particularly a kind of RCS conversion method based on the scattering center matrix.
Background technology
Target RCS (RCS) conversion method mainly contains conversion method, the conversion method based on being scattering into picture and the convolution extrapolation algorithm etc. that distribute based on the enveloping surface tangential field at present.Wherein, theoretical foundation based on enveloping surface tangential field distribution conversion method is Huygens' principle, and the precision of its RCS conversion is higher, but the method requires two stations RCS data volumes of measurement very huge, need to measure the two stations RCS data under different incidence angles and different scattering angle, the Project Realization cost is larger; The scattering data that conversion method based on being scattering into picture requires is less, can directly utilize single station RCS width phase measuring system commonly used to obtain, Project Realization is easier to, but but the method requires the combination that target Approximate Equivalent to be measured is one group of scattering point, the local effect that comes from high frequency region due to scattering point is approximate, and need supposition scattering point parameter not with variations such as observation angle, distances, the scope of application of this conversion method is subject to larger restriction; The convolution extrapolation algorithm is mainly based on the physical optics model, and the measurement data that its conversion requires is less, but has also adopted the local effect of high frequency region approximate, and there is larger limitation too in its scope of application.As can be seen here, the RCS conversion method of the easier Project Realization of latter two exists more approximate on principle, it is the hypothesis that the high frequency region local effect is approximate and the scattering point parameter does not change with angle, distance etc., without loss of generality, take the linear system with interaction of electromagnetic wave and metal surface target here describes as example.
Suppose that the irradiation source distance is
Figure 2013104368253100002DEST_PATH_IMAGE002
the time target surface induction source distribution be
Figure 2013104368253100002DEST_PATH_IMAGE004
, the target surface exposure field is distributed as
Figure 2013104368253100002DEST_PATH_IMAGE006
, between target scattering induction source and exposure field, relation can be write as matrix form,
Figure 2013104368253100002DEST_PATH_IMAGE008
（1a）
Figure 2013104368253100002DEST_PATH_IMAGE010
（1b）
In formula:
Figure 2013104368253100002DEST_PATH_IMAGE012
for matrix of coefficients, when discrete point is enough close, this entry of a matrix element is only relevant to how much physical features of target and irradiation wave frequency etc., and irrelevant with irradiation distance, angle etc.Visible, generally the target surface induction source distributes not only relevant to the local irradiation field, also relevant to the exposure field at other position.
And existing based on being scattering into the RCS conversion methods such as picture, convolution extrapolation, it is mainly approximate according to the local effect of high frequency region, has ignored matrix of coefficients
Figure DEST_PATH_IMAGE012A
in non-diagonal regions element, be about to that formula (1) is approximate to be reduced to
Figure 2013104368253100002DEST_PATH_IMAGE014
（2）
In formula,
Figure 2013104368253100002DEST_PATH_IMAGE016
mean column vector;
Figure 2013104368253100002DEST_PATH_IMAGE018
the scattering point number that means target Equivalent; mean that the local induction source of target that each scattering point is corresponding distributes;
Figure 2013104368253100002DEST_PATH_IMAGE022
mean that local exposure field distributes;
Figure 2013104368253100002DEST_PATH_IMAGE024
mean local coefficient's matrix that each scattering point is corresponding,
=
Figure 2013104368253100002DEST_PATH_IMAGE028
（3）
This is its first approximate condition.Simultaneously, because scattered field is mainly distributed and determined by the target induction source, the scattering point coefficient also will be distributed and be determined by local induction source, according to formula (3), the scattering point coefficient will be subject to the impact that exposure field distributes: when irradiation source close together, target regional area size that scattering point is corresponding do not meet far field condition, the local irradiation field will be difficult to its average magnitude approximate representation, and the scattering point coefficient will change with observed range; When the observation angle scope is larger, the local irradiation field also will be difficult to its average magnitude approximate representation, and the scattering point coefficient will change with observation angle; This is second approximate condition.
As can be seen here, existing RCS conversion method is difficult to the practical situations of applicable Electromagnetic Scattering of Complex Target, needs the further more effective target RCS conversion method of research.
Summary of the invention
There is larger limitation in order to solve existing RCS conversion approximation technique, can't be applicable to the problem of the situations such as each locally coupled strong, scattering point parameter variation of target, the invention provides a kind of RCS conversion method based on the scattering center matrix.
Technical program of the present invention lies in providing a kind of RCS conversion method based on the scattering center matrix, the method comprises the steps:
A, by mapping, obtain the target geometric data, set up the rough geometric model of target, the roughly profile that wherein at least needs to comprise each parts of target and mutual geometric relationship;
B, according to the rough geometric model of target, adopt method of moment or iterative physical optics method to obtain the approximation relation between its induction source and exposure field, set up RCS that target to be measured comprises undetermined coefficient and the scattering center transition matrix between exposure field;
C, the observed ray of need being carried out to the RCS conversion, for being fixed into firing angle, obtain lower pair of station scattering properties data of the different scattering angle of target to be measured with the observation center position;
D, employing least square method are determined the undetermined coefficient in transition matrix;
E, utilize the scattering center matrix to carry out the target RCS conversion.
Alternatively, in steps A, further comprise:
When setting up the rough geometric model of target, according to position relationship and observed range between each parts of target, according to the far-field approximation condition, carry out the piecemeal division, wherein via formula 4, obtain each block size:
Figure DEST_PATH_IMAGE030
In formula,
Figure DEST_PATH_IMAGE032
for irradiation source and target component spacing;
Figure DEST_PATH_IMAGE034
for irradiating electromagnetic wavelength.
RCS conversion method as claimed in claim 2, is characterized in that,
In step B, further comprise:
The rough geometric model of target is carried out to mesh generation, obtain the induction source of each grid and the piecemeal approximation relation matrix between exposure field, formula 5:
Figure DEST_PATH_IMAGE036
In formula,
Figure DEST_PATH_IMAGE038
for target piecemeal sum;
Figure DEST_PATH_IMAGE040
mean the
Figure DEST_PATH_IMAGE042
the corresponding induction source of piecemeal and the
Figure DEST_PATH_IMAGE044
piecemeal relational matrix between the corresponding exposure field of piecemeal;
The piecemeal relational matrix of realistic objective is expressed as to formula 6:
By this piecemeal relational matrix element aggregation, and build the scattering center matrix, described scattering center entry of a matrix element is relation between corresponding scattering function and exposure field, with formula 7a and formula 7b, means:
Figure DEST_PATH_IMAGE048
Figure 2013104368253100002DEST_PATH_IMAGE050
In formula, for the scattering center matrix element;
Figure DEST_PATH_IMAGE054
for the scattered field observation station to target's center's distance;
Figure DEST_PATH_IMAGE056
for scattering angle;
Figure DEST_PATH_IMAGE032A
for irradiation source to target's center's distance;
Figure 2013104368253100002DEST_PATH_IMAGE059
for incident angle;
Figure 2013104368253100002DEST_PATH_IMAGE061
variation function for the scattering center matrix element;
Figure 2013104368253100002DEST_PATH_IMAGE063
undetermined coefficient for the scattering center matrix element;
Figure 2013104368253100002DEST_PATH_IMAGE065
be
Figure 2013104368253100002DEST_PATH_IMAGE067
the target surface area that individual induction source is corresponding;
Figure 2013104368253100002DEST_PATH_IMAGE069
for electromagnetic wavelength;
Figure 2013104368253100002DEST_PATH_IMAGE071
for local coefficient's matrix
Figure 2013104368253100002DEST_PATH_IMAGE073
element;
Figure DEST_PATH_IMAGE075
be
Figure DEST_PATH_IMAGE042A
the piecemeal grid is assembled the distance vector at center to it;
Figure 2013104368253100002DEST_PATH_IMAGE078
be the piecemeal grid is assembled the distance vector at center to it;
Figure DEST_PATH_IMAGE081
,
Figure DEST_PATH_IMAGE083
be respectively scattering wave vector and incident wave vector.
In step C, according to the observed range of setting, and near the set angle scope described incident direction chosen is during as the scattering direction, and the two station of the target of choosing Scattering Amplitude phase data sets are designated as
Figure DEST_PATH_IMAGE085
,
Figure DEST_PATH_IMAGE087
for scattering angle,
Figure DEST_PATH_IMAGE089
,
Figure DEST_PATH_IMAGE091
for choosing the scattering angle number, , and calculate the now variation function of each element of scattering center matrix according to described formula 7b.
In step D, the two stations scattering data collection that utilizes the scattering piecemeal relational matrix of target to be measured and choose, set up the system of linear equations that solves undetermined coefficient, and formula 8:
Figure DEST_PATH_IMAGE095
In formula, in formula 7b
Figure DEST_PATH_IMAGE099
's ,
Figure DEST_PATH_IMAGE002A
,
Figure DEST_PATH_IMAGE059A
shorthand way during for fixed value;
The recycling least square method is estimated undetermined coefficient, and formula 9:
Figure DEST_PATH_IMAGE103
In formula,
Figure DEST_PATH_IMAGE105
Figure DEST_PATH_IMAGE107
Figure DEST_PATH_IMAGE109
In above formula, mean column vector.
In step e, utilize the undetermined coefficient solve, calculate the variation function of each element of scattering center matrix under the observed range of new settings and angular range, and the conversion of according to formula 7a, 7b and following formula 10, calculating RCS:
Figure DEST_PATH_IMAGE111
Figure DEST_PATH_IMAGE113
Compared with prior art, the present invention has the following advantages:
The present invention mainly utilizes relation between the rough geometric model of target and induction source and exposure field, the scattering center transition matrix that comprises undetermined coefficient between establishing target RCS and exposure field, two stations scattering data that recycling is measured is determined the undetermined coefficient in matrix, and then the conversion of realize target RCS, have following advantage: (1), the inventive method have counted the coupling between target component, the stronger situation applicable to coupling between parts, and be not limited to high frequency region; (2) the inventive method has counted target geometric properties information roughly, and its applicable angular range is larger; (3) data volume that the inventive method requires is less, only requires the two station of target RCS point audio data under limited incident angle, is easy to Project Realization.
The accompanying drawing explanation
Fig. 1 is the process flow diagram of the method for the invention;
Fig. 2 is prolate dihedral angle piecemeal schematic diagram;
Fig. 3 is for the prolate dihedral angle of conversion two stations RCS data closely;
Fig. 4 is the remote single station of the dihedral angle RCS data that conversion is obtained.
Embodiment
By describing technology contents of the present invention, structural attitude in detail, being reached purpose and effect, below in conjunction with embodiment and coordinate accompanying drawing to be described in detail.
As shown in Figure 1, the RCS conversion method based on the scattering center matrix of the present invention comprises the steps:
(A), the rough profile of target is obtained and partitioning pretreatment
Obtain the roughly profile of each parts of target and mutual geometry site by mapping, ignore the details such as target surface roughness, deformation, set up rough object module; According to position relationship and observed range between each parts of target, according to the far-field approximation condition, carry out the piecemeal division again, each block size is got
Figure DEST_PATH_IMAGE030A
（4）
In formula,
Figure DEST_PATH_IMAGE032AA
for irradiation source and target component spacing; for irradiating electromagnetic wavelength.From (4) formula, the desirable smaller value of block size when the adjacent component close together, and adjacent component desirable higher value of block size when distant; Simultaneously, material behavior that as far as possible guarantees each piecemeal during division etc. is close.
(B), the scattering center matrix of RCS conversion builds
Utilize the rough geometric model of target, by theoretical modeling pattern, the target geometric model is carried out to mesh generation, then adopt the electromagnetic calculations such as method of moment, iterative physical optics to obtain the induction source of its each grid and the piecemeal approximation relation matrix between exposure field,
Figure DEST_PATH_IMAGE036A
（5）
In formula,
Figure DEST_PATH_IMAGE038A
for target piecemeal sum;
Figure DEST_PATH_IMAGE040A
mean the the corresponding induction source of piecemeal and the
Figure DEST_PATH_IMAGE044AA
piecemeal relational matrix between the corresponding exposure field of piecemeal, if with undetermined coefficient
Figure DEST_PATH_IMAGE063A
characterize the difference of realistic objective and the simple model of building, the piecemeal relational matrix of realistic objective can be written as
Figure DEST_PATH_IMAGE046A
（6）
Now, if by this piecemeal relational matrix element aggregation, can build the scattering center matrix according to the scattering center concept, matrix element is relation between corresponding scattering function and exposure field, can be written as
Figure DEST_PATH_IMAGE048A
（7a）
Figure DEST_PATH_IMAGE116
（7b）
In formula,
Figure DEST_PATH_IMAGE052A
for the scattering center matrix element; for the scattered field observation station to target's center's distance;
Figure DEST_PATH_IMAGE056A
for scattering angle;
Figure DEST_PATH_IMAGE032AAA
for irradiation source to target's center's distance;
Figure DEST_PATH_IMAGE059AA
for incident angle;
Figure DEST_PATH_IMAGE061A
variation function for the scattering center matrix element; undetermined coefficient for the scattering center matrix element;
Figure DEST_PATH_IMAGE065A
be the target surface area that individual induction source is corresponding;
Figure DEST_PATH_IMAGE069A
for electromagnetic wavelength;
Figure DEST_PATH_IMAGE071A
for local coefficient's matrix
Figure DEST_PATH_IMAGE073A
element; be
Figure DEST_PATH_IMAGE042AAA
the piecemeal grid is assembled the distance vector at center to it;
Figure DEST_PATH_IMAGE078A
be
Figure DEST_PATH_IMAGE044AAA
the piecemeal grid is assembled the distance vector at center to it; ,
Figure DEST_PATH_IMAGE083A
be respectively scattering wave vector and incident wave vector.
(C), two stations scattering data collection is chosen
The observed ray that conversion requires according to RCS, choose a certain observed range, elect the observation center position required as incident direction, the scattering direction choose near incident direction for example, than polarizers of big angle scope (being the angular range more than 50 ° or 50 ° in a preferred example), usually need comprise the main scattering direction of target to be measured, the two station of the target of choosing by these conditions Scattering Amplitude phase data set is designated as ,
Figure DEST_PATH_IMAGE087A
for scattering angle,
Figure DEST_PATH_IMAGE089A
, for choosing scattering angle number, requirement
Figure DEST_PATH_IMAGE093A
,
Figure DEST_PATH_IMAGE038AA
for target piecemeal sum.
(D), undetermined coefficient solves
The two stations scattering data collection that utilizes the scattering partitioned matrix relation of target to be measured and choose, set up the system of linear equations that solves undetermined coefficient,
Figure DEST_PATH_IMAGE118
（8）
In formula, because two stations scattering data of choosing is concentrated ,
Figure DEST_PATH_IMAGE002AA
,
Figure DEST_PATH_IMAGE059AAA
for fixed value, here by formula 7b be abbreviated as
Figure DEST_PATH_IMAGE097A
.
The recycling least square method is estimated undetermined coefficient,
Figure DEST_PATH_IMAGE103A
（9）
Figure DEST_PATH_IMAGE105A
Figure DEST_PATH_IMAGE107A
Figure DEST_PATH_IMAGE109A
In formula, mean column vector.
(E), RCS conversion
The undetermined coefficient that utilization solves, utilize formula (7) to realize under other observed range, view angle RCS(, conversion σ),
Figure DEST_PATH_IMAGE111A
（10）
Figure DEST_PATH_IMAGE113A
The prolate dihedral angle target of take is example, and the embodiment of the RCS transition matrix method based on the scattering center matrix of the present invention is described.Here launch and receive and all adopt vertical polarization, suppose that wave frequency is 10GHz, the RCS of known 0.8m distance target distributes, and requires conversion to obtain the target RCS at 1km place.
Step 1: the rough profile of target is obtained and partitioning pretreatment
First set up the rough geometric shape model of prolate dihedral angle target, as shown in Figure 2, Fig. 2 is the piecemeal schematic diagram of overlooking prolate dihedral angle under condition, here with the two sides of line segment approximate representation dihedral angle, the long length of side on its two sides is 0.36m, along the narrow length of side of direction of visual lines, is 0.02m.Simultaneously, according to observed range, far field condition and target signature, it is divided into to 5 zones, every length of side presses 0.09m, 0.09m and 0.18m divides,
Figure DEST_PATH_IMAGE120
,
Figure DEST_PATH_IMAGE122
,
Figure DEST_PATH_IMAGE124
.Wherein, l 1 with l 1 'be piecemeal 1 together; l 2 'for piecemeal 2; l 3 'for piecemeal 3; l 2 for piecemeal 4; l 3 for piecemeal 5.
Step 2: the scattering center matrix builds
Adopt the iterative physical optics method to calculate the relational matrix between target induction source and exposure field, then build the concrete scattering center matrix form of this target according to target characteristic ,
Figure DEST_PATH_IMAGE128
（11）
Step 3: two stations scattering data collection is chosen
Here two stations scattering properties data set of choosing is as follows: observed range is 0.8m, and incident angle is fixed as
Figure DEST_PATH_IMAGE130
, scattering angle in scope, press
Figure DEST_PATH_IMAGE134
interval is chosen.As shown in Figure 3.Calculate the now variation function of each element of scattering center matrix according to formula (7b) again
Figure DEST_PATH_IMAGE136
.
Step 4: undetermined coefficient solves
Two stations scattering data collection that utilization is chosen, adopt least square method to solve the undetermined coefficient in formula (11)
Figure DEST_PATH_IMAGE138
, its matrix form result is as follows:
Figure DEST_PATH_IMAGE140
（12）
Step 5: RCS conversion
According to the observed range that will change RCS and angular range, again calculate the variation function of each element of scattering center matrix , utilize formula (7) and formula (10) to calculate target lower single station RCS distribution at a distance, result as shown in Figure 4, its
Figure DEST_PATH_IMAGE144
in scope, result is less with the RCS error of directly calculating, and the validity of the inventive method has been described.
Although content of the present invention has been done detailed introduction by above preferred embodiment, will be appreciated that above-mentioned description should not be considered to limitation of the present invention.After those skilled in the art have read foregoing, for multiple modification of the present invention with to substitute will be all apparent.Therefore, protection scope of the present invention should be limited to the appended claims.

Claims (6)

1. the RCS conversion method based on the scattering center matrix, is characterized in that, the method comprises the steps:
A, by mapping, obtain the target geometric data, set up the rough geometric model of target, the roughly profile that wherein at least needs to comprise each parts of target and mutual geometric relationship;
B, according to the rough geometric model of target, adopt method of moment or iterative physical optics method to obtain the approximation relation between its induction source and exposure field, set up RCS that target to be measured comprises undetermined coefficient and the scattering center transition matrix between exposure field;
C, the observed ray of need being carried out to the RCS conversion, for being fixed into firing angle, obtain lower pair of station scattering properties data of the different scattering angle of target to be measured with the observation center position;
D, employing least square method are determined the undetermined coefficient in transition matrix;
E, utilize the scattering center matrix to carry out the target RCS conversion.
2. RCS conversion method as claimed in claim 1, is characterized in that,
In steps A, further comprise:
When setting up the rough geometric model of target, according to position relationship and observed range between each parts of target, according to the far-field approximation condition, carry out the piecemeal division, wherein via formula 4, obtain each block size:
Figure 2013104368253100001DEST_PATH_IMAGE002
In formula,
Figure 2013104368253100001DEST_PATH_IMAGE004
for irradiation source and target component spacing;
Figure 2013104368253100001DEST_PATH_IMAGE006
for irradiating electromagnetic wavelength.
3. RCS conversion method as claimed in claim 2, is characterized in that,
In step B, further comprise:
The rough geometric model of target is carried out to mesh generation, obtain the induction source of each grid and the piecemeal approximation relation matrix between exposure field, formula 5:
In formula, for target piecemeal sum; mean the
Figure DEST_PATH_IMAGE014
the corresponding induction source of piecemeal and the
Figure DEST_PATH_IMAGE016
piecemeal relational matrix between the corresponding exposure field of piecemeal;
The piecemeal relational matrix of realistic objective is expressed as to formula 6:
Figure DEST_PATH_IMAGE018
By this piecemeal relational matrix element aggregation, and build the scattering center matrix, described scattering center entry of a matrix element is relation between corresponding scattering function and exposure field, with formula 7a and formula 7b, means:
Figure DEST_PATH_IMAGE020
Figure DEST_PATH_IMAGE022
In formula,
Figure DEST_PATH_IMAGE024
for the scattering center matrix element;
Figure DEST_PATH_IMAGE026
for the scattered field observation station to target's center's distance;
Figure DEST_PATH_IMAGE028
for scattering angle;
Figure DEST_PATH_IMAGE004A
for irradiation source to target's center's distance;
Figure 2013104368253100001DEST_PATH_IMAGE031
for incident angle;
Figure 2013104368253100001DEST_PATH_IMAGE033
variation function for the scattering center matrix element;
Figure DEST_PATH_IMAGE035
undetermined coefficient for the scattering center matrix element;
Figure DEST_PATH_IMAGE037
be
Figure DEST_PATH_IMAGE039
the target surface area that individual induction source is corresponding;
Figure DEST_PATH_IMAGE041
for electromagnetic wavelength;
Figure DEST_PATH_IMAGE043
for local coefficient's matrix
Figure DEST_PATH_IMAGE045
element;
Figure DEST_PATH_IMAGE047
be
Figure DEST_PATH_IMAGE014A
the piecemeal grid is assembled the distance vector at center to it;
Figure DEST_PATH_IMAGE050
be
Figure DEST_PATH_IMAGE016A
the piecemeal grid is assembled the distance vector at center to it; ,
Figure DEST_PATH_IMAGE055
be respectively scattering wave vector and incident wave vector.
4. RCS conversion method as claimed in claim 3, is characterized in that,
In step C, according to the observed range of setting, and near the set angle scope described incident direction chosen is during as the scattering direction, and the two station of the target of choosing Scattering Amplitude phase data sets are designated as
Figure DEST_PATH_IMAGE057
,
Figure DEST_PATH_IMAGE059
for scattering angle,
Figure DEST_PATH_IMAGE061
, for choosing the scattering angle number,
Figure DEST_PATH_IMAGE065
, and calculate the now variation function of each element of scattering center matrix according to described formula 7b.
5. RCS conversion method as claimed in claim 4, is characterized in that,
In step D, the two stations scattering data collection that utilizes the scattering piecemeal relational matrix of target to be measured and choose, set up the system of linear equations that solves undetermined coefficient, and formula 8:
Figure DEST_PATH_IMAGE067
In formula,
Figure DEST_PATH_IMAGE069
in formula 7b
Figure DEST_PATH_IMAGE071
's ,
Figure DEST_PATH_IMAGE073
,
Figure DEST_PATH_IMAGE031A
shorthand way during for fixed value;
The recycling least square method is estimated undetermined coefficient, and formula 9:
Figure DEST_PATH_IMAGE076
In formula,
Figure DEST_PATH_IMAGE080
Figure DEST_PATH_IMAGE082
In above formula, mean column vector.
6. RCS conversion method as claimed in claim 5, is characterized in that,
In step e, utilize the undetermined coefficient solve, calculate the variation function of each element of scattering center matrix under the observed range of new settings and angular range, and the conversion of according to formula 7a, 7b and following formula 10, calculating RCS:
Figure DEST_PATH_IMAGE086
。
CN201310436825.3A 2013-09-24 2013-09-24 A kind of RCS conversion method based on scattering center matrix Expired - Fee Related CN103487791B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310436825.3A CN103487791B (en) 2013-09-24 2013-09-24 A kind of RCS conversion method based on scattering center matrix

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310436825.3A CN103487791B (en) 2013-09-24 2013-09-24 A kind of RCS conversion method based on scattering center matrix

Publications (2)

Publication Number Publication Date
CN103487791A true CN103487791A (en) 2014-01-01
CN103487791B CN103487791B (en) 2015-12-23

Family

ID=49828161

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310436825.3A Expired - Fee Related CN103487791B (en) 2013-09-24 2013-09-24 A kind of RCS conversion method based on scattering center matrix

Country Status (1)

Country Link
CN (1) CN103487791B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103760544A (en) * 2014-01-14 2014-04-30 北京环境特性研究所 Scattering center extraction method and system for radar target
CN104215949A (en) * 2014-08-26 2014-12-17 上海无线电设备研究所 Target RCS (radar cross section) calibration method based on automatic gain control voltage correction
CN106096267A (en) * 2016-06-08 2016-11-09 上海无线电设备研究所 A kind of EM scattering of cavities characteristic quick calculation method
CN106291496A (en) * 2016-07-26 2017-01-04 上海无线电设备研究所 Moving targets based on four path models and ground combination scattering modeling method
CN109100692A (en) * 2018-06-19 2018-12-28 西安电子科技大学 Rough surface and multiple target combination scattering emulation modes based on iterative physical optics

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101216556A (en) * 2007-12-27 2008-07-09 复旦大学 Electrically Large complex target and rugged face background composite electromagnetic scattering numerical value emulation method

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101216556A (en) * 2007-12-27 2008-07-09 复旦大学 Electrically Large complex target and rugged face background composite electromagnetic scattering numerical value emulation method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
梁子长 等: "广义RCS及近场电磁散射建模应用", 《上海航天》 *
梁子长 等: "角反射器的线性非均匀电磁波散射特性", 《制导与引信》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103760544A (en) * 2014-01-14 2014-04-30 北京环境特性研究所 Scattering center extraction method and system for radar target
CN103760544B (en) * 2014-01-14 2016-05-11 北京环境特性研究所 The scattering centers extraction method and system of radar target
CN104215949A (en) * 2014-08-26 2014-12-17 上海无线电设备研究所 Target RCS (radar cross section) calibration method based on automatic gain control voltage correction
CN106096267A (en) * 2016-06-08 2016-11-09 上海无线电设备研究所 A kind of EM scattering of cavities characteristic quick calculation method
CN106096267B (en) * 2016-06-08 2018-12-11 上海无线电设备研究所 A kind of EM scattering of cavities characteristic quick calculation method
CN106291496A (en) * 2016-07-26 2017-01-04 上海无线电设备研究所 Moving targets based on four path models and ground combination scattering modeling method
CN106291496B (en) * 2016-07-26 2019-03-08 上海无线电设备研究所 Moving target and ground combination scattering modeling method based on four path models
CN109100692A (en) * 2018-06-19 2018-12-28 西安电子科技大学 Rough surface and multiple target combination scattering emulation modes based on iterative physical optics

Also Published As

Publication number Publication date
CN103487791B (en) 2015-12-23

Similar Documents

Publication Publication Date Title
CN101692257B (en) Method for registering complex curved surface
CN103487791A (en) RCS conversion method based on scattering center matrix
CN102865862B (en) Measurement method of non-horizontalization free established station
CN102831751B (en) Road high-dangerous slope monitoring method based on double-camera imaging technology
CN103177193A (en) Rotationally symmetric target electromagnetic scattering rapid calculation method of thin-medium-coated metal
CN104502911A (en) Wall parameter estimation method of through-wall imaging radar
CN107632964A (en) A kind of plane GEOMAGNETIC FIELD downward continuation recurrence cosine transform method
Wu et al. The modern high frequency methods for solving electromagnetic scattering problems
CN107391880A (en) A kind of accurate efficient radome electrical property energy analysis method
CN102708257A (en) Method for fast determining tolerance of active phased array antenna structure
CN109215072A (en) Chaff cloud RCS acquisition methods based on tanh-sinh integration method
CN107180409A (en) A kind of generalized cylinder projecting method for bending matrix type object three-dimensional point cloud
CN103592627B (en) Determine method and the device of vibroseis combination center
Zimiao et al. A simple and precise calibration method for binocular vision
CN201576308U (en) Registration system for complex curved surface
CN103777229A (en) VSP observation system design method facing objective layer
CN104778293B (en) The volume of inhomogeneous dielectric object electromagnetic scattering divides Nystrom analysis methods
CN108051792B (en) Distribution description method of scattering source coupled between target and rough surface
CN104750960A (en) Method for rapidly extracting electrical property parameter of metal truss-type radome
CN113052918B (en) Method, device, medium and equipment for evaluating calibration error of opposite-pole binocular camera
CN104680586A (en) Method for fitting on ellipsoidal surface in spatial arbitrary position based on minimum area
Wu et al. Noniterative calibration of a camera lens with radial distortion
CN104636551A (en) Equivalent electromagnetic parameter deduction method for net-shaped reflector antenna metal mesh
CN103808646A (en) Image analyzer and method for measuring line fractal dimension of granular material
CN104914431B (en) A kind of focusing time delay computational methods for the non-homogeneous wall of known multilayer

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20151223

Termination date: 20210924

CF01 Termination of patent right due to non-payment of annual fee