CN105301355B - The frequency power signal detection method and system modulated according to SIN function - Google Patents
The frequency power signal detection method and system modulated according to SIN function Download PDFInfo
- Publication number
- CN105301355B CN105301355B CN201510600695.1A CN201510600695A CN105301355B CN 105301355 B CN105301355 B CN 105301355B CN 201510600695 A CN201510600695 A CN 201510600695A CN 105301355 B CN105301355 B CN 105301355B
- Authority
- CN
- China
- Prior art keywords
- frequency
- sequence
- phase
- real
- sin function
- 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
- Inverter Devices (AREA)
Abstract
The present invention relates to a kind of frequency power signal detection method and system modulated according to SIN function.The present invention is modulated by SIN function for the larger signal sequence of initial phase excursion, obtains the SIN function modulation sequence of zero initial phase or initial phase near zero.The SIN function modulation sequence avoids the influence of the larger problem of initial phase excursion, SIN function modulation sequence carries the larger signal sequence all phase difference information of numerical value simultaneously, can significantly improve the degree of accuracy of frequency power signal calculating, improve anti-harmonic wave and noise jamming.The frequency accuracy for the electric power signal that the present invention obtains can reach 10‑10Magnitude, the degree of accuracy that frequency calculates are higher.
Description
Technical field
The present invention relates to technical field of electric power, is examined more particularly to a kind of frequency power signal modulated according to SIN function
Survey method, the frequency power signal detecting system modulated according to SIN function.
Background technology
The measurement of the sine parameter of power system includes frequency measurement, phase measurement, amplitude measurement etc..Fourier transformation is
The basic skills of sine parameter measurement is realized, is widely used in power system.But with sine parameter e measurement technology
Development, the problem of Fourier transformation is present also more aobvious protrusion, it is difficult to further meet power system to sine parameter high precision
Spend the requirement calculated.
In terms of power system sine parameter measurement, there are some improved measurement method of parameters, such as zero hands over method, based on filter
The mensuration of ripple, the mensuration based on Wavelet Transform, based on neutral net, based on DFT (Discrete Fourier
Transform, discrete Fourier transform) conversion mensuration etc..The specified power frequency of operation of power networks belongs near 50Hz (hertz)
The relatively low sinusoidal frequency of frequency.Due to the complexity that the limitation and signal of actual signal treatment technology are formed, such as signal discrete
Data quantization ambient noise caused by sampling influences, and signal sequence blocks caused spectrum leakage problem and is objectively difficult to avoid that,
The influence of any initial phase problem of signal, the influence of direct current and subharmonic and subharmonic problem in signal etc., these algorithms
Measurement accuracy it is relatively low, and anti-harmonic wave and noise jamming are poor.
The content of the invention
Based on this, it is necessary to detected in view of the above-mentioned problems, providing a kind of frequency power signal modulated according to SIN function
Method and system, it is possible to increase the degree of accuracy, the anti-harmonic wave of raising and the noise jamming that frequency calculates.
In order to solve the above technical problems, the present invention adopts the following technical scheme that：
A kind of frequency power signal detection method modulated according to SIN function, including step：
According to the lower limit of frequency power signal scope, preset sample frequency and default integer signal period number, obtain preliminary
Sample sequence length；
The electric power signal is tentatively sampled according to the preliminary sample sequence length, obtains the electric power signal
Preliminary sample sequence；
Line frequency preliminary survey is entered to the preliminary sample sequence, the preliminary frequency of the electric power signal is obtained, according to described first
Synchronizing frequency determines reference frequency；
According to the preset sample frequency and the reference frequency, the unit period sequence length of the electric power signal is obtained
Degree；
According to the default integer signal period number and the unit period sequence length, predetermined sequence length is obtained；
According to the predetermined sequence length, positive sequence is obtained from the preliminary sample sequence；
The positive sequence is reversely exported, obtains the anti-pleat sequence of the positive sequence；
The positive sequence and the anti-pleat sequence are subtracted each other, obtain the SIN function modulation sequence of zero initial phase；
The SIN function modulation sequence is truncated, obtains truncated sequence；
The cosine function of the reference frequency and the SIN function of the reference frequency are adjusted with the SIN function respectively
Sequence processed is multiplied, and obtains the first real frequency sequence vector and the first empty frequency sequence vector；
By the cosine function of the reference frequency and the SIN function of the reference frequency respectively with the truncated sequence phase
Multiply, obtain the second real frequency sequence vector and the second empty frequency sequence vector；
Digital notch is carried out to the described first real frequency sequence vector and the first empty frequency sequence vector respectively, obtains first
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,
Obtain the first real frequently vectorial integrated value and the first empty vectorial integrated value of frequency；
Digital notch is carried out to the described second real frequency sequence vector and the second empty frequency sequence vector respectively, obtains second
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,
Obtain the 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 frequency vector are integrated
Value is converted to first phase；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 cut-off phase transition rule, the first phase and the second phase are converted into the sine
The initial phase of FUNCTION MODULATION sequence；
According to default all phase difference transformation rule, by the initial phase of the SIN function modulation sequence be converted to it is described just
The all phase of string FUNCTION MODULATION sequence is poor；
According to default frequency power signal transformation rule, all phase difference of the SIN function modulation sequence is converted to
The frequency of the electric power signal.
A kind of frequency power signal detecting system modulated according to SIN function, including：
Preliminary sample sequence length determination modul, for the lower limit according to frequency power signal scope, preset sample frequency
With default integer signal period number, preliminary sample sequence length is obtained；
Preliminary sample sequence acquisition module, for being carried out just to the electric power signal according to the preliminary sample sequence length
Step sampling, obtain the preliminary sample sequence of the electric power signal；
Reference frequency determining module, for entering line frequency preliminary survey to the preliminary sample sequence, obtain the electric power signal
Preliminary frequency, reference frequency is determined according to the preliminary frequency；
Unit period sequence length determining module, for according to the preset sample frequency and the reference frequency, obtaining
The unit period sequence length of the electric power signal；
Predetermined sequence length determination modul, for according to the default integer signal period number and the unit period sequence
Length, obtain predetermined sequence length；
Positive retrieval module, for according to the predetermined sequence length, being obtained just from the preliminary sample sequence
To sequence；
Anti- pleat retrieval module, for the positive sequence reversely to be exported, obtain the anti-pleat sequence of the positive sequence
Row；
SIN function modulation sequence determining module, for the positive sequence and the anti-pleat sequence to be subtracted each other, obtain zero
The SIN function modulation sequence of initial phase；
Truncated sequence acquisition module, for the SIN function modulation sequence to be truncated, obtain truncated sequence；
Primary vector sequence determining module, for by the sine of the cosine function of the reference frequency and the reference frequency
Function is multiplied with the SIN function modulation sequence respectively, obtains the first real frequency sequence vector and the first empty frequency sequence vector；
Secondary vector sequence determining module, for by the sine of the cosine function of the reference frequency and the reference frequency
Function is multiplied with the truncated sequence respectively, obtains the second real frequency sequence vector and the second empty frequency sequence vector；
First trap sequence determining module, for respectively to the described first real frequency sequence vector and the first empty frequency vector
Sequence carries out digital notch, obtains the first real frequently vectorial trap sequence and the first empty vectorial trap sequence of frequency；
First integral value determining module, the described first real vectorial trap sequence of frequency and the first empty frequency vector are fallen into respectively
Wave train carries out integral operation, obtains the first real frequently vectorial integrated value and the first empty vectorial integrated value of frequency；
Second trap sequence determining module, for respectively to the described second real frequency sequence vector and the second empty frequency vector
Sequence carries out digital notch, obtains the second real frequently vectorial trap sequence and the second empty vectorial trap sequence of frequency；
Second integral value determining 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, obtain the second real frequently vectorial integrated value and the second empty vectorial integrated value of frequency；
Phase determination module, for regular according to default phase transition, by the described first empty frequently vectorial integrated value and institute
State the first real vectorial integrated value of frequency and be converted to first phase；By the described second empty vectorial integrated value of frequency and the described second real frequency vector
Integrated value is converted to second phase；
Initial phase determining module, for according to default cut-off phase transition rule, by the first phase and described the
Two phase is converted to the initial phase of the SIN function modulation sequence；
All phase difference determining module, for according to default all phase difference transformation rule, the SIN function to be modulated into sequence
The all phase that the initial phase of row is converted to the SIN function modulation sequence is poor；
Frequency detection module, for according to default frequency power signal transformation rule, the SIN function to be modulated into sequence
The all phase difference of row is converted to the frequency of the electric power signal.
The frequency power signal detection method and system that the present invention modulates according to SIN function, for initial phase excursion
Larger signal sequence, is modulated by SIN function, obtains the SIN function modulation sequence of zero initial phase or initial phase near zero
Row.The SIN function modulation sequence avoids the influence of the larger problem of initial phase excursion, while SIN function modulation sequence
Row carry the larger signal sequence all phase difference information of numerical value, can significantly improve frequency power signal calculating the degree of accuracy,
Improve anti-harmonic wave and noise jamming.The frequency accuracy for the electric power signal that the present invention obtains can reach 10-10Magnitude, frequency
The degree of accuracy of calculating is higher.
Brief description of the drawings
Fig. 1 is the schematic flow sheet for the frequency power signal detection method embodiment that the present invention modulates according to SIN function；
Fig. 2 is the schematic diagram of positive sequence of the invention and anti-pleat sequence length；
Fig. 3 is the experimental result schematic diagram that the frequency power signal obtained using the inventive method detects relative error；
Fig. 4 is the structural representation for the frequency power signal detecting system embodiment that the present invention modulates according to SIN function；
Fig. 5 is the structural representation of phase determination module embodiment of the present invention.
Embodiment
In order to which the technical scheme for being better understood from the technical problem to be solved in the present invention, taking and the technology reached are imitated
Fruit, the embodiment of the present invention is described in detail below in conjunction with the accompanying drawings.
As shown in figure 1, a kind of frequency power signal detection method modulated according to SIN function, including step：
S101, according to the lower limit of frequency power signal scope, preset sample frequency and default integer signal period number, obtain
Preliminary sample sequence length；
S102, according to the preliminary sample sequence length electric power signal is tentatively sampled, obtain the electric power
The preliminary sample sequence of signal；
S103, line frequency preliminary survey is entered to the preliminary sample sequence, the preliminary frequency of the electric power signal is obtained, according to institute
State preliminary frequency and determine reference frequency；
S104, according to the preset sample frequency and the reference frequency, obtain the unit period sequence of the electric power signal
Row length；
S105, according to the default integer signal period number and the unit period sequence length, obtain predetermined sequence length
Degree；
S106, according to the predetermined sequence length, obtain positive sequence from the preliminary sample sequence；
S107, the positive sequence reversely exported, obtain the anti-pleat sequence of the positive sequence；
S108, the positive sequence and the anti-pleat sequence subtracted each other, obtain the SIN function modulation sequence of zero initial phase；
S109, the SIN function modulation sequence truncated, obtain truncated sequence；
S110, by the cosine function of the reference frequency and the SIN function of the reference frequency respectively with the sinusoidal letter
Number modulation sequence is multiplied, and obtains the first real frequency sequence vector and the first empty frequency sequence vector；
S111, by the cosine function of the reference frequency and the SIN function of the reference frequency respectively with the truncation sequence
Row are multiplied, and obtain the second real frequency sequence vector and the second empty frequency sequence vector；
S112, digital notch is carried out to the described first real frequency sequence vector and the first empty frequency sequence vector respectively, obtained
To the first real frequently vectorial trap sequence and the first empty vectorial trap sequence of frequency；
S113, the described first real vectorial trap sequence of frequency and the first empty vectorial trap sequence of frequency are integrated respectively
Computing, obtain the first real frequently vectorial integrated value and the first empty vectorial integrated value of frequency；
S114, digital notch is carried out to the described second real frequency sequence vector and the second empty frequency sequence vector respectively, obtained
To the second real frequently vectorial trap sequence and the second empty vectorial trap sequence of frequency；
S115, the described second real vectorial trap sequence of frequency and the second empty vectorial trap sequence of frequency are integrated respectively
Computing, obtain the second real frequently vectorial integrated value and the second empty vectorial integrated value of frequency；
S116, according to default phase transition rule, will the described first empty vectorial integrated value of frequency with described first it is real frequently to
Amount integrated value is converted to first phase；The described second empty vectorial integrated value of frequency and the described second real vectorial integrated value of frequency are converted to
Second phase；
S117, according to it is default cut-off phase transition rule, the first phase and the second phase are converted into institute
State the initial phase of SIN function modulation sequence；
S118, according to default all phase difference transformation rule, the initial phase of the SIN function modulation sequence is converted to
The all phase of the SIN function modulation sequence is poor；
S119, according to default frequency power signal transformation rule, all phase of the SIN function modulation sequence is poor
Be converted to the frequency of the electric power signal.
Electric power signal is the sinusoidal signal based on a kind of fundamental wave component.Sinusoidal signal extensively criticizes string function signal and cosine letter
Number signal.Frequency power signal scope is generally 45Hz~55Hz.So lower limit f of frequency power signal scopeminIt can be taken as
45Hz.Default integer signal period number C2πIt can be configured according to being actually needed, for example, can be by the default integer signal period
Number C2πIt is arranged to 11 etc..In one embodiment, preliminary sample sequence length N can be determined according to formula (1)start：
Wherein, (int) represents to round；fnFor preset sample frequency, unit Hz；NstartUnit dimensionless, C2πList
Position dimensionless, fminUnit be Hz.
After obtaining preliminary sample sequence length, electric power signal is tentatively sampled according to preliminary sample sequence length.Example
Such as, the electric power signal is single fundamental frequency SIN function signal, sinusoidal to single fundamental frequency according to preliminary sample sequence length
Function signal is tentatively sampled, the preliminary sample sequence X of the electric power signal of acquisitionstart(n) it is formula (2)：
Wherein, A is signal amplitude, and unit can be v；ω is signal frequency, unit rad/s；TnFor sampling interval, list
Position is s；fnFor preset sample frequency, unit Hz；N is series of discrete number, unit dimensionless；For signal initial phase, unit is
rad；NstartFor preliminary sample sequence length, unit dimensionless.
After obtaining preliminary sample sequence, method, the algorithm based on filtering can be handed over, based on Wavelet Transformation Algorithm, base by zero
Algorithm in neutral net, the frequency algorithm based on DFT transform or frequency algorithm based on phase difference etc. are to preliminary sample sequence
Enter line frequency preliminary survey, obtain preliminary frequency ωo, preliminary frequency ωoUnit be rad/s.In one embodiment, can should
Preliminary frequency is as reference frequency, i.e. reference frequency ωs=ωo.Reference frequency ωsUnit is rad/s.
After obtaining reference frequency, in one embodiment, unit period sequence length N can be determined according to formula (3)2π：
Wherein, (int) is round numbers；fnFor preset sample frequency, unit Hz；fsFor the reference frequency of Hz units, ωs
For the reference frequency of rad/s units；N2πUnit dimensionless.Unit period sequence length integer was present in 1 sampling interval
Error.
After obtaining unit period sequence length, in one embodiment, predetermined sequence length N is determined according to formula (4)：
N=(int) (C2πN2π) (4)
Wherein, (int) is round numbers；N unit dimensionless.From formula (4) as can be seen that predetermined sequence length N and integer
Signal period number C2πIt is corresponding.Predetermined sequence length N can the periodic sequence length N in units of2π11 times.Due to error be present, institute
It is about to state the signal period integer that predetermined sequence length is included.
After obtaining default series length N, according to predetermined sequence length, positive sequence is obtained from the preliminary sample sequence
Row.In one embodiment, the preliminary sample sequence obtained based on formula (2), the positive sequence X that the present invention obtainsi(n) it is formula
(5)：
Wherein, Xstart(n) it is preliminary sample sequence；A is signal amplitude, and unit can be v；ω is signal frequency, unit
For rad/s；TnFor sampling interval, unit s；N is series of discrete number, unit dimensionless；For signal initial phase, unit is
rad；N is positive sequence length, that is to say predetermined sequence length, unit dimensionless.
The positive sequence obtained based on formula (5), anti-pleat sequence X of the invention-i(- n) is formula (6)：
Wherein, β is anti-pleat sequence initial phase, and anti-pleat sequence initial phase is the cut-off phase of positive sequence, i.e., described electric power
The cut-off phase of signal, unit rad；N is anti-pleat sequence length, as shown in Fig. 2 anti-pleat sequence length and positive sequence length
It is identical, unit dimensionless.
The anti-pleat sequence that positive sequence and formula (6) based on formula (5) obtain, the sinusoidal letter for zero initial phase that the present invention obtains
Number modulation sequence Xcos(n) it is formula (7)：
Wherein, AsinFor the amplitude of SIN function modulation sequence, unit v；For at the beginning of SIN function modulation sequence
Phase, unit rad；N is the series of discrete value of dimensionless unit；N is the SIN function modulation sequence length of dimensionless unit,
SIN function modulation sequence length is identical with positive sequence length.Allow positive sequence initial phaseWith anti-pleat sequence initial phase β
Excursion in 0~± 0.375 π rad.
It is reference frequency error correspondingly to preset integer signal period number to have one of error, reason due to predetermined sequence length
Caused error, the two of reason is the integer error of predetermined sequence length.If described error is zero, SIN function is adjusted
Sequence initial phase processedIt is zero, on the contrary SIN function modulation sequence initial phaseNear zero.The first phase
PositionError compared with null value and be proportional relation between the error of the default integer signal period number.
After obtaining SIN function modulation sequence, the SIN function modulation sequence is truncated, as shown in formula (8), obtained
Truncated sequence is obtained, wherein, the length of truncation can be the unit period sequence length N2π0.25 times：
Wherein, Xs(n) it is truncated sequence；N is the SIN function modulation sequence length of dimensionless unit, that is to say default sequence
Row length；NsFor the truncated sequence length of dimensionless unit；MsFor the shortening value of relative SIN function modulation sequence length, unit
Dimensionless；N2πFor the unit period sequence length of dimensionless unit.
When not considering to be mixed interfering frequency composition, obtained the first real frequency sequence vector R1(n) and the first empty frequency is vectorial
Sequence I1(n) it is formula (9)：
When not considering to be mixed interfering frequency composition, obtained the second real frequency sequence vector R2(n) and the second empty frequency is vectorial
Sequence I2(n) it is formula (10)：
Mixing interfering frequency is included in real frequency sequence vector and empty frequency sequence vector.When in input signal include direct current into
Point, subharmonic composition and during subharmonic composition, the mixing interfering frequency will be more complicated, and these mixing interfering frequencies will be tight
Ghost image rings the degree of accuracy that frequency calculates.Although window function and integral operation have good decay to mixing interfering frequency in itself
Effect, but there is no specific aim, it is impossible to the inhibitory action of depth is produced to the mixing interfering frequency of complexity, it is impossible to meet sinusoidal ginseng
Several high accuracies, which calculates, to be needed.
In order to targetedly suppress to be mixed the influence of interfering frequency, the present invention uses a kind of digital trap.Preferable feelings
Under condition, the null Frequency point of digital trap just corresponds to the mixing interfering frequency point, to the mixing interfering frequency tool
There is complete inhibitory action.In one embodiment, digital notch specifically uses digital averaging filtering algorithm, i.e., connects several
Continuous centrifugal pump is added, and then takes its arithmetic mean of instantaneous value to be exported as this filter value.Digital notch needs to set digital notch to join
Number, the digital notch parameter refer to the length N that several continuous centrifugal pumps are addedD.In digital notch parameter NDValue is unit
At 1.5 times of periodic sequence length, it can suppress to being mixed interfering frequency caused by 1/3 subharmonic.And NDValue is
At 2 times of unit period sequence length, can to direct current, 1/2 by several times, 1 time, 2 times, 3 times, 4 times, mix caused by 5 subharmonic etc.
Frequency interfering frequency is suppressed.Therefore, digital notch is made up of the digital trap of 2 kinds of parameters, it is contemplated that physical presence is missed
The factors such as difference, 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 is formed, and totally six grades of arithmetic mean of instantaneous value digital notch are formed.In one embodiment, six grades of arithmetic mean of instantaneous value numerals
Trap can be formula (11)：
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 the unit period sequence length of the reference frequency,
Notch parameter ND2Value is 2 times of the unit period sequence length of the reference frequency, and six grades of arithmetic mean of instantaneous value digital notch need
To use 10.5 times of unit period sequence lengths.
Under the premise of the mixing interfering frequency composition is completely suppressed, the described first real frequently vectorial trap sequence and institute
It is (12) to state the first empty vectorial trap sequence of frequency：
Wherein, RD1(n) it is the first real vectorial trap sequence of frequency；ID1(n) it is the first empty vectorial trap sequence of frequency；K (Ω) is
Dimensionless gain of the digital notch in frequency difference Ω；α (Ω) is rad unit phase shift of the digital filtering in frequency difference Ω；ND1Join for trap
Number 1；ND2For notch parameter 2.
Equally, under the premise of mixing interfering frequency composition is completely suppressed, the vectorial trap sequence of the second real frequency and
The described second empty vectorial trap sequence of frequency is formula (13)：
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；K (Ω) is dimensionless gain of the digital filtering in frequency difference Ω；α (Ω) is phase shift of the digital filtering in frequency difference Ω, and unit is
rad；ND1For notch parameter 1；ND2For notch parameter 2.
In one embodiment, the present invention carries out integral operation using existing integrator in the prior art.Using integration
Device carries out integral operation to the described first real vectorial trap sequence of frequency and the first empty vectorial trap sequence of frequency, and it is real to obtain first
Frequently vectorial integrated value and the first empty vectorial integrated value of frequency are formula (14)：
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, L1 are 0.5 times of unit period sequence length.
Equally, using integrator respectively to the described second real vectorial trap sequence of frequency and the second empty vectorial trap sequence of frequency
Row carry out integral operation, and it is formula (15) to obtain the second real frequently vectorial integrated value and the second empty vectorial integrated value of frequency：
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, L2 are 0.25 times of unit period sequence length.
In one embodiment, step S120 can include：
Determine the first ratio of the described first real vectorial integrated value of frequency and the described first empty vectorial integrated value of frequency；
Using the arctan function value of first ratio as the first phase；
Determine the second ratio of the described second real vectorial integrated value of frequency and the described second empty vectorial integrated value of frequency；
Using the arctan function value of second ratio as the second phase.
As shown in formula (16) and formula (17), obtained first phase and second phase are respectively：
Wherein, PH1For first phase, unit rad；R1For the first real vectorial integrated value of frequency；I1For the first empty frequency vector product
Score value；PH2For second phase, unit rad；R2For the second real vectorial integrated value of frequency；I2For the second empty vectorial integrated value of frequency；L1
For the integral and calculating length 1 of the dimensionless unit；L2 is the integral and calculating length 2 of the dimensionless unit.
It should be noted that above-mentioned first ratio, the second ratio, just to distinguish each ratio, reduced value is not suitable
Sequence is limited.
After obtaining first phase and second phase, in one embodiment, it can determine that SIN function is adjusted according to formula (18)
The initial phase of sequence processed
The all phase difference of SIN function modulation sequence is signal sinusoidal frequency and SIN function modulation sequence time span
Product.All phase difference △ PH and initial phaseRelation be：Initial phaseSubstantially represent the predetermined sequence length
The error amount of the corresponding integer signal period number, if error amount is zero, all phase difference is 2 π integral multiples, otherwise institute
State the non-2 π integral multiples of all phase difference, initial phaseOpposite number again just reflects the mistake of 2 π integral multiples and non-2 π integral multiples
Difference.
In one embodiment, all phase difference △ PH of SIN function modulation sequence can be determined according to formula (19)：
Wherein, 2 π C2πAll phase for 2 π integral multiples is poor；Work as C2πWhen=11,2 π C2π=22 π, then △ PH scope is in 22 π
Left and right, unit rad.
The relation for being signal sinusoidal frequency with the product of SIN function modulation sequence time span according to all phase difference, can be with
Draw：Frequency power signal is the ratio of all phase difference and SIN function modulation sequence time span.It is so true according to formula (20)
Determine the frequencies omega of electric power signal：
Wherein, ω is the frequency of the electric power signal of rad/s units；TnN is the modulation sequence time span of s units.
In order to be better understood from the technique effect that technical solution of the present invention reaches, carried out with reference to a specific embodiment
Explanation.
Assuming that the experimental signal provided is formula (21)：
For signal fundamental frequency excursion in 45Hz-55Hz, round numbers signal period number is about 11, the change of signal initial phase
The π rad of scope 0~± 0.375, the sample frequency of signal are 10kHz, the discrete data quantization digit 24bit of signal, frequency preliminary survey
Relative error<| ± 0.25% |, obtain frequency power signal detection absolute relative error | Ferr(f) | with signal fundamental frequency f
The experimental result of variation characteristic is as shown in Figure 3.According to Fig. 3 as can be seen that the experimental signal frequency detecting that provides of the present invention is accurate
Degree is 10-10Magnitude, the degree of accuracy are higher.
Based on same inventive concept, the present invention also provides a kind of frequency power signal modulated according to SIN function and detects system
Embodiments of systems of the invention are described in detail by system below in conjunction with the accompanying drawings.
As shown in figure 4, a kind of frequency power signal detecting system modulated according to SIN function, including：
Preliminary sample sequence length determination modul 101, for the lower limit according to frequency power signal scope, default sampling frequency
Rate and default integer signal period number, obtain preliminary sample sequence length；
Preliminary sample sequence acquisition module 102, for being entered according to the preliminary sample sequence length to the electric power signal
The preliminary sampling of row, obtain the preliminary sample sequence of the electric power signal；
Reference frequency determining module 103, for entering line frequency preliminary survey to the preliminary sample sequence, obtain the electric power letter
Number preliminary frequency, reference frequency is determined according to the preliminary frequency；
Unit period sequence length determining module 104, for according to the preset sample frequency and the reference frequency, obtaining
To the unit period sequence length of the electric power signal；
Predetermined sequence length determination modul 105, for according to the default integer signal period number and the unit period
Sequence length, obtain predetermined sequence length；
Positive retrieval module 106, for according to the predetermined sequence length, being obtained from the preliminary sample sequence
Positive sequence；
Anti- pleat retrieval module 107, for the positive sequence reversely to be exported, obtain the anti-pleat of the positive sequence
Sequence；
SIN function modulation sequence determining module 108, for the positive sequence and the anti-pleat sequence to be subtracted each other, obtain
The SIN function modulation sequence of zero initial phase；
Truncated sequence acquisition module 109, for the SIN function modulation sequence to be truncated, obtain truncated sequence；
Primary vector sequence determining module 110, for by the cosine function of the reference frequency and the reference frequency
SIN function is multiplied with the SIN function modulation sequence respectively, obtains the first real frequency sequence vector and the first empty vectorial sequence of frequency
Row；
Secondary vector sequence determining module 111, for by the cosine function of the reference frequency and the reference frequency
SIN function is multiplied with the truncated sequence respectively, obtains the second real frequency sequence vector and the second empty frequency sequence vector；
First trap sequence determining module 112, for respectively to the described first real frequency sequence vector and the first empty frequency
Sequence vector carries out digital notch, obtains the first real frequently vectorial trap sequence and the first empty vectorial trap sequence of frequency；
First integral value determining module 113, respectively the vectorial trap sequence of frequency real to described first and the first empty frequency to
Measure trap sequence and carry out integral operation, obtain the first real frequently vectorial integrated value and the first empty vectorial integrated value of frequency；
Second trap sequence determining module 114, for respectively to the described second real frequency sequence vector and the second empty frequency
Sequence vector carries out digital notch, obtains the second real frequently vectorial trap sequence and the second empty vectorial trap sequence of frequency；
Second integral value determining module 115, for respectively to the described second real vectorial trap sequence of frequency and second void
Frequently vectorial trap sequence carries out integral operation, obtains the second real frequently vectorial integrated value and the second empty vectorial integrated value of frequency；
Phase determination module 116, for according to default phase transition rule, will the vectorial integrated value of the described first empty frequency and
The described 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 to
Amount integrated value is converted to second phase；
Initial phase determining module 117, for according to default cut-off phase transition rule, by the first phase and described
Second phase is converted to the initial phase of the SIN function modulation sequence；
All phase difference determining module 118, for according to default all phase difference transformation rule, the SIN function to be modulated
The all phase that the initial phase of sequence is converted to the SIN function modulation sequence is poor；
Frequency detection module 119, for according to default frequency power signal transformation rule, the SIN function to be modulated
The all phase difference of sequence is converted to the frequency of the electric power signal.
Electric power signal is the sinusoidal signal based on a kind of fundamental wave component.Sinusoidal signal extensively criticizes string function signal and cosine letter
Number signal.Frequency power signal scope is generally 45Hz~55Hz.So lower limit f of frequency power signal scopeminIt can be taken as
45Hz.Default integer signal period number C2πIt can be configured according to being actually needed, for example, can be by the default integer signal period
Number C2πIt is arranged to 11.In one embodiment, the preliminary sample sequence length determination modul 101 can basisIt is determined that preliminary sample sequence length Nstart, wherein (int) represents to round, fnFor default sampling frequency
Rate.
After the preliminary sample sequence length determination modul 101 obtains preliminary sample sequence length, preliminary sample sequence obtains
Modulus block 102 is tentatively sampled according to preliminary sample sequence length to electric power signal.For example, the electric power signal is single
Fundamental frequency SIN function signal, preliminary sample sequence acquisition module 102 is according to preliminary sample sequence length to single fundamental wave frequency
Rate SIN function signal is tentatively sampled, the preliminary sample sequence X of the electric power signal of acquisitionstart(n) it is：Wherein, A is signal amplitude；ω is signal frequency；Between sampling
Every；fnFor preset sample frequency；N=0,1,2,3 ... .., Nstart- 1, it is series of discrete number；For signal initial phase；Nstart
For preliminary sample sequence length.
After preliminary sample sequence acquisition module 102 obtains preliminary sample sequence, reference frequency determining module 103 can pass through
Zero hands over method, the algorithm based on filtering, the algorithm based on Wavelet Transformation Algorithm, based on neutral net, the frequency based on DFT transform to calculate
Method or frequency algorithm based on phase difference etc. enter line frequency preliminary survey to preliminary sample sequence, obtain preliminary frequency ωo.In a reality
Apply in example, reference frequency determining module 103 can be using the preliminary frequency as reference frequency, i.e. reference frequency ωs=ωo。
After reference frequency determining module 103 obtains reference frequency, in one embodiment, unit period sequence length determines
Module 104 can basisDetermine unit period sequence length N2π.Wherein (int) is round numbers, fnFor
Preset sample frequency,fsFor the reference frequency of Hz units, ωsFor the reference frequency of rad/s units.Unit week
The error in 1 sampling interval be present in phase sequence length integer.
After unit period sequence length determining module 104 obtains unit period sequence length, in one embodiment, preset
Sequence length determining module 105 is according to N=(int) (C2πN2π) determine predetermined sequence length N.Wherein (int) is round numbers, C2π
For default integer signal period number, N2πFor unit periodic sequence length.Predetermined sequence length N and integer are can be seen that from the formula
Signal period number C2πIt is corresponding.Predetermined sequence length N can the periodic sequence length N in units of2π11 times.Due to error be present, institute
It is about to state the signal period integer that predetermined sequence length is included.
After predetermined sequence length determination modul 105 obtains default series length N, positive retrieval module 106 is according to pre-
If sequence length, positive sequence is obtained from the preliminary sample sequence.In one embodiment, based on above-mentioned preliminary sampling sequence
The preliminary sample sequence that row acquisition module 102 obtains, the positive sequence X that positive retrieval module 106 obtainsi(n) it is：Wherein, Xstart(n) it is preliminary sample sequence；A is signal amplitude；
ω is signal frequency；TnFor the sampling interval；N=0,1,2,3 ... .., N-1,For series of discrete number；N≤Nstart；For signal
Initial phase；N is positive sequence length, that is to say predetermined sequence length.
The positive sequence obtained based on above-mentioned positive retrieval module 106, anti-pleat retrieval module 107 obtains anti-
Pleat sequence X-i(- n) is：X-i(- n)=Xi(X-n)=Asin(-ωTnn+β).Wherein n=0,1,2,3 ... .., N-1；β is anti-pleat
Sequence initial phase, anti-pleat sequence initial phase are the cut-off phases of positive sequence, i.e., the cut-off phase of described electric power signal；N is anti-
Pleat sequence length, as shown in Fig. 2 anti-pleat sequence length is identical with positive sequence length.
What the positive sequence and anti-pleat retrieval module 107 obtained based on above-mentioned positive retrieval module 106 was obtained
Anti- pleat sequence, the SIN function modulation sequence X for zero initial phase that SIN function modulation sequence determining module 108 obtainscos(n) it is：
Wherein, AsinFor the amplitude of SIN function modulation sequence, unit v；For at the beginning of SIN function modulation sequence
Phase, unit rad；N is the series of discrete value of dimensionless unit；N is the SIN function modulation sequence length of dimensionless unit,
SIN function modulation sequence length is identical with positive sequence length.Allow positive sequence initial phaseWith anti-pleat sequence initial phase β
Excursion in 0~± 0.375 π rad.
It is reference frequency error correspondingly to preset integer signal period number to have one of error, reason due to predetermined sequence length
Caused error, the two of reason is the integer error of predetermined sequence length.If described error is zero, SIN function is adjusted
Sequence initial phase processedIt is zero, on the contrary SIN function modulation sequence initial phaseNear zero.It is described first
PhaseError compared with null value and be proportional relation between the error of the default integer signal period number.
After obtaining SIN function modulation sequence, truncated sequence acquisition module 109 carries out the SIN function modulation sequence
Truncate, be shown below, obtain truncated sequence：
Xs(n)=Xsin(n)
Ns=N-Ms
Ms=0.25N2π
Wherein, Xs(n) it is truncated sequence；N is the SIN function modulation sequence length of dimensionless unit, that is to say default sequence
Row length；NsFor the truncated sequence length of dimensionless unit；MsFor the shortening value of relative SIN function modulation sequence length, unit
Dimensionless；N2πFor the unit period sequence length of dimensionless unit.
When not considering to be mixed interfering frequency composition, the first real frequency that primary vector sequence determining module 110 obtains is vectorial
Sequence R1And the first empty frequency sequence vector I (n)1(n) it is：
When not considering to be mixed interfering frequency composition, the second real frequency that secondary vector sequence determining module 111 obtains is vectorial
Sequence R2And the second empty frequency sequence vector I (n)2(n) it is：
Mixing interfering frequency is included in real frequency sequence vector and empty frequency sequence vector.When in input signal include direct current into
Point, subharmonic composition and during subharmonic composition, the mixing interfering frequency will be more complicated, and these mixing interfering frequencies will be tight
Ghost image rings the degree of accuracy that frequency calculates.Although window function and integral operation have good decay to mixing interfering frequency in itself
Effect, but there is no specific aim, it is impossible to the inhibitory action of depth is produced to the mixing interfering frequency of complexity, it is impossible to meet sinusoidal ginseng
Several high accuracies, which calculates, to be needed.
In order to targetedly suppress to be mixed the influence of interfering frequency, the present invention uses a kind of digital trap.Preferable feelings
Under condition, the null Frequency point of digital trap just corresponds to the mixing interfering frequency point, to the mixing interfering frequency tool
There is complete inhibitory action.In one embodiment, digital notch specifically uses digital averaging filtering algorithm, i.e., connects several
Continuous centrifugal pump is added, and then takes its arithmetic mean of instantaneous value to be exported as this filter value.Digital notch needs to set digital notch to join
Number, the digital notch parameter refer to the length N that several continuous centrifugal pumps are addedD.In digital notch parameter NDValue is unit week
At 1.5 times of phase sequence length, it can suppress to being mixed interfering frequency caused by 1/3 subharmonic.And NDValue is single
At 2 times of bit period sequence length, can to direct current, 1/2 by several times, 1 time, 2 times, 3 times, 4 times, mixing caused by 5 subharmonic etc.
Interfering frequency is suppressed.Therefore, digital notch is made up of the digital trap of 2 kinds of parameters, it is contemplated that physical presence error
Deng factor, influenceed in order to which depth suppresses mixing interfering frequency, the digital trap of every kind of parameter is by the series of parameter identical three
Word trap is formed, and totally six grades of arithmetic mean of instantaneous value digital notch are formed.In one embodiment, six grades of arithmetic mean of instantaneous value numerals are fallen into
Waves are：
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 the unit period sequence length of the reference frequency,
Notch parameter ND2Value is 2 times of the unit period sequence length of the reference frequency, and six grades of arithmetic mean of instantaneous value digital notch need
To use 10.5 times of unit period sequence lengths.
Under the premise of the mixing interfering frequency composition is completely suppressed, the first trap sequence determining module 112 obtains
The described first real vectorial trap sequence of frequency and the first empty vectorial trap sequence of frequency be：
Wherein, RD1(n) it is the first real vectorial trap sequence of frequency；ID1(n) it is the first empty vectorial trap sequence of frequency；K (Ω) is
Dimensionless gain of the digital notch in frequency difference Ω；α (Ω) is rad unit phase shift of the digital filtering in frequency difference Ω；ND1Join for trap
Number 1；ND2For notch parameter 2.
Equally, under the premise of mixing interfering frequency composition is completely suppressed, the second trap sequence determining module 114
To the described second real vectorial trap sequence of frequency and the second empty vectorial trap sequence of frequency be：
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；K (Ω) is dimensionless gain of the digital filtering in frequency difference Ω；α (Ω) is phase shift of the digital filtering in frequency difference Ω, and unit is
rad；ND1For notch parameter 1；ND2For notch parameter 2.
In one embodiment, the present invention carries out integral operation using existing integrator in the prior art.First integral
Value determining module 113 is entered using integrator to the described first real vectorial trap sequence of frequency and the first empty vectorial trap sequence of frequency
Row integral operation, obtains the first real frequently vectorial integrated value and the first empty vectorial integrated value of frequency is：
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, L1 are 0.5 times of unit period sequence length.
Equally, second integral value determining module 115 using integrator respectively the vectorial trap sequence of frequency real to described second and
The described second empty vectorial trap sequence of frequency carries out integral operation, obtains the second real frequently vectorial integrated value and the second empty frequency vector integration
It is worth and is：
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, L2 are 0.25 times of unit period sequence length.
As shown in figure 5, in one embodiment, the phase determination module 116 can include：
First ratio determining unit 1161, for determining the described first real vectorial integrated value of frequency and the described first empty frequency vector
First ratio of integrated value；
First phase determining unit 1162, for using the arctan function value of first ratio as first phase
Position；
Second ratio determining unit 1163, for determining the described second real vectorial integrated value of frequency and the described second empty frequency vector
Second ratio of integrated value；
Second phase determining unit 1164, for using the arctan function value of second ratio as second phase
Position.
It is shown below, the first phase and second phase that the phase determination module 116 obtains are respectively：
Wherein, PH1For first phase, unit rad；R1For the first real vectorial integrated value of frequency；I1For the first empty frequency vector product
Score value；PH2For second phase, unit rad；R2For the second real vectorial integrated value of frequency；I2For the second empty vectorial integrated value of frequency；L1
For the integral and calculating length 1 of the dimensionless unit；L2 is the integral and calculating length 2 of the dimensionless unit.
It should be noted that above-mentioned first ratio, the second ratio, just to distinguish each ratio, reduced value is not suitable
Sequence is limited.
After phase determination module 116 obtains first phase and second phase, in one embodiment, the initial phase determines
The basis of module 117Determine the first phase of the SIN function modulation sequence
PositionWherein, PH1For first phase, PH2For second phase, N is predetermined sequence length, NsFor truncated sequence length.
According to the initial phase of zero phase modulation sequence and default integer signal period number, all phase difference determining module
118 can basisDetermine all phase difference △ of the SIN function modulation sequence
PH, wherein, C2πTo preset integer signal period number,For the initial phase of the SIN function modulation sequence.
In one embodiment, the basis of frequency detection module 119Determine the electric power signal
Frequencies omega, wherein, N is predetermined sequence length,fnFor preset sample frequency.
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 (4)
1. a kind of frequency power signal detection method modulated according to SIN function, it is characterised in that including step：
According to the lower limit of frequency power signal scope, preset sample frequency and default integer signal period number, tentatively sampled
Sequence length；
The electric power signal is tentatively sampled according to the preliminary sample sequence length, obtains the preliminary of the electric power signal
Sample sequence；
Line frequency preliminary survey is entered to the preliminary sample sequence, obtains the preliminary frequency of the electric power signal, according to the just cadence
Rate determines reference frequency；
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, predetermined sequence length is obtained；
According to the predetermined sequence length, positive sequence is obtained from the preliminary sample sequence；
The positive sequence is reversely exported, obtains the anti-pleat sequence of the positive sequence；
The positive sequence and the anti-pleat sequence are subtracted each other, obtain the SIN function modulation sequence of zero initial phase；
The SIN function modulation sequence is truncated, obtains truncated sequence；
The cosine function of the reference frequency and the SIN function of the reference frequency are modulated into sequence with the SIN function respectively
Row are multiplied, and obtain the first real frequency sequence vector and the first empty frequency sequence vector；
The cosine function of the reference frequency and the SIN function of the reference frequency are multiplied with the truncated sequence respectively, obtained
To the second real frequency sequence vector and the second empty frequency sequence vector；
Digital notch is carried out to the described first real frequency sequence vector and the first empty frequency sequence vector respectively, obtains the first real frequency
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, obtained
First real frequently vectorial integrated value and the first empty vectorial integrated value of frequency；
Digital notch is carried out to the described second real frequency sequence vector and the second empty frequency sequence vector respectively, obtains the second real frequency
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, obtained
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
It is changed to first phase；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 cut-off phase transition rule, the first phase and the second phase are converted into the SIN function
The initial phase of modulation sequence；According toDetermine the first phase of the SIN function modulation sequence
PositionWherein, PH1For first phase, PH2For second phase, N is predetermined sequence length, NsFor truncated sequence length
According to default all phase difference transformation rule, the initial phase of the SIN function modulation sequence is converted into the sinusoidal letter
The all phase of number modulation sequence is poor；According toDetermine the SIN function modulation sequence
All phase difference △ PH, wherein, C2πTo preset integer signal period number,For the first phase of the SIN function modulation sequence
Position；
According to default frequency power signal transformation rule, all phase difference of the SIN function modulation sequence is converted to described
The frequency of electric power signal；Wherein, the frequency power signal is the ratio of all phase difference and SIN function modulation sequence time span
Value；
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
It is changed to first phase；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：
Determine the first ratio of the described first real vectorial integrated value of frequency and the described first empty vectorial integrated value of frequency；
Using the arctan function value of first ratio as the first phase；
Determine the second ratio of the described second real vectorial integrated value of frequency and the described second empty vectorial integrated value of frequency；
Using the arctan function value of second ratio as the second phase.
2. the frequency power signal detection method according to claim 1 modulated according to SIN function, it is characterised in that root
According toThe frequencies omega of the electric power signal is determined, wherein, N is predetermined sequence length,fnIt is default
Sample frequency.
A kind of 3. frequency power signal detecting system modulated according to SIN function, it is characterised in that including：
Preliminary sample sequence length determination modul, for the lower limit according to frequency power signal scope, preset sample frequency and pre-
If integer signal period number, preliminary sample sequence length is obtained；
Preliminary sample sequence acquisition module, for tentatively being adopted to the electric power signal according to the preliminary sample sequence length
Sample, obtain the preliminary sample sequence of the electric power signal；
Reference frequency determining module, for entering line frequency preliminary survey to the preliminary sample sequence, obtain the first of the electric power signal
Synchronizing frequency, reference frequency is determined according to the preliminary frequency；
Unit period sequence length determining module, for according to the preset sample frequency and the reference frequency, obtaining described
The unit period sequence length of electric power signal；
Predetermined sequence length determination modul, for according to the default integer signal period number and unit period sequence length
Degree, obtains predetermined sequence length；
Positive retrieval module, for according to the predetermined sequence length, positive sequence to be obtained from the preliminary sample sequence
Row；
Anti- pleat retrieval module, for the positive sequence reversely to be exported, obtain the anti-pleat sequence of the positive sequence；
SIN function modulation sequence determining module, for the positive sequence and the anti-pleat sequence to be subtracted each other, obtain zero first phase
The SIN function modulation sequence of position；
Truncated sequence acquisition module, for the SIN function modulation sequence to be truncated, obtain truncated sequence；
Primary vector sequence determining module, for by the SIN function of the cosine function of the reference frequency and the reference frequency
It is multiplied respectively with the SIN function modulation sequence, obtains the first real frequency sequence vector and the first empty frequency sequence vector；
Secondary vector sequence determining module, for by the SIN function of the cosine function of the reference frequency and the reference frequency
It is multiplied respectively with the truncated sequence, obtains the second real frequency sequence vector and the second empty frequency sequence vector；
First trap sequence determining module, for respectively to the described first real frequency sequence vector and the first empty frequency sequence vector
Digital notch is carried out, obtains the first real frequently vectorial trap sequence and the first empty vectorial trap sequence of frequency；
First integral value determining module, respectively to the described first real vectorial trap sequence of frequency and the first empty vectorial trap sequence of frequency
Row carry out integral operation, obtain the first real frequently vectorial integrated value and the first empty vectorial integrated value of frequency；
Second trap sequence determining module, for respectively to the described second real frequency sequence vector and the second empty frequency sequence vector
Digital notch is carried out, obtains the second real frequently vectorial trap sequence and the second empty vectorial trap sequence of frequency；
Second integral value determining module, for being fallen into respectively to the described second real vectorial trap sequence of frequency and the second empty frequency vector
Wave train carries out integral operation, obtains the second real frequently vectorial integrated value and the second empty vectorial integrated value of frequency；
Phase determination 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；The described second empty vectorial integrated value of frequency and the described second real frequency vector are integrated
Value is converted to second phase；
Initial phase determining module, for regular according to default cut-off phase transition, by the first phase and second phase
Position is converted to the initial phase of the SIN function modulation sequence；
All phase difference determining module, for according to default all phase difference transformation rule, by the SIN function modulation sequence
The all phase that initial phase is converted to the SIN function modulation sequence is poor；
Frequency detection module, for according to default frequency power signal transformation rule, by the SIN function modulation sequence
All phase difference is converted to the frequency of the electric power signal；Wherein, the frequency power signal is that all phase difference is adjusted with SIN function
The ratio of sequence time length processed；
The all phase difference determining module according toDetermine the SIN function modulation sequence
The all phase difference △ PH of row, wherein, C2πTo preset integer signal period number,For the first phase of the SIN function modulation sequence
Position；
The phase determination module includes：
First ratio determining unit, for determining the described first real vectorial integrated value of frequency and the described first empty vectorial integrated value of frequency
First ratio；
First phase determining unit, for using the arctan function value of first ratio as the first phase；
Second ratio determining unit, for determining the described second real vectorial integrated value of frequency and the described second empty vectorial integrated value of frequency
Second ratio；
Second phase determining unit, for using the arctan function value of second ratio as the second phase；
The initial phase determining module according toDetermine the SIN function modulation sequence
Initial phaseWherein, PH1For first phase, PH2For second phase, N is predetermined sequence length, NsFor truncated sequence length.
4. the frequency power signal detecting system according to claim 3 modulated according to SIN function, it is characterised in that institute
State frequency detection module according toThe frequencies omega of the electric power signal is determined, wherein, N is predetermined sequence length,fnFor preset sample frequency.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510600695.1A CN105301355B (en) | 2015-09-18 | 2015-09-18 | The frequency power signal detection method and system modulated according to SIN function |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510600695.1A CN105301355B (en) | 2015-09-18 | 2015-09-18 | The frequency power signal detection method and system modulated according to SIN function |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105301355A CN105301355A (en) | 2016-02-03 |
| CN105301355B true CN105301355B (en) | 2018-04-03 |
Family
ID=55198853
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510600695.1A Active CN105301355B (en) | 2015-09-18 | 2015-09-18 | The frequency power signal detection method and system modulated according to SIN function |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105301355B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004515141A (en) * | 2000-11-21 | 2004-05-20 | ユビネティクス リミテッド | Method and apparatus for estimating signal phase |
| CN104502700A (en) * | 2014-12-29 | 2015-04-08 | 广东电网有限责任公司电力科学研究院 | Sinusoidal parameter measurement method and system of power signal |
| CN104502698A (en) * | 2014-12-10 | 2015-04-08 | 广东电网有限责任公司电力科学研究院 | Method and system for measuring frequency of electric power signal |
| CN104502701A (en) * | 2014-12-10 | 2015-04-08 | 广东电网有限责任公司电力科学研究院 | Method and system for detecting electric power signal frequency based on phase modulation |
| CN104502706A (en) * | 2014-12-29 | 2015-04-08 | 广东电网有限责任公司电力科学研究院 | Harmonic amplitude measurement method and system of power signal |
-
2015
- 2015-09-18 CN CN201510600695.1A patent/CN105301355B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004515141A (en) * | 2000-11-21 | 2004-05-20 | ユビネティクス リミテッド | Method and apparatus for estimating signal phase |
| CN104502698A (en) * | 2014-12-10 | 2015-04-08 | 广东电网有限责任公司电力科学研究院 | Method and system for measuring frequency of electric power signal |
| CN104502701A (en) * | 2014-12-10 | 2015-04-08 | 广东电网有限责任公司电力科学研究院 | Method and system for detecting electric power signal frequency based on phase modulation |
| CN104502700A (en) * | 2014-12-29 | 2015-04-08 | 广东电网有限责任公司电力科学研究院 | Sinusoidal parameter measurement method and system of power signal |
| CN104502706A (en) * | 2014-12-29 | 2015-04-08 | 广东电网有限责任公司电力科学研究院 | Harmonic amplitude measurement method and system of power signal |
Non-Patent Citations (1)
| Title |
|---|
| Mean square error analysis of unbiased modified plain gradient algorithm for second-order adaptive IIR notch filter;R.Punchalard;《SignalProcessing》;20120506;第2815-2820页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105301355A (en) | 2016-02-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104635045B (en) | Frequency power signal detection method and system based on phase-modulation | |
| CN104635044B (en) | Power signal frequency detection method and system based on amplitude modulation | |
| CN104502700B (en) | Sinusoidal parameter measurement method and system of power signal | |
| CN104502706B (en) | Harmonic amplitude measurement method and system of power signal | |
| CN104535836B (en) | Electric signal fundamental frequency measuring method and system | |
| CN104459321B (en) | Power signal base wave phase measurement method and system | |
| CN104502675B (en) | Fundamental wave amplitude method and system of power signal | |
| CN105301355B (en) | The frequency power signal detection method and system modulated according to SIN function | |
| CN105158564B (en) | The electric power signal all phase difference detection method and system modulated according to SIN function | |
| CN105319442B (en) | The frequency power signal detection method and system modulated according to cosine function | |
| CN105203844B (en) | The zero initial phase modulator approach and system of electric power signal | |
| CN105067885B (en) | Electric power signal is converted to the method and system of zero initial phase signal sequence | |
| CN105182077B (en) | The electric power signal all phase difference detection method and system modulated according to cosine function | |
| CN105137172B (en) | Electric power signal is converted to the method and system of zero initial phase cosine signal sequence | |
| CN105445551B (en) | The cut-off method for detecting phases and system of sinusoidal signal | |
| CN105203843B (en) | The average first phase position detecting method and system of electric power signal | |
| CN105137188B (en) | The method and system of quadrature signal sequence is extracted from electric power signal | |
| CN105425031B (en) | The orthogonal frequency multiplication sequence method of zero initial phase benchmark and system are obtained from electric power signal | |
| CN105044462B (en) | The cut-off method for detecting phases and system of electric power signal | |
| CN105092967B (en) | The frequency detecting method and system of electric power signal | |
| CN104991104A (en) | Power signal amplitude detection method and system | |
| CN105425040B (en) | Based on zero initial phase benchmark cosine function frequency multiplication sequence frequency measuring method and system | |
| CN105548692B (en) | The method and system of frequency measurement is carried out according to zero initial phase benchmark orthogonal sequence | |
| CN105353217B (en) | The all phase difference detection method and system of sinusoidal signal | |
| CN105158558B (en) | The frequency detecting method and system of electric power signal |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |