CN109782219A - A kind of signal accumulation survey phase method based on correlation - Google Patents
A kind of signal accumulation survey phase method based on correlation Download PDFInfo
- Publication number
- CN109782219A CN109782219A CN201811596460.XA CN201811596460A CN109782219A CN 109782219 A CN109782219 A CN 109782219A CN 201811596460 A CN201811596460 A CN 201811596460A CN 109782219 A CN109782219 A CN 109782219A
- Authority
- CN
- China
- Prior art keywords
- signal
- phase
- correlation
- measured
- measured signal
- 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
Links
Abstract
The present invention provides a kind of, and the signal accumulation based on correlation surveys phase method, is tested phase signals to two-way and is tentatively intercepted according to signal length complete cycle, measured signal A is divided into K sections of new signal, and K sections of new signal head and the tail are realigned and added up；Measured signal B is reprocessed, the signal A and signal B after adding up obtain the cross-correlation function of measured signal A and B both as two overall signals, and to the overall signal and reference signal progress discrete sampling after adding up；It obtains surveying phase result using the two cross-correlation function.The present invention provides signal phase measurement result from the angle of time domain, and calculation amount is smaller, and it is possible to prevente effectively from the problem of spectral leakage；In terms of noise suppressed, the present invention can effectively inhibit noise, to keep measurement result more accurate by way of Signal averaging under the premise of guarantee signal enough energy.
Description
Technical field
The present invention relates to a kind of survey phase methods.
Background technique
Phase measurement all has a wide range of applications in satellite navigation and signal analysis field, and the precision of phase is surveyed
Amount is always the hot research problem of academia and engineering circles.Currently, being mentioned to the research of phase detection techniques gradual perfection
Many Phase Detecting Devices are gone out, the measuring instrument of the various M&S of domestic and international development and production has all been pushed to city
, very important effect has been played in applications such as aerospace detection, production and processing, engineering technology, geodesic surveys.Modern phase
The development of position measuring technique mainly experienced following several stages:
First stage is method that early stage generally uses and poor, method of comparison and balancing method etc., its main feature is that measurement method is former
Reason is easily understood, but measurement accuracy is relatively low；
Second stage is to use various dedicated integrated circuits, MCU etc. to constitute phase difference detection system, is further mentioned
High phase difference detection precision；
Phase III is using computer technology, and emphasis is to make the adaptation of phase measuring equipment to the design of detection program
Property is stronger, and error is smaller.At the same time, the research of various phase detection algorithms is increasingly paid attention to.
Currently, the research hotspot of phase measurement mainly with it is following two aspect based on: first, it is assumed that phase information is
A part of sinusoidal (or cosine) signal, it is present in voltage or current signal, therefore how to reject voltage or electricity
Influence of the frequency variation of signal to phase measurement accuracy is flowed just as the critical issue in phase detection techniques research；Secondly,
To measure phase, not only need to guarantee that the frequency of two-way measured signal is essentially equal, it is often more important that exclude due to two-way
The influences of the factors to measurement result such as the amplitude of measured signal is inconsistent.
Summary of the invention
For overcome the deficiencies in the prior art, the present invention provides a kind of signal accumulation survey phase method based on correlation,
Under the premise of two paths of signals frequency is consistent, solve due to the inconsistent brought survey phase error interference of the amplitude of two-way measured signal
Problem (core manifestation is that there are white Gaussian noises for signal).
The technical solution adopted by the present invention to solve the technical problems the following steps are included:
(1) phase signals are tested to two-way tentatively to be intercepted according to signal length complete cycle；
(2) measured signal A starting point and modelled signal accumulation length are selected, measured signal A is divided into K sections of new letter
Number, K sections of new signal head and the tail are realigned and added up；Repeat this step process measured signal B；
(3) will it is cumulative after signal A and signal B as two overall signals, and to after cumulative overall signal and
Reference signal carries out discrete sampling, obtains the cross-correlation function of measured signal A and B both；It is obtained using the two cross-correlation function
To survey phase result.
The beneficial effects of the present invention are: reducing noise bring in such a way that signal accumulation is superimposed surveys phase error, and mention
Gone out by sampled point at correlation peak deviate in the way of calculate phase, compared to previous conventional method, the present invention is not necessarily to
Fourier transformation is carried out, provides signal phase measurement result from the angle of time domain, calculation amount is smaller, and it is possible to prevente effectively from
The problem of spectral leakage；In terms of noise suppressed, the present invention is by way of Signal averaging, before guaranteeing the enough energy of signal
Noise can effectively be inhibited by putting, to keep measurement result more accurate.
Detailed description of the invention
Fig. 1 is the basic flow chart that the present invention is accumulated by signal cutout and correlation is utilized to survey phase algorithm；
Fig. 2 is two-way measured signal pretreatment process figure；
Fig. 3 is the process flow diagram that phase signals are truncated and add up；
Fig. 4 is to obtain the step flow chart surveyed and be mutually worth according to two paths of signals correlation；
Fig. 5 is the effect contrast figure that phase result and conventional method are surveyed after improving in the case where signal-to-noise ratio is 0dB；
Fig. 6 is the effect contrast figure that phase result and conventional method are surveyed after improving in the case where signal-to-noise ratio is -10dB；
Fig. 7 is the effect contrast figure that phase result and conventional method are surveyed after improving in the case where signal-to-noise ratio is -20dB；
Fig. 8 is the effect contrast figure that phase result and conventional method are surveyed after improving in the case where signal-to-noise ratio is -30dB；
Fig. 9 is the effect contrast figure that phase result and conventional method are surveyed after improving in the case where signal-to-noise ratio is -35dB.
Specific embodiment
Present invention will be further explained below with reference to the attached drawings and examples, and the present invention includes but are not limited to following implementations
Example.
The defect that the present invention is influenced for measurement result in conventional phase measurement method vulnerable to stability of waveform, proposes one kind
Signal accumulation based on correlation surveys phase method, is modified from algorithm level to conventional measurement phase method, the essence of measurement result
Degree and stability are all significantly enhanced, and have very strong application value.
Specific implementation step is as follows:
Step 1: tentatively intercepting measured signal A, and concrete mode is the integer section week choosing coincidence measurement and requiring
Measured signal in phase, and rest part is cast out.Identical operation is carried out to measured signal B simultaneously, the difference is that due to
The measured signal A for needing to add up to measured signal A, therefore being typically selected to have in more number of cycles, and measured signal B
It can choose the less period.Significantly, since needing to carry out digitized processing to signal, therefore in order to guarantee to cut
The number of winning the confidence is number of cycles, it is necessary to set sample frequency to the integral multiple of signal frequency.Two-way measured signal is passed through later
Bandpass filter.Since handled signal is usually single frequency sinusoidal signal, which should have narrow-band characteristic.Filter
Noise can effectively be inhibited while wave, in addition, should use two paths of signals the filter of same structure, purpose exists
In elimination by filter bring phase offset.Measured signal A is indicated are as follows:
xr(t)=cos (2 π f0t+φ0)
Wherein f0Indicate signal frequency, φ0It indicates first phase, then can be indicated after bandpass filter are as follows:
x'r(t)=cos (2 π f0t+φ0+Δφ)
Wherein Δ φ indicates the additional phase delay after bandpass filter.Similarly, measured signal B is indicated are as follows:
x0(t)=cos (2 π f0t)
After identical bandpass filter, measured signal B can be indicated are as follows:
x'0(t)=cos (2 π f0t+Δφ)
Step 2: measured signal A is truncated and is added up.Truncation start-stop point is set according to noisy signal frequency, by cumulative effect
It answers it is found that choose complete cycle most appropriate, sets 100 for accumulation period in the present invention.By the signal alignment after truncation and tire out
Add, obtains new cumulative signal.
Due to that can be indicated after measured signal A filtering are as follows:
x'r(t)=cos (2 π f0t+Δφ+φ0)
Understand from the angle of random process, it is believed that the signal is a wide stationary random process, and noise is being propagated through
Meet gaussian random distribution in journey, measured signal A can be indicated again are as follows:
x'r(t)=cos (2 π f0t+Δφ+φ0)+N(t)
Wherein N (t) indicates Gaussian random process, and N (t)~(0, σ2)。
Then in t1,T+t1,…,KT+t1Moment, signal meet following distribution:
x′r(t1)~N (cos (2 π f0t1+Δφ+φ0), σ2)
x'r(T+t1)~N (cos (2 π f0t1+Δφ+φ0), σ2)
x'r((K-1)T+t1)~N (cos (2 π f0t1+Δφ+φ0), σ2)
It is hereby understood that signal again cumulative after truncation meets:
As can be seen that adding up with signal, noise is opposite to be reduced, this is also the original that signal is intercepted and added up
Reason.It is worth noting that, generally using add up complete cycle in real process.Then in the same way to measured signal
B is handled, and accumulation period can opposite reduction.
Step 3: same treatment in step 2 is carried out to signal B, and calculates two paths of signals relevance function.It enables:
Measured signal B uses x' after the processing of step 20To indicate.
In order to calculate the phase difference of the two, its correlation function can be calculated, for either signal x'c(t) and another delay
For the signal x' of τ0(t+ τ), its cross-correlation function in limited time interval T are as follows:
The value range of τ is-T to T in formula.If τ is variable, then correlation functionIt is exactly the function of delay, τ, by its property
Matter is it is found that work as x'c(t)、x'0When (t+ τ) is identical,With maximum value.
To x'c(t)、x'0(t+ τ) carries out sliding-model control, with fsSample frequency sampled, the place after being sampled
Manage result x'c(n) and x'0(n), wherein n=0,1,2 ... N, and obtain correlation function:
Remember that the corresponding sampled point of maximum related value is lmax, then available two-way sinusoidal signal phase difference are as follows:
ΔφIt surveys=(lmax-N)·2πf0/fs
The available signal accumulation based on correlation surveys phase result through the above steps.The processing method mainly includes two
A part, first is to carry out truncation accumulation process to signal, which plays good inhibiting effect to noise；Second is to adopt
With the survey phase method based on signal correlation, which is only handled, algorithm in the time domain without carrying out Fourier transformation
It is accurate and quick, there is very strong application value.
The above examples are only used to illustrate the technical scheme of the present invention, rather than its limitations；Although with reference to the foregoing embodiments
Invention is explained in detail, those skilled in the art should understand that: it still can be to aforementioned implementation institute
The technical solution of record is modified or equivalent replacement of some of the technical features；And these are modified or replaceed,
It does not separate the essence of the corresponding technical solution, and the present invention implements the spirit and scope of technical solution.
Claims (1)
1. a kind of signal accumulation based on correlation surveys phase method, it is characterised in that include the following steps:
(1) phase signals are tested to two-way tentatively to be intercepted according to signal length complete cycle；
(2) measured signal A starting point and modelled signal accumulation length are selected, measured signal A is divided into K sections of new signal, by K
Section new signal head and the tail are realigned and are added up；Repeat this step process measured signal B；
(3) will it is cumulative after signal A and signal B as two overall signals, and to the overall signal and benchmark after cumulative
Signal carries out discrete sampling, obtains the cross-correlation function of measured signal A and B both；It is surveyed using the two cross-correlation function
Phase result.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811596460.XA CN109782219A (en) | 2018-12-26 | 2018-12-26 | A kind of signal accumulation survey phase method based on correlation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811596460.XA CN109782219A (en) | 2018-12-26 | 2018-12-26 | A kind of signal accumulation survey phase method based on correlation |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109782219A true CN109782219A (en) | 2019-05-21 |
Family
ID=66497674
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811596460.XA Pending CN109782219A (en) | 2018-12-26 | 2018-12-26 | A kind of signal accumulation survey phase method based on correlation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109782219A (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1774689A (en) * | 2003-04-16 | 2006-05-17 | 索尼爱立信移动通讯股份有限公司 | Direct digital frequency synthesizer for cellular wireless communication systems based on fast frequency-hopped spread spectrum technology |
| CN104280748A (en) * | 2013-07-12 | 2015-01-14 | 曲卫 | Pre-accumulation post-correlation batch processing and differential coherence navigation signal capture method |
| CN104898172A (en) * | 2015-05-19 | 2015-09-09 | 吉林大学 | Cross-correlation-based nuclear magnetic resonance full wave signal noise filtering method |
| CN106842263A (en) * | 2017-01-21 | 2017-06-13 | 中国人民解放军信息工程大学 | Based on the satellite navigation method that navigation signal space time processing is combined with vector tracking |
| WO2018225250A1 (en) * | 2017-06-09 | 2018-12-13 | 三菱電機株式会社 | Radar device |
-
2018
- 2018-12-26 CN CN201811596460.XA patent/CN109782219A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1774689A (en) * | 2003-04-16 | 2006-05-17 | 索尼爱立信移动通讯股份有限公司 | Direct digital frequency synthesizer for cellular wireless communication systems based on fast frequency-hopped spread spectrum technology |
| CN104280748A (en) * | 2013-07-12 | 2015-01-14 | 曲卫 | Pre-accumulation post-correlation batch processing and differential coherence navigation signal capture method |
| CN104898172A (en) * | 2015-05-19 | 2015-09-09 | 吉林大学 | Cross-correlation-based nuclear magnetic resonance full wave signal noise filtering method |
| CN106842263A (en) * | 2017-01-21 | 2017-06-13 | 中国人民解放军信息工程大学 | Based on the satellite navigation method that navigation signal space time processing is combined with vector tracking |
| WO2018225250A1 (en) * | 2017-06-09 | 2018-12-13 | 三菱電機株式会社 | Radar device |
Non-Patent Citations (2)
| Title |
|---|
| 史芳芳: "基于DSP的相位差测量系统的研究", 《中国优秀博硕士学位论文全文数据库（硕士） 工程科技Ⅱ辑》 * |
| 徐成发 等: "基于互相关的快速角度估计算法", 《电子与信息学报》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Brennan et al. | On the relationship between time and frequency domain methods in time delay estimation for leak detection in water distribution pipes | |
| Zhang et al. | Analysis of white noise on power frequency estimation by DFT-based frequency shifting and filtering algorithm | |
| CN108548957B (en) | Dual-spectrum analysis method based on combination of cyclic modulation spectrum and piecewise cross correlation | |
| CN108318737A (en) | A kind of high accuracy number phase detecting method | |
| CN107305223B (en) | Improved phase difference frequency estimation method | |
| Sun et al. | Digital recursive demodulator based on Kalman filter | |
| CN114460527B (en) | Correlation degree continuation Hilbert phase-shifting electronic transformer calibrator source tracing method and system | |
| Ma et al. | Harmonic and interharmonic analysis of mixed dense frequency signals | |
| Chen et al. | Robust precise time difference estimation based on digital zero-crossing detection algorithm | |
| CN109030944B (en) | Real sinusoidal signal phase difference estimation method based on frequency shift and DFT | |
| CN109444539B (en) | Synchronous phasor measurement method based on Clark transformation | |
| CN113155210A (en) | Electromagnetic flow measurement system signal processing method based on improved Kalman filtering | |
| KR20080037136A (en) | The method of power frequency estimation using the difference between the gain and cosine and sine filter | |
| CN109782219A (en) | A kind of signal accumulation survey phase method based on correlation | |
| CN104677486A (en) | Aero-engine vibration signal phase measurement method based on revolving speed pulse reconstruction | |
| Huang et al. | A novel high-accuracy digitalized measuring phase method | |
| CN106093905B (en) | Precision instrumentation radar frequency deviation classification frequency measuring method based on apFFT phase differences | |
| CN110133738B (en) | IpDFT-based frequency estimation method for free induction attenuation signal of proton magnetometer | |
| CN111551785B (en) | Frequency and harmonic detection method based on unscented Kalman filter | |
| Chen et al. | A phase difference measurement method based on the extended kalman filter for Coriolis mass flowmeters | |
| Yang et al. | GPS multipath signal estimation based on all phase DFT | |
| CN112595889B (en) | under-Nyquist sampling and parameter measuring method for non-ideal multi-exponential decay sinusoidal signal | |
| Zaiqun | Study on detection of weak signal by phase lock-in amplification | |
| Xu et al. | An iterative three-point interpolation algorithm for one/multiple damped real-valued sinusoids | |
| Hu et al. | Research and implementation of amplitude ratio calculating system for coriolis mass flowmeter |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190521 |