CN109900952A - A kind of transient signal rapidly extracting processing method - Google Patents
A kind of transient signal rapidly extracting processing method Download PDFInfo
- Publication number
- CN109900952A CN109900952A CN201910188507.7A CN201910188507A CN109900952A CN 109900952 A CN109900952 A CN 109900952A CN 201910188507 A CN201910188507 A CN 201910188507A CN 109900952 A CN109900952 A CN 109900952A
- Authority
- CN
- China
- Prior art keywords
- data
- uproar
- mean filter
- filter
- transient 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
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Data Exchanges In Wide-Area Networks (AREA)
Abstract
The invention discloses a kind of transient signal rapidly extracting processing method, data acquisition unit first carries out continuous acquisition to signal, and data processing unit carries out subpackage reading, and calculates the first bag data bottom and make an uproar.Then data processing unit carries out high-order mean filter and low order mean filter to data respectively, then subtracts low order mean filter data with high-order mean filter data to obtain final filter result.It is finally made an uproar and is compared the bottom of with point by point with final filter result, if making an uproar the bottom of less than or equal to, be judged as no triggering, continue data acquisition;It makes an uproar the bottom of if more than, is then judged as triggering, and the data of triggering front and back certain period of time are subjected to packing storage, completes transient signal and extract.The present invention can be improved compared with transient signal extraction rate under the conditions of high sampling rate, reduce data processing algorithm step, improves transient signal and extracts handling capacity.
Description
Technical field
The invention belongs to digital processing field more particularly to a kind of transient signal rapidly extracting processing methods.
Background technique
Direction with power grid towards high voltage, Gao Zhineng, high reliability is developed, to transient overvoltage measuring device performance
Requirement it is higher and higher, power grid transient overvoltage not only in decision systems power equipment dielectric level, also threaten installations and facilities
Safe operation, be likely to result in insulation breakdown under overvoltage effect containing defective equipment, or even cause explosion accident.This
So that being extremely important to the measurement of power grid transient overvoltage and analysis.
Transient signal, which extracts, at present mainly filters out power frequency component and high-frequency noise by bandpass filter, then carries out transient state letter
It number extracts, but this method needs specialized hardware to carry out Digital Signal Processing, it can not be accurate based on PC and data collecting card framework hardware
Carry out the extraction of high sampling rate transient signal.
Summary of the invention
To solve the above-mentioned problems, the present invention proposes that a kind of transient signal rapid extracting method method, first data acquire
Unit carries out continuous acquisition to signal, and data processing unit carries out subpackage reading, and calculates the first bag data bottom and make an uproar.Then data
Processing unit carries out high-order mean filter and low order mean filter to data respectively, then subtracts low order with high-order mean filter data
Mean filter data are to obtain final filter result.Finally made an uproar and be compared the bottom of with point by point with final filter result, if be less than or
It makes an uproar the bottom of equal to, is then judged as no triggering, continue data acquisition;It makes an uproar the bottom of if more than, is then judged as triggering, and will triggering front and back
The data of certain period of time carry out packing storage, complete transient signal and extract.
Further, in data acquisition phase: the data that data collecting card DMA is uploaded are read in data processing unit subpackage,
Every bag data is total to N number of sampled point.
Further, the calculation method that bottom is made an uproar are as follows: the maximum sampled value of the first bag data subtracts minimum sampled value, i.e.,
Vmax=max (Vi) i ∈ [1, N] (formula 1);
Vmin=min (Vi) i ∈ [1, N] (formula 2);
VIt makes an uproar at bottom=Vmax-Vmin(formula 3);
Wherein, VmaxFor the first bag data maximum value, VminFor the first bag data minimum value, VIt makes an uproar at bottomIt makes an uproar the bottom of for.
Further, in filtering calculation stages:
(1) high-order mean filter carries out m sampled point mean filter, and low order mean filter carries out n sampled point mean value filter
Wave, wherein m > > n;
(2) calculate the cumulative of high-order mean filter and:
(3) the initial filter value of high-order mean filter: V is calculatedH_Init=SumH_Init/ m (formula 5);
(4) calculate the cumulative of low order mean filter and:
(5) the initial filter value of low order mean filter: V is calculatedL_Init=SumL_Init/ n (formula 7);
(6) difference of filter result twice: V is calculatedresult=VH_Init-VL_Init(formula 8).
Further, the filtering algorithm of the m/2+1 point is as follows:
The beneficial effects of the present invention are: the present invention can complete the rapidly extracting to transient voltage signal, improve transient state letter
Number extract sensitivity and handling capacity, reduce data processing amount, reduce false triggering rate.
Detailed description of the invention
Fig. 1 is flow chart of the invention;
Fig. 2 is the original waveform figure before filtering;
Fig. 3 is the waveform diagram of high-grade filting result;
Fig. 4 is the waveform diagram of low order filter result;
Fig. 5 is the waveform diagram of final filter result.
Specific embodiment
For a clearer understanding of the technical characteristics, objects and effects of the present invention, this hair of Detailed description of the invention is now compareed
Bright specific embodiment.
The present invention proposes a kind of transient signal rapidly extracting processing method, as shown in Figure 1, including that data acquisition, bottom are made an uproar meter
Calculation, filtering calculating, triggering, which are compared, stores this 5 steps with data, specific as follows:
S1. data acquire: data acquisition unit carries out continuous acquisition to signal, and data processing unit carries out subpackage reading;
S2. bottom is made an uproar calculating: being made an uproar at the bottom that data processing unit calculates the first bag data;
S3. filtering calculates: carrying out high-order mean filter and low order mean filter respectively to data, then with high-order mean filter
Data subtract low order mean filter data to obtain final filter result;
S4. triggering is compared: it is made an uproar and is compared the bottom of with point by point with final filter result, made an uproar the bottom of if more than, be then judged as triggering,
Execute step S5;If making an uproar the bottom of less than or equal to, it is judged as no triggering, continues data acquisition, execute step S1;
S5. data store: the data of triggering front and back certain period of time being carried out packing storage, transient signal is completed and extracts.
In a specific embodiment of the present invention, data collecting card carries out 40MS/s continuous sampling, and every bag data is set as 4M
Sampled point is uploaded to the end PC by DMA, and original waveform is as shown in Figure 2.PC termination receive data after, calculate bottom make an uproar benchmark electricity
It is flat, the low order mean filter of 10001 sampled point high-order mean filters and 3 points is then carried out respectively, and filtering starts Shi Xianji
It calculates preceding 10001 point sampling values and while divided by 10001, obtains the high-order mean filter of the 5000th point as a result, such as Fig. 3 institute
Show.Then calculate 4999,5000,5001 point sampling values and, while divided by 3, obtain the low order mean filter knot of the 5000th point
Fruit, as shown in figure 4, high-order mean filter result, which subtracts low order mean filter result, can be obtained the 5000th point of filter result.
When carrying out the 5001st filter result calculating, it is only necessary to by preceding 10001 points and the 1st point is subtracted, in addition the 10002nd point, so
High-grade filting can be obtained divided by 10001 afterwards as a result, preceding 4999,5000,5001 and subtracting at 4999 points plus 5002, then divided by
3 can be obtained low order filter result, and filter result can be obtained by then subtracting each other.The data processing of a point is directed in the above method
3 subtractions are only needed, transient signal extraction can be completed in 2 sub-additions and 2 divisions.Its final filter effect is as shown in Figure 5.
The above is only a preferred embodiment of the present invention, it should be understood that the present invention is not limited to described herein
Form should not be regarded as an exclusion of other examples, and can be used for other combinations, modifications, and environments, and can be at this
In the text contemplated scope, modifications can be made through the above teachings or related fields of technology or knowledge.And those skilled in the art institute into
Capable modifications and changes do not depart from the spirit and scope of the present invention, then all should be in the protection scope of appended claims of the present invention
It is interior.
Claims (5)
1. a kind of transient signal rapidly extracting processing method, which comprises the following steps:
S1. data acquire: data acquisition unit carries out continuous acquisition to signal, and data processing unit carries out subpackage reading;
S2. bottom is made an uproar calculating: being made an uproar at the bottom that data processing unit calculates the first bag data;
S3. filtering calculates: carrying out high-order mean filter and low order mean filter respectively to data, then with high-order mean filter data
Low order mean filter data are subtracted to obtain final filter result;
S4. triggering is compared: being made an uproar and is compared the bottom of with point by point with final filter result, made an uproar the bottom of if more than, be then judged as triggering, executed
Step S5;If making an uproar the bottom of less than or equal to, it is judged as no triggering, continues data acquisition, execute step S1;
S5. data store: the data of triggering front and back certain period of time being carried out packing storage, transient signal is completed and extracts.
2. a kind of transient signal rapidly extracting processing method according to claim 1, which is characterized in that the step S1 tool
Body are as follows: the data that data collecting card DMA is uploaded are read in data processing unit subpackage, and every bag data is total to N number of sampled point.
3. a kind of transient signal rapidly extracting processing method according to claim 1, which is characterized in that the step S2
In, calculation method that bottom is made an uproar are as follows: the maximum sampled value of the first bag data subtracts minimum sampled value, i.e.,
Vmax=max (Vi) i ∈ [1, N] (formula 1);
Vmin=min (Vi) i ∈ [1, N] (formula 2);
VIt makes an uproar at bottom=Vmax-Vmin(formula 3);
Wherein, VmaxFor the first bag data maximum value, VminFor the first bag data minimum value, VIt makes an uproar at bottomIt makes an uproar the bottom of for.
4. a kind of transient signal rapidly extracting processing method according to claim 1, which is characterized in that the step S3 packet
Include following sub-step:
S31. high-order mean filter carries out m sampled point mean filter, and low order mean filter carries out n sampled point mean filter,
Wherein m > > n;
S32. calculate the cumulative of high-order mean filter and:
S33. the initial filter value of high-order mean filter is calculated:
VH_Init=SumH_Init/ m (formula 5);
S34. calculate the cumulative of low order mean filter and:
S35. the initial filter value of low order mean filter is calculated:
VL_Init=SumL_Init/ n (formula 7);
S36. the difference of filter result twice is calculated:
Vresult=VH_Init-VL_Init(formula 8).
5. a kind of transient signal rapidly extracting processing method according to claim 4, which is characterized in that step S32-S36
What is completed is that the m/2 point filtering calculates, and the filtering algorithm of m/2+1 point is as follows:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910188507.7A CN109900952B (en) | 2019-03-13 | 2019-03-13 | Transient signal rapid extraction processing method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910188507.7A CN109900952B (en) | 2019-03-13 | 2019-03-13 | Transient signal rapid extraction processing method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109900952A true CN109900952A (en) | 2019-06-18 |
CN109900952B CN109900952B (en) | 2020-11-24 |
Family
ID=66952159
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910188507.7A Active CN109900952B (en) | 2019-03-13 | 2019-03-13 | Transient signal rapid extraction processing method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109900952B (en) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101832756A (en) * | 2009-03-10 | 2010-09-15 | 深圳迈瑞生物医疗电子股份有限公司 | Method and device for measuring displacement of targets in images and carrying out strain and strain rate imaging |
CN102497220A (en) * | 2011-12-13 | 2012-06-13 | 重庆大学 | Calling system based on flexible alternating current signal technology and communication method for calling system |
CN103493382A (en) * | 2012-03-09 | 2014-01-01 | 华为技术有限公司 | Method, device, apparatus and system for cancelling multicarrier transmitting interference |
EP3035530A1 (en) * | 2014-12-15 | 2016-06-22 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Circuit for comparing a voltage to a threshold |
CN106549651A (en) * | 2016-09-14 | 2017-03-29 | 芯海科技(深圳)股份有限公司 | A kind of highly-precise filtering method of quick foundation |
CN106685435A (en) * | 2016-12-14 | 2017-05-17 | 芯海科技(深圳)股份有限公司 | Method for improving effective-signal variation rapidly in low-SNR (Signal to Noise Ratio) signals |
CN107526032A (en) * | 2017-09-19 | 2017-12-29 | 海盐众信电子有限公司 | A kind of magnetic for being used to produce relay keeps machine examination machine and detection method |
CN108429228A (en) * | 2018-04-02 | 2018-08-21 | 中铁第四勘察设计院集团有限公司 | Split-phase Over-voltage Analysis method is crossed in a kind of electric railway |
CN109067393A (en) * | 2018-08-28 | 2018-12-21 | 南方电网科学研究院有限责任公司 | Phase locking method, device and equipment of power system |
CN109239447A (en) * | 2018-08-21 | 2019-01-18 | 上海赛璞乐电力科技有限公司 | A kind of power distribution network transient overvoltage monitoring system |
-
2019
- 2019-03-13 CN CN201910188507.7A patent/CN109900952B/en active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101832756A (en) * | 2009-03-10 | 2010-09-15 | 深圳迈瑞生物医疗电子股份有限公司 | Method and device for measuring displacement of targets in images and carrying out strain and strain rate imaging |
CN102497220A (en) * | 2011-12-13 | 2012-06-13 | 重庆大学 | Calling system based on flexible alternating current signal technology and communication method for calling system |
CN103493382A (en) * | 2012-03-09 | 2014-01-01 | 华为技术有限公司 | Method, device, apparatus and system for cancelling multicarrier transmitting interference |
EP3035530A1 (en) * | 2014-12-15 | 2016-06-22 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Circuit for comparing a voltage to a threshold |
CN106549651A (en) * | 2016-09-14 | 2017-03-29 | 芯海科技(深圳)股份有限公司 | A kind of highly-precise filtering method of quick foundation |
CN106685435A (en) * | 2016-12-14 | 2017-05-17 | 芯海科技(深圳)股份有限公司 | Method for improving effective-signal variation rapidly in low-SNR (Signal to Noise Ratio) signals |
CN107526032A (en) * | 2017-09-19 | 2017-12-29 | 海盐众信电子有限公司 | A kind of magnetic for being used to produce relay keeps machine examination machine and detection method |
CN108429228A (en) * | 2018-04-02 | 2018-08-21 | 中铁第四勘察设计院集团有限公司 | Split-phase Over-voltage Analysis method is crossed in a kind of electric railway |
CN109239447A (en) * | 2018-08-21 | 2019-01-18 | 上海赛璞乐电力科技有限公司 | A kind of power distribution network transient overvoltage monitoring system |
CN109067393A (en) * | 2018-08-28 | 2018-12-21 | 南方电网科学研究院有限责任公司 | Phase locking method, device and equipment of power system |
Non-Patent Citations (1)
Title |
---|
冯宇: "电流互感器暂态特性试验用标准信号的参数估计", 《电网技术》 * |
Also Published As
Publication number | Publication date |
---|---|
CN109900952B (en) | 2020-11-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lopes et al. | A traveling-wave detection method based on Park's transformation for fault locators | |
CN104977502B (en) | A kind of extra high voltage direct current transmission line internal fault external fault recognition methods | |
CN106066437B (en) | A kind of on-line monitoring device of metallic oxide arrester and its monitoring method | |
CN101373197B (en) | Method and apparatus for on-line monitoring over voltage | |
CN108805107B (en) | Method for identifying partial discharge defects in GIS based on PRPS signal | |
CN203133168U (en) | Power harmonic detector | |
CN109298287A (en) | Small current neutral grounding Distribution Network Failure kind identification method based on fault detector data | |
Dashtdar et al. | Fault location in the transmission network using a discrete wavelet transform | |
CN107561393A (en) | A kind of arrester early defect live testing system and method based on total current harmonic wave vector | |
CN107395157A (en) | Grounded screen potential difference filtering method based on wavelet transformation and weighted moving average | |
CN106680585B (en) | Harmonic wave/m-Acetyl chlorophosphonazo detection method | |
Bains et al. | Enhanced phasor estimation technique for fault location in series-compensated lines | |
CN111046791A (en) | Current signal filtering and denoising method based on generalized S transform containing variable factors | |
Han et al. | Fault type classification in transmission line using STFT | |
CN103091561B (en) | Device obtaining direct current signals from alternative current and direct current superposition signals and method thereof | |
Yusuff et al. | Stationary wavelet transform and single differentiator based decaying DC-offset filtering in post fault measurements | |
CN109900952A (en) | A kind of transient signal rapidly extracting processing method | |
CN109142867A (en) | Based on harmonic detecting method, the equipment for improving sliding window discrete Fourier transform | |
CN108872747B (en) | Surge protector resistive current extraction device and method based on correlation coefficient | |
CN107064634B (en) | The detection method of Harmonious Waves in Power Systems | |
CN116125208A (en) | Power distribution network single-phase earth fault positioning method based on data acquisition equipment | |
CN207408505U (en) | A kind of arrester early defect live testing system based on total current harmonic wave vector | |
Zou et al. | Mathematical morphology based phase selection scheme in digital relaying | |
CN106058810B (en) | A kind of excitation flow recognition method based on power failure component criterion | |
CN213986638U (en) | High-frequency grounding impedance test analyzer based on digital filtering technology |
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 |