CN108363076B - System and method for estimating instantaneous frequency of man-made interference based on short-time Renyi entropy - Google Patents
System and method for estimating instantaneous frequency of man-made interference based on short-time Renyi entropy Download PDFInfo
- Publication number
- CN108363076B CN108363076B CN201810037988.7A CN201810037988A CN108363076B CN 108363076 B CN108363076 B CN 108363076B CN 201810037988 A CN201810037988 A CN 201810037988A CN 108363076 B CN108363076 B CN 108363076B
- Authority
- CN
- China
- Prior art keywords
- frequency
- time
- instantaneous
- interference
- sweep
- 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
Images
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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/21—Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Noise Elimination (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
Abstract
The invention discloses an artificial interference instantaneous frequency estimation system and method based on a short-time Renyi entropy.
Description
Technical Field
The invention belongs to the technical field of navigation receiver interference instantaneous frequency estimation, and particularly relates to a system and a method for estimating artificial interference instantaneous frequency based on short-time Renyi entropy, which are used for estimating instantaneous frequency of artificial interference with sawtooth type abrupt change frequency characteristics.
Background
Global Navigation Satellite Systems (GNSS) provide high precision, all-weather, and global coverage navigation and timing functions that are widely used in more and more civilian applications. In recent years, however, the effective range of the human jammers, which is appearing in large numbers, can reach several tens of meters or even thousands of meters, and is becoming a major problem in the field of civil navigation, for example, causing many damaging events affecting commercial aviation and wireless communication systems.
Most jammers broadcast chirp signals in the GNSS band, obtained by a saw-tooth shaped input voltage controlled oscillator. The instantaneous frequency of the interference signal is equal to the input voltage of the voltage-controlled oscillator, so that the sawtooth-type sudden change frequency characteristic is presented, namely the instantaneous frequency of the interference signal is suddenly changed from the ending frequency to the starting frequency directly near a sudden change point.
Due to the property of overall broadband and instantaneous narrow-band, a time-frequency distribution-based method is usually adopted to process such GNSS receiver jamming signals to obtain jamming frequency sweeping characteristics. And calculating the time-frequency distribution of the received signals by using a time-frequency analysis tool, and obtaining an instantaneous frequency estimation value by detecting time-frequency ridge lines, namely frequency points corresponding to the peak values of all time slices. When the instantaneous frequency of the interference signal changes linearly, the ridge line detection method based on time-frequency distribution can provide higher instantaneous frequency estimation accuracy.
Because the GNSS receiver artificial interference signal has the sudden change frequency characteristic, obvious cross terms exist in the time-frequency distribution of the received signal near the sudden change point of the instantaneous frequency, the time-frequency resolution ratio is deteriorated, and the instantaneous frequency estimated value obtained by the ridge line detection algorithm based on the time-frequency distribution is obviously deviated from the true value.
Disclosure of Invention
In view of this, the invention provides a short-time Renyi entropy-based artificial interference instantaneous frequency estimation system and method, which can perform high-precision estimation on artificial interference instantaneous frequency with sawtooth-shaped frequency mutation characteristics, and an obtained instantaneous frequency estimation value is close to a true value.
The invention provides an artificial interference instantaneous frequency estimation system based on a short-time Renyi entropy, which comprises a time-frequency analysis module, a conventional instantaneous frequency estimation module, an instantaneous component number estimation module, an interference frequency sawtooth function model parameter estimation module and an improved instantaneous frequency estimation module;
the time-frequency analysis module is used for performing time-frequency analysis on the received signals and the reference signals to obtain time-frequency distribution of the received signals and the reference signals, sending the time-frequency distribution of the received signals to the conventional instantaneous frequency estimation module, and sending the time-frequency distribution of the received signals and the reference signals to the instantaneous component number estimation module; the reference signal is a signal with slow frequency spectrum change;
the conventional instantaneous frequency estimation module detects the time-frequency distribution of the received receiving signal to obtain a conventional instantaneous frequency estimation value of the interference signal, and sends the conventional instantaneous frequency estimation value of the interference signal to the interference frequency sawtooth function model parameter estimation module;
the instantaneous component number estimation module obtains short-time Renyi entropies of the received signals and the reference signals under each moment respectively through shifting a time-frequency shielding window according to the time-frequency distribution of the received signals and the reference signals; comparing the short-time Renyi entropies of the received signal and the reference signal at the same time to obtain an instantaneous component number estimated value of the received signal at the moment; comparing the instantaneous component number estimated values of the received signals at each time to obtain the time of the local peak of the instantaneous component number estimated values of the received signals, and arranging the time according to the time sequence to obtain a local peak time sequence P1,P2…PE-1,PEThe local peak time sequence P1,P2…PE-1,PESending the interference frequency sawtooth function model parameter to an interference frequency sawtooth function model parameter estimation module; wherein E is the total number of local peaks;
an interference frequency sawtooth function model parameter estimation module is used for estimating the time sequence P according to the conventional instantaneous frequency estimation value and the local peak value of the received interference signal1,P2…PE-1,PECarrying out model parameter estimation on the interference frequency sawtooth function model, and sending the parameter estimation value of the interference frequency sawtooth function model to an improved instantaneous frequency estimation module; the model parameters of the interference frequency sawtooth function model comprise a frequency sweep period, a frequency sweep starting moment, a frequency sweep starting frequency and a frequency sweep ending frequency;
the interference frequency sawtooth function model parameter estimation value is obtained in the following mode:
the local peak time sequence P1,P2…PE-1,PEThe average value of the time intervals of two adjacent moments is used as the frequency sweep period estimation value
Scanning frequency starting time corresponding to E local peak value timeIs used as the estimated value of the start time of the frequency sweep
Wherein the jth local peak time PjThe corresponding sweep start time is
j＝1,2,3......E；
Using the jth local peak time PjCorresponding interval
Fitting to obtain a ramp-up frequency straight line fup[t]Wherein D is less than or equal to
A constant of (d); according to a ramp-up frequency straight line fup[t]Obtaining
And fup[Pj-1]；
I.e. the jth sweep start frequency
fup[Pj-1]I.e. the jth sweep end frequency
Taking the average value of the sweep frequency starting frequencies corresponding to the E local peak moments as a sweep frequency starting frequency estimated value
Taking the average value of the sweep frequency end frequency corresponding to the E local peak moments as a sweep frequency end frequency estimated value
And the improved instantaneous frequency estimation module substitutes the received parameter estimation value of the interference frequency sawtooth function model into the interference frequency sawtooth function model to obtain an improved instantaneous frequency estimation value at the current moment.
Wherein the interference frequency sawtooth function model is:
wherein the content of the first and second substances,
fosc[n]for the nth modified instantaneous frequency estimate, mod () is a rounding function, where N is 0, …, NTF-1, N is the signal sampling instant, NTFIs the length of the time-frequency observation window, NSWFor the sweep period, fstartFor sweeping the starting frequency, fendIs the sweep end frequency.
The invention also provides a short-time Renyi entropy-based artificial interference instantaneous frequency estimation method, which adopts the short-time Renyi entropy-based artificial interference instantaneous frequency estimation system to estimate and comprises the following steps:
step 1, selecting a signal with slow frequency spectrum change as a reference signal;
according to the time-frequency distribution of the received signals and the reference signals, the short-time Renyi entropies of the received signals and the reference signals at all the moments corresponding to the received signals and the reference signals are obtained by shifting a time-frequency shielding window;
step 2, comparing the short-time Renyi entropies of the received signal and the reference signal at the same time to obtain an instantaneous component number estimated value of the received signal at the moment;
step 3, comparing the instantaneous component number estimated value of the received signal at each time to obtain the time of the local peak value of the instantaneous component number estimated value of the received signal, and arranging according to the time sequence to obtain a local peak value time sequence P1,P2…PE-1,PE；
Step 4, utilizing the local peak value time sequence P1,P2…PE-1,PEObtaining an interference frequency sawtooth function model parameter estimation value by using a conventional instantaneous frequency estimation value;
the interference frequency sawtooth function model parameter estimation value comprises a sweep frequency period estimation value
Starting time of frequency sweep
Estimation value of frequency sweep starting frequency
And end-of-sweep frequency estimate
The conventional instantaneous frequency estimation value is obtained by carrying out time-frequency ridge line detection on the time-frequency distribution of the received signals;
step 5, substituting the parameter estimation value of the interference frequency sawtooth function model into the corresponding parameter of the interference frequency sawtooth function model to obtain an improved instantaneous frequency estimation value fosc[i]；
The interference frequency sawtooth function model is as follows:
wherein the content of the first and second substances,
fosc[n]for the nth modified instantaneous frequency estimate, mod () is a rounding function, where N is 0, …, NTF-1, N is the signal sampling instant, NTFIs the length of the time-frequency observation window, NSWFor the sweep period, fstartFor sweeping the starting frequency, fendIs the sweep end frequency.
Has the advantages that:
(1) the artificial interference instantaneous frequency estimation system based on the short-time Renyi entropy can directly calculate the number of the obtained instantaneous components based on the short-time Renyi entropy, and the instantaneous frequency mutation point information of the interference signal extracted from the instantaneous component number; the improved instantaneous frequency estimated value is calculated by utilizing the interference frequency sawtooth function model parameter estimated value, so that a larger estimation error near an instantaneous frequency mutation point is reduced.
(2) The method combines the conventional instantaneous frequency estimation value obtained by the conventional time-frequency ridge line detection method, calculates the interference frequency sawtooth function model parameter estimation value by using the conventional instantaneous frequency estimation value, calculates the improved instantaneous frequency estimation value by using the interference frequency sawtooth function model parameter estimation value, reduces larger estimation error near the instantaneous frequency catastrophe point, and further realizes high-precision artificial interference instantaneous frequency estimation.
Drawings
FIG. 1 is a flow chart of an instantaneous frequency estimation method according to the present invention.
Detailed Description
The invention is described in detail below by way of example with reference to the accompanying drawings.
The invention provides an artificial interference instantaneous frequency estimation system and method based on a short-time Renyi entropy.
The specific principle is as follows: the time-frequency distribution of the received signals is multiplied by a two-dimensional time-frequency shielding window to separate the time-frequency distribution of the received signals in a short time period and the time-frequency distribution outside the time period is set to be zero. And calculating the short-time Renyi entropy of the received signal at the central moment q of the time-frequency shielding window by utilizing the shielded time-frequency distribution of the received signal. Calculating the short-time Renyi entropy of a reference signal at the central moment q of a time-frequency shielding window by adopting the same process as that of receiving the signal, wherein the reference signal is generally a signal with slow frequency spectrum change; combining the short-time Renyi entropy of the received signal and the reference signal (the Renyi entropy is used to represent the signal complexity), the number n (q) of transient components of the received signal at the time q at the center of the time-frequency mask window can be estimated.
The interference signal in a time-frequency observation window can be regarded as several chirp components with continuous short-time support. At the temporal frequency discontinuity point, the time-frequency distribution within the time-frequency mask window contains one upcoming chirp component and one upcoming chirp component, which results in a sudden increase in the number of temporal components. Therefore, the corresponding relation exists between the moment of the local peak value of the instantaneous component number and the moment of the instantaneous frequency catastrophe point, the instantaneous component number of each moment can be obtained by utilizing the short-time Renyi entropy of each moment of the received signal and the reference signal, the instantaneous component number estimated value of the received signal at each moment is compared, the moment of the local peak value of the instantaneous component number estimated value of the received signal is obtained, and the instantaneous frequency catastrophe point information is further obtained. And calculating an interference frequency sawtooth function model parameter estimation value by combining instantaneous frequency catastrophe point information and an instantaneous frequency estimation value obtained by a conventional time-frequency ridge line detection method, thereby obtaining an improved instantaneous frequency estimation value.
The system comprises a time frequency analysis module, a conventional instantaneous frequency estimation module, an instantaneous component number estimation module, a sawtooth function model parameter estimation module and an improved instantaneous frequency estimation module.
The time-frequency analysis module is used for performing time-frequency analysis on the received signals and the reference signals to obtain time-frequency distribution of the received signals and the reference signals, sending the time-frequency distribution of the received signals to the conventional instantaneous frequency estimation module, and sending the time-frequency distribution of the received signals and the reference signals to the instantaneous component number estimation module; the reference signal is a signal with slow frequency spectrum change; the reference signal is generally chip signal or fixed frequency cosine signal;
the conventional instantaneous frequency estimation module is used for carrying out time-frequency ridge line detection on the time-frequency distribution of the received signals, detecting frequency points corresponding to the peak values of all time slices in the time-frequency distribution of the received signals, obtaining a conventional instantaneous frequency estimation value of the interference signals and sending the conventional instantaneous frequency estimation value of the interference signals to the sawtooth function model parameter estimation module;
the conventional instantaneous frequency estimate of the interference signal is
Where N is the signal sampling time, NTFIs the length of the time-frequency observation window;
the instantaneous component number estimation module obtains short-time Renyi entropies of the received signals and the reference signals under each moment respectively through shifting a time-frequency shielding window according to the time-frequency distribution of the received signals and the reference signals; comparing the short-time Renyi entropies of the received signal and the reference signal at the same time to obtain an instantaneous component number estimated value of the received signal at the moment; comparing the instantaneous component number estimated values of the received signals at each time to obtain the time of the local peak of the instantaneous component number estimated values of the received signals, and arranging the time according to the time sequence to obtain a local peak time sequence P1,P2…PE-1,PEThe local peak time sequence P1,P2…PE-1,PESending the interference frequency sawtooth function model parameter to an interference frequency sawtooth function model parameter estimation module; wherein E is the total number of local peaks;
an interference frequency sawtooth function model parameter estimation module is used for estimating the time sequence P according to the conventional instantaneous frequency estimation value and the local peak value of the received interference signal1,P2…PE-1,PECarrying out model parameter estimation on the interference frequency sawtooth function model, and sending the model parameter estimation value of the interference frequency sawtooth function model to an improved instantaneous frequency estimation module; the model parameters of the interference frequency sawtooth function model comprise a frequency sweep period, a frequency sweep starting moment, a frequency sweep starting frequency and a frequency sweep ending frequency;
the interference frequency sawtooth function model parameter estimation value is obtained in the following mode:
calculating a local peak time sequence P1,P2…PE-1,PEObtaining E-1 time intervals of the time intervals of two adjacent local peak values, averaging to obtain the estimated value of the sweep frequency period
By using
Calculating the starting moment of the frequency sweep: the jth local peak time PjThe corresponding estimated value of the sweep start time is as follows:
wherein j is 1,2,3.. E;
by using
Obtaining sweep frequency starting time corresponding to E local peak value time
And averaging to obtain the estimated value of the start time of the frequency sweep
Selecting the jth sweep period
In the middle part of, i.e. the interval
Wherein D is less than or equal to
A constant of (d);
using the jth local peak time PjCorresponding interval
Fitting to obtain a ramp-up frequency straight line fup[t]Calculating a slope rising frequency straight line fup[t]To obtain a slope rising frequency straight line fup[t]The expression of (1);
according to fup[t]Is expressed to obtain
And fup[Pj-1]；
I.e. the jth sweep start frequency
fup[Pj-1]I.e. the jth sweep end frequency
Taking the average value of the sweep frequency starting frequencies corresponding to the E local peak moments as a sweep frequency starting frequency estimated value
Taking the average value of the sweep frequency end frequency corresponding to the E local peak moments as a sweep frequency end frequency estimated value
The improved instantaneous frequency estimation module substitutes the received interference frequency sawtooth function model parameter estimation value into the corresponding parameter of the interference frequency sawtooth function model to obtain the improved instantaneous frequency estimation value f at the current momentosc[n]；
The interference frequency sawtooth function model is as follows:
wherein the content of the first and second substances,
fosc[n]for the nth modified instantaneous frequency estimate, mod () is a rounding function, where N is 0, …, NTF-1, N is the signal sampling instant, NTFIs the length of the time-frequency observation window, NSWFor the sweep period, fstartFor sweeping the starting frequency, fendIs the sweep end frequency.
The instantaneous frequency estimation method of the man-made interference instantaneous frequency estimation system based on the short-time Renyi entropy comprises the following steps:
step 1, multiplying the time-frequency distribution of the received signals by a time-frequency shielding window, separating the time-frequency distribution of the received signals in a short time period, and setting the time-frequency distribution outside the time period to be zero; calculating a short-time Renyi entropy of a time-frequency shielding window center moment q by utilizing the separated time-frequency distribution of the received signals, and obtaining the short-time Renyi entropy of the received signals at the moment, wherein the Renyi entropy is used for representing the signal complexity;
selecting a cosine signal with fixed frequency as a reference signal, and processing the reference signal by adopting the same process as that of receiving the signal to obtain a short-time Renyi entropy of the reference signal at each moment;
step 2, comparing the short-time Renyi entropies of the received signal and the reference signal at the same time to obtain an instantaneous component number estimated value of the received signal at the moment, and recording the instantaneous component number estimated value as N (q);
step 3, comparing the instantaneous component number estimated value of the received signal at each time to obtain the time of the local peak value of the instantaneous component number estimated value of the received signal, and arranging according to the time sequence to obtain a local peak value time sequence P1,P2…PE-1,PE；
Since the interference signal in a time-frequency observation window can be regarded as several chirp components with continuous short-time support. At the temporal frequency break point, the time-frequency distribution within the time-frequency mask window contains one upcoming chirp component and one upcoming chirp component, which results in a sudden increase in n (q). Therefore, the corresponding relation between the N (q) local peak value and the instantaneous frequency catastrophe point can be utilized to obtain the instantaneous frequency catastrophe point information, namely the local peak value moment of the invention.
Step 4, utilizing all instantaneous frequency catastrophe point information, namely local peak value time sequence P1,P2…PE-1,PECalculating an interference frequency sawtooth function model parameter estimation value according to an instantaneous frequency estimation value obtained by a conventional time-frequency ridge line detection method;
and 5, substituting the parameter estimation value of the interference frequency sawtooth function model into the corresponding parameter of the interference frequency sawtooth function model to obtain an improved instantaneous frequency estimation value, wherein the improved instantaneous frequency estimation value is an instantaneous frequency estimation value with higher precision.
The interference frequency sawtooth function model is as follows:
wherein the content of the first and second substances,
fosc[n]for the nth modified instantaneous frequency estimate, mod () is a rounding function, where N is 0, …, NTF-1, N is the signal sampling instant, NTFIs the length of the time-frequency observation window, NSWFor the sweep period, fstartFor sweeping the starting frequency, fendIs the sweep end frequency.
In summary, the above description is only a preferred embodiment of the present invention, and is not intended to limit the scope of the present invention. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (2)
1. An artificial interference instantaneous frequency estimation system based on a short-time Renyi entropy is characterized by comprising a time-frequency analysis module, a conventional instantaneous frequency estimation module, an instantaneous component number estimation module, an interference frequency sawtooth function model parameter estimation module and an improved instantaneous frequency estimation module;
the time-frequency analysis module is used for performing time-frequency analysis on the received signals and the reference signals to obtain time-frequency distribution of the received signals and the reference signals, sending the time-frequency distribution of the received signals to the conventional instantaneous frequency estimation module, and sending the time-frequency distribution of the received signals and the reference signals to the instantaneous component number estimation module; the reference signal is a signal with slow frequency spectrum change;
the conventional instantaneous frequency estimation module detects the time-frequency distribution of the received receiving signal to obtain a conventional instantaneous frequency estimation value of the interference signal, and sends the conventional instantaneous frequency estimation value of the interference signal to the interference frequency sawtooth function model parameter estimation module;
the instantaneous component number estimation module obtains short-time Renyi entropies of the received signals and the reference signals under each moment respectively through shifting a time-frequency shielding window according to the time-frequency distribution of the received signals and the reference signals; comparing the short-time Renyi entropies of the received signal and the reference signal at the same time to obtain an instantaneous component number estimated value of the received signal at the moment; comparing the instantaneous component number estimated values of the received signals at each time to obtain the time of the local peak of the instantaneous component number estimated values of the received signals, and arranging the time according to the time sequence to obtain a local peak time sequence P1,P2…PE-1,PEThe local peak time sequence P1,P2…PE-1,PESending the interference frequency sawtooth function model parameter to an interference frequency sawtooth function model parameter estimation module; wherein E is the total number of local peaks;
an interference frequency sawtooth function model parameter estimation module is used for estimating the time sequence P according to the conventional instantaneous frequency estimation value and the local peak value of the received interference signal1,P2…PE-1,PECarrying out model parameter estimation on the interference frequency sawtooth function model, and sending the parameter estimation value of the interference frequency sawtooth function model to an improved instantaneous frequency estimation module; the model parameters of the interference frequency sawtooth function model comprise a frequency sweep period, a frequency sweep starting time, a frequency sweep starting frequency and a frequency sweep ending frequency；
The interference frequency sawtooth function model parameter estimation value is obtained in the following mode:
the local peak time sequence P1,P2…PE-1,PEThe average value of the time intervals of two adjacent moments is used as the frequency sweep period estimation value
Taking the average value of the sweep frequency starting time corresponding to the E local peak time as the estimated value of the sweep frequency starting time
Wherein the jth local peak time PjThe corresponding sweep start time is
Using the jth local peak time PjCorresponding interval
Fitting to obtain a ramp-up frequency straight line fup[t]Wherein D is less than or equal to
A constant of (d); according to a ramp-up frequency straight line fup[t]Obtaining
And fup[Pj-1]；
I.e. the jth sweep start frequency
fup[Pj-1]I.e. the jth sweep end frequency
Taking the average value of the sweep frequency starting frequencies corresponding to the E local peak moments as a sweep frequency starting frequency estimated value
Taking the average value of the sweep frequency end frequency corresponding to the E local peak moments as a sweep frequency end frequency estimated value
The improved instantaneous frequency estimation module substitutes the received parameter estimation value of the interference frequency sawtooth function model into the interference frequency sawtooth function model to obtain an improved instantaneous frequency estimation value at the current moment;
the interference frequency sawtooth function model is as follows:
wherein the content of the first and second substances,
fosc[n]for the nth modified instantaneous frequency estimate, mod () is a rounding function, where N is 0, …, NTF-1, N is the signal sampling instant, NTFIs the length of the time-frequency observation window, NSWFor the sweep period, fstartFor sweeping the starting frequency, fendIs the sweep end frequency.
2. An artificial interference instantaneous frequency estimation method based on short-time Renyi entropy is characterized in that the estimation is carried out by adopting the artificial interference instantaneous frequency estimation system based on short-time Renyi entropy as claimed in any claim 1, and the method comprises the following steps:
step 1, selecting a signal with slow frequency spectrum change as a reference signal;
according to the time-frequency distribution of the received signals and the reference signals, the short-time Renyi entropies of the received signals and the reference signals at all the moments corresponding to the received signals and the reference signals are obtained by shifting a time-frequency shielding window;
step 2, comparing the short-time Renyi entropies of the received signal and the reference signal at the same time to obtain an instantaneous component number estimated value of the received signal at the moment;
step 3, comparing the instantaneous component number estimated value of the received signal at each time to obtain the time of the local peak value of the instantaneous component number estimated value of the received signal, and arranging according to the time sequence to obtain a local peak value time sequence P1,P2…PE-1,PE；
Step 4, utilizing the local peak value time sequence P1,P2…PE-1,PEObtaining an interference frequency sawtooth function model parameter estimation value by using a conventional instantaneous frequency estimation value;
the interference frequency sawtooth function model parameter estimation value comprises a sweep frequency period estimation value
Starting time of frequency sweep
Estimation value of frequency sweep starting frequency
And end-of-sweep frequency estimate
The conventional instantaneous frequency estimation value is obtained by carrying out time-frequency ridge line detection on the time-frequency distribution of the received signals;
step 5, a sawtooth function of the interference frequencySubstituting the model parameter estimated value into the corresponding parameter of the interference frequency sawtooth function model to obtain an improved instantaneous frequency estimated value fosc[i]；
The interference frequency sawtooth function model is as follows:
wherein the content of the first and second substances,
fosc[n]for the nth modified instantaneous frequency estimate, mod () is a rounding function, where N is 0, …, NTF-1, N is the signal sampling instant, NTFIs the length of the time-frequency observation window, NSWFor the sweep period, fstartFor sweeping the starting frequency, fendIs the sweep end frequency.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810037988.7A CN108363076B (en) | 2018-01-16 | 2018-01-16 | System and method for estimating instantaneous frequency of man-made interference based on short-time Renyi entropy |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810037988.7A CN108363076B (en) | 2018-01-16 | 2018-01-16 | System and method for estimating instantaneous frequency of man-made interference based on short-time Renyi entropy |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108363076A CN108363076A (en) | 2018-08-03 |
| CN108363076B true CN108363076B (en) | 2021-09-10 |
Family
ID=63006128
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810037988.7A Active CN108363076B (en) | 2018-01-16 | 2018-01-16 | System and method for estimating instantaneous frequency of man-made interference based on short-time Renyi entropy |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108363076B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102053242A (en) * | 2010-10-26 | 2011-05-11 | 北京理工大学 | Signal parameter real-time precision measurement method of synthetic aperture radar |
| CN102809686A (en) * | 2012-07-05 | 2012-12-05 | 北京邮电大学 | Instantaneous frequency detection method and system for frequency sweeping source |
| WO2014203854A1 (en) * | 2013-06-17 | 2014-12-24 | 独立行政法人産業技術総合研究所 | Method and device for measurement of propagation delay characteristic in multipath propagation environment, and external audio perception device |
| CN105320969A (en) * | 2015-11-20 | 2016-02-10 | 北京理工大学 | A heart rate variability feature classification method based on multi-scale Renyi entropy |
| WO2016123484A1 (en) * | 2015-01-29 | 2016-08-04 | Worcester Polytechnic Institute | Motion and noise artifact detection and reconstruction algorithms for photoplethysmogram and equivalent signals |
| CN106054159A (en) * | 2016-05-12 | 2016-10-26 | 北京航空航天大学 | Instantaneous frequency extraction method of Doppler signals |
| WO2017147255A1 (en) * | 2016-02-25 | 2017-08-31 | Trellisware Technologies, Inc. | Estimation and mitigation of swept-tone interferers in frequency-hopped systems |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20040107566A (en) * | 2003-06-05 | 2004-12-23 | 삼성전자주식회사 | estimater for frequency offset of Multi-carrier transmission receiving system and a method thereof |
| CN101496325B (en) * | 2006-07-25 | 2012-07-04 | 富士通株式会社 | Interference noise estimating method in multicarrier communication system, reception processing method, interference noise estimating device and receiver |
| CN103425992B (en) * | 2012-05-16 | 2017-03-08 | 中国科学院声学研究所 | A kind of seafloor sediment classification method and system based on synthetic aperture sonar picture |
-
2018
- 2018-01-16 CN CN201810037988.7A patent/CN108363076B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102053242A (en) * | 2010-10-26 | 2011-05-11 | 北京理工大学 | Signal parameter real-time precision measurement method of synthetic aperture radar |
| CN102809686A (en) * | 2012-07-05 | 2012-12-05 | 北京邮电大学 | Instantaneous frequency detection method and system for frequency sweeping source |
| WO2014203854A1 (en) * | 2013-06-17 | 2014-12-24 | 独立行政法人産業技術総合研究所 | Method and device for measurement of propagation delay characteristic in multipath propagation environment, and external audio perception device |
| WO2016123484A1 (en) * | 2015-01-29 | 2016-08-04 | Worcester Polytechnic Institute | Motion and noise artifact detection and reconstruction algorithms for photoplethysmogram and equivalent signals |
| CN105320969A (en) * | 2015-11-20 | 2016-02-10 | 北京理工大学 | A heart rate variability feature classification method based on multi-scale Renyi entropy |
| WO2017147255A1 (en) * | 2016-02-25 | 2017-08-31 | Trellisware Technologies, Inc. | Estimation and mitigation of swept-tone interferers in frequency-hopped systems |
| CN106054159A (en) * | 2016-05-12 | 2016-10-26 | 北京航空航天大学 | Instantaneous frequency extraction method of Doppler signals |
Non-Patent Citations (1)
| Title |
|---|
| "Estimating the number of components of a multicomponent nonstationary signal using the short-term time-frequency Renyi entropy";Victor Sucic等;《EURASIP journal on advances in signal processing》;20111231;第1-11页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108363076A (en) | 2018-08-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP1611464B1 (en) | System and method for spectral corrected lightning detection | |
| EP3256876B1 (en) | Radio receiver for determining location of a signal source | |
| JP2017003553A (en) | Phase calibration of stepwise chirp signal for synthetic aperture radar | |
| CN111190198B (en) | Satellite receiver of GBAS test equipment and pseudo-range generation method and system thereof | |
| CN109975770B (en) | Separation method and device of time-frequency overlapped multi-component linear frequency modulation signals | |
| KR101135459B1 (en) | Spread spectrum signal receiver, Method for multipath super-resolution thereof, and Recording medium thereof | |
| CN103217669B (en) | Sub-range profile offset deviation-based satellite borne SAR (Synthetic Aperture Radar) ionosphere calibration method | |
| US20200064457A1 (en) | Radar range accuracy improvement method | |
| CN108363076B (en) | System and method for estimating instantaneous frequency of man-made interference based on short-time Renyi entropy | |
| US7336736B2 (en) | Method of detecting and processing pulsed signals in a radio signal | |
| CN110673139B (en) | Target detection radar system pulse compression coefficient generation method and device | |
| JP2008249541A (en) | Pulse signal detector | |
| Sun et al. | Interference Detection Based on Time-Frequency Analysis for GNSS | |
| CN108241144B (en) | FMCW radar waveform modulation method and device | |
| Kang et al. | A study on pulsed-LFM and pulsed-NLFM waveforms for radar systems | |
| Sun et al. | Design of an ultrawideband ionosonde | |
| CN110988832B (en) | Software-defined frequency modulation continuous wave radar system and transmitting signal modulation and echo signal processing method thereof | |
| CN113552597A (en) | Satellite navigation receiver pulse interference monitoring and suppressing method | |
| CN109557537B (en) | On-orbit monitoring and altimetry data correction method for frequency reference drift of radar altimeter | |
| Guangyan et al. | Unambiguous reconstruction and imaging of nonuniform sampling SAR signals | |
| KR100727653B1 (en) | Method for Fast Signal Acquisition in GPS Receiver and Dual Carrier Correlator Therefor | |
| Stein et al. | Suppression of range ambiguities in synthetic aperture radar systems | |
| Ivchenko et al. | Plasma line observations from the EISCAT Svalbard Radar during the International Polar Year | |
| Moreira | The residual error of the reflectivity displacement method of SAR motion extraction | |
| EP2395367A1 (en) | Improvements in or relating to differential carrier phase determination |
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 |