CN102707252B - Method for removing time base flutter of high-speed sampling oscilloscope - Google Patents
Method for removing time base flutter of high-speed sampling oscilloscope Download PDFInfo
- Publication number
- CN102707252B CN102707252B CN201210166086.6A CN201210166086A CN102707252B CN 102707252 B CN102707252 B CN 102707252B CN 201210166086 A CN201210166086 A CN 201210166086A CN 102707252 B CN102707252 B CN 102707252B
- Authority
- CN
- China
- Prior art keywords
- signal
- time base
- imf component
- sampling
- base flutter
- 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
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Complex Calculations (AREA)
Abstract
The invention provides a method for removing time base flutter of a high-speed sampling oscilloscope, which relates to the data processing field of instruments and meters. The method is characterized in that an EEMD (Ensemble Empirical Mode Decomposition) method is used to decompose in the data processing process, IMF (Intrinsic Mode Function) components are subjected to frequency-domain analysis to screen the useful signal, and the noise signals are removed to reconstruct the new signals, thereby realizing the effect for removing the time base flutter.
Description
Technical field
The present invention relates to a kind of method of removing the time base flutter of high-speed sampling oscillograph, belong to the data processing field of instrument and meter.
Background technology
High-speed sampling oscillograph for 10GHz and above bandwidth there will be time base flutter in the process of carrying out signal measurement and calibration.Time base flutter is because reality sampling is constantly constantly unequal caused with ideal sampling.In actual sampling process, sample interval is not often accurately to equal desirable sample time, but and between ideal value, has a time deviation.In precision measurement, high deadbeat sampling oscilloscope is during to sample of signal, and measurement result more or less all exists time base flutter, and it has a strong impact on test result.For example, in Nose-to-Nose collimation technique, an important step is exactly the removal to time base flutter.
Summary of the invention
The invention provides a kind of data processing method of high-speed sampling oscillograph being removed the time base flutter in signal source sampling based on population mean empirical mode decomposition (Ensemble Empirical Mode Decomposition, EEMD) method.First the signal with high-speed sampling oscillograph, a signal source being produced carries out the sampling of N point, samples M time.Wherein, N represents sampling number, and M represents the sampling number for N point.The fewer time base flutter of sampling number is more obvious, uses EEMD method only need to sample and obtains one group of data that comprise N numerical value for 1 time.Wherein when N > 1000, time base flutter is obvious.The N obtaining numerical value is imported to computing machine, carry out data processing.Method provided by the invention, is characterised in that in data handling procedure, by EEMD method, obtains intrinsic mode function (Intrinsic Mode Function, IMF), filter out useful signal, noise signal is removed, and then the signal that makes new advances of reconstruct, reach the effect of removing time base flutter.As shown in Figure 1, concrete steps are a kind of method flow diagram of removing the time base flutter of high-speed sampling oscillograph:
Step 1, with the signal that high-speed sampling oscillograph produces signal source, once sample, obtain original signal s (t), in s (t), have the data of N numerical value.The N obtaining numerical value is imported to computing machine and carry out data processing, wherein N is sampling number, and t is the sampling time.
Step 2, by computer drawing signal s to be analyzed (t) time-domain diagram.Signal waveform is a discontinuous band shape, is considered as time base flutter and causes the original signal high frequency noise that superposeed.
Step 3, to sampling original signal s (t) carry out population mean empirical mode decomposition (Ensemble Empirical Mode Decomposition, EEMD).The length of original signal s (t) is N, through EEMD, decomposes and obtains i intrinsic mode function (Intrinsic Mode Function, IMF) component c
j(t) and one remaining component r (t), j=1 wherein, 2 ..., i.C
j(t) instantaneous amplitude of domain space while being, the average tendency that r (t) is signal; In formula, j represents the sign subscript of IMF component, and i represents the number of intrinsic mode function IMF component, and t represents the sampling time;
Step 4, to obtaining obtaining i IMF component c in step 3
j(t) carry out respectively FFT conversion, by abscissa value get 10 times take 10 as end logarithm, ordinate value get 20 times take 10 as end logarithm, try to achieve each IMF component c
j(t) corresponding amplitude versus frequency characte C
j(w), wherein w represent to get 10 times take 10 after the end logarithm frequency;
Step 5, the amplitude versus frequency characte corresponding according to each IMF component, the C that maximum amplitude is greater than to 0
j(w) corresponding IMF component c
j(t) be considered as non-noise signal, the C that maximum amplitude is less than to 0
j(w) corresponding IMF component c
j(t) be considered as raw data composition and unmatched white noise completely in EEMD process.
Step 6, the raw data composition obtaining in step 5 and unmatched white noise are completely rejected from s (t), by the IMF component c of non-noise
j(t) carry out matching, reconstruct becomes new signal.Be the signal of removing after time base flutter.A kind of method of removing the time base flutter of high-speed sampling oscillograph according to claim 1, is characterized in that: the formula that in described step 3, EEMD decomposes is:
In formula:
S (t) is original signal;
C
j(t) be intrinsic mode function IMF;
R (t) is remaining component;
J is the sign subscript of IMF component;
T represents the sampling time.A kind of method of removing the time base flutter of high-speed sampling oscillograph according to claim 1, is characterized in that: in described step 6 to each IMF component c
j(t) carrying out fitting formula is: y (t)=∑
jc
j(t) in formula:
Y (t) is reconstruction signal;
C
j(t) be IMF component.
Described high-speed sampling oscillograph is that bandwidth is at the above high-speed sampling oscillograph of 10GHz.
The present invention can obtain following beneficial effect:
When sampling to signal source, high-speed sampling oscillograph there is time base flutter, sampled signal waveform is a discontinuous band shape, and after above-mentioned steps is processed, reconstruction signal waveform is a continuous curve, and retained original signal useful feature, the sampled result error obtaining obviously reduces.Remove a method for high-speed sampling oscillograph time base flutter, there is debounce successful, simple to operate, be easy to the advantage that program realizes.The method has good removal effect to time base flutter, more accurate to data processing after the sampling of high-speed sampling ripple device.
Accompanying drawing explanation
Fig. 1, a kind of method flow diagram of removing the time base flutter of high-speed sampling oscillograph
Fig. 2, system hardware block diagram
Fig. 3, time base debounce before waveform
Fig. 4, be respectively front 5 IMF components
Fig. 5, be respectively rear 5 IMF components
Fig. 6, be respectively front 5 amplitude-versus-frequency curves that IMF component is corresponding
Fig. 7, be respectively rear 5 amplitude-versus-frequency curves that IMF component is corresponding
Fig. 8, the waveform after matching again
Embodiment
Below in conjunction with accompanying drawing, preferred embodiment is elaborated.Should be emphasized that, following explanation is only exemplary, rather than in order to limit the scope of the invention and to apply.
Below in conjunction with accompanying drawing, provide most preferred embodiment of the present invention:
(1) with connector, Picosecond4016 signal source is connected with Aglient86100C broadband 70GHz sampling oscilloscope 86118A.By signal source, produced the negative step pulse signal that is nominally 5ps, high-speed sampling oscillograph carries out primary sample to signal source, and sampled data is imported to computing machine, and hardware connects as shown in Figure 2.
(2) by computer drawing signal time-domain diagram to be analyzed.Signal waveform is a discontinuous band shape, is considered as time base flutter and causes the original signal high frequency noise that superposeed, as shown in Figure 3.
(3) sampled data is carried out to population mean empirical mode decomposition (Ensemble Empirical Mode Decomposition, EEMD), sampled data length is 1350.Through EEMD, decompose and obtain 10 intrinsic mode functions (Intrinsic Mode Function, IMF) component c
j(t) and one remaining component r (t), j=1 wherein, 2 ..., 10.The instantaneous amplitude of domain space when each IMF component is, the average tendency that r (t) is signal.The formula that wherein EEMD decomposes is:
In formula:
S (t) is original signal;
C
j(t) be intrinsic mode function IMF;
R (t) is remaining component.
In EEMD process, add original signal standard deviation than the white noise that is 0.1, set sample number is 1000, decomposes and obtains as Fig. 4 waveform and Fig. 5 waveform.
(4) then to 10 IMF component c
j(t) carry out respectively FFT conversion, by abscissa value get 10 times take 10 as end logarithm, ordinate value get 20 times take 10 as end logarithm, try to achieve each IMF component c
j(t) corresponding amplitude versus frequency characte C
j(w), wherein w represent to get 10 times take 10 after the end logarithm frequency, obtain as Fig. 6 waveform and Fig. 7 waveform.
(5) amplitude versus frequency characte corresponding according to each IMF component again, the C that maximum amplitude is greater than to 0
j(w) corresponding IMF component c
j(t) be considered as non-noise signal, the C that maximum amplitude is less than to 0
j(w) corresponding IMF component c
j(t) be considered as the set of raw data composition and unmatched white noise completely in EEMD process, wherein the value of the corresponding j value of non-noise signal is 5,6 ..., 10.
(6) each IMF component c of non-noise will be considered as
j(t) carry out matching, reconstruct becomes new signal.Be the signal waveform of removing after time base flutter, as shown in Figure 8.Wherein fitting formula is:
In formula:
Y (t) is reconstruction signal;
C
j(t) be IMF component.
When sampling to signal source, high-speed sampling oscillograph there is time base flutter, sampled signal waveform is a discontinuous band shape, and after above-mentioned steps is processed, reconstruction signal waveform is a continuous curve, and retained original signal useful feature, the sampled result error obtaining obviously reduces.Remove a method for high-speed sampling oscillograph time base flutter, there is debounce successful, simple to operate, be easy to the advantage that program realizes.The method has good removal effect to time base flutter, more accurate to data processing after the sampling of high-speed sampling ripple device.
The above; be only the present invention's embodiment preferably, but protection scope of the present invention is not limited to this, is anyly familiar with in technical scope that those skilled in the art disclose in the present invention; the variation that can expect easily or replacement, within all should being encompassed in protection scope of the present invention.Therefore, protection scope of the present invention should be as the criterion with the protection domain of claim.
Claims (4)
1. a method of removing the time base flutter of high-speed sampling oscillograph, it is the data processing method of for high-speed sampling oscillograph, the time base flutter of signal source sampling being removed based on population mean empirical mode decomposition method EEMD, this method is in for data handling procedure, by EEMD method, obtain intrinsic mode function, filter out useful signal, noise signal is removed, and then the signal that makes new advances of reconstruct, reach the effect of removing time base flutter; It is characterized in that: its concrete steps are:
Step 1, with the signal that high-speed sampling oscillograph produces signal source, once sample and obtain original signal s (t), in s (t), have the data of N numerical value; The N obtaining numerical value is imported to computing machine and carry out data processing, wherein N is sampling number, and t is the sampling time;
Step 2, by computer drawing, treat the time-domain diagram of original signal s (t); Signal waveform is a discontinuous band shape, i.e. time base flutter causes the original signal high frequency noise that superposeed;
Step 3, sampling original signal s (t) is carried out to population mean empirical modal EEMD decompose, the length of original signal s (t) is N; Through EEMD, decompose and obtain i intrinsic mode function IMF component c
j(t) and one remaining component r (t), j=1 wherein, 2 ..., i; c
j(t) instantaneous amplitude of domain space while being, the average tendency of r (t) signal; In formula, j represents the sign subscript of IMF component, and i represents the number of intrinsic mode function IMF component, and t represents the sampling time;
Step 4, to obtaining obtaining i IMF component c in step 3
j(ε) carry out respectively FFT conversion, by abscissa value get 10 times take 10 as end logarithm, ordinate value get 20 times take 10 as end logarithm, try to achieve each IMF component c
j(t) corresponding amplitude versus frequency characte c
j(t), wherein w represent to get 10 times take 10 after the end logarithm frequency;
Step 5, the amplitude versus frequency characte corresponding according to each IMF component, the c that maximum amplitude is greater than to 0
j(w) corresponding IMF component c
j(t) be considered as non-noise signal, the c that maximum amplitude is less than to 0
j(w) corresponding IMF component c
j(t) be considered as raw data composition and unmatched white noise completely in EEMD process;
Step 6, the raw data composition obtaining in step 5 and unmatched white noise are completely rejected from s (t), by the IMF component c of non-noise
j(t) carry out matching, reconstruct becomes new signal; Be the signal of removing after time base flutter.
2. a kind of method of removing the time base flutter of high-speed sampling oscillograph according to claim 1, is characterized in that: the formula that in described step 3, EEMD decomposes is:
in formula:
S (t) is original signal;
C
j(t) be intrinsic mode function IMF;
R (t) is remaining component;
J is the sign subscript of IMF component;
T represents the sampling time.
3. a kind of method of removing the time base flutter of high-speed sampling oscillograph according to claim 1, is characterized in that: in described step 6 to each IMF component c
j(t) carrying out fitting formula is:
In formula:
Y (t) is reconstruction signal;
C
j(t) be IMF component.
4. a kind of method of removing the time base flutter of high-speed sampling oscillograph according to claim 1, is characterized in that: described high-speed sampling oscillograph is that bandwidth is oscillograph more than 10GHz.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210166086.6A CN102707252B (en) | 2012-05-24 | 2012-05-24 | Method for removing time base flutter of high-speed sampling oscilloscope |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210166086.6A CN102707252B (en) | 2012-05-24 | 2012-05-24 | Method for removing time base flutter of high-speed sampling oscilloscope |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102707252A CN102707252A (en) | 2012-10-03 |
CN102707252B true CN102707252B (en) | 2014-10-22 |
Family
ID=46900182
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210166086.6A Expired - Fee Related CN102707252B (en) | 2012-05-24 | 2012-05-24 | Method for removing time base flutter of high-speed sampling oscilloscope |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102707252B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104679698B (en) * | 2015-03-09 | 2017-10-31 | 电子科技大学 | The decomposition method of dither signal in high-speed serial communication |
CN105068131A (en) * | 2015-08-03 | 2015-11-18 | 中国科学院电子学研究所 | Aeromagnetic data leveling method |
CN107246891A (en) * | 2017-06-28 | 2017-10-13 | 中海艾普油气测试(天津)有限公司 | One kind test HTHP underground data signal extraction method |
CN109581018B (en) * | 2019-01-08 | 2022-02-01 | 北京工业大学 | Time base error synchronous compensation method and device of broadband sampling oscilloscope |
CN113447873B (en) * | 2021-07-12 | 2022-07-19 | 北京无线电计量测试研究所 | Sampling oscilloscope complex frequency response calibration device and method |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100465944C (en) * | 2007-04-13 | 2009-03-04 | 北京工业大学 | Time-base dither method for compensated oscilloscope |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7941298B2 (en) * | 2006-09-07 | 2011-05-10 | DynaDx Corporation | Noise-assisted data analysis method, system and program product therefor |
TWI439874B (en) * | 2009-01-10 | 2014-06-01 | Ind Tech Res Inst | Signal processing method for hierarchical empirical mode decomposition and apparatus therefor |
-
2012
- 2012-05-24 CN CN201210166086.6A patent/CN102707252B/en not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100465944C (en) * | 2007-04-13 | 2009-03-04 | 北京工业大学 | Time-base dither method for compensated oscilloscope |
Non-Patent Citations (6)
Title |
---|
System Health Management Conference》.2010,正文第1-7页. * |
Wei Guo等.Enhancing the Ability of Ensemble Empirical Mode Decomposition in Machine Fault Diagnosis.《2010Progno stics& * |
Wei Guo等.Enhancing the Ability of Ensemble Empirical Mode Decomposition in Machine Fault Diagnosis.《2010Progno stics&System Health Management Conference》.2010,正文第1-7页. |
孙圣和等.消除等效取样时基抖动影响的反卷积方法.《计量学报》.1989,第10卷(第2期),第91-96页. * |
陈东方.采用EMD方法消除瞬态散射回波中的高斯白噪声干扰.《电子学报》.2004,第32卷(第3期),第496-498页. * |
龚鹏伟等.取样示波器时基抖动的修正.《宇航计测技术》.2011,第31卷(第1期),第17-20页. * |
Also Published As
Publication number | Publication date |
---|---|
CN102707252A (en) | 2012-10-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102707252B (en) | Method for removing time base flutter of high-speed sampling oscilloscope | |
Peng et al. | Sparse signal decomposition method based on multi-scale chirplet and its application to the fault diagnosis of gearboxes | |
US20220343898A1 (en) | Speech recognition method and apparatus, and computer-readable storage medium | |
EP2209117A1 (en) | Method for determining unbiased signal amplitude estimates after cepstral variance modification | |
Xie et al. | Fast-varying AM–FM components extraction based on an adaptive STFT | |
CN104215459A (en) | Bearing fault diagnosis method | |
CN112595782B (en) | Ultrasonic transverse wave take-off point identification method and system based on EEMD algorithm | |
Melgoza et al. | Comparing radar receiver pulse deinterleaving performance of differing window functions for bandpass FIR filter design | |
Zhidong et al. | A new method for processing end effect in empirical mode decomposition | |
CN105280178A (en) | audio signal processing device and audio signal processing method thereof | |
Faisal et al. | Suppression of false-terms in wigner-ville distribution using time and frequency windowing | |
EP2690448A3 (en) | System for improving probability of transient event detection | |
Liu et al. | A method for blind source separation of multichannel electromagnetic radiation in the field | |
Gao et al. | A Method Using EEMD and L-Kurtosis to detect faults in roller bearings | |
JP2010002394A (en) | Feeble signal analytical apparatus, feeble signal analytical method, and feeble signal analytical program | |
Zhang et al. | Application of Hilbert-Huang transform to extract arrival time of ultrasonic lamb waves | |
Belega et al. | Choice of the window used in the interpolated discrete Fourier transform method | |
CN117330939B (en) | ATE-based switched capacitor filter adjustment and measurement method, system and storage medium | |
US20240019470A1 (en) | Detector circuit, signal processing circuit, and measurement instrument | |
Anqi et al. | Random vibration measured signal induction method research | |
US20230408551A1 (en) | Oscilloscope having a principal component analyzer | |
Liu et al. | Rotor local rubbing fault feature analysis based on Teager-Huang transform | |
Shankar et al. | EMD based Hilbert energy spectrum analysis of adventitious lung wave signals | |
CN116165486A (en) | Method and system for recovering time domain waveform of partial discharge pulse electric field | |
CN205454097U (en) | Audio frequency test analytic system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
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: 20141022 Termination date: 20170524 |