CN106872773A - A kind of the multiple-pulse Precision Method of Freuqency Measurement and device of single carrier frequency pulse signal - Google Patents

A kind of the multiple-pulse Precision Method of Freuqency Measurement and device of single carrier frequency pulse signal Download PDF

Info

Publication number
CN106872773A
CN106872773A CN201710278139.6A CN201710278139A CN106872773A CN 106872773 A CN106872773 A CN 106872773A CN 201710278139 A CN201710278139 A CN 201710278139A CN 106872773 A CN106872773 A CN 106872773A
Authority
CN
China
Prior art keywords
signal
sample
pulse
frequency
complex
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
CN201710278139.6A
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.)
CETC 2 Research Institute
Southwest China Research Institute Electronic Equipment
Original Assignee
CETC 2 Research Institute
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 CETC 2 Research Institute filed Critical CETC 2 Research Institute
Priority to CN201710278139.6A priority Critical patent/CN106872773A/en
Publication of CN106872773A publication Critical patent/CN106872773A/en
Pending legal-status Critical Current

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
    • G01R23/12—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage by converting frequency into phase shift
    • 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)
  • Mathematical Physics (AREA)
  • Measuring Frequencies, Analyzing Spectra (AREA)

Abstract

The invention discloses a kind of the multiple-pulse Precision Method of Freuqency Measurement and device of single carrier frequency pulse signal, it is related to the Precision Method of Freuqency Measurement of pulse signal.Technical key point includes：Step 1：The real number list carrier frequency pulse string signal sample that will be received carries out orthogonal transformation；The 1st Complex pulse sample of signal is processed using Kay algorithms, obtains frequency iteration initial value fsp；Step 2：Fourier transformation is carried out to each Complex pulse sample of signal, the phase difference of the peak point after each Complex pulse sample of signal Fourier transformation is asked for；Step 3：Time difference between calculating two neighboring Complex pulse sample of signal and carrying out the starting sample point of Fourier transformation；Step 4：Be iterated be calculated f (2), f (3) ..., f (K)；F (K) is multiple-pulse frequency precise measurements of the real number list carrier frequency pulse string signal sample etc..

Description

A kind of the multiple-pulse Precision Method of Freuqency Measurement and device of single carrier frequency pulse signal
Technical field
It is more particularly to a kind of to single carrier frequency pulse signal the present invention relates to the Precision Method of Freuqency Measurement of pulse signal Multiple-pulse Precision Method of Freuqency Measurement.
Background technology
Steven Kay were in 1989《IEEE Transactions on Acoustics.Speech.And Signal Processing》On delivered entitled《A Fast and Accurate Single Frequency Estimator》Opinion Text, in white Gaussian noise environment using phase weighting method enters line frequency accurate the paper proposes one kind to single CF signal The algorithm of measurement, here it is current widely used Kay algorithms.Single carrier frequency continuous wave is believed using phase weighting method in paper Number phase weighting treatment is carried out, continuous wave signal frequency on variance least meaning can be obtained under conditions of high signal/noise ratio Unbiased measured value.
When carrying out single CF signal frequency measurement using Kay algorithms, the precision of its frequency measurement adds with line phase is entered The sample of signal length of power is relevant, within the specific limits as the sample of signal of phase weighting is more long, what Kay algorithms can be obtained Frequency measurement accuracy is higher.Therefore, for the frequency accurate measurement requirement of single carrier frequency continuous wave signal, Kay algorithms can pass through Increase the sample of signal length of phase weighting to obtain required frequency precise measurements.However, when being directed to pulse signal, due to Its pulsewidth is fixed, and the sample of signal length for carrying out phase weighting can not possibly exceed pulsewidth, therefore Kay algorithms can reach most Good frequency measurement accuracy has been restricted, it is impossible to further lift frequency measurement accuracy.In order to solve Kay algorithms to list Carrier frequency pulse signal cannot further lift the problem of frequency measurement accuracy when carrying out frequency measurement, the present invention proposes a kind of right The multiple-pulse Precision Method of Freuqency Measurement of single carrier frequency pulse signal, the method can further be lifted on the basis of Kay algorithms The precision of frequency measurement.
The content of the invention
Carry out further lifting frequency measurement accuracy when single carrier frequency pulse signal frequency is measured to solve Kay algorithms Problem, the present invention on the basis of pulse frequency measurement utilize multiple-pulse iteration, realize to single carrier frequency pulse signal High accuracy frequency accurate measurement.
The multiple-pulse Precision Method of Freuqency Measurement of single carrier frequency pulse signal that the present invention is provided, including：
Step 1：The real number list carrier frequency pulse string signal sample that will be received carries out orthogonal transformation, obtains K Complex pulse Sample of signal；The 1st Complex pulse sample of signal is processed using Kay algorithms, obtains frequency iteration initial value fsp；Real number Single carrier frequency pulse string signal sample includes K real number list carrier frequency pulse signal；K is the integer more than 1；
Step 2：Fourier transformation is carried out to each Complex pulse sample of signal, each Complex pulse signal sample is asked for The phase value of the peak point after this Fourier transformation, so as to obtain K phase value φ (k), k=1,2 ... ..., K；By adjacent two Individual phase value subtracts each other φ (k)-φ (k-1) respectively, and k takes 2 herein ... ..., K, so as to obtain K-1 phase difference φ (2), Δ φ (3) ... ..., Δ φ (K)；
Step 3：Calculate two neighboring Complex pulse sample of signal carry out between the starting sample point of Fourier transformation when Between difference t0(k)-t0(k-1), t0K () represents that the starting sample point of k-th Fourier transformation of Complex pulse sample of signal is corresponding Time, herein k take 2 ... ..., K obtains Δ t (2), Δ t (3) ... ..., Δ t (K)；
Step 4：Iterative calculation below equation K-1 times successively, obtain f (2), f (3) ..., f (K)：
Φ (n)=2 π f (n) Δ t (n+1)-Δ φ (n+1)；
Φ ' (n)=2m (n) π+Δ φ (n+1)；
Wherein, n take 1 successively, 2,3 ..., K-1；F (1)=fsp, [] represents and asks for integer according to the principle for rounding up；
F (K) is the multiple-pulse frequency precise measurements of the real number list carrier frequency pulse string signal sample.
Preferably, step 1 is further included：
Step 11：The phase value of whole sampled points is asked in the 1st Complex pulse sample of signal, by two neighboring sampling The phase value of point subtracts each other and obtainsWherein i=1,2 ... ..., L, represent adjacent two in the 1st Complex pulse sample of signal The order label of the phase difference of individual sampled point, L represents the sum of the two neighboring sampled point phase difference of Complex pulse sample of signal；
Step 12：According to formulaAsk for L phase weighting coefficients, wherein i= 1,2 ... ..., L；Computing formulaCan obtain the 1st simple venation of Complex pulse sample of signal Rush frequency measurement fsp, wherein tsIt is the sampling interval duration of Complex pulse sample of signal, i=1,2 ... ..., L.
The multiple-pulse frequency device for accurately measuring of a kind of single carrier frequency pulse signal that the present invention is provided, including：
Frequency iteration calculation of initial value unit, the real number list carrier frequency pulse string signal sample for that will receive carries out orthogonal Conversion, obtains K Complex pulse sample of signal；The 1st Complex pulse sample of signal is processed using Kay algorithms, is obtained Frequency iteration initial value fsp；Real number list carrier frequency pulse string signal sample includes K real number list carrier frequency pulse signal；K is more than 1 Integer；
Fourier transformation peak point phase difference calculating unit, for carrying out Fourier to each Complex pulse sample of signal Conversion, asks for the phase value of the peak point after each Complex pulse sample of signal Fourier transformation, so as to obtain K phase value φ (k), k=1,2 ... ..., K；Two neighboring phase value is subtracted each other into φ (k)-φ (k-1) respectively, k takes 2 herein ... ..., K, from And K-1 phase difference φ (2) is obtained, Δ φ (3) ... ..., Δ φ (K)；
Fourier transformation starting sample point time difference calculating unit, enters for calculating two neighboring Complex pulse sample of signal Time difference t between the starting sample point of row Fourier transformation0(k)-t0(k-1), t0K () represents k-th Complex pulse signal sample The starting sample point corresponding time of this Fourier transformation, herein k take 2 ... ..., K obtains Δ t (2), Δ t (3) ... ..., Δt(K)；
Multiple-pulse frequency precise measurements iterate to calculate unit, for iterating to calculate below equation K-1 times successively, obtain f (2)、f(3)、…、f(K)：
Φ (n)=2 π f (n) Δ t (n+1)-Δ φ (n+1)；
Φ ' (n)=2m (n) π+Δ φ (n+1)；
Wherein, n take 1 successively, 2,3 ..., K-1；F (1)=fsp, [] represents and asks for integer according to the principle for rounding up；
F (K) is the multiple-pulse frequency precise measurements of the real number list carrier frequency pulse string signal sample.
Preferably, frequency iteration calculation of initial value unit is further included：
Signal phase difference computing unit, the phase for asking for whole sampled points in the 1st Complex pulse sample of signal Value, the phase value of two neighboring sampled point is subtracted each other and is obtainedWherein i=1,2 ... ..., L, represent the 1st Complex pulse The order label of the phase difference of two neighboring sampled point in sample of signal, L represents the two neighboring sampling of Complex pulse sample of signal The sum of point phase difference；
Phase weighting unit, for according to formulaAsk for L phase weighting system Number, wherein i=1,2 ... ..., L；Computing formulaCan obtain the 1st Complex pulse signal The pulse frequency measurement f of samplesp, wherein tsIt is the sampling interval duration of Complex pulse sample of signal, i=1,2 ... ..., L。
The present invention, using multiple-pulse iteration method, solves to be directed on the basis of Kay algorithm pulse phase weightings Kay algorithms cannot further lift the problem of frequency measurement accuracy during single carrier frequency pulse signal.For single carrier frequency pulse signal pair As when, the inventive method needs more pulse signal samples compared to Kay algorithms, but can increase substantially frequency measurement essence Degree.
It is 1us, SNR=15dB in single carrier frequency pulse signal width, the present invention uses 6 pulse signal samples, enters respectively 1000 Monte Carlo experiments of row, result of the test is as shown in accompanying drawing 2 and accompanying drawing 3；The pulse frequency error measurement of Kay algorithms is equal Between -200Hz~200Hz, frequency error measurement mean square deviation is between 300Hz~450Hz for value；Under similarity condition, the present invention The multiple-pulse frequency error measurement average of method between -0.4Hz~0.4Hz, frequency error measurement mean square deviation 0Hz~ Between 1.5Hz.
In summary, when carrying out frequency measurement to single carrier frequency pulse signal, the inventive method compares Kay algorithms, although needing The signal pulse number wanted is more, but can but obtain precision frequency precise measurements higher, and the inventive method is adapted to application In the situation higher to frequency measurement accuracy requirement.
Brief description of the drawings
Examples of the present invention will be described by way of reference to the accompanying drawings, wherein：
Fig. 1 is the inventive method FB(flow block).
Fig. 2 is the error mean and mean square deviation of Kay algorithm pulse frequency measurements.
Fig. 3 is the error mean and mean square deviation of the inventive method frequency measurement.
Specific embodiment
All features disclosed in this specification, or disclosed all methods or during the step of, except mutually exclusive Feature and/or step beyond, can combine by any way.
Any feature disclosed in this specification, unless specifically stated otherwise, can be equivalent or with similar purpose by other Alternative features are replaced.I.e., unless specifically stated otherwise, each feature is an example in a series of equivalent or similar characteristics .
Referring to Fig. 1, the inventive method includes：
Orthogonal transformation step：The real number list carrier frequency arteries and veins being made up of K real number list carrier frequency pulse signal first to receiving Rushing string signal sample carries out orthogonal transformation, so that the real number burst signal sample is converted into Complex pulse string signal sample This, is easy to extraction and treatment of the subsequent step to sample of signal phase.
Complex pulse string signal sample phase difference asks for step：Complex pulse string signal from after orthogonal transformation step process Whole Complex pulse sample of signal are extracted in sample, altogether comprising K Complex pulse letter in the Complex pulse string signal sample Number sample.The phase value of whole sampled points is asked in the 1st Complex pulse sample of signal, by whole two neighboring sampled points Phase value subtract each other, the phase difference of as two neighboring sampled pointWherein i=1,2 ... ..., L represent the 1st plural number The order label of the phase difference of two neighboring sampled point in pulse signal sample, L represents that Complex pulse sample of signal is two neighboring The sum of sampled point phase difference.
Phase weighting step：The number L of the phase difference of two neighboring sampled point is in 1st Complex pulse sample of signal The number of phase weighting coefficients, according to formulaWhole L phase weighting coefficients are asked for, Wherein i=1,2 ... ..., L represent the order label of phase weighting coefficients, and the label is corresponded with the order label of phase difference.Will be preceding Whole L phase differences of the 1st Complex pulse sample of signal that one step is obtainedWith order label identical phase weighting Coefficient w (i) correspondence is multiplied, and L product is all added up and divided by 2 π ts, i.e., according to formula Processed, obtained the 1st pulse frequency measurement f of Complex pulse sample of signalsp, wherein tsIt is Complex pulse signal sample This sampling interval duration, i=1,2 ... ..., L, fspAlso it is the pulse frequency measurement obtained using Kay algorithms, subsequently by it As the initial value that frequency iteration is calculated.
The phase difference calculating step of Fourier transformation peak point：By each Complex pulse after orthogonal transformation step process Sample of signal carries out Fourier transformation, asks for the phase value of peak point after Fourier transformation, so as to obtain k-th Complex pulse letter Number Fourier transformation peak point phase value φ (k), k=1,2 ... ..., K.Will be adjacent in whole K Complex pulse sample of signal Two Complex pulse sample of signal Fourier transformation peak point phases are subtracted each other respectively, i.e. Δ φ (2)=φ (2)-φ (1), Δ φ (3)=φ (3)-φ (2) ..., Δ φ (K)=φ (K)-φ (K-1), so as to obtain the phase difference of peak point after Fourier transformation Δ φ (2), Δ φ (3) ... ..., Δ φ (K).
The time difference calculation procedure of Fourier transformation starting sample point：Calculate the Complex pulse signal subtracted each other in previous step Sample carries out the time difference between the starting sample point of Fourier transformation, that is, calculate Δ t (2)=t0(2)-t0(1), Δ t (3)= t0(3)-t0(2) ..., Δ t (K)=t0(K)-t0(K-1), t0K () represents k-th Fourier transformation of Complex pulse sample of signal The starting sample point corresponding time, herein k take 2 ... ..., K obtains K-1 time difference Δ t (2), Δ t (3) ... ..., Δ t (K)。
Iterative step first：Frequency of the pulse frequency measurement that phase weighting step is obtained as frequency accurate measurement Rate iteration initial value f (1)=fsp, the Δ φ (2) and Δ t (2) for then obtaining the first two steps accurately survey as frequency The Fourier transformation peak phase difference initial value and time difference initial value of amount, they are all brought into formula Φ (1)=2 π f first (1) Δ t (2)-Δ φ (2) andIn, the 2nd Complex pulse letter that wherein Φ (1) the 1st iteration of expression is obtained The phase difference value for subtracting the small several times of 2 π of the phase difference comprising error, m number between sample and the 1st Complex pulse sample of signal (1) phase difference between the 2nd Complex pulse sample of signal and the 1st Complex pulse sample of signal that the 1st iteration obtain is represented Take 2 π as the integral multiple periodic quantity in cycle, wherein, [] represents and asks for integer according to the principle for rounding up.
M (1) and Δ φ (2) are brought into formula Φ ' (1)=2m (1) π+Δ φ (2), wherein Φ ' (1) is represented and included 2 π 2nd Complex pulse sample of signal Fourier transformation starting sample point of decimal multiple and the 1st Complex pulse sample of signal Fu In phase difference between leaf transformation starting sample point.Bring Φ ' (1) and Δ t (2) into formula againIn, so that Ask for the frequency measurement f (2) that iteration first is obtained.
Multiple-pulse iteration step：Complex pulse string signal sample includes K Complex pulse sample of signal altogether, therefore altogether Carry out K-1 iteration.First after iterative step, next in nth iteration (wherein, n=2,3 ... ..., K-1), will above Iterative step parameter f (n), Δ t (n), Δ φ (n) that obtain substitute into equation below respectively：Φ (n)=2 π f (n) Δ t (n+ 1)-Δφ(n+1)、Φ ' (n)=2m (n) π+Δ φ (n+1) andIn repetition Process is stated, after completing whole K-1 iteration, frequency measurement f (K) can be obtained.Frequency measurement f (K) is what is received The multiple-pulse frequency precise measurements of real number list carrier frequency pulse string signal sample, the result is exported and respective application is carried out.
Present invention also offers a kind of and one-to-one floppy disk system of above method step.
In order to be expressly recited the implementation process and validity of the inventive method, an embodiment is herein proposed：From list The real number train of pulse of carrier frequency is processed as sample of signal according to the inventive method flow shown in Fig. 1, real number list carrier frequency arteries and veins Rushing string signal sample parameter is：Pulsewidth be 1us, the pulse repetition period be 60us, signal to noise ratio be 15dB, sample rate be 2GHz, in Frequency of heart is 100MHz, and burst signal sample includes 6 pulse signal samples.
First, to frequency for single carrier frequency real number burst signal sample of 100MHz carries out orthogonal transformation, this step terminates After obtain Complex pulse string signal sample.
The Complex pulse string signal sample extraction obtained from previous step goes out whole K=6 Complex pulse sample of signal, and Each Complex pulse sample of signal is processed respectively.Each pulse is 1us, and sample rate is 2GHz, therefore each is multiple Rapid pulse rushes sample of signal and has 2000 sampled points, and now all the total L=1999 of the phase difference of two neighboring sampled point is individual.It is right 6 Complex pulse sample of signal ask for phase and ask for the phase difference of two neighboring sampled point in sequence respectively, each plural number Pulse signal sample can obtain 1999 phase differences of two neighboring sampled point, be expressed asWherein i represents phase difference label, meets i=1, 2 ... ..., 1999, leftover bits and pieces is designated as the serial number of Complex pulse signal.During actual calculating, only needed in subsequent step of the invention Use the 1st phase difference of Complex pulse sample of signal.
L=1999 is brought into formulaIn, wherein, i=1, in 2 ... ..., L. So as to try to achieve all 1999 phase weighting coefficients w (1), w (2) ... ..., w (1999).In the 1st Complex pulse sample of signal Middle identical phase weighting coefficients w (i) that will have been tried to achieve respectively and the corresponding Complex pulse numbered with same sequence The phase difference of the two neighboring sampled point of sample of signalIt is brought into respectivelyIn, wherein ts=1/Fs, sample rate Fs=2GHz.The pulse frequency measurement f of the utilization Complex pulse sample of signal for now obtainingsp= 100.000469MHz.The value is the frequency measurement of single carrier frequency pulse sample of signal that Kay algorithms are obtained, now Kay algorithms Frequency error measurement be 469Hz.
1024 point quick Fourier conversion are carried out respectively to 6 Complex pulse sample of signal, the peak point after conversion is asked for Phase simultaneously asks for peak point phase difference respectively, and phase difference is respectively Δ φ (2), Δ φ (3), Δ φ (4), Δ φ (5), Δ φ (6)。
Time difference between calculating different Complex pulse sample of signal and carrying out the starting sample point of Fourier transformation, respectively Δt(2),Δt(3),Δt(4),Δt(5),Δt(6)。
F (1)=f that step before is asked forsp=100.000469MHz, Δ φ (2), Δ t (2) difference iteration initial values, Respectively substitute into formula Φ (1)=2 π f (1) Δ t (2)-Δ φ (2) andIn, then Φ (1) and m (1) is substituted into To in formula Φ ' (1)=2m (1) π+Δ φ (2), then itself and Δ t (2) are brought intoIn, can obtain first Frequency measurement f (the 2)=99.9999978MHz of iteration.
K-1=5 iteration is carried out altogether, and the frequency measurement for obtaining is respectively：F (3)=99.9999996MHz, f (4)=99.9999996MHz, f (5)=99.9999996MHz, f (6)=99.9999996MHz.Therefore, whole iteration terminate The frequency measurement f (6)=99.9999996 for obtaining afterwards is the frequency measurement that the inventive method is obtained.In summary, originally The frequency error measurement of inventive method is 0.4Hz, and contrast understands the inventive method compared to Kay algorithms in terms of frequency error measurement There is greatly lifting.
F (6)=99.9999996MHz is the real number list carrier frequency pulse string signal sample for receiving that the inventive method is obtained This multiple-pulse frequency precise measurements, the result is exported and respective application is carried out.
A frequency values are taken as real number list carrier frequency every 10MHz in the range of 100MHz~900MHz according to as above method The frequency of burst signal sample, other conditions are constant, and being utilized respectively Kay algorithms and the inventive method carries out pulse frequency survey Amount and multiple-pulse frequency accurate measurement, carry out 1000 Monte Carlo experiments respectively, can obtain Kay algorithms and the inventive method Frequency measurement error mean and mean square deviation, respectively as shown in accompanying drawing 2 and accompanying drawing 3.It is visible in figure, obtained using Kay algorithms The frequency measurement error mean for arriving in the range of -200Hz~200Hz, the error mean square difference of frequency measurement 300Hz~ In the range of 450Hz；The frequency measurement error mean obtained using the inventive method in the range of -0.4Hz~0.4Hz, frequency The error mean square difference of measured value is in the range of 0Hz~1.5Hz.
Result above shows：When carrying out single carrier frequency pulse signal frequency accurate measurement, though the inventive method compares Kay algorithms The pulse signal sample for so using is more, but can be significantly to lift frequency measurement accuracy, so as to realize to single carrier frequency pulse The purpose of the multiple-pulse frequency accurate measurement of signal.
The invention is not limited in foregoing specific embodiment.The present invention is expanded to and any in this manual disclosed New feature or any new combination, and disclose any new method or process the step of or any new combination.

Claims (4)

1. a kind of multiple-pulse Precision Method of Freuqency Measurement of single carrier frequency pulse signal, it is characterised in that including：
Step 1：The real number list carrier frequency pulse string signal sample that will be received carries out orthogonal transformation, obtains K Complex pulse signal Sample；The 1st Complex pulse sample of signal is processed using Kay algorithms, obtains frequency iteration initial value fsp；Real number list is carried Frequency burst signal sample includes K real number list carrier frequency pulse signal；K is the integer more than 1；
Step 2：Fourier transformation is carried out to each Complex pulse sample of signal, each Complex pulse sample of signal Fu is asked for In peak point after leaf transformation phase value, so as to obtain K phase value φ (k), k=1,2 ... ..., K；By two neighboring phase Place value subtracts each other φ (k)-φ (k-1) respectively, and k takes 2 herein ... ..., K, so as to obtain K-1 phase difference φ (2), Δ φ (3) ... ..., Δ φ (K)；
Step 3：Time difference t between calculating two neighboring Complex pulse sample of signal and carrying out the starting sample point of Fourier transformation0 (k)-t0(k-1), t0K () represents the k-th starting sample point corresponding time of the Fourier transformation of Complex pulse sample of signal, K takes 2 herein ... ..., K, obtains Δ t (2), Δ t (3) ... ..., Δ t (K)；
Step 4：Iterative calculation below equation K-1 times successively, obtain f (2), f (3) ..., f (K)：
Φ (n)=2 π f (n) Δ t (n+1)-Δ φ (n+1)；
m ( n ) = [ Φ ( n ) 2 π ] ;
Φ ' (n)=2m (n) π+Δ φ (n+1)；
f ( n + 1 ) = Φ ′ ( n ) 2 π Δ t ( n + 1 ) ;
Wherein, n take 1 successively, 2,3 ..., K-1；F (1)=fsp, [] represents and asks for integer according to the principle for rounding up；
F (K) is the multiple-pulse frequency precise measurements of the real number list carrier frequency pulse string signal sample.
2. a kind of multiple-pulse Precision Method of Freuqency Measurement of single carrier frequency pulse signal according to claim 1, its feature exists In step 1 is further included：
Step 11：The phase value of whole sampled points is asked in the 1st Complex pulse sample of signal, by two neighboring sampled point Phase value subtracts each other and obtainsWherein i=1,2 ... ..., L, it is two neighboring in the 1st Complex pulse sample of signal of expression to adopt The order label of the phase difference of sampling point, L represents the sum of the two neighboring sampled point phase difference of Complex pulse sample of signal；
Step 12：According to formulaL phase weighting coefficients, wherein i=1 are asked for, 2 ... ..., L；Computing formulaCan obtain the 1st pulse of Complex pulse sample of signal Frequency measurement fsp, wherein tsIt is the sampling interval duration of Complex pulse sample of signal, i=1,2 ... ..., L.
3. a kind of multiple-pulse frequency device for accurately measuring of single carrier frequency pulse signal, it is characterised in that including：
Frequency iteration calculation of initial value unit, the real number list carrier frequency pulse string signal sample for that will receive carries out positive alternation Change, obtain K Complex pulse sample of signal；The 1st Complex pulse sample of signal is processed using Kay algorithms, obtains frequency Rate iteration initial value fsp；Real number list carrier frequency pulse string signal sample includes K real number list carrier frequency pulse signal；K is whole more than 1 Number；
Fourier transformation peak point phase difference calculating unit, for carrying out Fourier's change to each Complex pulse sample of signal Change, ask for the phase value of the peak point after each Complex pulse sample of signal Fourier transformation, so as to obtain K phase value φ (k), k=1,2 ... ..., K；Two neighboring phase value is subtracted each other into φ (k)-φ (k-1) respectively, k takes 2 herein ... ..., K, so that Obtain K-1 phase difference φ (2), Δ φ (3) ... ..., Δ φ (K)；
Fourier transformation starting sample point time difference calculating unit, Fu is carried out for calculating two neighboring Complex pulse sample of signal In leaf transformation starting sample point between time difference t0(k)-t0(k-1), t0K-th Complex pulse sample of signal of (k) expression The starting sample point corresponding time of Fourier transformation, herein k take 2 ... ..., K obtains Δ t (2), Δ t (3) ... ..., Δ t (K)；
Multiple-pulse frequency precise measurements iterate to calculate unit, for iterating to calculate below equation K-1 times successively, obtain f (2), f (3)、…、f(K)：
Φ (n)=2 π f (n) Δ t (n+1)-Δ φ (n+1)；
m ( n ) = [ Φ ( n ) 2 π ] ;
Φ ' (n)=2m (n) π+Δ φ (n+1)；
f ( n + 1 ) = Φ ′ ( n ) 2 π Δ t ( n + 1 ) ;
Wherein, n take 1 successively, 2,3 ..., K-1；F (1)=fsp, [] represents and asks for integer according to the principle for rounding up；
F (K) is the multiple-pulse frequency precise measurements of the real number list carrier frequency pulse string signal sample.
4. a kind of multiple-pulse frequency device for accurately measuring of single carrier frequency pulse signal according to claim 3, its feature exists In frequency iteration calculation of initial value unit is further included：
Signal phase difference computing unit, the phase value for asking for whole sampled points in the 1st Complex pulse sample of signal will The phase value of two neighboring sampled point subtracts each other and obtainsWherein i=1,2 ... ..., L, represent the 1st Complex pulse signal The order label of the phase difference of two neighboring sampled point in sample, L represents the two neighboring sampled point phase of Complex pulse sample of signal The sum of potential difference；
Phase weighting unit, for according to formulaL phase weighting coefficients are asked for, its Middle i=1,2 ... ..., L；Computing formulaI.e. available 1st Complex pulse sample of signal Pulse frequency measurement fsp, wherein tsIt is the sampling interval duration of Complex pulse sample of signal, i=1,2 ... ..., L.
CN201710278139.6A 2017-04-25 2017-04-25 A kind of the multiple-pulse Precision Method of Freuqency Measurement and device of single carrier frequency pulse signal Pending CN106872773A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710278139.6A CN106872773A (en) 2017-04-25 2017-04-25 A kind of the multiple-pulse Precision Method of Freuqency Measurement and device of single carrier frequency pulse signal

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710278139.6A CN106872773A (en) 2017-04-25 2017-04-25 A kind of the multiple-pulse Precision Method of Freuqency Measurement and device of single carrier frequency pulse signal

Publications (1)

Publication Number Publication Date
CN106872773A true CN106872773A (en) 2017-06-20

Family

ID=59161688

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710278139.6A Pending CN106872773A (en) 2017-04-25 2017-04-25 A kind of the multiple-pulse Precision Method of Freuqency Measurement and device of single carrier frequency pulse signal

Country Status (1)

Country Link
CN (1) CN106872773A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109375128A (en) * 2018-10-15 2019-02-22 合肥本源量子计算科技有限责任公司 A kind of carrier frequency pulse signal frequency calibration method
CN111896922A (en) * 2020-08-06 2020-11-06 南京理工雷鹰电子科技有限公司 Method for measuring emission frequency of magnetron
CN112162152A (en) * 2020-08-31 2021-01-01 南京亿杰明信息技术有限公司 Sine wave coherent pulse train signal frequency estimation method based on phase straight line fitting
CN113824494A (en) * 2021-10-29 2021-12-21 北京邮电大学 Linear light sampling-based peak value extraction method and system

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5706202A (en) * 1995-03-08 1998-01-06 Anritsu Corporation Frequency spectrum analyzing apparatus and transmitter characteristics measuring apparatus using the same
EP1239455A2 (en) * 2001-03-09 2002-09-11 Alcatel Method and system for implementing a Fourier transformation which is adapted to the transfer function of human sensory organs, and systems for noise reduction and speech recognition based thereon
CN101057150A (en) * 2004-09-13 2007-10-17 菲尔特罗尼克公开有限公司 Method and apparatus for measuring the frequency of a received signal
CN102680826A (en) * 2012-05-18 2012-09-19 中国电子科技集团公司第四十一研究所 Method for realizing test of embedded local-oscillator frequency converter by utilizing vector network analyzer
JP2015036635A (en) * 2013-08-12 2015-02-23 株式会社近計システム Frequency measuring method and frequency measuring device using fourier analysis
CN104407212A (en) * 2014-10-23 2015-03-11 电子科技大学 Real-time and high-precision coherent pulse signal frequency measurement method
CN104502701A (en) * 2014-12-10 2015-04-08 广东电网有限责任公司电力科学研究院 Method and system for detecting electric power signal frequency based on phase modulation
CN104635045A (en) * 2015-02-05 2015-05-20 广东电网有限责任公司电力科学研究院 Power signal frequency detection method and system based on phase modulation
CN204989407U (en) * 2015-06-30 2016-01-20 绵阳市维博电子有限责任公司 Phase sensitivity track signal detecting system
CN105445547A (en) * 2015-07-09 2016-03-30 深圳市科润宝实业有限公司 Phase detection method and system of sinusoidal signals
CN105785123A (en) * 2016-03-22 2016-07-20 电子科技大学 Radar signal frequency calculation method based on apFFT phase difference

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5706202A (en) * 1995-03-08 1998-01-06 Anritsu Corporation Frequency spectrum analyzing apparatus and transmitter characteristics measuring apparatus using the same
EP1239455A2 (en) * 2001-03-09 2002-09-11 Alcatel Method and system for implementing a Fourier transformation which is adapted to the transfer function of human sensory organs, and systems for noise reduction and speech recognition based thereon
CN101057150A (en) * 2004-09-13 2007-10-17 菲尔特罗尼克公开有限公司 Method and apparatus for measuring the frequency of a received signal
CN102680826A (en) * 2012-05-18 2012-09-19 中国电子科技集团公司第四十一研究所 Method for realizing test of embedded local-oscillator frequency converter by utilizing vector network analyzer
JP2015036635A (en) * 2013-08-12 2015-02-23 株式会社近計システム Frequency measuring method and frequency measuring device using fourier analysis
CN104407212A (en) * 2014-10-23 2015-03-11 电子科技大学 Real-time and high-precision coherent pulse signal frequency measurement method
CN104502701A (en) * 2014-12-10 2015-04-08 广东电网有限责任公司电力科学研究院 Method and system for detecting electric power signal frequency based on phase modulation
CN104635045A (en) * 2015-02-05 2015-05-20 广东电网有限责任公司电力科学研究院 Power signal frequency detection method and system based on phase modulation
CN204989407U (en) * 2015-06-30 2016-01-20 绵阳市维博电子有限责任公司 Phase sensitivity track signal detecting system
CN105445547A (en) * 2015-07-09 2016-03-30 深圳市科润宝实业有限公司 Phase detection method and system of sinusoidal signals
CN105785123A (en) * 2016-03-22 2016-07-20 电子科技大学 Radar signal frequency calculation method based on apFFT phase difference

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
STEVEN KAY: "A Fast and accurate Single Frequency Estimator", 《IEEE TRANSACTIONS ON ACOUSTICS. SPEECH. AND SIGNAL PROCESSING》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109375128A (en) * 2018-10-15 2019-02-22 合肥本源量子计算科技有限责任公司 A kind of carrier frequency pulse signal frequency calibration method
CN109375128B (en) * 2018-10-15 2020-09-15 合肥本源量子计算科技有限责任公司 Carrier frequency pulse signal frequency calibration method
CN111896922A (en) * 2020-08-06 2020-11-06 南京理工雷鹰电子科技有限公司 Method for measuring emission frequency of magnetron
CN111896922B (en) * 2020-08-06 2023-08-18 南京雷鹰电子科技有限公司 Method for measuring emission frequency of magnetron
CN112162152A (en) * 2020-08-31 2021-01-01 南京亿杰明信息技术有限公司 Sine wave coherent pulse train signal frequency estimation method based on phase straight line fitting
CN112162152B (en) * 2020-08-31 2024-01-26 南京亿杰明信息技术有限公司 Sine wave coherent pulse train signal frequency estimation method based on phase straight line fitting
CN113824494A (en) * 2021-10-29 2021-12-21 北京邮电大学 Linear light sampling-based peak value extraction method and system
CN113824494B (en) * 2021-10-29 2022-07-01 北京邮电大学 Linear light sampling-based peak value extraction method and system

Similar Documents

Publication Publication Date Title
Eichstädt et al. Deconvolution filters for the analysis of dynamic measurement processes: a tutorial
CN106872773A (en) A kind of the multiple-pulse Precision Method of Freuqency Measurement and device of single carrier frequency pulse signal
Chen et al. Digital fractional order Savitzky-Golay differentiator
CN107102255B (en) Single ADC acquisition channel dynamic characteristic test method
CN107085140B (en) Nonequilibrium system frequency estimating methods based on improved SmartDFT algorithm
CN101403774B (en) Harmonic wave analysis method based on non-synchronous sampling
CN108414833B (en) Accurate estimation method for signal component frequency
CN104142425B (en) Phase matching method for sinusoidal signal frequency estimation
CN104218973A (en) Frequency hopping signal parameter estimation method based on Myriad filtering
CN105024951B (en) A kind of power delay spectrum PDP methods of estimation and device
CN105137373B (en) A kind of denoising method of exponential signal
CN111159888A (en) Covariance matrix sparse iteration time delay estimation method based on cross-correlation function
CN110705041A (en) Linear structure working modal parameter identification method based on EASI
CN107766293B (en) Signal spectrum analysis method and system when part of sampling data is missing regularly
CN105300386B (en) A kind of frequency domain weighting of X-ray pulsar photon sequence compares phase method
Zygarlicki et al. Prony’s method with reduced sampling-numerical aspects
CN110320400B (en) Voltage flicker envelope parameter extraction method for quasi-synchronous sampling and improved energy operator
CN108459992A (en) A kind of power distribution network synchronous phasor measuring method based on Prony algorithms
Matusiak et al. Noniterative method for frequency estimation based on interpolated DFT with low-order harmonics elimination
CN112014811B (en) Fine estimation method for radar carrier frequency
CN108020761A (en) A kind of Denoising of Partial Discharge
CN103441975B (en) A kind of Coded Signals parameter estimation method based on power spectrum
Pálfi et al. Efficient execution of ADC test with sine fitting with verification of excitation signal parameter settings
CN113155384A (en) Sensor arrangement method for reducing uncertainty of structural damping ratio identification
Goos et al. Continuous time frequency domain LPV state space identification via periodic time-varying input-output modeling

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20170620