CN105548706B - Obtain the method and system of the reduction frequency sine sequence of function of electric power signal - Google Patents
Obtain the method and system of the reduction frequency sine sequence of function of electric power signal Download PDFInfo
- Publication number
- CN105548706B CN105548706B CN201510891736.7A CN201510891736A CN105548706B CN 105548706 B CN105548706 B CN 105548706B CN 201510891736 A CN201510891736 A CN 201510891736A CN 105548706 B CN105548706 B CN 105548706B
- Authority
- CN
- China
- Prior art keywords
- sequence
- electric power
- power signal
- frequency
- 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/16—Spectrum analysis; Fourier analysis
- G01R23/165—Spectrum analysis; Fourier analysis using filters
- G01R23/167—Spectrum analysis; Fourier analysis using filters with digital filters
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Testing Relating To Insulation (AREA)
- Inverter Devices (AREA)
Abstract
The invention discloses it is a kind of obtain electric power signal reduction the frequency sine sequence of function method and system, the method includes：Obtain electric power signal preliminary sequence length and preliminary sequence；To preliminary sequence into line frequency preliminary survey, reference frequency is set；Electric power signal unit period sequence length and predetermined sequence length are obtained, first positive/anti-pleat sequence is further obtained；Averagely initial phase is obtained for the first time according to first positive/anti-pleat sequence；Further according to first averagely initial phase and the first phase fiducial value of ± π/4 and new starting point obtain again positive/anti-pleat sequence, be averaged initial phase again；According to positive/anti-pleat sequence again, be averaged initial phase again, obtains remaining/SIN function modulation sequence；First/second multiplication sequence is obtained according to default fine tuning frequency, remaining/SIN function modulation sequence, finally being obtained according to first/second multiplication sequence reduces the frequency sine sequence of function.The present invention improves the accuracy that sine parameter calculates, and is suitble to application.
Description
Technical field
The present invention relates to technical field of electric power, more particularly to a kind of reduction frequency sine function sequence obtaining electric power signal
The method and system of row.
Background technology
Frequency measurement, phase measurement, the amplitude measurement etc. of electric system are the measurement of sine parameter in itself.Quickly
Fourier transformation and discrete Fourier transform are the basic skills realized sine parameter and measured, and have widely answer in the power system
With.But the above method is blocked in the non-integer of signal sampling process, causes spectrum leakage, and spectrum leakage will produce phase
The error answered.
In terms of electric system sine parameter measurement, there are many measurement methods, are repaiied such as adding window mouth functional based method, using interpolation
Correction method etc. reduces the influence of spectrum leakage problem.But above-mentioned measurement method measurement accuracy is low, is being not suitable for high accuracy just
The measurement of string parameter.
Invention content
Based on the above situation, the present invention proposes a kind of method for the reduction frequency sine sequence of function obtaining electric power signal
And system, improve the accuracy that sine parameter measures.
To achieve the goals above, the embodiment of technical solution of the present invention is：
A method of the reduction frequency sine sequence of function obtaining electric power signal includes the following steps：
According to the lower limit of frequency power signal range, preset sample frequency and default integer signal period number, described in acquisition
The preliminary sequence length of electric power signal；
The electric power signal is sampled according to the preliminary sequence length, obtains the preliminary sequence of the electric power signal
Row；
To the preliminary sequence into line frequency preliminary survey, the preliminary frequency of the electric power signal is generated, and according to described preliminary
Frequency sets the reference frequency of the electric power signal；
According to the preset sample frequency and the reference frequency, the unit period sequence for obtaining the electric power signal is long
Degree；
According to the default integer signal period number and the unit period sequence length, the pre- of the electric power signal is obtained
If sequence length；
According to default starting point and the predetermined sequence length, the first of the electric power signal is obtained from the preliminary sequence
Secondary forward direction sequence；
The first anti-pleat sequence of the electric power signal is obtained according to the first positive sequence；
The first positive phase of the electric power signal is obtained according to the first positive sequence, and according to the first anti-pleat sequence
Row obtain the first antiphase of the electric power signal；
The initial phase that is averaged for the first time of the electric power signal is obtained according to the first positive phase and the first antiphase；
The averagely initial phase for the first time is compared with ± π/4, obtains the first phase bit comparison compared with ± π/4
Value, and according to the first phase fiducial value and the default starting point, obtain new starting point；
According to the new starting point and the predetermined sequence length, the electric power signal is obtained from the preliminary sequence
Positive sequence again, and obtain according to the sequence positive again the pleat sequence anti-again of the electric power signal；
The positive phase again of the electric power signal is obtained according to the sequence positive again, and according to the pleat sequence anti-again
Row obtain the antiphase again of the electric power signal；
The initial phase that is averaged again of the electric power signal is obtained according to the positive phase again and the antiphase again；
The sequence of forward direction again is added with the pleat sequence anti-again, and according to result after being added and it is described again
Average initial phase, obtains the cosine function modulation sequence of the electric power signal；
The sequence positive again and the pleat sequence anti-again are subtracted each other, and according to after subtracting each other result and it is described again
Average initial phase, obtains the SIN function modulation sequence of the electric power signal；
The discrete sine function of default fine tuning frequency is multiplied with the cosine function modulation sequence and obtains the electric power letter
Number the first multiplication sequence, the discrete cosine function of the default fine tuning frequency is multiplied with the SIN function modulation sequence and obtains
Obtain the second multiplication sequence of the electric power signal；
Second multiplication sequence and first multiplication sequence are subtracted each other, obtaining the reduction frequency of the electric power signal just
The string sequence of function.
A kind of system for the reduction frequency sine sequence of function obtaining electric power signal, including：
Preliminary sequence length modules are used for according to the lower limit of frequency power signal range, preset sample frequency and preset whole
Number signal period number, obtains the preliminary sequence length of the electric power signal；
Preliminary sequence module, for being sampled to the electric power signal according to the preliminary sequence length, described in acquisition
The preliminary sequence of electric power signal；
Frequency preliminary survey module, for, into line frequency preliminary survey, generating the first cadence of the electric power signal to the preliminary sequence
Rate, and set according to the preliminary frequency reference frequency of the electric power signal；
Unit period sequence length module, for according to the preset sample frequency and the reference frequency, described in acquisition
The unit period sequence length of electric power signal；
Predetermined sequence length modules, for long according to the default integer signal period number and the unit period sequence
Degree, obtains the predetermined sequence length of the electric power signal；
First forward direction block, presets starting point and the predetermined sequence length, from the preliminary sequence for basis
The middle first positive sequence for obtaining the electric power signal；
First anti-pleat block, the first anti-pleat sequence for obtaining the electric power signal according to the first positive sequence
Row；
First phase module, the first positive phase for obtaining the electric power signal according to the first positive sequence, and
The first antiphase of the electric power signal is obtained according to the first anti-pleat sequence；
Be averaged initial phase module for the first time, for obtaining the electric power according to the first positive phase and the first antiphase
The initial phase that is averaged for the first time of signal；
Phase bits comparison module obtains and ± π/4 for the averagely initial phase for the first time to be compared with ± π/4
The first phase fiducial value compared, and according to the first phase fiducial value and the default starting point, obtain new starting point；
Block again is used for according to the new starting point and the predetermined sequence length, from the preliminary sequence
The sequence positive again of the electric power signal is obtained, and the anti-again of the electric power signal is obtained according to the sequence positive again
Pleat sequence；
Phase module again, the positive phase again for obtaining the electric power signal according to the sequence of forward direction again, and
The antiphase again of the electric power signal is obtained according to the pleat sequence anti-again；
Be averaged initial phase module again, for positive phase and the antiphase again to obtain the electric power again according to
The initial phase that is averaged again of signal；
Cosine function modulation sequence module, for the sequence of forward direction again to be added with the pleat sequence anti-again, and
According to result after being added and the initial phase that is averaged again, the cosine function modulation sequence of the electric power signal is obtained；
SIN function modulation sequence module, for the sequence positive again and the pleat sequence anti-again to be subtracted each other, and
According to after subtracting each other result and it is described again be averaged initial phase, obtain the SIN function modulation sequence of the electric power signal；
Multiplication sequence module is used for the discrete sine function of default fine tuning frequency and the cosine function modulation sequence phase
Multiply the first multiplication sequence for obtaining the electric power signal, by the discrete cosine function of the default fine tuning frequency and the sinusoidal letter
Number modulation sequence, which is multiplied, obtains the second multiplication sequence of the electric power signal；
Frequency sine sequence of function module is reduced, is used for second multiplication sequence and the first multiplication sequence phase
Subtract, obtains the reduction frequency sine sequence of function of the electric power signal.
Compared with prior art, beneficial effects of the present invention are：The present invention obtains the reduction frequency sine letter of electric power signal
The method and system of Number Sequence, according to the lower limit of frequency power signal range, preset sample frequency and preliminary design integer signal period
Number obtains preliminary sequence length, and is tentatively sampled to electric power signal, obtains preliminary sequence；To preliminary sequence into line frequency
Preliminary survey generates preliminary frequency, sets reference frequency；According to preset sample frequency and reference frequency, it is long to obtain unit period sequence
Degree；According to default integer signal period number and unit period sequence length, predetermined sequence length is obtained；It is obtained from preliminary sequence
First forward direction sequence, further obtains first anti-pleat sequence；First positive is obtained according to first positive sequence, first anti-pleat sequence
Position, first antiphase further obtain for the first time averagely initial phase；It is compared with ± π/4 and obtains first phase fiducial value, into one
Step obtains new starting point；According to new starting point and predetermined sequence length, positive sequence, again anti-pleat sequence again are obtained, into one
Step obtain again positive phase, again antiphase, again be averaged initial phase；According to positive phase again, again antiphase, average again
Initial phase obtains cosine function modulation sequence, SIN function modulation sequence；Sequence is modulated according to default fine tuning frequency, cosine function
Row, SIN function modulation sequence obtain reduce the frequency sine sequence of function, by reduce signal sequence fundamental frequency, it can be achieved that
The number of cycles of fundamental signal is blocked, with the amplitude maximum of fundamental wave component in electric power signal, solves the integer multiples of fundamental signal
Phase truncated problem improves the accuracy of sine parameter calculating, is suitble to practical application.
Description of the drawings
Fig. 1 is the method flow diagram for the reduction frequency sine sequence of function that electric power signal is obtained in one embodiment；
Fig. 2 is first positive sequence and first anti-pleat sequence diagram in one embodiment；
Fig. 3 is the system structure diagram for the reduction frequency sine sequence of function that electric power signal is obtained in one embodiment.
Specific implementation mode
To make the objectives, technical solutions, and advantages of the present invention clearer, below in conjunction with attached drawing to the present invention make into
One step it is described in detail.
Although the step in the present invention is arranged with label, it is not used to limit the precedence of step, unless
Based on the execution of the order or certain step that specify step needs other steps, otherwise the relative rank of step is
It is adjustable.
In one embodiment obtain electric power signal reduction the frequency sine sequence of function method, as shown in Figure 1, include with
Lower step：
Step S101：According to the lower limit of frequency power signal range, preset sample frequency and default integer signal period number,
Obtain the preliminary sequence length of the electric power signal；
Step S102：The electric power signal is sampled according to the preliminary sequence length, obtains the electric power signal
Preliminary sequence；
Step S103：To the preliminary sequence into line frequency preliminary survey, the preliminary frequency of the electric power signal is generated, and according to
The preliminary frequency sets the reference frequency of the electric power signal；
Step S104：According to the preset sample frequency and the reference frequency, the unit week of the electric power signal is obtained
Phase sequence length；
Step S105：According to the default integer signal period number and the unit period sequence length, the electricity is obtained
The predetermined sequence length of force signal；
Step S106：According to default starting point and the predetermined sequence length, the electricity is obtained from the preliminary sequence
The first positive sequence of force signal；
Step S107：The first anti-pleat sequence of the electric power signal is obtained according to the first positive sequence；
Step S108：The first positive phase of the electric power signal is obtained according to the first positive sequence, and according to described
First anti-pleat sequence obtains the first antiphase of the electric power signal；
Step S109：The first average of the electric power signal is obtained according to the first positive phase and the first antiphase
Initial phase；
Step S110：The averagely initial phase for the first time is compared with ± π/4, is obtained first compared with ± π/4
Secondary phase fiducial value, and according to the first phase fiducial value and the default starting point, obtain new starting point；
Step S111：According to the new starting point and the predetermined sequence length, obtained from the preliminary sequence described
The sequence positive again of electric power signal, and obtain according to the sequence positive again the pleat sequence anti-again of the electric power signal；
Step S112：The positive phase again of the electric power signal is obtained according to the sequence positive again, and according to described
Anti- pleat sequence obtains the antiphase again of the electric power signal again；
Step S113：The average again of the electric power signal is obtained according to the positive phase again and the antiphase again
Initial phase；
Step S114：Positive sequence is added with the pleat sequence anti-again again by described in, and according to result after being added
With the initial phase that is averaged again, the cosine function modulation sequence of the electric power signal is obtained；
Step S115：The sequence positive again and the pleat sequence anti-again are subtracted each other, and according to the result after subtracting each other
With the initial phase that is averaged again, the SIN function modulation sequence of the electric power signal is obtained；
Step S116：The discrete sine function of default fine tuning frequency is multiplied with the cosine function modulation sequence and obtains institute
The discrete cosine function of the default fine tuning frequency is modulated sequence by the first multiplication sequence for stating electric power signal with the SIN function
Row, which are multiplied, obtains the second multiplication sequence of the electric power signal；
Step S117：Second multiplication sequence and first multiplication sequence are subtracted each other, the electric power signal is obtained
Reduce the frequency sine sequence of function.
It is evidenced from the above discussion that the present invention, which obtains, reduces the frequency sine sequence of function, the accurate of sine parameter calculating is improved
Degree, actual application value are high.
Wherein, for step S101, believed according to the lower limit of frequency power signal range, preset sample frequency and default integer
Number periodicity, obtains the preliminary sequence length of the electric power signal；
The electric power signal is the sinusoidal signal based on a kind of fundamental wave component.Sinusoidal signal extensively criticizes string function signal and remaining
String function signal.
In one embodiment, power system frequency range takes frequency power signal lower limit f in 45Hz-55HzminFor
45Hz；The default integer signal period number C is set according to actual needs2π, in one embodiment, take C2πIt is 13.
In one embodiment, it is formula (1) to obtain the preliminary sequence length：
In formula, NstartFor preliminary sequence length；(int) rounding is indicated；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 electric power signal is sampled according to the preliminary sequence length, obtains the electric power
The preliminary sequence of signal；
In one embodiment, the electric power signal is the cosine function signal of single fundamental frequency, obtains the electric power letter
Number preliminary sequence be formula (2)：
Wherein, Xstart(n) it is the preliminary sequence；A is signal amplitude, unit v；ωiFor signal frequency, T is between sampling
Every the time, f is preset sample frequency, and unit Hz, n are series of discrete number,For preliminary sequence initial phase, NstartFor preliminary sequence
Length.
The preliminary frequency of the electric power signal is generated to the preliminary sequence into line frequency preliminary survey for step S103, and
The reference frequency of the electric power signal is set according to the preliminary frequency；
Can by zero friendship method, the algorithm based on filtering, the algorithm based on Wavelet Transformation Algorithm, based on neural network, be based on
The frequency algorithm of DFT transform or based on the frequency algorithm of phase difference to the preliminary sequence into line frequency preliminary survey, obtain it is described just
Synchronizing frequency.
In one embodiment, it is formula (3) to generate the preliminary frequency：
ωo(3)；
Wherein, ωoFor preliminary frequency；
Preferably, it is formula (4) that the reference frequency, which is equal to the preliminary frequency,：
ωs=ωo(4)；
Wherein, ωsFor reference frequency, ωoFor preliminary frequency.
The list of the electric power signal is obtained according to the preset sample frequency and the reference frequency for step S104
Bit period sequence length；
In one embodiment, the unit period sequence length for obtaining the electric power signal is formula (5)：
In formula, N2πFor the unit period 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.
Institute is obtained according to the default integer signal period number and the unit period sequence length for step S105
State the predetermined sequence length of electric power signal；
In one embodiment, it is formula (6) to obtain the predetermined sequence length：
N=(int) [(C2π-1)N2π] (6)；
Wherein, N is predetermined sequence length, and (int) is round numbers, N2πFor the unit period sequence length, C2πIt is default
Integer signal period number.
Institute is obtained from the preliminary sequence according to default starting point and the predetermined sequence length for step S106
State the first positive sequence of electric power signal；
In one embodiment, 0.5 times that starting point is the unit period sequence length is preset；
In one embodiment, it is formula (7) to obtain the first positive sequence：
Wherein, Xstart(n) it is preliminary sequence, X+start(n) it is first positive sequence, PstartDefault starting point, N2πFor institute
Unit period sequence length is stated, (int) is round numbers, and A is signal amplitude, unit v, ωiFor signal frequency, T is the sampling interval
Time, n are series of discrete number,For first positive sequence initial phase, N predetermined sequence length.
First forward direction sequence pattern expression, shown in Fig. 2.
For step S107, the first anti-pleat sequence of the electric power signal is obtained according to the first positive sequence；
In one embodiment, it is formula (8) to obtain first anti-pleat sequence：
X-start(- n)=X+start(N-n)=Acos (- ωi Tn+β1)
(8)；
N=0,1,2 ..., N-1
Wherein, X-start(- n) is first anti-pleat sequence, X+start(n) it is 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 first anti-pleat sequence initial phase, and N presets sequence
Row length.
The avatars of the first anti-pleat sequence, as shown in Figure 2.
For step S108, the first positive phase of the electric power signal is obtained according to the first positive sequence, and according to
The first anti-pleat sequence obtains the first antiphase of the electric power signal；
In one embodiment, first positive phase and first antiphase be based on it is orthogonal mixing and integral and calculating as a result,
In the mixing interfering frequency for not considering orthogonal mixing, orthogonal mixing is expressed as formula (9), and integral and calculating is expressed as formula (10)：
Wherein, R+start(n) it is first positive real frequency mixing sequence, I+start(n) it is that first weakened body resistance frequency is mixed sequence, R-start
(- n) is first anti-real frequency mixing sequence, I-start(- n) is first anti-empty frequency mixing sequence, cos (ωs) or cos (- ω Tns
Tn it is) 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, ωiFor signal frequency, T is sampling interval duration, and n is series of discrete
Number,For first positive sequence initial phase, β 1 is first positive sequence initial phase, and N is predetermined sequence length.
In formula, R+startFirst positive real frequency integrated value, unit dimensionless, I+startFor first weakened body resistance frequency integrated value, unit without
Dimension, R-startFor first anti-real frequency integrated value, unit dimensionless, I-startIt is immeasurable for first anti-empty frequency mixing integrated value, unit
Guiding principle, Ω 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
Row length,For first positive sequence initial phase, β 1 is first positive sequence initial phase.
In one embodiment, the expression formula for obtaining first positive phase and first antiphase is (11)：
In formula, PH+startFor first positive phase, PH-startFor first antiphase, R+startFirst positive real frequency integrated value, unit
Dimensionless, I+startFor first weakened body resistance frequency integrated value, unit dimensionless, R-startFor first anti-real frequency integrated value, unit dimensionless,
I-startFor first anti-empty frequency mixing integrated value, unit dimensionless, Ω is signal frequency ωiWith reference frequency ωsFrequency difference, T is
Sampling interval duration, N are predetermined sequence length,For first positive sequence initial phase, β 1 is first anti-pleat sequence initial phase.
For step S109, the first of the electric power signal is obtained according to the first positive phase and the first antiphase
Average initial phase；
In one embodiment, the expression formula for obtaining first averagely initial phase is (12)：
In formula, PHstart-avgFor the initial phase that is averaged for the first time, PH+startFor first positive phase, PH-startFor first antiphase,For first positive sequence initial phase, β 1 is first anti-pleat sequence initial phase.
For step S110, the averagely initial phase for the first time is compared with ± π/4, is obtained compared with ± π/4
First phase fiducial value, and according to the first phase fiducial value and the default starting point, obtain new starting point；
In one embodiment, by initial phase and the PH of being averaged for the first timestart-avgIt is compared as formula (13) with ± π/4：
In formula, Δ PHcomFor first phase fiducial value, unit rad, PHstart-avgFor the initial phase that is averaged for the first time.
In one embodiment, it is formula (14) to obtain the new starting point：
In formula, PnewFor new starting point, unit dimensionless, PstartTo preset starting point, Δ PHcomFor first phase bit comparison
Value, unit rad, N2πFor unit periodic sequence length, (int) is round numbers.
Step S111 is obtained according to the new starting point and the predetermined sequence length from the preliminary sequence
The sequence positive again of the electric power signal, and obtain according to the sequence positive again the pleat sequence anti-again of the electric power signal
Row；
In one embodiment, positive sequence and anti-pleat sequence is formula (15) again again：
In formula, X+end(n) it is positive sequence again, X-end(- n) is anti-pleat sequence again, PnewFor new starting point, unit without
Dimension,For positive sequence initial phase again, β 2 is anti-pleat sequence initial phase again, ωiFor signal frequency, T is the sampling interval
Time, n are series of discrete number, and N is predetermined sequence length.
For step S112, the positive phase again of the electric power signal is obtained according to the sequence positive again, and according to
The pleat sequence anti-again obtains the antiphase again of the electric power signal；
In one embodiment, again positive phase and again antiphase be based on it is orthogonal mixing and digital filtering calculate knot
Fruit.The digital filtering is made of 6 grades of rectangular window arithmetic mean filters of 2 kinds of filtering parameters.
In the mixing interfering frequency for not considering orthogonal mixing, orthogonal mixing is expressed as formula (16), and the 6 of 2 kinds of filtering parameters
Grade rectangular window arithmetic mean filter filtering calculation expression is formula (17)：
In formula, R+end(n) it is that positive real frequency is mixed sequence, I again+end(n) it is that weakened body resistance frequency is mixed sequence, R again-end(-n)
For anti-real frequency mixing sequence again, I-end(- n) is that anti-empty frequency is mixed sequence, cos (ω agains) or cos (- ω TnsTn) it is
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, ωiFor signal frequency, T is sampling interval duration, and n is series of discrete number,
For positive sequence initial phase again, β 2 is that anti-pleat sequence initial phase, N are predetermined sequence length again.
In formula, R+endFor positive real frequency digital filtering final value again, unit dimensionless；I+endFor weakened body resistance frequency digital filtering again
Final value, unit dimensionless；R-endFor anti-digital filtering final value again, unit dimensionless；I-endIt is whole for anti-empty frequency digital filtering again
Value, unit dimensionless；Ω is signal frequency ωiWith reference frequency ωsFrequency difference；K (Ω) is amplitude of the digital filtering in frequency difference Ω
Gain, unit dimensionless；T is sampling interval duration；For positive sequence initial phase again；β 2 is anti-pleat sequence first phase again
Position；ND1For filtering parameter 1, i.e., to ND1A continuous centrifugal pump is added, and then takes its arithmetic mean of instantaneous value defeated as this filter value
Go out；ND2For filtering parameter 2, i.e., to ND2A continuous centrifugal pump is added, and then takes its arithmetic mean of instantaneous value defeated as this filter value
Go out；NDSequence length is used for digital filtering, is the summation of 6 grades of rectangular window arithmetic mean filter filtering parameters in quantity, it is small
In 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,
The mixing interfering frequency that 1/3 subharmonic of purpose pair generates carries out depth inhibition；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 ND2Expression formula is formula (18)：
ND1=(int) (1.5N2π)
(18)；
ND2=2N2π
In formula, 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, again positive phase and again antiphase expression formula be (19)：
In formula, PH+endFor positive phase again, PH-endFor antiphase again, R+endFor positive real frequency integrated value again, unit without
Dimension, I+endFor weakened body resistance frequency integrated value again, unit dimensionless, R-endFor anti-real frequency integrated value again, unit dimensionless, I-end
For anti-empty frequency mixing integrated value again, unit dimensionless, Ω is signal frequency ωiWith reference frequency ωsFrequency difference, T be sampling between
Every the time, NDSequence length is used for digital filtering,For positive sequence initial phase again, β 2 is anti-pleat sequence first phase again
Position.
For step S113, the electric power signal is obtained again according to the positive phase again and the antiphase again
Average initial phase；
In one embodiment, the expression formula for obtaining the initial phase that is averaged again is (20)：
In formula, PHend-avgFor the initial phase that is averaged again, PH+endFor positive phase again, PH-endFor antiphase again,For
Positive sequence initial phase again, β 2 are anti-pleat sequence initial phase again.
For step S114, positive sequence is added with the pleat sequence anti-again again by described in, and according to after being added
As a result with the initial phase that is averaged again, the cosine function modulation sequence of the electric power signal is obtained；
In one embodiment, it is (21) to obtain cosine function modulation sequence expression formula：
In formula, Xcos(n) it is cosine function modulation sequence；A is cosine function modulation sequence amplitude, unit v；For cosine function modulation sequence initial phase, ωiFor signal frequency, T is sampling interval duration, and n is series of discrete
Number, N are predetermined sequence length,For positive sequence initial phase again, β 2 is anti-pleat sequence initial phase again.
For step S115, the sequence positive again and the pleat sequence anti-again are subtracted each other, and according to subtracting each other after
As a result with the initial phase that is averaged again, the SIN function modulation sequence of the electric power signal is obtained；
In one embodiment, it is (22) to obtain SIN function modulation sequence expression formula：
In formula, Xsin(n) be SIN function modulation sequence, A be SIN function modulation sequence amplitude, unit v,For cosine function modulation sequence initial phase, ωiFor signal frequency, T is sampling interval duration, and n is series of discrete
Number, N are predetermined sequence length,For positive sequence initial phase again, β 2 is anti-pleat sequence initial phase again.
For step S116, the discrete sine function of default fine tuning frequency is multiplied with the cosine function modulation sequence and is obtained
The first multiplication sequence for obtaining the electric power signal, by the discrete cosine function of the default fine tuning frequency and the SIN function tune
Sequence processed, which is multiplied, obtains the second multiplication sequence of the electric power signal；
In one embodiment, the default fine tuning frequency is the positive real number less than or equal to actual signal frequency 1%, unit
Rad/s is expressed as formula (23)：
Ωset
(23)
Ωset≤0.01ωi
In formula, ΩsetTo finely tune frequency, unit rad/s, Ωset≤0.01ωi。
The discrete sine function of the fine tuning frequency is multiplied with the cosine function modulation sequence and obtains the first multiplication sequence
It is classified as formula (24)：
The discrete cosine function of the fine tuning frequency is multiplied with the SIN function modulation sequence and obtains the second multiplication sequence
It is classified as formula (25)：
In formula, X1 (n) is the first multiplication sequence, and X2 (n) is the second multiplication sequence, sin (ΩsetTn) it is fine tuning frequency
The discrete sine function of rate, cos (ΩsetTn) it is the discrete cosine function for finely tuning frequency.
For step S117, second multiplication sequence and first multiplication sequence are subtracted each other, obtain the electric power letter
Number the reduction frequency sine sequence of function；
In one embodiment, it is formula (26) that obtaining, which reduces the frequency sine sequence of function,：
In formula, Xsin+f(n) it is to reduce the frequency sine sequence of function, sequence frequency reduces Ωset。
The system that the reduction frequency sine sequence of function of electric power signal is obtained in one embodiment, as shown in figure 3, including：
Preliminary sequence length modules 301 are used for according to the lower limit of frequency power signal range, preset sample frequency and preset
Integer signal period number obtains the preliminary sequence length of the electric power signal；
Preliminary sequence module 302 obtains institute for being sampled to the electric power signal according to the preliminary sequence length
State the preliminary sequence of electric power signal；
Frequency preliminary survey module 303, for, into line frequency preliminary survey, generating the preliminary of the electric power signal to the preliminary sequence
Frequency, and set according to the preliminary frequency reference frequency of the electric power signal；
Unit period sequence length module 304, for according to the preset sample frequency and the reference frequency, obtaining institute
State the unit period sequence length of electric power signal；
Predetermined sequence length modules 305, for according to the default integer signal period number and the unit period sequence
Length obtains the predetermined sequence length of the electric power signal；
First forward direction block 306, presets starting point and the predetermined sequence length, from the preliminary sequence for basis
The first positive sequence of the electric power signal is obtained in row；
First anti-pleat block 307, for obtaining the first anti-of the electric power signal according to the first positive sequence
Pleat sequence；
First phase module 308, the first positive phase for obtaining the electric power signal according to the first positive sequence,
And the first antiphase of the electric power signal is obtained according to the first anti-pleat sequence；
Be averaged initial phase module 309 for the first time, described in being obtained according to the first positive phase and the first antiphase
The initial phase that is averaged for the first time of electric power signal；
Phase bits comparison module 310, for will it is described for the first time averagely initial phase be compared with ± π/4, acquisition and it is described ±
The first phase fiducial value that π/4 are compared, and according to the first phase fiducial value and the default starting point, obtain new starting
Point；
Block 311 again are used for according to the new starting point and the predetermined sequence length, from the preliminary sequence
The middle sequence positive again for obtaining the electric power signal, and the electric power signal is obtained again according to the sequence positive again
Anti- pleat sequence；
Phase module 312 again, the positive phase again for obtaining the electric power signal according to the sequence of forward direction again,
And the antiphase again of the electric power signal is obtained according to the pleat sequence anti-again；
Be averaged initial phase module 313 again, described in positive phase and the antiphase again obtain again according to
The initial phase that is averaged again of electric power signal；
Cosine function modulation sequence module 314, for the sequence of forward direction again to be added with the pleat sequence anti-again,
And according to result after being added and the initial phase that is averaged again, obtain the cosine function modulation sequence of the electric power signal；
SIN function modulation sequence module 315, for the sequence positive again and the pleat sequence anti-again to be subtracted each other,
And according to after subtracting each other result and it is described again be averaged initial phase, obtain the SIN function modulation sequence of the electric power signal；
Multiplication sequence module 316, for the discrete sine function of default fine tuning frequency to be modulated sequence with the cosine function
Row, which are multiplied, obtains the first multiplication sequence of the electric power signal, by the discrete cosine function of the default fine tuning frequency with it is described just
String FUNCTION MODULATION sequence, which is multiplied, obtains the second multiplication sequence of the electric power signal；
Frequency sine sequence of function module 317 is reduced, is used for second multiplication sequence and first multiplication sequence
Subtract each other, obtains the reduction frequency sine sequence of function of the electric power signal.
It is evidenced from the above discussion that the present invention, which obtains, reduces the frequency sine sequence of function, the accurate of sine parameter calculating is improved
Degree meets actual needs.
Wherein, preliminary sequence length modules 301 according to the lower limit of frequency power signal range, preset sample frequency and are preset
Integer signal period number obtains the preliminary sequence length of the electric power signal；
The electric power signal is the sinusoidal signal based on a kind of fundamental wave component.Sinusoidal signal extensively criticizes string function signal and remaining
String function signal.
In one embodiment, power system frequency range takes frequency power signal lower limit f in 45Hz-55HzminFor
45Hz；The default integer signal period number C is set according to actual needs2π, in one embodiment, take C2πIt is 13.
In one embodiment, it is formula (1) to obtain the preliminary sequence length：
In formula, NstartFor preliminary sequence length；(int) rounding is indicated；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.
Preliminary sequence module 302 samples the electric power signal according to the preliminary sequence length, obtains the electricity
The preliminary sequence of force signal；In one embodiment, the electric power signal is the cosine function signal of single fundamental frequency, obtains institute
The preliminary sequence for stating electric power signal is formula (2)：
Wherein, Xstart(n) it is the preliminary sequence；A is signal amplitude, unit v；ωiFor signal frequency, T is between sampling
Every the time, f is preset sample frequency, and unit Hz, n are series of discrete number,For preliminary sequence initial phase, NstartFor preliminary sequence
Length.
Frequency preliminary survey module 303 can hand over method, the algorithm based on filtering by zero, be based on Wavelet Transformation Algorithm, based on nerve
The algorithm of network, the frequency algorithm based on DFT transform or based on the frequency algorithm of phase difference to the preliminary sequence into line frequency
Preliminary survey obtains the preliminary frequency.
In one embodiment, it is formula (3) to generate the preliminary frequency：
ωo(3)；
Preferably, it is formula (4) that the reference frequency, which is equal to the preliminary frequency,：
ωs=ωo(4)；
Wherein, ωsFor reference frequency, ωoFor preliminary frequency.
Unit period sequence length module 304 obtains the electricity according to the preset sample frequency and the reference frequency
The unit period sequence length of force signal；
In one embodiment, the unit period sequence length for obtaining the electric power signal is formula (5)：
In formula, N2πFor the unit period 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 unit period sequence length integer.
Predetermined sequence length modules 305 according to the default integer signal period number and the unit period sequence length,
Obtain the predetermined sequence length of the electric power signal；
In one embodiment, it is formula (6) to obtain the predetermined sequence length：
N=(int) [(C2π-1)N2π] (6)；
Wherein, N is predetermined sequence length, and (int) is round numbers, N2πFor the unit period sequence length, C2πIt is default
Integer signal period number.
First forward direction block 306 is according to default starting point and the predetermined sequence length, from the preliminary sequence
Obtain the first positive sequence of the electric power signal；
In one embodiment, 0.5 times that starting point is the unit period sequence length is preset；
In one embodiment, it is formula (7) to obtain the first positive sequence：
Wherein, Xstart(n) it is preliminary sequence, X+start(n) it is first positive sequence, PstartDefault starting point, N2πFor institute
Unit period sequence length is stated, (int) is round numbers, and A is signal amplitude, unit v, ωiFor signal frequency, T is the sampling interval
Time, n are series of discrete number,For first positive sequence initial phase, N predetermined sequence length.First forward direction sequence pattern table
It reaches, shown in Fig. 2.
First anti-pleat block 307 obtains the first anti-pleat sequence of the electric power signal according to the first positive sequence
Row；
In one embodiment, it is formula (8) to obtain first anti-pleat sequence：
X-start(- n)=X+start(N-n)=Acos (- ωi Tn+β1)
(8)；
N=0,1,2 ..., N-1
Wherein, X-start(- n) is first anti-pleat sequence, X+start(n) it is 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 first anti-pleat sequence initial phase, and N presets sequence
Row length.The avatars of the first anti-pleat sequence, as shown in Figure 2.
First phase module 308 obtains the first positive phase of the electric power signal, and root according to the first positive sequence
The first antiphase of the electric power signal is obtained according to the first anti-pleat sequence；
In one embodiment, first positive phase and first antiphase be based on it is orthogonal mixing and integral and calculating as a result,
In the mixing interfering frequency for not considering orthogonal mixing, orthogonal mixing is expressed as formula (9), and integral and calculating is expressed as formula (10)：
Wherein, R+start(n) it is first positive real frequency mixing sequence, I+start(n) it is that first weakened body resistance frequency is mixed sequence, R-start
(- n) is first anti-real frequency mixing sequence, I-start(- n) is first anti-empty frequency mixing sequence, cos (ωs) or cos (- ω Tns
Tn it is) 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, ωiFor signal frequency, T is sampling interval duration, and n is series of discrete
Number,For first positive sequence initial phase, β 1 is first positive sequence initial phase, and N is predetermined sequence length.
In formula, R+startFirst positive real frequency integrated value, unit dimensionless, I+startFor first weakened body resistance frequency integrated value, unit without
Dimension, R-startFor first anti-real frequency integrated value, unit dimensionless, I-startIt is immeasurable for first anti-empty frequency mixing integrated value, unit
Guiding principle, Ω 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
Row length,For first positive sequence initial phase, β 1 is first positive sequence initial phase.
In one embodiment, the expression formula for obtaining first positive phase and first antiphase is (11)：
In formula, PH+startFor first positive phase, PH-startFor first antiphase, R+startFirst positive real frequency integrated value, unit
Dimensionless, I+startFor first weakened body resistance frequency integrated value, unit dimensionless, R-startFor first anti-real frequency integrated value, unit dimensionless,
I-startFor first anti-empty frequency mixing integrated value, unit dimensionless, Ω is signal frequency ωiWith reference frequency ωsFrequency difference, T is
Sampling interval duration, N are predetermined sequence length,For first positive sequence initial phase, β 1 is first anti-pleat sequence initial phase.
In one embodiment, it is formula (12) that averagely initial phase module 309, which obtains for the first time averagely initial phase, for the first time：
In formula, PHstart-avgFor the initial phase that is averaged for the first time, PH+startFor first positive phase, PH-startFor first antiphase,For first positive sequence initial phase, β 1 is first anti-pleat sequence initial phase.
The averagely initial phase for the first time is compared by phase bits comparison module 310 with ± π/4, is obtained and is compared with ± the π/4
Compared with first phase fiducial value, and according to the first phase fiducial value and the default starting point, obtain new starting point；
In one embodiment, by initial phase and the PH of being averaged for the first timestart-avgIt is compared as formula (13) with ± π/4：
In formula, Δ PHcomFor first phase fiducial value, unit rad, PHstart-avgFor the initial phase that is averaged for the first time.
In one embodiment, it is formula (14) to obtain the new starting point：
In formula, PnewFor new starting point, unit dimensionless, PstartTo preset starting point, Δ PHcomFor first phase bit comparison
Value, unit rad, N2πFor unit periodic sequence length, (int) is round numbers.
Block 311 is obtained according to the new starting point and the predetermined sequence length from the preliminary sequence again
The sequence positive again of the electric power signal is obtained, and obtains the pleat anti-again of the electric power signal according to the sequence positive again
Sequence；
In one embodiment, positive sequence and anti-pleat sequence is formula (15) again again：
In formula, X+end(n) it is positive sequence again, X-end(- n) is anti-pleat sequence again, PnewFor new starting point, unit without
Dimension,For positive sequence initial phase again, β 2 is anti-pleat sequence initial phase again, ωiFor signal frequency, T is the sampling interval
Time, n are series of discrete number, and N is predetermined sequence length.
Phase module 312 obtains the positive phase again of the electric power signal, and root according to the sequence positive again again
The antiphase again of the electric power signal is obtained according to the pleat sequence anti-again；
In one embodiment, again positive phase and again antiphase be based on it is orthogonal mixing and digital filtering calculate knot
Fruit.The digital filtering is made of 6 grades of rectangular window arithmetic mean filters of 2 kinds of filtering parameters.
In the mixing interfering frequency for not considering orthogonal mixing, orthogonal mixing is expressed as formula (16), and the 6 of 2 kinds of filtering parameters
Grade rectangular window arithmetic mean filter filtering calculation expression is formula (17)：
In formula, R+end(n) it is that positive real frequency is mixed sequence, I again+end(n) it is that weakened body resistance frequency is mixed sequence, R again-end(-n)
For anti-real frequency mixing sequence again, I-end(- n) is that anti-empty frequency is mixed sequence, cos (ω agains) or cos (- ω TnsTn) it is
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, ωiFor signal frequency, T is sampling interval duration, and n is series of discrete number,
For positive sequence initial phase again, β 2 is that anti-pleat sequence initial phase, N are predetermined sequence length again.
In formula, R+endFor positive real frequency digital filtering final value again, unit dimensionless；I+endFor weakened body resistance frequency digital filtering again
Final value, unit dimensionless；R-endFor anti-digital filtering final value again, unit dimensionless；I-endIt is whole for anti-empty frequency digital filtering again
Value, unit dimensionless；Ω is signal frequency ωiWith reference frequency ωsFrequency difference；K (Ω) is amplitude of the digital filtering in frequency difference Ω
Gain, unit dimensionless；T is sampling interval duration；For positive sequence initial phase again；β 2 is anti-pleat sequence first phase again
Position；ND1For filtering parameter 1, i.e., to ND1A continuous centrifugal pump is added, and then takes its arithmetic mean of instantaneous value defeated as this filter value
Go out；ND2For filtering parameter 2, i.e., to ND2A continuous centrifugal pump is added, and then takes its arithmetic mean of instantaneous value defeated as this filter value
Go out；NDSequence length is used for digital filtering, is the summation of 6 grades of rectangular window arithmetic mean filter filtering parameters in quantity, it is small
In 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,
The mixing interfering frequency that 1/3 subharmonic of purpose pair generates carries out depth inhibition；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 ND2Expression formula is formula (18)：
ND1=(int) (1.5N2π)
(18)；
ND2=2N2π
In formula, 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, again positive phase and again antiphase expression formula be (19)：
In formula, PH+endFor positive phase again, PH-endFor antiphase again, R+endFor positive real frequency integrated value again, unit without
Dimension, I+endFor weakened body resistance frequency integrated value again, unit dimensionless, R-endFor anti-real frequency integrated value again, unit dimensionless, I-end
For anti-empty frequency mixing integrated value again, unit dimensionless, Ω is signal frequency ωiWith reference frequency ωsFrequency difference, T be sampling between
Every the time, NDSequence length is used for digital filtering,For positive sequence initial phase again, β 2 is anti-pleat sequence first phase again
Position.
The initial phase module 313 that is averaged again obtains the electric power according to the positive phase again and the antiphase again and believes
Number again be averaged initial phase；
In one embodiment, the expression formula for obtaining the initial phase that is averaged again is (20)：
In formula, PHend-avgFor the initial phase that is averaged again, PH+endFor positive phase again, PH-endFor antiphase again,For
Positive sequence initial phase again, β 2 are anti-pleat sequence initial phase again.
Positive sequence is added cosine function modulation sequence module 314 with the pleat sequence anti-again again by described in, and root
According to result after being added and the initial phase that is averaged again, the cosine function modulation sequence of the electric power signal is obtained；At one
In embodiment, it is (21) to obtain cosine function modulation sequence expression formula：
In formula, Xcos(n) it is cosine function modulation sequence；A is cosine function modulation sequence amplitude, unit v；For cosine function modulation sequence initial phase, ωiFor signal frequency, T is sampling interval duration, and n is series of discrete
Number, N are predetermined sequence length,For positive sequence initial phase again, β 2 is anti-pleat sequence initial phase again.
SIN function modulation sequence module 315 subtracts each other the sequence positive again and the pleat sequence anti-again, and root
According to after subtracting each other result and it is described again be averaged initial phase, obtain the SIN function modulation sequence of the electric power signal；
In one embodiment, it is (22) to obtain SIN function modulation sequence expression formula：
In formula, Xsin(n) be SIN function modulation sequence, A be SIN function modulation sequence amplitude, unit v,For cosine function modulation sequence initial phase, ωiFor signal frequency, T is sampling interval duration, and n is series of discrete
Number, N are predetermined sequence length,For positive sequence initial phase again, β 2 is anti-pleat sequence initial phase again.
The discrete sine function of default fine tuning frequency is multiplied by multiplication sequence module 316 with the cosine function modulation sequence
The first multiplication sequence for obtaining the electric power signal, by the discrete cosine function of the default fine tuning frequency and the SIN function
Modulation sequence, which is multiplied, obtains the second multiplication sequence of the electric power signal；
In one embodiment, the default fine tuning frequency is the positive real number less than or equal to actual signal frequency 1%, unit
Rad/s is expressed as formula (23)：
Ωset
(23)
Ωset≤0.01ωi
In formula, ΩsetTo finely tune frequency, unit rad/s, Ωset≤0.01ωi。
The discrete sine function of the fine tuning frequency is multiplied with the cosine function modulation sequence and obtains the first multiplication sequence
It is classified as formula (24)：
The discrete cosine function of the fine tuning frequency is multiplied with the SIN function modulation sequence and obtains the second multiplication sequence
It is classified as formula (25)：
In formula, X1 (n) is the first multiplication sequence, and X2 (n) is the second multiplication sequence, sin (ΩsetTn) it is fine tuning frequency
The discrete sine function of rate, cos (ΩsetTn) it is the discrete cosine function for finely tuning frequency.
Frequency sine sequence of function module 317 is reduced to subtract each other second multiplication sequence and first multiplication sequence,
Obtain the reduction frequency sine sequence of function of the electric power signal；
In one embodiment, it is formula (26) that obtaining, which reduces the frequency sine sequence of function,：
In formula, Xsin+f(n) it is to reduce the frequency sine sequence of function, sequence frequency reduces Ωset。
Each technical characteristic of embodiment described above can be combined arbitrarily, to keep 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, it is all considered to be the range of this specification record.
Several embodiments of the invention above described embodiment only expresses, the description thereof 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 coming 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
Range.Therefore, the protection domain of patent of the present invention should be determined by the appended claims.
Claims (6)
1. a kind of method for the reduction frequency sine sequence of function obtaining electric power signal, which is characterized in that include the following steps：
According to the lower limit of frequency power signal range, preset sample frequency and default integer signal period number, the electric power is obtained
The preliminary sequence length of signal；
The 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 generated, and according to the preliminary frequency
Set the reference frequency of the electric power signal；
According to the preset sample frequency and the reference frequency, the unit period sequence length of the electric power signal is obtained；
According to the default integer signal period number and the unit period sequence length, the default sequence of the electric power signal is obtained
Row length；
According to default starting point and the predetermined sequence length, obtained from the preliminary sequence electric power signal it is first just
To sequence；
The first anti-pleat sequence of the electric power signal is obtained according to the first positive sequence；
The first positive phase of the electric power signal is obtained according to the first positive sequence, and is obtained according to the first anti-pleat sequence
Obtain the first antiphase of the electric power signal；
The initial phase that is averaged for the first time of the electric power signal is obtained according to the first positive phase and the first antiphase；
The averagely initial phase for the first time is compared with ± π/4, obtains the first phase fiducial value compared with ± π/4,
And according to the first phase fiducial value and the default starting point, obtain new starting point；
According to the new starting point and the predetermined sequence length, the electric power signal is obtained from the preliminary sequence again
Positive sequence, and obtain according to the sequence positive again the pleat sequence anti-again of the electric power signal；
The positive phase again of the electric power signal is obtained according to the sequence positive again, and is obtained according to the pleat sequence anti-again
Obtain the antiphase again of the electric power signal；
The initial phase that is averaged again of the electric power signal is obtained according to the positive phase again and the antiphase again；
The sequence positive again is added with the pleat sequence anti-again, and described is averaged according to result after being added and again
Initial phase, according to expression formulaObtain the cosine function modulation sequence of the electric power signal
Xcos(n), wherein n=0,1,2 ..., N-1, X+end(n) it is the sequence positive again, X-end(- n) is the pleat sequence anti-again
Row, PHend-avgFor the initial phase that is averaged again, n is series of discrete number, and N is the predetermined sequence length；
The sequence positive again and the pleat sequence anti-again are subtracted each other, and according to result after subtracting each other and described average again
Initial phase, according to expression formulaObtain the SIN function modulation sequence of the electric power signal
Xsin(n), wherein n=0,1,2 ..., N-1, X+end(n) it is the sequence positive again, X-end(- n) is the pleat sequence anti-again
Row, PHend-avgFor the initial phase that is averaged again, n is series of discrete number, and N is the predetermined sequence length；
The discrete sine function of default fine tuning frequency is multiplied with the cosine function modulation sequence and obtains the electric power signal
The discrete cosine function of the default fine tuning frequency is multiplied with the SIN function modulation sequence and obtains institute by the first multiplication sequence
State the second multiplication sequence of electric power signal；
Second multiplication sequence and first multiplication sequence are subtracted each other, the reduction frequency sine letter of the electric power signal is obtained
Number Sequence.
2. the method for the reduction frequency sine sequence of function according to claim 1 for obtaining electric power signal, which is characterized in that
The electric power signal is the cosine function signal of single fundamental frequency, according to expression formula
Obtain the preliminary sequence Xstart(n), wherein A is signal amplitude, ωiFor signal frequency
Rate,For the initial phase of the preliminary sequence, T is sampling interval duration, and f is the preset sample frequency, and n is series of discrete number,
NstartFor the preliminary sequence length.
3. the method for the reduction frequency sine sequence of function according to claim 1 for obtaining electric power signal, which is characterized in that
According to expression formula X1 (n)=Xcos(n)sin(ΩsetTn) obtain the first multiplication sequence X1 (n), according to expression formula X2 (n)=
Xsin(n)cos(ΩsetTn the second multiplication sequence X2 (n) the wherein n=0,1,2 ..., N-1, X) are obtainedcos(n) it is described remaining
String FUNCTION MODULATION sequence, Xsin(n) it is the SIN function modulation sequence, sin (ΩsetTn) be it is described it is default fine tuning frequency from
Dissipate SIN function, cos (ΩsetTn) it is the default discrete cosine function for finely tuning frequency, ΩsetFor the default fine tuning frequency
Rate, T are sampling interval duration, and n is series of discrete number, and N is the predetermined sequence length.
4. a kind of system for the reduction frequency sine sequence of function obtaining electric power signal, which is characterized in that including：
Preliminary sequence length modules, for being believed according to the lower limit of frequency power signal range, preset sample frequency and default integer
Number periodicity, obtains the preliminary sequence length of the electric power signal；
Preliminary sequence module obtains the electric power for being sampled to the electric power signal according to the preliminary sequence length
The preliminary sequence of signal；
Frequency preliminary survey module, for, into line frequency preliminary survey, generating the preliminary frequency of the electric power signal to the preliminary sequence, and
The reference frequency of the electric power signal is set according to the preliminary frequency；
Unit period sequence length module, for according to the preset sample frequency and the reference frequency, obtaining the electric power
The unit period sequence length of signal；
Predetermined sequence length modules, for according to the default integer signal period number and the unit period sequence length, obtaining
Obtain the predetermined sequence length of the electric power signal；
First forward direction block, for according to starting point and the predetermined sequence length is preset, being obtained from the preliminary sequence
Obtain the first positive sequence of the electric power signal；
First anti-pleat block, the first anti-pleat sequence for obtaining the electric power signal according to the first positive sequence；
First phase module, the first positive phase for obtaining the electric power signal according to the first positive sequence, and according to
The first anti-pleat sequence obtains the first antiphase of the electric power signal；
Be averaged initial phase module for the first time, for obtaining the electric power signal according to the first positive phase and the first antiphase
For the first time averagely initial phase；
Phase bits comparison module is obtained for the averagely initial phase for the first time to be compared with ± π/4 compared with ± π/4
First phase fiducial value, and according to the first phase fiducial value and the default starting point, obtain new starting point；
Block again, for according to the new starting point and the predetermined sequence length, being obtained from the preliminary sequence
The sequence positive again of the electric power signal, and obtain according to the sequence positive again the pleat sequence anti-again of the electric power signal
Row；
Phase module again, the positive phase again for obtaining the electric power signal according to the sequence of forward direction again, and according to
The pleat sequence anti-again obtains the antiphase again of the electric power signal；
Be averaged initial phase module again, for positive phase and the antiphase again to obtain the electric power signal again according to
Again be averaged initial phase；
Cosine function modulation sequence module, for the sequence of forward direction again to be added with the pleat sequence anti-again, and according to
Result after being added and the initial phase that is averaged again, according to expression formula Obtain institute
State the cosine function modulation sequence X of electric power signalcos(n), wherein n=0,1,2 ..., N-1, X+end(n) it is described positive again
Sequence, X-end(- n) is the pleat sequence anti-again, PHend-avgFor the initial phase that is averaged again, n is series of discrete number, and N is
The predetermined sequence length；
SIN function modulation sequence module, for the sequence positive again and the pleat sequence anti-again to be subtracted each other, and according to
Result after subtracting each other and the initial phase that is averaged again, according to expression formula Obtain institute
State the SIN function modulation sequence X of electric power signalsin(n), wherein n=0,1,2 ..., N-1, X+end(n) it is described positive again
Sequence, X-end(- n) is the pleat sequence anti-again, PHend-avgFor the initial phase that is averaged again, n is series of discrete number, and N is
The predetermined sequence length；
Multiplication sequence module is obtained for the discrete sine function of default fine tuning frequency to be multiplied with the cosine function modulation sequence
The first multiplication sequence for obtaining the electric power signal, by the discrete cosine function of the default fine tuning frequency and the SIN function tune
Sequence processed, which is multiplied, obtains the second multiplication sequence of the electric power signal；
Frequency sine sequence of function module is reduced to obtain for subtracting each other second multiplication sequence and first multiplication sequence
Obtain the reduction frequency sine sequence of function of the electric power signal.
5. the system of the reduction frequency sine sequence of function according to claim 4 for obtaining electric power signal, which is characterized in that
The electric power signal is the cosine function signal of single fundamental frequency, and the preliminary sequence module is according to expression formulaObtain the preliminary sequence Xstart(n), wherein
A is signal amplitude, ωiFor signal frequency,For the initial phase of the preliminary sequence, T is sampling interval duration, and f is described default
Sample frequency, n are series of discrete number, NstartFor the preliminary sequence length.
6. the system of the reduction frequency sine sequence of function according to claim 4 for obtaining electric power signal, which is characterized in that
The multiplication sequence module is according to expression formula X1 (n)=Xcos(n)sin(ΩsetTn the first multiplication sequence X1 (n)) is obtained,
According to expression formula X2 (n)=Xsin(n)cos(ΩsetTn) obtain the second multiplication sequence X2 (n), wherein n=0,1,2 ...,
N-1, Xcos(n) it is the cosine function modulation sequence, Xsin(n) it is the SIN function modulation sequence, sin (ΩsetTn it is) institute
State the discrete sine function of default fine tuning frequency, cos (ΩsetTn) it is the default discrete cosine function for finely tuning frequency, Ωset
For the default fine tuning frequency, T is sampling interval duration, and n is series of discrete number, and N is the predetermined sequence length.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510891736.7A CN105548706B (en) | 2015-12-04 | 2015-12-04 | Obtain the method and system of the reduction frequency sine sequence of function of electric power signal |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510891736.7A CN105548706B (en) | 2015-12-04 | 2015-12-04 | Obtain the method and system of the reduction frequency sine sequence of function of electric power signal |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105548706A CN105548706A (en) | 2016-05-04 |
| CN105548706B true CN105548706B (en) | 2018-09-04 |
Family
ID=55828027
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510891736.7A Active CN105548706B (en) | 2015-12-04 | 2015-12-04 | Obtain the method and system of the reduction frequency sine sequence of function of electric power signal |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105548706B (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 |
| CN104502700A (en) * | 2014-12-29 | 2015-04-08 | 广东电网有限责任公司电力科学研究院 | Sinusoidal parameter measurement method and system of power signal |
| CN104502706A (en) * | 2014-12-29 | 2015-04-08 | 广东电网有限责任公司电力科学研究院 | Harmonic amplitude 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 CN201510891736.7A patent/CN105548706B/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 |
| CN104502700A (en) * | 2014-12-29 | 2015-04-08 | 广东电网有限责任公司电力科学研究院 | Sinusoidal parameter measurement method and system of power signal |
| CN104502706A (en) * | 2014-12-29 | 2015-04-08 | 广东电网有限责任公司电力科学研究院 | Harmonic amplitude 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 |
|---|---|
| CN105548706A (en) | 2016-05-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105067880B (en) | The method and system of orthogonal modulation is carried out to electric power signal | |
| CN105548706B (en) | Obtain the method and system of the reduction frequency sine sequence of function of electric power signal | |
| CN105425036B (en) | Obtain the method and system of the reduction frequency cosine function sequence of electric power signal | |
| CN105548691B (en) | Obtain the method and system of the raising frequency sine sequence of function of electric power signal | |
| CN105548699B (en) | Obtain the method and system of the raising frequency reference SIN function sequence of electric power signal | |
| CN105467212B (en) | Obtain the method and system of the raising frequency cosine function sequence of electric power signal | |
| CN105548700B (en) | Obtain the method and system of the raising frequency orthogonal sequence of electric power signal | |
| CN105425035B (en) | Obtain the method and system for dropping low-frequency zero initial phase reference sinusoidal sequence of function | |
| CN105467211B (en) | Obtain the method and system of the raising frequency reference cosine function sequence of electric power signal | |
| CN105548701B (en) | Obtain the method and system of the reduction frequency orthogonal sequence of electric power signal | |
| CN105548705B (en) | Obtain the method and system of the raising frequency reference orthogonal sequence of electric power signal | |
| CN105548689B (en) | Obtain the method and system of the reduction frequency reference orthogonal sequence of electric power signal | |
| CN105548704B (en) | The method and system of fine setting frequency reference sequence is obtained from electric power signal | |
| CN105425031B (en) | The orthogonal frequency multiplication sequence method of zero initial phase benchmark and system are obtained from electric power signal | |
| CN105548702B (en) | Obtain the method and system for dropping low-frequency zero initial phase benchmark cosine function sequence | |
| CN105548698B (en) | Frequency power signal method for trimming and system | |
| CN105403768B (en) | The method and system of cosine function times frequency sequence is obtained from electric power signal | |
| CN105425032B (en) | The method and system of SIN function times frequency sequence is obtained from electric power signal | |
| CN105372490B (en) | The method and system of frequency measurement based on zero orthogonal times of frequency sequence of initial phase benchmark | |
| CN105548693B (en) | Obtain the method and system of zero initial phase reference sinusoidal function times frequency sequence | |
| CN105388359B (en) | The method and system of orthogonal times of frequency sequence is obtained from electric power signal | |
| CN105548688B (en) | Frequency measurement method and system are carried out according to zero initial phase benchmark cosine function sequence | |
| CN105548692B (en) | The method and system of frequency measurement is carried out according to zero initial phase benchmark orthogonal sequence | |
| CN105301356B (en) | Obtain the method and system of zero initial phase benchmark cosine function times frequency sequence | |
| CN105425040B (en) | Based on zero initial phase benchmark cosine function frequency multiplication sequence frequency measuring method and system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |