CN105203842B - The method and system of SIN function reference signal is extracted from electric power signal - Google Patents
The method and system of SIN function reference signal is extracted from electric power signal Download PDFInfo
- Publication number
- CN105203842B CN105203842B CN201510599156.0A CN201510599156A CN105203842B CN 105203842 B CN105203842 B CN 105203842B CN 201510599156 A CN201510599156 A CN 201510599156A CN 105203842 B CN105203842 B CN 105203842B
- Authority
- CN
- China
- Prior art keywords
- sequence
- frequency
- electric power
- power signal
- 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
Links
Landscapes
- Electrotherapy Devices (AREA)
- Inverter Devices (AREA)
Abstract
The invention discloses a kind of method and system that SIN function reference signal is extracted from electric power signal, methods described includes：According to frequency power signal, preset sample frequency and integer signal period number, predetermined sequence length is obtained；Electric power signal is sampled, obtains positive, anti-pleat sequence, obtains zero initial phase SIN function modulation sequence；Frequency power signal cosine, SIN function are multiplied with positive sequence respectively, the real frequency of generation first, empty frequency sequence vector, then be multiplied respectively with anti-pleat sequence, the real frequency of generation second, empty frequency sequence；By carrying out digital notch, integration, phase transition to above-mentioned sequence vector respectively, first, second phase is obtained, then is converted into the average initial phase of positive sequence, amplitude is obtained and recovers sequence；Point on the basis of the central point for recovering sequence by amplitude, obtains SIN function reference signal sequence.The present invention can extract SIN function reference signal sequence from Physical Experiment signal, determine its zero initial phase datum mark, be adapted to application.
Description
Technical field
The present invention relates to technical field of electric power, and SIN function reference signal is extracted from electric power signal more particularly to one kind
Method and system.
Background technology
In power science experimental study, the determination of the initial phase datum mark of sinusoidal signal zero is to testing and research has in itself
Important meaning, in Simulating Test Study, the problem of any is not present in the determination of zero initial phase datum mark.But in physics
In experimental study, for discrete acquisitions data, the determination of Physical Experiment signal zero initial phase datum mark is not an easy thing, zero
The result of the uncertain influence Physical Experiment research of initial phase datum mark.
The content of the invention
Based on the above situation, the present invention propose a kind of method that SIN function reference signal is extracted from electric power signal and
System, SIN function reference signal sequence can be extracted from Physical Experiment signal, determine its zero initial phase datum mark, met actual
Need.
To achieve these goals, the embodiment of technical solution of the present invention is：
A kind of method that SIN function reference signal is extracted from electric power signal, comprises the following steps：
According to the frequency of electric power signal, preset sample frequency and default integer signal period number, the electric power signal is obtained
Predetermined sequence length, the predetermined sequence length is odd number；
The electric power signal is sampled according to the predetermined sequence length, obtains the positive sequence of the electric power signal
Row；
The positive sequence is reversely exported, obtains anti-pleat sequence corresponding with the positive sequence；
The positive sequence and the anti-pleat sequence are subtracted each other, the SIN function for obtaining the initial phase of electric power signal zero is adjusted
Sequence processed；
Respectively by the SIN function of the cosine function of the frequency power signal and the frequency power signal with it is described just
It is multiplied to sequence, generates the first real frequency sequence vector and the first empty frequency sequence vector；
Respectively by the SIN function of the cosine function of the frequency power signal and the frequency power signal with it is described anti-
Pleat sequence is multiplied, and generates the second real frequency sequence vector and the second empty frequency sequence vector；
Digital notch, generation first are carried out to the described first real frequency sequence vector and the first empty frequency sequence vector respectively
Real frequently vectorial trap sequence and the first empty vectorial trap sequence of frequency；
Integral operation is carried out to the described first real vectorial trap sequence of frequency and the first empty vectorial trap sequence of frequency respectively,
The vectorial integrated value of the real frequency of generation first and the first empty vectorial integrated value of frequency；
Digital notch, generation second are carried out to the described second real frequency sequence vector and the second empty frequency sequence vector respectively
Real frequently vectorial trap sequence and the second empty vectorial trap sequence of frequency；
Integral operation is carried out to the described second real vectorial trap sequence of frequency and the second empty vectorial trap sequence of frequency respectively,
The vectorial integrated value of the real frequency of generation second and the second empty vectorial integrated value of frequency；
According to default phase transition rule, the described first empty vectorial integrated value of frequency and the described first real frequency vector are integrated
Value is converted to first phase, and the described second empty vectorial integrated value of frequency and the described second real vectorial integrated value of frequency are converted into the second phase
Position；
According to default sequence average initial phase transformation rule, the first phase and the second phase are converted into institute
State the average initial phase of positive sequence；
According to the average initial phase of the positive sequence, by the SIN function modulation sequence of the initial phase of electric power signal zero
Amplitude revert to the amplitude of the electric power signal, obtain amplitude and recover sequence；
Point on the basis of the central point for recovering sequence by the amplitude, obtain the SIN function reference signal of the electric power signal
Sequence.
A kind of system that SIN function reference signal is extracted from electric power signal, including：
Predetermined sequence length modules, for the frequency according to electric power signal, preset sample frequency and default integer signal week
Issue, obtains the predetermined sequence length of the electric power signal, and the predetermined sequence length is odd number；
Positive block, for being sampled according to the predetermined sequence length to the electric power signal, described in acquisition
The positive sequence of electric power signal；
Anti- pleat block, for the positive sequence reversely to be exported, obtain anti-pleat corresponding with the positive sequence
Sequence；
SIN function modulation sequence module, for the positive sequence and the anti-pleat sequence to be subtracted each other, obtain the electricity
The SIN function modulation sequence of the initial phase of force signal zero；
Primary vector sequence generating module, for respectively believing the cosine function of the frequency power signal and the electric power
The SIN function of number frequency is multiplied with the positive sequence, generates the first real frequency sequence vector and the first empty frequency sequence vector；
Secondary vector sequence generating module, for respectively believing the cosine function of the frequency power signal and the electric power
The SIN function of number frequency is multiplied with the anti-pleat sequence, generates the second real frequency sequence vector and the second empty frequency sequence vector；
Primary vector trap sequence generating module, for respectively to the described first real frequency sequence vector and the first empty frequency
Sequence vector carries out digital notch, the vectorial trap sequence of the real frequency of generation first and the first empty vectorial trap sequence of frequency；
Primary vector integrated value generation module, for respectively to the described first real vectorial trap sequence of frequency and first void
Frequently vectorial trap sequence carries out integral operation, the vectorial integrated value of the real frequency of generation first and the first empty vectorial integrated value of frequency；
Secondary vector trap sequence generating module, for respectively to the described second real frequency sequence vector and the second empty frequency
Sequence vector carries out digital notch, the vectorial trap sequence of the real frequency of generation second and the second empty vectorial trap sequence of frequency；
Secondary vector integrated value generation module, for respectively to the described second real vectorial trap sequence of frequency and second void
Frequently vectorial trap sequence carries out integral operation, the vectorial integrated value of the real frequency of generation second and the second empty vectorial integrated value of frequency；
Phase module, will the vectorial integrated value of the described first empty frequency and described for according to default phase transition rule
The one real vectorial integrated value of frequency is converted to first phase, and the described second empty vectorial integrated value of frequency and the described second real frequency vector are integrated
Value is converted to second phase；
Average initial phase module, for according to default sequence average initial phase transformation rule, by the first phase and
The second phase is converted to the average initial phase of the positive sequence；
Amplitude recovers block, for the average initial phase according to the positive sequence, by the beginning of the electric power signal zero
The amplitude of the SIN function modulation sequence of phase reverts to the amplitude of the electric power signal, obtains amplitude and recovers sequence；
SIN function reference signal module, for sequence is recovered by the amplitude central point on the basis of point, described in acquisition
The SIN function reference signal sequence of electric power signal.
Compared with prior art, beneficial effects of the present invention are：The present invention extracts SIN function benchmark from electric power signal
The method and system of signal, according to the frequency of electric power signal, preset sample frequency and default integer signal period number, preset
Sequence length；Electric power signal is sampled according to predetermined sequence length, obtains positive sequence；Positive sequence is reversely exported,
Obtain anti-pleat sequence；Positive sequence and anti-pleat sequence are subtracted each other, obtain the SIN function modulation sequence of zero initial phase；Respectively by electricity
The cosine function and SIN function of force signal frequency are multiplied with positive sequence, generate the first real frequency sequence vector and the first empty frequency to
Sequence is measured, then is multiplied respectively with anti-pleat sequence, generates the second real frequency sequence vector and the second empty frequency sequence vector；By right respectively
The sequence vector of above-mentioned generation carries out digital notch, integration, phase transition, and then obtains the first phase and second phase；
First phase and second phase are converted to the average initial phase of positive sequence again；According to the average initial phase of positive sequence, obtain
Obtain amplitude and recover sequence；Point on the basis of the central point for recovering sequence by amplitude, obtain the higher SIN function benchmark letter of accuracy
Number sequence, determines its zero initial phase datum mark, is adapted to practical application.
Brief description of the drawings
Fig. 1 is the method flow diagram for extracting SIN function reference signal in one embodiment from electric power signal；
Fig. 2 is positive sequence and anti-pleat sequence diagram in one embodiment；
Fig. 3 is zero initial phase datum mark schematic diagram of SIN function reference signal sequence in one embodiment；
Fig. 4 is the system structure diagram for extracting SIN function reference signal in one embodiment from electric power signal.
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.
The method for extracting SIN function reference signal in one embodiment from electric power signal, as shown in figure 1, including following
Step：
Step S101：According to the frequency of electric power signal, preset sample frequency and default integer signal period number, described in acquisition
The predetermined sequence length of electric power signal, the predetermined sequence length are odd number；
Step S102：The electric power signal is sampled according to the predetermined sequence length, obtains the electric power signal
Positive sequence；
Step S103：The positive sequence is reversely exported, obtains anti-pleat sequence corresponding with the positive sequence；
Step S104：The positive sequence and the anti-pleat sequence are subtracted each other, obtain the initial phase of electric power signal zero
SIN function modulation sequence；
Step S105：Respectively by the cosine function of the frequency power signal and the SIN function of the frequency power signal
It is multiplied with the positive sequence, generates the first real frequency sequence vector and the first empty frequency sequence vector；
Step S106：Respectively by the cosine function of the frequency power signal and the SIN function of the frequency power signal
It is multiplied with the anti-pleat sequence, generates the second real frequency sequence vector and the second empty frequency sequence vector；
Step S107：Are carried out by numeral and is fallen into for the described first real frequency sequence vector and the first empty frequency sequence vector respectively
Ripple, the vectorial trap sequence of the real frequency of generation first and the first empty vectorial trap sequence of frequency；
Step S108：The described first real vectorial trap sequence of frequency and the first empty vectorial trap sequence of frequency are carried out respectively
Integral operation, the vectorial integrated value of the real frequency of generation first and the first empty vectorial integrated value of frequency；
Step S109：Are carried out by numeral and is fallen into for the described second real frequency sequence vector and the second empty frequency sequence vector respectively
Ripple, the vectorial trap sequence of the real frequency of generation second and the second empty vectorial trap sequence of frequency；
Step S110：The described second real vectorial trap sequence of frequency and the second empty vectorial trap sequence of frequency are carried out respectively
Integral operation, the vectorial integrated value of the real frequency of generation second and the second empty vectorial integrated value of frequency；
Step S111：It is according to default phase transition rule, the described first empty vectorial integrated value of frequency and described first is real
Frequently vectorial integrated value is converted to first phase, and the described second empty vectorial integrated value of frequency and the described second real vectorial integrated value of frequency are turned
It is changed to second phase；
Step S112：According to default sequence average initial phase transformation rule, by the first phase and second phase
Position is converted to the average initial phase of the positive sequence；
Step S113：According to the average initial phase of the positive sequence, by the sinusoidal letter of the initial phase of electric power signal zero
The amplitude of number modulation sequence reverts to the amplitude of the electric power signal, obtains amplitude and recovers sequence；
Step S114：Point on the basis of the central point for recovering sequence by the amplitude, obtain the sinusoidal letter of the electric power signal
Base calibration signal sequence.
It is evidenced from the above discussion that the present invention extracts accurately SIN function reference signal sequence in Physical Experiment signal,
There is important use on power science experimental study.
Wherein, for step S101, according to the frequency of electric power signal, preset sample frequency and default integer signal period
Number, the predetermined sequence length of electric power signal is obtained, and require that default integer signal period number is even number and predetermined sequence length
For odd number, predetermined sequence length correspondingly presets integer signal period number in principle；
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.
The frequency of usual Physical Experiment electric power signal is known, and in one embodiment, the signal frequency is expressed as
Formula (1)：
ω (1)；
In formula, ω is the signal frequency of rad/s units；
In one embodiment, it is formula (2) to obtain predetermined sequence length：
In formula, N is the predetermined sequence length of dimensionless unit；(int) round numbers is represented；C2πTo preset the integer signal period
Number；N2πFor the signal unit periodic sequence length of dimensionless unit；fnFor Hz unit preset sample frequencies；ω is rad/s units
Signal frequency；
In one embodiment, integer signal period number C is preset2π=12, preset sample frequency fn=10KHz.
For step S102, the electric power signal is sampled according to predetermined sequence length, obtains the electric power signal
Positive sequence；
In one embodiment, the electric power signal is the SIN function signal of single fundamental frequency, obtains the positive sequence
It is classified as formula (3)：
Wherein, X+(n) it is positive sequence；A is signal amplitude, unit v；ω is the signal frequency of rad/s units；TnIt is mono- for s
The sampling interval of position；fnFor the sample frequency of Hz units；N is the series of discrete number of immeasurable unit；At the beginning of signal for rad units
Phase；N is the positive sequence length of dimensionless unit, equal to predetermined sequence length.
The avatars of the positive sequence, as shown in Figure 2.
For step S103, the positive sequence is reversely exported, obtains anti-pleat sequence corresponding with the positive sequence；
In one embodiment, relatively positive sequence, it is formula (4) to obtain anti-pleat sequence：
X-(- n)=X+(N-n)=Asin (- ω Tnn+β)
(4)；
N=0,1,2,3 ..., N-1
In formula, X-(- n) is anti-pleat sequence；β is the anti-pleat sequence initial phase of rad units, in relation, anti-pleat sequence first phase
Position is the cut-off phase of positive sequence；N is the series of discrete number of immeasurable unit；N is the anti-pleat sequence length of dimensionless unit, is equal to
Positive sequence length.
The avatars of the anti-pleat sequence, as shown in Figure 2.
For step S104, the positive sequence and the anti-pleat sequence are subtracted each other, obtain the first phase of electric power signal zero
The SIN function modulation sequence of position, abbreviation modulation sequence；
In one embodiment, the SIN function modulation sequence for obtaining zero initial phase is formula (5)：
In formula, Xsin(n) be zero initial phase SIN function modulation sequence；AsinFor above-mentioned modulation sequence amplitude, unit v；For above-mentioned modulation sequence initial phase, unit rad；For the positive sequence initial phase of rad/s units；β is rad/s
The anti-pleat sequence initial phase of unit.
Allow positive sequence initial phaseExcursion with anti-pleat sequence initial phase β is in 0~± 0.375 π rad.
Error be present because the predetermined sequence length corresponds to integer signal period number, mainly described predetermined sequence length
Integer error.If described error is zero, above-mentioned modulation sequence initial phaseIt is zero, on the contrary initial phaseNear zero.The initial phaseError and the mistake of the integer signal period number compared with null value
It is proportional relation between difference.
For step S105, respectively by the cosine function of the frequency power signal and the sine of the frequency power signal
Function is multiplied with the positive sequence, generates the first real frequency sequence vector and the first empty frequency sequence vector；
In one embodiment, when not considering the mixing interfering frequency composition, it is formula (6) to obtain being mixed sequence：
In formula, R1(n) it is the first real frequency sequence vector；I1(n) it is the first empty frequency sequence vector；N be dimensionless unit just
To sequence length.
For step S106, respectively by the cosine function of the frequency power signal and the sine of the frequency power signal
Function is multiplied with the anti-pleat sequence, generates the second real frequency sequence vector and the second empty frequency sequence vector；
In one embodiment, when not considering the mixing interfering frequency composition, it is formula (7) to obtain being mixed sequence：
In formula, R2(n) it is the second real frequency sequence vector；I2(n) it is the second empty frequency sequence vector；N is the anti-of dimensionless unit
Pleat sequence length.
For step S107, numeral is carried out to the described first real frequency sequence vector and the first empty frequency sequence vector respectively
Trap, the vectorial trap sequence of the real frequency of generation first and the first empty vectorial trap sequence of frequency；
Mixing interfering frequency is included in the real frequency sequence vector and the empty frequency sequence vector.When also existing in input signal
When flip-flop, subharmonic composition and subharmonic composition, the mixing interfering frequency will be more complicated, these mixing interference frequencies
Rate has a strong impact on accuracy in computation.Although window function and integral operation there is good decay to make mixing interfering frequency in itself
With, but there is no specific aim, it is impossible to the inhibitory action of depth is produced to the complicated mixing interfering frequency, it is impossible to meet sinusoidal
The high accuracy of parameter, which calculates, to be needed.
In order to targetedly suppress the influence of the mixing interfering frequency, using a kind of digital trap, ideal situation
Under, the null Frequency point of digital trap just corresponds to the mixing interfering frequency point, has to the mixing interfering frequency
Complete inhibitory action.Preferably, digital notch specifically uses arithmetic average notch algorithm, i.e., by several continuous centrifugal pump phases
Add, then take its arithmetic mean of instantaneous value to be exported as this trap value.Digital notch needs to set digital notch parameter, the numeral
Notch parameter refers to the length N that several continuous centrifugal pumps are addedD.In digital notch parameter NDValue is signal unit periodic sequence
1.5 times of length, it can suppress to being mixed interfering frequency caused by 1/3 subharmonic.And NDValue is signal unit week
2 times of phase sequence length, can to direct current, 1/2 by several times, 1 time, 2 times, 3 times, 4 times, mixing interference frequency caused by 5 subharmonic etc.
Rate is suppressed.Therefore, digital notch is made up of the digital trap of 2 kinds of parameters, it is contemplated that the factor such as physical presence error,
Influenceed in order to which depth suppresses mixing interfering frequency, the digital trap of every kind of parameter is by parameter identical three-level digital notch group
Into totally six grades of arithmetic mean of instantaneous value digital notch are formed.
In one embodiment, six grades of arithmetic mean of instantaneous value digital notch formulas can be formula (8)：
To X (n) n=0,1,2,3 ..., N-1
To XD (n) n=0,1,2,3 ..., N-3ND1-3ND2-1
Wherein, X (n) is digital notch list entries, sequence length N；XD(n) it is digital notch output sequence, sequence length
Spend N-3ND1-3ND2；ND1For notch parameter 1, i.e., continuous centrifugal pump is added quantity；ND2For notch parameter 2, i.e., continuous centrifugal pump phase
Addend amount.
In one embodiment, notch parameter ND1Value is 1.5 times of signal unit periodic sequence length, notch parameter
ND2Value is 2 times of signal unit periodic sequence length, and six grades of arithmetic mean of instantaneous value digital notch need to use 10.5 times of signal lists
Bit period sequence length.
In one embodiment, under the premise of the mixing interfering frequency composition is completely suppressed, the described first real frequency
Vectorial trap sequence and the first empty vectorial trap sequence of frequency are (9)：
To R1(n)I1(n) n=0,1,2,3 ..., N-1
To RD1(n)ID1(n) n=0,1,2,3 ..., N-3ND1-3ND2-1
Wherein, RD1(n) it is the described first real vectorial trap sequence of frequency；ID1(n) it is the described first empty vectorial trap sequence of frequency
Row.
For step S108, respectively to the described first real vectorial trap sequence of frequency and the first empty vectorial trap sequence of frequency
Carry out integral operation, the vectorial integrated value of the real frequency of generation first and the first empty vectorial integrated value of frequency；
In one embodiment, the vectorial integrated value of the real frequency of generation first and the first empty vectorial integrated value of frequency, are formula (10)：
Wherein, R1For the first real vectorial integrated value of frequency；I1For the first empty vectorial integrated value of frequency；L1 is integral and calculating length 1,
Unit dimensionless.
In one embodiment, L1 is 1.5 times of signal unit periodic sequence length.
For step S109, numeral is carried out to the described second real frequency sequence vector and the second empty frequency sequence vector respectively
Trap, the vectorial trap sequence of the real frequency of generation second and the second empty vectorial trap sequence of frequency；
In one embodiment, under the premise of the mixing interfering frequency composition is completely suppressed, the described second real frequency
Vectorial trap sequence and the second empty vectorial trap sequence of frequency are formula (11)：
To R2(n)I2(n) n=0,1,2,3 ..., N-1
To RD2(n)ID2(n) n=0,1,2,3 ..., N-3ND1-3ND2-1
Wherein, RD2(n) it is the described second real vectorial trap sequence of frequency；ID2(n) it is the described second empty vectorial trap sequence of frequency
Row.
For step S110, respectively to the described second real vectorial trap sequence of frequency and the second empty vectorial trap sequence of frequency
Carry out integral operation, the vectorial integrated value of the real frequency of generation second and the second empty vectorial integrated value of frequency；
In one embodiment, the vectorial integrated value of the real frequency of generation second and the second empty vectorial integrated value of frequency, are formula (12)：
Wherein, R2For the second real vectorial integrated value of frequency；I2For the second empty vectorial integrated value of frequency；L2 is integral and calculating length 2,
Unit dimensionless.
In one embodiment, L2 is 1.5 times of signal unit periodic sequence length.
, will the vectorial integrated value of the described first empty frequency and described the according to default phase transition rule for step S111
The one real vectorial integrated value of frequency is converted to first phase, and the described second empty vectorial integrated value of frequency and the described second real frequency vector are integrated
Value is converted to second phase；
In one embodiment, can by below equation (13)-(14) by the described first empty vectorial integrated value of frequency with it is described
The first real vectorial integrated value of frequency is converted to first phase, by the described second empty vectorial integrated value of frequency and the described second real frequency vector product
Score value is converted to second phase：
Wherein, PH1For first phase, unit rad；R1For the first real vectorial integrated value of frequency；I1For the first empty frequency vector integration
Value, PH2For second phase, unit rad；R2For the second real vectorial integrated value of frequency；I2For the second empty vectorial integrated value of frequency.
In one embodiment, according to default phase transition rule, will the described first empty frequently vectorial integrated value with it is described
The first real vectorial integrated value of frequency is converted to first phase, by the described second empty vectorial integrated value of frequency and the described second real frequency vector product
Score value, which is converted to the step of second phase, to be included：
Obtain the first ratio and described the of the described first real vectorial integrated value of frequency and the described first empty vectorial integrated value of frequency
Second ratio of the two real vectorial integrated values of frequency and the described second empty vectorial integrated value of frequency；
The arctan function value of first ratio is obtained, generates the first phase, obtains the anti-of second ratio
Tan value, generate the second phase.
For step S112, according to default sequence average initial phase transformation rule, by the first phase and described
Two phase is converted to the average initial phase of the positive sequence；
In one embodiment, the first phase and the second phase can be converted to by following formula (15)
The average initial phase of the positive sequence：
Wherein, PH0For the average initial phase of positive sequence, unit rad；PH1For first phase, unit rad；PH2For second
Phase, unit rad.
The average initial phase of the positive sequence is the flat of the positive sequence initial phase and the anti-pleat sequence initial phase
Average.
For step S113, according to the average initial phase of the positive sequence, by the initial phase of electric power signal zero just
The amplitude of string FUNCTION MODULATION sequence reverts to the amplitude of the electric power signal, obtains amplitude and recovers sequence；
In one embodiment, can be by following formula (16), by the SIN function of the initial phase of electric power signal zero
The amplitude of modulation sequence reverts to the amplitude of the electric power signal：
In formula, Xout(n) sequence is recovered for the amplitude, N is the positive sequence length, the positive sequence length etc.
In the predetermined sequence length.
For step S114, point on the basis of the central point for recovering sequence by the amplitude, the electric power signal is being obtained just
String function reference signal sequence；
In one embodiment, it is formula (17), formula (18) to obtain SIN function reference signal sequence：
In formula, Xsin+(n) it is SIN function forward reference signal sequence, SIN function forward reference signal sequence initial phase
It is zero；N is the series of discrete value of dimensionless unit, and N is the positive sequence length, and the positive sequence length is equal to described pre-
If sequence length.
In formula, Xsin-(- n) is the reverse reference signal sequence of SIN function, the reverse reference signal sequence first phase of SIN function
Position is zero；N is the series of discrete value of dimensionless unit, and N is the positive sequence length, and the length of the positive sequence is equal to institute
State predetermined sequence length.
The avatars of the SIN function zero initial phase datum mark of reference signal sequence forward or backwards, such as Fig. 3 institutes
Show.
The system for extracting SIN function reference signal in one embodiment from electric power signal, as shown in figure 4, including：
Predetermined sequence length modules 401, for the frequency according to electric power signal, preset sample frequency and default integer signal
Periodicity, obtains the predetermined sequence length of the electric power signal, and the predetermined sequence length is odd number；
Positive block 402, for being sampled according to the predetermined sequence length to the electric power signal, obtain institute
State the positive sequence of electric power signal；
Anti- pleat block 403, for the positive sequence reversely to be exported, obtain corresponding with the positive sequence anti-
Pleat sequence；
SIN function modulation sequence module 404, for the positive sequence and the anti-pleat sequence to be subtracted each other, described in acquisition
The SIN function modulation sequence of the initial phase of electric power signal zero；
Primary vector sequence generating module 405, for respectively by the cosine function of the frequency power signal and the electricity
The SIN function of force signal frequency is multiplied with the positive sequence, generates the first real frequency sequence vector and the first empty vectorial sequence of frequency
Row；
Secondary vector sequence generating module 406, for respectively by the cosine function of the frequency power signal and the electricity
The SIN function of force signal frequency is multiplied with the anti-pleat sequence, generates the second real frequency sequence vector and the second empty vectorial sequence of frequency
Row；
Primary vector trap sequence generating module 407, for respectively to the described first real frequency sequence vector and described first
Empty frequency sequence vector carries out digital notch, the vectorial trap sequence of the real frequency of generation first and the first empty vectorial trap sequence of frequency；
Primary vector integrated value generation module 408, for the vectorial trap sequence of frequency real to described first respectively and described the
The one empty vectorial trap sequence of frequency carries out integral operation, the vectorial integrated value of the real frequency of generation first and the first empty vectorial integrated value of frequency；
Secondary vector trap sequence generating module 409, for respectively to the described second real frequency sequence vector and described second
Empty frequency sequence vector carries out digital notch, the vectorial trap sequence of the real frequency of generation second and the second empty vectorial trap sequence of frequency；
Secondary vector integrated value generation module 410, for the vectorial trap sequence of frequency real to described second respectively and described the
The two empty vectorial trap sequences of frequency carry out integral operation, the vectorial integrated value of the real frequency of generation second and the second empty vectorial integrated value of frequency；
Phase module 411, for according to default phase transition rule, will the described first empty frequently vectorial integrated value with it is described
The first real vectorial integrated value of frequency is converted to first phase, by the described second empty vectorial integrated value of frequency and the described second real frequency vector product
Score value is converted to second phase；
Average initial phase module 412, for according to default sequence average initial phase transformation rule, by the first phase
The average initial phase of the positive sequence is converted to the second phase；
Amplitude recovers block 413, for the average initial phase according to the positive sequence, by the electric power signal zero
The amplitude of the SIN function modulation sequence of initial phase reverts to the amplitude of the electric power signal, obtains amplitude and recovers sequence；
SIN function reference signal module 414, for the point on the basis of the central point of amplitude recovery sequence, obtain institute
State the SIN function reference signal sequence of electric power signal.
Based on the system of the present embodiment shown in Fig. 4, a specific course of work can be discussed further below：
Predetermined sequence length modules 401 are according to the frequency of electric power signal, preset sample frequency and default integer signal first
Periodicity, obtain predetermined sequence length；Positive block 402 samples according to predetermined sequence length to electric power signal, obtains
Obtain positive sequence；Anti- pleat block 403 reversely exports positive sequence, obtains anti-pleat sequence；SIN function modulation sequence mould
Block 404 subtracts each other positive sequence and anti-pleat sequence, obtains the SIN function modulation sequence of zero initial phase；Primary vector sequence generates
The cosine function of frequency power signal and SIN function are multiplied by module 405 with positive sequence respectively, the real frequency vector of generation first
Sequence and the first empty frequency sequence vector；Secondary vector sequence generating module 406 respectively by the cosine function of frequency power signal and
SIN function is multiplied with anti-pleat sequence, generates the second real frequency sequence vector and the second empty frequency sequence vector；Primary vector trap sequence
Column-generation module 407 carries out digital notch to the first real frequency sequence vector and the first empty frequency sequence vector respectively, and generation first is real
Frequently vectorial trap sequence and the first empty vectorial trap sequence of frequency；Primary vector integrated value generation module 408 is respectively to the first real frequency
Vectorial trap sequence and the first empty frequently vectorial trap sequence carry out integral operation, the vectorial integrated value of the real frequency of generation first and the first void
Frequently vectorial integrated value；Secondary vector trap sequence generating module 409 is respectively to the second real frequency sequence vector and the second empty frequency vector
Sequence carries out digital notch, the vectorial trap sequence of the real frequency of generation second and the second empty vectorial trap sequence of frequency；Secondary vector integrates
Being worth generation module 410, frequently vectorial trap sequence real to second and the second empty vectorial trap sequence of frequency carry out integral operation respectively, raw
Into the second real frequently vectorial integrated value and the second empty vectorial integrated value of frequency；Phase module 411, will according to default phase transition rule
The first empty vectorial integrated value of frequency is converted to first phase with the first real vectorial integrated value of frequency, by the second empty vectorial integrated value of frequency and the
The two real vectorial integrated values of frequency are converted to second phase；Average initial phase module 412 is changed according to default sequence average initial phase
Rule, first phase and second phase are converted to the average initial phase of positive sequence；Amplitude recovers block 413 according to just
To the average initial phase of sequence, the amplitude of the SIN function modulation sequence of zero initial phase is reverted to the amplitude of electric power signal, obtained
Obtain amplitude and recover sequence；Point on the basis of the central point that SIN function reference signal module 414 recovers sequence by amplitude, obtain sinusoidal
Function reference signal sequence.
It is evidenced from the above discussion that the present invention can extract SIN function reference signal sequence from Physical Experiment signal, it is determined that
Its zero initial phase datum mark, meets to be actually needed.
Wherein, in one embodiment, the predetermined sequence length modules 401 according to the frequency of electric power signal, default adopt
Sample frequency and default integer signal period number, the predetermined sequence length of electric power signal is obtained, and require default integer signal week
Issue is even number and predetermined sequence length is odd number, and predetermined sequence length correspondingly presets integer signal period number in principle；
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.
The frequency of usual Physical Experiment electric power signal is known, and in one embodiment, the signal frequency is expressed as
Formula (1)：
ω (1)；
In formula, ω is the signal frequency of rad/s units；
In one embodiment, it is formula (2) to obtain predetermined sequence length：
In formula, N is the predetermined sequence length of dimensionless unit；(int) round numbers is represented；C2πTo preset the integer signal period
Number；N2πFor the signal unit periodic sequence length of dimensionless unit；fnFor Hz unit preset sample frequencies；ω is rad/s units
Signal frequency；
In one embodiment, integer signal period number C is preset2π=12, preset sample frequency fn=10KHz.
The positive block 402 samples according to predetermined sequence length to the electric power signal, obtains the electricity
The positive sequence of force signal；
In one embodiment, the electric power signal is the SIN function signal of single fundamental frequency, obtains the positive sequence
It is classified as formula (3)：
Wherein, X+(n) it is positive sequence；A is signal amplitude, unit v；ω is the signal frequency of rad/s units；TnIt is mono- for s
The sampling interval of position；fnFor the sample frequency of Hz units；N is the series of discrete number of immeasurable unit；At the beginning of signal for rad units
Phase；N is the positive sequence length of dimensionless unit, equal to predetermined sequence length.
The avatars of the positive sequence, as shown in Figure 2.
The anti-pleat block 403 reversely exports the positive sequence, obtains corresponding with the positive sequence anti-
Pleat sequence；
In one embodiment, relatively positive sequence, it is formula (4) to obtain anti-pleat sequence：
X-(- n)=X+(N-n)=Asin (- ω Tnn+β)
(4)；
N=0,1,2,3 ..., N-1
In formula, X-(- n) is anti-pleat sequence；β is the anti-pleat sequence initial phase of rad units, in relation, anti-pleat sequence first phase
Position is the cut-off phase of positive sequence；N is the series of discrete number of immeasurable unit；N is the anti-pleat sequence length of dimensionless unit, is equal to
Positive sequence length.
The avatars of the anti-pleat sequence, as shown in Figure 2.
The SIN function modulation sequence module 404 subtracts each other the positive sequence and the anti-pleat sequence, described in acquisition
The SIN function modulation sequence of the initial phase of electric power signal zero, abbreviation modulation sequence；
In one embodiment, the SIN function modulation sequence for obtaining zero initial phase is formula (5)：
In formula, Xsin(n) be zero initial phase SIN function modulation sequence；AsinFor above-mentioned modulation sequence amplitude, unit v；For above-mentioned modulation sequence initial phase, unit rad；For the positive sequence initial phase of rad/s units；β is rad/s
The anti-pleat sequence initial phase of unit.
Allow positive sequence initial phaseExcursion with anti-pleat sequence initial phase β is in 0~± 0.375 π rad.
Error be present because the predetermined sequence length corresponds to integer signal period number, mainly described predetermined sequence length
Integer error.If described error is zero, above-mentioned modulation sequence initial phaseIt is zero, on the contrary initial phaseNear zero.The initial phaseError and the mistake of the integer signal period number compared with null value
It is proportional relation between difference.
The primary vector sequence generating module 405 is respectively by the cosine function of the frequency power signal and the electric power
The SIN function of signal frequency is multiplied with the positive sequence, generates the first real frequency sequence vector and the first empty frequency sequence vector；
In one embodiment, when not considering the mixing interfering frequency composition, it is formula (6) to obtain being mixed sequence：
In formula, R1(n) it is the first real frequency sequence vector；I1(n) it is the first empty frequency sequence vector；N be dimensionless unit just
To sequence length.
The secondary vector sequence generating module 406 is respectively by the cosine function of the frequency power signal and the electric power
The SIN function of signal frequency is multiplied with the anti-pleat sequence, generates the second real frequency sequence vector and the second empty frequency sequence vector；
In one embodiment, when not considering the mixing interfering frequency composition, it is formula (7) to obtain being mixed sequence：
In formula, R2(n) it is the second real frequency sequence vector；I2(n) it is the second empty frequency sequence vector；N is the anti-of dimensionless unit
Pleat sequence length.
The primary vector trap sequence generating module 407 is respectively to the described first real frequency sequence vector and first void
Frequency sequence vector carries out digital notch, the vectorial trap sequence of the real frequency of generation first and the first empty vectorial trap sequence of frequency；
Mixing interfering frequency is included in the real frequency sequence vector and the empty frequency sequence vector.When also existing in input signal
When flip-flop, subharmonic composition and subharmonic composition, the mixing interfering frequency will be more complicated, these mixing interference frequencies
Rate has a strong impact on accuracy in computation.Although window function and integral operation there is good decay to make mixing interfering frequency in itself
With, but there is no specific aim, it is impossible to the inhibitory action of depth is produced to the complicated mixing interfering frequency, it is impossible to meet sinusoidal
The high accuracy of parameter, which calculates, to be needed.
In order to targetedly suppress the influence of the mixing interfering frequency, using a kind of digital trap, ideal situation
Under, the null Frequency point of digital trap just corresponds to the mixing interfering frequency point, has to the mixing interfering frequency
Complete inhibitory action.Preferably, digital notch specifically uses arithmetic average notch algorithm, i.e., by several continuous centrifugal pump phases
Add, then take its arithmetic mean of instantaneous value to be exported as this trap value.Digital notch needs to set digital notch parameter, the numeral
Notch parameter refers to the length N that several continuous centrifugal pumps are addedD.In digital notch parameter NDValue is signal unit periodic sequence
1.5 times of length, it can suppress to being mixed interfering frequency caused by 1/3 subharmonic.And NDValue is signal unit week
2 times of phase sequence length, can to direct current, 1/2 by several times, 1 time, 2 times, 3 times, 4 times, mixing interference frequency caused by 5 subharmonic etc.
Rate is suppressed.Therefore, digital notch is made up of the digital trap of 2 kinds of parameters, it is contemplated that the factor such as physical presence error,
Influenceed in order to which depth suppresses mixing interfering frequency, the digital trap of every kind of parameter is by parameter identical three-level digital notch group
Into totally six grades of arithmetic mean of instantaneous value digital notch are formed.
In one embodiment, six grades of arithmetic mean of instantaneous value digital notch formulas can be formula (8)：
To X (n) n=0,1,2,3 ..., N-1
To XD (n) n=0,1,2,3 ..., N-3ND1-3ND2-1
Wherein, X (n) is digital notch list entries, sequence length N；XD(n) it is digital notch output sequence, sequence length
Spend N-3ND1-3ND2；ND1For notch parameter 1, i.e., continuous centrifugal pump is added quantity；ND2For notch parameter 2, i.e., continuous centrifugal pump phase
Addend amount.
In one embodiment, notch parameter ND1Value is 1.5 times of signal unit periodic sequence length, notch parameter
ND2Value is 2 times of signal unit periodic sequence length, and six grades of arithmetic mean of instantaneous value digital notch need to use 10.5 times of signal lists
Bit period sequence length.
In one embodiment, under the premise of the mixing interfering frequency composition is completely suppressed, the described first real frequency
Vectorial trap sequence and the first empty vectorial trap sequence of frequency are (9)：
To R1(n)I1(n) n=0,1,2,3 ..., N-1
To RD1(n)ID1(n) n=0,1,2,3 ..., N-3ND1-3ND2-1
Wherein, RD1(n) it is the described first real vectorial trap sequence of frequency；ID1(n) it is the described first empty vectorial trap sequence of frequency
Row.
The primary vector integrated value generation module 408 is respectively to the vectorial trap sequence of the described first reality frequency and described first
The empty vectorial trap sequence of frequency carries out integral operation, the vectorial integrated value of the real frequency of generation first and the first empty vectorial integrated value of frequency；
In one embodiment, the vectorial integrated value of the real frequency of generation first and the first empty vectorial integrated value of frequency, are formula (10)：
Wherein, R1For the first real vectorial integrated value of frequency；I1For the first empty vectorial integrated value of frequency；L1 is integral and calculating length 1,
Unit dimensionless.
In one embodiment, L1 is 1.5 times of signal unit periodic sequence length.
The secondary vector trap sequence generating module 409 is respectively to the described second real frequency sequence vector and second void
Frequency sequence vector carries out digital notch, the vectorial trap sequence of the real frequency of generation second and the second empty vectorial trap sequence of frequency；
In one embodiment, under the premise of the mixing interfering frequency composition is completely suppressed, the described second real frequency
Vectorial trap sequence and the second empty vectorial trap sequence of frequency are formula (11)：
To R2(n)I2(n) n=0,1,2,3 ..., N-1
To RD2(n)ID2(n) n=0,1,2,3 ..., N-3ND1-3ND2-1
Wherein, RD2(n) it is the described second real vectorial trap sequence of frequency；ID2(n) it is the described second empty vectorial trap sequence of frequency
Row.
The secondary vector integrated value generation module 410 is respectively to the vectorial trap sequence of the described second reality frequency and described second
The empty vectorial trap sequence of frequency carries out integral operation, the vectorial integrated value of the real frequency of generation second and the second empty vectorial integrated value of frequency；
In one embodiment, the vectorial integrated value of the real frequency of generation second and the second empty vectorial integrated value of frequency, are formula (12)：
Wherein, R2For the second real vectorial integrated value of frequency；I2For the second empty vectorial integrated value of frequency；L2 is integral and calculating length 2,
Unit dimensionless.
In one embodiment, L2 is 1.5 times of signal unit periodic sequence length.
The phase module 411 according to default phase transition rule, will the described first empty frequently vectorial integrated value with it is described
The first real vectorial integrated value of frequency is converted to first phase, by the described second empty vectorial integrated value of frequency and the described second real frequency vector product
Score value is converted to second phase；
In one embodiment, can by below equation (13)-(14) by the described first empty vectorial integrated value of frequency with it is described
The first real vectorial integrated value of frequency is converted to first phase, by the described second empty vectorial integrated value of frequency and the described second real frequency vector product
Score value is converted to second phase：
Wherein, PH1For first phase, unit rad；R1For the first real vectorial integrated value of frequency；I1For the first empty frequency vector integration
Value, PH2For second phase, unit rad；R2For the second real vectorial integrated value of frequency；I2For the second empty vectorial integrated value of frequency.
In one embodiment, according to default phase transition rule, will the described first empty frequently vectorial integrated value with it is described
The first real vectorial integrated value of frequency is converted to first phase, by the described second empty vectorial integrated value of frequency and the described second real frequency vector product
Score value, which is converted to the step of second phase, to be included：
Obtain the first ratio and described the of the described first real vectorial integrated value of frequency and the described first empty vectorial integrated value of frequency
Second ratio of the two real vectorial integrated values of frequency and the described second empty vectorial integrated value of frequency；
The arctan function value of first ratio is obtained, generates the first phase, obtains the anti-of second ratio
Tan value, generate the second phase.
The average initial phase module 412 is according to default sequence average initial phase transformation rule, by the first phase
The average initial phase of the positive sequence is converted to the second phase；
In one embodiment, the first phase and the second phase can be converted to by following formula (15)
The average initial phase of the positive sequence：
Wherein, PH0For the average initial phase of positive sequence, unit rad；PH1For first phase, unit rad；PH2For second
Phase, unit rad.
The average initial phase of the positive sequence is the flat of the positive sequence initial phase and the anti-pleat sequence initial phase
Average.
The amplitude recovers average initial phase of the block 413 according to the positive sequence, by the electric power signal zero
The amplitude of the SIN function modulation sequence of initial phase reverts to the amplitude of the electric power signal, obtains amplitude and recovers sequence；
In one embodiment, can be by following formula (16), by the SIN function of the initial phase of electric power signal zero
The amplitude of modulation sequence reverts to the amplitude of the electric power signal：
In formula, Xout(n) sequence is recovered for the amplitude, N is the positive sequence length, the positive sequence length etc.
In the predetermined sequence length.
Point on the basis of the central point that the SIN function reference signal module 414 recovers sequence by the amplitude, obtain institute
State the SIN function reference signal sequence of electric power signal；
In one embodiment, it is formula (17), formula (18) to obtain SIN function reference signal sequence：
In formula, Xsin+(n) it is SIN function forward reference signal sequence, SIN function forward reference signal sequence initial phase
It is zero；N is the series of discrete value of dimensionless unit, and N is the positive sequence length, and the positive sequence length is equal to described pre-
If sequence length.
In formula, Xsin-(- n) is the reverse reference signal sequence of SIN function, the reverse reference signal sequence first phase of SIN function
Position is zero；N is the series of discrete value of dimensionless unit, and N is the positive sequence length, and the length of the positive sequence is equal to institute
State predetermined sequence length.
The avatars of the SIN function zero initial phase datum mark of reference signal sequence forward or backwards, such as Fig. 3 institutes
Show.
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 (8)
- A kind of 1. method that SIN function reference signal is extracted from electric power signal, it is characterised in that comprise the following steps：According to the frequency of electric power signal, preset sample frequency and default integer signal period number, the pre- of the electric power signal is obtained If sequence length, the predetermined sequence length is odd number；The electric power signal is sampled according to the predetermined sequence length, obtains the positive sequence of the electric power signal；The positive sequence is reversely exported, obtains anti-pleat sequence corresponding with the positive sequence；The positive sequence and the anti-pleat sequence are subtracted each other, obtain the SIN function modulation sequence of the initial phase of electric power signal zero Row；Respectively by the SIN function of the cosine function of the frequency power signal and the frequency power signal and the positive sequence Row are multiplied, and generate the first real frequency sequence vector and the first empty frequency sequence vector；Respectively by the SIN function of the cosine function of the frequency power signal and the frequency power signal and the anti-pleat sequence Row are multiplied, and generate the second real frequency sequence vector and the second empty frequency sequence vector；Digital notch, the real frequency of generation first are carried out to the described first real frequency sequence vector and the first empty frequency sequence vector respectively Vectorial trap sequence and the first empty vectorial trap sequence of frequency；Integral operation, generation are carried out to the described first real vectorial trap sequence of frequency and the first empty vectorial trap sequence of frequency respectively First real frequently vectorial integrated value and the first empty vectorial integrated value of frequency；Digital notch, the real frequency of generation second are carried out to the described second real frequency sequence vector and the second empty frequency sequence vector respectively Vectorial trap sequence and the second empty vectorial trap sequence of frequency；Integral operation, generation are carried out to the described second real vectorial trap sequence of frequency and the second empty vectorial trap sequence of frequency respectively Second real frequently vectorial integrated value and the second empty vectorial integrated value of frequency；According to default phase transition rule, the described first empty vectorial integrated value of frequency and the described first real vectorial integrated value of frequency are turned First phase is changed to, the described second empty vectorial integrated value of frequency and the described second real vectorial integrated value of frequency are converted into second phase；According to default sequence average initial phase transformation rule, by the first phase and the second phase be converted to it is described just To the average initial phase of sequence；Wherein, the average initial phase of the positive sequence is the first phase and the second phase Average value；According to the average initial phase of the positive sequence, by the width of the SIN function modulation sequence of the initial phase of electric power signal zero Value reverts to the amplitude of the electric power signal, obtains amplitude and recovers sequence；Point on the basis of the central point for recovering sequence by the amplitude, obtain the SIN function reference signal sequence of the electric power signal Row；According to default phase transition rule, the described first empty vectorial integrated value of frequency and the described first real vectorial integrated value of frequency are turned First phase is changed to, the described second empty vectorial integrated value of frequency and the described second real vectorial integrated value of frequency are converted into second phase Step includes：Obtain the first ratio of the described first real vectorial integrated value of frequency and the described first empty vectorial integrated value of frequency and described second real Frequently the second ratio of vectorial integrated value and the described second empty vectorial integrated value of frequency；The arctan function value of first ratio is obtained, generates the first phase, obtains the arc tangent of second ratio Functional value, generate the second phase.
- 2. the method according to claim 1 that SIN function reference signal is extracted from electric power signal, it is characterised in that institute The SIN function signal that electric power signal is single fundamental frequency is stated, according to expression formulaObtain it is described just To sequence X+(n), whereinA is the amplitude of the electric power signal, and ω is the electric power The frequency of signal,For the initial phase of the electric power signal, TnFor sampling interval, fnFor the preset sample frequency, N is described Positive sequence length, the positive sequence length are equal to the predetermined sequence length.
- 3. the method according to claim 2 that SIN function reference signal is extracted from electric power signal, it is characterised in that root According to expression formulaObtain the amplitude and recover sequence Xout(n), whereinXsin(n) believe for the electric power The SIN function modulation sequence of number zero initial phase, PH0For the average initial phase of the positive sequence, A is the electric power signal Amplitude,For the positive sequence initial phase, β is the anti-pleat sequence initial phase, and ω is the frequency of the electric power signal,N be the electric power signal series of discrete number, TnFor sampling interval, fnTo be described pre- If sample frequency, N is the positive sequence length, and the positive sequence length is equal to the predetermined sequence length, PH1To be described First phase, PH2For the second phase.
- 4. the method according to claim 2 that SIN function reference signal is extracted from electric power signal, it is characterised in that root According to expression formulaObtain the SIN function forward reference signal sequence X of the electric power signalsin+ (n), according to expression formulaObtain the reverse reference signal sequence of SIN function of the electric power signal Arrange Xsin-(- n), wherein Xout(n) sequence is recovered for the amplitude,N is the electric power signal Series of discrete number, N are the positive sequence length, and the positive sequence length is equal to the predetermined sequence length.
- A kind of 5. system that SIN function reference signal is extracted from electric power signal, it is characterised in that including：Predetermined sequence length modules, for the frequency according to electric power signal, preset sample frequency and default integer signal period number, The predetermined sequence length of the electric power signal is obtained, the predetermined sequence length is odd number；Positive block, for being sampled according to the predetermined sequence length to the electric power signal, obtain the electric power The positive sequence of signal；Anti- pleat block, for the positive sequence reversely to be exported, obtain anti-pleat sequence corresponding with the positive sequence；SIN function modulation sequence module, for the positive sequence and the anti-pleat sequence to be subtracted each other, obtain the electric power letter The SIN function modulation sequence of number zero initial phase；Primary vector sequence generating module, for respectively by the cosine function of the frequency power signal and electric power signal frequency The SIN function of rate is multiplied with the positive sequence, generates the first real frequency sequence vector and the first empty frequency sequence vector；Secondary vector sequence generating module, for respectively by the cosine function of the frequency power signal and electric power signal frequency The SIN function of rate is multiplied with the anti-pleat sequence, generates the second real frequency sequence vector and the second empty frequency sequence vector；Primary vector trap sequence generating module, for respectively to the described first real frequency sequence vector and the first empty frequency vector Sequence carries out digital notch, the vectorial trap sequence of the real frequency of generation first and the first empty vectorial trap sequence of frequency；Primary vector integrated value generation module, for the vectorial trap sequence of frequency real to described first respectively and the first empty frequency to Measure trap sequence and carry out integral operation, the vectorial integrated value of the real frequency of generation first and the first empty vectorial integrated value of frequency；Secondary vector trap sequence generating module, for respectively to the described second real frequency sequence vector and the second empty frequency vector Sequence carries out digital notch, the vectorial trap sequence of the real frequency of generation second and the second empty vectorial trap sequence of frequency；Secondary vector integrated value generation module, for the vectorial trap sequence of frequency real to described second respectively and the second empty frequency to Measure trap sequence and carry out integral operation, the vectorial integrated value of the real frequency of generation second and the second empty vectorial integrated value of frequency；Phase module, it is for regular according to default phase transition, the described first empty vectorial integrated value of frequency and described first is real Frequently vectorial integrated value is converted to first phase, and the described second empty vectorial integrated value of frequency and the described second real vectorial integrated value of frequency are turned It is changed to second phase；Average initial phase module, for according to default sequence average initial phase transformation rule, by the first phase and described Second phase is converted to the average initial phase of the positive sequence；Wherein, the average initial phase of the positive sequence is described the The average value of one phase and the second phase；Amplitude recovers block, for the average initial phase according to the positive sequence, by the initial phase of electric power signal zero The amplitude of SIN function modulation sequence revert to the amplitude of the electric power signal, obtain amplitude and recover sequence；SIN function reference signal module, for the point on the basis of the central point of amplitude recovery sequence, obtain the electric power The SIN function reference signal sequence of signal；The phase module obtains the first ratio of the described first real vectorial integrated value of frequency and the described first empty vectorial integrated value of frequency With the second ratio of the described second real vectorial integrated value of frequency and the described second empty vectorial integrated value of frequency；Obtain first ratio Arctan function value, the first phase is generated, obtain the arctan function value of second ratio, generate second phase Position.
- 6. the system according to claim 5 that SIN function reference signal is extracted from electric power signal, it is characterised in that the electricity Force signal is the SIN function signal of single fundamental frequency, and the positive block is according to expression formula Obtain the positive sequence X+(n), wherein A is the amplitude of the electric power signal, and ω is The frequency of the electric power signal,For the initial phase of the electric power signal, TnFor sampling interval, fnFor the preset sample frequency, N is the positive sequence length, and the positive sequence length is equal to the predetermined sequence length.
- 7. the system according to claim 6 that SIN function reference signal is extracted from electric power signal, it is characterised in that institute State amplitude and recover block according to expression formulaObtain the amplitude and recover sequence Xout(n), whereinXsin(n) believe for the electric power The SIN function modulation sequence of number zero initial phase, PH0For the average initial phase of the positive sequence, A is the electric power signal Amplitude,For the positive sequence initial phase, β is the anti-pleat sequence initial phase, and ω is the frequency of the electric power signal,N be the electric power signal series of discrete number, TnFor sampling interval, fnTo be described pre- If sample frequency, N is the positive sequence length, and the positive sequence length is equal to the predetermined sequence length, PH1To be described First phase, PH2For the second phase.
- 8. the system according to claim 6 that SIN function reference signal is extracted from electric power signal, it is characterised in that institute SIN function reference signal module is stated according to expression formulaObtain the sinusoidal letter of the electric power signal Number forward reference signal sequence Xsin+(n), according to expression formulaObtain the electric power signal The reverse reference signal sequence X of SIN functionsin- (- n), wherein Xout(n) sequence is recovered for the amplitude,N is the series of discrete number of the electric power signal, and N is the positive sequence length, the forward direction Sequence length is equal to the predetermined sequence length.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510599156.0A CN105203842B (en) | 2015-09-18 | 2015-09-18 | The method and system of SIN function reference signal is extracted from electric power signal |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510599156.0A CN105203842B (en) | 2015-09-18 | 2015-09-18 | The method and system of SIN function reference signal is extracted from electric power signal |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105203842A CN105203842A (en) | 2015-12-30 |
| CN105203842B true CN105203842B (en) | 2018-03-20 |
Family
ID=54951625
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510599156.0A Active CN105203842B (en) | 2015-09-18 | 2015-09-18 | The method and system of SIN function reference signal is extracted from electric power signal |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105203842B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH10300801A (en) * | 1997-04-25 | 1998-11-13 | Hitachi Ltd | Phase detecting filter circuit |
| CN1382995A (en) * | 2002-04-24 | 2002-12-04 | 清华大学 | Digital measuring method of frequency and phase |
| CN102135568A (en) * | 2011-02-28 | 2011-07-27 | 无锡辐导微电子有限公司 | Period detection method and phase difference detection method for alternating current signal |
| CN104459320A (en) * | 2014-12-29 | 2015-03-25 | 广东电网有限责任公司电力科学研究院 | Power signal harmonic wave phase measurement method and system |
| CN104655929A (en) * | 2015-01-04 | 2015-05-27 | 中国科学院物理研究所 | Measuring method for digital time frequency of time domain signal and corresponding target identification method |
| CN104914308A (en) * | 2014-12-03 | 2015-09-16 | 许继集团有限公司 | Signal phase measurement method based on two DFT plural spectral lines |
-
2015
- 2015-09-18 CN CN201510599156.0A patent/CN105203842B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH10300801A (en) * | 1997-04-25 | 1998-11-13 | Hitachi Ltd | Phase detecting filter circuit |
| CN1382995A (en) * | 2002-04-24 | 2002-12-04 | 清华大学 | Digital measuring method of frequency and phase |
| CN102135568A (en) * | 2011-02-28 | 2011-07-27 | 无锡辐导微电子有限公司 | Period detection method and phase difference detection method for alternating current signal |
| CN104914308A (en) * | 2014-12-03 | 2015-09-16 | 许继集团有限公司 | Signal phase measurement method based on two DFT plural spectral lines |
| CN104459320A (en) * | 2014-12-29 | 2015-03-25 | 广东电网有限责任公司电力科学研究院 | Power signal harmonic wave phase measurement method and system |
| CN104655929A (en) * | 2015-01-04 | 2015-05-27 | 中国科学院物理研究所 | Measuring method for digital time frequency of time domain signal and corresponding target identification method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105203842A (en) | 2015-12-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105067880B (en) | The method and system of orthogonal modulation is carried out to electric power signal | |
| CN104535836B (en) | Electric signal fundamental frequency measuring method and system | |
| EP1729141B1 (en) | Method for estimating the rotor time constant of an induction machine | |
| CN104330622B (en) | The frequency measurement method and system of sine wave signal in power system | |
| CN104459321B (en) | Power signal base wave phase measurement method and system | |
| CN104459320B (en) | The harmonic phase measuring method and system of electric power signal | |
| CN102508022B (en) | Method for detecting power grid frequency by using optimal multiplier Newton algorithm | |
| CN105203842B (en) | The method and system of SIN function reference signal is extracted from electric power signal | |
| CN105116264B (en) | The method and system of cosine function reference signal is extracted from electric power signal | |
| CN105067885B (en) | Electric power signal is converted to the method and system of zero initial phase signal sequence | |
| CN104330644A (en) | Method and system for detecting amplitude value of sine wave signal in power system | |
| CN102832931A (en) | Phase demodulation method, phase demodulation device and phase-locked loop based on incomplete period grid voltage signal | |
| CN105137172B (en) | Electric power signal is converted to the method and system of zero initial phase cosine signal sequence | |
| CN104991104B (en) | The amplitude detection method and system of electric power signal | |
| CN105158560B (en) | The method and system of orthogonal reference signals sequence is extracted from electric power signal | |
| CN108362940A (en) | A kind of dynamic harmonic frequency extraction method of half cycles distorted signal | |
| CN105092967B (en) | The frequency detecting method and system of electric power signal | |
| CN105203843B (en) | The average first phase position detecting method and system of electric power signal | |
| CN105044462B (en) | The cut-off method for detecting phases and system of electric power signal | |
| CN105319442B (en) | The frequency power signal detection method and system modulated according to cosine function | |
| CN105445551B (en) | The cut-off method for detecting phases and system of sinusoidal signal | |
| CN105425031B (en) | The orthogonal frequency multiplication sequence method of zero initial phase benchmark and system are obtained from electric power signal | |
| CN104977467B (en) | The first phase position detecting method and system of electric power signal | |
| CN105445547A (en) | Phase detection method and system of sinusoidal signals | |
| CN105353217B (en) | The all phase difference detection method and system of sinusoidal signal |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |