CN112014811A - Method for finely estimating radar carrier frequency - Google Patents
Method for finely estimating radar carrier frequency Download PDFInfo
- Publication number
- CN112014811A CN112014811A CN202010829967.6A CN202010829967A CN112014811A CN 112014811 A CN112014811 A CN 112014811A CN 202010829967 A CN202010829967 A CN 202010829967A CN 112014811 A CN112014811 A CN 112014811A
- Authority
- CN
- China
- Prior art keywords
- fft
- carrier frequency
- formula
- estimation
- frequency
- 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
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000001228 spectrum Methods 0.000 claims abstract description 44
- 238000005070 sampling Methods 0.000 claims description 6
- 235000007164 Oryza sativa Nutrition 0.000 claims description 4
- 235000009566 rice Nutrition 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 239000006185 dispersion Substances 0.000 claims description 3
- 240000007594 Oryza sativa Species 0.000 claims 1
- 241000209094 Oryza Species 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Mathematical Physics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a fine estimation method of radar carrier frequency, which belongs to the field of fine estimation of carrier frequency, and comprises the following steps of firstly, extracting a maximum value in a radar carrier signal FFT output amplitude spectrum sequence and a maximum value in a value adjacent to the maximum value; then, obtaining an estimation threshold value according to the statistical characteristics of the FFT output sequence, and judging whether the two values exceed the estimation threshold value so as to judge whether the radar carrier wave is successfully captured; finally, according to the constructed discrimination function linearly related to the carrier frequency, the compensation quantity of the FFT process for carrier frequency estimation is calculated, so that the estimation precision of the radar carrier frequency is improved; the invention has simple structure and small calculation amount, can greatly improve the estimation precision of the carrier frequency by adding the carrier frequency fine estimation module after the carrier frequency estimation structure in the original FFT process, and has higher engineering use value.
Description
Technical Field
The invention belongs to the field of fine estimation of radar carrier frequency, and particularly relates to a fine estimation method of radar carrier frequency.
Background
The carrier frequency of radar is an important parameter in radar frequency domain parameters, and the existing radar carrier frequency estimation generally adopts an FFT frequency estimator, but the frequency estimation precision of FFT is related to the data length.
Increasing the frequency estimation accuracy by increasing the length of the FFT processing sequence can lead to the counting of the system
Computational resources are consumed exponentially.
Disclosure of Invention
In order to solve the technical problems of the background art, the invention aims to improve the estimation precision of the radar carrier frequency according to the FFT output amplitude spectrum characteristic, and provides a fine estimation method of the radar carrier frequency on the basis of not increasing the calculated amount and not changing the original structure of the FFT estimated carrier frequency.
In order to achieve the purpose, the invention provides the technical scheme that: radar carrier frequency essence
A fine estimation method comprising the steps of:
step one, extracting a maximum value M in an FFT output amplitude spectrum sequence S (k) of a radar carrier signalFFTAnd the maximum value S of the values adjacent theretoubm;
Step two, respectively judging the maximum value MFFTAnd the most adjacent value to itLarge value SubmAnd an estimated threshold value YtIf both are larger than the estimated threshold value YtIf yes, continuing to estimate, otherwise, ending;
step three, obtaining the compensation value of the FFT process to the carrier frequency estimation according to the formula (1),
wherein, Deltaf' is the compensation value of the FFT process to the carrier frequency estimation, Deltaf is the carrier frequency after frequency reduction, TLIs the size of the rectangular window;
and step four, obtaining a radar carrier frequency fine estimation value according to the compensation value.
Further, the first step includes:
the structure of the radar carrier signal s (t) after frequency reduction is shown in formula (2),
where A is the amplitude of the carrier signal, i is the imaginary unit, Δ f is the down-converted carrier frequency,is the carrier initial phase;
after the radar carrier signal s (t) subjected to frequency reduction is dispersed and cut, as shown in formula (3),
where (·) is a Dirac Delta function, TsIs the sampling period, rect () is a rectangular window function, TLIs the size of the rectangular window;
the formula (3) is abbreviated as:
wherein d (t) is a sampling pulse signal; w (t) is a rectangular window signal;
the spectrum of S (n) is obtained according to the time domain multiplication frequency domain convolution principle as shown in the following formula,
wherein, is a convolution symbol;is the frequency spectrum of a discrete finite long time domain signal s (n); (f) is the frequency spectrum of the continuous infinite time domain signal S (t); d (f) is the frequency spectrum of the continuous infinite time domain signal d (t); w (f) is the frequency spectrum of a continuous infinite time domain signal w (t);
each term in the formula (5) is developed, the amplitude spectrum sequence of the DTFT output is shown as the following formula,
where sinc (·) is a sinc function;
according to FFT principle, continuous infinite length frequency spectrum obtained by DTFTPerforming dispersion and truncation to obtain an amplitude spectrum sequence output by the FFT:
as can be seen from equation (7), the maximum value of the FFT output magnitude spectrum sequence can be represented by equation (8),
in the FFT output amplitude spectrum sequence, one value among the points adjacent to the maximum value point of the FFT output is only smaller than the maximum value, as shown in equation (9),
where round (. cndot.) is an integer function, MFFTRepresenting the maximum value, S, of the FFT output amplitude spectrum sequenceubmRepresents and MFFTThe larger of the neighboring points.
The estimated threshold value Y in the second steptThe method comprises the following steps:
noise mixed in the FFT amplitude spectrum belongs to Rayleigh (Rayleigh) distribution, and FFT amplitude output belongs to Rice (Rice) distribution, so that a calculation formula (11) for estimating a threshold value can be deduced according to a formula (10);
in the formula (f)n(. is a probability density function of the Rayleigh distribution, PfaIs the false alarm rate, Yt is the estimated threshold, and σ is the noise statistical standard deviation.
In the third step, a compensation value of the FFT process for carrier frequency estimation is obtained according to formula (1), further comprising:
a linear discrimination function D related to the down-converted carrier frequency Deltaf is constructed according to the formula (9)isc(. cndot.) represented by formula (12).
In the formula, | · | is an absolute value symbol.
If SubmAnd MFFTAdjacent positions in the FFT output amplitude sequence, the FFT estimation compensation formula can be obtained by formula (1),
where Δ f' is a compensation value for the carrier frequency estimation in the FFT process.
Optionally, the step four further includes:
if the maximum value S among the adjacent valuesubmAt a maximum M at the position of the FFT output amplitude sequenceFFTTo the left of the position, the estimated carrier frequency estimateCan be expressed by formula (13); if on the right, it can be expressed by equation (14),
in the formula (f)FFTFor the FFT process to estimate the carrier frequency,a fine estimate of the compensated carrier frequency.
Compared with the prior art, the invention has the beneficial effects that: the fine estimation method comprises the steps of firstly, constructing an identification function linearly related to carrier frequency according to two values (the maximum value and the maximum value in the adjacent values) in an FFT output amplitude spectrum sequence; then, designing an estimation threshold value according to the statistical characteristics of the FFT output sequence, and judging whether the radar carrier is captured or not by judging whether the two values exceed the estimation threshold value or not; and finally, calculating to obtain the compensation quantity of the FFT process to the carrier frequency estimation according to the constructed linear correlation discrimination function, thereby improving the estimation precision of the radar carrier frequency. The invention has simple structure and small calculation amount, can greatly improve the estimation precision of the carrier frequency by adding the carrier frequency fine estimation module after the carrier frequency estimation structure in the original FFT process, and has higher engineering use value.
Drawings
FIG. 1 is a flow chart of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
As used in this application and the appended claims, the terms "a," "an," "the," and/or "the" are not intended to be inclusive in the singular, but rather are intended to be inclusive in the plural unless the context clearly dictates otherwise. In general, the terms "comprises" and "comprising" merely indicate that steps and elements are included which are explicitly identified, that the steps and elements do not form an exclusive list, and that a method or apparatus may include other steps or elements.
The relative arrangement of the components and steps, the numerical expressions, and numerical values set forth in these embodiments do not limit the scope of the present application unless specifically stated otherwise. Meanwhile, it should be understood that the sizes of the respective portions shown in the drawings are not drawn in an actual proportional relationship for the convenience of description. Techniques, methods, and apparatus known to those of ordinary skill in the relevant art may not be discussed in detail but are intended to be part of the specification where appropriate. In all examples shown and discussed herein, any particular value should be construed as merely illustrative, and not limiting. Thus, other examples of the exemplary embodiments may have different values. It should be noted that: like reference numbers and letters refer to like items in the following figures, and thus, once an item is defined in one figure, further discussion thereof is not required in subsequent figures.
Referring to fig. 1, the present embodiment provides a method for fine estimation of a radar carrier frequency, including the following steps:
step 1: two values (maximum value M) in an FFT (Fast Fourier Transform) output amplitude spectrum sequence S (k) of a radar carrier signal are extractedFFTAnd the maximum value S of the values adjacent theretoubm);
Specifically, the structure of the radar carrier signal s (t) after frequency reduction is shown in formula (1),
where A is the amplitude of the carrier signal, i is the imaginary unit, Δ f is the down-converted carrier frequency,is the carrier initial phase.
After the radar carrier signal s (t) subjected to frequency reduction is dispersed and cut, as shown in formula (2),
where (·) is a Dirac Delta function, TsIs the sampling period, rect () is a rectangular window function, TLIs the size of a rectangular window.
The formula (2) is abbreviated as:
wherein d (t) is a sampling pulse signal; w (t) is a rectangular window signal;
the spectrum of S (n) is obtained according to the time domain multiplication frequency domain convolution principle as shown in the following formula,
in the formula, denotes a convolution symbol.
The amplitude spectrum sequence of the DTFT (Discrete-time Fourier Transform) output by expanding each term in the formula (4) is shown as the following formula,
wherein, is a convolution symbol;is the frequency spectrum of a discrete finite long time domain signal s (n); (f) is the frequency spectrum of the continuous infinite time domain signal S (t); d (f) is the frequency spectrum of the continuous infinite time domain signal d (t); w (f) is the frequency spectrum of a continuous infinite time domain signal w (t);
according to FFT principle, continuous infinite length frequency spectrum obtained by DTFTPerforming dispersion and truncation to obtain an amplitude spectrum sequence output by the FFT:
as can be seen from equation (6), the maximum value of the FFT output magnitude spectrum sequence can be represented by equation (7),
in the FFT output magnitude spectrum sequence, one of the points (or left point, or right point) adjacent to the FFT output maximum value point is only smaller than the maximum value, which can be expressed by equation (8),
where round (. cndot.) is an integer function, MFFTRepresenting the maximum value, S, of the FFT output amplitude spectrum sequenceubmRepresents and MFFTThe larger of the neighboring points.
Step 2: for the two values (M)FFTAnd Subm) Performing an analysis by comparison with an estimated threshold value YtComparing and judging MFFTAnd SubmWhether a condition that can be used to improve the estimation accuracy of the FFT on the frequency is satisfied or not is used to attenuate the influence of noise on the method;
wherein the threshold value Y is estimatedtThe calculation formula (10) for estimating the threshold value can be derived from the following calculation, in which the noise mixed in the FFT magnitude spectrum belongs to the rayleigh distribution and the FFT magnitude output belongs to the rice distribution, according to the formula (9).
In the formula (f)n(. is a probability density function of the Rayleigh distribution, PfaIs the false alarm rate, Yt is the estimated threshold, and σ is the noise statistical standard deviation.
Will be obtained by outputting the maximum M in the amplitude spectrum sequence S (k) by FFTFFTAnd the maximum value S of the values adjacent theretoubmRespectively comparing with the estimated threshold value Yt, if both are greater than Yt, then M is indicatedFFTAnd SubmCan be used to improve the estimation accuracy of FFT to frequency, and then proceed the subsequent steps, if M isFFTIf the value is less than or equal to Yt, the carrier wave is not detected, the calculation is finished, and if S is less than or equal to Yt, the calculation is finishedubmIf it is less than or equal to Yt, then S is indicatedubmCan not be used for improving the frequency estimation precision of FFT, and the calculation is finished。
And step 3: if M is as aboveFFTAnd SubmBoth values are greater than the threshold value Yt, then the two values (M) are usedFFTAnd Subm) And calculating to obtain the compensation quantity of the FFT process to the carrier frequency estimation, thereby improving the estimation precision of the radar carrier frequency.
Specifically, a linear discrimination function D associated with the down-converted carrier frequency Δ f can be constructed according to equation (8)isc(. cndot.) represented by formula (9).
In the formula, | · | is an absolute value symbol.
If S isubmAnd MFFTAdjacent positions in the FFT output amplitude sequence, the FFT estimation compensation formula can be obtained by formula (11),
where Δ f' is a compensation value for the carrier frequency estimation in the FFT process.
And 4, step 4: judgment SubmWhether or not at M at the position of the FFT output amplitude sequenceFFTOn the left-hand side of the position,
if on the right, it can be expressed by equation (14),
in the formula (f)FFTFor the FFT process to estimate the carrier frequency,a fine estimate of the compensated carrier frequency.
The invention has simple structure and small calculation amount, can greatly improve the estimation precision of the carrier frequency by adding the carrier frequency fine estimation module after the carrier frequency estimation structure in the original FFT process, and has higher engineering use value.
The above description is only an embodiment of the present invention, and not intended to limit the scope of the present invention, and all modifications of equivalent structures and equivalent processes performed by the present specification and drawings, or directly or indirectly applied to other related technical fields, are included in the scope of the present invention.
Claims (5)
1. A method for fine estimation of a radar carrier frequency, comprising:
extracting maximum value M in FFT output amplitude spectrum sequence S (k) of radar carrier signalFFTAnd the maximum value S of the values adjacent theretoubm;
Respectively judging the maximum values MFFTAnd the maximum value S of the values adjacent to itubmAnd an estimated threshold value YtIf both are larger than the estimated threshold value YtIf yes, continuing to estimate, otherwise, ending;
the compensation value of the FFT process for the carrier frequency estimation is obtained according to the formula (1),
wherein, Deltaf' is the compensation value of the FFT process to the carrier frequency estimation, Deltaf is the carrier frequency after frequency reduction, TLIs the size of the rectangular window; mFFTOutputting a maximum value in the amplitude spectrum sequence for the FFT; submIs the maximum value among the adjacent values of the maximum values;
and obtaining a carrier frequency fine estimation value of the radar carrier according to the compensation value.
2. The fine estimation method according to claim 1, wherein the maximum value M in the extracted radar carrier signal FFT output magnitude spectrum sequence s (k)FFTAnd the maximum value S of the values adjacent theretoubmFurther comprising:
the structure of the radar carrier signal s (t) after frequency reduction is shown in formula (2),
where A is the amplitude of the carrier signal, i is the imaginary unit, Δ f is the down-converted carrier frequency,is the carrier initial phase;
after the radar carrier signal s (t) subjected to frequency reduction is dispersed and cut, as shown in formula (3),
where (·) is a Dirac Delta function, TsIs the sampling period, rect () is a rectangular window function, TLIs the size of the rectangular window;
the formula (3) is abbreviated as:
wherein d (t) is a sampling pulse signal; w (t) is a rectangular window signal;
according to the principle of time-domain multiplied frequency-domain convolution, s (n) is transformed into the frequency domain by a discrete-time fourier transform, as shown in equation (5):
wherein, is a convolution symbol;is the frequency spectrum of a discrete finite long time domain signal s (n); (f) is the frequency spectrum of the continuous infinite time domain signal S (t); d (f) is the frequency spectrum of the continuous infinite time domain signal d (t); w (f) is the frequency spectrum of a continuous infinite time domain signal w (t);
each term in the formula (5) is expanded to obtain, the amplitude spectrum sequence of the DTFT output is shown as the following formula,
where sinc (·) is a sinc function;
according to FFT principle, continuous infinite length frequency spectrum obtained by DTFTPerforming dispersion and truncation to obtain an amplitude spectrum sequence output by the FFT:
as can be seen from equation (7), the maximum value of the FFT output magnitude spectrum sequence can be represented by equation (8),
in the FFT output amplitude spectrum sequence, one value among the points adjacent to the maximum value point of the FFT output is only smaller than the maximum value, as shown in equation (9),
where round (. cndot.) is an integer function, MFFTRepresenting the maximum value, S, of the FFT output amplitude spectrum sequenceubmRepresents and MFFTThe larger of the neighboring points.
3. The fine estimation method according to claim 1, characterized in that the estimation threshold value YtThe method comprises the following steps:
the noise mixed in the FFT amplitude spectrum belongs to Rayleigh distribution, and the FFT amplitude output belongs to Rice distribution, so that a calculation formula (11) for estimating the threshold value can be deduced according to a formula (10);
in the formula (f)n(. is a probability density function of the Rayleigh distribution, PfaIs the false alarm rate, Yt is the estimated threshold, and σ is the noise statistical standard deviation.
4. The fine estimation method according to claim 2, wherein the deriving the compensation value for the carrier frequency estimation by the FFT process according to equation (1) further comprises:
a linear discrimination function D related to the down-converted carrier frequency Deltaf is constructed according to the formula (9)isc(. cndot.) represented by formula (12).
In the formula, | · | is an absolute value symbol.
If SubmAnd MFFTAdjacent positions in the FFT output amplitude sequence, the FFT estimation compensation formula can be obtained by formula (1),
where Δ f' is a compensation value for the carrier frequency estimation in the FFT process.
5. The fine estimation method of claim 4, wherein the deriving a radar carrier frequency fine estimate from the compensation value further comprises:
if the maximum value S among the adjacent valuesubmAt a maximum M at the position of the FFT output amplitude sequenceFFTTo the left of the position, the estimated carrier frequency estimateCan be expressed by formula (13); if on the right, it can be expressed by equation (14),
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010829967.6A CN112014811B (en) | 2020-08-18 | 2020-08-18 | Fine estimation method for radar carrier frequency |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010829967.6A CN112014811B (en) | 2020-08-18 | 2020-08-18 | Fine estimation method for radar carrier frequency |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112014811A true CN112014811A (en) | 2020-12-01 |
CN112014811B CN112014811B (en) | 2023-11-07 |
Family
ID=73504872
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010829967.6A Active CN112014811B (en) | 2020-08-18 | 2020-08-18 | Fine estimation method for radar carrier frequency |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112014811B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117518108A (en) * | 2024-01-04 | 2024-02-06 | 江苏方测建筑工程技术有限公司 | High-precision measurement and compensation method for carrier frequency difference between airborne double-station radars |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1626549A2 (en) * | 2004-08-11 | 2006-02-15 | Rohde & Schwarz GmbH & Co. KG | Method and apparatus for estimating carrier frequency using a block-wise coarse estimation |
CN101666874A (en) * | 2009-09-07 | 2010-03-10 | 北京科技大学 | Measuring device for accurate correction by FMCW radar distance measurement and compensation correction method |
CN106018956A (en) * | 2016-08-10 | 2016-10-12 | 北京妙微科技有限公司 | Power system frequency calculation method of windowing spectral line interpolation |
CN106154257A (en) * | 2016-06-21 | 2016-11-23 | 成都雷尼尔科技有限公司 | Precision instrumentation radar secondary frequency measuring method based on FFT Yu apFFT |
CN107064629A (en) * | 2017-06-07 | 2017-08-18 | 东南大学 | It is a kind of that simple signal frequency estimating methods are integrated based on the segmentation that frequency relative bias is estimated |
CN110691051A (en) * | 2019-09-29 | 2020-01-14 | 天津大学 | GMSK signal frequency offset estimation algorithm based on FFT |
CN111537969A (en) * | 2020-05-25 | 2020-08-14 | 电子科技大学 | Broadband instantaneous frequency measuring and tracking device |
-
2020
- 2020-08-18 CN CN202010829967.6A patent/CN112014811B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1626549A2 (en) * | 2004-08-11 | 2006-02-15 | Rohde & Schwarz GmbH & Co. KG | Method and apparatus for estimating carrier frequency using a block-wise coarse estimation |
CN101666874A (en) * | 2009-09-07 | 2010-03-10 | 北京科技大学 | Measuring device for accurate correction by FMCW radar distance measurement and compensation correction method |
CN106154257A (en) * | 2016-06-21 | 2016-11-23 | 成都雷尼尔科技有限公司 | Precision instrumentation radar secondary frequency measuring method based on FFT Yu apFFT |
CN106018956A (en) * | 2016-08-10 | 2016-10-12 | 北京妙微科技有限公司 | Power system frequency calculation method of windowing spectral line interpolation |
CN107064629A (en) * | 2017-06-07 | 2017-08-18 | 东南大学 | It is a kind of that simple signal frequency estimating methods are integrated based on the segmentation that frequency relative bias is estimated |
CN110691051A (en) * | 2019-09-29 | 2020-01-14 | 天津大学 | GMSK signal frequency offset estimation algorithm based on FFT |
CN111537969A (en) * | 2020-05-25 | 2020-08-14 | 电子科技大学 | Broadband instantaneous frequency measuring and tracking device |
Non-Patent Citations (3)
Title |
---|
孙宏军;王小威;: "基于幅值-相角判据的修正Rife正弦波频率估计算法", 天津大学学报(自然科学与工程技术版), no. 08 * |
李小捷;许录平;郑文松;: "GPS信号的码相位与多普勒频率联合捕获", 宇航学报, no. 12 * |
齐国清: "利用FFT相位差校正信号频率和初相估计的误差分析", 数据采集与处理, no. 01 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117518108A (en) * | 2024-01-04 | 2024-02-06 | 江苏方测建筑工程技术有限公司 | High-precision measurement and compensation method for carrier frequency difference between airborne double-station radars |
CN117518108B (en) * | 2024-01-04 | 2024-04-16 | 江苏方测建筑工程技术有限公司 | High-precision measurement and compensation method for carrier frequency difference between airborne double-station radars |
Also Published As
Publication number | Publication date |
---|---|
CN112014811B (en) | 2023-11-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106845010B (en) | Low-frequency oscillation dominant mode identification method based on improved SVD noise reduction and Prony | |
CN106597408B (en) | High-order PPS signal parameter estimation method based on time-frequency analysis and instantaneous frequency curve fitting | |
Raykar et al. | Fast computation of kernel estimators | |
CN108344895B (en) | Signal detection method and device | |
CN105785324B (en) | Linear frequency-modulated parameter estimating method based on MGCSTFT | |
CN102222911A (en) | Power system interharmonic estimation method based on auto-regression (AR) model and Kalman filtering | |
CN103684350B (en) | A kind of particle filter method | |
CN104901909B (en) | The method for parameter estimation of chirp signals under a kind of α non-Gaussian noises | |
Ozkurt et al. | Determination of wavelet ridges of nonstationary signals by singular value decomposition | |
CN106209703A (en) | A kind of Frequency Hopping Signal Blind Parameter Estimation and device | |
Kuang et al. | Efficient and adaptive signal denoising based on multistage singular spectrum analysis | |
Dehkordi et al. | Sound source localization using compressive sensing-based feature extraction and spatial sparsity | |
CN112014811B (en) | Fine estimation method for radar carrier frequency | |
KR20160116440A (en) | SNR Extimation Apparatus and Method of Voice Recognition System | |
Najeeb et al. | Review of parameter estimation techniques for time-varying autoregressive models of biomedical signals | |
CN113468474B (en) | Power grid frequency estimation method based on root Mini-Norm | |
Eappen et al. | Semiparametric algorithm for processing MST radar data | |
CN115378776A (en) | MFSK modulation identification method based on cyclic spectrum parameters | |
Morgan | On level-crossing excursions of Gaussian low-pass random processes | |
Romulus | A comparison between instantaneous frequency estimation methods of frequency modulated signals covered with gaussian noise | |
Yin et al. | Noise reduction method for vibration signals 2D time‐frequency distribution using anisotropic diffusion equation | |
Song | Random signal frequency identification based on AR model spectral estimation | |
CN115795233B (en) | Method, device and storage medium for extracting time delay term of frequency domain data | |
Chen et al. | An improved method of wavelet threshold denoising for satellite signal | |
CN106980043A (en) | A kind of improvement phase difference correction method based on Hanning window |
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 |