CN106771600B - Adaptive unbiased phase difference estimation method under the conditions of frequency mismatch - Google Patents
Adaptive unbiased phase difference estimation method under the conditions of frequency mismatch Download PDFInfo
- Publication number
- CN106771600B CN106771600B CN201710030146.4A CN201710030146A CN106771600B CN 106771600 B CN106771600 B CN 106771600B CN 201710030146 A CN201710030146 A CN 201710030146A CN 106771600 B CN106771600 B CN 106771600B
- Authority
- CN
- China
- Prior art keywords
- phase difference
- frequency
- signal
- way
- cosine 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R25/00—Arrangements for measuring phase angle between a voltage and a current or between voltages or currents
- G01R25/005—Circuits for comparing several input signals and for indicating the result of this comparison, e.g. equal, different, greater, smaller, or for passing one of the input signals as output signal
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Phase Differences (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
Adaptive unbiased phase difference estimation method the present invention relates to field of signal processing, under the conditions of especially a kind of frequency mismatch.Applicable object of the invention is phase difference measurement of any two-way with frequency cosine signal.For the present invention the following steps are included: firstly, obtaining the frequency estimation of two paths of signals, building frequency signal estimates two-way cosine signal, and generates two-way error signal；Then, ART network is carried out to two-way cosine signal parameter using two-way error signal, obtains estimates of parameters；Finally, calculating the phase difference of two-way cosine signal using obtained two-way cosine signal estimates of parameters.Method for measuring phase difference structure of the present invention is simple, and strong operability has better noiseproof feature, can obtain unbiased phase difference estimation as a result, can effectively overcome the influence caused by phase difference measurement of signal frequency mismatch.
Description
Technical field
The present invention relates to field of signal processing, especially a kind of two-way that is directed to is with frequency cosine signal, in frequency mismatch condition
Under adaptive unbiased phase difference estimation method.
Background technique
In electric system, intelligence instrument, fault diagnosis, the fields such as telecommunications have more wide the estimation method of phase difference
It is general using for example, the coriolis mass flowmeters being used widely in Oil-Gas Storage and Transportation Engineering, measured fluid
Mass flow is dependent on the Phase Difference Measuring Precision of two-way cosine signal；For another example, the determination of power factor (PF) in the power system,
It is dependent on the estimated value of voltage and current and phase difference.Therefore, realization, which seems to the high-acruracy survey of phase difference, is even more important.It is existing
Some phase difference estimation methods, are mainly the following method:
(1) (bibliography [1]: Tu Yaqing, Li Ming, Shen Tingao wait to count to the phase difference estimation method based on DFT and DTFT
And DFT and DTFT phase difference measurement error analysis [J] the vibration and impact of negative frequency, 2015,34 (20): 85-91. is with reference to text
It offers [2]: Shen Tingao, Tu Yaqing, Haitao Zhang, waits sliding DTFT high-precision phase position difference measurements algorithm [J] of based on rectangle double window
Central South University's journal (natural science edition), 2015,46 (2): 554-560.):
Such phase difference estimation method precision is counted with sampled data, signal noise and frequency mismatch are related.It counts more
Precision is higher, but causes calculation amount bigger than normal, and real-time is not high, practical application.And influenced by signal noise and frequency mismatch,
Lead to this method there are influence factors excessive, calculates complicated problem.
(2) (the bibliography [3]: Tu Yaqing, Shen Ting chela, Li Ming, Zhang Hai of the method for measuring phase difference based on digital correlation
Non-integer-period signal phase difference Measurement Algorithm [J] the Chinese journal of scientific instrument of great waves based on multiple cross-correlation, 2014,35 (7):
1578-1585. bibliography [4]: Shen Y, Tu Y.Correlation theory-based signal processing
Method for CMF signals [J] .Measurement Science & Technology, 2016,27 (6):
065006.):
Such correlation method can obtain the accurate measurement of phase difference under conditions of signal carries out stringent integer-period sampled,
But due to the influence of signal noise and frequency mismatch, tends not to realize the integer-period sampled of signal, cause to adopt in non-integer-period
Phase Difference Measuring Precision declines under the conditions of sample, and there are biggish evaluated errors.
(3) (bibliography [5]: Yang Huiyue, Tu Yaqing, Haitao Zhang wait a kind of to the adaptive phase difference method based on LMS
The adaptive unbiased esti-mator method of vibration signal phase difference and application [J] based on LMS are vibrated and are impacted, 2016,35 (10): 55-
59.):
Such method seeks phase information by Hilbert transform, although the frequency information of undesired signal, not by
Frequency mismatch influences, but Hilbert transform pairs are weaker in the rejection ability of noise, causes the noise immunity of total algorithm insufficient, easily
Affected by noise, phase difference estimation precision is to be improved.
Summary of the invention
The purpose of the present invention is to propose to a kind of noise immunity is strong, estimated accuracy is high, and the good phase difference estimation method of applicability disappears
Influence except frequency estimation mismatch to phase difference measurement is solved suitable for any two-way with the phase difference measurement of frequency cosine signal
Main problem existing for existing method for measuring phase difference expands its application range.
The present invention proposes the adaptive unbiased phase difference estimation method under the conditions of a kind of frequency mismatch:
The present invention the following steps are included:
It is illustrated so that two-way is with frequency cosine signal as an example.
1) first via cosine signal x1(n) and the second road cosine signal x2(n) frequency is all ω0, frequency mismatch value is Δ
ω, it is assumed that estimated value ω=ω of given frequency0Δ ω constructs frequency signal cos ω n and sin ω n, and adaptive calculate is arranged
Step size mu=0.01 of method；
2) using frequency signal cos ω n and sin the ω n of building, if a1' (n), b1' (n) be respectively x1(n) parameter a in1
And b1Estimated value, then the estimated value of first via cosine signal be s1' (n)=a1’(n)cosωn-b1'(n)sinωn；
3) also with frequency signal cos ω n and sin the ω n of building, if a2' (n), b2' (n) be respectively x2(n) ginseng in
Number a2And b2Estimated value, then the estimated value of the second road cosine signal be s2' (n)=a2’(n)cosωn-b2’(n)sinωn
4) two-way cosine signal and its estimated value are utilized, error signal e is constructed1(n)=x1(n)-s1' (n) and e2(n)=x2
(n)-s2'(n)；
5) a is estimated respectively using adaptive iteration calculation1' (n+1), b1' (n+1), a2' (n+1) and b2' (n+1),
Iterative calculation method are as follows:
a1' (n+1)=a1’(n)+μe1(n)cosωn；
b1' (n+1)=b1’(n)-μe1(n)sinωn；
a2' (n+1)=a2’(n)+μe2(n)cosωn；
b2' (n+1)=b2’(n)-μe2(n)sinωn；
6) according to a of calculating1' (n+1), b1' (n+1), a2' (n+1) and b2' (n+1) value, calculate two-way known to frequency
With the phase difference of frequency cosine signal, then phase difference are as follows:
Detailed description of the invention
Below according to attached drawing and specific implementation, the present invention is further elaborated.By taking two-way cosine signal containing white noise as an example
It is illustrated.
Fig. 1 is method implementation flow chart.
In figure: 1 indicates to construct frequency signal using the frequency values obtained；2 indicate using frequency signal to first via cosine
Signal is estimated；3 indicate to estimate the second road cosine signal using frequency signal；4 indicate to produce into first via error letter
Number；5 indicate to generate the second tunnel error signal；6 indicate using first via error signal to parameter a1And b1Estimated；7 indicate benefit
With the second tunnel error signal to parameter a2And b2Estimated；8 indicate to calculate two-way with the phase difference of frequency cosine signal.
Specific embodiment
A specific embodiment of the invention is as follows:
Step 1: building frequency signal cos ω n and sin ω n, is arranged algorithm step-size μ=0.01；
Step 2: calculating first via cosine signal x1(n) estimated value s1' (n)=a1’(n)cosωn-b1’(n)sinω
n；
Step 3: calculating the second road cosine signal x2(n) estimated value s2' (n)=a2’(n)cosωn-b2’(n)sinω
n；
Step 3: building error signal e1(n)=x1(n)-s1' (n) and e2(n)=x2(n)-s2'(n)；
Step 4: estimating a respectively using adaptive iteration calculation1' (n+1), b1' (n+1), a2' (n+1) and b2’(n
+1)；
Step 5: according to a of calculating1' (n+1), b1' (n+1), a2' (n+1) and b2' (n+1) value, it calculates known to frequency
Phase difference of the two-way with frequency cosine signal.
Claims (1)
1. the adaptive unbiased phase difference estimation method under the conditions of frequency mismatch, applicable object is two-way known to frequency with more than frequency
String signal；
Method includes the following steps:
Step 1: constructing frequency signal cos ω n and sin ω n, algorithm step-size μ=0.01 using the signal frequency value ω obtained；
Step 2: calculating first via cosine signal x1(n) estimated value s1' (n)=a1’(n)cosωn-b1' (n) sin ω n,
Middle a1' (n), b1' (n) be respectively x1(n) parameter a in1And b1Estimated value；
Step 3: calculating the second road cosine signal x2(n) estimated value s2' (n)=a2’(n)cosωn-b2' (n) sin ω n,
Middle a2' (n), b2' (n) be respectively x2(n) parameter a in2And b2Estimated value；
Step 4: building error signal e1(n)=x1(n)-s1' (n) and e2(n)=x2(n)-s2'(n)；
Step 5: estimating a respectively using adaptive iteration calculation1' (n+1), b1' (n+1), a2' (n+1) and b2' (n+1),
Iterative calculation method are as follows:
a1' (n+1)=a1’(n)+μe1(n)cosωn；
b1' (n+1)=b1’(n)-μe1(n)sinωn；
a2' (n+1)=a2’(n)+μe2(n)cosωn；
b2' (n+1)=b2’(n)-μe2(n)sinωn；
Step 6: according to a of calculating1' (n+1), b1' (n+1), a2' (n+1) and b2' (n+1) value, calculate two-way known to frequency
With the phase difference of frequency cosine signal, then phase difference are as follows:
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710030146.4A CN106771600B (en) | 2017-01-07 | 2017-01-07 | Adaptive unbiased phase difference estimation method under the conditions of frequency mismatch |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710030146.4A CN106771600B (en) | 2017-01-07 | 2017-01-07 | Adaptive unbiased phase difference estimation method under the conditions of frequency mismatch |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106771600A CN106771600A (en) | 2017-05-31 |
| CN106771600B true CN106771600B (en) | 2019-07-30 |
Family
ID=58945949
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710030146.4A Expired - Fee Related CN106771600B (en) | 2017-01-07 | 2017-01-07 | Adaptive unbiased phase difference estimation method under the conditions of frequency mismatch |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106771600B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109030944B (en) * | 2018-07-16 | 2020-11-03 | 中国人民解放军陆军勤务学院 | Real sinusoidal signal phase difference estimation method based on frequency shift and DFT |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101018083A (en) * | 2007-02-14 | 2007-08-15 | 哈尔滨工业大学 | Dopla frequency shift compensation method in the MPSK mobile communication system |
| CN101825660A (en) * | 2010-05-05 | 2010-09-08 | 天津大学 | High-efficiency measurement method for sinusoidal signal frequency in undersampling and implementation device |
| CN102735937A (en) * | 2012-06-29 | 2012-10-17 | 中国科学院半导体研究所 | Method for measuring signal phase difference |
| CN103297190A (en) * | 2013-05-11 | 2013-09-11 | 哈尔滨工业大学深圳研究生院 | Deep-space communication oriented code-aided carrier phase synchronization system and deep-space communication oriented code-aided carrier phase synchronization method |
| CN103926463A (en) * | 2014-02-16 | 2014-07-16 | 中国人民解放军后勤工程学院 | Phase difference measurement method related to data continuation |
| CN103983261A (en) * | 2014-04-28 | 2014-08-13 | 北京大学 | A fiber optic gyroscope based on vector space analysis and a signal processing method thereof |
| CN104483549A (en) * | 2014-12-13 | 2015-04-01 | 中国人民解放军后勤工程学院 | Phase difference measuring method based on data extension and Hilbert transform |
| CN105092969A (en) * | 2015-07-23 | 2015-11-25 | 中国人民解放军后勤工程学院 | Phase frequency coupling method for phase difference estimation |
| CN105738696A (en) * | 2016-04-18 | 2016-07-06 | 天津大学 | Frequency estimation method and device for all-phase time-shift phase difference |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005030978A (en) * | 2003-07-09 | 2005-02-03 | Sony Corp | Phase difference measuring apparatus, phase difference measurement method, and testing apparatus |
-
2017
- 2017-01-07 CN CN201710030146.4A patent/CN106771600B/en not_active Expired - Fee Related
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101018083A (en) * | 2007-02-14 | 2007-08-15 | 哈尔滨工业大学 | Dopla frequency shift compensation method in the MPSK mobile communication system |
| CN101825660A (en) * | 2010-05-05 | 2010-09-08 | 天津大学 | High-efficiency measurement method for sinusoidal signal frequency in undersampling and implementation device |
| CN102735937A (en) * | 2012-06-29 | 2012-10-17 | 中国科学院半导体研究所 | Method for measuring signal phase difference |
| CN103297190A (en) * | 2013-05-11 | 2013-09-11 | 哈尔滨工业大学深圳研究生院 | Deep-space communication oriented code-aided carrier phase synchronization system and deep-space communication oriented code-aided carrier phase synchronization method |
| CN103926463A (en) * | 2014-02-16 | 2014-07-16 | 中国人民解放军后勤工程学院 | Phase difference measurement method related to data continuation |
| CN103983261A (en) * | 2014-04-28 | 2014-08-13 | 北京大学 | A fiber optic gyroscope based on vector space analysis and a signal processing method thereof |
| CN104483549A (en) * | 2014-12-13 | 2015-04-01 | 中国人民解放军后勤工程学院 | Phase difference measuring method based on data extension and Hilbert transform |
| CN105092969A (en) * | 2015-07-23 | 2015-11-25 | 中国人民解放军后勤工程学院 | Phase frequency coupling method for phase difference estimation |
| CN105738696A (en) * | 2016-04-18 | 2016-07-06 | 天津大学 | Frequency estimation method and device for all-phase time-shift phase difference |
Non-Patent Citations (3)
| Title |
|---|
| 一种基于LMS的振动信号相位差自适应无偏估计方法及应用;杨辉跃 等;《振动与冲击》;20160528;第35卷(第10期);全文 |
| 广义正态信号处理理论及在通信中应用的研究;唐洪;《中国优秀博硕士学位论文全文数据库 (博士) 信息科技辑》;20060815(第08期);全文 |
| 插值FFT和滑动DTFT的科氏流量计信号处理方法;易鹏 等;《计算机工程与应用》;20121023;第49卷(第5期);全文 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106771600A (en) | 2017-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105841762B (en) | The flow metering method and system of ultrasonic water meter | |
| Wang et al. | The gravity matching area selection criteria for underwater gravity-aided navigation application based on the comprehensive characteristic parameter | |
| CN105588883B (en) | Three-dimensional rock mechanical parameter acquisition method and system | |
| CN107305223B (en) | Improved phase difference frequency estimation method | |
| JP2014228549A (en) | Method and apparatus for evaluating satellite positioning quality | |
| CN109407501B (en) | Time interval measuring method based on relevant signal processing | |
| CN103926463B (en) | The method for measuring phase difference that a kind of Data extension formula is relevant | |
| CN106771600B (en) | Adaptive unbiased phase difference estimation method under the conditions of frequency mismatch | |
| CN109030944B (en) | Real sinusoidal signal phase difference estimation method based on frequency shift and DFT | |
| CN102193029B (en) | Method for measuring short-term frequency stability of unconventional sampling time | |
| Xu et al. | A new robust filtering method of GNSS/MINS integrated system for land vehicle navigation | |
| CN103793614B (en) | A kind of mutation filtering method | |
| Ventura et al. | Estimation of urban noise with the assimilation of observations crowdsensed by the mobile application Ambiciti | |
| CN105092969B (en) | A kind of phase frequency matching process of phase difference estimation | |
| CN105242225B (en) | A kind of calibration system and its calibration method of dynamic phasor measurement device | |
| WO2020220579A1 (en) | Method for quickly obtaining integer solution of satellite positioning in high reliability | |
| US9673862B1 (en) | System and method of analyzing crosstalk without measuring aggressor signal | |
| Jian et al. | Acoustic velocity measurement in seawater based on phase difference of signal | |
| Liu et al. | Signal processing of Coriolis mass flowmeters under gas-liquid two-phase flow conditions | |
| CN107064546B (en) | A kind of fluid velocity estimation method based on non-circular signal frequency algorithm for estimating | |
| CN105846822A (en) | Method and system for measuring mismatch error of time-interleaved analog-to-digital conversion system | |
| Hammes et al. | Semi-parametric geolocation estimation in NLOS environments | |
| CN112255647A (en) | Optimization method for trajectory measurement information post fusion processing | |
| Kang et al. | Performance analysis on digital phase difference measurement techniques for interferometer direction finder | |
| CN106218348B (en) | Leaf spring grasswort spacing detection method and apparatus |
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 | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: 401311 Logistics Engineering College, No. 20 North Road, University Town, Shapingba District, Chongqing Applicant after: PLA Military Service College Address before: 401311 Logistics Engineering College, No. 20 North Road, University Town, Shapingba District, Chongqing Applicant before: Rare Service Engineering College of Chinese PLA. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190730 Termination date: 20200107 |