CN108362939A - A kind of frequency domain parameter measurement method of linear FM signal - Google Patents
A kind of frequency domain parameter measurement method of linear FM signal Download PDFInfo
- Publication number
- CN108362939A CN108362939A CN201810098772.1A CN201810098772A CN108362939A CN 108362939 A CN108362939 A CN 108362939A CN 201810098772 A CN201810098772 A CN 201810098772A CN 108362939 A CN108362939 A CN 108362939A
- Authority
- CN
- China
- Prior art keywords
- spectrum
- signal
- frequency
- data
- data segment
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/02—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Resistance Or Impedance (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
The invention discloses a kind of frequency domain parameter measurement methods of linear FM signal, include the following steps：S1. time domain data obtains：Time domain is carried out to signal and draws window, signal is divided into multiple data segments；S2. data convert：Discrete Fourier transform is carried out to each data segment, obtains the corresponding subband spectrum of each data segment；S3. frequency spectrum interpolation：Interpolation processing is carried out to each subband spectrum, obtains the spectrum signal after the corresponding points extension of each data segment；S4. frequency spectrum time delay：Increase time domain delayed data on the spectrum signal after points extend；S5. frequency spectrum integrates：All spectrum signals with time domain delayed data are subjected to frequency spectrum synthesis, restore the complete frequency spectrum of original signal；S6. frequency domain parameter measures：The complete frequency spectrum obtained using recovery carries out frequency domain parameter measurement.The present invention recovers original signal spectrum using short-term spectrum synthesis, has fitted the more spectrum informations of original signal, improves frequency resolution, and then improve bandwidth, the measurement accuracy of initial frequency.
Description
Technical field
The present invention relates to frequency domain parameter measurements, are specifically related to a kind of frequency domain parameter measurement method of linear FM signal.
Background technology
During to spatial domain signal detection, to obtain the frequency domain parameter in intercepted signal pulse descriptive word in time, such as
Bandwidth, initial frequency or centre frequency, Short Time Fourier Transform (STFT) are a kind of common processing methods；By existing to signal
Signal is divided into time shorter data segment, then carries out DFT (discrete Fourier transform) to each short signal and obtain by time domain sliding window
To subband spectrum, its corresponding frequency domain information is analyzed, but for modulated signal such as linear FM signal (LFM), respectively in short-term
The frequency spectrum of data segment can not accurately reflect the frequency domain information of original complete signal, and Fig. 1 is to be carried out to signal in narrow data window
The schematic diagram (short signal frequency spectrum) of spectrum analysis, Fig. 2 is original signal spectrum, and two figures of comparison can be seen that by each section of frequency spectrum
Only reflect the partial frequency spectrum information of original signal, and the measurement of each subband spectrum parameter is influenced by frequency resolution etc., frequency
Rate measurement error is larger, and the parameters such as the bandwidth of original signal, centre frequency are further obtained inherently by each subband measurement parameter
Bring large error.
Invention content
It is an object of the invention to overcome the deficiencies of the prior art and provide a kind of measurements of the frequency domain parameter of linear FM signal
Method recovers original signal spectrum using short-term spectrum synthesis, has fitted the more spectrum informations of original signal, improved frequency
Rate resolution ratio, and then improve bandwidth, the measurement accuracy of initial frequency.
The purpose of the present invention is achieved through the following technical solutions：A kind of frequency domain parameter measurement of linear FM signal
Method includes the following steps：
S1. time domain data obtains：Time domain is carried out to signal and draws window, signal is divided into the data segment x (n) of multiple time domains, n
=0,1 ..., N-1, each data segment include N number of discrete data；
S2. data convert：Discrete Fourier transform is carried out to each data segment, obtains the corresponding subband of each data segment
Frequency spectrum X (k)：
S3. frequency spectrum interpolation：Interpolation processing is carried out to the corresponding subband spectrum X (k) of each data segment, obtains each data segment pair
Spectrum signal X ' (k) after the points extension answered, the spectrum signal X ' (k) is the frequency-region signal for including r*N data, and r is
Positive integer；
S4. frequency spectrum time delay：Using the time domain delay property of Fourier transformation, in the corresponding spectrum signal of each data segment
Increase time domain delayed data on X ' (k)；
S5. frequency spectrum integrates：All spectrum signal X ' (k) with time domain delayed data are subjected to frequency spectrum synthesis, restore former
The complete frequency spectrum of beginning signal；
S6. frequency domain parameter measures：The complete frequency spectrum obtained using recovery carries out frequency domain parameter measurement to original signal.
Wherein, the step S3 includes following sub-step：
S301. it is positive integer to set and mend N number of zero, the r of (r-1) * in x (n), obtains the new sequence x ' (n) of r*N points, then by from
Fourier transformation is dissipated, derives the corresponding spectrum signal X ' (k) of x ' (n)：
It is found that X ' (k) is the linear weighted function of X (k)；
S302. according to derivation result, interpolation is carried out to X (k) in frequency domain, it is described to obtain required spectrum signal X ' (k)
Spectrum signal X ' (k) is the frequency-region signal for including r*N data.
Wherein, the step S4 includes following sub-step：For each data segment, it is delayed using the time domain of Fourier transformation
Property will increase time domain delayed data, obtain on corresponding spectrum signal X ' (k)：
DFT { } indicates that discrete Fourier transform, t indicate that t-th of data segment that time-domain window is chosen, q indicate time-domain window
The time interval of data segment is chosen, (t-1) q is the delayed data of each data segment；In formula,
As increase the spectrum signal after delayed data.
The beneficial effects of the invention are as follows：The present invention recovers original signal spectrum using short-term spectrum synthesis, fits
The more spectrum informations of original signal improve frequency resolution, and then improve bandwidth, the measurement accuracy of initial frequency.
Description of the drawings
Fig. 1 is the schematic diagram of signal spectral analysis in narrow data window；
Fig. 2 is original signal spectrum schematic diagram；
Fig. 3 is flow chart of the method for the present invention；
Fig. 4 is the spectrum recovery schematic diagram in the embodiment of the present invention.
Specific implementation mode
Technical scheme of the present invention is described in further detail below in conjunction with the accompanying drawings, but protection scope of the present invention is not limited to
It is as described below.
As shown in figure 3, a kind of frequency domain parameter measurement method of linear FM signal, includes the following steps：
S1. time domain data obtains：Time domain is carried out to signal and draws window, signal is divided into the data segment x (n) of multiple time domains, n
=0,1 ..., N-1, each data segment include N number of discrete data；
S2. data convert：Discrete Fourier transform is carried out to each data segment, obtains the corresponding subband of each data segment
Frequency spectrum X (k)：
S3. frequency spectrum interpolation：Interpolation processing is carried out to the corresponding subband spectrum X (k) of each data segment, obtains each data segment pair
Spectrum signal X ' (k) after the points extension answered, the spectrum signal X ' (k) is the frequency-region signal for including r*N data, and r is
Positive integer；
S4. frequency spectrum time delay：Using the time domain delay property of Fourier transformation, in the corresponding spectrum signal of each data segment
Increase time domain delayed data on X ' (k)；
S5. frequency spectrum integrates：All spectrum signal X ' (k) with time domain delayed data are subjected to frequency spectrum synthesis, fitting is extensive
The complete frequency spectrum of multiple original signal；
S6. frequency domain parameter measures：The complete frequency spectrum obtained using recovery carries out frequency domain parameter measurement to original signal.
Wherein, the step S3 includes following sub-step：
S301. it is positive integer to set and mend N number of zero, the r of (r-1) * in x (n), obtains the new sequence x ' (n) of r*N points, then by from
Fourier transformation is dissipated, derives the corresponding spectrum signal X ' (k) of x ' (n)：
It is found that in k=m, X ' (k) is the linear weighted function of X (k), can be further simplified as
In k ≠ m, X ' (k) is still the linear weighted function of X (k), and only weights are complex.
S302. according to derivation result, interpolation is carried out to X (k) in frequency domain, it is described to obtain required spectrum signal X ' (k)
Spectrum signal X ' (k) is the frequency-region signal for including r*N data.
The step S4 includes following sub-step：
Each data segment will be increased using the time domain delay property of Fourier transformation on corresponding spectrum signal X ' (k)
Added-time domain delayed data, obtains：
DFT { } indicates that discrete Fourier transform, t indicate that t-th of data segment that time-domain window is chosen, q indicate time-domain window
The time interval of data segment is chosen, (t-1) q is the delayed data of each data segment；In formula,
As increase the spectrum signal after delayed data.
It is that the LFM signal synthesis that a bandwidth is 40MHz restores using the present processes in embodiments herein
The complete signal frequency spectrum gone out, obtains that the results are shown in Figure 4, and comparison diagram 2 is it is found that the frequency spectrum recovered approximate with original signal spectrum one
It causes, it is possible thereby to accurately measure bandwidth, initial frequency parameter.
To sum up, the present invention recovers original signal spectrum using short-term spectrum synthesis, and it is more to have fitted original signal
Spectrum information improves frequency resolution, and then improves bandwidth, the measurement accuracy of initial frequency.
Claims (3)
1. a kind of frequency domain parameter measurement method of linear FM signal, it is characterised in that：Include the following steps：
S1. time domain data obtains：Time domain is carried out to signal and draws window, signal is divided into data segment x (n), the n=0 of multiple time domains,
1 ..., N-1, each data segment include N number of discrete data；
S2. data convert：Discrete Fourier transform is carried out to each data segment, obtains the corresponding subband spectrum of each data segment
X(k)：
S3. frequency spectrum interpolation：Interpolation processing is carried out to the corresponding subband spectrum X (k) of each data segment, it is corresponding to obtain each data segment
Spectrum signal X ' (k) after points extension, the spectrum signal X ' (k) are the frequency-region signal for including r*N data, and r is just whole
Number；
S4. frequency spectrum time delay：Using the time domain delay property of Fourier transformation, in the corresponding spectrum signal X ' of each data segment
(k) increase time domain delayed data on；
S5. frequency spectrum integrates：All spectrum signal X ' (k) with time domain delayed data are subjected to frequency spectrum synthesis, restore original letter
Number complete frequency spectrum；
S6. frequency domain parameter measures：The complete frequency spectrum obtained using recovery carries out frequency domain parameter measurement to original signal.
2. a kind of frequency domain parameter measurement method of linear FM signal according to claim 1, it is characterised in that：The step
Rapid S3 includes following sub-step：
S301. it is positive integer to set and mend N number of zero, the r of (r-1) * in x (n), obtains the new sequence x ' (n) of r*N points, then passes through discrete Fu
In leaf transformation, derive the corresponding spectrum signal X ' (k) of x ' (n)：
It is found that X ' (k) is the linear weighted function of X (k)；
S302. according to derivation result, interpolation is carried out to X (k) in frequency domain, to obtain required spectrum signal X ' (k), the frequency spectrum
Signal X ' (k) is the frequency-region signal for including r*N data.
3. a kind of frequency domain parameter measurement method of linear FM signal according to claim 1, it is characterised in that：The step
Rapid S4 includes following sub-step：
For each data segment, using the time domain delay property of Fourier transformation, when by increasing on corresponding spectrum signal X ' (k)
Domain delayed data, obtains：
DFT { } indicates that discrete Fourier transform, t indicate that t-th of data segment that time-domain window is chosen, q indicate that time-domain window is chosen
The time interval of data segment, (t-1) q are the delayed data of each data segment；In formula,
As increase the spectrum signal after delayed data.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810098772.1A CN108362939B (en) | 2018-01-31 | 2018-01-31 | Frequency domain parameter measuring method of linear frequency modulation signal |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810098772.1A CN108362939B (en) | 2018-01-31 | 2018-01-31 | Frequency domain parameter measuring method of linear frequency modulation signal |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108362939A true CN108362939A (en) | 2018-08-03 |
| CN108362939B CN108362939B (en) | 2020-06-23 |
Family
ID=63007713
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810098772.1A Active CN108362939B (en) | 2018-01-31 | 2018-01-31 | Frequency domain parameter measuring method of linear frequency modulation signal |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108362939B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110007140A (en) * | 2019-04-16 | 2019-07-12 | 北京遥感设备研究所 | A kind of radar equipment LFM signal chirp rate test method and system |
| CN117118536A (en) * | 2023-10-25 | 2023-11-24 | 南京派格测控科技有限公司 | Method, device, equipment and storage medium for determining frequency modulation stability |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2003012416A1 (en) * | 2001-08-01 | 2003-02-13 | Liposcience, Inc. | Method of determining presence and concentration of lipoprotein x in blood plasma and serum |
| CN101566684A (en) * | 2009-06-08 | 2009-10-28 | 武汉大学 | Radar digital signal processing method and device thereof |
| CN101776747A (en) * | 2010-01-15 | 2010-07-14 | 山东大学 | Adaptive resampling-based radar signal compression method |
| CN101833035A (en) * | 2010-04-19 | 2010-09-15 | 天津大学 | Linear frequency-modulated parameter estimating method and implementing device thereof |
| CN102055411A (en) * | 2010-12-22 | 2011-05-11 | 成都凯腾四方数字广播电视设备有限公司 | Power amplifier linearization correcting circuit and method based on multi-channel feedback |
| CN102308546A (en) * | 2011-07-26 | 2012-01-04 | 华为技术有限公司 | Receiving method and device of multiple carrier optical signal |
| CN102680948A (en) * | 2012-05-15 | 2012-09-19 | 东南大学 | Method for estimating modulation frequency and starting frequency of linear frequency-modulated signal |
| CN103226190A (en) * | 2013-03-19 | 2013-07-31 | 南京大学 | True delay-based wideband phased array wave beam formation method |
| CN103532656A (en) * | 2013-08-08 | 2014-01-22 | 北京理工大学 | Broadband linear frequency-modulated (LFM) signal multi-decoy interference method based on fractional Fourier domain channelization |
| CN103869162A (en) * | 2014-03-05 | 2014-06-18 | 湖南大学 | Dynamic signal phasor measurement method based on time domain quasi-synchronization |
-
2018
- 2018-01-31 CN CN201810098772.1A patent/CN108362939B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2003012416A1 (en) * | 2001-08-01 | 2003-02-13 | Liposcience, Inc. | Method of determining presence and concentration of lipoprotein x in blood plasma and serum |
| CN101566684A (en) * | 2009-06-08 | 2009-10-28 | 武汉大学 | Radar digital signal processing method and device thereof |
| CN101776747A (en) * | 2010-01-15 | 2010-07-14 | 山东大学 | Adaptive resampling-based radar signal compression method |
| CN101833035A (en) * | 2010-04-19 | 2010-09-15 | 天津大学 | Linear frequency-modulated parameter estimating method and implementing device thereof |
| CN102055411A (en) * | 2010-12-22 | 2011-05-11 | 成都凯腾四方数字广播电视设备有限公司 | Power amplifier linearization correcting circuit and method based on multi-channel feedback |
| CN102308546A (en) * | 2011-07-26 | 2012-01-04 | 华为技术有限公司 | Receiving method and device of multiple carrier optical signal |
| CN102680948A (en) * | 2012-05-15 | 2012-09-19 | 东南大学 | Method for estimating modulation frequency and starting frequency of linear frequency-modulated signal |
| CN103226190A (en) * | 2013-03-19 | 2013-07-31 | 南京大学 | True delay-based wideband phased array wave beam formation method |
| CN103532656A (en) * | 2013-08-08 | 2014-01-22 | 北京理工大学 | Broadband linear frequency-modulated (LFM) signal multi-decoy interference method based on fractional Fourier domain channelization |
| CN103869162A (en) * | 2014-03-05 | 2014-06-18 | 湖南大学 | Dynamic signal phasor measurement method based on time domain quasi-synchronization |
Non-Patent Citations (3)
| Title |
|---|
| KAIHUI DING: "Parameters estimation of LFM signal based on fractional order cross spectrum", 《PROCEEDINGS OF 2011 INTERNATIONAL CONFERENCE ON ELECTRONIC & MECHANICAL ENGINEERING AND INFORMATION TECHNOLOGY》 * |
| YIN ZHIPING: "LFM signal detection using the origin moment of fractional spectrum", 《2008 9TH INTERNATIONAL CONFERENCE ON SIGNAL PROCESSING》 * |
| 黄慧慧: "连续波背景下多信号处理方法", 《数据采集与处理》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110007140A (en) * | 2019-04-16 | 2019-07-12 | 北京遥感设备研究所 | A kind of radar equipment LFM signal chirp rate test method and system |
| CN117118536A (en) * | 2023-10-25 | 2023-11-24 | 南京派格测控科技有限公司 | Method, device, equipment and storage medium for determining frequency modulation stability |
| CN117118536B (en) * | 2023-10-25 | 2023-12-19 | 南京派格测控科技有限公司 | Method, device, equipment and storage medium for determining frequency modulation stability |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108362939B (en) | 2020-06-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105678273B (en) | The starting point detection algorithm of radio-frequency fingerprint identification technology transient signal | |
| Liu et al. | Deep learning and recognition of radar jamming based on CNN | |
| Wang et al. | A reducing iteration orthogonal matching pursuit algorithm for compressive sensing | |
| CN110133632B (en) | Composite modulation signal identification method based on CWD time-frequency analysis | |
| CN103837740A (en) | High-precision digital instantaneous frequency measurement method and device | |
| KR101294681B1 (en) | Apparatus and method for processing weather signal | |
| KR20080093989A (en) | Detection of presence of television signals embedded in noise using cyclostationary toolbox | |
| CN109343020A (en) | A kind of S-transformation Time-Frequency Analysis Method based on improvement window function | |
| US11449697B2 (en) | Combined higher order statistics and artificial intelligence signal analysis | |
| CN108267657B (en) | Power quality disturbance detection method and system based on S transformation | |
| CN112985639A (en) | Distributed optical fiber sensing method based on chirp removal and time domain sub-chirp pulse extraction | |
| CN108362939A (en) | A kind of frequency domain parameter measurement method of linear FM signal | |
| CN107918146B (en) | A kind of Weak Signal Detection Method based on non-linear extruding S time-frequency conversions | |
| CN110687595B (en) | Seismic data processing method based on time resampling and synchronous extrusion transformation | |
| WO2020078291A1 (en) | Method and device for improving phase measurement accuracy | |
| CN114545342A (en) | Radar pulse signal parameter measuring method using multi-channel reconnaissance receiver | |
| US8175829B2 (en) | Analyzer for signal anomalies | |
| CN117113004A (en) | Linear frequency modulation signal parameter estimation method under low signal-to-noise ratio condition | |
| Razzaq et al. | Instantaneous Frequency Estimation for Frequency-Modulated Signals under Gaussian and Symmetric α-Stable Noise | |
| Sayeed et al. | Optimal quadratic detection and estimation using generalized joint signal representations | |
| CN115345216A (en) | FMCW radar interference elimination method fusing prior information | |
| CN109490857B (en) | Method and system for determining frequency modulation nonlinearity of LFM pulse signal of radar equipment | |
| Boudraa | Instantaneous frequency estimation of fm signals by ψb-energy operator | |
| CN109490625A (en) | A kind of harmonic signal analysis method based on sliding window and Semidefinite Programming | |
| CN110083865B (en) | Anti-pulse noise known signal detection method and detection circuit |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |