CN102707157A - Single-frequency impulse signal parameter estimating method based on power spectrum - Google Patents
Single-frequency impulse signal parameter estimating method based on power spectrum Download PDFInfo
- Publication number
- CN102707157A CN102707157A CN2012101550368A CN201210155036A CN102707157A CN 102707157 A CN102707157 A CN 102707157A CN 2012101550368 A CN2012101550368 A CN 2012101550368A CN 201210155036 A CN201210155036 A CN 201210155036A CN 102707157 A CN102707157 A CN 102707157A
- Authority
- CN
- China
- Prior art keywords
- signal
- spectrum
- power
- pure
- tau
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The invention discloses a single-frequency impulse signal parameter estimating method based on a power spectrum. The method comprises the following steps: extracting the single-frequency impulse signal power spectrum whose value is to be estimated, mapping a frequency axis of the power spectrum to a time axis to obtain the power time signal; performing the Fourier transform on the power time signal to obtain the magnitude spectrum and the phase spectrum of the power time signal respectively; calculating the impulse width according to the relationship between the impulse width and the magnitude spectrum of the single-frequency impulse signal; calculating a carrier frequency according to the relationship between the carrier frequency and the phase spectrum of the single-frequency impulse signal to estimate the value of the single-frequency impulse signal. According to the single-frequency impulse signal parameter estimating method disclosed by invention, the problem in the prior art that the single-frequency impulse signal parameter cannot be exactly estimated under the environment having low signal to noise ratio is solved; the estimating method is applicable to the cases of spectral peak wrap, oblique change and submersed zero crossing point of a main lobe of the single-frequency impulse signal power spectrum, resulting from noise influences.
Description
Technical field
The present invention relates to signal processing technology, particularly a kind of pure-tone polse signal parameter estimation method based on power spectrum.
Background technology
The valuation of signal parameter always is the focus of communications field research and application, all is extremely important for military electronic countermeasure and the supervision of civilian frequency spectrum.Valuation is meant measured signal analysis, obtains the useful parameter of signal, so that confirm signal purposes and restoring signal.For pure-tone polse signal parameter estimation algorithm, simple parameters valuation scheme is the trend and the main flow of modulation parameter valuation relatively at present, but the noiseproof feature of algorithm is poor.
The pure-tone polse signal is a signal form commonly used in radar and the sonar, is meant that carrier frequency is to immobilize in each pulse of sending, and has advantage convenient, fast, that be easy to realize.Can be expressed as:
Wherein: A is a signal intensity, and τ is a pulse width, f
cBe carrier frequency,
It is the signal initial phase.
Power spectrum signal is expressed as:
P(f)=A
2τ
2sinc
2[τ(f-f
c)]
Pure-tone polse signal parameter commonly used has carrier frequency, pulse width.The signal carrier frequency can characterization signal the radiation frequency range, can reflect function of radar and purposes, to the bearing accuracy that improves the sonar tracking measuring system with to improve tracking effect significant.It generally not in time the place change and change, be a comparatively fixing sorting parameter.Pulse width is the time width between end-of-pulsing time and time of arrival, and it is an important parameter describing pulse characteristics.
Pure-tone polse signal estimation method commonly used at present comprises instantaneous correlation method and zero crossing detection method, but all has the shortcoming of noiseproof feature difference, is not suitable under complex environment, estimating.
Summary of the invention
The object of the invention is: overcome the problem that can not estimate signal parameter under the low signal-to-noise ratio environment exactly that existing pure-tone polse signal estimation technology exists, be applicable to that the pure-tone polse power spectrum signal is owing to the spectrum peak that noise effect causes twists, tiltedly becomes and the covered situation of main lobe zero crossing.
Technical scheme of the present invention is: the invention provides a kind of pure-tone polse signal parameter estimation method based on power spectrum, it is characterized in that this method comprises:
Wherein: with said power spectrum P (f)=A
2τ
2Sinc
2[τ (f-f
c)] frequency axis be mapped to time shaft according to the principle of f=t, obtain said power time signal P (t), A is said pure-tone polse signal intensity, expression formula is as follows:
P(t)=A
2τ
2sinc
2[τ(t-f
c)]
Said power time signal P (t) is carried out Fourier transform obtain F (x), x is the quadrature component after the conversion,
F (x) expression formula is represented as follows:
Because the frequency shift property F (x) of Fourier transform is about longitudinal axis symmetry, the part F ' that gets x >=0 (x) analyzes, and expression formula is as follows:
Obtain the amplitude spectrum of said power time signal | F ' (x) |, expression formula is as follows:
Wherein:, calculate the first zero crossing x of said power time signal amplitude spectrum according to the expression formula of said power time signal amplitude spectrum
0, said first zero crossing x
0Be the pulse width τ=x of said pure-tone polse signal
0
Wherein:, calculate the horizontal ordinate x of first maximum value in the phase spectrum of said power time signal according to the expression formula of said power time signal phase spectrum
1Horizontal ordinate x with last maximum value
N, N is the maximum value number, the horizontal ordinate x of said two maximum value
1And x
NDifference be Δ x=x
N-x
1, the carrier frequency f of said pure-tone polse signal
cExpression formula as follows:
The invention has the beneficial effects as follows: the present invention is through extracting pure-tone polse power spectrum signal to be valuated, and the frequency axis of said power spectrum is mapped to time shaft, obtains the power time signal; Said power time signal is carried out Fourier transform, obtain the amplitude spectrum and the phase spectrum of said power time signal respectively; According to the pulse width of said pure-tone polse signal and the relation between the said amplitude spectrum, calculate said pulse width; According to the carrier frequency of said pure-tone polse signal and the relation between the said phase spectrum, calculate said carrier frequency, realize the valuation of pure-tone polse signal.The invention solves the pure-tone polse signal parameter problem that can not estimate exactly in the prior art under the low signal-to-noise ratio environment, be applicable to that the pure-tone polse power spectrum signal is owing to the spectrum peak that noise effect causes twists, tiltedly becomes and the covered situation of main lobe zero crossing.
Description of drawings
Fig. 1 is the process flow diagram of the pure-tone polse signal parameter valuation based on power spectrum provided by the invention;
Fig. 2 is the synoptic diagram of embodiment of the invention power time signal;
Fig. 3 is the amplitude spectrum of embodiment of the invention power time signal;
Fig. 4 is the phase spectrum of embodiment of the invention power time signal;
Fig. 5 is the pure-tone polse signal pulse width valuation graph of errors of embodiment of the invention signal to noise ratio (S/N ratio)-10dB to 10dB;
Fig. 6 is the pure-tone polse signal carrier frequency valuation graph of errors of embodiment of the invention signal to noise ratio (S/N ratio)-10dB to 10dB.
Embodiment
Below will describe with reference to Fig. 1-6 specific embodiments of the invention.
As shown in Figure 1, the embodiment of the invention carries out comprising the following steps: based on the pure-tone polse signal parameter estimation method of power spectrum
Wherein, definite method of the amplitude spectrum of the time signal of power described in the step 1 and phase spectrum specifically comprises:
With said power spectrum P (f)=A
2τ
2Sinc
2[τ (f-f
c)] frequency axis be mapped to time shaft according to the principle of f=t, obtain said power time signal P (t), A is said pure-tone polse signal intensity, expression formula is as follows:
P(t)=A
2τ
2sinc
2[τ(t-f
c)]
Said power time signal P (t) is carried out Fourier transform obtain F (x), x is the quadrature component after the conversion,
F (x) expression formula is represented as follows:
Because the frequency shift property F (x) of Fourier transform is about longitudinal axis symmetry, the part F ' that gets x >=0 (x) analyzes,, expression formula is as follows:
Obtain the amplitude spectrum of said power time signal | F ' (x) |, expression formula is as follows:
Wherein, the method for the said pure-tone polse signal pulse width of calculating specifically comprises in the step 2:
According to the expression formula of said power time signal amplitude spectrum, calculate the first zero crossing x of said power time signal amplitude spectrum
0, said first zero crossing x
0Be the pulse width τ=x of said pure-tone polse signal
0
Need to prove:
Confirm that the method that concerns between the amplitude spectrum of pulse width and said power time signal of pure-tone polse signal specifically comprises:
The amplitude spectrum expression formula of power time signal is as follows:
According to rectangular function
Can know, when x>=τ | F ' (x) |=0, obtain the first zero crossing x of amplitude spectrum thus
0Pulse width τ=x for the pure-tone polse signal
0
Wherein, the method for the said pure-tone polse signal carrier frequency of calculating specifically comprises in the step 3:
According to the expression formula of said power time signal phase spectrum, calculate the horizontal ordinate x of first maximum value in the phase spectrum of said power time signal
1Horizontal ordinate x with last maximum value
N, N is the maximum value number, the horizontal ordinate x of said two maximum value
1And x
NDifference be Δ x=x
N-x
1, the carrier frequency f of said pure-tone polse signal
cExpression formula as follows:
Need to prove:
Confirm that the method that concerns between the phase spectrum of carrier frequency and said power time signal of pure-tone polse signal specifically comprises:
The phase spectrum expression formula of power time signal is as follows:
The frequency that can draw phase spectrum according to the expression formula of phase spectrum is the carrier frequency of signal, obtains the carrier frequency f of pure-tone polse signal thus
cExpression formula as follows:
The present invention is that embodiment is elaborated to the inventive method with the pure-tone polse signal.
The embodiment of the invention can be applied to military electronic countermeasure and civilian frequency spectrum supervision field; Under the low signal-to-noise ratio environment, estimate signal parameter exactly; Can reflect function of radar and purposes, to the bearing accuracy that improves the sonar tracking measuring system with to improve tracking effect significant.
The following content of the present invention is represented the pure-tone polse signal with abbreviation CW, and SNR representes signal to noise ratio (S/N ratio) with abbreviation.Wherein, carrier frequency f
c=500MHz, the power time signal that pulse width τ=the CW signal of 1 μ s obtains is as shown in Figure 2, and the power time signal is carried out Fourier transform, and the amplitude spectrum that obtains the power time signal respectively is as shown in Figure 3 as shown in Figure 4 with phase spectrum.
The present invention estimates the CW signal parameter under the MATLAB6.5 environment.
CW signal carrier frequency f
c=500MHz, pulse width τ=1 μ s, SNR variation range-10dB-10dB.The present invention to the parameter estimation result of CW signal like table 1, shown in 2; Pulse width evaluated error curve is as shown in Figure 5, and carrier frequency evaluated error curve is as shown in Figure 6.Can know by above result, when SNR at-10dB-10dB, pulse width valuation error is less than 0.019, carrier frequency valuation error is 0, analyzes thus and can know that the inventive method still can be carried out parameter estimation to the CW signal under the low signal-to-noise ratio situation, and precision is higher.
Table 1 signal to noise ratio (S/N ratio) changes parameter estimation result down
Table 2 signal to noise ratio (S/N ratio) changes parameter estimation result down
Claims (1)
1. pure-tone polse signal parameter estimation method based on power spectrum is characterized in that this method comprises:
Step 1, extraction pure-tone polse power spectrum signal P (f) to be valuated; The frequency axis of said power spectrum is mapped to time shaft; Receive the intensity and the time relation of power spectrum; Be called power time signal P (t); Said power time signal is carried out Fourier transform, obtains the amplitude spectrum of said power time signal respectively | F ' (x) | and phase spectrum
Wherein: with said power spectrum P (f)=A
2τ
2Sinc
2[τ (f-f
c)] frequency axis be mapped to time shaft according to the principle of f=t, obtain said power time signal P (t), A is said pure-tone polse signal intensity, expression formula is as follows:
P(t)=A
2τ
2sinc
2[τ(t-f
c)]
Said power time signal P (t) is carried out Fourier transform obtain F (x), x is the quadrature component after the conversion,
F (x) expression formula is represented as follows:
Because the frequency shift property F (x) of Fourier transform is about longitudinal axis symmetry, the part F ' that gets x >=0 (x) analyzes, and expression formula is as follows:
Obtain the amplitude spectrum of said power time signal | F ' (x) |, expression formula is as follows:
Step 2, according to the pulse width τ of said pure-tone polse signal and the amplitude spectrum of said power time signal | F ' (x) | between relation, calculate the pulse width τ of said pure-tone polse signal;
Wherein:, calculate the first zero crossing x of said power time signal amplitude spectrum according to the expression formula of said power time signal amplitude spectrum
0, said first zero crossing x
0Be the pulse width τ=x of said pure-tone polse signal
0
Step 3, according to the carrier frequency f of said pure-tone polse signal
cPhase spectrum with said power time signal
Between relation, calculate the carrier frequency f of said pure-tone polse signal
c
Wherein:, calculate the horizontal ordinate x of first maximum value in the phase spectrum of said power time signal according to the expression formula of said power time signal phase spectrum
1Horizontal ordinate x with last maximum value
N, N is the maximum value number, the horizontal ordinate x of said two maximum value
1And x
NDifference be Δ x=x
N-x
1, the carrier frequency f of said pure-tone polse signal
cExpression formula as follows:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210155036.8A CN102707157B (en) | 2012-05-18 | 2012-05-18 | A kind of pure-tone polse signal parameter estimation method based on power spectrum |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210155036.8A CN102707157B (en) | 2012-05-18 | 2012-05-18 | A kind of pure-tone polse signal parameter estimation method based on power spectrum |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102707157A true CN102707157A (en) | 2012-10-03 |
CN102707157B CN102707157B (en) | 2015-10-28 |
Family
ID=46900090
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210155036.8A Expired - Fee Related CN102707157B (en) | 2012-05-18 | 2012-05-18 | A kind of pure-tone polse signal parameter estimation method based on power spectrum |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102707157B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103063909A (en) * | 2012-12-18 | 2013-04-24 | 天津理工大学 | Linear frequency modulation signal parameter valuation method based on power spectrum |
CN103441975A (en) * | 2013-08-30 | 2013-12-11 | 天津理工大学 | Two-phase coding signal parameter estimation method based on power spectrum |
CN108469602A (en) * | 2018-03-29 | 2018-08-31 | 东南大学 | A kind of pulse signal type automatic distinguishing method based on spectrum signature extraction |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20010037189A1 (en) * | 2000-01-20 | 2001-11-01 | Dan Onu | Method of estimating phase noise spectral density and jitter in a periodic signal |
CN1395400A (en) * | 2002-07-05 | 2003-02-05 | 清华大学 | Automatic identification method of conventional digital signal modulation mode having no need of a priori knowledge |
US20090207897A1 (en) * | 2008-02-14 | 2009-08-20 | Advantest Corporation | Measuring apparatus, test apparatus, recording medium, program and electronic device |
CN101743731A (en) * | 2007-11-05 | 2010-06-16 | 英国广播公司 | Method and apparatus for signal discovery |
-
2012
- 2012-05-18 CN CN201210155036.8A patent/CN102707157B/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20010037189A1 (en) * | 2000-01-20 | 2001-11-01 | Dan Onu | Method of estimating phase noise spectral density and jitter in a periodic signal |
CN1395400A (en) * | 2002-07-05 | 2003-02-05 | 清华大学 | Automatic identification method of conventional digital signal modulation mode having no need of a priori knowledge |
CN101743731A (en) * | 2007-11-05 | 2010-06-16 | 英国广播公司 | Method and apparatus for signal discovery |
US20090207897A1 (en) * | 2008-02-14 | 2009-08-20 | Advantest Corporation | Measuring apparatus, test apparatus, recording medium, program and electronic device |
Non-Patent Citations (3)
Title |
---|
ZHANG J,WANG H P: "A new approach of pitch detection based on morphology filter and wavelet transform", 《PROC. OF THE SECOND INTERNATIONAL WORKSHOP ON EDUCATION TECHNOLOGY AND COMPUTER SCIENCE》 * |
刘柳等: "基于能量分离算法的多分量线性调频信号频率估计", 《天津理工大学学报》 * |
张佳明等: "非参数模型谱估计算法性能比较", 《天津理工大学学报》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103063909A (en) * | 2012-12-18 | 2013-04-24 | 天津理工大学 | Linear frequency modulation signal parameter valuation method based on power spectrum |
CN103063909B (en) * | 2012-12-18 | 2015-10-28 | 天津理工大学 | A kind of linear frequency-modulated parameter estimation method based on power spectrum |
CN103441975A (en) * | 2013-08-30 | 2013-12-11 | 天津理工大学 | Two-phase coding signal parameter estimation method based on power spectrum |
CN103441975B (en) * | 2013-08-30 | 2016-12-28 | 天津理工大学 | A kind of Coded Signals parameter estimation method based on power spectrum |
CN108469602A (en) * | 2018-03-29 | 2018-08-31 | 东南大学 | A kind of pulse signal type automatic distinguishing method based on spectrum signature extraction |
Also Published As
Publication number | Publication date |
---|---|
CN102707157B (en) | 2015-10-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105556332B (en) | The dispersion correction of fmcw radar in pipeline or pipe | |
CN103364783B (en) | Moving target radial velocity non-fuzzy estimation method based on single-channel SAR (synthetic aperture radar) | |
CN107942322B (en) | Multipath radar time-based target distance estimation method | |
US20130127655A1 (en) | Radar device | |
Li et al. | ISAR imaging of nonuniformly rotating target based on the multicomponent CPS model under low SNR environment | |
CN102778674B (en) | Chirp pulse time delay estimation method for non-uniform sampling | |
CN107612587A (en) | A kind of method for parameter estimation for being used for Frequency Hopping Signal in frequency hopping non-cooperative communication | |
CN101702701B (en) | Method for estimating and compensating frequency deviation under very low signal to noise ratio | |
CN103063909B (en) | A kind of linear frequency-modulated parameter estimation method based on power spectrum | |
CN103941089A (en) | Method for estimating sinusoidal signal frequency based on DFT | |
CN103630878A (en) | Frequency band fusion method and device based on adaptive iteration scheme | |
CN106054159A (en) | Instantaneous frequency extraction method of Doppler signals | |
CN103278807A (en) | Time delay estimation method for double-channel under-sampling line sweeping frequency pulse signal | |
CN102707157A (en) | Single-frequency impulse signal parameter estimating method based on power spectrum | |
CN105445767B (en) | BOC signal parameter blind estimating methods based on average Ambiguity Function | |
CN114895248A (en) | Sinusoidal frequency modulation signal parameter estimation method, system and medium | |
CN106330342A (en) | Water sound communication doppler factor estimation method with low computation complexity | |
CN111007473A (en) | High-speed weak target detection method based on distance frequency domain autocorrelation function | |
CN102540153A (en) | Array amplitude and phase error correcting method based on interstation direct wave interference | |
CN106154257A (en) | Precision instrumentation radar secondary frequency measuring method based on FFT Yu apFFT | |
CN103605139B (en) | Carrier frequency and phase estimation method and carrier frequency and phase estimation system applicable to GNSS (global navigation satellite system) receiver | |
Wu et al. | Frequency estimation algorithm for ranging of millimeter wave LFMCW radar | |
CN104914439A (en) | Ultrasonic ranging-based double-phase measuring method | |
CN104320360B (en) | A kind of linear FM signal delay time estimation method based on Fourier Transform of Fractional Order | |
CN103308895A (en) | Overlapped spectrum based low-order phase error estimation 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: 20151028 Termination date: 20200518 |