CN115684718A - Frequency estimation method - Google Patents

Frequency estimation method Download PDF

Info

Publication number
CN115684718A
CN115684718A CN202211201107.3A CN202211201107A CN115684718A CN 115684718 A CN115684718 A CN 115684718A CN 202211201107 A CN202211201107 A CN 202211201107A CN 115684718 A CN115684718 A CN 115684718A
Authority
CN
China
Prior art keywords
frequency
signals
signal
carrying
fft
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202211201107.3A
Other languages
Chinese (zh)
Inventor
王赟
张东良
李鹏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Jiuhang Electronics Co ltd
Original Assignee
Shanghai Jiuhang Electronics Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Jiuhang Electronics Co ltd filed Critical Shanghai Jiuhang Electronics Co ltd
Priority to CN202211201107.3A priority Critical patent/CN115684718A/en
Publication of CN115684718A publication Critical patent/CN115684718A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses a frequency estimation method, which comprises the following steps: 1) A radar echo simulator receives a radio frequency emission signal emitted by a radar; 2) Carrying out down-conversion treatment on the signal in the step 1); 3) Performing analog-to-digital conversion processing on the signals in the step 2); 4) Carrying out digital IQ decomposition processing on the signals in the step 3); 5) Performing FFT on the signals in the step 4); 6) Solving the frequency f1 and f2 of the peak value and the secondary peak value of the FFT according to the signal in the step 5); 7) Generating signals with frequencies from f1 to f2, and respectively carrying out correlation processing on the signals in the step 4); 8) And (3) the frequency of the signal maximum value in the step 7) is the estimated frequency. The digital frequency measurement technology is optimized and improved, the calculation times can be reduced, the calculation time is shortened, the purpose of rapidly measuring the frequency is achieved, and meanwhile, the cost of the system is reduced.

Description

Frequency estimation method
Technical Field
The invention relates to the technical field of digital frequency measurement, in particular to a frequency estimation method.
Background
Frequency estimation of signals is a classical and important signal processing problem, and has extremely wide application in the fields of science and engineering. In an electric power system, frequency is one of the most important power quality parameters, the change of the frequency reflects the dynamic balance between power generation and load, and the dynamic process monitoring of the system can be realized through the accurate estimation of the frequency and the change of the frequency. In addition, in specific applications such as measurement benchmark testing and standard tracing, the frequency estimation accuracy is often required to be higher, the existing frequency estimation method is generally calculated based on hardware, and in order to meet the frequency accuracy, the resource overhead is increased, and the purpose of measuring the frequency at low cost cannot be achieved.
Disclosure of Invention
The present invention is directed to a frequency estimation method to solve the above problems.
The technical problem solved by the invention can be realized by adopting the following technical scheme:
a method of frequency estimation comprising the steps of:
1) A radar echo simulator receives a radio frequency emission signal emitted by a radar;
2) Carrying out down-conversion treatment on the signals in the step 1);
3) Performing analog-to-digital conversion processing on the signals in the step 2);
4) Carrying out digital IQ decomposition processing on the signals in the step 3);
5) Performing FFT on the signal in the step 4);
6) Solving the frequency f1 and f2 of the peak value and the secondary peak value of the FFT of the signal in the step 5);
7) Generating signals with frequencies from f1 to f2, and respectively carrying out correlation processing on the signals in the step 4);
8) And (3) the frequency of the signal maximum value in the step 7) is the estimated frequency.
Due to the adoption of the technical scheme, the digital frequency measurement technology is optimized and improved, so that the calculation times can be reduced, the calculation time can be shortened, the aim of rapidly measuring the frequency can be fulfilled, and the cost of the system can be reduced.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a system block diagram of one embodiment of the present invention.
FIG. 2 is a signal spectrum diagram of an embodiment of the present invention when the sampling frequency fs is 1GHz, the signal is 35MHz, and the number of FFT points is 1024 points.
Fig. 3 is a signal spectrum diagram of fig. 2 after amplification in the vicinity of the peak.
Fig. 4 is a signal spectrum diagram of fig. 2 after amplification in the vicinity of the secondary peak.
Fig. 5 is a schematic diagram of the correlation between the peak and sub-peak values of fig. 3 and 4 and the original signal after correlation processing.
Fig. 6 is a partially enlarged view of fig. 5.
FIG. 7 is a signal spectrum diagram of an embodiment of the present invention when the sampling frequency fs is 1GHz, the signal is 35MHz, and the number of FFT points is 65536 points.
Fig. 8 is a graph of the signal spectrum of fig. 7 after amplification near the peak.
Detailed Description
In order to make the technical means, the creation characteristics, the achievement purposes and the effects of the invention easy to understand, the invention is further explained below.
Referring to fig. 1, a frequency estimation method includes the steps of:
1) A radar echo simulator receives a radio frequency emission signal emitted by a radar;
2) Carrying out down-conversion treatment on the signals in the step 1);
3) Performing analog-to-digital conversion processing on the signals in the step 2);
4) Carrying out digital IQ decomposition processing on the signals in the step 3);
5) Performing FFT on the signals in the step 4);
6) Solving the frequency f1 and f2 of the peak value and the secondary peak value of the FFT of the signal in the step 5);
7) Generating signals with frequencies from f1 to f2, and respectively carrying out relevant processing on the signals in the step 4);
8) And (4) determining the frequency of the maximum value of the signal in the step 7) as the estimated frequency.
For a signal with unknown frequency, an FFT can be performed to find a peak value to determine the frequency, if the frequency of the signal is to be accurately estimated, the number of FFT points is required to be sufficient, because the frequency resolution of the FFT is F = fs/(N-1), fs is a sampling frequency, and N is the number of FFT points, when fs is constant, the larger the N is, the higher the frequency resolution is, the more accurate the frequency estimation is, but the larger the number of FFT points is, the larger the calculation frequency is, the higher the system overhead is, and thus, the frequency measurement in a short time cannot be achieved.
For example, when fs is 1GHz, a signal is 35mhz, and fft points are 1024 points, the spectrum of the signal is shown in fig. 2. The vicinity of the peak is enlarged as shown in fig. 3. It can be seen from fig. 3 that the frequency of the signal peak is 35.68MHz when the number of FFT points is 1024 points, but is different from the real frequency of 35MHz by 0.68 MHz.
The method provided by the invention is implemented as follows, for the above example, the known peak value is 35.68MHz, the frequency of the secondary peak value is 34.7MHz (as shown in fig. 4), signals with the frequency from 34.71MHz to 35.67MHz are generated, the interval is 0.01MHz, the signals are respectively correlated with the original signals, the frequency of the peak value after correlation processing is searched and is the estimated frequency, as shown in fig. 5, and the frequency is amplified and is shown in fig. 6. It can be seen that the estimated frequency for the embodiment of the present invention is 35MHz.
Of course, the number of FFT points may be changed according to actual conditions, for example, when fs is 1GHz, the signal is 35mhz, and the number of FFT points is 65536 points, the spectrum of the signal is as shown in fig. 7, and after the vicinity of the peak is amplified as shown in fig. 8, it can be seen from fig. 8 that the frequency of the peak of the signal is 35.01MHz when the number of FFT points is 65536 points, but the frequency is different from the real frequency of 35MHz by 0.01 MHz. Similarly, the estimated frequency of the embodiment of the present invention can be estimated to be 35MHz by the above method.
The digital frequency measurement technology is optimized and improved, the calculation times can be reduced, the calculation time is shortened, the purpose of rapidly measuring the frequency is achieved, and meanwhile, the cost of the system is reduced.
The foregoing shows and describes the general principles and broad features of the present invention and advantages thereof. It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, which are described in the specification and illustrated only to illustrate the principle of the present invention, but that various changes and modifications may be made therein without departing from the spirit and scope of the present invention, which fall within the scope of the invention as claimed. The scope of the invention is defined by the appended claims and equivalents thereof.

Claims (1)

1. A method of frequency estimation, comprising the steps of:
1) A radar echo simulator receives a radio frequency emission signal emitted by a radar;
2) Carrying out down-conversion treatment on the signals in the step 1);
3) Performing analog-to-digital conversion processing on the signals in the step 2);
4) Carrying out digital IQ decomposition processing on the signals in the step 3);
5) Performing FFT on the signals in the step 4);
6) Solving the frequency f1 and f2 of the peak value and the secondary peak value of the FFT of the signal in the step 5);
7) Generating signals with frequencies from f1 to f2, and respectively carrying out correlation processing on the signals in the step 4);
8) And (3) the frequency of the signal maximum value in the step 7) is the estimated frequency.
CN202211201107.3A 2022-09-29 2022-09-29 Frequency estimation method Pending CN115684718A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211201107.3A CN115684718A (en) 2022-09-29 2022-09-29 Frequency estimation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211201107.3A CN115684718A (en) 2022-09-29 2022-09-29 Frequency estimation method

Publications (1)

Publication Number Publication Date
CN115684718A true CN115684718A (en) 2023-02-03

Family

ID=85064468

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211201107.3A Pending CN115684718A (en) 2022-09-29 2022-09-29 Frequency estimation method

Country Status (1)

Country Link
CN (1) CN115684718A (en)

Similar Documents

Publication Publication Date Title
JP3836328B2 (en) Integrated circuit having AD converter
CN108351405A (en) The circuit and method of the chirp signal linearity and phase noise for determining fmcw radar
CN111610503B (en) Linear frequency modulation signal parameter estimation method based on improved LVD
CN112462356B (en) Method for improving radar ranging accuracy based on Gaussian interpolation
Jinzhi et al. A novel single tone frequency estimation by interpolation using DFT samples with zero-padding
CN115684718A (en) Frequency estimation method
Sato et al. Revisit to accurate ADC testing with incoherent sampling using proper sinusoidal signal and sampling frequencies
CN105282086B (en) The double-tone orthogonal measuring method of frequency change system group delay
CN113406386B (en) Signal frequency accurate estimation method based on digital down-conversion
Attivissimo et al. A study on nonlinear averagings to perform the characterization of power spectral density estimation algorithms
CN112883787B (en) Short sample low-frequency sinusoidal signal parameter estimation method based on spectrum matching
JP5282914B2 (en) Radio identification device
CN112737613B (en) Automatic determination method for frequency range of broadband signal
Artyushenko et al. Correlation Functions of the Noise Modulation Function under the Influence of Stationary and Slow Multiplicative Noise
Nunzi et al. A procedure for highly reproducible measurements of ADC spectral parameters
Belega et al. An effective procedure for the estimation of harmonic parameters of distorted sine-waves
CN115882977A (en) Spectral noise processing method and device, storage medium and electronic equipment
Niranjan et al. FPGA based Identification of Frequency and Phase Modulated Signals by Time Domain Digital Techniques for ELINT Systems.
Keller A New Concept for a Wideband FFT-Based EMI Receiver
CN110940857A (en) Spectrum parameter detection method
Belega et al. Choice of the window used in the interpolated discrete Fourier transform method
Zhivomirov et al. A method for single-tone frequency estimation
CN114740274B (en) Pulse width measuring method, device and system for laser ranging signal
JPWO2002033427A1 (en) Noise measuring device, method, and recording medium
Rahim Evolution and trends of EMI receiver

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