CN105548705B - Obtain the method and system of the raising frequency reference orthogonal sequence of electric power signal - Google Patents

Obtain the method and system of the raising frequency reference orthogonal sequence of electric power signal Download PDF

Info

Publication number
CN105548705B
CN105548705B CN201510891667.XA CN201510891667A CN105548705B CN 105548705 B CN105548705 B CN 105548705B CN 201510891667 A CN201510891667 A CN 201510891667A CN 105548705 B CN105548705 B CN 105548705B
Authority
CN
China
Prior art keywords
sequence
electric power
frequency
power signal
initial phase
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
Application number
CN201510891667.XA
Other languages
Chinese (zh)
Other versions
CN105548705A (en
Inventor
朱亚清
李军
潘凤萍
吴乐
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Electric Power Research Institute of Guangdong Power Grid Co Ltd
Original Assignee
Electric Power Research Institute of Guangdong Power Grid Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Electric Power Research Institute of Guangdong Power Grid Co Ltd filed Critical Electric Power Research Institute of Guangdong Power Grid Co Ltd
Priority to CN201510891667.XA priority Critical patent/CN105548705B/en
Publication of CN105548705A publication Critical patent/CN105548705A/en
Application granted granted Critical
Publication of CN105548705B publication Critical patent/CN105548705B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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

Abstract

The invention discloses a kind of method and system for the raising frequency reference orthogonal sequence for obtaining electric power signal, methods described includes：Obtain electric power signal preliminary sequence length and preliminary sequence；Enter line frequency preliminary survey to preliminary sequence, set reference frequency；Electric power signal unit period sequence length and predetermined sequence length are obtained, further obtains first positive/anti-pleat sequence；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, obtain zero initial phase benchmark it is remaining/SIN function modulation sequence, the multiplication sequence of first/second/the 3rd/the 4th is further obtained, is finally obtained according to the multiplication sequence of first/second/the 3rd/the 4th and carries high-frequency benchmark orthogonal sequence.The present invention improves the degree of accuracy that sine parameter calculates, and is adapted to application.

Description

Obtain the method and system of the raising frequency reference orthogonal sequence of electric power signal
Technical field
The present invention relates to technical field of electric power, more particularly to a kind of orthogonal sequence of raising frequency reference for obtaining electric power signal The method and system of row.
Background technology
The frequency measurement of power system, phase measurement, amplitude measurement etc. are the measurement of sine parameter in itself.Quickly Fourier transformation and discrete Fourier transform are to realize the basic skills of sine parameter measurement, widely should be had in power system With.But the above method blocks in the non-integer of signal sampling process, cause spectrum leakage, and spectrum leakage can produce phase The error answered.
In terms of power system sine parameter measurement, there are many measuring methods, repaiied such as adding window mouth functional based method, using interpolation Correction method etc., reduce the influence of spectrum leakage problem.But above-mentioned measuring method measurement accuracy is low, it is being not suitable for high accuracy just The measurement of string parameter.
The content of the invention
Based on the above situation, the present invention proposes a kind of method for the raising frequency reference orthogonal sequence for obtaining electric power signal And system, improve the degree of accuracy that sine parameter measures.
To achieve these goals, the embodiment of technical solution of the present invention is：
A kind of method for the raising frequency reference orthogonal sequence for obtaining electric power signal, comprises the following steps：
According to the lower limit of frequency power signal scope, 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；
Line frequency preliminary survey is entered to the preliminary sequence, generates the preliminary frequency of the electric power signal, and according to described preliminary The reference frequency of electric power signal described in frequency setting；
According to the preset sample frequency and the reference frequency, the unit period sequence for obtaining the electric power signal is grown 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, the predetermined sequence length is odd number；
According to default starting point and the predetermined sequence length, the first of the electric power signal is obtained from the preliminary sequence Secondary positive 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；
By the averagely initial phase for the first time compared with ± π/4, the first phase bit comparison compared with ± π/4 is obtained 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 according to the pleat sequence anti-again of the positive sequence acquisition again 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；
Will described in again positive sequence be added with the pleat sequence anti-again, and according to the result after addition and it is described again Average initial phase, obtain 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 the result after subtracting each other and it is described again Average initial phase, obtain the SIN function modulation sequence of the electric power signal；
Exported from the cosine function modulation sequence central point, obtain zero initial phase benchmark cosine letter of the electric power signal Number modulation sequence；Exported from the SIN function modulation sequence central point, obtaining zero initial phase benchmark of the electric power signal just String FUNCTION MODULATION sequence；
The discrete sine function of default fine setting frequency is multiplied with the zero initial phase benchmark cosine function modulation sequence and obtained The first multiplication sequence of the electric power signal is obtained, by the discrete cosine function of the default fine setting frequency and the zero initial phase base Quasi sinusoids modulation sequence, which is multiplied, obtains the second multiplication sequence of the electric power signal, by the discrete sine of default fine setting frequency Function is multiplied with the zero initial phase reference sinusoidal FUNCTION MODULATION sequence obtains the 3rd multiplication sequence of the electric power signal, by institute The discrete cosine function for stating default fine setting frequency is multiplied the acquisition electricity with the zero initial phase benchmark cosine function modulation sequence 4th multiplication sequence of force signal；
First multiplication sequence is added with second multiplication sequence, it is high-frequency to obtain carrying for the electric power signal The zero initial phase reference sinusoidal sequence of function；4th multiplication sequence is subtracted each other with the 3rd multiplication sequence, obtains the electricity Force signal puies forward high-frequency zero initial phase benchmark cosine function sequence；
Put forward the high-frequency zero initial phase reference sinusoidal sequence of function according to described and described carry high-frequency zero initial phase base The quasi- cosine function sequence acquisition electric power signal carries high-frequency benchmark orthogonal sequence.
A kind of system for the raising frequency reference orthogonal sequence for obtaining electric power signal, including：
Preliminary sequence length modules, for the lower limit according to frequency power signal scope, preset sample frequency and preset whole Number signal period number, obtain the preliminary sequence length of the electric power signal；
Preliminary sequence module, for being sampled according to the preliminary sequence length to the electric power signal, described in acquisition The preliminary sequence of electric power signal；
Frequency preliminary survey module, for entering line frequency preliminary survey to the preliminary sequence, generate the first cadence of the electric power signal Rate, and according to the reference frequency of the preliminary frequency setting 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 according to the default integer signal period number and unit period sequence length Degree, obtains the predetermined sequence length of the electric power signal, the predetermined sequence length is odd number；
First positive block, for the default starting point of basis and the predetermined sequence length, from the preliminary sequence The middle first positive sequence for obtaining the electric power signal；
First anti-pleat block, for obtaining the first anti-pleat sequence of the electric power signal according to the first positive sequence Row；
First phase module, for obtaining the first positive phase of 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, for the averagely initial phase for the first time compared with ± π/4, to be obtained and ± π/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, 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, for the positive phase again according to the positive sequence acquisition again electric power signal, 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 will described in again positive sequence be added with the pleat sequence anti-again, and According to the result after addition 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 the result after subtracting each other and the initial phase that is averaged again, the SIN function modulation sequence of the electric power signal is obtained；
Reference function modulation sequence module, for being exported from the cosine function modulation sequence central point, obtain the electricity Zero initial phase benchmark cosine function modulation sequence of force signal；Exported from the SIN function modulation sequence central point, obtain institute State zero initial phase reference sinusoidal FUNCTION MODULATION sequence of electric power signal；
Multiplication sequence module, for by discrete sine function and the zero initial phase benchmark cosine letter of default fine setting frequency Number modulation sequence, which is multiplied, obtains the first multiplication sequence of the electric power signal, by the discrete cosine function of the default fine setting frequency It is multiplied with the zero initial phase reference sinusoidal FUNCTION MODULATION sequence and obtains the second multiplication sequence of the electric power signal, will be default micro- The discrete sine function of frequency modulation rate is multiplied with the zero initial phase reference sinusoidal FUNCTION MODULATION sequence obtains the electric power signal 3rd multiplication sequence, the discrete cosine function of the default fine setting frequency is modulated into sequence with the zero initial phase benchmark cosine function Row, which are multiplied, obtains the 4th multiplication sequence of the electric power signal；
Frequency reference sequence of function module is improved, first multiplication sequence is added with second multiplication sequence, obtained Obtain the electric power signal puies forward the high-frequency zero initial phase reference sinusoidal sequence of function；By the 4th multiplication sequence and described the Three multiplication sequences subtract each other, and obtain the electric power signal puies forward high-frequency zero initial phase benchmark cosine function sequence；
High-frequency benchmark orthogonal sequence module is put forward, for carrying high-frequency zero initial phase reference sinusoidal function according to Sequence and it is described carry that high-frequency zero initial phase benchmark cosine function sequence obtains the electric power signal carry high-frequency benchmark Orthogonal sequence.
Compared with prior art, beneficial effects of the present invention are：The present invention is obtaining the raising frequency reference of electric power signal just The method and system of sequence is handed over, according to the lower limit of frequency power signal scope, preset sample frequency and preliminary design integer signal period Number, preliminary sequence length is obtained, and electric power signal is tentatively sampled, obtain preliminary sequence；Line frequency is entered to preliminary sequence Preliminary survey, preliminary frequency is generated, set reference frequency；According to preset sample frequency and reference frequency, unit period sequence length is obtained Degree；According to default integer signal period number and unit period sequence length, predetermined sequence length is obtained；Obtained from preliminary sequence First positive sequence, further obtain 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；First phase fiducial value is obtained compared with ± π/4, enters one Step obtains new starting point；According to new starting point and predetermined sequence length, positive sequence, again anti-pleat sequence again are obtained, enters 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；Export, obtain from cosine function modulation sequence central point Zero initial phase benchmark cosine function modulation sequence, exported from SIN function modulation sequence central point, obtaining zero initial phase benchmark just String FUNCTION MODULATION sequence；According to default fine setting frequency, zero initial phase benchmark cosine function modulation sequence, zero initial phase reference sinusoidal FUNCTION MODULATION sequence obtains and carries high-frequency benchmark orthogonal sequence, by improving the fundamental frequency of signal sequence, can be achieved to base The number of cycles of ripple signal is blocked, and with the amplitude maximum of fundamental wave component in electric power signal, the number of cycles for solving fundamental signal is cut Disconnected problem, improves the degree of accuracy of sine parameter calculating, is adapted to practical application.
Brief description of the drawings
Fig. 1 is the method flow diagram for the raising frequency reference SIN function sequence 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 zero initial phase datum mark schematic diagram in one embodiment；
Fig. 4 is the system architecture signal for the raising frequency reference SIN function sequence that electric power signal is obtained in one embodiment Figure.
Embodiment
In order that the object, technical solutions and advantages of the present invention are clearer, the present invention is made below in conjunction with accompanying drawing 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 It specify that the order of step or based on the execution of certain step needs other steps, otherwise the relative rank of step is It is adjustable.
In one embodiment obtain electric power signal raising frequency reference orthogonal sequence method, as shown in figure 1, including with Lower step：
Step S101：According to the lower limit of frequency power signal scope, 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：Line frequency preliminary survey is entered to the preliminary sequence, generates the preliminary frequency of the electric power signal, 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 for obtaining the electric power signal is all 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, the predetermined sequence length are odd number；
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：By the averagely initial phase for the first time compared with ± π/4, obtain 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 according to the pleat sequence anti-again of the positive sequence acquisition again 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 the result after addition 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：Exported from the cosine function modulation sequence central point, obtain zero initial phase of the electric power signal Benchmark cosine function modulation sequence；Exported from the SIN function modulation sequence central point, at the beginning of obtain the electric power signal zero Phase reference SIN function modulation sequence；
Step S117：By the discrete sine function of default fine setting frequency and zero initial phase benchmark cosine function modulation sequence Row, which are multiplied, obtains the first multiplication sequence of the electric power signal, by the discrete cosine function and described zero of the default fine setting frequency Initial phase reference sinusoidal FUNCTION MODULATION sequence, which is multiplied, obtains the second multiplication sequence of the electric power signal, by default fine setting frequency Discrete sine function is multiplied with the zero initial phase reference sinusoidal FUNCTION MODULATION sequence obtains the 3rd multiplication of the electric power signal Sequence, the discrete cosine function of the default fine setting frequency is multiplied with the zero initial phase benchmark cosine function modulation sequence and obtained Obtain the 4th multiplication sequence of the electric power signal；
Step S118：First multiplication sequence is added with second multiplication sequence, obtains the electric power signal Put forward the high-frequency zero initial phase reference sinusoidal sequence of function；4th multiplication sequence is subtracted each other with the 3rd multiplication sequence, Obtain the electric power signal puies forward high-frequency zero initial phase benchmark cosine function sequence；
Step S119：Put forward the high-frequency zero initial phase reference sinusoidal sequence of function according to described and described carry high-frequency zero The initial phase benchmark cosine function sequence acquisition electric power signal carries high-frequency benchmark orthogonal sequence.
It is evidenced from the above discussion that the present invention, which obtains, carries high-frequency benchmark orthogonal sequence, the standard that sine parameter calculates is improved Exactness, actual application value are high.
Wherein, the electric power signal is the sinusoidal signal based on a kind of fundamental wave component.Sinusoidal signal extensively refers to SIN function letter Number and cosine function signal.
In one embodiment, power system frequency scope is in 45Hz-55Hz, power taking force signal lower-frequency limit fminFor 45Hz；According to being actually needed, the default integer signal period number C is set2π, in one embodiment, take C2πFor 13.At one In embodiment, for step S101, it is formula (1) to obtain the preliminary sequence length：
In formula, NstartFor preliminary sequence length；(int) represent to round；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.
In one embodiment, the electric power signal is the cosine function signal of single fundamental frequency, for step S102, is obtained The preliminary sequence for obtaining the 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.
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 at the beginning of entering line frequency to the preliminary sequence based on the frequency algorithm of phase difference Survey, obtain the preliminary frequency.
In one embodiment, it is formula (3) to generate the preliminary frequency：ωo(3)；
Preferably, the reference frequency is formula (4) equal to the preliminary frequency：ωs=ωo(4)；
Wherein, ωsFor reference frequency, ωoFor preliminary frequency.
In one embodiment, for step S104, 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.The error in 1 sampling interval be present in unit period sequence length integer.
In one embodiment, for step S105, it is formula (6) to obtain the predetermined sequence length：
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.
In one embodiment, 0.5 times that starting point is the unit period sequence length is preset；
In one embodiment, for step S106, 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 positive sequence pattern table Reach, shown in Fig. 2.
In one embodiment, for step S107, it is formula (8) to obtain first anti-pleat sequence：
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.
In one embodiment, first positive phase and first antiphase are the results based on orthogonal mixing and integral and calculating, When not considering the mixing interfering frequency of 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) sequence, R are mixed for first weakened body resistance frequency-start (- n) is that first anti-real frequency is mixed sequence, I-start(- n) is that first anti-empty frequency is mixed 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 just 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, for step S108, 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 just 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, for step S109, 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.
In one embodiment, for step S110, by initial phase and the PH of being averaged for the first timestart-avgCarried out with ± π/4 Compare for formula (13)：
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.
In one embodiment, for step S111, 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.
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 up of 6 grades of rectangular window arithmetic mean filters of 2 kinds of filtering parameters.
When not considering the mixing interfering frequency of orthogonal mixing, orthogonal mixing is expressed as formula (16), and the 6 of 2 kinds of filtering parameters Level rectangular window arithmetic mean filter filtering calculation expression is formula (17)：
In formula, R+end(n) it is positive real frequency mixing sequence again, I+end(n) sequence, R are mixed for weakened body resistance frequency again-end(-n) For anti-real frequency mixing sequence again, I-end(- n) is mixed sequence, cos (ω for anti-empty frequency agains) 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, Ω is signal Frequencies omegaiWith reference frequency ωsFrequency difference, ωiFor signal frequency, T is sampling interval duration, and n is series of discrete number,For again Secondary positive sequence initial phase, β 2 is anti-pleat sequence initial phase, N are predetermined sequence length again.
In formula, R+endFor just 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 ND1Individual 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 ND2Individual 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, Purpose carries out depth suppression to being mixed interfering frequency caused by 1/3 subharmonic；Filtering parameter ND2Value is the reference frequency 2 times of unit period sequence length, purpose carries out depth to being mixed interfering frequency caused by direct current, 1/2 gradation, subharmonic etc. Suppress.6 grades of rectangular window arithmetic mean filters filtering of 2 kinds of filtering parameters, which calculates, needs to use signal period sequence length 10.5 times.
Filtering parameter ND1With filtering parameter ND2Expression formula is formula (18)：
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, for step S112, positive phase and antiphase expression formula is (19) again again：
In formula, PH+endFor positive phase again, PH-endFor antiphase again, R+endFor just 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 for 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.
In one embodiment, for step S113, 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.
In one embodiment, for step S114, 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 is predetermined sequence length,For positive sequence initial phase again, β 2 is anti-pleat sequence initial phase again.
In one embodiment, for step S115, it is (22) to obtain SIN function modulation sequence expression formula：
In formula, Xsin(n) it is SIN function modulation sequence, A is 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 is predetermined sequence length,For positive sequence initial phase again, β 2 is anti-pleat sequence initial phase again.
In one embodiment, for step S116, obtaining zero initial phase benchmark cosine function modulation sequence expression formula is (23)：
In formula, X0cos(n) it is zero initial phase benchmark cosine function modulation sequence, A is cosine function modulation sequence amplitude, list Position v, ωiSignal frequency, T are sampling interval duration, and n is series of discrete number, and N is predetermined sequence length.
In one embodiment, for step S116, obtaining zero initial phase reference sinusoidal FUNCTION MODULATION sequence expression formula is (24)：
In formula, X0sin(n) it is zero initial phase reference sinusoidal FUNCTION MODULATION sequence, A is cosine function modulation sequence amplitude, list Position v, ωiSignal frequency, T are sampling interval duration, and n is series of discrete number, and N is predetermined sequence length.Zero initial phase datum mark Avatars, shown in Fig. 3.
In one embodiment, the default frequency of finely tuning is the arithmetic number less than or equal to actual signal frequency 1%, unit Rad/s, it is expressed as formula (25)：
In formula, ΩsetTo finely tune frequency, unit rad/s, Ωset≤0.01ωi。
For step S117, the first multiplication sequence of acquisition is formula (26)：
For step S117, the second multiplication sequence of acquisition is formula (27)：
For step S117, the 3rd multiplication sequence of acquisition is formula (28)：
For step S117, the 4th multiplication sequence of acquisition is formula (29)：
In formula, X1 (n) is the first multiplication sequence, and X2 (n) is the second multiplication sequence, and X3 (n) is the 3rd multiplication sequence, X4 (n) it is the 4th multiplication sequence, sin (ΩsetTn it is) discrete sine function of the fine setting frequency, cos (ΩsetTn) to be described micro- The discrete cosine function of frequency modulation rate.
In one embodiment, put forward the high-frequency zero initial phase reference sinusoidal sequence of function for step S118, acquisition and be Formula (30)：
In one embodiment, put forward high-frequency zero initial phase benchmark cosine function sequence for step S118, acquisition and be Formula (31)：
In formula, X0sin+f(n) it is the raising frequency zero initial phase reference sinusoidal sequence of function, X0cos+f(n) it is raising frequency zero The initial phase reference sinusoidal sequence of function, sequence frequency improve Ωset。
In one embodiment, for step S119, it is formula (32) that acquisition, which carries high-frequency benchmark orthogonal sequence,：
In formula, sequence X 0sin+fAnd sequence X 0 (n)cos+f(n) it is orthogonal each other.
The system that the raising frequency reference orthogonal sequence of electric power signal is obtained in one embodiment, as shown in figure 4, including：
Preliminary sequence length modules 401, for the lower limit according to frequency power signal scope, preset sample frequency and preset Integer signal period number, obtain the preliminary sequence length of the electric power signal；
Preliminary sequence module 402, for being sampled according to the preliminary sequence length to the electric power signal, obtain institute State the preliminary sequence of electric power signal；
Frequency preliminary survey module 403, for entering line frequency preliminary survey to the preliminary sequence, generate the preliminary of the electric power signal Frequency, and according to the reference frequency of the preliminary frequency setting electric power signal；
Unit period sequence length module 404, 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 405, 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, and the predetermined sequence length is odd number；
First positive block 406, for the default starting point of basis and the predetermined sequence length, from the preliminary sequence The first positive sequence of the electric power signal is obtained in row；
First anti-pleat block 407, for obtaining the first anti-of the electric power signal according to the first positive sequence Pleat sequence；
First phase module 408, for obtaining the first positive phase of 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 409 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 410, for by it is described for the first time averagely initial phase with ± π/4 compared with, obtain 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 411 again, 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 412 again, for the positive phase again according to the positive sequence acquisition again electric power signal, And the antiphase again of the electric power signal is obtained according to the pleat sequence anti-again；
Be averaged initial phase module 413 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 414, for will described in again positive sequence be added with the pleat sequence anti-again, And according to the result after addition and the initial phase that is averaged again, obtain the cosine function modulation sequence of the electric power signal；
SIN function modulation sequence module 415, for the sequence positive again and the pleat sequence anti-again to be subtracted each other, And according to the result after subtracting each other and the initial phase that is averaged again, obtain the SIN function modulation sequence of the electric power signal；
Reference function modulation sequence module 416, for being exported from the cosine function modulation sequence central point, described in acquisition Zero initial phase benchmark cosine function modulation sequence of electric power signal；Export, obtain from the SIN function modulation sequence central point Zero initial phase reference sinusoidal FUNCTION MODULATION sequence of the electric power signal；
Multiplication sequence module 417, for by more than discrete sine function and the zero initial phase benchmark of default fine setting frequency String FUNCTION MODULATION sequence, which is multiplied, obtains the first multiplication sequence of the electric power signal, by the discrete cosine of the default fine setting frequency Function is multiplied with the zero initial phase reference sinusoidal FUNCTION MODULATION sequence obtains the second multiplication sequence of the electric power signal, will be pre- If the discrete sine function of fine setting frequency is multiplied with the zero initial phase reference sinusoidal FUNCTION MODULATION sequence obtains the electric power letter Number the 3rd multiplication sequence, discrete cosine function and the zero initial phase benchmark cosine function of the default fine setting frequency are adjusted Sequence processed, which is multiplied, obtains the 4th multiplication sequence of the electric power signal；
Frequency reference sequence of function module 418 is improved, first multiplication sequence is added with second multiplication sequence, Obtain the electric power signal puies forward the high-frequency zero initial phase reference sinusoidal sequence of function；By the 4th multiplication sequence with it is described 3rd multiplication sequence subtracts each other, and obtain the electric power signal puies forward high-frequency zero initial phase benchmark cosine function sequence；
High-frequency benchmark orthogonal sequence module 419 is put forward, for carrying high-frequency zero initial phase reference sinusoidal according to The sequence of function and described put forward high-frequency zero initial phase benchmark cosine function sequence to obtain carrying for the electric power signal high-frequency Benchmark orthogonal sequence.
It is evidenced from the above discussion that the present invention, which obtains, carries high-frequency benchmark orthogonal sequence, the standard that sine parameter calculates is improved Exactness, meet to be actually needed.
In one embodiment, power system frequency scope in one embodiment, takes C in 45Hz-55Hz2πFor 13. In one embodiment, it is formula (1) that preliminary sequence length modules 401, which obtain the preliminary sequence length,：
In formula, NstartFor preliminary sequence length；(int) represent to round；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.
In one embodiment, the electric power signal be single fundamental frequency cosine function signal, preliminary sequence module 402 The preliminary sequence for obtaining the 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.
In one embodiment, it is formula (3) that frequency preliminary survey module 403, which generates the preliminary frequency,：
ωo(3)；
Preferably, the reference frequency is formula (4) equal to the preliminary frequency：ωs=ωo(4)；
Wherein, ωsFor reference frequency, ωoFor preliminary frequency.
In one embodiment, unit period sequence length module 404 obtains the unit period sequence of the electric power signal Length 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.The error in 1 sampling interval be present in unit period sequence length integer.
In one embodiment, it is formula (6) that predetermined sequence length modules 405, which obtain the predetermined sequence length,：
Wherein, N is predetermined sequence length, and (int) is round numbers, N2πFor the unit period sequence length.
In one embodiment, 0.5 times that starting point is the unit period sequence length is preset；
In one embodiment, it is formula (7) that first positive block 406, which obtains 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, T is sampling interval duration, and n is series of discrete number,It is pre- for first positive sequence initial phase, N If sequence length.First positive sequence pattern expression, shown in Fig. 2.
In one embodiment, it is formula (8) that first anti-pleat block 407, which obtains first anti-pleat sequence,：
Wherein, X-start(- n) is first anti-pleat sequence, X+start(n) it is first positive sequence, β 1 is first anti-pleat sequence Initial phase, N predetermined sequence length.The avatars of the first anti-pleat sequence, as shown in Figure 2.
In one embodiment, first positive phase and first antiphase are the results based on orthogonal mixing and integral and calculating, When not considering the mixing interfering frequency of 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) sequence, R are mixed for first weakened body resistance frequency-start (- n) is that first anti-real frequency is mixed sequence, I-start(- n) is that first anti-empty frequency is mixed sequence, cos (ωs) or cos (- ω TnsTn) For the discrete cosine function of reference frequency, sin (ωs) or sin (- ω TnsTn) it is for the discrete sine function of reference frequency, Ω Signal frequency ωiWith reference frequency ωsFrequency difference, N is predetermined sequence length.
In formula, R+startFirst just 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, T are sampling interval duration, and n is series of discrete number, and N is predetermined sequence length,For first positive sequence initial phase, β 1 is First positive sequence initial phase.
In one embodiment, first phase module 408 obtains first positive phase and the expression formula of first antiphase is (11)：
In formula, PH+startFor first positive phase, PH-startFor first antiphase, R+startFirst just 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, T is sampling interval duration, and N is predetermined sequence length,To be first Secondary positive sequence initial phase, β 1 are first anti-pleat sequence initial phase.
In one embodiment, the expression formula that averagely initial phase module 409 obtains first averagely initial phase for the first time is (12)：
PHstart-avgFor the initial phase that is averaged for the first time, PH+startFor first positive phase, PH-startFor first antiphase.
In one embodiment, phase bits comparison module 410 by it is described for the first time averagely initial phase and PHstart-avgEnter with ± π/4 Row is relatively formula (13)：
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.
In one embodiment, block 411 obtains positive sequence again again and anti-pleat sequence is formula (15) 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.
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 up of 6 grades of rectangular window arithmetic mean filters of 2 kinds of filtering parameters.
When not considering the mixing interfering frequency of orthogonal mixing, orthogonal mixing is expressed as formula (16), and the 6 of 2 kinds of filtering parameters Level rectangular window arithmetic mean filter filtering calculation expression is formula (17)：
In formula, R+end(n) it is positive real frequency mixing sequence again, I+end(n) sequence, R are mixed for weakened body resistance frequency again-end(-n) For anti-real frequency mixing sequence again, I-end(- n) is mixed sequence, cos (ω for anti-empty frequency agains) 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.
In formula, R+endFor just 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；K (Ω) is digital filtering in frequency difference Ω amplitude gain, unit dimensionless；ND1For filtering parameter 1, i.e., pair ND1Individual continuous centrifugal pump is added, and then takes its arithmetic mean of instantaneous value to be exported as this filter value；ND2For filtering parameter 2, i.e., to ND2 Individual continuous centrifugal pump is added, and then takes its arithmetic mean of instantaneous value to be exported as this filter value；NDGrown for digital filtering using sequence Spend, be the summation of 6 grades of rectangular window arithmetic mean filter filtering parameters in quantity, 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 suppression to being mixed interfering frequency caused by 1/3 subharmonic；Filtering parameter ND2Value is the reference frequency 2 times of unit period sequence length, purpose carries out depth to being mixed interfering frequency caused by direct current, 1/2 gradation, subharmonic etc. Suppress.6 grades of rectangular window arithmetic mean filters filtering of 2 kinds of filtering parameters, which calculates, needs to use signal period sequence length 10.5 times.
Filtering parameter ND1With filtering parameter ND2Expression formula is formula (18)：
In one embodiment, phase module 412 obtains positive phase again again and antiphase expression formula is (19) again：
In formula, PH+endFor positive phase again, PH-endFor antiphase again, R+endFor just 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, NDSequence length is used for digital filtering,For positive sequence again Initial phase, β 2 are anti-pleat sequence initial phase again.
In one embodiment, the expression formula that averagely initial phase module 413 obtains the initial phase that is averaged again again is (20)：
PHend-avgFor the initial phase that is averaged again, PH+endFor positive phase again, PH-endFor antiphase again
In one embodiment, the acquisition of cosine function modulation sequence module 414 cosine function modulation sequence expression formula is (21)：
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.
In one embodiment, the acquisition of SIN function modulation sequence module 415 SIN function modulation sequence expression formula is (22)：
In formula, Xsin(n) it is SIN function modulation sequence, A is SIN function modulation sequence amplitude, unit v,For cosine function modulation sequence initial phase.
In one embodiment, reference function modulation sequence module 416 obtains zero initial phase benchmark cosine function modulation sequence Row expression formula is (23)：
In formula, X0cos(n) it is zero initial phase benchmark cosine function modulation sequence, A is cosine function modulation sequence amplitude, list Position v.
In one embodiment, reference function modulation sequence module 416 obtains zero initial phase reference sinusoidal FUNCTION MODULATION sequence Row expression formula is (24)：
In formula, X0sin(n) it is zero initial phase reference sinusoidal FUNCTION MODULATION sequence, A is cosine function modulation sequence amplitude, list Position v.Zero initial phase datum mark avatars, shown in Fig. 3.
In one embodiment, the default frequency of finely tuning is the arithmetic number less than or equal to actual signal frequency 1%, unit Rad/s, it is expressed as formula (25)：
In formula, ΩsetTo finely tune frequency, unit rad/s, Ωset≤0.01ωi。
It is formula (26) that multiplication sequence module 417, which obtains the first multiplication sequence,：
It is formula (27) that multiplication sequence module 417, which obtains the second multiplication sequence,：
It is formula (28) that multiplication sequence module 417, which obtains the 3rd multiplication sequence,：
The multiplication sequence of multiplication sequence module 417 the 4th is formula (29)：
In one embodiment, improve the acquisition of frequency reference sequence of function module 418 and carry high-frequency zero initial phase benchmark SIN function sequence is formula (30)：
In one embodiment, improve the acquisition of frequency reference sequence of function module 418 and carry high-frequency zero initial phase benchmark Cosine function sequence is formula (31)：
In formula, X0sin+f(n) it is the raising frequency zero initial phase reference sinusoidal sequence of function, X0cos+f(n) it is raising frequency zero The initial phase reference sinusoidal sequence of function, sequence frequency improve Ωset。
In one embodiment, propose the high-frequency acquisition of benchmark orthogonal sequence module 419 and put forward the orthogonal sequence of high-frequency benchmark It is classified as formula (32)：
In formula, sequence X 0sin+fAnd sequence X 0 (n)cos+f(n) it is orthogonal each other.
Each technical characteristic of embodiment described above can be combined arbitrarily, to make description succinct, not to above-mentioned reality Apply 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 its description is more specific and detailed, but simultaneously Can not therefore it be construed as limiting the scope of the patent.It should be pointed out that come for one of ordinary skill in the art Say, 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 (6)

  1. A kind of 1. method for the raising frequency reference orthogonal sequence for obtaining electric power signal, it is characterised in that comprise the following steps：
    According to the lower limit of frequency power signal scope, 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；
    Line frequency preliminary survey is entered to the preliminary sequence, generates the preliminary frequency of the electric power signal, 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, the predetermined sequence length are odd number；
    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；
    By the averagely initial phase for the first time compared with ± π/4, the first phase fiducial value compared with ± π/4 is obtained, 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 according to the pleat sequence anti-again of the positive sequence acquisition again electric power signal；
    The positive phase again of the electric power signal is obtained according to the sequence positive again, and 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 positive sequence again is added with the pleat sequence anti-again, and according to the result after addition and it is described be averaged again it is first Phase, the cosine function modulation sequence of the electric power signal is obtained, wherein, according to expression formula Obtain the cosine function modulation sequence Xcos(n), n=0,1,2 ..., N-1, X+end(n) it is the sequence positive again, X-end (- n) pleat sequence anti-again, PH for described inend-avgFor the initial phase that is averaged again, n is series of discrete number, and N is described default Sequence length；
    The sequence positive again and the pleat sequence anti-again are subtracted each other, and according to the result after subtracting each other and it is described again it is average just Phase, the SIN function modulation sequence of the electric power signal is obtained, wherein, according to expression formula Obtain the SIN function modulation sequence Xsin(n), n=0,1,2 ..., N-1, X+end(n) it is the sequence positive again, X-end (- n) pleat sequence anti-again, PH for described inend-avgFor the initial phase that is averaged again, n is series of discrete number, and N is described default Sequence length；
    Exported from the cosine function modulation sequence central point, the zero initial phase benchmark cosine function for obtaining the electric power signal is adjusted Sequence processed；Exported from the SIN function modulation sequence central point, obtain zero initial phase reference sinusoidal letter of the electric power signal Number modulation sequence；
    The discrete sine function of default fine setting frequency is multiplied with the zero initial phase benchmark cosine function modulation sequence and obtains institute The first multiplication sequence of electric power signal is stated, by discrete cosine function and the zero initial phase benchmark of the default fine setting frequency just String FUNCTION MODULATION sequence, which is multiplied, obtains the second multiplication sequence of the electric power signal, by the discrete sine function of default fine setting frequency It is multiplied with the zero initial phase reference sinusoidal FUNCTION MODULATION sequence and obtains the 3rd multiplication sequence of the electric power signal, will be described pre- If the discrete cosine function of fine setting frequency is multiplied with the zero initial phase benchmark cosine function modulation sequence obtains the electric power letter Number the 4th multiplication sequence；
    First multiplication sequence is added with second multiplication sequence, obtains carrying at the beginning of high-frequency zero for the electric power signal Phase reference SIN function sequence；4th multiplication sequence is subtracted each other with the 3rd multiplication sequence, obtains the electric power letter Number put forward high-frequency zero initial phase benchmark cosine function sequence；
    Put forward the high-frequency zero initial phase reference sinusoidal sequence of function according to described and described carry more than high-frequency zero initial phase benchmark The string sequence of function acquisition electric power signal carries high-frequency benchmark orthogonal sequence.
  2. 2. the method for the raising frequency reference orthogonal sequence according to claim 1 for obtaining electric power signal, it is characterised in that The electric power signal is the cosine function signal of single fundamental frequency, according to expression formulaObtain institute State preliminary sequence Xstart(n), whereinN=0,1,2,3 ..., Nstart- 1, A are signal amplitude, ωiFor signal frequency,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. 3. the method for the raising frequency reference orthogonal sequence according to claim 1 for obtaining electric power signal, it is characterised in that According to expression formula X1 (n)=X0cos(n)sin(ΩsetTn the first multiplication sequence X1 (n)) is obtained, according to expression formula X2 (n) =X0sin(n)cos(ΩsetTn the second multiplication sequence X2 (n)) is obtained, according to expression formula X3 (n)=X0sin(n)sin (ΩsetTn the 3rd multiplication sequence X3 (n)) is obtained, according to expression formula X4 (n)=X0cos(n)cos(ΩsetTn) described in obtaining 4th multiplication sequence X4 (n), whereinX0cos(n) modulated for the zero initial phase benchmark cosine function Sequence, X0sin(n) it is the zero initial phase reference sinusoidal FUNCTION MODULATION sequence, sin (ΩsetTn) it is the default fine setting frequency Discrete sine function, cos (ΩsetTn it is) the discrete cosine function of the default fine setting frequency, ΩsetFor the default fine setting Frequency, T are sampling interval duration, and n is series of discrete number, and N is the predetermined sequence length.
  4. A kind of 4. system for the raising frequency reference orthogonal sequence for obtaining electric power signal, it is characterised in that including：
    Preliminary sequence length modules, believe for the lower limit according to frequency power signal scope, preset sample frequency and default integer Number periodicity, obtain the preliminary sequence length of the electric power signal；
    Preliminary sequence module, for being sampled according to the preliminary sequence length to the electric power signal, obtain the electric power The preliminary sequence of signal；
    Frequency preliminary survey module, for entering line frequency preliminary survey to the preliminary sequence, the preliminary frequency of the electric power signal is generated, 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 The predetermined sequence length of the electric power signal is obtained, the predetermined sequence length is odd number；
    First positive block, for according to default starting point and the predetermined sequence length, being obtained from the preliminary sequence Obtain the first positive sequence of the electric power signal；
    First anti-pleat block, for obtaining the first anti-pleat sequence of the electric power signal according to the first positive sequence；
    First phase module, for the first positive phase according to the first positive sequence acquisition electric power signal, 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, for the averagely initial phase for the first time compared with ± π/4, to be 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；
    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 according to the pleat sequence anti-again of the positive sequence acquisition again electric power signal Row；
    Phase module again, for the positive phase again according to the positive sequence acquisition again electric power signal, 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 will described in again positive sequence be added with the pleat sequence anti-again, and according to Result and the initial phase that is averaged again after addition, the cosine function modulation sequence of the electric power signal is obtained, wherein, it is described Cosine function modulation sequence module is according to expression formulaObtain the cosine function modulation sequence Arrange Xcos(n), 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；
    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 and the initial phase that is averaged again after subtracting each other, the SIN function modulation sequence of the electric power signal is obtained, wherein, it is described SIN function modulation sequence module is according to expression formulaObtain the SIN function modulation sequence Arrange Xsin(n), 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；
    Reference function modulation sequence module, for being exported from the cosine function modulation sequence central point, obtain the electric power letter Number zero initial phase benchmark cosine function modulation sequence；Exported from the SIN function modulation sequence central point, obtain the electricity Zero initial phase reference sinusoidal FUNCTION MODULATION sequence of force signal；
    Multiplication sequence module, for the discrete sine function of default fine setting frequency to be adjusted with the zero initial phase benchmark cosine function Sequence processed, which is multiplied, obtains the first multiplication sequence of the electric power signal, by the discrete cosine function of the default fine setting frequency and institute The second multiplication sequence that the multiplication of zero initial phase reference sinusoidal FUNCTION MODULATION sequence obtains the electric power signal is stated, by default fine setting frequency The discrete sine function of rate is multiplied with the zero initial phase reference sinusoidal FUNCTION MODULATION sequence obtains the 3rd of the electric power signal Multiplication sequence, by the discrete cosine function of the default fine setting frequency and the zero initial phase benchmark cosine function modulation sequence phase Multiply the 4th multiplication sequence for obtaining the electric power signal；
    Frequency reference sequence of function module is improved, first multiplication sequence is added with second multiplication sequence, obtains institute That states electric power signal puies forward the high-frequency zero initial phase reference sinusoidal sequence of function；4th multiplication sequence is multiplied with the described 3rd Method sequence is subtracted each other, and obtain the electric power signal puies forward high-frequency zero initial phase benchmark cosine function sequence；
    High-frequency benchmark orthogonal sequence module is put forward, for putting forward the high-frequency zero initial phase reference sinusoidal sequence of function according to With it is described carry high-frequency zero initial phase benchmark cosine function sequence obtain the electric power signal to carry high-frequency benchmark orthogonal Sequence.
  5. 5. the system of the raising frequency reference orthogonal sequence according to claim 4 for obtaining electric power signal, it is characterised 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), whereinN=0,1,2,3 ..., Nstart- 1, A is signal amplitude, ωiFor signal frequency,For the initial phase of the preliminary sequence, T is sampling interval duration, and f is described pre- If sample frequency, n is series of discrete number, NstartFor the preliminary sequence length.
  6. 6. the system of the raising frequency reference orthogonal sequence according to claim 4 for obtaining electric power signal, it is characterised in that The multiplication sequence module is according to expression formula X1 (n)=X0cos(n)sin(ΩsetTn the first multiplication sequence X1 (n)) is obtained, According to expression formula X2 (n)=X0sin(n)cos(ΩsetTn the second multiplication sequence X2 (n)) is obtained, according to expression formula X3 (n) =X0sin(n)sin(ΩsetTn the 3rd multiplication sequence X3 (n)) is obtained, according to expression formula X4 (n)=X0cos(n)cos (ΩsetTn the 4th multiplication sequence X4 (n)) is obtained, whereinX0cos(n) it is zero initial phase Benchmark cosine function modulation sequence, X0sin(n) it is the zero initial phase reference sinusoidal FUNCTION MODULATION sequence, sin (ΩsetTn) it is The discrete sine function of the default fine setting frequency, cos (ΩsetTn it is) the discrete cosine function of the default fine setting frequency, ΩsetFor the default fine setting frequency, T is sampling interval duration, and n is series of discrete number, and N is the predetermined sequence length.
CN201510891667.XA 2015-12-04 2015-12-04 Obtain the method and system of the raising frequency reference orthogonal sequence of electric power signal Active CN105548705B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510891667.XA CN105548705B (en) 2015-12-04 2015-12-04 Obtain the method and system of the raising frequency reference orthogonal sequence of electric power signal

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510891667.XA CN105548705B (en) 2015-12-04 2015-12-04 Obtain the method and system of the raising frequency reference orthogonal sequence of electric power signal

Publications (2)

Publication Number Publication Date
CN105548705A CN105548705A (en) 2016-05-04
CN105548705B true CN105548705B (en) 2018-01-30

Family

ID=55828026

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510891667.XA Active CN105548705B (en) 2015-12-04 2015-12-04 Obtain the method and system of the raising frequency reference orthogonal sequence of electric power signal

Country Status (1)

Country Link
CN (1) CN105548705B (en)

Citations (9)

* Cited by examiner, † Cited by third party
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
CN104635045A (en) * 2015-02-05 2015-05-20 广东电网有限责任公司电力科学研究院 Power signal frequency detection method and system based on phase modulation
CN104635044A (en) * 2015-02-05 2015-05-20 广东电网有限责任公司电力科学研究院 Power signal frequency detection method and system based on amplitude modulation
CN105067880A (en) * 2015-09-18 2015-11-18 广东电网有限责任公司电力科学研究院 Method and system for orthogonal modulation of power signals

Patent Citations (9)

* Cited by examiner, † Cited by third party
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
CN104635045A (en) * 2015-02-05 2015-05-20 广东电网有限责任公司电力科学研究院 Power signal frequency detection method and system based on phase modulation
CN104635044A (en) * 2015-02-05 2015-05-20 广东电网有限责任公司电力科学研究院 Power signal frequency detection method and system based on amplitude 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
CN105548705A (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
CN105548699B (en) Obtain the method and system of the raising frequency reference SIN function 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
CN105548700B (en) Obtain the method and system of the raising frequency orthogonal sequence of electric power signal
CN105548701B (en) Obtain the method and system of the reduction frequency orthogonal sequence of electric power signal
CN105548702B (en) Obtain the method and system for dropping low-frequency zero initial phase benchmark cosine function sequence
CN105467212B (en) Obtain the method and system of the raising frequency cosine function sequence of electric power signal
CN105467211B (en) Obtain the method and system of the raising frequency reference cosine function sequence of electric power signal
CN105425035B (en) Obtain the method and system for dropping low-frequency zero initial phase reference sinusoidal sequence of function
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
CN105092970B (en) Obtain the method and system of the initial phase datum mark of electric power signal sequence SIN function zero
CN105548706B (en) Obtain the method and system of the reduction frequency sine sequence of function of electric power signal
CN105425031B (en) The orthogonal frequency multiplication sequence method of zero initial phase benchmark and system are obtained from electric power signal
CN105067885B (en) Electric power signal is converted to the method and system of zero initial phase signal sequence
CN105548698B (en) Frequency power signal method for trimming and system
CN105203843B (en) The average first phase position detecting method and system of electric power signal
CN105548687B (en) The method and system of any initial phase orthogonal sequence is obtained from electric power signal
CN105548703B (en) The method and system of any initial phase SIN function sequence is obtained from electric power signal
CN105445551B (en) The cut-off method for detecting phases and system of sinusoidal signal
CN105548693B (en) Obtain the method and system of zero initial phase reference sinusoidal function times frequency sequence
CN105372490B (en) The method and system of frequency measurement based on zero orthogonal times of frequency sequence of initial phase benchmark
CN105301356B (en) Obtain the method and system of zero initial phase benchmark cosine function times frequency sequence

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