CN105548694B - Frequency measurement method and system are carried out according to the zero initial phase reference sinusoidal sequence of function - Google Patents
Frequency measurement method and system are carried out according to the zero initial phase reference sinusoidal sequence of function Download PDFInfo
- Publication number
- CN105548694B CN105548694B CN201510891522.XA CN201510891522A CN105548694B CN 105548694 B CN105548694 B CN 105548694B CN 201510891522 A CN201510891522 A CN 201510891522A CN 105548694 B CN105548694 B CN 105548694B
- Authority
- CN
- China
- Prior art keywords
- sequence
- frequency
- initial phase
- phase
- positive
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/02—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
Frequency measurement method and system are carried out according to the zero initial phase reference sinusoidal sequence of function the present invention relates to a kind of, it obtains real frequency mixing sequence and empty frequency is mixed sequence, then digital filtering and integration are carried out to the real frequency mixing sequence and the empty frequency mixing sequence, obtain real frequency integrated value and empty frequency integrated value, according to the real frequency integrated value and empty frequency integrated value, the phase of zero initial phase reference sinusoidal FUNCTION MODULATION sequence is obtained, the frequency of electric power signal is obtained according to the phase.The frequency of electric power signal is obtained by mode of the present invention, is remarkably improved the accuracy of frequency power signal measurement.
Description
Technical field
The present invention relates to technical field of power systems, more particularly to a kind of according to the zero initial phase reference sinusoidal sequence of function
Carry out frequency measurement method and system.
Background technology
Frequency measurement, phase measurement, amplitude measurement of electric system etc. are the measurement of sine parameter in itself.Fu Li
Leaf transformation is the basic skills for realizing sine parameter measurement, is widely used in electric system.But as sine parameter is surveyed
The development of amount technology, Fourier transformation there are the problem of also more aobvious protrusion, be difficult to further to meet electric system and sine joined
The requirement that number high accuracy calculates.
In terms of electric system sine parameter measurement, there are some improved measurement method of parameters, such as zero hands over method, based on filter
The mensuration of ripple, the mensuration based on Wavelet Transform, based on neutral net, based on DFT (Discrete Fourier
Transform, discrete Fourier transform) conversion mensuration etc..Since the specified power frequency of operation of power networks is near 50Hz (hertz),
Belong to the relatively low sinusoidal frequency of frequency, and there is interference in actual electric power signal, such as harmonic wave interference, electric load are small
In the range of the similar white noise acoustic jamming etc. that generates of random fluctuation, under interference environment, the measurement of these algorithm generally existings is accurate
The problem of degree is not high.
The content of the invention
Based on this, it is necessary in view of the above-mentioned problems, provide it is a kind of according to the zero initial phase reference sinusoidal sequence of function into line frequency
Rate measuring method and system can improve the accuracy of frequency power signal measurement.
In order to solve the above technical problems, the present invention adopts the following technical scheme that：
It is a kind of that frequency measurement method is carried out according to the zero initial phase reference sinusoidal sequence of function, including step：
According to the lower limit of frequency power signal scope, preset sample frequency and default integer signal period number, obtain preliminary
Sequence length；
Electric power signal is sampled according to the preliminary sequence length, obtains the preliminary sequence of the electric power signal；
To the preliminary sequence into line frequency preliminary survey, the preliminary frequency of the electric power signal is obtained, according to the just cadence
Rate obtains reference frequency；
According to the preset sample frequency and the reference frequency, the unit period sequence for obtaining the electric power signal is long
Degree；
The default integer signal period number with the unit period sequence length is multiplied, it is long to obtain preprocessing sequence
Degree；
According to the preprocessing sequence length, preprocessing sequence is obtained from the preliminary sequence of the electric power signal；
Comb filtering processing is carried out to the preprocessing sequence, obtains comb filtering sequence, wherein comb filtering sequence is long
It spends and is carrying out comb filtering treated residue length for the preprocessing sequence；
The ratio integer of the comb filtering sequence length and the unit period sequence length is determined, according to the ratio
Integer and the unit period sequence length obtain predetermined sequence length, wherein the ratio integer is odd number, the default sequence
Row length is odd number；
According to the predetermined sequence length and default starting point, the first positive sequence is obtained from the comb filtering sequence
Row obtain the first anti-pleat sequence according to the described first positive sequence；
First positive phase is obtained according to the described first positive sequence, the first reverse phase is obtained according to the described first anti-pleat sequence
Position；
First average initial phase is obtained according to first positive phase and first antiphase；
According to the described first average initial phase and preset phase value, obtain phase fiducial value, according to the phase fiducial value,
The default starting point and the unit period sequence length, obtain new starting point；
According to the predetermined sequence length and the new starting point, the second positive sequence is obtained from the comb filtering sequence
Row obtain the second anti-pleat sequence according to the second positive sequence；
Second positive phase is obtained according to the described second positive sequence, the second reverse phase is obtained according to the described second anti-pleat sequence
Position；
Second average initial phase is obtained according to second positive phase and second antiphase；
Described second positive sequence with the described second anti-pleat sequence is subtracted each other, obtains difference sequence, according to the difference sequence with
The sine function of described second average initial phase, obtains SIN function modulation sequence；
It is exported from the SIN function modulation sequence central point, obtains zero initial phase reference sinusoidal FUNCTION MODULATION sequence；
The zero initial phase reference sinusoidal FUNCTION MODULATION sequence with the discrete cosine function of reference frequency is multiplied, is obtained real
Frequency mixing sequence, the zero initial phase reference sinusoidal FUNCTION MODULATION sequence with the discrete sine function of reference frequency is multiplied, is obtained
Obtain empty frequency mixing sequence；
Digital filtering is carried out to the real frequency mixing sequence and the empty frequency mixing sequence respectively, obtains real frequency filtered sequence
With imaginary frequency filtering sequence；
The real frequency filtered sequence and the imaginary frequency filtering sequence are integrated respectively, obtain real frequency integrated value and empty frequency
Integrated value；
According to the real frequency integrated value and empty frequency integrated value, the phase of zero initial phase reference sinusoidal FUNCTION MODULATION sequence of acquisition
Position obtains the frequency of electric power signal according to the phase.
It is a kind of that frequency measuring system is carried out according to the zero initial phase reference sinusoidal sequence of function, including：
Preliminary sequence length determination modul, for lower limit, the preset sample frequency and pre- according to frequency power signal scope
If integer signal period number, preliminary sequence length is obtained；
Preliminary sequence acquisition module, for being sampled according to the preliminary sequence length to electric power signal, described in acquisition
The preliminary sequence of electric power signal；
Reference frequency determining module, for, into line frequency preliminary survey, obtaining the first of the electric power signal to the preliminary sequence
Synchronizing frequency obtains reference frequency according to the preliminary frequency；
Unit period sequence length determining module, for according to the preset sample frequency and the reference frequency, obtaining
The unit period sequence length of the electric power signal；
Preprocessing sequence length determination modul, for by the default integer signal period number and the unit period sequence
Length is multiplied, and obtains preprocessing sequence length；
Preprocessing sequence acquisition module, for according to the preprocessing sequence length, from the preliminary sequence of the electric power signal
Preprocessing sequence is obtained in row；
Comb filtering sequence determining module for carrying out comb filtering processing to the preprocessing sequence, obtains pectination filter
Wave train, wherein comb filtering sequence length are carrying out comb filtering treated residue length for the preprocessing sequence；
Predetermined sequence length determination modul, for determining that the comb filtering sequence length and the unit period sequence are long
The ratio integer of degree obtains predetermined sequence length, wherein described according to the ratio integer and the unit period sequence length
Ratio integer is odd number, and the predetermined sequence length is odd number；
First ray acquisition module, for according to the predetermined sequence length and default starting point, from the comb filtering
The first positive sequence is obtained in sequence, the first anti-pleat sequence is obtained according to the described first positive sequence；
First positive and negative phase determination module, for obtaining the first positive phase according to the described first positive sequence, according to described
First anti-pleat sequence obtains the first antiphase；
First average initial phase determining module, for obtaining first according to first positive phase and first antiphase
Average initial phase；
New starting point determining module, for according to the described first average initial phase and preset phase value, obtaining phase bit comparison
Value, according to the phase fiducial value, the default starting point and the unit period sequence length, obtains new starting point；
Second retrieval module, for according to the predetermined sequence length and the new starting point, being filtered from the pectination
The second positive sequence is obtained in wave train, the second anti-pleat sequence is obtained according to the second positive sequence；
Second positive and negative phase determination module, for obtaining the second positive phase according to the described second positive sequence, according to described
Second anti-pleat sequence obtains the second antiphase；
Second average initial phase determining module, for obtaining second according to second positive phase and second antiphase
Average initial phase；
SIN function modulation sequence determining module, for by the described second positive sequence and the described second anti-pleat sequence phase
Subtract, obtain difference sequence, according to the difference sequence and the sine function of the second average initial phase, obtain SIN function modulation
Sequence；
Zero initial phase modulation sequence acquisition module for being exported from the SIN function modulation sequence central point, obtains zero
Initial phase reference sinusoidal FUNCTION MODULATION sequence；
Be mixed sequence determining module, for by the zero initial phase reference sinusoidal FUNCTION MODULATION sequence and reference frequency from
It dissipates cosine function to be multiplied, real frequency mixing sequence is obtained, by the zero initial phase reference sinusoidal FUNCTION MODULATION sequence and reference frequency
Discrete sine function be multiplied, obtain empty frequency mixing sequence；
Filtered sequence determining module, for carrying out number to the real frequency mixing sequence and the empty frequency mixing sequence respectively
Filtering obtains real frequency filtered sequence and imaginary frequency filtering sequence；
Integrated value determining module, for being integrated respectively to the real frequency filtered sequence and the imaginary frequency filtering sequence,
Obtain real frequency integrated value and empty frequency integrated value；
Frequency power signal determining module, for according to the real frequency integrated value and empty frequency integrated value, obtaining zero initial phase
The phase of reference sinusoidal FUNCTION MODULATION sequence obtains the frequency of electric power signal according to the phase.
The present invention carries out frequency measurement method and system according to the zero initial phase reference sinusoidal sequence of function, obtains real frequency and is mixed
Sequence and empty frequency are mixed sequence, then carry out digital filtering and product to the real frequency mixing sequence and the empty frequency mixing sequence
Point, real frequency integrated value and empty frequency integrated value are obtained, according to the real frequency integrated value and empty frequency integrated value, obtains zero initial phase benchmark
The phase of SIN function modulation sequence obtains the frequency of electric power signal according to the phase.Electric power is obtained by mode of the present invention
The frequency of signal is remarkably improved the accuracy of frequency power signal measurement.
Description of the drawings
Fig. 1 is that the present invention shows according to the flow of zero initial phase reference sinusoidal sequence of function progress frequency measurement method embodiment
It is intended to；
Fig. 2 is schematic diagram of the comb filtering of the present invention processing in frequency domain amplitude versus frequency characte；
Fig. 3 is comb filtering sequence of the present invention, the first positive sequence, the schematic diagram of the first anti-pleat sequence；
Fig. 4 is the pictorial diagram of zero initial phase datum mark of the invention；
Fig. 5 is that the present invention shows according to the structure of zero initial phase reference sinusoidal sequence of function progress frequency measuring system embodiment
It is intended to.
Specific embodiment
It is below in conjunction with the accompanying drawings and preferably real further to illustrate the effect of the technological means of the invention taken and acquirement
Example is applied, to technical scheme, carries out clear and complete description.
As shown in Figure 1, a kind of carry out frequency measurement method according to the zero initial phase reference sinusoidal sequence of function, including step：
S101, according to the lower limit of frequency power signal scope, preset sample frequency and default integer signal period number, obtain
Preliminary sequence length；
S102, electric power signal is sampled according to the preliminary sequence length, obtains the preliminary sequence of the electric power signal
Row；
S103, to the preliminary sequence into line frequency preliminary survey, obtain the preliminary frequency of the electric power signal, according to it is described just
Synchronizing frequency obtains reference frequency；
S104, according to the preset sample frequency and the reference frequency, obtain the unit period sequence of the electric power signal
Row length；
S105, the default integer signal period number with the unit period sequence length is multiplied, obtains pretreatment sequence
Row length；
S106, according to the preprocessing sequence length, obtain preprocessing sequence from the preliminary sequence of the electric power signal；
S107, comb filtering processing is carried out to the preprocessing sequence, obtains comb filtering sequence, wherein comb filtering sequence
Row length is carrying out comb filtering treated residue length for the preprocessing sequence；
S108, the ratio integer for determining the comb filtering sequence length and the unit period sequence length, according to institute
It states ratio integer and the unit period sequence length obtains predetermined sequence length, wherein the ratio integer is odd number, it is described
Predetermined sequence length is odd number；
S109, according to the predetermined sequence length and default starting point, obtaining first just from the comb filtering sequence
To sequence, the first anti-pleat sequence is obtained according to the described first positive sequence；
S110, the first positive phase is obtained according to the described first positive sequence, first is obtained according to the described first anti-pleat sequence
Antiphase；
S111, the first average initial phase is obtained according to first positive phase and first antiphase；
S112, be averaged initial phase and preset phase value according to described first, phase fiducial value is obtained, according to the phase ratio
Compared with value, the default starting point and the unit period sequence length, new starting point is obtained；
S113, according to the predetermined sequence length and the new starting point, obtain second from the comb filtering sequence
Positive sequence obtains the second anti-pleat sequence according to the second positive sequence；
S114, the second positive phase is obtained according to the described second positive sequence, second is obtained according to the described second anti-pleat sequence
Antiphase；
S115, the second average initial phase is obtained according to second positive phase and second antiphase；
S116, the described second positive sequence with the described second anti-pleat sequence is subtracted each other, obtains difference sequence, according to the poor sequence
The sine function of row and the described second average initial phase, obtains SIN function modulation sequence；
S117, exported from the SIN function modulation sequence central point, obtain zero initial phase reference sinusoidal FUNCTION MODULATION sequence
Row；
S118, the zero initial phase reference sinusoidal FUNCTION MODULATION sequence is multiplied with the discrete cosine function of reference frequency,
Real frequency mixing sequence is obtained, by the discrete sine function phase of the zero initial phase reference sinusoidal FUNCTION MODULATION sequence and reference frequency
Multiply, obtain empty frequency mixing sequence；
S119, digital filtering is carried out to the real frequency mixing sequence and the empty frequency mixing sequence respectively, obtains real frequency and filter
Wave train and imaginary frequency filtering sequence；
S120, the real frequency filtered sequence and the imaginary frequency filtering sequence are integrated respectively, obtains real frequency integrated value
With empty frequency integrated value；
S121, according to the real frequency integrated value and empty frequency integrated value, obtain zero initial phase reference sinusoidal FUNCTION MODULATION sequence
Phase, according to the phase obtain electric power signal frequency.
Actual electric power signal is a kind of sinusoidal signal based on fundamental wave component, and unless otherwise specified, electric power signal refers both to
Fundamental signal, frequency power signal refer both to fundamental frequency.
For step S101, power system frequency scope is generally in 45Hz (hertz) -55Hz, so frequency power signal model
The lower limit f enclosedminIt can be taken as 45Hz.Default integer signal period number C2πIt can set according to actual needs, for example, by C2πIt takes
For 17.Preliminary sequence length can be calculated according to formula (1)：
Wherein, NstartFor preliminary sequence length；(int) rounding is represented；C2πTo preset integer signal period number；fminFor electricity
The lower limit of force signal frequency range, unit Hz；F is preset sample frequency, unit Hz.
For step S102, the cosine function signal representation electric power signal of single fundamental frequency can be used, then preliminary sequence
For formula (2)：
Wherein, Xstart(n) it is preliminary sequence；A is signal amplitude, unit v；ωiFor signal frequency, when T is the sampling interval
Between, f is preset sample frequency, and unit Hz, n are series of discrete number,For the initial phase of preliminary sequence, NstartIt is long for preliminary sequence
Degree.
For step S103, can by zero hand over method, the algorithm based on filtering, based on Wavelet Transformation Algorithm, based on nerve net
The algorithm of network, the frequency algorithm based on DFT transform or based on the frequency algorithm of phase difference to preliminary sequence into line frequency preliminary survey, obtain
Take preliminary frequency ωo.In one embodiment, the reference frequency ωs=ωo。
For step S104, in one embodiment, the unit period sequence length of the electric power signal calculates, and is formula
(3)：
Wherein, N2πFor unit periodic sequence length；(int) it is round numbers；F is preset sample frequency, unit Hz；ωsFor
Reference frequency.There are the errors in 1 sampling interval for the unit period sequence length integer.
For step S105, the preprocessing sequence length is formula (4)：
Nset=C2πN2π (4)
Wherein, NsetFor preprocessing sequence length；C2πTo preset integer signal period number；N2πIt is long for unit periodic sequence
Degree.
For step S106, in one embodiment, preprocessing sequence is obtained, be formula (5)：
Wherein, Xset(n) it is preprocessing sequence, Xstart(n) it is preliminary sequence, NsetFor preprocessing sequence length.
There are subharmonic, even-order harmonic, odd harmonics etc., in optical mixing process, generated mixing in electric power signal
Interfering frequency seriously affects the accuracy in computation of sine parameter, and carrying out comb filtering processing can be from source to mixing interfering frequency production
Raw factor is effectively inhibited.
For step S107,2 different list entries of starting point are subtracted each other, the frequency domain amplitude-frequency filtering that can obtain pectination is special
Property, referred to as comb filtering processing.The interval of 2 list entries is defined as comb filtering parameter, the processing of single-stage comb filtering
It is expressed as formula (6)：
Wherein, XL(n) it is single-stage comb filtering output sequence；Xset(n) it is preprocessing sequence；Xset(NL+ n) it is from NLIt rises
The preprocessing sequence of beginning；NLFor 2 train intervals or single-stage comb filtering parameter；NsetFor preprocessing sequence length；N2πFor unit
Periodic sequence length.Comb filtering parameter NLValue is unit periodic sequence length N2π0.5 times, even-order harmonic can be carried out
Inhibit and decay to subharmonic.
In one embodiment, can comb filtering processing be carried out to the preprocessing sequence by comb filter.By
, in order to improve comb-filter effects, 8 can be passed through there are integer error there are error, comb filtering parameter in reference frequency
Grade comb filter carries out comb filtering processing, is expressed as formula (7)：
Wherein, X8L(n) it is 8 grades of comb filter or comb filtering sequence；Filter[8,NL,Xset(n)] 8 represent comb in
Shape filtering series is 8, NLFor single-stage comb filtering parameter, Xset(n) it is preprocessing sequence；KL(ωi) exist for comb filtering processing
Signal frequency ωiDimensionless amplitude gain, T is sampling interval duration；For comb filtering sequence initial phase；NsetTo locate in advance
Manage sequence length.
Comb filtering processing needs to use unit period sequence length N2π4 times of sequence lengths.Comb filtering processing is in frequency
Domain amplitude versus frequency characte is as shown in Figure 2.
For step S108, in one embodiment, ratio integer is set to odd number, determines the comb filtering sequence length
It is formula (8) with the ratio integer of the unit period sequence length：
Wherein, k be the ratio integer, NsetFor preprocessing sequence length, NLFor single-stage comb filtering parameter, N2πFor list
Bit period sequence length.
In one embodiment, predetermined sequence length is arranged to odd number, and predetermined sequence length computation is formula (9)：
Wherein, N be predetermined sequence length, k be the ratio integer, N2πFor unit periodic sequence length.
For step S109, in one embodiment, default starting point can be the 0.5 of the unit period sequence length
Times, the described first positive sequence, is formula (10)：
Wherein, X8L(n) it is 8 grades of comb filtering sequences, X+start(n) it is the first positive sequence, PstartTo preset starting point,
N2πFor unit periodic sequence length, (int) is round numbers, and A is signal amplitude, unit v, ωiFor signal frequency, KL(ωi) it is comb
Shape filtering process is in signal frequency ωiDimensionless amplitude gain, T is sampling interval duration, and n is series of discrete number,For
One positive sequence initial phase, N are predetermined sequence length.
In one embodiment, the described first anti-pleat sequence is formula (11)：
Wherein, X-start(- n) be the first anti-pleat sequence, X+start(n) it is the first positive sequence, A is signal amplitude, unit
V, ωiFor signal frequency, T is sampling interval duration, and n is series of discrete number, and β 1 is the first anti-pleat sequence initial phase, and N is default
Sequence length.Comb filtering sequence, the first positive sequence and the first anti-pleat sequence pattern expression are as shown in Figure 3.
For step S110, the calculating of the first positive phase and the first antiphase is the knot based on orthogonal mixing and integral and calculating
Fruit.
When without considering the mixing interfering frequency of orthogonal mixing, orthogonal mixing is expressed as formula (12), and integral and calculating is expressed as
Formula (13)：
Wherein, R+start(n) it is the first positive real frequency mixing sequence, I+start(n) sequence, R are mixed for the first weakened body resistance frequency-start
(- n) is mixed sequence, I for the first anti-real frequency-start(- n) is mixed sequence, cos (ω for the first anti-empty frequencys) or cos (- ω TnsTn)
For the discrete cosine function of reference frequency, sin (ωs) or sin (- ω TnsTn it is) discrete sine function of reference frequency, Ω is
Signal frequency ωiWith reference frequency ωsFrequency difference, T is sampling interval duration, and n is series of discrete number, KL(ωi) it is comb filtering
Processing is in signal frequency ωiDimensionless amplitude gain,For the first positive sequence initial phase, β 1 is at the beginning of the first positive sequence
Phase, N are predetermined sequence length.
Wherein, R+startFor the first positive real frequency integrated value, unit dimensionless, I+startFor the first weakened body resistance frequency integrated value, unit
Dimensionless, R-startFor the first anti-real frequency integrated value, unit dimensionless, I-startFor the first anti-empty frequency mixing integrated value, unit without
Dimension, Ω are signal frequency ωiWith reference frequency ωsFrequency difference, T is sampling interval duration, and n is series of discrete number, and N is default
Sequence length,For the first positive sequence initial phase, β 1 is the first anti-pleat sequence initial phase, and N is predetermined sequence length.
In one embodiment, the calculating of the first positive phase and the first antiphase is expressed as formula (14)：
Wherein, PH+startFor the first positive phase, PH-startFor the first antiphase, R+startIt is single for the first positive real frequency integrated value
Position dimensionless, I+startFor the first weakened body resistance frequency integrated value, unit dimensionless, R-startFor first, anti-reality frequency integrated value, unit are immeasurable
Guiding principle, I-startFor the first anti-empty frequency mixing integrated value, unit dimensionless, Ω is signal frequency ωiWith reference frequency ωsFrequency difference, T
For sampling interval duration, N is predetermined sequence length,For the first positive sequence initial phase, β 1 is the first anti-pleat sequence first phase
Position.
For step S111, in one embodiment, the first average initial phase computational methods are expressed as formula (15)：
Wherein, PHstart-avgFor the first average initial phase, PH+startFor the first positive phase, PH-startFor the first antiphase,For the first positive sequence initial phase, β 1 is the first anti-pleat sequence initial phase.
For step S112, in one embodiment, the preset phase value can be ± π/4；It is flat according to described first
The step of equal initial phase and preset phase value, acquisition phase fiducial value, can include：
If the described first average initial phase is more than or equal to 0 and is less than or equal to pi/2, the described first average first phase is subtracted according to π/4
Position obtains phase fiducial value；
If the described first average initial phase is more than or equal to-pi/2 is less than or equal to 0, according to-π/4 subtract described first it is average just
Phase obtains phase fiducial value.
Specially formula (16)：
Wherein, △ PHcomFor phase fiducial value, unit rad, PHstart-avgFor the first average initial phase.
In one embodiment, the new starting point calculates, and is formula (17)：
Wherein, PnewFor new starting point, unit dimensionless, PstartTo preset starting point, △ PHcomIt is single for phase fiducial value
Position rad, N2πFor unit periodic sequence length, (int) is round numbers.
For step S113, the second positive sequence and the second anti-pleat sequence are formula (18)：
Wherein, X8L(n) it is 8 grades of comb filtering sequences, X+end(n) it is the second positive sequence, X-end(- n) is the second anti-pleat
Sequence, PnewFor new starting point, unit dimensionless, KL(ωi) handled for comb filtering in signal frequency ωiDimensionless amplitude increase
Benefit,For the second positive sequence initial phase, β 2 is the second anti-pleat sequence initial phase, ωiFor signal frequency, when T is the sampling interval
Between, n is series of discrete number, and N is predetermined sequence length.
For step S114, the computational methods of the second positive phase and the second antiphase are based on orthogonal mixing and number filter
The result that ripple calculates.The digital filtering is made of 6 grades of rectangular window arithmetic mean filters of 2 kinds of filtering parameters.
When without considering the mixing interfering frequency of orthogonal mixing, orthogonal mixing is expressed as formula (19), and the 6 of 2 kinds of filtering parameters
Grade rectangular window arithmetic mean filter filtering calculation expression is formula (20)：
Wherein, R+end(n) it is the second positive real frequency mixing sequence, I+end(n) sequence, R are mixed for the second weakened body resistance frequency-end(-n)
For the second anti-real frequency mixing sequence, I-end(- n) is mixed sequence, cos (ω for the second anti-empty frequencys) or cos (- ω TnsTn it is) ginseng
Examine the discrete cosine function of frequency, sin (ωs) or sin (- ω TnsTn it is) discrete sine function of reference frequency, KL(ωi) be
Comb filtering processing is in signal frequency ωiDimensionless amplitude gain, Ω be signal frequency ωiWith reference frequency ωsFrequency difference,
ωiFor signal frequency, T is sampling interval duration, and n is series of discrete number,For the first positive sequence initial phase, β 1 is first
Anti- pleat sequence initial phase, N are predetermined sequence length.
Wherein, R+endFor the second positive real frequency digital filtering final value, unit dimensionless；I+endFor the second weakened body resistance frequency digital filtering
Final value, unit dimensionless；R-endFor the second anti-digital filtering final value, unit dimensionless；I-endIt is whole for the second anti-empty frequency digital filtering
Value, unit dimensionless；KL(ωi) handled for comb filtering in signal frequency ωiDimensionless amplitude gain, Ω is signal frequency
ωiWith reference frequency ωsFrequency difference；K (Ω) is digital filtering in the amplitude gain of frequency difference Ω, unit dimensionless；T is between sampling
Every the time；For the second positive sequence initial phase；β 2 is the second anti-pleat sequence initial phase；ND1For filtering parameter 1, i.e., to ND1It is a
Continuous centrifugal pump is added, its arithmetic mean of instantaneous value is then taken to be exported as this filter value；ND2For filtering parameter 2, i.e., to ND2A company
Continuous centrifugal pump is added, its arithmetic mean of instantaneous value is then taken to be exported as this filter value；NDSequence length, number are used for digital filtering
It is the summation of 6 grades of rectangular window arithmetic mean filter filtering parameters in amount, less than or equal to predetermined sequence length N.
In one embodiment, filtering parameter ND1Value is 1.5 times of the unit period sequence length of the reference frequency,
Purpose carries out depth inhibition to the mixing interfering frequency that 1/3 subharmonic generates；Filtering parameter ND2Value is the reference frequency
2 times of unit period sequence length, purpose carries out depth to the mixing interfering frequency of the generations such as direct current, 1/2 gradation, subharmonic
Inhibit.6 grades of rectangular window arithmetic mean filters filtering of 2 kinds of filtering parameters, which calculates, to be needed using signal period sequence length
10.5 times.
Filtering parameter ND1With filtering parameter ND2It is calculated as formula (21)：
Wherein, ND1For digital filter parameters 1, unit dimensionless, (int) is round numbers, ND2It is single for digital filter parameters 2
Position dimensionless, N2πFor unit periodic sequence length.
In one embodiment, the computational methods of the second positive phase and the second antiphase are expressed as formula (22)：
Wherein, PH+endFor the second positive phase, PH-endFor the second antiphase, R+endFor the second positive real frequency integrated value, unit without
Dimension, I+endFor the second weakened body resistance frequency integrated value, unit dimensionless, R-endFor the second anti-real frequency integrated value, unit dimensionless, I-end
For the second anti-empty frequency mixing integrated value, unit dimensionless, Ω is signal frequency ωiWith reference frequency ωsFrequency difference, T for sampling between
Every the time, NDSequence length is used for digital filtering,For the second positive sequence initial phase, β 2 is the second anti-pleat sequence first phase
Position.
For step S115, the second average initial phase computational methods are expressed as formula (23)：
Wherein, PHend-avgFor the second average initial phase, PH+endFor the second positive phase, PH-endFor the second antiphase,For
Second positive sequence initial phase, β 2 are the second positive sequence initial phase.
For step S116, SIN function modulation sequence is expressed as formula (24)：
Wherein, Xsin(n) it is SIN function modulation sequence, AKL(ωi) for SIN function modulation sequence amplitude, unit v,For cosine function modulation sequence initial phase, KL(ωi) handled for comb filtering in signal frequency ωiDimensionless
Amplitude gain, ωiSignal frequency, T are sampling interval duration, and n is series of discrete number, and N is predetermined sequence length,For second
Positive sequence initial phase, β 2 are the second anti-pleat sequence initial phase.
For step S117, in one embodiment, zero initial phase reference sinusoidal FUNCTION MODULATION sequence is obtained, is expressed as formula
(25)：
Wherein, X0sin(n) it is zero initial phase reference sinusoidal FUNCTION MODULATION sequence, AKL(ωi) it is SIN function modulation sequence
Amplitude, unit v, ωiFor signal frequency, KL(ωi) handled for comb filtering in signal frequency ωiDimensionless amplitude gain,
T is sampling interval duration, and n is series of discrete number, and N is predetermined sequence length.Zero initial phase datum mark avatars such as Fig. 4 institutes
Show.
For step S118, in one embodiment, obtain real frequency mixing sequence and empty frequency is mixed sequence, be expressed as formula
(26)：
Wherein, R (n) is mixed sequence for real frequency, and I (n) is mixed sequence, X0 for empty frequencysin(n) it is zero initial phase reference sinusoidal
FUNCTION MODULATION sequence, cos (ωsTn it is) the discrete cosine function of reference frequency, sin (ωsTn it is) discrete sine of reference frequency
Function, U are public amplitude, and unit v, Ω are signal frequency ωiWith reference frequency ωsFrequency difference, T is sampling interval duration, and n is
Series of discrete number, 0.5 (N-1) are mixing sequence length.Its middle frequency difference Ω is partly useful component, frequency (ωi+ωs) be partly
It is mixed interference component.
For step S119, the filtering series of digital filtering is set, in one embodiment, rectangular window arithmetic can be passed through
Average value filtering algorithm carries out digital filtering to the real frequency mixing sequence and the empty frequency mixing sequence respectively.
In one embodiment, the filtering series of digital filtering could be provided as 5 etc., and digital filter parameters are the unit cycle
1 times of sequence length, it is therefore an objective to which depth inhibition is carried out to subharmonic mixing interfering frequency.Without considering mixing interfering frequency, real frequency
Filtered sequence and imaginary frequency filtering sequence are expressed as formula (27)：
Wherein, RD(n) it is real frequency filtered sequence；ID(n) it is imaginary frequency filtering sequence；R (n) is mixed sequence for real frequency；I (n) is
Empty frequency mixing sequence；U is public amplitude, unit v；NDFor the filtering parameter of rectangular window arithmetic mean of instantaneous value filtering algorithm, i.e., to NDIt is a
Continuous centrifugal pump is added, its arithmetic mean of instantaneous value is then taken to be exported as this filter value；T is sampling interval duration；N2πFor unit
Periodic sequence length；KD(Ω) is digital filtering in the amplitude gain of frequency difference Ω, unit dimensionless；α (Ω) exists for digital filtering
The phase shift of frequency difference Ω；M is digital filtering output sequence length；0.5 (N-1) is digital filtering list entries length；N is default sequence
Row length.
For step S120, in one embodiment, the digital filtering output sequence length is integration lengths, is obtained real
Frequency integrated value and empty frequency integrated value, are expressed as formula (28)：
Wherein, RDFor real frequency integrated value；IDFor empty frequency integrated value；RD(n) it is real frequency filtered sequence；ID(n) it is imaginary frequency filtering
Sequence；T is sampling interval duration, and M is integration lengths；NDFor the filtering parameter of rectangular window arithmetic mean of instantaneous value filtering algorithm；0.5
(N-1) it is digital filtering list entries length；N is predetermined sequence length.
For step S121, in one embodiment, the phase of zero initial phase reference sinusoidal FUNCTION MODULATION sequence, table are obtained
Up to for formula (29)：
Wherein, PH is the phase of zero initial phase reference sinusoidal FUNCTION MODULATION sequence, IDFor empty frequency integrated value, RDFor real frequency product
Score value, M are integration lengths；NDFor the filtering parameter of rectangular window arithmetic mean of instantaneous value filtering algorithm.
In one embodiment, according to the phase of zero initial phase reference sinusoidal FUNCTION MODULATION sequence, electric power signal frequency is obtained
Rate is expressed as formula (30)：
Wherein, ωiFor the frequency of electric power signal, PH is the phase of zero initial phase reference sinusoidal FUNCTION MODULATION sequence, and T is to adopt
Sample interval time, M are integration lengths, NDFor the filtering parameter of rectangular window arithmetic mean of instantaneous value filtering algorithm, ωsFor reference frequency.
Based on same inventive concept, the present invention also provides it is a kind of according to the zero initial phase reference sinusoidal sequence of function into line frequency
Measuring system is below in conjunction with the accompanying drawings described in detail embodiments of systems of the invention.
As shown in figure 5, a kind of carry out frequency measuring system according to the zero initial phase reference sinusoidal sequence of function, including：
Preliminary sequence length determination modul 101, for according to frequency power signal scope lower limit, preset sample frequency and
Default integer signal period number, obtains preliminary sequence length；
Preliminary sequence acquisition module 102 for being sampled according to the preliminary sequence length to electric power signal, obtains institute
State the preliminary sequence of electric power signal；
Reference frequency determining module 103, for, into line frequency preliminary survey, obtaining the electric power signal to the preliminary sequence
Preliminary frequency obtains reference frequency according to the preliminary frequency；
Unit period sequence length determining module 104, for according to the preset sample frequency and the reference frequency, obtaining
Obtain the unit period sequence length of the electric power signal；
Preprocessing sequence length determination modul 105, for by the default integer signal period number and the unit period
Sequence length is multiplied, and obtains preprocessing sequence length；
Preprocessing sequence acquisition module 106, for according to the preprocessing sequence length, from the preliminary of the electric power signal
Preprocessing sequence is obtained in sequence；
Comb filtering sequence determining module 107 for carrying out comb filtering processing to the preprocessing sequence, obtains pectination
Filtered sequence, wherein comb filtering sequence length are carrying out comb filtering treated residue length for the preprocessing sequence；
Predetermined sequence length determination modul 108, for determining the comb filtering sequence length and the unit period sequence
The ratio integer of row length obtains predetermined sequence length according to the ratio integer and the unit period sequence length, wherein
The ratio integer is odd number, and the predetermined sequence length is odd number；
First ray acquisition module 109, for according to the predetermined sequence length and default starting point, being filtered from the pectination
The first positive sequence is obtained in wave train, the first anti-pleat sequence is obtained according to the described first positive sequence；
First positive and negative phase determination module 110, for obtaining the first positive phase according to the described first positive sequence, according to institute
It states the first anti-pleat sequence and obtains the first antiphase；
First average initial phase determining module 111, for according to first positive phase and first antiphase acquisition
First average initial phase；
New starting point determining module 112, for according to the described first average initial phase and preset phase value, obtaining phase ratio
Compared with value, according to the phase fiducial value, the default starting point and the unit period sequence length, new starting point is obtained；
Second retrieval module 113, for according to the predetermined sequence length and the new starting point, from the pectination
The second positive sequence is obtained in filtered sequence, the second anti-pleat sequence is obtained according to the second positive sequence；
Second positive and negative phase determination module 114, for obtaining the second positive phase according to the described second positive sequence, according to institute
It states the second anti-pleat sequence and obtains the second antiphase；
Second average initial phase determining module 115, for according to second positive phase and second antiphase acquisition
Second average initial phase；
SIN function modulation sequence determining module 116, for by the described second positive sequence and the described second anti-pleat sequence
Subtract each other, obtain difference sequence, according to the difference sequence and the sine function of the second average initial phase, obtain SIN function tune
Sequence processed；
Zero initial phase modulation sequence acquisition module 117 for being exported from the SIN function modulation sequence central point, obtains
Zero initial phase reference sinusoidal FUNCTION MODULATION sequence；
Sequence determining module 118 is mixed, for by the zero initial phase reference sinusoidal FUNCTION MODULATION sequence and reference frequency
Discrete cosine function be multiplied, obtain real frequency mixing sequence, will the zero initial phase reference sinusoidal FUNCTION MODULATION sequence with referring to
The discrete sine function of frequency is multiplied, and obtains empty frequency mixing sequence；
Filtered sequence determining module 119, for being carried out respectively to the real frequency mixing sequence and the empty frequency mixing sequence
Digital filtering obtains real frequency filtered sequence and imaginary frequency filtering sequence；
Integrated value determining module 120, for being accumulated respectively to the real frequency filtered sequence and the imaginary frequency filtering sequence
Point, obtain real frequency integrated value and empty frequency integrated value；
Frequency power signal determining module 121, for according to the real frequency integrated value and empty frequency integrated value, obtaining zero first phase
The phase of position reference sinusoidal FUNCTION MODULATION sequence obtains the frequency of electric power signal according to the phase.
In one embodiment, the mixing sequence determining module 118 can be according to expression formula R (n)=X0sin(n)cos
(ωsTn real frequency mixing sequence R (n)) is obtained, wherein, X0sin(n) it is zero initial phase reference sinusoidal FUNCTION MODULATION sequence, cos (ωsTn it is) the discrete cosine function of reference frequency, ωsFor reference frequency, T is sampling interval duration, and n is series of discrete number,
The mixing sequence determining module 118 can be according to expression formula I (n)=X0sin(n)sin(ωsTn empty frequency) is obtained
Sequence I (n) is mixed, wherein, X0sin(n) it is zero initial phase reference sinusoidal FUNCTION MODULATION sequence, sin (ωsTn) it is reference frequency
Discrete sine function, ωsFor reference frequency, T is sampling interval duration, and n is series of discrete number,
In one embodiment, the filtered sequence determining module 119 can be calculated by rectangular window arithmetic average value filtering
Method carries out digital filtering to the real frequency mixing sequence and the empty frequency mixing sequence respectively.
In one embodiment, the frequency power signal determining module 121 can be according to expression formula
The phase PH of zero initial phase reference sinusoidal FUNCTION MODULATION sequence is obtained, wherein, RDFor real frequency integrated value, IDFor empty frequency integrated value.
In one embodiment, the frequency power signal determining module 121 can be according to expression formulaObtain the frequencies omega of electric power signali, wherein, PH is zero initial phase reference sinusoidal FUNCTION MODULATION
The phase of sequence, T are sampling interval duration, and M is integration lengths, NDJoin for the filtering of rectangular window arithmetic mean of instantaneous value filtering algorithm
Number, ωsFor reference frequency.
The other technical characteristics of present system are identical with the method for the present invention, and it will not be described here.
Each technical characteristic of embodiment described above can be combined arbitrarily, to make description succinct, not to above-mentioned reality
It applies all possible combination of each technical characteristic in example to be all described, as long as however, the combination of these technical characteristics is not deposited
In contradiction, the scope that this specification is recorded all is considered to be.
Embodiment described above only expresses the several embodiments of the present invention, and description is more specific and detailed, but simultaneously
It cannot therefore be construed as limiting the scope of the patent.It should be pointed out that come for those of ordinary skill in the art
It says, without departing from the inventive concept of the premise, various modifications and improvements can be made, these belong to the protection of the present invention
Scope.Therefore, the protection domain of patent of the present invention should be determined by the appended claims.
Claims (10)
1. a kind of carry out frequency measurement method according to the zero initial phase reference sinusoidal sequence of function, which is characterized in that including step：
According to the lower limit of frequency power signal scope, preset sample frequency and default integer signal period number, preliminary sequence is obtained
Length；
Electric power signal is sampled according to the preliminary sequence length, obtains the preliminary sequence of the electric power signal；
To the preliminary sequence into line frequency preliminary survey, the preliminary frequency of the electric power signal is obtained, is obtained according to the preliminary frequency
To reference frequency；
According to the preset sample frequency and the reference frequency, the unit period sequence length of the electric power signal is obtained；
The default integer signal period number with the unit period sequence length is multiplied, obtains preprocessing sequence length；
According to the preprocessing sequence length, preprocessing sequence is obtained from the preliminary sequence of the electric power signal；
Comb filtering processing is carried out to the preprocessing sequence, obtains comb filtering sequence, wherein comb filtering sequence length is
The preprocessing sequence is carrying out comb filtering treated residue length；
The ratio integer of the comb filtering sequence length and the unit period sequence length is determined, according to the ratio integer
Predetermined sequence length is obtained with the unit period sequence length, wherein the ratio integer is odd number, the predetermined sequence is long
It spends for odd number；
According to the predetermined sequence length and default starting point, the first positive sequence, root are obtained from the comb filtering sequence
The first anti-pleat sequence is obtained according to the described first positive sequence；
First positive phase is obtained according to the described first positive sequence, the first antiphase is obtained according to the described first anti-pleat sequence；
First average initial phase is obtained according to first positive phase and first antiphase；
According to the described first average initial phase and preset phase value, phase fiducial value is obtained, according to the phase fiducial value, described
Default starting point and the unit period sequence length, obtain new starting point；
According to the predetermined sequence length and the new starting point, the second positive sequence is obtained from the comb filtering sequence,
Second anti-pleat sequence is obtained according to the second positive sequence；
Second positive phase is obtained according to the described second positive sequence, the second antiphase is obtained according to the described second anti-pleat sequence；
Second average initial phase is obtained according to second positive phase and second antiphase；
Described second positive sequence with the described second anti-pleat sequence is subtracted each other, obtains difference sequence, according to the difference sequence and described
The sine function of second average initial phase, obtains SIN function modulation sequence；
According to expression formulaObtain SIN function modulation sequence Xsin(n), wherein,
X+end(n) it is the second positive sequence, X-end(- n) be the second anti-pleat sequence, PHend-avgFor the second average initial phase；
It is exported from the SIN function modulation sequence central point, obtains zero initial phase reference sinusoidal FUNCTION MODULATION sequence；
The zero initial phase reference sinusoidal FUNCTION MODULATION sequence with the discrete cosine function of reference frequency is multiplied, it is mixed to obtain real frequency
The zero initial phase reference sinusoidal FUNCTION MODULATION sequence is multiplied by frequency sequence with the discrete sine function of reference frequency, obtains empty
Frequency mixing sequence；
Digital filtering is carried out to the real frequency mixing sequence and the empty frequency mixing sequence respectively, obtains real frequency filtered sequence and void
Frequency filtered sequence；
The real frequency filtered sequence and the imaginary frequency filtering sequence are integrated respectively, real frequency integrated value is obtained and empty frequency integrates
Value；
According to the real frequency integrated value and empty frequency integrated value, the phase of zero initial phase reference sinusoidal FUNCTION MODULATION sequence of acquisition, root
The frequency of electric power signal is obtained according to the phase.
2. according to claim 1 carry out frequency measurement method, feature according to the zero initial phase reference sinusoidal sequence of function
It is, according to expression formulaThe phase PH of zero initial phase reference sinusoidal FUNCTION MODULATION sequence is obtained, wherein,
RDFor real frequency integrated value, IDFor empty frequency integrated value.
3. according to claim 1 carry out frequency measurement method, feature according to the zero initial phase reference sinusoidal sequence of function
It is, the real frequency mixing sequence and the empty frequency mixing sequence is carried out respectively by rectangular window arithmetic mean of instantaneous value filtering algorithm
Digital filtering.
4. according to claim 3 carry out frequency measurement method, feature according to the zero initial phase reference sinusoidal sequence of function
It is, according to expression formulaObtain the frequencies omega of electric power signali, wherein, PH is zero initial phase base
The phase of quasi sinusoids modulation sequence, T are sampling interval duration, and M is integration lengths, NDFor rectangular window arithmetic average value filtering
The filtering parameter of algorithm, ωsFor reference frequency.
5. according to any one of claims 1 to 4, carry out frequency measurement according to the zero initial phase reference sinusoidal sequence of function
Method, it is characterised in that：
According to expression formula R (n)=X0sin(n)cos(ωsTn real frequency mixing sequence R (n)) is obtained, wherein, X0sin(n) it is zero first phase
Position reference sinusoidal FUNCTION MODULATION sequence, cos (ωsTn it is) the discrete cosine function of reference frequency, ωsFor reference frequency, T is to adopt
Sample interval time, n are series of discrete number,
According to expression formula I (n)=X0sin(n)sin(ωsTn empty frequency mixing sequence I (n)) is obtained, wherein, X0sin(n) it is zero first phase
Position reference sinusoidal FUNCTION MODULATION sequence, sin (ωsTn it is) discrete sine function of reference frequency, ωsFor reference frequency, T is to adopt
Sample interval time, n are series of discrete number,
6. a kind of carry out frequency measuring system according to the zero initial phase reference sinusoidal sequence of function, which is characterized in that including：
Preliminary sequence length determination modul, for lower limit, the preset sample frequency and default whole according to frequency power signal scope
Number signal period number, obtains preliminary sequence length；
Preliminary sequence acquisition module for being sampled according to the preliminary sequence length to electric power signal, obtains the electric power
The preliminary sequence of signal；
Reference frequency determining module, for, into line frequency preliminary survey, obtaining the first cadence of the electric power signal to the preliminary sequence
Rate obtains reference frequency according to the preliminary frequency；
Unit period sequence length determining module, for according to the preset sample frequency and the reference frequency, described in acquisition
The unit period sequence length of electric power signal；
Preprocessing sequence length determination modul, for by the default integer signal period number and the unit period sequence length
It is multiplied, obtains preprocessing sequence length；
Preprocessing sequence acquisition module, for according to the preprocessing sequence length, from the preliminary sequence of the electric power signal
Obtain preprocessing sequence；
Comb filtering sequence determining module for carrying out comb filtering processing to the preprocessing sequence, obtains comb filtering sequence
Row, wherein comb filtering sequence length are carrying out comb filtering treated residue length for the preprocessing sequence；
Predetermined sequence length determination modul, for determining the comb filtering sequence length and the unit period sequence length
Ratio integer obtains predetermined sequence length, wherein the ratio according to the ratio integer and the unit period sequence length
Integer is odd number, and the predetermined sequence length is odd number；
First ray acquisition module, for according to the predetermined sequence length and default starting point, from the comb filtering sequence
It is middle to obtain the first positive sequence, the first anti-pleat sequence is obtained according to the described first positive sequence；
First positive and negative phase determination module, for obtaining the first positive phase according to the described first positive sequence, according to described first
Anti- pleat sequence obtains the first antiphase；
First average initial phase determining module, it is average for obtaining first according to first positive phase and first antiphase
Initial phase；
New starting point determining module, for according to the described first average initial phase and preset phase value, obtaining phase fiducial value, root
According to the phase fiducial value, the default starting point and the unit period sequence length, new starting point is obtained；
Second retrieval module, for according to the predetermined sequence length and the new starting point, from the comb filtering sequence
The second positive sequence is obtained in row, the second anti-pleat sequence is obtained according to the second positive sequence；
Second positive and negative phase determination module, for obtaining the second positive phase according to the described second positive sequence, according to described second
Anti- pleat sequence obtains the second antiphase；
Second average initial phase determining module, it is average for obtaining second according to second positive phase and second antiphase
Initial phase；
SIN function modulation sequence determining module for the described second positive sequence to be subtracted each other with the described second anti-pleat sequence, obtains
To difference sequence, according to the difference sequence and the sine function of the second average initial phase, SIN function modulation sequence is obtained；
According to expression formulaObtain SIN function modulation sequence Xsin(n), wherein,
X+end(n) it is the second positive sequence, X-end(- n) be the second anti-pleat sequence, PHend-avgFor the second average initial phase；
Zero initial phase modulation sequence acquisition module for being exported from the SIN function modulation sequence central point, obtains zero first phase
Position reference sinusoidal FUNCTION MODULATION sequence；
Sequence determining module is mixed, for by the discrete remaining of the zero initial phase reference sinusoidal FUNCTION MODULATION sequence and reference frequency
String function is multiplied, and obtains real frequency mixing sequence, by the zero initial phase reference sinusoidal FUNCTION MODULATION sequence and reference frequency from
It dissipates SIN function to be multiplied, obtains empty frequency mixing sequence；
Filtered sequence determining module, for carrying out digital filter to the real frequency mixing sequence and the empty frequency mixing sequence respectively
Ripple obtains real frequency filtered sequence and imaginary frequency filtering sequence；
Integrated value determining module for being integrated respectively to the real frequency filtered sequence and the imaginary frequency filtering sequence, obtains
Real frequency integrated value and empty frequency integrated value；
Frequency power signal determining module, for according to the real frequency integrated value and empty frequency integrated value, obtaining zero initial phase benchmark
The phase of SIN function modulation sequence obtains the frequency of electric power signal according to the phase.
7. according to claim 6 carry out frequency measuring system, feature according to the zero initial phase reference sinusoidal sequence of function
It is, the frequency power signal determining module is according to expression formulaObtain zero initial phase reference sinusoidal letter
The phase PH of number modulation sequence, wherein, RDFor real frequency integrated value, IDFor empty frequency integrated value.
8. according to claim 6 carry out frequency measuring system, feature according to the zero initial phase reference sinusoidal sequence of function
Be, the filtered sequence determining module by rectangular window arithmetic mean of instantaneous value filtering algorithm respectively to the real frequency mixing sequence and
The empty frequency mixing sequence carries out digital filtering.
9. according to claim 8 carry out frequency measuring system, feature according to the zero initial phase reference sinusoidal sequence of function
It is, the frequency power signal determining module is according to expression formulaObtain the frequency of electric power signal
Rate ωi, wherein, PH is the phase of zero initial phase reference sinusoidal FUNCTION MODULATION sequence, and T is sampling interval duration, and M grows for integration
Degree, NDFor the filtering parameter of rectangular window arithmetic mean of instantaneous value filtering algorithm, ωsFor reference frequency.
10. frequency measurement is carried out according to the zero initial phase reference sinusoidal sequence of function according to claim 6 to 9 any one
System, it is characterised in that：
The mixing sequence determining module is according to expression formula R (n)=X0sin(n)cos(ωsTn real frequency mixing sequence R (n)) is obtained,
Wherein, X0sin(n) it is zero initial phase reference sinusoidal FUNCTION MODULATION sequence, cos (ωsTn it is) the discrete cosine letter of reference frequency
Number, ωsFor reference frequency, T is sampling interval duration, and n is series of discrete number,
The mixing sequence determining module is according to expression formula I (n)=X0sin(n)sin(ωsTn empty frequency mixing sequence I (n)) is obtained,
Wherein, X0sin(n) it is zero initial phase reference sinusoidal FUNCTION MODULATION sequence, sin (ωsTn it is) the discrete sine letter of reference frequency
Number, ωsFor reference frequency, T is sampling interval duration, and n is series of discrete number,
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510891522.XA CN105548694B (en) | 2015-12-04 | 2015-12-04 | Frequency measurement method and system are carried out according to the zero initial phase reference sinusoidal sequence of function |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510891522.XA CN105548694B (en) | 2015-12-04 | 2015-12-04 | Frequency measurement method and system are carried out according to the zero initial phase reference sinusoidal sequence of function |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105548694A CN105548694A (en) | 2016-05-04 |
| CN105548694B true CN105548694B (en) | 2018-06-01 |
Family
ID=55828015
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510891522.XA Active CN105548694B (en) | 2015-12-04 | 2015-12-04 | Frequency measurement method and system are carried out according to the zero initial phase reference sinusoidal sequence of function |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105548694B (en) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3441871A (en) * | 1967-01-05 | 1969-04-29 | Gen Electric Co Ltd | Oscillator system selectively operable as frequency measuring device or automatic frequency control |
| US3675124A (en) * | 1970-10-28 | 1972-07-04 | Sperry Rand Corp | Apparatus for measuring frequency modulation noise signals and for calibrating same |
| DE3327876A1 (en) * | 1982-08-02 | 1984-03-01 | Raytheon Co., 02173 Lexington, Mass. | FREQUENCY MEASURING DEVICE |
| JP2015025726A (en) * | 2013-07-26 | 2015-02-05 | 東芝三菱電機産業システム株式会社 | Frequency detecting device, frequency detecting method, and electric power converter |
| CN104502706A (en) * | 2014-12-29 | 2015-04-08 | 广东电网有限责任公司电力科学研究院 | Harmonic amplitude measurement method and system of power signal |
| CN104502700A (en) * | 2014-12-29 | 2015-04-08 | 广东电网有限责任公司电力科学研究院 | Sinusoidal parameter measurement method and system of power signal |
| CN104635044A (en) * | 2015-02-05 | 2015-05-20 | 广东电网有限责任公司电力科学研究院 | Power signal frequency detection method and system based on amplitude modulation |
| CN104635045A (en) * | 2015-02-05 | 2015-05-20 | 广东电网有限责任公司电力科学研究院 | Power signal frequency detection method and system based on phase modulation |
| CN105067880A (en) * | 2015-09-18 | 2015-11-18 | 广东电网有限责任公司电力科学研究院 | Method and system for orthogonal modulation of power signals |
-
2015
- 2015-12-04 CN CN201510891522.XA patent/CN105548694B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3441871A (en) * | 1967-01-05 | 1969-04-29 | Gen Electric Co Ltd | Oscillator system selectively operable as frequency measuring device or automatic frequency control |
| US3675124A (en) * | 1970-10-28 | 1972-07-04 | Sperry Rand Corp | Apparatus for measuring frequency modulation noise signals and for calibrating same |
| DE3327876A1 (en) * | 1982-08-02 | 1984-03-01 | Raytheon Co., 02173 Lexington, Mass. | FREQUENCY MEASURING DEVICE |
| JP2015025726A (en) * | 2013-07-26 | 2015-02-05 | 東芝三菱電機産業システム株式会社 | Frequency detecting device, frequency detecting method, and electric power converter |
| CN104502706A (en) * | 2014-12-29 | 2015-04-08 | 广东电网有限责任公司电力科学研究院 | Harmonic amplitude measurement method and system of power signal |
| CN104502700A (en) * | 2014-12-29 | 2015-04-08 | 广东电网有限责任公司电力科学研究院 | Sinusoidal parameter measurement method and system of power signal |
| CN104635044A (en) * | 2015-02-05 | 2015-05-20 | 广东电网有限责任公司电力科学研究院 | Power signal frequency detection method and system based on amplitude modulation |
| CN104635045A (en) * | 2015-02-05 | 2015-05-20 | 广东电网有限责任公司电力科学研究院 | Power signal frequency detection method and system based on phase modulation |
| CN105067880A (en) * | 2015-09-18 | 2015-11-18 | 广东电网有限责任公司电力科学研究院 | Method and system for orthogonal modulation of power signals |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105548694A (en) | 2016-05-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104635044B (en) | Power signal frequency detection method and system based on amplitude modulation | |
| CN105067880B (en) | The method and system of orthogonal modulation is carried out to electric power signal | |
| CN104635045B (en) | Frequency power signal detection method and system based on phase-modulation | |
| CN104502700B (en) | Sinusoidal parameter measurement method and system of power signal | |
| CN105548694B (en) | Frequency measurement method and system are carried out according to the zero initial phase reference sinusoidal sequence of function | |
| CN105425031B (en) | The orthogonal frequency multiplication sequence method of zero initial phase benchmark and system are obtained from electric power signal | |
| CN105548690B (en) | Based on zero initial phase reference sinusoidal function frequency multiplication sequence frequency measuring method and system | |
| CN105548688B (en) | Frequency measurement method and system are carried out according to zero initial phase benchmark cosine function sequence | |
| CN105548693B (en) | Obtain the method and system of zero initial phase reference sinusoidal function times frequency sequence | |
| CN105425040B (en) | Based on zero initial phase benchmark cosine function frequency multiplication sequence frequency measuring method and system | |
| CN105548692B (en) | The method and system of frequency measurement is carried out according to zero initial phase benchmark orthogonal sequence | |
| CN105372490B (en) | The method and system of frequency measurement based on zero orthogonal times of frequency sequence of initial phase benchmark | |
| CN105388359B (en) | The method and system of orthogonal times of frequency sequence is obtained from electric power signal | |
| CN105403768B (en) | The method and system of cosine function times frequency sequence is obtained from electric power signal | |
| CN105301356B (en) | Obtain the method and system of zero initial phase benchmark cosine function times frequency sequence | |
| CN105425032B (en) | The method and system of SIN function times frequency sequence is obtained from electric power signal | |
| CN105067885B (en) | Electric power signal is converted to the method and system of zero initial phase signal sequence | |
| CN105548699B (en) | Obtain the method and system of the raising frequency reference SIN function sequence of electric power signal | |
| CN105203843B (en) | The average first phase position detecting method and system of electric power signal | |
| CN105425035B (en) | Obtain the method and system for dropping low-frequency zero initial phase reference sinusoidal sequence of function | |
| CN105548704B (en) | The method and system of fine setting frequency reference sequence is obtained from electric power signal | |
| CN105548691B (en) | Obtain the method and system of the raising frequency sine sequence of function of electric power signal | |
| CN105158560B (en) | The method and system of orthogonal reference signals sequence is extracted from electric power signal | |
| CN105548689B (en) | Obtain the method and system of the reduction frequency reference orthogonal sequence of electric power signal | |
| CN105467211B (en) | Obtain the method and system of the raising frequency reference cosine function sequence of electric power signal |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |