CN108710026A - Frequency stability measurement method based on high-precision phase frequency analysis and system - Google Patents
Frequency stability measurement method based on high-precision phase frequency analysis and system Download PDFInfo
- Publication number
- CN108710026A CN108710026A CN201810614999.7A CN201810614999A CN108710026A CN 108710026 A CN108710026 A CN 108710026A CN 201810614999 A CN201810614999 A CN 201810614999A CN 108710026 A CN108710026 A CN 108710026A
- Authority
- CN
- China
- Prior art keywords
- frequency
- digital
- phase
- signal
- precision
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/02—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage
Abstract
The invention discloses a kind of frequency stability measurement method and system based on high-precision phase frequency analysis, the measurement method uses full digital implementation, the analog frequency signal that measured frequency source and derived reference signal are directly treated using multi-channel high-speed ADC executes sampling and digitized processing, double mixing frequency measurements are realized followed by digital quadrature down conversion technology, and then pass through high-precision phase frequency analysis, phase difference calculates and frequency differential is calculated high-precision frequency differential and measures sequence, calculate overlapping Allan variance, to obtain the frequency stability measurement of frequency source to be measured.The measurement method proposed avoids the shortcomings of simulation and the design difficulty of Semi-digital implementation are big, the development cycle is long；Using high-precision phase frequency parser relative to current digital phase discriminator measurement method, reduce the requirement of the quantization resolution to measuring system ADC devices, hardware cost is reduced, and provides and uses higher rate ADC devices to improve the ability that system can measure frequency range.
Description
Technical field
The present invention relates to a kind of frequency stability measurement method and system based on high-precision phase frequency analysis, belong to electronics
Field of measuring technique.
Background technology
High stable frequency signal source, such as high stable crystal oscillator, atomic clock, be widely used in aerospace,
The numerous areas such as navigator fix, national defense and military, communication.Frequency stability is to evaluate the important indicator of frequency signal source performance, instead
The output frequency value for having reflected frequency source in a period of time is maintained at the ability of a numerical value, is quantified usually using Allan variance
Description.Currently, the second stability of crystal oscillator is up to 10-12, the second stability of atomic clock is up to 10-15, therefore, research
High-precision frequency stability measurement method is highly important.
The highest frequency stability measurement method of measurement accuracy is double mixing frequency measurement methods at present.Realize that double mixing frequency measurement methods can
Using simulated implementation scheme, Semi-digital implementation and full digital implementation.Simulation and Semi-digital implementation pass through mould
Quasi- circuit realizes double mixing, the former usually directly carries out phase measurement using time-interval counter to simulation difference frequency signal；Afterwards
Low speed analog-digital converter (ADC) obtains digital difference frequency signal to simulation difference frequency signal sampling and digitlization in person's use, and then makes
Realize that phase frequency measures with digital signal processing method.Full digital implementation is using high-speed ADC directly to measured source and reference
The analog frequency signal in source executes sampling and digitized processing, generates digital difference frequency signal by digital down-conversion technology, then
Phase frequency is completed by digital signal processing method to measure, most common method is that digital phase discriminator surveys phase method.Currently, most
Measurement method depends heavily on the performance of analog device and circuit using simulation or Semi-digital implementation.When waiting for measured frequency
Source stability is high (such as atomic clock), it is especially desirable to while when measuring the frequency stability in multiple sources to be measured, it is necessary to it makes multiple high
Performance and analog down module as consistent as possible, research and development difficulty is big, the period is long.Full digital implementation passes through multichannel
Analog signal to be measured is synchronously converted to digital signal by ADC, then realizes double mixing frequency measurement methods by digital down-conversion technology,
Not only avoid the difficulty of design high performance analog circuits, additionally it is possible to rapidly expand out stable and completely the same multiple of performance
Measurement channel.However, the measurement result of most common digital phase discriminator survey phase method is influenced seriously by the quantization error of ADC at present,
It is required that the quantization resolution for the ADC chips that measuring system hardware uses is not less than 12, system is caused to be difficult with sample frequency
Higher ADC chips, limit system can measured frequency range, be unable to measure frequency signal source output very high frequency band signal
(such as 100MHz signals of atomic clock output).In addition, the phase measurement sequence of digital phase discriminator output must also pass through solution
After a variety of data processing methods such as winding, filtering, data fitting, sufficiently high measurement accuracy could be provided.
Invention content
For problems of the prior art and disadvantage, the present invention provides a kind of frequency based on high-precision phase frequency analysis
Stable measurement method and system.The measurement method uses full digital implementation, directly right using multi-channel high-speed ADC
The analog frequency signal of frequency source and derived reference signal to be measured executes sampling and digitized processing, followed by becoming under digital quadrature
Frequency technology realizes double mixing frequency measurements, and then height is calculated by high-precision phase frequency analysis, phase difference calculating and frequency differential
The frequency differential of precision measures sequence, overlapping Allan variance is calculated, to obtain the frequency stability measurement of frequency source to be measured.Institute
The measurement method of proposition avoids the shortcomings of simulation and the design difficulty of Semi-digital implementation are big, the development cycle is long；It uses
FFT+CZT joint spectrums interpolator arithmetic is reduced relative to current digital phase discriminator measurement method to measuring system ADC devices
Quantization resolution requirement, reduce hardware cost, and provide use higher rate ADC devices with improve system can measure
The ability of frequency range.
The present invention uses following technical scheme to solve above-mentioned technical problem：
On the one hand, the present invention provides a kind of frequency stability measurement method based on high-precision phase frequency analysis, including following
Specific steps：
S1. analog-to-digital conversion is carried out respectively by analog frequency signal to be measured and with reference to analog frequency signal, obtain two railway digitals
Sampled signal；
After the two-way digital sampled signal that S2.S1 is obtained carries out digital quadrature down conversion processing respectively, two railway digitals are obtained
Difference frequency signal；
After the two-way number difference frequency signal that S3.S2 is obtained carries out high-precision phase frequency analyzing processing respectively, two-way frequency is obtained
Phase measurement；
The two-way frequency and phase position measurement result that S4.S3 is obtained carries out phase difference calculating respectively, obtains two-way frequency measurement
As a result sequence；
S5. frequency differential calculating is carried out to the two-way frequency measurement sequence that S4 is obtained, obtains analog frequency letter to be measured
Number and with reference to analog frequency signal frequency differential measurement result sequence；
S6. it carries out overlapping Allan variance to the frequency differential measurement result sequence that S5 is obtained to calculate, obtains simulation frequency to be measured
The opposite frequency stability measurement result for referring to analog frequency signal of rate signal.
As the further technical solution of the present invention, divided by the digital quadrature down conversion unit of two independent synchronizations in S2
It is other that digital quadrature down conversion processing is carried out to two-way digital sampled signal.
As the further technical solution of the present invention, analog frequency signal to be measured in S1 and with reference to analog frequency signal simultaneously
Row, synchronous progress analog-to-digital conversion process.
As the further technical solution of the present invention, the high-precision phase frequency analyzing processing in S3 uses Fast Fourier Transform (FFT)
FFT and CZT joint spectrum interpolator arithmetics：First, the FFT for calculating digital difference frequency signal obtains the low resolution of digital difference frequency signal
Rate frequency spectrum, to obtain coarse frequency estimation result；Then, using coarse frequency estimation result as amplification center, high score is carried out by CZT
The local spectrum of resolution amplifies, and obtains high-precision signal spectrum as a result, to obtain high-precision frequency and phase measurement knot
Fruit.
On the other hand, the present invention also provides a kind of frequency stability measurement system based on high-precision phase frequency analysis, the surveys
Amount system includes frequency source to be measured, reference frequency source, the first and second analog-digital converters, the first and second digital quadrature down conversions
Unit, the first and second high-precision phase frequency analytic units, the first and second phase difference computing units, frequency differential computing unit
And overlapping Allan variance computing unit, wherein frequency source to be measured, the first analog-digital converter, the first digital quadrature down conversion list
Member, the first high-precision phase frequency analytic unit and first phase difference computing unit are sequentially connected, and reference frequency source, the second modulus turn
Parallel operation, the second digital quadrature down conversion unit, the second high-precision phase frequency analytic unit and second phase Difference Calculation unit are successively
Connection, the first and second phase difference computing units are connect with frequency differential computing unit respectively, frequency differential computing unit
It is connected with Chong Die Allan variance computing unit.
As the further technical solution of the present invention, the first and second digital quadrature down conversion units by local frequency and
Accurately adjustable digital controlled oscillator and multiplier realize digital quadrature down conversion to initial phase, and then are reduced by down-sampled technology
The sample frequency of output digit signals, the digital difference frequency signal exported to parallel real-time synchronization.
As the further technical solution of the present invention, the first and second analog-digital converters are parallel, synchronously complete to be measured
Frequency source and the analog frequency signal of reference frequency source output execute sampling and digitized processing.
As the further technical solution of the present invention, the first and second digital quadrature down conversion units are the two of independent synchronization
A digital quadrature down conversion unit.
As the further technical solution of the present invention, the first and second high-precision phase frequency analytic units are all made of in quick Fu
Leaf transformation FFT carries out high-precision phase frequency analyzing processing with CZT joint spectrum interpolator arithmetics, specially：First, digital difference frequency is calculated
The FFT of signal obtains the low resolution frequency spectrum of digital difference frequency signal, to obtain coarse frequency estimation result；Then, with coarse frequency
Estimated result is amplification center, carries out high-resolution local spectrum amplification by CZT, obtains high-precision signal spectrum knot
Fruit, to obtain high-precision frequency and phase measurement.
The present invention has the following technical effects using above technical scheme is compared with the prior art：
1) full digital scheme is used to realize double mixing frequency measurement methods, compared to simulation and Semi-digital implementation, hardware system
It is small to research and develop difficulty, it is easy to accomplish, performance is stablized, and can rapidly expand Measurement channel quantity, maintenance and upgrade is simple；
2) the high-precision phase frequency parser used is calculated by FFT+CZT joint spectrum interpolator arithmetic combination phase differences
Method surveys phase method instead of common digital phase discriminator, significantly reduces quantization resolution of the measuring system to ADC chips
Requirement, to reduce hardware cost；
3) high-precision phase frequency parser allows measuring system to measure needs according to specific, using low quantization resolution, height
The single core multichannel ADC of sample frequency carrys out the frequency measurement range of expansion system, and high quantization resolution ratio, high sample frequency
Multi-core multichannel ADC is difficult to obtain；
4) Allan variance is replaced using overlapping Allan variance, avoided when average time, τ was larger, Allan variance calculating makes
The shortcomings that data points are reduced, and confidence level reduces.
Description of the drawings
Fig. 1 is the solution principle schematic diagram of the frequency stability measurement method based on high-precision phase frequency analysis；
Fig. 2 is the principle schematic for the FFT+CZT joint spectrum interpolator arithmetics that high-precision phase frequency analytic unit uses.
Fig. 3 is that system bottom in embodiment is made an uproar measurement result.
Specific implementation mode
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation describes, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
With reference to the accompanying drawings and detailed description, the measurement method of the present invention is made to be better understood.
The present invention provides a kind of frequency stability measurement methods based on high-precision phase frequency analysis, which is characterized in that packet
Include step in detail below：
S1. analog-to-digital conversion is carried out respectively by analog frequency signal to be measured and with reference to analog frequency signal, obtain two railway digitals
Sampled signal；
After the two-way digital sampled signal that S2.S1 is obtained carries out digital quadrature down conversion processing respectively, two railway digitals are obtained
Difference frequency signal；
After the two-way number difference frequency signal that S3.S2 is obtained carries out high-precision phase frequency analyzing processing respectively, two-way frequency is obtained
Phase measurement；
The two-way frequency and phase position measurement result that S4.S3 is obtained carries out phase difference calculating respectively, obtains two-way frequency measurement
As a result sequence；
S5. frequency differential calculating is carried out to the two-way frequency measurement sequence that S4 is obtained, obtains analog frequency letter to be measured
Number and with reference to analog frequency signal frequency differential measurement result sequence；
S6. it carries out overlapping Allan variance to the frequency differential measurement result sequence that S5 is obtained to calculate, obtains simulation frequency to be measured
The opposite frequency stability measurement result for referring to analog frequency signal of rate signal.
The present invention also provides a kind of frequency stability measurement systems based on high-precision phase frequency analysis, as shown in Figure 1, should
Measuring system includes：Frequency source 100, reference frequency source 110, analog-digital converter 101 and 111, digital quadrature down conversion list to be measured
Member 102 and 112, high-precision phase frequency analytic unit 103 and 113, phase difference computing unit 104 and 114, frequency differential calculate single
Member 120 and overlapping Allan variance computing unit 121.The analog-digital converter 101 and it is 111 parallel, synchronously complete to be measured
The two-way analog frequency signal that source 100 and reference source 110 export executes sampling and digitlization.The digital quadrature down conversion unit
102 and 112, it is realized under digital quadrature by local frequency and the accurate adjustable orthogonal digital controlled oscillator of initial phase and multiplier
Frequency conversion, and then reduce by down-sampled technology the sample frequency of output signal, it is parallel, in real time, synchronously output digital difference frequency letter
Number.The high-precision phase frequency analytic unit 103 and 113 is parallel, real-time, accurate by FFT+CZT joint spectrum interpolator arithmetics
Ground measurement frequency and phase.The phase difference computing unit 104 and 114 utilizes high-precision phase frequency analytic unit 103 and 113
Frequency plot result calculate phase difference, export higher precision and with high number of significant digit frequency measurement sequence.Institute
It states frequency differential computing unit 120 and calculates and export the frequency differential result sequence that two channel frequences measure the correspondence moment of sequence
Row.The overlapping Allan variance unit 121 calculates simultaneously according to the frequency differential result sequence that frequency differential computing unit 120 exports
Output overlapping Allan variance measurement result.
Technical scheme of the present invention is further elaborated below by specific embodiment：
(1) frequency source 100 and reference frequency source 110 output nominal frequency to be measured is f0Analog frequency signal.Wherein,
The analog frequency signal x to be measured that frequency source 100 to be measured exportsmea(t), the reference analog frequency signal x of reference frequency source outputref
(t) analog-digital converter 101 and 111 is respectively enterd, respectively obtains and exports to be measured and reference number sampled signal xmea(n) and xref
(n)。
In this example, the 100MHz analog frequencies signal that the same frequency source exports is divided through power divider one and is
Two, it is regarded as the mould analog frequency signal from frequency source 100 to be measured and reference frequency source 110 respectively, is respectively connected to analog-to-digital conversion
Device 101 and 111.
In this example, measuring system use analog-digital converter ADC have four-way, sample frequency highest 1.2GHz,
Quantization resolution is 8.The analog-digital converter 101 and 111 is respectively two in 4 cores of the chip.Sample frequency fs
=1152MHz, measurable nominal frequency range are f0≤500MHz。
(2) according to the nominal frequency f of input signal0, the numerical control of the digital quadrature down conversion unit 102 and 112 is set
The local frequency f of oscillatorLO=f0-fb, wherein fbIt is the frequency of digital difference frequency signal.
Digital quadrature down conversion unit 102 and 112 real-time parallels according to fLOWith identical initial phase θ0, according to formula
(1) and (2) are to xmea(n) and xref(n) digital quadrature down conversion processing is executed：
WhereinWithIt is x respectivelymea(n) and xref(n) corresponding two-way digital quadrature down conversion output
Complex signal.
Then, digital quadrature down conversion unit 102 and 112 is right respectivelyWithExecute anti-aliasing low pass filtered
Wave processing and D times of down-sampled operation, by the sample frequency of complex digital signal from fsIt is reduced to fbase=fs/ D, wherein D are just whole
Number exports to be measured and reference number difference frequency signalWithWherein nD=nD.
Digital quadrature down conversion can be effectively prevented from the problems such as local oscillator is uneven in simulation quadrature frequency conversion, and can guarantee
The synchronism and consistency of multiple intermodules are easy to extend multichannel measuring system.Down-sampled technology can not lose information
Under the premise of realize data compression, greatly reduce data volume and follow-up signal processing operand.
In this example, digital quadrature down conversion unit is realized by hardware program language in programmable device FPGA chip
On, use the cascade realization structure of multi-stage digital down coversion.The frequency f of digital difference frequency signalb=100Hz, local frequency are set
It is set to fLO=99.999900MHz.After down-sampled, the sample frequency of digital difference frequency signal is fbase=1000Hz.Extracting multiple is
D=1152000 uses total 9 grades anti-aliasing decimation filters and realizes.
(3) the high-precision phase frequency analytic unit 103 and 113, the minimum value τ in the average time τ specified by user0Make
It, will for the time of integrationWithBe divided into a data frame, wherein N=τ per N number of complex sampling point0×
fbase.The serial number of data frame is set as i, wherein i=1,2 ..., I.The measurement total duration T that I is specified by user, by I=T/ τ0It calculates
It obtains.
High-precision phase frequency analytic unit 103 and 113 is respectively to digital difference frequency signalWithEvery number
FFT+CZT joint spectrum interpolator arithmetics are executed according to frame, its principle is as shown in Figure 2, and then the frequency of the frequency source to be measured of the i-th data frame
Rate phase measurementWith the frequency and phase position measurement result of reference frequency source
Specific steps are as follows：
A. assume x (nD) be some digital difference frequency signal i-th of data frame digital signal sequences, to x (nD) L
Complex digital signal sequence executes the plural FFT of L points, obtains x (nD) L point plural number FFT result sequence XsFFT(l)：
Wherein, l indicates the frequency domain serial number of plural FFT result sequence.
And then calculate XFFT(l) the Mo Zhi | at each l;XFFT(l)|：
Wherein, Re { } and Im { } indicates to take real and imaginary parts to plural number respectively.
B. Gen Ju |XFFT(l)|Frequency domain serial number l corresponding to maximum valuemax, calculate the initial position l of CZTstart=lmax-
0.5.And then x (n are calculated as followsD) M point CZT, obtain M point plural number CZT result sequence XsCZT(m)：
Wherein, m indicates the frequency domain serial number of CZT result sequences.And then calculate XCZT(m) the Mo Zhi | of each complex result;XCZT
(m)|。
C. according to lmaxAnd |XCZT(m)|Frequency domain serial number m corresponding to maximum valuemax, it is calculated as follows to obtain digital difference frequency letter
Number x (nD) High Precision Frequency result
Wherein, the frequency resolution Δ f of FFTFFT=fbaseThe frequency resolution Δ f of/L, CZTCZT=fbase/(L×M)。
D. reference number difference frequency signal x (n as the following formulaD) high-precision phase position analysis result
Wherein, the argument of arg { } operator representation calculated complex.
Using high-precision phase frequency analysis (FFT+CZT joint spectrums interpolator arithmetic) into line frequency and phase measurement, Ke Yiyou
Influence of the noise that quantization error in the reduction ADC digital and Digital Signal Processing of effect introduces to measurement result precision.
In this example, minimum average B configuration time τ0=1 second, the sampling number N=1000 points that each data frame packet contains measured total
Duration T=40000 seconds, data frame sum I=40000.FFT points L=2048, the frequency of CZT amplification factors M=1024, FFT
Resolution ax fFFT=fbaseThe frequency resolution Δ f of/L ≈ 0.488Hz, CZTCZT=fbase/(L×M)≈0.477mHz
(4) the phase difference computing unit 104 and 114 is provided according to high-precision phase frequency analytic unit 103 and 113
The frequency and phase position measurement result sequence of frequency source to be measuredWith the frequency and phase position measurement result of reference frequency source
SequenceThe frequency measurement sequence with higher precision and number of significant digit is calculated by formula (8) and (9)With
Wherein, i=1,2 ..., I-1, operatorIndicate downward rounding.
(5) the frequency differential computing unit 120 is by source signal to be measuredWith reference source signal's
Each result correspondence is subtracted each other, and is obtained differential frequency and is measured sequenceWherein i=1,2 ..., I-
1。
Since the signal path that source 100 to be measured and reference source 110 pass through is almost the same, especially digital signal processing path
It is completely the same, after difference processing, noise (the especially system work clock of measuring system itself can be effectively canceled out
Unstability) and influence, obtain precise frequencies stable measurement result.
In this example, differential frequency result sequence F all the way is obtainedb[i], physical significance is the measurement system in embodiment
It makes an uproar at the system bottom of system.
(6) the overlapping Allan variance calculates 121, according to Fb[i]And y[i]=Fb[i]/f0Calculate fractional frequency result sequences y
[i], wherein being f0It is the nominal frequency of signal.And then according to formula
() calculates overlapping Allan variance：
Wherein, τ=s τ0The average time specified for user.
Although Allan variance is the classical measurement in frequency stability measurement calculates Allan when average time τ increases
The fractional frequency data F that variance usesb[i]Points tail off, and the confidence level of result reduces.Using overlapping Allan variance
(Overlapping Allan Variance) can overcome this disadvantage of Allan variance, improve the confidence level of assessment result.
In the present embodiment, average time be τ=1,10,100,1000,3600,7200 seconds.Fig. 3 gives a practical survey
The system bottom of amount is made an uproar result.It can be seen from the figure that under 1 second average time, system bottom makes an uproar only 9.8 × 10-15；1000
System bottom makes an uproar only 3.6 × 10 under second average time-17, the measurement of highly stable frequency signal source (such as hydrogen atomic clock) can be met
It is required that.
Frequency signal source and reference frequency signal to be measured source samples and digitizes through multipath A/D converter in the present invention
Afterwards, it is concurrently handled by multi-channel digital quadrature frequency conversion unit and exports digital difference frequency signal；Digital difference frequency signal is sent into high
Precision phase frequency analytic unit and phase difference computing unit, output high-precision, the frequency measurement sequence with high number of significant digit；Frequently
After rate Difference Calculation unit calculates the differential frequency sequence of frequency source and reference frequency source to be measured, calculated by overlapping Allan variance single
Member calculates overlapping Allan variance.This method has the following advantages：Full digital implementation, relative analog implementation and half
Word implementation reduces the hardware research and development difficulty of measuring system, improves performance and system stability, it is easy to accomplish, it safeguards
And upgrading；Kilosecond grade can be realized in the ADC that high-precision phase frequency analysis and phase difference computational methods only use 8 quantization resolutions
10-17The systematic survey bottom of magnitude is made an uproar, and 12 must not be less than to ADC quantization resolutions by breaching existing full digital implementation
Limitation, so as to use the ADC devices of higher sample frequency to extend to broader frequency measurement range；Using overlapping Allan
Variance replaces Allan variance, the confidence level of measurement result when improving long time-average measurement.
The above, the only specific implementation mode in the present invention, but scope of protection of the present invention is not limited thereto, appoints
What is familiar with the people of the technology within the technical scope disclosed by the invention, it will be appreciated that expects transforms or replaces, and should all cover
Within the scope of the present invention, therefore, the scope of protection of the invention shall be subject to the scope of protection specified in the patent claim.
Claims (9)
1. the frequency stability measurement method based on high-precision phase frequency analysis, which is characterized in that including step in detail below：
S1. analog-to-digital conversion is carried out respectively by analog frequency signal to be measured and with reference to analog frequency signal, obtain two-way digital sample
Signal；
After the two-way digital sampled signal that S2.S1 is obtained carries out digital quadrature down conversion processing respectively, two-way number difference frequency is obtained
Signal；
After the two-way number difference frequency signal that S3.S2 is obtained carries out high-precision phase frequency analyzing processing respectively, two-way frequency plot is obtained
Measurement result；
The two-way frequency and phase position measurement result that S4.S3 is obtained carries out phase difference calculating respectively, obtains two-way frequency measurement
Sequence；
S5. the two-way frequency measurement sequence obtained to S4 carries out frequency differential calculating, obtain analog frequency signal to be measured and
With reference to the frequency differential measurement result sequence of analog frequency signal；
S6. it carries out overlapping Allan variance to the frequency differential measurement result sequence that S5 is obtained to calculate, obtains analog frequency letter to be measured
Number opposite frequency stability measurement result for referring to analog frequency signal.
2. the frequency stability measurement method according to claim 1 based on high-precision phase frequency analysis, which is characterized in that S2
In respectively two-way digital sampled signal under digital quadrature become by the digital quadrature down conversion units of two independent synchronizations
Frequency is handled.
3. the frequency stability measurement method according to claim 1 based on high-precision phase frequency analysis, which is characterized in that S1
In analog frequency signal to be measured and with reference to analog frequency signal parallel, synchronous progress analog-to-digital conversion process.
4. the frequency stability measurement method according to claim 1 based on high-precision phase frequency analysis, which is characterized in that S3
In high-precision phase frequency analyzing processing using Fast Fourier Transform (FFT) FFT and CZT joint spectrum interpolator arithmetics：First, number is calculated
The FFT of word difference frequency signal obtains the low resolution frequency spectrum of digital difference frequency signal, to obtain coarse frequency estimation result；Then, with
Coarse frequency estimation result is amplification center, and high-resolution local spectrum amplification is carried out by CZT, obtains high-precision signal frequency
Spectrum is as a result, to obtain high-precision frequency and phase measurement.
5. the frequency stability measurement system based on high-precision phase frequency analysis, which is characterized in that the measuring system includes waiting for frequency measurement
Rate source, reference frequency source, the first and second analog-digital converters, the first and second digital quadrature down conversion units, first and second
High-precision phase frequency analytic unit, the first and second phase difference computing units, frequency differential computing unit and overlapping Allan side
Poor computing unit, wherein frequency source to be measured, the first analog-digital converter, the first digital quadrature down conversion unit, the first high-precision phase
Frequency analysis unit is sequentially connected with first phase difference computing unit, and reference frequency source, the second analog-digital converter, the second number are just
Down-converter unit, the second high-precision phase frequency analytic unit and second phase Difference Calculation unit is handed over to be sequentially connected, first and second
Phase difference computing unit is connect with frequency differential computing unit respectively, frequency differential computing unit and Chong Die Allan variance meter
Calculate unit connection.
6. the frequency stability measurement system according to claim 5 based on high-precision phase frequency analysis, which is characterized in that the
One and second digital quadrature down conversion unit pass through local frequency and initial phase accurately adjustable digital controlled oscillator and multiplier
It realizes digital quadrature down conversion, and then reduces the sample frequency of output digit signals, parallel real-time synchronization by down-sampled technology
The digital difference frequency signal of ground output.
7. the frequency stability measurement system according to claim 5 based on high-precision phase frequency analysis, which is characterized in that the
One and second analog-digital converter it is parallel, synchronously complete to treat measured frequency source and the analog frequency signal of reference frequency source output is held
Row sampling and digitized processing.
8. the frequency stability measurement system according to claim 5 based on high-precision phase frequency analysis, which is characterized in that the
One and second digital quadrature down conversion unit be independent synchronization two digital quadrature down conversion units.
9. the frequency stability measurement system according to claim 5 based on high-precision phase frequency analysis, which is characterized in that the
One and second high-precision phase frequency analytic unit be all made of Fast Fourier Transform (FFT) FFT and CZT joint spectrum interpolator arithmetics carry out it is high
Precision phase frequency analyzing processing, specially：First, the FFT for calculating digital difference frequency signal obtains the low resolution of digital difference frequency signal
Frequency spectrum, to obtain coarse frequency estimation result；Then, using coarse frequency estimation result as amplification center, high-resolution is carried out by CZT
The local spectrum of rate amplifies, and obtains high-precision signal spectrum as a result, to obtain high-precision frequency and phase measurement.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810614999.7A CN108710026B (en) | 2018-06-14 | 2018-06-14 | Frequency stability measuring method and system based on high-precision phase frequency analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810614999.7A CN108710026B (en) | 2018-06-14 | 2018-06-14 | Frequency stability measuring method and system based on high-precision phase frequency analysis |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108710026A true CN108710026A (en) | 2018-10-26 |
| CN108710026B CN108710026B (en) | 2020-05-05 |
Family
ID=63872509
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810614999.7A Active CN108710026B (en) | 2018-06-14 | 2018-06-14 | Frequency stability measuring method and system based on high-precision phase frequency analysis |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108710026B (en) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109521269A (en) * | 2018-11-09 | 2019-03-26 | 中电科仪器仪表有限公司 | A kind of am signals digitlization frequency measuring method |
| CN110007150A (en) * | 2019-03-28 | 2019-07-12 | 河北远东通信系统工程有限公司 | A kind of linear ratio phase method of Direct Digital Phase Processing |
| CN110007149A (en) * | 2019-03-28 | 2019-07-12 | 河北远东通信系统工程有限公司 | A kind of linear ratio phase method of digitlization phase shift auxiliary |
| CN111367157A (en) * | 2020-04-17 | 2020-07-03 | 中国计量科学研究院 | Multipath phase comparison measuring system and method |
| CN111487500A (en) * | 2020-06-08 | 2020-08-04 | 上海航天测控通信研究所 | System and method for testing high-stability crystal oscillator of satellite-borne atomic clock |
| CN112014640A (en) * | 2020-09-03 | 2020-12-01 | 中国计量科学研究院 | Multi-channel frequency standard comparison test system and working method thereof |
| CN113406386A (en) * | 2021-06-23 | 2021-09-17 | 中国电子科技集团公司第二十九研究所 | Signal frequency accurate estimation method based on digital down-conversion |
| CN114070679A (en) * | 2021-10-25 | 2022-02-18 | 中国电子科技集团公司第二十九研究所 | Frequency-phase characteristic analysis method for intelligent pulse classification |
| CN116506277A (en) * | 2023-04-27 | 2023-07-28 | 南京海汇装备科技有限公司 | Multichannel intelligent test analysis processing system and method |
| CN117491744A (en) * | 2024-01-02 | 2024-02-02 | 石家庄博铭仪器有限公司 | Comprehensive time-frequency parameter measurement method for wide-range arbitrary frequency point |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101021555A (en) * | 2007-03-13 | 2007-08-22 | 熊猫电子集团有限公司 | Frequency marker calibrating system based on GPS frequency standard source |
| KR20110066590A (en) * | 2009-12-11 | 2011-06-17 | 한국표준과학연구원 | Apparatus for measuring frequency stability of oscillator and method for measuring using the same |
| CN102411091A (en) * | 2011-07-27 | 2012-04-11 | 江汉大学 | Device for detecting multi-channel signal stability and detection method |
| CN103293376A (en) * | 2013-05-31 | 2013-09-11 | 江汉大学 | Frequency stability measuring method and device |
| CN103472299A (en) * | 2013-08-14 | 2013-12-25 | 西安电子科技大学 | High-resolution transient frequency stability measuring method |
| CN107247181A (en) * | 2017-04-24 | 2017-10-13 | 西安电子科技大学 | A kind of digitization frequencies stable measurement method of total reponse time |
-
2018
- 2018-06-14 CN CN201810614999.7A patent/CN108710026B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101021555A (en) * | 2007-03-13 | 2007-08-22 | 熊猫电子集团有限公司 | Frequency marker calibrating system based on GPS frequency standard source |
| KR20110066590A (en) * | 2009-12-11 | 2011-06-17 | 한국표준과학연구원 | Apparatus for measuring frequency stability of oscillator and method for measuring using the same |
| CN102411091A (en) * | 2011-07-27 | 2012-04-11 | 江汉大学 | Device for detecting multi-channel signal stability and detection method |
| CN103293376A (en) * | 2013-05-31 | 2013-09-11 | 江汉大学 | Frequency stability measuring method and device |
| CN103472299A (en) * | 2013-08-14 | 2013-12-25 | 西安电子科技大学 | High-resolution transient frequency stability measuring method |
| CN107247181A (en) * | 2017-04-24 | 2017-10-13 | 西安电子科技大学 | A kind of digitization frequencies stable measurement method of total reponse time |
Non-Patent Citations (1)
| Title |
|---|
| 许龙飞等: "一种全面响应时间的频率稳定度测量", 《西安电子科技大学学报（自然科学版）》 * |
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109521269B (en) * | 2018-11-09 | 2021-09-28 | 中电科思仪科技股份有限公司 | Amplitude modulation signal digital frequency measurement method |
| CN109521269A (en) * | 2018-11-09 | 2019-03-26 | 中电科仪器仪表有限公司 | A kind of am signals digitlization frequency measuring method |
| CN110007150A (en) * | 2019-03-28 | 2019-07-12 | 河北远东通信系统工程有限公司 | A kind of linear ratio phase method of Direct Digital Phase Processing |
| CN110007149A (en) * | 2019-03-28 | 2019-07-12 | 河北远东通信系统工程有限公司 | A kind of linear ratio phase method of digitlization phase shift auxiliary |
| CN111367157A (en) * | 2020-04-17 | 2020-07-03 | 中国计量科学研究院 | Multipath phase comparison measuring system and method |
| CN111367157B (en) * | 2020-04-17 | 2022-02-11 | 中国计量科学研究院 | Multipath phase comparison measuring system and method |
| CN111487500A (en) * | 2020-06-08 | 2020-08-04 | 上海航天测控通信研究所 | System and method for testing high-stability crystal oscillator of satellite-borne atomic clock |
| CN111487500B (en) * | 2020-06-08 | 2022-07-05 | 上海航天测控通信研究所 | System and method for testing high-stability crystal oscillator of satellite-borne atomic clock |
| CN112014640A (en) * | 2020-09-03 | 2020-12-01 | 中国计量科学研究院 | Multi-channel frequency standard comparison test system and working method thereof |
| CN113406386A (en) * | 2021-06-23 | 2021-09-17 | 中国电子科技集团公司第二十九研究所 | Signal frequency accurate estimation method based on digital down-conversion |
| CN113406386B (en) * | 2021-06-23 | 2023-04-25 | 中国电子科技集团公司第二十九研究所 | Signal frequency accurate estimation method based on digital down-conversion |
| CN114070679A (en) * | 2021-10-25 | 2022-02-18 | 中国电子科技集团公司第二十九研究所 | Frequency-phase characteristic analysis method for intelligent pulse classification |
| CN114070679B (en) * | 2021-10-25 | 2023-05-23 | 中国电子科技集团公司第二十九研究所 | Pulse intelligent classification-oriented frequency-phase characteristic analysis method |
| CN116506277A (en) * | 2023-04-27 | 2023-07-28 | 南京海汇装备科技有限公司 | Multichannel intelligent test analysis processing system and method |
| CN116506277B (en) * | 2023-04-27 | 2024-01-09 | 南京海汇装备科技有限公司 | Multichannel intelligent test analysis processing system and method |
| CN117491744A (en) * | 2024-01-02 | 2024-02-02 | 石家庄博铭仪器有限公司 | Comprehensive time-frequency parameter measurement method for wide-range arbitrary frequency point |
| CN117491744B (en) * | 2024-01-02 | 2024-03-29 | 石家庄博铭仪器有限公司 | Comprehensive time-frequency parameter measurement method for wide-range arbitrary frequency point |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108710026B (en) | 2020-05-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108710026A (en) | Frequency stability measurement method based on high-precision phase frequency analysis and system | |
| CN104007316B (en) | A kind of High Precision Frequency method under lack sampling speed and measuring instrument thereof | |
| CN105656485B (en) | A kind of multichannel time-interleaved AD C measures calibration method and device | |
| CN104237622B (en) | The method of sampling and wideband voltage/calibration of power device based on software based frequency tracking | |
| CN102857225A (en) | Mismatch error calibration method for multi-channel high-speed parallel alternate sampling system | |
| CN103457603B (en) | A kind of method that ADC dynamic parameter is tested based on average frequency spectrum | |
| CN101915874A (en) | Harmonic wave detection method based on Fourier transformation | |
| CN101713795B (en) | Method of digitalized measuring frequency in restriction of sampling rate | |
| Shaw et al. | A Kalman-filter spectral envelope preprocessor | |
| CN103983849A (en) | Real-time high-accuracy power harmonic analysis method | |
| CN112014640A (en) | Multi-channel frequency standard comparison test system and working method thereof | |
| CN109581018A (en) | Wide-band sampling oscillograph time base error synchroballistic method and device | |
| CN105182069A (en) | High resolution group quantization phase processing method under pilot frequency architecture | |
| CN109813962A (en) | Frequency conversion system group delay measurement method and system based on Hilbert transform | |
| CN103941093B (en) | A kind of two-way DFT symmetrical compensations Method for Phase Difference Measurement and its device | |
| CN105807128B (en) | Multicycle strategy is measured to the method and system of alternating voltage for digital-to-analogue conversion | |
| JPH06242158A (en) | Transmission line constant measuring device and improving method for its measurement precision | |
| RU88157U1 (en) | INFORMATION-MEASURING SYSTEM FOR ELECTRIC ENERGY QUALITY CONTROL | |
| CN114624513A (en) | Anti-harmonic interference phase detection method and device for periodic signal | |
| CN109191543B (en) | AC sampling common-section data generation method | |
| CN106452693A (en) | Clock phase jittering measurement method based on double-frequency-point noise floor energy analysis | |
| CN113422658A (en) | Correction method and system for sampling time sequence asynchronism between channels | |
| CN117491744B (en) | Comprehensive time-frequency parameter measurement method for wide-range arbitrary frequency point | |
| RU2363005C1 (en) | Method of spectral analysis of polyharmonic signals and device to this end | |
| CN110865577B (en) | Digital sampling method and device for alternating current resistance calibration |
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 |