CN105158737B - A kind of frequency-stepped signal time-frequency detection method - Google Patents
A kind of frequency-stepped signal time-frequency detection method Download PDFInfo
- Publication number
- CN105158737B CN105158737B CN201510587124.9A CN201510587124A CN105158737B CN 105158737 B CN105158737 B CN 105158737B CN 201510587124 A CN201510587124 A CN 201510587124A CN 105158737 B CN105158737 B CN 105158737B
- Authority
- CN
- China
- Prior art keywords
- signal
- frequency
- stepped
- new
- detection method
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
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/021—Auxiliary means for detecting or identifying radar signals or the like, e.g. radar jamming signals
Abstract
The invention provides a kind of frequency-stepped signal time-frequency detection method, comprise the following steps:Auto-correlation processing is carried out to frequency-stepped signal x (t), recurrence interval estimate T between arteries and veins is obtainedr;Radon STFT conversion is carried out to frequency-stepped signal x (t), local linear FM signal is obtainedInitial frequency estimation f0With chirp rate estimate k;Local linear FM signalFrequency-stepped signal x (t) is multiplied, and obtains thick demodulated signal xde(t);To thick demodulated signal xde(t) it is filtered acquisition signal xfilter(t);To signal xfilter(t) reduced, obtained new signal xnew(t);With TrFor window, to new signal xnew(t) the time-frequency detection that Short Time Fourier Transform obtains frequency-stepped signal is carried out;The present invention can detect frequency-stepped signal under low signal-to-noise ratio, and the accuracy detected is higher.
Description
Technical field
The present invention relates to field of signal processing, and in particular to a kind of frequency-stepped signal time-frequency detection method.
Background technology
Frequency-stepped signal is a kind of typical Low probability intercept signal, is had a wide range of applications in field of radar, in letter
It is in uniform stepping Energy distribution in number frequency band range;And linear FM signal is in the straight line of certain slope in the range of signal band
Energy distribution, therefore, linear FM signal can realize the approximate match with frequency-stepped signal to a certain extent;Especially for
Frequency interval is smaller, the longer frequency-stepped signal of pulse train, and from the point of view of time-frequency domain, linear FM signal can be matched more
It is good;It therefore, it can carry out frequency-stepped signal with linear FM signal just demodulation, the energy compression of frequency-stepped signal existed
In one less frequency band;How can be from compared with detecting that it is important that frequency-stepped signal has under low signal-to-noise ratio under non-cooperation
Meaning;Under low signal-to-noise ratio frequency-stepped signal can be detected there is presently no a kind of detection method under non-cooperation.
The content of the invention
The present invention provides frequency-stepped signal time-frequency detection method under a kind of non-condition for cooperation.
A kind of frequency-stepped signal time-frequency detection method, comprises the following steps:
Auto-correlation processing is carried out to frequency-stepped signal x (t), recurrence interval estimate T between arteries and veins is obtainedr;
Radon-STFT conversion is carried out to frequency-stepped signal x (t), local linear FM signal is obtained's
Initial frequency estimation f0With chirp rate estimate k;
Local linear FM signalIt is multiplied with frequency-stepped signal x (t), obtains thick demodulated signal xde(t);
To thick demodulated signal xde(t) it is filtered acquisition signal xfilter(t);
To signal xfilter(t) reduced, obtained new signal xnew(t);
With TrFor window, to new signal xnew(t) the time-frequency inspection that Short Time Fourier Transform obtains frequency-stepped signal is carried out
Survey.
Preferably, the auto-correlation processing is according to the interval between the main peak and secondary peak of auto-correlation function, obtain between arteries and veins
Recurrence interval estimate Tr。
Preferably, the peak that the Radon-STFT conversion is converted according to Radon-STFT obtains original frequency and estimated
Evaluation f0With chirp rate estimate k.
Preferably, the filter bandwidht that the filtering is used is no more than kTr。
Preferably, the computational methods of the reduction are as follows:
The beneficial effects of the invention are as follows:
The present invention can detect frequency-stepped signal under low signal-to-noise ratio, and the accuracy detected is higher.
Brief description of the drawings
Fig. 1 is flow chart of the present invention.
Embodiment
The present invention will be further described with specific embodiment below in conjunction with the accompanying drawings.
A kind of frequency-stepped signal time-frequency detection method, comprises the following steps:
Auto-correlation processing is carried out to frequency-stepped signal x (t), recurrence interval estimate T between arteries and veins is obtainedr;
Radon-STFT conversion is carried out to frequency-stepped signal x (t), local linear FM signal is obtained's
Initial frequency estimation f0With chirp rate estimate k;
Local linear FM signalIt is multiplied with frequency-stepped signal x (t), obtains thick demodulated signal xde(t);
To thick demodulated signal xde(t) it is filtered acquisition signal xfilter(t);
To signal xfilter(t) reduced, obtained new signal xnew(t);
With TrFor window, to new signal xnew(t) the time-frequency inspection that Short Time Fourier Transform obtains frequency-stepped signal is carried out
Survey.
The time-domain expression of frequency-stepped signal is:
In formula:TrFor recurrence interval estimate between arteries and veins, TpFor pulse width, fAFor carrier frequency initial frequency estimate, Δ f is
The stepping-in amount of carrier frequency, rect (t) is rectangular function, and N is u in the total length of frequency stepped pulse trains sequence, formulai(t) it is pulse bag
Network function;
Auto-correlation processing is carried out to frequency-stepped signal:
χ (τ)=∫ x (t) x* (t- τ) dt
By searching for the interval between the main peak of auto-correlation function and secondary peak, pulse recurrence interval estimate T is obtainedr;
It can be handled by other auto-correlation processing methods.
Radon-STFT conversion is done to frequency-stepped signal, if the STFT time-frequency representations for receiving signal are STFTx(t, ω),
Then the Radon-STFT of signal is transformed to:
By searching for the peak of Radon-STFT conversion, local linear FM signal is obtainedIt is initial
Frequency estimation f0With chirp rate estimate k.
Linear FM signalIt is multiplied with frequency-stepped signal x (t), obtains and the stepping frequency modulation received is believed
Number thick demodulation xde(t),
Bandwidth is carried out to thick demodulated signal and is no more than kTrFiltering, filter out-of-band noise, to lift signal to noise ratio, filtered
Signal x after ripplefilter(t)。
Filtered signal is reduced, new signal x is obtainednew(t):
With TrSTFT (Short Time Fourier Transform) conversion is done to the signal after reduction for window length, time-frequency detection STFT is obtainedXnew
(t,ω)。
Presently preferred embodiments of the present invention is the foregoing is only, is not intended to limit the invention, all essences in the present invention
Any modification, equivalent and improvement made within refreshing and principle etc., should be included within the scope of the present invention.
Claims (5)
1. a kind of frequency-stepped signal time-frequency detection method, it is characterised in that comprise the following steps:
Auto-correlation processing is carried out to frequency-stepped signal x (t), recurrence interval estimate T between arteries and veins is obtainedr;
Radon-STFT conversion is carried out to frequency-stepped signal x (t), local linear FM signal is obtainedIt is initial
Frequency estimation f0With chirp rate estimate k;
Local linear FM signalIt is multiplied with frequency-stepped signal x (t), obtains thick demodulated signal xde(t);
To thick demodulated signal xde(t) it is filtered acquisition signal xfilter(t);
To signal xfilter(t) reduced, obtained new signal xnew(t);
With TrFor window, to new signal xnew(t) the time-frequency detection that Short Time Fourier Transform obtains frequency-stepped signal is carried out.
2. a kind of frequency-stepped signal time-frequency detection method according to claim 1, it is characterised in that at the auto-correlation
The interval between the main peak and secondary peak according to auto-correlation function is managed, recurrence interval estimate T between arteries and veins is obtainedr。
3. a kind of frequency-stepped signal time-frequency detection method according to claim 1, it is characterised in that the Radon-
STFT conversion converts peak according to Radon-STFT and obtains local linear FM signalOriginal frequency estimate
Evaluation f0With chirp rate estimate k.
4. a kind of frequency-stepped signal time-frequency detection method according to claim 1, it is characterised in that the filtering is used
Filter bandwidht be no more than kTr。
5. a kind of frequency-stepped signal time-frequency detection method according to claim 1, it is characterised in that the meter of the reduction
Calculation method is as follows:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510587124.9A CN105158737B (en) | 2015-09-14 | 2015-09-14 | A kind of frequency-stepped signal time-frequency detection method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510587124.9A CN105158737B (en) | 2015-09-14 | 2015-09-14 | A kind of frequency-stepped signal time-frequency detection method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105158737A CN105158737A (en) | 2015-12-16 |
CN105158737B true CN105158737B (en) | 2017-07-11 |
Family
ID=54799655
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510587124.9A Active CN105158737B (en) | 2015-09-14 | 2015-09-14 | A kind of frequency-stepped signal time-frequency detection method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105158737B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105911349B (en) * | 2016-05-31 | 2019-01-11 | 清华大学 | Based on the Linear chirp basic parameter evaluation method and device for resetting time-frequency spectrum |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101984363B (en) * | 2010-10-27 | 2013-03-06 | 南京航空航天大学 | Ultrahigh-resolution synthetic aperture radar (SAR) imaging method based on frequency-stepped system |
CN102866391B (en) * | 2012-09-05 | 2014-08-27 | 中北大学 | Multi-target detection method based on short-time Fourier transform and fractional Fourier transform |
-
2015
- 2015-09-14 CN CN201510587124.9A patent/CN105158737B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN105158737A (en) | 2015-12-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102680948B (en) | Method for estimating modulation frequency and starting frequency of linear frequency-modulated signal | |
CN101944926B (en) | Compressed sampling based estimating method of arrival time of pulse ultra-wide band signal | |
CN102546499B (en) | Fractional-order channelized receiving method of real linear frequency modulation (LFM) signal | |
CN104268883A (en) | Time-frequency spectrum curve extracting method based on edge detection | |
CN109412644B (en) | Doppler frequency estimation method for direct sequence spread spectrum MSK signal | |
CN105743612A (en) | Method for real-time blind demodulation of large-frequency-shift short-time burst signal | |
CN106054159A (en) | Instantaneous frequency extraction method of Doppler signals | |
CN107171994B (en) | Radio fuze signal identification and reconstruction system and method | |
CN104215954B (en) | Networking method for stepping variable-frequency radar communication integration | |
CN105429719A (en) | Strong interference signal detection method based on power spectrum and multiple dimensioned wavelet transformation analysis | |
JP5752481B2 (en) | Passive radar device, guidance device, and radio wave detection method | |
CN102866429B (en) | Method for determining groundwater occurrence | |
CN103063909A (en) | Linear frequency modulation signal parameter valuation method based on power spectrum | |
CN105158737B (en) | A kind of frequency-stepped signal time-frequency detection method | |
CN100349007C (en) | Method for eliminating single frequency interference in signal of seismological record | |
CN102307055A (en) | DSSS (direct sequence spread spectrum) frequency domain interference detection method | |
CN105445767A (en) | BOC signal parameter blind estimation method based on average ambiguity function | |
CN106330342B (en) | A kind of underwater sound communication Doppler factor estimation method of low computation complexity | |
CN106772271A (en) | A kind of radar signal TOA methods of estimation based on matching auto-correlation algorithm | |
CN103944845A (en) | Related detection method for full-frequency-spectrum carrier modulation | |
CN102185811B (en) | Carrier frequency estimation method | |
CN106501787A (en) | Coded Signals method for parameter estimation based on smooth pseudo derivative feedback | |
CN104320360B (en) | A kind of linear FM signal delay time estimation method based on Fourier Transform of Fractional Order | |
CN105812119A (en) | High-precision ultra wide band signal propagation time delay measurement method | |
CN104181508A (en) | Threatening radar signal detection method based on compressed sensing |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20190528 Address after: 610000 9-storey 928, No. 56 Tianhui Middle Street, Chengdu High-tech Zone, Sichuan Province Patentee after: Chengdu bang bang technology Co., Ltd. Address before: 610225 24 section 1 Xuefu Road, Southwest Airport Economic Development Zone, Chengdu, Sichuan Patentee before: Chengdu information engineering university |