CN102879681A - Flickering real-time calculation method - Google Patents
Flickering real-time calculation method Download PDFInfo
- Publication number
- CN102879681A CN102879681A CN2012103635816A CN201210363581A CN102879681A CN 102879681 A CN102879681 A CN 102879681A CN 2012103635816 A CN2012103635816 A CN 2012103635816A CN 201210363581 A CN201210363581 A CN 201210363581A CN 102879681 A CN102879681 A CN 102879681A
- Authority
- CN
- China
- Prior art keywords
- flickering
- array
- measured data
- real
- integers
- 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
- Testing, Inspecting, Measuring Of Stereoscopic Televisions And Televisions (AREA)
- Complex Calculations (AREA)
Abstract
The invention discloses a flickering real-time calculation method which is used in flickering measurement and calculation of an electric energy quality on-line monitoring device. The flickering real-time calculation method comprises the following steps of: measuring flickering transient visibility data, and setting the data to be floating point numbers which are more than zero to form a measured data floating point number set; multiplying each floating point number in the measured data floating point number set with a coefficient to obtain a group of measured data integer sets; setting a group of statistical numbers; zeroing every element in the statistical number group, traversing the measured data integer sets, and once an integer is obtained, adding 1 to corresponding variables in the statistical number group; and after a traversing process is finished, obtaining the probability of the flickering transient visibility data. By the quick flickering calculation method, real-time calculation of flickering measurement can be finished; flickering indexes can be mastered in real time; and the power supply quality of a power grid can be monitored in real time.
Description
Technical field
The present invention relates to a kind of electric energy quality on-line monitoring device flickering computing method.
Background technology
Flickering is an important technical indicator in the electric energy quality monitoring, and it is a kind of phenomenon that voltage fluctuation causes, but therefore reflecting voltage fluctuation situation, gets more and more for the measurement of flickering, the research of calculating.Probability distribution for the instantaneous visual sense degree in the flickering calculating is calculated, and method commonly used is piecewise linear interpolation method and sort method method at present.Flickering measurement method according to the IEC recommendation, calculated an instantaneous flicker in per 10 minutes, if the sample frequency according to 800Hz, the instantaneous visual sense degrees of data of accumulation will reach 480000 points in 10 minutes, using above conventional method is to be difficult to finish at short notice calculating, the real-time grasp to the flickering index can't be satisfied, the real-time quality monitoring to grid supply quality can not be satisfied.
Summary of the invention
In order to solve the real time problems of the calculating of instantaneous visual sense degrees of data in the existing flickering computing technique, the invention provides a kind of accuracy that both satisfied, satisfied again the flickering computing method of real-time.
In order to achieve the above object, the technical solution used in the present invention is:
A kind of flickering real-time computing technique is characterized in that may further comprise the steps:
(1), measure the instantaneous visual sense degrees of data of flickering, the floating number that these data are greater than zero is referred to as measured data floating-point manifold.
(2), each floating number that this measured data floating number is concentrated multiply by a coefficient, obtains one group of set of integers and is referred to as the measured data set of integers, and the probability statistics of measured data floating-point manifold in fact just are equal to probability statistics to the measured data set of integers.
(3), set one group of integer array and be referred to as to add up array.
(4), when the measured data set of integers is carried out probability statistics, will add up first each the variable zero clearing in the array, then travel through the measured data set of integers, and whenever obtain a round values, the relevant variable of then adding up in the array adds 1.
(5), ergodic process finishes, the numerical value of then adding up each variable of array has just represented the number of times that respective integer occurs in the measured data set of integers in whole measured data set of integers, obtain the probability of the instantaneous visual sense degrees of data of flickering.
Aforesaid a kind of flickering real-time computing technique is characterized in that may further comprise the steps:
(1), establishing the data that need statistics is greater than zero floating number { X
nBe referred to as measured data floating-point manifold, and n={0 wherein, 1,2.。。, N}, N are the data point number that needs statistics, determine { X
nMeasuring accuracy namely need the scale M that keeps.
(2), with measured data floating-point manifold { X
nMultiply by 10
MBe transformed into integer { Y
nBe referred to as the measured data set of integers, to { X
nProbability statistics in fact be exactly to { Y
nProbability statistics.
(3), to { Y
nWhen carrying out probability statistics, estimate first { Y
nMaximal value Y
Max, according to this maximal value Y
Max, set an integer array sort[count] and be called the statistics array, add up the big or small count of array〉Y
Max
(4), at every turn add up beginning, will add up array sort[count first] zero clearing, then from n=0 to n=N, travel through { Y
n, whenever obtain a Y
n, then add up accordingly the variable sort[Y in the array
n] add 1.
(5), this process finishes, and then adds up array sort[count] in each variable sort[y] (0≤y≤Y
Max) just represented numerical value y at { Y
nThe middle number of times m that occurs, m/N is that numerical value y is at { Y just
nIn the probability that occurs, such as y at { X
nIn corresponding value be x, then m/N namely x at { X
nThe middle probability that occurs.
The invention has the beneficial effects as follows:
The present invention can finish the calculating of flickering data in a short period of time, and is not subjected to the impact of statistics element number, can realize that the real-time statistics of large capacity flicker measurement data is calculated.Can satisfy the real-time grasp to the flickering index, satisfy the real-time quality monitoring to grid supply quality.
Embodiment
Below embodiments of the invention are done concrete introduction:
(1), establishes x (n)=f (n), 0≤n≤N, 0≤x (n)≤X
Max
X (n) is a discrete floating number ordered series of numbers, represents the results of intermediate calculations in the detecting instrument testing process, whenever obtains N point and will do a probability statistics calculating to x (n).
Suppose the maximal value X of this sequence
Max＜99, the scale that needs to keep is 2.
(2), establish y (n)=x (n) * 10
2, 0≤n≤N.Y (n) is defined as set of integers.
0≤y (n)≤9999 then.So just obtained with x (n) deserved set of integers y (n).
(3), set one group of integer array sort[count], count=10000.
(4), in the actual computation, install every once sampling and just calculate a some x (n).
When n=0, sort[count] the array zero clearing.
y(0)＝x(0)*10
2，sort[y(0)]=sort[y(0)]+1；
Whenever obtain later on an x (n), just obtain a y (n), corresponding sort[y (n)]=sort[y (n)]+1.
When n=N-1, carry out probability statistics.
The probability that occurs in x (n) as obtaining floating number x becomes x into integer first, gets y=x*10
2
Then this probable value is: sort[y]/N.
Equally, can carry out probability statistics to other numerical value.
Above demonstration and described ultimate principle of the present invention, principal character and advantage.The technician of the industry should understand; the present invention is not restricted to the described embodiments; that describes in above-described embodiment and the instructions just illustrates principle of the present invention; without departing from the spirit and scope of the present invention; the present invention also has various changes and modifications, and these changes and improvements all fall in the claimed scope of the invention.The claimed scope of the present invention is by appending claims and equivalent circle thereof.
Claims (2)
1. a flickering real-time computing technique is characterized in that, may further comprise the steps:
(1), measure the instantaneous visual sense degrees of data of flickering, these data are made as floating number greater than zero, be called measured data floating-point manifold;
(2), each floating number that this measured data floating number is concentrated multiply by a coefficient, obtains one group of set of integers, is referred to as the measured data set of integers;
(3), set one group of integer array and be referred to as to add up array;
(4), will add up each the element zero clearing in the array, then travel through the measured data set of integers, whenever obtain a round values, the relevant variable of then adding up in the array adds 1;
(5), above-mentioned ergodic process finishes, the numerical value of then adding up each variable of array has just represented the number of times that respective integer occurs in the measured data set of integers in whole measured data set of integers, obtain the probability of the instantaneous visual sense degrees of data of flickering.
2. a kind of flickering real-time computing technique according to claim 1 is characterized in that, may further comprise the steps:
(1), establishing the instantaneous visual sense degrees of data of flickering is greater than zero floating number { X
nBe referred to as measured data floating-point manifold, and n={0 wherein, 1,2 ..., N}, N are the data point number that needs statistics, determine { X
nMeasuring accuracy namely need the scale M that keeps;
(2), with measured data floating-point manifold { X
nMultiply by 10
MBe transformed into integer { Y
nBe referred to as the measured data set of integers;
(3), to { Y
nWhen carrying out probability statistics, estimate first { Y
nMaximal value Y
Max, according to this maximal value Y
Max, set one group of integer array sort[count], be called the statistics array, the big or small count of statistics array〉Y
Max
(4), at every turn add up beginning, will add up array sort[count first] zero clearing, then from n=0 to n=N, travel through { Y
n, whenever obtain a Y
n, then add up accordingly the variable sort[Y in the array
n] add 1;
(5), said process finishes, and then adds up array sort[count] in each variable sort[y] (0≤y≤Y
Max) just represented numerical value y at { Y
nThe middle number of times m that occurs, m/N is that numerical value y is at { Y just
nIn the probability that occurs, such as y at { X
nIn corresponding value be x, then m/N namely x at { X
nThe middle probability that occurs.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210363581.6A CN102879681B (en) | 2012-09-26 | 2012-09-26 | A kind of flickering real-time computing technique |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210363581.6A CN102879681B (en) | 2012-09-26 | 2012-09-26 | A kind of flickering real-time computing technique |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102879681A true CN102879681A (en) | 2013-01-16 |
| CN102879681B CN102879681B (en) | 2015-08-26 |
Family
ID=47481077
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210363581.6A Active CN102879681B (en) | 2012-09-26 | 2012-09-26 | A kind of flickering real-time computing technique |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102879681B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103543317A (en) * | 2013-09-26 | 2014-01-29 | 光一科技股份有限公司 | Method for measuring IEC (International Electrical Commission) standard voltage flicker |
| CN103713223A (en) * | 2014-01-03 | 2014-04-09 | 烟台东方威思顿电气有限公司 | Method for flickering measuring of self-adaption range low storage data volume |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060044933A1 (en) * | 2000-08-25 | 2006-03-02 | Micro Technology, Inc. | Burst read addressing in a non-volatile memory device |
| CN102288807A (en) * | 2011-05-14 | 2011-12-21 | 苏州大学 | Method for measuring electric network voltage flicker |
| CN102565879A (en) * | 2011-08-04 | 2012-07-11 | 南京信息工程大学 | Processing method of wind data |
-
2012
- 2012-09-26 CN CN201210363581.6A patent/CN102879681B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060044933A1 (en) * | 2000-08-25 | 2006-03-02 | Micro Technology, Inc. | Burst read addressing in a non-volatile memory device |
| CN102288807A (en) * | 2011-05-14 | 2011-12-21 | 苏州大学 | Method for measuring electric network voltage flicker |
| CN102565879A (en) * | 2011-08-04 | 2012-07-11 | 南京信息工程大学 | Processing method of wind data |
Non-Patent Citations (4)
| Title |
|---|
| 余健明等: "电压闪变的评价与工程计算方法研究", 《西安理工大学学报》, vol. 13, no. 01, 31 January 1997 (1997-01-31) * |
| 刘诚等: "交流电弧炉的电压闪变及计算", 《工业加热》, vol. 41, no. 04, 30 April 2012 (2012-04-30) * |
| 史三省等: "基于FFT的电压波动与闪变测量算法", 《电力***及其自动化学报》, vol. 22, no. 06, 31 December 2010 (2010-12-31) * |
| 郭晓丽等: "短时间闪变严重度Pst的计算方法的研究", 《南通工学院学报（自然科学版）》, vol. 03, no. 03, 30 September 2004 (2004-09-30) * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103543317A (en) * | 2013-09-26 | 2014-01-29 | 光一科技股份有限公司 | Method for measuring IEC (International Electrical Commission) standard voltage flicker |
| CN103543317B (en) * | 2013-09-26 | 2016-06-01 | 光一科技股份有限公司 | For measuring the method for IEC standard voltage flickering |
| CN103713223A (en) * | 2014-01-03 | 2014-04-09 | 烟台东方威思顿电气有限公司 | Method for flickering measuring of self-adaption range low storage data volume |
| CN103713223B (en) * | 2014-01-03 | 2016-01-27 | 烟台东方威思顿电气股份有限公司 | A kind of low storage data volume flickering measurement method of self-adaptation range |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102879681B (en) | 2015-08-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106443285B (en) | Multiple-harmonic-source harmonic responsibility quantitative analysis method based on total least square method | |
| CN102662106B (en) | Method for measuring electric energy of harmonic power network | |
| CN103996084B (en) | Wind power probability Forecasting Methodology based on longitudinal moment Markov chain model | |
| CN103116064A (en) | Method and device for detecting voltage fluctuation and flicker based on energy operator and spectrum correction | |
| CN109284933B (en) | Electronic transformer state evaluation system and method based on mathematical statistics | |
| CN103728570B (en) | Battery-thermal-characteristic-based health state detection method | |
| CN101408577B (en) | Analysis method of low-frequency oscillation node contribution factor based on wide field measuring information | |
| WO2015061227A1 (en) | Performance analysis of power grid monitors | |
| CN106383501A (en) | Carbon emission process monitoring and control method | |
| CN106526312A (en) | Electric energy metering method based on R-V(III) window FFT double-peak interpolation | |
| CN105093166B (en) | A kind of field test method of electronic electric energy meter | |
| CN108490257A (en) | A kind of short time-window inter-harmonic wave measuring method based on spectrum curve fitting | |
| CN103472373B (en) | Method for evaluating external insulation equivalence of natural pollution insulator and artificial pollution insulator | |
| CN106405280A (en) | Intelligent transformer station online monitoring parameter trend early warning method | |
| CN105486945A (en) | Determination method for line loss abnormity of 10kV line | |
| CN105093163A (en) | Error debugging method for electronic electric energy meter | |
| CN102879681B (en) | A kind of flickering real-time computing technique | |
| CN109359822A (en) | Electronic type voltage transformer measuring state appraisal procedure and system | |
| CN108415880A (en) | A kind of line loss characteristic analysis method based on Sample Entropy and wavelet transformation | |
| CN107679017A (en) | A kind of self-healing cable status evaluation method and device | |
| CN107101800A (en) | Wind power measuring equipment and method for power transmission iron tower model | |
| CN103218755A (en) | Micro-grid evaluating method with inverse non-extensive entropy adopted | |
| Germer | High-precision ac measurements using the Monte Carlo method | |
| CN108969153A (en) | A kind of knee joint structural equation measuring device and its measurement method | |
| CN105158555A (en) | Method for measuring and verifying saved electric power and electric quantity in boosting operation of power distribution line |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20130116 Assignee: NANJING GUODIAN NANZI RURAL POWER GRID AUTOMATION ENGINEERING CO., LTD. Assignor: Nanjing Automation Co., Ltd., China Electronics Corp. Contract record no.: 2013320000276 Denomination of invention: Flickering real-time calculation method License type: Exclusive License Record date: 20130401 |
|
| LICC | Enforcement, change and cancellation of record of contracts on the licence for exploitation of a patent or utility model | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| EC01 | Cancellation of recordation of patent licensing contract |
Assignee: NANJING GUODIAN NANZI RURAL POWER GRID AUTOMATION ENGINEERING CO., LTD. Assignor: Nanjing Automation Co., Ltd., China Electronics Corp. Contract record no.: 2013320000276 Date of cancellation: 20181016 |
|
| EC01 | Cancellation of recordation of patent licensing contract |