CN103530529B - Magnanimity target characteristic data quickly reads and precise interpolation method - Google Patents
Magnanimity target characteristic data quickly reads and precise interpolation method Download PDFInfo
- Publication number
- CN103530529B CN103530529B CN201310530355.7A CN201310530355A CN103530529B CN 103530529 B CN103530529 B CN 103530529B CN 201310530355 A CN201310530355 A CN 201310530355A CN 103530529 B CN103530529 B CN 103530529B
- Authority
- CN
- China
- Prior art keywords
- interpolation
- formula
- data
- frequency
- pitch
- 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.)
- Expired - Fee Related
Links
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
The open magnanimity target characteristic data of the present invention quickly reads and precise interpolation method, comprises: 1, calculate the Bistatic RCS of target, requires to store the derivative information of its single order, second order for the RCS at frequency end simultaneously；2, determine according to dbjective state and read useful Bistatic RCS scope；3, carry out the incident angle of pitch, incident orientation angle respectively, scatter the angle of pitch and scatter azimuthal interpolation；4, Bistatic RCS is carried out data format organization；5, dual station data complete for tissue are sent into Graphics Processing Unit, carry out the interpolation that frequency is one-dimensional；6, judged whether the incident angle of pitch, incident orientation angle, scattered the angle of pitch and scatter azimuthal interpolation；7, the interpolated data produced in Graphics Processing Unit is passed back central processing unit.The present invention uses lagrange-interpolation, it is achieved magnanimity target characteristic data quickly reads, it is to avoid data read on a large scale, improves interpolation precision, saves double counting.
Description
Technical field
The present invention relates to a kind of electromagnetic characteristic of scattering data processing technique, be specifically related to a kind of magnanimity target property number
According to quickly reading and precise interpolation method.
Background technology
Electrically large sizes, the broadband radar target performance data of labyrinth target calculate the most quite time-consuming, so
Setting up target characteristic data storehouse is the most general a kind of solution.With Bistatic RCS (Radar Cross
Section, RCS) as a example by, under a certain characteristic frequency, the incident angle of pitch, incident orientation angle and the scattering angle of pitch, scattering side
Parallactic angle, with 2 ° of sampling intervals, uses real part and the imaginary part of double-precision floating point access dual station RCS, and desired volume is 90 × 181 × 91
×181×8×2/1024/1024/1024G≈4G.Consider the target of Ku wave band (12GHz-18 GHz, with 0.1GHz for interval)
These data are stored internal memory (being obtained the RCS of other frequencies by these data interpolatings), need about 240G by property database
Internal memory.No matter on the reading time, or in memory requirements, these are all unallowed, are also unscientific.
Interpolation mainly includes interpolation and extrapolation, interpolation technique be generally divided into Lagrange's interpolation, newton polynomial interopolation, three
Secondary spline interpolation and Ai Er meter Te interpolation.Wherein Lagrange's interpolation uses the most universal.For more accurate interpolation, generally
Wish that the information near interpolation point is The more the better.In magnanimity target characteristic data interpolation, we often can be according to electromagnetism meter
Of calculation method itself, is not expending internal memory and on the premise of the time, is obtaining the numerical information near more interpolation point, and this is also slotting
Data basis is provided on the accuracy of value technology.More conventional direct interpolation can have greatly improved in precision.
Summary of the invention
The present invention provides a kind of magnanimity target characteristic data quickly to read and precise interpolation method, it is to avoid data extensive
Read, improve interpolation precision, save substantial amounts of double counting.
For achieving the above object, the present invention provides a kind of magnanimity target characteristic data quickly to read and precise interpolation method,
Being characterized in, the method includes the steps of:
Step 1, the Bistatic RCS of calculating target, require simultaneously for the RCS at frequency end
Store the derivative information of its single order, second order；
Step 2, determine according to dbjective state and read useful Bistatic RCS scope；
Step 3, carry out the incident angle of pitch, incident orientation angle respectively, scatter the angle of pitch and scatter azimuthal interpolation；
Step 4, Bistatic RCS is carried out data format organization；
Step 5, the complete dual station data of tissue are sent into Graphics Processing Unit, first carry out the interpolation that frequency is one-dimensional；
Step 6, judge whether the incident angle of pitch, incident orientation angle, scatter the angle of pitch and scatter azimuthal interpolation,
The most then jump to step 7, if it is not, then jump to step 4；
Step 7, pass the interpolated data produced in Graphics Processing Unit back central processing unit.
In above-mentioned steps 1, Bistatic RCS and derivative information thereof use binary storage.
In above-mentioned steps 2, Bistatic RCS scope uses the location function reading in c/c++ to carry out useful data
Quick obtaining.
In above-mentioned steps 3, Interpolation uses Lagrange's interpolation or polynomial interopolation.
When using Lagrange's interpolation in above-mentioned steps 4, data format organization method comprises:
5-D Lagrange's interpolation is the direct extension of 1-D Lagrange's interpolation, and its form shows such as formula (1):
(1)
In formula,,、、WithExpression formula is similar,；
If being first grouped according to frequency by the useful data read, on each frequency, data type of organization is too；
According in formula (1)Definition, all scattering datas on each frequencyIt is same；With
Time, for the respective element in each frequency matrix, the weighted sum operation carried out is identical.
The interpolation method using Lagrange's interpolation in above-mentioned steps 5 comprises:
Use the specific form-polynomial interopolation of Lagrange's interpolation, carry out the interpolation that frequency is one-dimensional, such as formula (2):
(2)
Wherein,For frequency variable,It it is polynomial exponent number.Under a certain state, it is known that, then undetermined constantNeed to meet with the matrix relationship of following formula (3):
(3)
In order to carry out the interpolation near end points more accurately, such as formula (4), (5), the single order at end points, second order can be led
Numerical information is applied in the matrix of formula (3), it may be assumed that
(4)
(5)
By latter two equation in formula (4), (5) alternate form (3), formula (6) can be obtained:
(6)
The coefficient matrix of formula (6) is generalized circular matrix deformation.
Magnanimity target characteristic data of the present invention quickly reads the broadband radar target with precise interpolation method and prior art
Performance data acquisition methods is compared, and has an advantage in that, the invention discloses one lagrange-interpolation fast and accurately, real
Show the quick reading of magnanimity target characteristic data, it is to avoid the extensive reading of data, improved interpolation precision, save substantial amounts of heavy
Multiple calculating.
Accompanying drawing explanation
Fig. 1 is that magnanimity target characteristic data of the present invention quickly reads and the flow chart of precise interpolation method；
Fig. 2 be magnanimity target characteristic data of the present invention quickly read with precise interpolation method in use derivative information hold
The schematic diagram of point interpolation；
Fig. 3 is the data schematic illustration of tissue that magnanimity target characteristic data of the present invention quickly reads with precise interpolation method；
Fig. 4 is that the metal ball that radius is 0.3m using Lagrange's interpolation in the present invention is shone at 10GHz uniform plane wave
Dual station RCS contrast schematic diagram under penetrating；
Fig. 5 is that to use the length of side of Lagrange's interpolation to be the rectangular metal trihedral angle of 0.3m in the present invention equal at 10GHz
Mono-static RCS contrast schematic diagram under even plane wave illumination.
Detailed description of the invention
Below in conjunction with accompanying drawing, further illustrate the specific embodiment of the present invention.
As it is shown in figure 1, the open a kind of magnanimity target characteristic data of the present invention quickly reads and precise interpolation method, the method
Comprise the steps of
Step 1, employing physical optics (Physical Optics, PO) or Shooting and bouncing rays, calculate the list station or double of target
Stand RCS (RCS).RCS at frequency end is required to store simultaneously the derivative information of its single order, second order.
As a example by dual station RCS interpolation, comprise five variablees (the incident angle of pitch, incident orientation angle, the scattering angle of pitch, scattering
Azimuth and irradiation frequency), the most front four variablees are all circumferential cycle, but irradiation frequency is different.Therefore in frequency two
All can be lower than centre position in interpolation near end and extrapolation accuracy.In order to improve interpolation precision, PO method can used
Calculate scattering time pre-stored frequency end points near the single order of scatter echo, second dervative, as shown in Figure 2.Because using Gordon
Integration method can be with Analytical Solution PO integration, so the calculating of all-order derivative does not expend the time.
Reading to accelerate magnanimity target characteristic data, when data volume is big, dual station RCS and derivative information thereof use binary system
Storage.
Step 2, determine according to dbjective state (angle of pitch, azimuth) and read useful dual station RCS scope.Use c/c++
In location function reading carry out the quick obtaining of useful data and specify original position and end position, i.e. determine dual station RCS
Scope.Avoid the extensive reading of data, but read available information therein.
About dbjective state, such as, in Simulation for Electronic Countermeasures, the flight speed of target, acceleration and positional information are all
Assume that known.Azimuth angle theta and the angle of pitch φ of now target thus can be obtained according to observation position.General need to be obtained
Take (θ ± 10 °, φ ± 10 °) EM scattering echo information.
Step 3, carry out the incident angle of pitch, incident orientation angle respectively, scatter the angle of pitch and scatter azimuthal interpolation.Interpolation
Form can use Lagrange's interpolation and polynomial interopolation.
Step 4, dual station RCS read is carried out data format organization.
As a example by Lagrange's interpolation, 5-D Lagrange's interpolation is the direct extension of 1-D Lagrange's interpolation, its shape
Formula shows such as formula (1):
(1)
In formula,,、、WithExpression formula is similar,。
If being first grouped according to frequency by the useful data read, on each frequency, data type of organization is too.
As it is shown on figure 3, according in formula (1)Definition, all scattering datas on each frequency
It is same, so can save substantial amounts of double counting.Simultaneously for the respective element in each frequency matrix, carry out
Weighted sum operation is consistent.The most single data manipulation is particularly well-suited to Graphics Processing Unit (Graphic
Processing Unit, GPU) speed-up computation.Good interpolation acceleration can be carried out based on above 2.
Step 5, the complete dual station data of tissue are sent into GPU, first carry out the interpolation that frequency is one-dimensional.This interpolation concrete
Whether form is at end points according to frequency determines.Frequency at non-end points uses the polynomial interopolation form under ordinary meaning to insert
Value, is in the frequency at end points and uses the polynomial interopolation form interpolation of deformation.It is explained in detail below.
As a example by Lagrange's interpolation, use the specific form-polynomial interopolation of Lagrange's interpolation, carry out frequency one-dimensional
Interpolation, such as formula (2):
(2)
Wherein,For frequency variable,It it is polynomial exponent number.Under a certain state, it is known that, then undetermined constantNeed to meet with the matrix relationship of following formula (3):
(3)
In order to carry out the interpolation near end points more accurately, such as formula (4), (5), the single order at end points, second order can be led
Numerical information is applied in the matrix of above-mentioned formula (3), it may be assumed that
(4)
(5)
By latter two equation in formula (4), (5) alternate form (3), formula (6) can be obtained:
(6)
The coefficient matrix of formula (6) be vandermonde (Vandermonde) matrix deformation, coefficient solve can be the most succinct obtain
Take.
Step 6, judge whether the incident angle of pitch, incident orientation angle, scatter the angle of pitch and scatter azimuthal interpolation,
The most then jump to step 7, if it is not, then jump to step 4；
Step 7, the interpolated data produced at GPU is passed back CPU (Central Processing Unit,
CPU).
As shown in Figure 4, it is shown that use the metal ball that radius is 0.3m of Lagrange's interpolation at 10GHz uniform plane wave
Dual station RCS schematic diagram under Zhao Sheing.
As shown in Figure 5, it is shown that use the length of side of Lagrange's interpolation to be the rectangular metal trihedral angle of 0.3m at 10GHz
Mono-static RCS schematic diagram under uniform plane wave irradiation.
Although present disclosure has been made to be discussed in detail by above preferred embodiment, but it should be appreciated that above-mentioned
Description is not considered as limitation of the present invention.After those skilled in the art have read foregoing, for the present invention's
Multiple amendment and replacement all will be apparent from.Therefore, protection scope of the present invention should be limited to the appended claims.
Claims (5)
1. a magnanimity target characteristic data quickly reads and precise interpolation method, it is characterised in that the method comprises following step
Rapid:
Step 1, the Bistatic RCS of calculating target, require to store for the RCS at frequency end simultaneously
Its single order, the derivative information of second order；
Step 2, determine according to dbjective state and read useful Bistatic RCS scope；
Step 3, carry out the incident angle of pitch, incident orientation angle respectively, scatter the angle of pitch and scatter azimuthal interpolation；
Step 4, Bistatic RCS is carried out data format organization；
Step 5, the complete dual station data of tissue are sent into Graphics Processing Unit, first carry out the interpolation that frequency is one-dimensional；
Step 6, judge whether the incident angle of pitch, incident orientation angle, scatter the angle of pitch and scatter azimuthal interpolation, if
It is then to jump to step 7, if it is not, then jump to step 4；
Step 7, pass the interpolated data produced in Graphics Processing Unit back central processing unit;
When using Lagrange's interpolation in described step 4, data format organization method comprises:
5-D Lagrange's interpolation is the direct extension of 1-D Lagrange's interpolation, and its form shows such as formula (1):
(1)
In formula,,、、WithExpression formula is similar,；
If being first grouped according to frequency by the useful data read, on each frequency, data type of organization is too；
According in formula (1)Definition, all scattering datas on each frequencyIt is same；Meanwhile,
For the respective element in each frequency matrix, the weighted sum operation carried out is identical.
2. magnanimity target characteristic data as claimed in claim 1 quickly reads and precise interpolation method, it is characterised in that described
In step 1, Bistatic RCS and derivative information thereof use binary storage.
3. magnanimity target characteristic data as claimed in claim 1 quickly reads and precise interpolation method, it is characterised in that described
In step 2, Bistatic RCS scope uses the location function reading in c/c++ to carry out the quick obtaining of useful data.
4. magnanimity target characteristic data as claimed in claim 1 quickly reads and precise interpolation method, it is characterised in that described
In step 3, Interpolation uses Lagrange's interpolation or polynomial interopolation.
5. magnanimity target characteristic data as claimed in claim 4 quickly reads and precise interpolation method, it is characterised in that described
The interpolation method using Lagrange's interpolation in step 5 comprises:
Use the specific form-polynomial interopolation of Lagrange's interpolation, carry out the interpolation that frequency is one-dimensional, such as formula (2):
(2)
Wherein,For frequency variable,It it is polynomial exponent number；
Under a certain state, it is known that, then undetermined constantNeed to meet and close with the matrix of following formula (3)
System:
(3)
In order to carry out the interpolation near end points more accurately, can be if formula (4), (5) are by the single order at end points, second derivative values
Information application is in the matrix of formula (3), it may be assumed that
(4)
(5)
By latter two equation in formula (4), (5) alternate form (3), formula (6) can be obtained:
(6)
The coefficient matrix of formula (6) is generalized circular matrix deformation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310530355.7A CN103530529B (en) | 2013-11-01 | 2013-11-01 | Magnanimity target characteristic data quickly reads and precise interpolation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310530355.7A CN103530529B (en) | 2013-11-01 | 2013-11-01 | Magnanimity target characteristic data quickly reads and precise interpolation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103530529A CN103530529A (en) | 2014-01-22 |
| CN103530529B true CN103530529B (en) | 2016-10-26 |
Family
ID=49932534
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310530355.7A Expired - Fee Related CN103530529B (en) | 2013-11-01 | 2013-11-01 | Magnanimity target characteristic data quickly reads and precise interpolation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103530529B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106872972B (en) * | 2017-03-31 | 2019-05-24 | 北京环境特性研究所 | Near space Electromagnetic Scattering of Target data capture method based on sextuple interpolation |
| CN109816590B (en) * | 2018-12-26 | 2023-03-14 | 呈像科技(北京)有限公司 | Image extrapolation processing method |
| CN116540219B (en) * | 2023-07-04 | 2023-09-22 | 北醒（北京）光子科技有限公司 | Method and device for correcting radar emergent light angle, storage medium and electronic equipment |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101587500A (en) * | 2008-05-23 | 2009-11-25 | 中国科学院电子学研究所 | Computer emulation method for sea-surface imaging of bistatic synthetic aperture radar |
-
2013
- 2013-11-01 CN CN201310530355.7A patent/CN103530529B/en not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101587500A (en) * | 2008-05-23 | 2009-11-25 | 中国科学院电子学研究所 | Computer emulation method for sea-surface imaging of bistatic synthetic aperture radar |
Non-Patent Citations (4)
| Title |
|---|
| A Potential Integral Equations Method for Electromagnetic Scattering by Penetrable Bodies;Philippe De Doncker.;《IEEE transactions on antennas and propagation》;20010731;第49卷(第7期);1037-1042 * |
| Application of Asymptotic Waveform Evaluation to Hybrid FE-BI-MLFMA for Fast RCS Computation Over a Frequency Band;Bi-Yi Wu,et al.;《IEEE transactions on antennas and propagation》;20130531;第61卷(第5期);2597-2604 * |
| Hermite插值结合FDTD法快速计算三维目标宽角度RCS;王立峰，等.;《航空学报》;20080731;第29卷(第4期);924-930 * |
| Interpolation/Extrapolation of Radar Cross-Section (RCS) Data in the Frequency Domain Using the Cauchy Method;Jie Yang,et al.;《IEEE transactions on antennas and propagation》;20071031;第55卷(第10期);2844-2851 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103530529A (en) | 2014-01-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104038295B (en) | A kind of distortion array antenna scattering property analytical method based on mechanical-electric coupling | |
| Jenn et al. | Wind turbine radar cross section | |
| CN103530469B (en) | A kind of computing method of corner reflector radar cross section | |
| CN103530529B (en) | Magnanimity target characteristic data quickly reads and precise interpolation method | |
| CN106291473A (en) | Nested type aerial array method to set up | |
| CN112733364B (en) | Foil cloud scattering rapid calculation method based on impedance matrix partitioning | |
| Sanchez-Olivares et al. | EM modelling of monostatic RCS for different complex targets in the near-field range: Experimental evaluation for traffic applications | |
| Niu et al. | ARCS: Active radar cross section for multi-radiator problems in complex EM environments | |
| CN104915465A (en) | Metal target transient electromagnetic scattering analysis method based on delay laguerre polynomials | |
| Pienaar et al. | RCS validation of asymptotic techniques using measured data of an electrically large complex model airframe | |
| Pienaar et al. | RCS results for an electrically large realistic model airframe | |
| Rotgerink et al. | Modelling of effects of nose radomes on radar antenna performance | |
| Ji et al. | Fast shadowing test algorithm based on target division by cubes | |
| Preobrazhensky et al. | The use of intellectual technologies for construction of models of objects with necessary average values scattering characteristics | |
| Zhang et al. | Investigation of the Shielding and Attenuation Effects of a Dynamical High-Density Chaff Cloud on the Signal Based on Voxel Splitting | |
| Lee | Investigation of the effects of target feature variations on ballistic missile RCS | |
| Fromentin‐Denozière et al. | Comparative study of miscellaneous methods applied to a benchmark, inlet scattering problem | |
| Wang et al. | EM scattering of a target above canyon/valley environment based on composite rough surface modeling method and modified SBR-FBSSA algorithm | |
| CN106156431B (en) | Target conductor electromagnetic scattering emulation mode based on nurbs surface modeling | |
| Hua et al. | High-frequency method for terahertz radar cross section of conductive targets in free space | |
| Bian et al. | Accurate Geometric Modeling and Scattering Characteristics Analysis of Complex Targets | |
| Ergül et al. | Benchmark solutions of large problems for evaluating accuracy and efficiency of electromagnetics solvers | |
| Hao et al. | Intelligent Target Design Based on Complex Target Simulation | |
| He et al. | A Wide-Angle Alternate Direction Implicit–Parabolic Equation Method for Electromagnetic Scattering from Electrically Large Targets | |
| Oswald et al. | Radar cross sections of objects with simulated defects using the parallel FDTD method |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20161026 Termination date: 20201101 |