CN109917360A - A kind of irregular PRI estimation method of aliasing pulse - Google Patents
A kind of irregular PRI estimation method of aliasing pulse Download PDFInfo
- Publication number
- CN109917360A CN109917360A CN201910154517.9A CN201910154517A CN109917360A CN 109917360 A CN109917360 A CN 109917360A CN 201910154517 A CN201910154517 A CN 201910154517A CN 109917360 A CN109917360 A CN 109917360A
- Authority
- CN
- China
- Prior art keywords
- pri
- aliasing
- pulse
- irregular
- phase
- 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
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a kind of diversity impulse recurrence interval estimation methods of aliasing pulse.The specific steps of which are as follows: determining the numberical range of radar PRI, the corresponding discrete Fourier transform of each PRI is calculated；The remainder that aliasing pulse arrival time is counted in each PRI is distributed, and is made scalar product with corresponding DFT sequence and obtained aliasing pulse energy value under the PRI；Standardization is made to remainder distribution, energy value respectively by the remainder theory quantity of each PRI, obtains the aliasing phase-resolved partial discharge and standard energy value under each PRI；Energy threshold and phase threshold are set, judge aliasing pulse PRI numerical value that may be present and corresponding phase, makees after harmonic management the irregular PRI in the available aliasing pulse of comprehensive analysis again.The present invention can accurately estimate the irregular PRI hidden in aliasing pulse, and have completely anti-false ability and good anti-missing ability.
Description
Technical field
The present invention relates to radar signal fields, in particular to the irregular PRI estimation of aliasing pulse under complex electromagnetic environment.
Background technique
As information technology is in the fast development of military field, electromagnetic space at except sea, land and air day four-dimension battlefield space it
Outer quintuple space, electromagnetism power processed become the basis of information superiority.Complex electromagnetic environment on Information Battlefield causes to fight
Acquisition of information is unsmooth, and primary resolving ideas is the radar signal separated under complex electromagnetic environment, needs to obtain aliasing state thus
Under each portion's radar PRI (Pulse Repeat Interval, write a Chinese character in simplified form PRI herein).It is how accurate under complex electromagnetic environment
Identify various types, the PRI of different radars is current urgent problem.
Kofler proposes statistic histogram method first: measuring the reaching time-difference between adjacent pulse, then one by one with pulse
Numerical value is spaced as abscissa, the pulse spacing frequency occurs as ordinate, draws statistic histogram, analysis histogram determines
PRI.The algorithm can it is few in radiation source quantity, shake it is small in the case where accurately estimate PRI, but will appear in result PRI times
Number, i.e. harmonic wave.In order to inhibit the appearance of subharmonic, avoid subharmonic interference PRI identification as a result, Mardia propose it is accumulative
Poor histogram (CDIF) method: drawing multistage pulses interval histogram, in addition subharmonic detects, comprehensive analysis multistage histogram
As a result PRI is obtained.CDIF algorithm is suitable for the case where more multi radiation sources, and has better jitter immunity, but multistage histogram
Figure means the increase of algorithm execution time.In order to reduce operand,Propose sequence difference histogram (SDIF) calculation
Method: histograms at different levels are not accumulated, and obtain PRI by the way that suitable threshold check function is arranged.Do not influencing the anti-missing energy of anti-jitter
On the basis of power, SDIF algorithm ratio CDIF algorithm speed is faster.For the purposes of inhibiting harmonic wave, Nishiguchi proposes PRI change
It changes method: using estimate of autocorrelation PRI, and passing through the influence of phase factor counteracting subharmonic.PRI converter technique has completely
Inhibit the ability of harmonic wave, and there is good anti-missing ability, but vulnerable to effect of jitter.Nishiguchi analyzes PRI
The reason of transformation is vulnerable to effect of jitter, using the overlapping PRI case with shift time source, improves PRI converter technique, makes PRI
Converter technique is provided with good anti-jitter ability.Although the above classics PRI deinterleaving algorithm can be estimated in multiple radiation sources
PRI, but the signal of complex electromagnetic environment is intensive, broad categories degree far surpasses its upper limit.
In compressed sensing field, Vaidyanathan proposes nested period dictionary structure, and phase estimate problem is stated
It is solved at Optimal solution problem is sought.Experiment shows that nested period dictionary structure is sub all to estimating to hide from complex data
Phase has good result, wherein best based on Ramanujan summation and the Ramanujan summation dictionary effect of building.Xu Chengwei, which has studied, to be drawn
In Ma Nujin dictionary application to minor cycle radar pulse signal, and show the dictionary with good anti-jitter, it is anti-missing,
Anti- falseness ability.Although Ramanujan summation dictionary estimates that hiding subcycle works well, operand is limited its application.For
Ramanujan summation dictionary is designed to Ramanujan summation filter group by simplified calculating, Tenneti, is substituted dilute in compressed sensing
Reconstruct is dredged, greatly shortens Riming time of algorithm, and ensure that algorithm effect is constant.Nevertheless, compressed sensing cycle estimator is answered
It is still infeasible to use the estimation field radar PRI: first is that the Ramanujan summation dictionary data amount under large period is excessive, second is that algorithm is wanted
The aliasing signal data volume asked is much larger than its pulse arrival time (TOA) data volume.Moreover, the class of compressed sensing cycle estimator
Type does not include stagger cycle.
It is difficult to apply to the estimation field radar PRI to solve compressed sensing algorithm, and straight not under complex electromagnetic environment
The problem of estimating irregular PRI is connect, this paper presents period profile is obtained using compressed sensing principle, is obtained by statistics remainder
Phase distribution, in conjunction with the new method of period profile and phase distribution analysis estimation PRI numerical value, PRI type.Imitative legitimate reading
Show under complex electromagnetic environment, this method can accurately estimate conventional, shake, irregular PRI, and have good anti-tremble
Dynamic, anti-missing, anti-false ability, have the value for being applied to and estimating the field aliasing signal PRI under complex electromagnetic environment.
Summary of the invention
The purpose of the present invention is to propose to a kind of irregular PRI estimation methods of aliasing pulse.
Specific step is as follows:
Step 1: determining the numberical range of radar PRI, the corresponding DFT sequence of each PRI is calculated；
Step 2: remainder of the statistics aliasing pulse TOA in each PRI is distributed, and make scalar product with corresponding DFT sequence
Obtain the aliasing pulse energy value under the PRI；
Step 3: making standardization to remainder distribution, energy value respectively by the remainder theory quantity of each PRI, obtain
Aliasing phase-resolved partial discharge and standard energy value under each PRI；
Step 4: setting energy threshold and phase threshold, judge aliasing pulse PRI numerical value that may be present and corresponding phase
The irregular PRI in the available aliasing pulse of comprehensive analysis again is made after harmonic management in position.
DFT sequence in the step 1, the specific steps are as follows:
If radar PRI numberical range is [Pmin,Pmax], some PRI numerical value is q ∈ [Pmin,Pmax], and q ∈ Z, then numerical value q
Corresponding DFT sequence are as follows:
DqData representation be q tie up row vector.
Its remainder distribution and energy value at each PRI is obtained using aliasing pulse TOA in the step 2, it is specific to walk
It is rapid as follows:
If the TOA of aliasing pulse is [t1,t2,...,tn] ∈ Z, some PRI numerical value is q, EqFor the unit square of q × q dimension
Battle array, then remainder statistical result of the aliasing pulse at q are as follows:
Wherein, Eq(j :) indicate EqJth row vector, be remainder symbol, NqData representation be q tie up row vector.
To NqAnd DqAliasing pulse energy value at q is obtained as vector product:
Standardization is made to remainder distribution and energy value in the step 3, obtains the phase distribution and mark of each PRI
Quasi-energy value, the specific steps are as follows:
If some PRI numerical value is q, the TOA range of aliasing pulse is Δ t=tn-t1, then the remainder theory quantity under q be
Δ t/q, respectively to NqAnd pqMake standardization, obtains phase distribution R of the aliasing pulse at qqWith standard energy value
Rq=Nq×q/Δt
RqData representation be q tie up row vector.
Energy threshold and phase threshold are first set in the step 4, hidden by phase distribution and standard energy value
PRI and corresponding phase value, then comprehensive analysis obtain the irregular PRI in aliasing pulse, the specific steps are as follows:
Setting energy threshold is τp, judgment criteria energy valueWith τpSize relation, define energy threshold discriminant function:
Setting phase threshold is τr, judge phase distribution RqWith τrSize relation, define phase threshold discriminant function:
Wherein, Rq(i) phase distribution R is indicatedqThe i-th element,Data representation be indefinite length row to
Amount, if its occurrence is
PRI meets if it existsD | q, and m > 1 then show that the stagger cycle of the PRI is [r2-r1,...,rm-
rm-1,r1-rm+q]；D | q indicate d be q aliquot, i.e. q d=0；
If PRI meetsAnd m=1, then show that the PRI is general type；
If PRI meetsAndThen show that the PRI is wobble type.
Beneficial effect
The invention proposes a kind of irregular PRI of aliasing pulse to estimate that new method, advantage are: Method And Principle is simply easy to
It realizes；In the case where without any known parameter, can direct estimation go out in pulse hide unknown irregular PRI；In pulse aliasing
Under serious situation, which can also estimate the irregular PRI hidden in aliasing pulse；In addition, the algorithm can also be estimated
Conventional PRI in aliasing pulse and shake PRI out.
Detailed description of the invention
Fig. 1 is flow diagram of the invention；
Fig. 2 is the irregular PRI discrimination of single pulse with false rate variation tendency；
Fig. 3 is the irregular PRI discrimination of single pulse with miss rate variation tendency；
Fig. 4 is the irregular PRI discrimination of aliasing pulse with false rate variation tendency；
Fig. 5 is the irregular PRI discrimination of aliasing pulse with false rate variation tendency；
Fig. 6 is the energy profile of aliasing pulse；
Fig. 7 is the maximum phase distribution map of aliasing pulse；
Fig. 8 is phase distribution figure of the aliasing pulse when PRI is respectively 2232,3397,3700,4100
Specific embodiment
Below by specific embodiment, the present invention will be described in detail.
(1) numberical range for determining radar PRI calculates the corresponding DFT sequence of each PRI (Discrete Fourier
Transform writes a Chinese character in simplified form DFT herein)；
If radar PRI numberical range is [Pmin,Pmax], some PRI numerical value is q ∈ [Pmin,Pmax], and q ∈ Z, then numerical value q
Corresponding DFT sequence are as follows:
DqData representation be q tie up row vector.
(2) statistics aliasing pulse TOA (Time Of Arrival, write a Chinese character in simplified form TOA herein) is distributed in the remainder of each PRI, and
Make scalar product with corresponding DFT sequence and obtains aliasing pulse energy value under the PRI；
If the TOA of aliasing pulse is [t1,t2,...,tn] ∈ Z, some PRI numerical value is q, EqFor the unit square of q × q dimension
Battle array, then remainder statistical result of the aliasing pulse at q are as follows:
Wherein, Eq(j :) indicate EqJth row vector, be remainder symbol, NqData representation be q tie up row vector.
To NqAnd DqAliasing pulse energy value at q is obtained as vector product:
(3) standardization is made to remainder distribution, energy value respectively by the remainder theory quantity of each PRI, obtains each
Aliasing phase-resolved partial discharge and standard energy value under PRI；
If some PRI numerical value is q, the TOA range of aliasing pulse is Δ t=tn-t1, then the remainder theory quantity under q be
Δ t/q, respectively to NqAnd pqMake standardization, obtains phase distribution R of the aliasing pulse at qqWith standard energy value
Rq=Nq×q/Δt
RqData representation be q tie up row vector.
(4) energy threshold and phase threshold are set, judge aliasing pulse PRI numerical value that may be present and corresponding phase, is made
Irregular PRI in the available aliasing pulse of comprehensive analysis again after harmonic management.
Setting energy threshold is τp, judgment criteria energy valueWith τpSize relation, define energy threshold discriminant function:
Setting phase threshold is τr, judge phase distribution RqWith τrSize relation, define phase threshold discriminant function:
Wherein, Rq(i) phase distribution R is indicatedqThe i-th element,Data representation be indefinite length row to
Amount, if its occurrence is
PRI meets if it existsD | q, and m > 1 then show that the stagger cycle of the PRI is [r2-r1,...,rm-
rm-1,r1-rm+q].D | q indicate d be q aliquot, i.e. q d=0.
If PRI meetsAnd m=1, then show that the PRI is general type；
If PRI meetsAndThen show that the PRI is wobble type.
Example 1: the irregular PRI estimation performance of single pulse is analyzed
Analyze performance experimental verification in terms of two that context of methods estimates single irregular PRI: recognition success rate is with irregular arteries and veins
Rush the variation tendency of false rate；Recognition success rate with diversity impulse miss rate variation tendency.
Parameter setting in experiment: minimax period Pmin=60, Pmax=5000；Maximum non-uniform quantity is 5；Aliasing pulse
TAO range is Δ t=30 × Pmax, it is meant that every minimum pulse number of radar signals is 30；τ is arranged in threshold parameterp=0.6,
τr=0.75.
The setting of staggered PRI radar pulse signal: irregular PRI quantity is 2~5；Irregular PRI numerical value summation belongs to Pmin~Pmax；Arteries and veins
Rushing TOA range is Δ t；Diversity impulse falseness rate is 0%~90% (stepping numerical value is 10%), and miss rate is 0%~80% (step
It is 1%) into numerical value.
Experiment is tested 9000 times altogether, has been respectively obtained context of methods and has been identified as under different false rates and different miss rates
Power, then the influence of false rate and miss rate to context of methods is drawn respectively, as a result such as Fig. 2~Fig. 3.Fig. 2 is observed it is found that ginseng
Poor pulse falseness rate does not influence this paper algorithm, this, which also demonstrates this paper algorithm from some angle, can be used for estimating aliasing arteries and veins
Irregular PRI in punching；Observe Fig. 3 it is found that context of methods under single pulse have good anti-missing ability, although with
Diversity impulse miss rate increases, and the accuracy rate of the estimation irregular PRI of single pulse can be gradually reduced, even if being in missing 50%
When, accuracy rate is still up to 70%.Therefore, context of methods is with good performance in terms of estimating single irregular PRI.
Example 2: the irregular PRI of analysis aliasing pulse estimates performance
Analyze the performance also experimental verification in terms of two of the context of methods estimation irregular PRI of aliasing: recognition success rate is with aliasing
The variation tendency of pulse falseness rate；Recognition success rate with aliasing pulse miss rate variation tendency.
Parameter setting in experiment: minimax period Pmin=60, Pmax=5000；Maximum non-uniform quantity is 5；Aliasing pulse
TAO range is Δ t=30 × Pmax, it is meant that every minimum pulse number of radar signals is 30；τ is arranged in threshold parameterp=0.6,
τr=0.75.
The setting of aliasing radar pulse signal: two staggered PRI radar pulse signal aliasings, every radar diversity PRI quantity are 2
~5；Every radar diversity PRI numerical value summation belongs to Pmin~Pmax；Aliasing pulse TOA range is Δ t；Aliasing pulse falseness rate
For 0%~90% (stepping numerical value is 10%), miss rate is 0%~80% (stepping numerical value is 1%).
Experiment is tested 9000 times altogether, has been respectively obtained context of methods and has been identified as under different false rates and different miss rates
Power, then the influence of false rate and miss rate to this paper algorithm is drawn respectively, as a result such as Fig. 4~Fig. 5.Fig. 4 is observed it is found that mixed
Folded pulse falseness rate does not influence context of methods still；Observation Fig. 5 is it is found that context of methods still has under aliasing pulse
Good anti-missing ability, even accuracy rate still has 64.6% when lacking 50%, though the standard of relatively more single irregular PRI estimation
True rate has decline, but undeniably context of methods works well in terms of irregular PRI in estimation aliasing pulse.
Example 3: verifying context of methods can estimate different type PRI
Actual conditions are copied, the simulation parameter that aliasing radar signal is arranged is as shown in the table:
1 simulation parameter of table
Parameter setting in experiment: minimax period Pmin=60, Pmax=5000；Aliasing pulse TAO range is Δ t=
30×Pmax, it is meant that every minimum pulse number of radar signals is 30；τ is arranged in threshold parameterp=0.6, τr=0.75.
It is handled by context of methods, the Energy distribution of available above-mentioned aliasing pulse, maximum phase distribution, phase distribution
(Fig. 6~Fig. 8), it includes that PRI is as shown in table 2, and the simulation parameter in contrast table 1 shows the side this paper that analysis, which can obtain aliasing pulse,
Method can actually estimate the irregular PRI in aliasing pulse, and have good anti-miss rate.In addition, context of methods can also be estimated
Routine PRI and shake PRI are counted, and there is good anti-jitter ability.
2 PRI estimated result of table
| Initial phase | Modulation system | PRI |
| 827,1703 | It is irregular | 876us,1356us |
| 384,1327,2107,2981 | It is irregular | 943us,780us,874us,800us |
| 2257 | It is conventional | 3700us |
| 1601 | It is conventional | 4100us |
| Shake | 4774.5us | |
| Shake | 1775.5us |
Claims (5)
1. a kind of irregular PRI estimation method of aliasing pulse, it is characterised in that: the following steps are included:
Step 1: determining the numberical range of radar PRI, the corresponding DFT sequence of each PRI is calculated；
Step 2: remainder of the statistics aliasing pulse TOA in each PRI is distributed, and makees scalar product with corresponding DFT sequence and must be somebody's turn to do
Aliasing pulse energy value under PRI；
Step 3: making standardization to remainder distribution, energy value respectively by the remainder theory quantity of each PRI, obtain each
Aliasing phase-resolved partial discharge and standard energy value under PRI；
Step 4: setting energy threshold and phase threshold, judge aliasing pulse PRI numerical value that may be present and corresponding phase, make
Irregular PRI in the available aliasing pulse of comprehensive analysis again after harmonic management.
2. a kind of irregular PRI estimation method of aliasing pulse according to right 1, it is characterised in that: in the step 1
DFT sequence, the specific steps are as follows:
If radar PRI numberical range is [Pmin,Pmax], some PRI numerical value is q ∈ [Pmin,Pmax], and q ∈ Z, then numerical value q is corresponding
DFT sequence are as follows:
DqData representation be q tie up row vector.
3. a kind of irregular PRI estimation method of aliasing pulse according to right 1, it is characterised in that: sharp in the step 2
Its remainder distribution and energy value at each PRI is obtained with aliasing pulse TOA, the specific steps are as follows:
If the TOA of aliasing pulse is [t1,t2,...,tn] ∈ Z, some PRI numerical value is q, EqFor the unit matrix of q × q dimension, then mix
Folded remainder statistical result of the pulse at q are as follows:
Wherein, Eq(j :) indicate EqJth row vector, be remainder symbol, NqData representation be q tie up row vector.
To NqAnd DqAliasing pulse energy value at q is obtained as vector product:
4. a kind of irregular PRI estimation method of aliasing pulse according to right 1, it is characterised in that: right in the step 3
Remainder distribution and energy value make standardization, obtain the phase distribution and standard energy value of each PRI, the specific steps are as follows:
If some PRI numerical value is q, the TOA range of aliasing pulse is Δ t=tn-t1, then the remainder theory quantity under q is Δ t/q,
Respectively to NqAnd pqMake standardization, obtains phase distribution R of the aliasing pulse at qqWith standard energy value
Rq=Nq×q/Δt
RqData representation be q tie up row vector.
5. a kind of irregular PRI estimation method of aliasing pulse according to right 1, it is characterised in that: in the step 4 first
Energy threshold and phase threshold are set, obtained hiding PRI and corresponding phase value, then comprehensive point by phase distribution and standard energy value
Analysis obtains the irregular PRI in aliasing pulse, the specific steps are as follows:
Setting energy threshold is τp, judgment criteria energy valueWith τpSize relation, define energy threshold discriminant function:
Setting phase threshold is τr, judge phase distribution RqWith τrSize relation, define phase threshold discriminant function:
Wherein, Rq(i) phase distribution R is indicatedqThe i-th element,Data representation be indefinite length row vector, if
Its occurrence is
PRI meets if it existsD | q, and m > 1 then show that the stagger cycle of the PRI is [r2-r1,...,rm-rm-1,
r1-rm+q]；D | q indicate d be q aliquot, i.e. q d=0；
If PRI meetsAnd m=1, then show that the PRI is general type；
If PRI meetsAndThen show that the PRI is wobble type.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910154517.9A CN109917360A (en) | 2019-03-01 | 2019-03-01 | A kind of irregular PRI estimation method of aliasing pulse |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910154517.9A CN109917360A (en) | 2019-03-01 | 2019-03-01 | A kind of irregular PRI estimation method of aliasing pulse |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109917360A true CN109917360A (en) | 2019-06-21 |
Family
ID=66962892
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910154517.9A Pending CN109917360A (en) | 2019-03-01 | 2019-03-01 | A kind of irregular PRI estimation method of aliasing pulse |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109917360A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110334471A (en) * | 2019-07-17 | 2019-10-15 | 吉林工程技术师范学院 | The PRI subcycle extracting method of non-uniform series and arteries and veins group non-uniform series |
| CN110988836A (en) * | 2019-12-06 | 2020-04-10 | 航天恒星科技有限公司 | Method and system for measuring pulse arrival time |
| CN111257839A (en) * | 2020-03-30 | 2020-06-09 | 吉林大学 | Radar signal sorting method |
| CN116821658A (en) * | 2023-06-29 | 2023-09-29 | 中国船舶集团有限公司第七二三研究所 | Clock period fingerprint feature extraction method suitable for different repetition interval types |
Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4721958A (en) * | 1985-10-23 | 1988-01-26 | Trw Inc. | Real-time pulse processor |
| US5583505A (en) * | 1995-09-11 | 1996-12-10 | Lockheed Martin Corporation | Radar pulse detection and classification system |
| CN1437746A (en) * | 2000-02-29 | 2003-08-20 | 高通股份有限公司 | Method and apparatus for tracking the phase of a quasi-periodic signal |
| WO2011158233A1 (en) * | 2010-06-15 | 2011-12-22 | Elbit Systems Ew And Sigint - Elisra Ltd. | A method and a system for determining location properties of an emitter |
| CN102298141A (en) * | 2011-05-16 | 2011-12-28 | 中国人民解放军海军航空工程学院 | Airborne pulse doppler radar iterative solution range ambiguity method |
| JP2012083264A (en) * | 2010-10-13 | 2012-04-26 | Mitsubishi Electric Corp | Positioning method |
| CN103814273A (en) * | 2011-06-14 | 2014-05-21 | 南方创新国际股份有限公司 | Method and apparatus for identifying pulses in detector output data |
| CN104105449A (en) * | 2011-12-01 | 2014-10-15 | 毛伊图像公司 | Motion detection using ping-based and multiple aperture Doppler ultrasound |
| CN104360334A (en) * | 2014-11-21 | 2015-02-18 | 吉林大学 | MIMO radar target blind detection method based on characteristic values under correlated noise background |
| CN104897962A (en) * | 2015-06-19 | 2015-09-09 | 天津大学 | Single-frequency signal short sample high precision frequency measurement method and device based on relatively prime perception |
| CN105277923A (en) * | 2015-11-20 | 2016-01-27 | 清华大学 | Single channel radar signal sorting method |
| CN105891793A (en) * | 2016-03-30 | 2016-08-24 | 中国人民解放军国防科学技术大学 | Radar seeker space-time correlation coherent K distribution clutter modeling method |
| CN107085213A (en) * | 2017-05-19 | 2017-08-22 | 中国人民解放军63892部队 | The moving target ISAR imaging methods designed based on random Based on Modulated Step Frequency Waveform |
| US9791563B1 (en) * | 2014-01-08 | 2017-10-17 | National Technology & Engineering Solutions Of Sandia, Llc | Joint synthetic aperture radar plus ground moving target indicator from single-channel radar using compressive sensing |
| CN109031215A (en) * | 2018-06-27 | 2018-12-18 | 扬州大学 | Staggered PRI radar pulse repetition period estimation method under high impulse loss |
-
2019
- 2019-03-01 CN CN201910154517.9A patent/CN109917360A/en active Pending
Patent Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4721958A (en) * | 1985-10-23 | 1988-01-26 | Trw Inc. | Real-time pulse processor |
| US5583505A (en) * | 1995-09-11 | 1996-12-10 | Lockheed Martin Corporation | Radar pulse detection and classification system |
| CN1437746A (en) * | 2000-02-29 | 2003-08-20 | 高通股份有限公司 | Method and apparatus for tracking the phase of a quasi-periodic signal |
| WO2011158233A1 (en) * | 2010-06-15 | 2011-12-22 | Elbit Systems Ew And Sigint - Elisra Ltd. | A method and a system for determining location properties of an emitter |
| JP2012083264A (en) * | 2010-10-13 | 2012-04-26 | Mitsubishi Electric Corp | Positioning method |
| CN102298141A (en) * | 2011-05-16 | 2011-12-28 | 中国人民解放军海军航空工程学院 | Airborne pulse doppler radar iterative solution range ambiguity method |
| CN103814273A (en) * | 2011-06-14 | 2014-05-21 | 南方创新国际股份有限公司 | Method and apparatus for identifying pulses in detector output data |
| CN104105449A (en) * | 2011-12-01 | 2014-10-15 | 毛伊图像公司 | Motion detection using ping-based and multiple aperture Doppler ultrasound |
| US9791563B1 (en) * | 2014-01-08 | 2017-10-17 | National Technology & Engineering Solutions Of Sandia, Llc | Joint synthetic aperture radar plus ground moving target indicator from single-channel radar using compressive sensing |
| CN104360334A (en) * | 2014-11-21 | 2015-02-18 | 吉林大学 | MIMO radar target blind detection method based on characteristic values under correlated noise background |
| CN104897962A (en) * | 2015-06-19 | 2015-09-09 | 天津大学 | Single-frequency signal short sample high precision frequency measurement method and device based on relatively prime perception |
| CN105277923A (en) * | 2015-11-20 | 2016-01-27 | 清华大学 | Single channel radar signal sorting method |
| CN105891793A (en) * | 2016-03-30 | 2016-08-24 | 中国人民解放军国防科学技术大学 | Radar seeker space-time correlation coherent K distribution clutter modeling method |
| CN107085213A (en) * | 2017-05-19 | 2017-08-22 | 中国人民解放军63892部队 | The moving target ISAR imaging methods designed based on random Based on Modulated Step Frequency Waveform |
| CN109031215A (en) * | 2018-06-27 | 2018-12-18 | 扬州大学 | Staggered PRI radar pulse repetition period estimation method under high impulse loss |
Non-Patent Citations (6)
| Title |
|---|
| MEHMET BURAK KOCAMIS等: "Deinterleaving for radar warning receivers with missed pulse consideration", 《2016 EUROPEAN RADAR CONFERENCE (EURAD)》 * |
| 万鹏武: "TDOA被动定位关键技术研究与应用", 《中国博士学位论文全文数据库社会科学Ⅰ辑》 * |
| 张保群: "一种抗SDIF分选的脉冲重复间隔参差设计方法", 《兵器装备工程学报》 * |
| 苏焕程等: "一种基于余数周期的PRI精确估计算法", 《航天电子对抗》 * |
| 陈潜等: "基于随机脉冲重复间隔Radon-Fourier变换的相参积累", 《电子与信息学报》 * |
| 隋金坪等: "基于随机PRI压缩感知雷达的速度假目标识别方法", 《电子学报》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110334471A (en) * | 2019-07-17 | 2019-10-15 | 吉林工程技术师范学院 | The PRI subcycle extracting method of non-uniform series and arteries and veins group non-uniform series |
| CN110334471B (en) * | 2019-07-17 | 2023-06-27 | 长春电子科技学院 | PRI sub-period extraction method for ragged sequence and pulse group ragged sequence |
| CN110988836A (en) * | 2019-12-06 | 2020-04-10 | 航天恒星科技有限公司 | Method and system for measuring pulse arrival time |
| CN110988836B (en) * | 2019-12-06 | 2021-12-28 | 航天恒星科技有限公司 | Method and system for measuring pulse arrival time |
| CN111257839A (en) * | 2020-03-30 | 2020-06-09 | 吉林大学 | Radar signal sorting method |
| CN111257839B (en) * | 2020-03-30 | 2021-12-21 | 吉林大学 | Radar signal sorting method |
| CN116821658A (en) * | 2023-06-29 | 2023-09-29 | 中国船舶集团有限公司第七二三研究所 | Clock period fingerprint feature extraction method suitable for different repetition interval types |
| CN116821658B (en) * | 2023-06-29 | 2024-04-12 | 中国船舶集团有限公司第七二三研究所 | Clock period fingerprint feature extraction method suitable for different repetition interval types |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109917360A (en) | A kind of irregular PRI estimation method of aliasing pulse | |
| CN105137498B (en) | A kind of the underground objects detection identifying system and method for feature based fusion | |
| Guan et al. | Adaptive fractional Fourier transform-based detection algorithm for moving target in heavy sea clutter | |
| Antipov | Simulation of sea clutter returns | |
| CN106872958B (en) | Radar target self-adapting detecting method based on linear fusion | |
| CN109521404B (en) | FMCW radar-based vibration measurement accuracy evaluation method and system | |
| CN109270497A (en) | The multi-Dimensional parameters Pre-sorting method of radar pulse signal | |
| CN103439692A (en) | STAP method based on wide symmetrical characteristic of covariance matrix | |
| CN103399297A (en) | Machine learning based ultra-wideband NLOS (non line of sight) identification method | |
| CN108802722A (en) | It is a kind of based on tracking before the Faint target detection virtually composed | |
| CN102230961A (en) | Method for detecting spread target based on phase compensation processing | |
| CN109613527A (en) | A kind of the detection threshold generation method and device of moving target | |
| CN104714225A (en) | Dynamic programming tracking-before-detection method based on generalized likelihood ratios | |
| CN109633598A (en) | Array radar object detection method based on signature analysis | |
| CN109932695A (en) | A kind of method and device improving object recognition speed | |
| CN105022044A (en) | Denoising processing-based actual measurement sea clutter modeling method | |
| CN113109784A (en) | Radar pulse repetition interval estimation method based on blind source separation | |
| CN103777189A (en) | Radar weak target detecting method based on information geometry multiple autoregressive model | |
| CN107390194A (en) | A kind of radar target detection method based on the graceful divergence of full Donald Bragg | |
| CN111401168A (en) | Multi-layer radar feature extraction and selection method for unmanned aerial vehicle | |
| CN104881521A (en) | Label random-set filter distribution type fusion method | |
| CN111562570A (en) | Vehicle sensing method for automatic driving based on millimeter wave radar | |
| CN110161478A (en) | A kind of waveform design method based on clutter power spectrum density self-optimizing | |
| CN107861108B (en) | A kind of method for detecting targets at sea based on power spectrum box counting dimension | |
| CN107479037A (en) | A kind of PD radar clutters area judging method |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190621 |