CN105353187A - Real-time waveform reconstruction method based on three-point asynchronous sampling - Google Patents
Real-time waveform reconstruction method based on three-point asynchronous sampling Download PDFInfo
- Publication number
- CN105353187A CN105353187A CN201510782283.4A CN201510782283A CN105353187A CN 105353187 A CN105353187 A CN 105353187A CN 201510782283 A CN201510782283 A CN 201510782283A CN 105353187 A CN105353187 A CN 105353187A
- Authority
- CN
- China
- Prior art keywords
- sampling
- real
- waveform
- time
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R13/00—Arrangements for displaying electric variables or waveforms
Abstract
The invention provides a real-time waveform reconstruction method based on three-point asynchronous sampling. The real-time waveform reconstruction method comprises the following steps that equal interval sampling of signals is performed according to fixed sampling frequency; frequency components apart from fundamental wave are filtered by a digital filter; a waveform reconstruction parameter deltak is calculated by selecting three sampling points after filtering processing; and software resampling is performed on an original sampling sequence by utilizing the parameter, a time domain sampling location is changed in real time and a synchronous sampling sequence is obtained, and a signal waveform is reconstructed. According to the method, error caused by incomplete synchronous sampling can be avoided so that real-time performance is high, calculation is convenient and realization is easy.
Description
Technical field
The present invention relates to power system signal process field, be specifically related to a kind of real-time waveform reconstructing method based on 3 asynchronous-samplings.
Background technology
Along with the development of Power Electronic Technique, a large amount of non-linear equipment comes into operation in electric system, causes the quality of power supply to decline.The acquisition of real-time waveform data, improves quality of power supply tool and is of great significance the actual conditions understanding electric system.
Obtain real-time waveform data usually to need to sample to continuous wave, convert simulating signal to discrete Serial No., just obtain relevant information by digital signal processing.In order to maintain the synchronous of sampling, prior art adjusts the length in sampling period usually in real time according to the change of signal frequency, mainly comprises software synchronization and hardware synchronization two kinds of modes.
Frequency due to electric system is real-time change, and the sampling period determines to there is certain hysteresis quality according to Wave data before.When adjusting the sampling period by the mode of timer interruption, due to the restriction by timer minimum resolution, the timing value not necessarily integer of actual computation, there is certain quantization error in the sampling period simultaneously.
Summary of the invention
Fundamental purpose of the present invention is to provide a kind of real-time waveform reconstructing method based on 3 asynchronous-samplings, can be solved the synchronous error problem of real-time waveform data acquisition by the method.
Technical scheme of the present invention is as follows:
Based on a real-time waveform reconstructing method for 3 asynchronous-samplings, comprise the following steps:
(1) with a fixed frequency
equal interval sampling is carried out to signal, sampled point is saved as crude sampling sample.
(2) sampled signal is passed through the component beyond digital filter filtering first-harmonic.
(3) from the sampled signal after filtering process, get three sampled points calculate Waveform Reconstructing parameter
.
(4) according to parameter
software resampling is carried out to original non-same sample sequence, changes time-domain sampling position, obtain a synchronized sampling sequence, the real-time waveform of reconstruction signal.
Described method, the digital filter in step (2) is the logical FIR filter of band.
Described method, waveform reconstruction parameter in step (3)
for signal actual frequency
with reference frequency
ratio.
Described method, calculates Waveform Reconstructing parameter and comprises in step (3):
Choose three to be spaced apart
sampled value
,
,
, wherein
for the sampling number of one-period,
;
Determine Waveform Reconstructing parameter
, wherein
,
.
Described method, the software resampling in step (4) is by three non-synchronous sampling points
,
,
calculate each synchronously sampled data
(
), realize real-time waveform reconstruct, wherein
, INT [], for rounding operator, specifically comprises:
Two sample sequences get identical starting point, of reconfiguration waveform
individual data point and corresponding
the time interval between individual non-synchronous sampling point is
, wherein
for the fixing non-synchronous sampling cycle,
for the synchronized sampling cycle obtained according to signal actual frequency, the ratio between this time interval and fixing non-synchronous sampling cycle
.
The weighting coefficient of three non-synchronous sampling points is determined by the method for Quadratic interpolation:
，
，
。
Determine the synchronously sampled data value reconstructed:
.
Beneficial effect of the present invention: avoid the error that sampling not exclusively synchronously causes, do not need to change sample frequency, without the need to longer sample sequence, real-time synchronously sampled data can be obtained rapidly and accurately, be convenient to carry out real-time analysis to electric power signal in the electric system of frequency moment change, thus understand the actual conditions of electric system more exactly.
Accompanying drawing explanation
Fig. 1 is the process flow diagram of the real-time waveform reconstructing method realized in the present invention based on 3 asynchronous-samplings.
Embodiment
As shown in Figure 1, the real-time waveform reconstructing method based on 3 asynchronous-samplings of the present invention comprises the following steps:
(1) with a fixed frequency
carry out equal interval sampling to signal, obtain N+2 sampled point, sampled point is saved as crude sampling sample, N is the sampling number of one-period simultaneously, is taken as the integral multiple of 4.
(2) sampled signal is carried out filtering process by the logical FIR filter of band, extract fundamental signal.
(3) Waveform Reconstructing parameter is calculated
, concrete steps are as follows:
From the signal processed after filtering, choose three be spaced apart
sampled value
,
and
, wherein
.
Determine Waveform Reconstructing parameter
, wherein
,
.
(4) crude sampling sample is extracted, by three continuous print non-synchronous sampling points
,
,
calculate each synchronously sampled data
(
), complete the reconstruct of real time data waveform, wherein
, INT [] is for rounding operator, and concrete steps are as follows:
Two sample sequences get identical starting point, calculate the of reconfiguration waveform
individual data point and corresponding the
the time interval between individual non-synchronous sampling point and the ratio between the fixing non-synchronous sampling cycle:
, wherein
represent the fixing non-synchronous sampling cycle,
represent the synchronized sampling cycle obtained according to signal actual frequency.
By the method for Quadratic interpolation, determine that the weighting coefficient that three non-synchronous sampling points are corresponding is respectively:
，
，
。
Quadratic interpolation formula is:
, wherein,
for in reconstruct Wave data the
in individual point and non-synchronous sampling sequence
the time interval between individual sampled point, namely
.
Calculate each synchronously sampled data
, realize the reconstruct of real-time waveform.
In the embodiment of the present invention, step (3) determines that the Method And Principle of Waveform Reconstructing parameter is specific as follows.
The fundamental signal obtained after filtering can be described as:
, wherein,
for signal amplitude,
for signal frequency,
for the phasing degree of signal.Adopt fixed frequency to sample to continuous time signal, obtain discrete sampling sequence:
, wherein
,
for the sampling number of one-period,
for the ratio of signal actual frequency and reference frequency
, select power frequency as reference frequency.
Fixed sampling frequency in described step (1)
Two are spaced apart
sampled value sum can be expressed as:
。
According to two jiaos and, difference sine and cosine formula, can obtain:
。
when being taken as even number,
, wherein
,
.
Sine function is carried out Maclaurin series expansion, can obtain:
。
Retain progression Section 1 and obtain approximate value
.
The present invention retains first two, can obtain:
.
Therefore Waveform Reconstructing parameter
.
Claims (5)
1., based on a real-time waveform reconstructing method for 3 asynchronous-samplings, it is characterized in that, the method includes the steps of:
(1) with a fixed frequency
equal interval sampling is carried out to signal, sampled point is saved as crude sampling sample;
(2) sampled signal is passed through the component beyond digital filter filtering first-harmonic;
(3) from the sampled signal after filtering process, get three sampled points calculate Waveform Reconstructing parameter
;
(4) according to parameter
software resampling is carried out to original non-same sample sequence, changes time-domain sampling position, obtain a synchronized sampling sequence, the real-time waveform of reconstruction signal.
2. method according to claim 1, is characterized in that, the digital filter in step (2) is the logical FIR filter of band.
3. method according to claim 1, is characterized in that, the Waveform Reconstructing parameter in step (3)
for signal actual frequency
with reference frequency
ratio.
4. method according to claim 1, is characterized in that, calculates Waveform Reconstructing parameter and comprise described in step (3):
Choose three to be spaced apart
sampled value
,
,
, wherein
for the sampling number of one-period,
;
Determine Waveform Reconstructing parameter
, wherein
,
.
5. method according to claim 1, is characterized in that, the software resampling in step (4) is by three non-synchronous sampling points
,
,
calculate each synchronously sampled data
(
), realize real-time waveform reconstruct, wherein
, INT [], for rounding operator, specifically comprises:
Two sample sequences get identical starting point, of reconfiguration waveform
individual data point and corresponding
the time interval between individual non-synchronous sampling point is
, wherein
for the fixing non-synchronous sampling cycle,
for the synchronized sampling cycle obtained according to signal actual frequency, the ratio between this time interval and fixing non-synchronous sampling cycle
;
Weighting coefficient by the method determination non-synchronous sampling point of Quadratic interpolation:
,
,
;
Determine the synchronously sampled data value reconstructed:
.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510782283.4A CN105353187B (en) | 2015-11-16 | 2015-11-16 | A kind of real-time waveform reconstructing method based on 3 asynchronous-samplings |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510782283.4A CN105353187B (en) | 2015-11-16 | 2015-11-16 | A kind of real-time waveform reconstructing method based on 3 asynchronous-samplings |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105353187A true CN105353187A (en) | 2016-02-24 |
| CN105353187B CN105353187B (en) | 2018-06-26 |
Family
ID=55329189
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510782283.4A Active CN105353187B (en) | 2015-11-16 | 2015-11-16 | A kind of real-time waveform reconstructing method based on 3 asynchronous-samplings |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105353187B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106405210A (en) * | 2016-09-05 | 2017-02-15 | 深圳紫光继保测控技术有限公司 | Sampling data anomaly detection method and system |
| CN109582176A (en) * | 2018-11-30 | 2019-04-05 | 北京集创北方科技股份有限公司 | A kind of touch screen anti-noise method and device |
| CN110888020A (en) * | 2019-11-30 | 2020-03-17 | 国网辽宁省电力有限公司锦州供电公司 | Synchronous and asynchronous alternating current mixed sampling method based on threshold value |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4928251A (en) * | 1987-04-27 | 1990-05-22 | Hewlett-Packard Company | Method and apparatus for waveform reconstruction for sampled data system |
| CN101126772A (en) * | 2007-09-11 | 2008-02-20 | 电子科技大学 | Oscilloscope high speed signal reconstruction method |
| CN101241145A (en) * | 2008-03-11 | 2008-08-13 | 江苏绿扬电子仪器集团有限公司 | Method for waveform sampling and reconstitution for oscillographs |
-
2015
- 2015-11-16 CN CN201510782283.4A patent/CN105353187B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4928251A (en) * | 1987-04-27 | 1990-05-22 | Hewlett-Packard Company | Method and apparatus for waveform reconstruction for sampled data system |
| CN101126772A (en) * | 2007-09-11 | 2008-02-20 | 电子科技大学 | Oscilloscope high speed signal reconstruction method |
| CN101241145A (en) * | 2008-03-11 | 2008-08-13 | 江苏绿扬电子仪器集团有限公司 | Method for waveform sampling and reconstitution for oscillographs |
Non-Patent Citations (1)
| Title |
|---|
| AZIME CAN等: "ASYNCHRONOUS SAMPLING AND RECONSTRUCTION OF SPARSE SIGNALS", 《20TH EUROPEAN SIGNAL PROCESSING CONFERENCE》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106405210A (en) * | 2016-09-05 | 2017-02-15 | 深圳紫光继保测控技术有限公司 | Sampling data anomaly detection method and system |
| CN109582176A (en) * | 2018-11-30 | 2019-04-05 | 北京集创北方科技股份有限公司 | A kind of touch screen anti-noise method and device |
| CN109582176B (en) * | 2018-11-30 | 2021-12-24 | 北京集创北方科技股份有限公司 | Anti-noise method and device for touch screen |
| CN110888020A (en) * | 2019-11-30 | 2020-03-17 | 国网辽宁省电力有限公司锦州供电公司 | Synchronous and asynchronous alternating current mixed sampling method based on threshold value |
| CN110888020B (en) * | 2019-11-30 | 2022-03-01 | 国网辽宁省电力有限公司锦州供电公司 | Synchronous and asynchronous alternating current mixed sampling method based on threshold value |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105353187B (en) | 2018-06-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103196547B (en) | One realizes the synchronous Computed order tracking analytical approach of characteristic of rotating machines vibration signal | |
| CN105353187A (en) | Real-time waveform reconstruction method based on three-point asynchronous sampling | |
| CN104897962B (en) | The short sample high-precision frequency measurement method of simple signal and its device based on coprime perception | |
| CN106406174B (en) | A kind of multimode multichannel collecting synchronization system and working method | |
| CN203287435U (en) | A micro electrical network harmonic wave and inter-harmonic wave test apparatus based on an STM32F107VCT6 | |
| CN102539915A (en) | Method for accurately calculating power harmonic wave parameters through adopting time delay Fourier transform frequency measurement method | |
| CN103257271A (en) | Device and method for detecting micro grid harmonic wave and inter-harmonics based on STM32F107VCT6 | |
| CN103969508B (en) | A kind of Electric Power Harmonic Analysis method high-accuracy in real time and device | |
| CN109298362B (en) | Three-phase synchronization precision testing method, device and system for distribution line fault indicator | |
| CN103983849B (en) | A kind of Electric Power Harmonic Analysis method of real-time high-precision | |
| CN106911624A (en) | A kind of channel compensation calibration method and system | |
| CN103117972A (en) | Method and device for analyzing vector signals | |
| CN102236048A (en) | Method for measuring phasor frequency of electric system | |
| CN104135326B (en) | It is applicable to radio frequency testing device and the method for EDGE system 8PSK signal | |
| CN104717051B (en) | A kind of Interpolate estimation method in parallel demodulation bit synchronization | |
| CN103957009A (en) | Method for compensating for low-pass filter of compressed sampling system | |
| CN203132688U (en) | Device for realizing synchronous order-domain trace analysis on rotating mechanical vibration signals | |
| CN109632070B (en) | Digital array time domain quasi-synchronous calibration method based on Newton interpolation | |
| CN105119702A (en) | Timing synchronization method and device for signal processing | |
| CN105300437A (en) | VLBI baseband signal Decimal delay simulation method | |
| CN111625769A (en) | PMU-SCADA data time synchronization and fusion method based on Lagrange interpolation and cubic exponential smoothing | |
| CN104852744A (en) | Sinusoidal signal-based signal reconstruction method for acquiring sensing matrix under modulated wideband converter system | |
| JP2013135401A (en) | Resampling device, digital modulation signal generation device using the same, and resampling method | |
| Olkkonen et al. | Fractional delay filter based on the B-spline transform | |
| CN104091323A (en) | Weld tracking feature signal extracting method based on MP-NBFE method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: 310053 Zhejiang city of Hangzhou province high and New Technology Industrial Park, Jiang Dongxin Road No. 69 building 4 Building 801 room 8 Applicant after: Hangzhou Gaea Electric Co., Ltd. Address before: Hangzhou City, Zhejiang Province, Binjiang District Jiangnan Road 310053 No. 3880 Huarong Times Building 15 Jiahe Electric Company Applicant before: Hangzhou Gaea Electric Co., Ltd. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |