Nothing Special   »   [go: up one dir, main page]

CN110069984A - One-variable linear regression frequency measuring method, electronic equipment and storage medium - Google Patents

One-variable linear regression frequency measuring method, electronic equipment and storage medium Download PDF

Info

Publication number
CN110069984A
CN110069984A CN201910179979.6A CN201910179979A CN110069984A CN 110069984 A CN110069984 A CN 110069984A CN 201910179979 A CN201910179979 A CN 201910179979A CN 110069984 A CN110069984 A CN 110069984A
Authority
CN
China
Prior art keywords
sinusoidal signal
dead
best estimate
time
center position
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.)
Pending
Application number
CN201910179979.6A
Other languages
Chinese (zh)
Inventor
魏强
何立林
罗雄宇
陈世奎
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhuhai Xingnuo Energy Technology Co Ltd
Original Assignee
Zhuhai Xingnuo Energy Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhuhai Xingnuo Energy Technology Co Ltd filed Critical Zhuhai Xingnuo Energy Technology Co Ltd
Priority to CN201910179979.6A priority Critical patent/CN110069984A/en
Publication of CN110069984A publication Critical patent/CN110069984A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/02Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/08Feature extraction
    • G06F2218/10Feature extraction by analysing the shape of a waveform, e.g. extracting parameters relating to peaks

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring Frequencies, Analyzing Spectra (AREA)

Abstract

The invention discloses one-variable linear regression frequency measuring method, include the following steps: the sinusoidal signal for obtaining DC charging motor input voltage, the sinusoidal signal of dead-center position is subjected to linearization process;The sinusoidal signal of dead-center position is defined as about the linear function of time t and being superimposed for random parameter;The sinusoidal signal of dead-center position is converted into discrete digital signal, to obtain the sampling of sinusoidal signal value of several dead-center positions;Regression analysis is carried out to the sampling of sinusoidal signal value, to obtain the best estimate of parametric variable in the linear function about time t;The best time estimated value when sinusoidal signal zero crossing position is calculated according to the best estimate.The present invention, which avoids traditional zero-crossing examination method using the method for linear regression, will receive the unstable characteristic of the longer measurement result of filtering time caused by the reasons such as Harmonious Waves in Power Systems, noise and interference, frequency computational accuracy is improved, and more accurate to the zero crossing position acquisition of actual signal.

Description

One-variable linear regression frequency measuring method, electronic equipment and storage medium
Technical field
The present invention designs the accurate survey that a kind of one-variable linear regression Frequency Measurement Algorithm is applied to DC charging motor input source frequency Amount.
Background technique
Currently, the non-vehicle-mounted type DC charging motor in market all has the high feature of charge power, common integral type direct current Charger power is differed from 30kw to 210kw, and the more power of charger quantity at direct current charge station are big and concentrate, if more are filled Motor, which works at the same time, will cause network load increase, causes system active power insufficient, as a result generator speed is caused to decline, electricity Net frequency reduces;On the contrary, it is superfluous that it will cause system active power if network load reduces, as a result lead to generator speed Rise, mains frequency increases.If wanting to make mains frequency to revert to rated frequency operation at this time, generator prime machine need to increase or The input power of reduction, or network load is made to generate variation, i.e., charger, which has, changes output power according to mains frequency fluctuation Ability.
But currently existing scheme has the following deficiencies:
Mains frequency has many detection methods as important performance characteristic at present, and wherein cross zero detecting method is a kind of letter Single practical, very widely used frequency measurement method.Its principle be by judge the time difference of two zero crossing of sinusoidal signal come Obtain frequency.But reality is in the application, since signal contains harmonic wave, by reasons such as ambient noise interferences, causes to measure reality Border signal zero-crossing positional fluctuation limits the application of this method so that the frequency error actually calculated is bigger than normal.Currently on the market Most DC charging motors will not change output power in frequency fluctuation.
Summary of the invention
For overcome the deficiencies in the prior art, one of the objects of the present invention is to provide one-variable linear regression frequency measuring method, Its prior art that can solve has deviation to dead-center position accuracy, is unable to adjust output power problem.
The second object of the present invention is to provide a kind of electronic equipment, and the prior art that can solve is accurate to dead-center position Degree has deviation, is unable to adjust output power problem.
The third object of the present invention is to provide a kind of computer readable storage medium, and the prior art that can solve is to zero Point position accuracy has deviation, is unable to adjust output power problem.
An object of the present invention adopts the following technical scheme that realization:
One-variable linear regression frequency measuring method, which comprises the steps of:
Linear process step: obtaining the sinusoidal signal of DC charging motor input voltage, by the sinusoidal letter of dead-center position Number carry out linearization process;
Function superposition step: the sinusoidal signal of dead-center position is defined as about the linear function of time t and random The superposition of parameter;
Signal switch process: the sinusoidal signal of dead-center position is converted into discrete digital signal, to obtain several The sampling of sinusoidal signal value of dead-center position;
Regression analysis step: regression analysis is carried out to the sampling of sinusoidal signal value, to obtain the linear letter about time t The best estimate of parametric variable in number;
Time valuation step: the best time when sinusoidal signal zero crossing position is calculated according to the best estimate Estimated value.
Further, sinusoidal signal u (t)=at+b+ ε of dead-center position is defined in function superposition step, wherein at+b For the linear function about time t, a and b are parametric variable, and ε is random number.
Further, the sampling of sinusoidal signal value of several dead-center positions is (t in signal switch processi, ui), i=1, 2 ... ..., n.
Further, in regression analysis step, best estimate is obtained especially by following steps:
In regression analysis step, best estimate is obtained especially by following manner:
Pass through formula u (t)=at+b+ ε and sampling of sinusoidal signal value (ti, ui) be calculatedAnd
DefinitionWhenMoment best estimate can then be calculated WhereinFor the best estimate of a,For the best estimate of b,For the best estimate of u.
The second object of the present invention adopts the following technical scheme that realization:
A kind of electronic equipment including memory, processor and stores the meter that can be run on a memory and in processor Calculation machine program, the processor perform the steps of when executing the computer program
Linear process step: obtaining the sinusoidal signal of DC charging motor input voltage, by the sinusoidal letter of dead-center position Number carry out linearization process;
Function superposition step: the sinusoidal signal of dead-center position is defined as about the linear function of time t and random The superposition of parameter;
Signal switch process: the sinusoidal signal of dead-center position is converted into discrete digital signal, to obtain several The sampling of sinusoidal signal value of dead-center position;
Regression analysis step: regression analysis is carried out to the sampling of sinusoidal signal value, to obtain the linear letter about time t The best estimate of parametric variable in number;
Time valuation step: the best time when sinusoidal signal zero crossing position is calculated according to the best estimate Estimated value.
Further, sinusoidal signal u (t)=at+b+ ε of dead-center position is defined in function superposition step, wherein at+ B is the linear function about time t, and a and b are parametric variable, and ε is random number.
Further, the sampling of sinusoidal signal value of several dead-center positions is (t in signal switch processi, ui), i=1, 2 ... ..., n.
Further, in regression analysis step, best estimate is obtained especially by following manner:
Pass through formula u (t)=at+b+ ε and sampling of sinusoidal signal value (ti, ui) be calculatedAnd
DefinitionWhenMoment best estimate can then be calculated WhereinFor the best estimate of a,For the best estimate of b,For the best estimate of u.
The third object of the present invention adopts the following technical scheme that realization:
A kind of computer readable storage medium, is stored thereon with computer program, and the computer program is held by processor One-variable linear regression frequency measuring method described in the above-mentioned any one of row.
Compared with prior art, the beneficial effects of the present invention are
The present invention avoids traditional zero-crossing examination method using the method for linear regression and will receive Harmonious Waves in Power Systems, makes an uproar The unstable characteristic of the longer measurement result of filtering time caused by the reasons such as sound and interference, improves frequency computational accuracy, and It is more accurate to the zero crossing position acquisition of actual signal.
Detailed description of the invention
Fig. 1 is the flow chart of one-variable linear regression frequency measuring method of the invention;
Fig. 2 is the frequency waveform diagram that frequency is 50Hz;
Fig. 3 is charger downslope time figure.
Specific embodiment
In the following, being described further in conjunction with attached drawing and specific embodiment to the present invention, it should be noted that not Under the premise of conflicting, new implementation can be formed between various embodiments described below or between each technical characteristic in any combination Example.
Embodiment one
As shown in Figure 1, the so-called dead-center position of the present invention exists the present invention provides a kind of one-variable linear regression frequency measuring method It include zero crossings position in practical application.The present invention be based on sinusoidal signal zero crossings can approximate linearization it is substantially former The signal value of near zero-crossing point is carried out simple linear regression analysis, is somebody's turn to do by reason after detecting that signal crosses zero position roughly Intercept of the signal on time axial coordinate, as the signal zero-crossing moment, simple linear regression analysis are believed zero crossings It is number for statistical analysis, while analyzing the result signal statistics analysis method based on principle of least square method again, i.e., it is calculated Result out meets error sum of squares minimum relationship.Specifically comprise the following steps:
S1: obtaining the sinusoidal signal of DC charging motor input voltage, the sinusoidal signal of dead-center position is carried out linear Change processing;
SIN function is decomposed into the form of power series:When signal zero is attached When nearly x value very little, it is approximately sinx ≈ x that high-order term formula, which can be ignored,.If approximate expression is required to be less than or be waited with true value relative error In 0.1%, then the radian value range of x is x ∈ [- 0.0774,0.0774].
S2: the sinusoidal signal of dead-center position is defined as about the linear function of time t and being superimposed for random parameter;
In electric system, voltageMacroscopic view sees, system voltage waveform, amplitude and time variable T is sinusoidal functional relation, generally belongs to deterministic dependence, and zero crossings can local linearization;But due to electrically making an uproar The influence of sound and system interference, the parameter voltages U of the SIN functionm, frequency f and phaseIt is all to obey a fixed distribution function Stochastic variable causes system that certain randomness is presented in local features such as zero crossing times.Therefore, we can be detected letter Number u (t) is regarded as in the observation result of zero crossings and is formed by stacking by two parts: a part causes by the linear function of time t, It is denoted as at+b;Another part is to be denoted as the sinusoidal letter that u (t)=at+b+ ε, u (t) is dead-center position as caused by enchancement factor Number, at+b is the linear function about time t, and a and b are parametric variable, and ε is random number.Parametric variable a and b is mainly in formula It is determined by the mathematical expectation of the voltage of signal, frequency, phase three's stochastic variable;For random number ε, it is to cause signal mistake Crossing point jitter, or even repeat the principal element of zero passage.
S3: the sinusoidal signal of dead-center position is converted into discrete digital signal, to obtain several dead-center positions Sampling of sinusoidal signal value;
U (t) is analog signal, in digital detection system, is usually sampled as analog signal u (t) by analog-to-digital conversion Discrete digital signal, borrow statistics in data regression thought, using n sampled value near tested signal zero as Independent observation sample.The sampling of sinusoidal signal value of several dead-center positions is (ti, ui), i=1,2 ... ..., n carry out unitary line Property regression analysis, obtains the best estimate of formula parametric variable a, b.
S4: carrying out regression analysis to the sampling of sinusoidal signal value, is become with obtaining parameter in the linear function about time t The best estimate of amount;
Pass through formula u (t)=at+b+ ε and sampling of sinusoidal signal value (ti, ui) be calculatedAnd
WhereinFor the best estimate of a,For the best estimate of b.
S5: the best time estimated value when sinusoidal signal zero crossing position is calculated according to the best estimate.
DefinitionWhenIt can then be calculated For the best estimate of u
The as best estimate zero-acrross ing moment value of current demand signal zero crossing.
In actual application, 80 points of the every cycle sampling number of power network signal, one-variable linear regression frequency measurement calculates point used It can be even number that number n, which can be odd number,.Specific formula for calculation is as follows:
One-variable linear regression points n is the linear regression frequency measurement calculation formula of even number point are as follows:
One-variable linear regression counts n as surprise The linear regression frequency measurement calculation formula of several points are as follows: Taking n in practical applications is even number point convenient for programming realization.
Following table is the simulation result for using one-variable linear regression points to obtain for even number point.Observation linear regression frequency measurement obtains Signal different frequency (40Hz~60Hz) measurement result and precision.Essential record offset signal actual frequency or so is maximum Error.
Frequency is the frequency wave shape that 50Hz is calculated as shown in Fig. 2, output power changes according to frequency fluctuation, using frequency wave Dynamic 0.5Hz power decline 10%, charger electric current decline as shown in Figure 3 are completed in 1s.
Embodiment two
Embodiment two discloses a kind of electronic equipment, including memory, processor and storage are on a memory and can be The computer program of processor operation, the processor perform the steps of when executing the computer program
Linear process step: obtaining the sinusoidal signal of DC charging motor input voltage, by the sinusoidal letter of dead-center position Number carry out linearization process;
Function superposition step: the sinusoidal signal of dead-center position is defined as about the linear function of time t and random The superposition of parameter;
Signal switch process: the sinusoidal signal of dead-center position is converted into discrete digital signal, to obtain several The sampling of sinusoidal signal value of dead-center position;
Regression analysis step: regression analysis is carried out to the sampling of sinusoidal signal value, to obtain the linear letter about time t The best estimate of parametric variable in number;
Time valuation step: the best time when sinusoidal signal zero crossing position is calculated according to the best estimate Estimated value.
Sinusoidal signal u (t)=at+b+ ε of dead-center position is defined in function superposition step, wherein at+b is about the time The linear function of t, a and b are parametric variable, and ε is random number.The sinusoidal signal of several dead-center positions in signal switch process Sampled value is (ti, ui), i=1,2 ... ..., n.
In regression analysis step, best estimate is obtained especially by following manner:
Pass through formula u (t)=at+b+ ε and sampling of sinusoidal signal value (ti, ui) be calculatedAnd
DefinitionWhenMoment best estimate can then be calculated WhereinFor the best estimate of a,For the best estimate of b,For the best estimate of u.
Embodiment three
Embodiment three discloses a kind of readable computer storage medium, which is somebody's turn to do for storing program When program is executed by processor, the one-variable linear regression frequency measuring method of embodiment one is realized.
The above embodiment is only the preferred embodiment of the present invention, and the scope of protection of the present invention is not limited thereto, The variation and replacement for any unsubstantiality that those skilled in the art is done on the basis of the present invention belong to institute of the present invention Claimed range.

Claims (9)

1. one-variable linear regression frequency measuring method, which comprises the steps of:
Linear process step: obtain DC charging motor input voltage sinusoidal signal, by the sinusoidal signal of dead-center position into Row linearization process;
Function superposition step: the sinusoidal signal of dead-center position is defined as the linear function and random parameter about time t Superposition;
Signal switch process: the sinusoidal signal of dead-center position is converted into discrete digital signal, to obtain several zero points The sampling of sinusoidal signal value of position;
Regression analysis step: regression analysis is carried out to the sampling of sinusoidal signal value, to obtain in the linear function about time t The best estimate of parametric variable;
Time valuation step: the best time estimation when sinusoidal signal zero crossing position is calculated according to the best estimate Value.
2. one-variable linear regression frequency measuring method as described in claim 1, which is characterized in that define zero point in function superposition step The sinusoidal signal u (t) of position=at+b+ ε, wherein at+b is the linear function about time t, and a and b are parametric variable, and ε is Random number.
3. one-variable linear regression frequency measuring method as claimed in claim 2, which is characterized in that in signal switch process several The sampling of sinusoidal signal value of dead-center position is (ti, ui), i=1,2 ... ..., n.
4. one-variable linear regression frequency measuring method as claimed in claim 3, which is characterized in that in regression analysis step, obtain Best estimate is especially by following manner:
Pass through formula u (t)=at+b+ ε and sampling of sinusoidal signal value (ti, ui) be calculated And
DefinitionWhenMoment best estimate can then be calculated WhereinFor the best estimate of a,For the best estimate of b,For the best estimate of u.
5. a kind of electronic equipment including memory, processor and stores the calculating that can be run on a memory and in processor Machine program, which is characterized in that the processor performs the steps of when executing the computer program
Linear process step: obtain DC charging motor input voltage sinusoidal signal, by the sinusoidal signal of dead-center position into Row linearization process;
Function superposition step: the sinusoidal signal of dead-center position is defined as the linear function and random parameter about time t Superposition;
Signal switch process: the sinusoidal signal of dead-center position is converted into discrete digital signal, to obtain several zero points The sampling of sinusoidal signal value of position;
Regression analysis step: regression analysis is carried out to the sampling of sinusoidal signal value, to obtain in the linear function about time t The best estimate of parametric variable;
Time valuation step: the best time estimation when sinusoidal signal zero crossing position is calculated according to the best estimate Value.
6. electronic equipment as claimed in claim 5, which is characterized in that define the sinusoidal letter of dead-center position in function superposition step Number u (t)=at+b+ ε, wherein at+b, a are the linear function about time t, and a and b are parametric variable, and ε is random number.
7. electronic equipment as claimed in claim 6, which is characterized in that several dead-center positions are being just in signal switch process String signal sampled value is (ti, ui), i=1,2 ... ..., n.
8. electronic equipment as claimed in claim 7, which is characterized in that in regression analysis step, in regression analysis step, Best estimate is obtained especially by following manner:
Pass through formula u (t)=at+b+ ε and sampling of sinusoidal signal value (ti, ui) be calculated And
DefinitionWhenMoment best estimate can then be calculated WhereinFor the best estimate of a,For the best estimate of b,For the best estimate of u.
9. a kind of computer readable storage medium, is stored thereon with computer program, it is characterised in that: the computer program quilt The one-variable linear regression frequency measuring method as described in claim 1-4 any one is realized when processor executes.
CN201910179979.6A 2019-03-11 2019-03-11 One-variable linear regression frequency measuring method, electronic equipment and storage medium Pending CN110069984A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910179979.6A CN110069984A (en) 2019-03-11 2019-03-11 One-variable linear regression frequency measuring method, electronic equipment and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910179979.6A CN110069984A (en) 2019-03-11 2019-03-11 One-variable linear regression frequency measuring method, electronic equipment and storage medium

Publications (1)

Publication Number Publication Date
CN110069984A true CN110069984A (en) 2019-07-30

Family

ID=67366197

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910179979.6A Pending CN110069984A (en) 2019-03-11 2019-03-11 One-variable linear regression frequency measuring method, electronic equipment and storage medium

Country Status (1)

Country Link
CN (1) CN110069984A (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101871965A (en) * 2010-06-11 2010-10-27 威胜集团有限公司 Method for detecting zero crossing time, frequency and phase difference of power sinusoidal signals
CN102542170A (en) * 2012-01-16 2012-07-04 南京磐能电力科技股份有限公司 Method for acquiring slip frequency difference of electric system
CN105044456A (en) * 2015-07-21 2015-11-11 电子科技大学 Power grid instantaneous frequency measuring and tracking method based on orthogonal sub-band
CN107834678A (en) * 2017-08-17 2018-03-23 山东电工豪迈节能科技有限公司 Application of the one-variable linear regression Frequency Measurement Algorithm in DC charging motor
CN109374966A (en) * 2018-10-23 2019-02-22 国网重庆市电力公司电力科学研究院 A kind of mains frequency estimation method

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101871965A (en) * 2010-06-11 2010-10-27 威胜集团有限公司 Method for detecting zero crossing time, frequency and phase difference of power sinusoidal signals
CN102542170A (en) * 2012-01-16 2012-07-04 南京磐能电力科技股份有限公司 Method for acquiring slip frequency difference of electric system
CN105044456A (en) * 2015-07-21 2015-11-11 电子科技大学 Power grid instantaneous frequency measuring and tracking method based on orthogonal sub-band
CN107834678A (en) * 2017-08-17 2018-03-23 山东电工豪迈节能科技有限公司 Application of the one-variable linear regression Frequency Measurement Algorithm in DC charging motor
CN109374966A (en) * 2018-10-23 2019-02-22 国网重庆市电力公司电力科学研究院 A kind of mains frequency estimation method

Similar Documents

Publication Publication Date Title
Jain et al. Fast harmonic estimation of stationary and time-varying signals using EA-AWNN
CN104316894B (en) A kind of realistically displayed method of electric energy meter actual motion environment
CN102033161B (en) Frequency measuring method of alternating current signal
CN109633262A (en) Three phase harmonic electric energy gauging method, device based on composite window multiline FFT
CN109752584B (en) Method for measuring effective value of periodic signal
CN102981049B (en) Frequency detection method used for microgrid system
CN110045186B (en) System harmonic impedance estimation method and system at power system common connection point
WO2017113754A1 (en) Method for detecting dynamic reactive power of grid and reactive power detector
CN102253281A (en) Method for measuring power grid frequency by resisting harmonic interference and automatically changing phase based on software
Zhang et al. A non-intrusive identification method of harmonic source loads for industrial users
CN102095935A (en) Method for measuring instantaneous phase of alternating current electrical signal
CN109446643B (en) Method for establishing household appliance load harmonic model based on measured data
CN104808055B (en) A kind of digital measuring method of signal frequency
CN112505391B (en) Frequency-adaptive alternating current signal effective value acquisition method
CN103513133B (en) The measurement apparatus of MCR type SVC device dynamic response time and method
CN110069984A (en) One-variable linear regression frequency measuring method, electronic equipment and storage medium
CN110007129B (en) A kind of three-phase voltage real-time estimation method applied to dynamic electric energy metering
Liu et al. Asynchronous harmonic analysis based on out-of-sequence measurement for large-scale residential power network
CN111769595B (en) Wind driven generator network side converter equivalent capacitance solving method based on parameter identification
CN109541308B (en) Harmonic wave analysis method based on distributed quasi-synchronous power grid measurement system
CN107834678A (en) Application of the one-variable linear regression Frequency Measurement Algorithm in DC charging motor
CN114184838A (en) Power system harmonic detection method, system and medium based on SN mutual convolution window
Di Leo et al. Smart meters and demand response: A double challenge
CN106052562B (en) A kind of when grid displacement measurement method and signal processing system of adaptive dynamic than phase
Luo et al. A novel method to dertermine the harmonic responsibility for power distribution network including photovoltaic system

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20190730