CN108535542A - A kind of peak-seeking phase detecting method - Google Patents
A kind of peak-seeking phase detecting method Download PDFInfo
- Publication number
- CN108535542A CN108535542A CN201810332080.9A CN201810332080A CN108535542A CN 108535542 A CN108535542 A CN 108535542A CN 201810332080 A CN201810332080 A CN 201810332080A CN 108535542 A CN108535542 A CN 108535542A
- Authority
- CN
- China
- Prior art keywords
- signal
- peak
- point
- window function
- correlation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 230000003044 adaptive effect Effects 0.000 claims abstract description 13
- 238000001914 filtration Methods 0.000 claims description 9
- 230000006835 compression Effects 0.000 claims description 6
- 238000007906 compression Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 238000000638 solvent extraction Methods 0.000 claims description 3
- 230000007423 decrease Effects 0.000 claims description 2
- 230000011218 segmentation Effects 0.000 claims description 2
- 230000010363 phase shift Effects 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 3
- 239000010931 gold Substances 0.000 description 3
- 229910052737 gold Inorganic materials 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 210000003739 neck Anatomy 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
- Complex Calculations (AREA)
Abstract
A kind of peak-seeking phase detecting method, the phase changing capacity for detecting analog signal all the way, this method are by finding signal peak location determination signal phase, phase demodulation process：Analog signal is sampled and is converted into Serial No., find out the maximum point position in Serial No., the raised cosine adaptive using frequency is as window function, cross correlation filter is carried out to the signal near maximum point position, maximum point position range is constantly reduced using Fibonacci method simultaneously, again maximum point position, i.e. peak position is accurately positioned.Since the Fibonacci method of use can constantly reduce maximum point position range, the feature which has operand small, real-time, in addition, as a result of cross correlation filter, which has noise very high tolerance.
Description
Technical field
The present invention relates to signal disposal and analysis technical field more particularly to a kind of peak-seeking phase detecting methods.
Background technology
Phase acquisition method is for necks such as Laser Range Finding Based on Phase, laser interferometer, optical fibre gyro, laser Doppler measurings
Domain.Typical phase acquisition method has fourier transform method, phase shift method, correlation method etc..Fourier transform method, that is, Spectral Analysis Method,
Time-domain signal is transformed into frequency domain using discrete fourier variation principle, basis signal phase-frequency characteristic realizes phase detection, Fu Li
There is very serious spectral leakage defect in leaf transformation method, measurement accuracy is limited by frequency resolution.Phase shift method uses Hilbert
Transformation carries out signal 90 ° of phase shift,
Phase demodulating is carried out according to original signal and the signal after phase shift, the phase shifting accuracy of phase-shifting method will lead to phase detection
Error, while it is endless that Hilbert transform, which requires data length,.
The principle of correlation method is the remaining of zero propagation cross-correlation function value and its phase difference of the sinusoidal signal of two same frequencys
String value is directly proportional.Due to the correlation very little of noise and useful signal, there are good noise inhibiting ability, but phase in this way
Pass method needs to carry out convolution algorithm to entire signal and stencil function, and operand is big, and algorithm real-time is poor.
Invention content
The present invention provides a kind of peak-seeking phase detecting method, and to overcome, operand of the existing technology is big and algorithm real-time is poor
The problem of.
To achieve the goals above, technical solution provided by the invention is：
A kind of peak-seeking phase detecting method, comprises the steps of：
Step 1：Analog signal is sampled and is converted into Serial No.；
Step 2：Method used by signal peak position is according to a preliminary estimate：According to the signal period, using larger space pair
Signal is sampled, and the sampling interval is about the 1/10 of signal period, it is ensured that change amount signal is much larger than noise, adjacent by judging
The size of two sampled points determines the variable condition of signal, if signal is decline state by propradation variation, it can be inferred that
There are a maximum between two sampling intervals, and the maximum value between two sampling intervals is denoted as and estimates peak value Max；
Step 3：The adaptive raised cosine of frequency is template, centered on maximum point position, determines one greatly
Value point position range, carries out digital filtering, while reducing maximum position range using Fibonacci method using cross-correlation；
Step 4：Using the time difference of two adjacent peaks positions as signal period T, the variable quantity of peak position is Δ T,
It is corresponding to phase changing capacity
Further, the raised cosine that the cross-correlation digital filtering in the step 3 uses frequency adaptive is as window
Function, raised cosine windows are：
Wherein, M+1 is the window function period, and n is window function time series, is an array by raised cosine windows tabulation,
For the digital signal of different cycles, time compression is carried out to the array, realizes that the frequency of window function is adaptive, time compression pair
The relational expression is answered to be：
Wherein, n is raised cosine windows time series, and M+1 is the window function period, and T is the signal period, and t is to be risen after compressing
Cosine window function time series, to the computing cross-correlation of certain point, by a cycle signal using centered on the point as input
Signal x (t), window function are w (t), and computing cross-correlation is
Y (τ)=∫ w ' (t) x (t+ τ) dt
In formula, τ is delay time, corresponds to maximum position range, and computing cross-correlation is that the sliding of Serial No. is multiplied,
Output signal y (τ) is the cross correlation results of Serial No. within the scope of maximum position, determines the maximum position of output signal y (τ)
It sets, as the filtered peak position of signal.
Further, maximum point position range is split using Fibonacci method in the step 3, specifically
Partitioning scheme is：It is M to estimate peak position, and it is [M+ Δs, M- Δs] that peak position, which corresponds to section, with the left and right endpoint position of the range
It sets and carries out golden section on the basis of distinguishing, corresponding left golden section point and right golden section point can be obtained, correspond to left gold
Divide domain and right golden section domain, cross correlation filter is carried out to golden section point, selection filtering exports corresponding to high cut-point
Region carries out golden section again, and can all be reduced by often carrying out a golden section peak ranges by 0.618 times.
Compared with prior art, beneficial effects of the present invention are as follows：
1, this method only needs to carry out phase demodulation to limited signal segment, and phase is carried out to signal cutout using raised cosine windows
Operation is closed, reduces limit signal and carries out spectrum energy leakage caused by periodic extension when related operation, reduce to signal
The requirement of length.
2, when cross correlation filter, window function frequency is better closer to filter effect with signal frequency.This method is to raised cosine
Function is improved, it can be achieved that its frequency is adaptive, for the signal of different frequency, can all have good filter effect.
3, the present invention reduces peak position range using Fibonacci method, reduces the test number (TN) of cross correlation filter, in turn
The operand for reducing algorithm, improves the real-time of algorithm.
4, the present invention is filtered signal using cross-correlation, can accurately obtain peak point position, high to noise tolerance.
Description of the drawings
Fig. 1 is the functional block diagram of the peak-seeking phase discriminator built with phase detecting method of the present invention；
Fig. 2 is cross correlation filter physical model；
Fig. 3 is Fibonacci method implementation steps flow chart；
Fig. 4 is sampled data and phase detecting method handling result.
Specific implementation mode
To make the objectives, technical solutions, and advantages of the present invention clearer, below in conjunction with attached drawing to the present invention make into
It is described in detail to one step, it is clear that described embodiments are only a part of the embodiments of the present invention, rather than whole implementation
Example.Based on the embodiments of the present invention, obtained by those of ordinary skill in the art without making creative efforts
All other embodiment, shall fall within the protection scope of the present invention.
A kind of peak-seeking phase detecting method provided by the invention, samples analog signal and is converted into Serial No., found out
Maximum point position in Serial No., using the adaptive raised cosine of frequency as template, near maximum point position
Signal carries out cross correlation filter, while constantly reducing maximum point position range using Fibonacci method, is accurately positioned again
Maximum point position, and then obtain signal phase.
Referring to Fig. 1, the present invention is by analog-to-digital conversion, peak value, Frequency Estimation, peak ranges golden section, cross correlation filter, peak
Value point amendment step is constituted.It comprises the steps of：
Step 1：Analog signal is sampled and is converted into Serial No.；
Step 2：Maximum point position in Serial No. according to a preliminary estimate；
Step 3：The adaptive raised cosine of frequency is template, centered on maximum point position, determines one greatly
Value point position range, digital filtering is carried out using cross-correlation, to which maximum point position be accurately positioned, while using gold point
It cuts method and reduces maximum position range, reduce cross correlation filter operand；
Step 4：Using the time difference of two adjacent peaks positions as signal period T, the variable quantity of peak position is Δ T,
It is corresponding to phase changing capacity
In step 3, used cross correlation filter physical model, referring to Fig. 2, using the adaptive raised cosine of frequency
As window function, raised cosine windows are：
Wherein, M+1 is the window function period, and n is window function time series.
It is an array or pointer by raised cosine windows tabulation, for the pending serial number of different cycles, to the array
Time compression is carried out, realizes that the frequency of window function is adaptive.Time compresses correspondence formula：
Wherein, n is raised cosine windows time series, and M+1 is the window function period, and T is the signal period, and t is to be risen after compressing
Cosine window function time series, to the computing cross-correlation of certain point, by a cycle signal using centered on the point as input
Signal x (t), window function are w (t), and computing cross-correlation is
Y (τ)=∫ w ' (t) x (t+ τ) dt
In formula, τ is delay time, corresponds to maximum position range, and computing cross-correlation is that the sliding of Serial No. is multiplied,
Output signal y (τ) is the cross correlation results of Serial No. within the scope of maximum position, determines the maximum position of output signal y (τ)
It sets, as the filtered peak position of signal.
In step 3, used Fibonacci method implementation steps, referring to Fig. 3, specific partitioning scheme is：Estimate peak value position
Be set to M, it is [M+ Δs, M- Δs] that peak position, which corresponds to section, by the left and right endpoint location of the range respectively on the basis of carry out gold
Segmentation, can be obtained corresponding left golden section point and right golden section point, correspond to left golden section domain and right golden section domain,
Cross correlation filter is carried out to golden section point, selection filtering exports region corresponding to high cut-point and carries out golden section again,
Can all be reduced by often carrying out a golden section peak ranges by 0.618 times.
Referring to Fig. 4, sampled data and phase detecting method handling result figure, solid line plot for original noisy signal, the song
Point on line is labeled as the peak position point primarily determined.Dotted line plots mutual to be carried out to the signal near peak position
Filtered signal is closed, the point on the curve is labeled as identified peak position after filtering.It can be seen that：Adaptive raised cosine
Computing cross-correlation has good filter effect, improves the precision of peak detection.
Use above specific case is illustrated the present invention, is merely used to help understand the present invention, not limiting
The system present invention.For those skilled in the art, according to the thought of the present invention, can also make several simple
It deduces, deform or replaces.
Claims (3)
1. a kind of peak-seeking phase detecting method, which is characterized in that comprise the steps of：
Step 1：Analog signal is sampled and is converted into Serial No.；
Step 2：Method used by signal peak position is according to a preliminary estimate：According to the signal period, using larger space to signal
It is sampled, the sampling interval is about the 1/10 of signal period, it is ensured that change amount signal is much larger than noise, two neighboring by judging
The size of sampled point determines the variable condition of signal, if signal is decline state by propradation variation, it can be inferred that two
There are a maximum between sampling interval, and the maximum value between two sampling intervals is denoted as and estimates peak value Max；
Step 3：The adaptive raised cosine of frequency is that template determines a maximum point centered on maximum point position
Position range carries out digital filtering using cross-correlation, while reducing maximum position range using Fibonacci method；
Step 4：Using the time difference of two adjacent peaks positions as signal period T, the variable quantity of peak position is Δ T, corresponding
It is in phase changing capacity
2. peak-seeking phase detecting method as described in claim 1, which is characterized in that
Cross-correlation digital filtering in the step 3 is using the adaptive raised cosine of frequency as window function, raised cosine window
Function is：
Wherein, M+1 is the window function period, and n is window function time series, is an array by raised cosine windows tabulation, for
The digital signal of different cycles carries out time compression to the array, realizes that the frequency of window function is adaptive, and time compression is corresponding to close
It is that formula is：
Wherein, n is raised cosine windows time series, and M+1 is the window function period, and T is the signal period, and t is raised cosine after compression
Window function time series, to the computing cross-correlation of certain point, by a cycle signal using centered on the point as input signal x
(t), window function is w (t), and computing cross-correlation is
Y (τ)=∫ w*(t)x(t+τ)dt
In formula, τ is delay time, corresponds to maximum position range, and computing cross-correlation is that the sliding of Serial No. is multiplied, output
Signal y (τ) is the cross correlation results of Serial No. within the scope of maximum position, determines the maximum position of output signal y (τ),
The as filtered peak position of signal.
3. the peak-seeking phase detecting method as described in claim 1 and 2, which is characterized in that
Maximum point position range is split using Fibonacci method in the step 3, specific partitioning scheme is：In advance
It is M to estimate peak position, and it is [M+ Δs, M- Δs] that peak position, which corresponds to section, on the basis of the left and right endpoint location difference of the range
Golden section is carried out, corresponding left golden section point and right golden section point can be obtained, corresponds to left golden section domain and right Huang
Gold segmentation domain carries out cross correlation filter to golden section point, and selection filtering exports region corresponding to high cut-point and carries out again
Golden section, can all be reduced by often carrying out a golden section peak ranges by 0.618 times.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810332080.9A CN108535542B (en) | 2018-04-13 | 2018-04-13 | Peak-seeking phase discrimination method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810332080.9A CN108535542B (en) | 2018-04-13 | 2018-04-13 | Peak-seeking phase discrimination method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108535542A true CN108535542A (en) | 2018-09-14 |
| CN108535542B CN108535542B (en) | 2020-06-19 |
Family
ID=63480455
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810332080.9A Expired - Fee Related CN108535542B (en) | 2018-04-13 | 2018-04-13 | Peak-seeking phase discrimination method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108535542B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111623958A (en) * | 2020-05-18 | 2020-09-04 | 长春欧意光电技术有限公司 | Wavelet peak-peak value extraction method in interference signal |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1371501A (en) * | 1999-08-26 | 2002-09-25 | 艾利森公司 | Method and apparatus for locating peak of function using filter bank |
| CN101237524A (en) * | 2008-03-03 | 2008-08-06 | 中国科学院光电技术研究所 | Image noise removing method capable of reserving high-frequency information |
| CN102984099A (en) * | 2012-11-23 | 2013-03-20 | 陕西理工学院 | Algorithm for digital predistortion fraction time delay estimation and signal alignment and achievement thereof |
-
2018
- 2018-04-13 CN CN201810332080.9A patent/CN108535542B/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1371501A (en) * | 1999-08-26 | 2002-09-25 | 艾利森公司 | Method and apparatus for locating peak of function using filter bank |
| CN101237524A (en) * | 2008-03-03 | 2008-08-06 | 中国科学院光电技术研究所 | Image noise removing method capable of reserving high-frequency information |
| CN102984099A (en) * | 2012-11-23 | 2013-03-20 | 陕西理工学院 | Algorithm for digital predistortion fraction time delay estimation and signal alignment and achievement thereof |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111623958A (en) * | 2020-05-18 | 2020-09-04 | 长春欧意光电技术有限公司 | Wavelet peak-peak value extraction method in interference signal |
| CN111623958B (en) * | 2020-05-18 | 2021-11-12 | 长春欧意光电技术有限公司 | Wavelet peak-peak value extraction method in interference signal |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108535542B (en) | 2020-06-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108955734B (en) | Cavity length demodulation method of optical fiber F-P temperature/pressure composite sensor | |
| CN110133598B (en) | FrFT-based quick estimation method for parameters of linear frequency modulation signals | |
| CN102043091B (en) | Digitized high-precision phase detector | |
| CN109490862B (en) | Carrier frequency estimation method based on phase difference statistical spectrum | |
| CN106546817B (en) | A kind of Frequency Estimation and energy state postulate with feedback function | |
| CN108318737A (en) | A kind of high accuracy number phase detecting method | |
| CN109214371A (en) | A kind of sinusoidal signal filtering system based on fast Fourier transform | |
| CN101320060A (en) | Fast phase meter | |
| CN103823177B (en) | Performance of filter detection method based on THE DESIGN OF WINDOW FUNCTION and system | |
| CN111693136B (en) | Acoustic surface wave resonator frequency estimation algorithm adopting echo signal autocorrelation phase spectrum | |
| CN114460527B (en) | Correlation degree continuation Hilbert phase-shifting electronic transformer calibrator source tracing method and system | |
| CN102495280B (en) | Anti-noise wide-range frequency measurement method and phase locking frequency meter | |
| CN108535542A (en) | A kind of peak-seeking phase detecting method | |
| CN108181486B (en) | The processing method and processing device of acceleration signal | |
| CN108900445B (en) | Method and device for estimating signal symbol rate | |
| CN113156206A (en) | Time-frequency combined noise-containing signal parameter estimation new algorithm | |
| CN115265691B (en) | Coriolis flowmeter vibration frequency tracking method and system | |
| Nannan | The application of an improved wavelet threshold function in de-noising of heart sound signal | |
| CN109521269A (en) | A kind of am signals digitlization frequency measuring method | |
| CN105738698A (en) | Harmonic parameter estimation algorithm based on center frequency shift | |
| CN112883787B (en) | Short sample low-frequency sinusoidal signal parameter estimation method based on spectrum matching | |
| CN108020818A (en) | A kind of sine pulse signal distance measuring method under noise based on sliding DFT | |
| US9759751B1 (en) | Line cycle correlated spectral analysis for power measurement systems | |
| CN110582697B (en) | Apparatus and method for measuring and tracking the quantity or concentration of a component in a fluid | |
| CN107576842B (en) | Broadband synchronous sampling method |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200619 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |