CN116008925A - Improved target radar sectional area estimation algorithm - Google Patents
Improved target radar sectional area estimation algorithm Download PDFInfo
- Publication number
- CN116008925A CN116008925A CN202211617974.5A CN202211617974A CN116008925A CN 116008925 A CN116008925 A CN 116008925A CN 202211617974 A CN202211617974 A CN 202211617974A CN 116008925 A CN116008925 A CN 116008925A
- Authority
- CN
- China
- Prior art keywords
- rcs
- distance
- waveform
- pitch angle
- gain
- 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.)
- Pending
Links
- 238000013507 mapping Methods 0.000 claims abstract description 10
- 238000010586 diagram Methods 0.000 claims abstract description 9
- 238000012360 testing method Methods 0.000 claims description 20
- 238000000034 method Methods 0.000 claims description 9
- 238000005259 measurement Methods 0.000 claims description 7
- 230000001427 coherent effect Effects 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000009825 accumulation Methods 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 3
- 238000000691 measurement method Methods 0.000 claims 1
- 238000012545 processing Methods 0.000 abstract description 2
- 230000005684 electric field Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses an improved target radar cross-sectional area estimation algorithm, which particularly relates to the field of radar data processing, and specifically comprises the following steps: s1, according to a receiving and transmitting direction diagram of an antenna, obtaining a mapping relation between an azimuth angle, a pitch angle and an antenna gain through fitting; s2, under different waveform modes with known waveform parameters, obtaining RCS coefficients under each waveform mode according to echo parameters of a standard target with known RCS under the condition that azimuth angle and pitch angle are known; s3, for any primary trace, calculating to obtain the RCS of the primary trace according to the RCS coefficient, the distance, the receiving link gain, the echo amplitude, the antenna gain corresponding to the azimuth angle and the pitch angle and other information of the corresponding waveform mode.
Description
Technical Field
The invention relates to the field of radar data processing, in particular to an improved target radar cross-sectional area estimation algorithm.
Background
The key of the object classification recognition problem is to select characteristic statistics which can better distinguish objects. The radar cross-sectional area (Radar Cross Section, RCS) can characterize the scattering properties of the target, and the RCS of different targets are clearly different, so that the RCS is very suitable as an important feature in the field of target recognition.
The formal definition of RCS given in the handbook of radars is:
wherein E is 0 Representing the electric field intensity of an incident wave impinging on the target; e (E) s Representing the electric field intensity of scattered waves at the radar; a typical value of R is the distance of the radar to the target.
The existing similar achievements mainly comprise: patent publication 202111217829.3, a real-time estimation algorithm for radar cross-sectional area, describes an algorithm that first estimates the RCS of the primary trace detected by the radar and then updates the RCS of the trace using a recursive filter. However, the algorithm does not consider the antenna receiving and transmitting pattern gain and also does not consider the fluctuation of the RCS coefficient in different distance segments, so the invention provides an improved target radar cross-sectional area estimation algorithm.
Disclosure of Invention
The invention aims to solve the defects in the prior art and provides an improved target radar cross-sectional area estimation algorithm.
In order to achieve the above purpose, the present invention adopts the following technical scheme:
an improved target radar cross-sectional area estimation algorithm specifically comprises the following steps:
s1, according to a receiving and transmitting direction diagram of an antenna, obtaining a mapping relation between an azimuth angle, a pitch angle and an antenna gain through fitting;
s2, under different waveform modes with known waveform parameters, obtaining RCS coefficients under each waveform mode according to echo parameters of a standard target with known RCS under the condition that azimuth angle and pitch angle are known;
and S3, for any primary trace, calculating the RCS of the primary trace according to the RCS coefficient, the distance, the receiving link gain, the echo amplitude, the antenna gain corresponding to the azimuth angle and the pitch angle and other information corresponding to the waveform mode.
Preferably, the method for calculating the mapping relationship between the azimuth angle, the pitch angle and the antenna gain in the step S1 is as follows: according to the receiving and transmitting direction diagram of the antenna, obtaining the antenna gain F and the corresponding azimuth angle in the beam main lobe by fittingThe mapping relation between the pitch angle theta and the pitch angle theta is as follows
Preferably, the calculation formula of the RCS coefficient in step S2 is as follows:in the above formula, W i RCS coefficient, τ, for the ith waveform pattern i The pulse width of the transmission signal in the ith waveform mode is expressed as seconds s, N i The pulse number accumulated for the phase of the ith waveform mode is dimensionless and d r To test the reference target, the total gain of the receiving link, τ r To test the pulse width of the transmitted signal in waveform mode at the time of reference target, in seconds,N r pulse number accumulated for coherent waveform pattern in test reference target, +.>For testing the reference target, the receiving and transmitting gain of the antenna R r For testing the reference object, the distance of the reference object is expressed in meters, P r For testing the reference target, the echo power of the reference target is expressed as W, sigma r RCS, which is a reference target, is in square meters.
Preferably, the method for measuring the RCS coefficient in step S2 is as follows: and (3) measuring M times in different distance ranges by adopting different targets and different receiving link gains, and taking the average value of RCS coefficient measurement results, namely:if RCS coefficient W i Along with the larger fluctuation of the distance change, the distance can be divided into K distance segments, and the RCS coefficient in each distance segment is ensured to be not changed greatly along with the distance.
Preferably, the RCS calculating method in step S3 is as follows:in sigma 1 RCS of the target corresponding to the one-time trace is expressed as square meter, R 1 The unit of the distance of the one-time trace is meter, P 1 The echo power of the primary point trace is expressed as W, i is the subscript of the waveform mode corresponding to the primary point trace, and d is dimensionless i (R 1 ) At distance R for the ith waveform pattern 1 Receiving link gain at->And theta 1 The azimuth angle and the pitch angle corresponding to the primary track are respectively expressed in degrees +.>For the transceiving gain of the antenna at the corresponding azimuth and elevation angles,dimensionless, W ik RCS coefficients for the ith waveform pattern in the kth distance segment.
The beneficial effects of the invention are as follows:
the invention is suitable for pulse Doppler system radar, the RCS coefficient provided by the invention not only considers parameters such as distance, echo energy and the like in the traditional radar equation, but also considers parameters such as antenna pattern receiving and transmitting gain, pulse width of different waveform modes, coherent accumulation times and the like, in addition, because the RCS coefficient possibly fluctuates along with the distance, the distance is divided into different distance sections, the RCS coefficient fluctuation in each distance section is ensured not to be large, in addition, in the algorithm provided by the invention, a darkroom is not needed for measuring the RCS coefficient, the RCS coefficient of one waveform mode can be obtained through an external field test, and the RCS coefficient of other waveform modes can be obtained. Therefore, the method has the advantages of good practicability, high estimation accuracy and the like.
Drawings
FIG. 1 is an algorithm flow chart;
FIG. 2 is a sequence diagram of RCS coefficients for waveform mode 0;
FIG. 3 is a sequence diagram of RCS coefficients for waveform mode 1;
fig. 4 is an RCS actual measurement result of waveform mode 0 of the target of rcs= -20dBm 2;
fig. 5 is an actual measurement result of RCS in waveform mode 1 of the target of rcs= -20dBm 2.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments.
Examples
As shown in fig. 1-5, an improved target radar cross-sectional area estimation algorithm specifically includes the following steps:
s1, according to a receiving and transmitting direction diagram of an antenna, obtaining a mapping relation between an azimuth angle, a pitch angle and an antenna gain through fitting;
s2, under different waveform modes with known waveform parameters, obtaining RCS coefficients under each waveform mode according to echo parameters of a standard target with known RCS under the condition that azimuth angle and pitch angle are known;
and S3, for any primary trace, calculating the RCS of the primary trace according to the RCS coefficient, the distance, the receiving link gain, the echo amplitude, the antenna gain corresponding to the azimuth angle and the pitch angle and other information corresponding to the waveform mode.
Preferably, the method for calculating the mapping relationship between the azimuth angle, the pitch angle and the antenna gain in the step S1 is as follows: according to the receiving and transmitting direction diagram of the antenna, obtaining the antenna gain F and the corresponding azimuth angle in the beam main lobe by fittingThe mapping relation between the pitch angle theta and the pitch angle theta is as follows
Wherein, the calculation formula of the RCS coefficient in step S2 is as follows:in the above formula, W i RCS coefficient, τ, for the ith waveform pattern i The pulse width of the transmission signal in the ith waveform mode is expressed as seconds s, N i The pulse number accumulated for the phase of the ith waveform mode is dimensionless and d r To test the reference target, the total gain of the receiving link, τ r Pulse width of transmitted signal in seconds, N, of waveform pattern for testing reference target r Pulse number accumulated for coherent waveform pattern in test reference target, +.>For testing the reference target, the receiving and transmitting gain of the antenna R r For testing the reference object, the distance of the reference object is expressed in meters, P r For testing the reference target, the echo power of the reference target is expressed as W, sigma r R as reference targetCS, in square meters.
The method for measuring the RCS coefficient in the step S2 is as follows: and (3) measuring M times in different distance ranges by adopting different targets and different receiving link gains, and taking the average value of RCS coefficient measurement results, namely:if RCS coefficient W i Along with the larger fluctuation of the distance change, the distance can be divided into K distance segments, and the RCS coefficient in each distance segment is ensured to be not changed greatly along with the distance.
The RCS calculating method in the step S3 is as follows:in sigma 1 RCS of the target corresponding to the one-time trace is expressed as square meter, R 1 The unit of the distance of the one-time trace is meter, P 1 The echo power of the primary point trace is expressed as W, i is the subscript of the waveform mode corresponding to the primary point trace, and d is dimensionless i (R 1 ) At distance R for the ith waveform pattern 1 Receiving link gain at->And theta 1 The azimuth angle and the pitch angle corresponding to the primary track are respectively expressed in degrees +.>The receiving and transmitting gain of the antenna at the corresponding azimuth angle and pitch angle is dimensionless, W ik RCS coefficients for the ith waveform pattern in the kth distance segment.
The invention is suitable for pulse Doppler system radar, the RCS coefficient provided by the invention not only considers parameters such as distance, echo energy and the like in the traditional radar equation, but also considers parameters such as antenna pattern receiving and transmitting gain, pulse width of different waveform modes, coherent accumulation times and the like, in addition, because the RCS coefficient possibly fluctuates along with the distance, the distance is divided into different distance sections, the RCS coefficient fluctuation in each distance section is ensured not to be large, in addition, in the algorithm provided by the invention, a darkroom is not needed for measuring the RCS coefficient, the RCS coefficient of one waveform mode can be obtained through an external field test, and the RCS coefficient of other waveform modes can be obtained. Therefore, the method has the advantages of good practicability, high estimation accuracy and the like;
the actual measurement results of the RCS coefficients of two different waveform modes are shown in the following table, and the waveform parameters adopted by the two waveform modes and the RCS of the reference target are shown in the following table. The RCS coefficients for waveform mode 0 and waveform mode 1 are shown in fig. 2 and 3, respectively, and specific experimental data are as follows:
the RCS actual measurement results of waveform pattern 0 and waveform pattern 1 of the target RCS with rcs= -20dBm2 are shown in fig. 4 and 5, respectively.
The foregoing is only a preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any person skilled in the art, who is within the scope of the present invention, should make equivalent substitutions or modifications according to the technical scheme of the present invention and the inventive concept thereof, and should be covered by the scope of the present invention.
Claims (5)
1. An improved target radar cross-sectional area estimation algorithm is characterized by comprising the following steps:
s1, according to a receiving and transmitting direction diagram of an antenna, obtaining a mapping relation between an azimuth angle, a pitch angle and an antenna gain through fitting;
s2, under different waveform modes with known waveform parameters, obtaining RCS coefficients under each waveform mode according to echo parameters of a standard target with known RCS under the condition that azimuth angle and pitch angle are known;
and S3, for any primary trace, calculating the RCS of the primary trace according to the RCS coefficient, the distance, the receiving link gain, the echo amplitude, the antenna gain corresponding to the azimuth angle and the pitch angle and other information corresponding to the waveform mode.
2. The improved target radar cross-sectional area estimation algorithm according to claim 1, wherein the method for calculating the mapping relationship between the azimuth angle, the pitch angle and the antenna gain in step S1 is as follows: according to the receiving and transmitting direction diagram of the antenna, obtaining the antenna gain F and the corresponding azimuth angle in the beam main lobe by fittingThe mapping relation between the pitch angle theta and the pitch angle theta is as follows
3. The improved target radar cross-sectional area estimation algorithm according to claim 1, wherein the calculation formula of the RCS coefficient in step S2 is:in the above formula, W i RCS coefficient, τ, for the ith waveform pattern i The pulse width of the transmission signal in the ith waveform mode is expressed as seconds s, N i The pulse number accumulated for the phase of the ith waveform mode is dimensionless and d r To test the reference target, the total gain of the receiving link, τ r Pulse width of transmitted signal in seconds, N, of waveform pattern for testing reference target r To test the number of pulses of coherent accumulation of waveform patterns at the reference target,for testing the reference target, the receiving and transmitting gain of the antenna R r For testing the reference object, the distance of the reference object is expressed in meters, P r For testing the reference target, the echo power of the reference target is expressed as W, sigma r RCS, which is a reference target, is in square meters.
4. The improved target radar cross-sectional area estimation algorithm according to claim 1, wherein the RCS coefficient measurement method in step S2 is as follows: and (3) measuring M times in different distance ranges by adopting different targets and different receiving link gains, and taking the average value of RCS coefficient measurement results, namely:if RCS coefficient W i Along with the larger fluctuation of the distance change, the distance can be divided into K distance segments, and the RCS coefficient in each distance segment is ensured to be not changed greatly along with the distance.
5. The improved target radar cross-sectional area estimation algorithm according to claim 1, wherein the RCS calculation method in step S3 is as follows:in sigma 1 RCS of the target corresponding to the one-time trace is expressed as square meter, R 1 The unit of the distance of the one-time trace is meter, P 1 The echo power of the primary point trace is expressed as W, i is the subscript of the waveform mode corresponding to the primary point trace, and d is dimensionless i (R 1 ) At distance R for the ith waveform pattern 1 Receiving link gain at->And theta 1 The azimuth angle and the pitch angle corresponding to the primary track are respectively expressed in degrees +.>The receiving and transmitting gain of the antenna at the corresponding azimuth angle and pitch angle is dimensionless, W ik RCS coefficients for the ith waveform pattern in the kth distance segment. />
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211617974.5A CN116008925A (en) | 2022-12-15 | 2022-12-15 | Improved target radar sectional area estimation algorithm |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211617974.5A CN116008925A (en) | 2022-12-15 | 2022-12-15 | Improved target radar sectional area estimation algorithm |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116008925A true CN116008925A (en) | 2023-04-25 |
Family
ID=86027857
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211617974.5A Pending CN116008925A (en) | 2022-12-15 | 2022-12-15 | Improved target radar sectional area estimation algorithm |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116008925A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116299297A (en) * | 2023-05-05 | 2023-06-23 | 南京隼眼电子科技有限公司 | Method and device for calculating radar scattering cross section area and storage medium |
-
2022
- 2022-12-15 CN CN202211617974.5A patent/CN116008925A/en active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116299297A (en) * | 2023-05-05 | 2023-06-23 | 南京隼眼电子科技有限公司 | Method and device for calculating radar scattering cross section area and storage medium |
CN116299297B (en) * | 2023-05-05 | 2024-01-12 | 南京隼眼电子科技有限公司 | Method and device for calculating radar scattering cross section area and storage medium |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107861117B (en) | Multi-target parameter measuring method suitable for continuous wave perimeter surveillance radar | |
JP2010535344A (en) | Ground reflected wave mitigation using parametric time-domain technique | |
KR102013205B1 (en) | Simulation Apparatus and Method for Radar Signal Processing | |
US7336222B2 (en) | System and method for measuring characteristics of a continuous medium and/or localized targets using multiple sensors | |
CN111220955B (en) | Airborne weather radar ground clutter suppression method based on vertical array covariance matrix characteristic decomposition | |
EP2406658B1 (en) | Method of detecting and compensating for precipitation in sodar systems | |
CN112462356A (en) | Method for improving radar ranging precision based on Gaussian interpolation | |
CN116008925A (en) | Improved target radar sectional area estimation algorithm | |
CN116973864A (en) | Method for carrying out RCS measurement on target and realizing target classification under motion condition by ship-borne radar | |
CN106959442B (en) | Method for extracting first-order echo spectrum of ground wave radar in strong interference environment based on multi-domain information | |
CN100585429C (en) | A Passive Channel Correction Method Based on Nonlinear Antenna Array | |
CN113239506A (en) | Tower deformation and wire galloping risk assessment model based on typhoon weather | |
JP3702347B2 (en) | Signal processing method and program and apparatus for wind profiler | |
Vizinho et al. | Evaluation of the use of the modified-covariance method in HF radar ocean measurement | |
Tian et al. | Wave-height map extraction from compact HF surface-wave radar network | |
JP5196959B2 (en) | Radar equipment | |
CN116125406A (en) | Method for evaluating performance of space-based surveillance radar based on SNR estimation and track-pointing track report | |
Jiang et al. | A Novel Airborne Radar Sea Clutter Simulation Method Based on Full Link Factors | |
KR101990078B1 (en) | Simulation Apparatus for Radar Signal Processing | |
Navarro et al. | A calibration-free methodology to estimate significant wave height using time-sequences of x-band marine radar images | |
Tian et al. | Quality control of compact high-frequency radar-retrieved wave data | |
Oyedokun | Sea clutter simulation | |
Huang et al. | Ship detection and tracking using multi-frequency HFSWR | |
JP2000310678A (en) | Wave observation system with radar | |
Chufeng et al. | Signal processing techniques in high-resolution RCS measurement system |
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 |