CN105067880B - The method and system of orthogonal modulation is carried out to electric power signal - Google Patents

The method and system of orthogonal modulation is carried out to electric power signal Download PDF

Info

Publication number
CN105067880B
CN105067880B CN201510600888.7A CN201510600888A CN105067880B CN 105067880 B CN105067880 B CN 105067880B CN 201510600888 A CN201510600888 A CN 201510600888A CN 105067880 B CN105067880 B CN 105067880B
Authority
CN
China
Prior art keywords
sequence
frequency
phase
electric power
preliminary
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510600888.7A
Other languages
Chinese (zh)
Other versions
CN105067880A (en
Inventor
李军
任娟娟
欧阳春明
陈华忠
湖康涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Electric Power Research Institute of Guangdong Power Grid Co Ltd
Original Assignee
Electric Power Research Institute of Guangdong Power Grid Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Electric Power Research Institute of Guangdong Power Grid Co Ltd filed Critical Electric Power Research Institute of Guangdong Power Grid Co Ltd
Priority to CN201510600888.7A priority Critical patent/CN105067880B/en
Publication of CN105067880A publication Critical patent/CN105067880A/en
Application granted granted Critical
Publication of CN105067880B publication Critical patent/CN105067880B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The invention discloses a kind of method and system that orthogonal modulation is carried out to electric power signal, methods described includes：According to frequency power signal lower limit, preset sample frequency and preliminary design integer signal period number, preliminary sequence length, preliminary sequence are obtained；Reference frequency is obtained according to preliminary sequence, further obtains unit period sequence length, predetermined sequence length；According to predetermined sequence starting point and predetermined sequence length, the first positive, anti-pleat sequence is obtained；According to first positive, anti-the first real frequency of pleat sequence generation, empty frequency sequence vector and the second real frequency, empty frequency sequence vector；It is again the first positive sequence average initial phase by first, second obtained phase transition；New predetermined sequence starting point is obtained according to average initial phase；Further according to the second positive, the anti-pleat sequence obtained from preliminary sequence, cosine function modulation sequence, SIN function modulation sequence are obtained.The present invention obtains the consistent orthogonal modulation sequences of amplitude height, improves the degree of accuracy and anti-interference that sine parameter calculates.

Description

The method and system of orthogonal modulation is carried out to electric power signal
Technical field
The present invention relates to technical field of electric power, more particularly to a kind of method to electric power signal progress orthogonal modulation and is System.
Background technology
The frequency measurement of power system, phase measurement, amplitude measurement etc. are the measurement of sine parameter in itself.Fu Li Leaf transformation is the basic skills for realizing sine parameter measurement, is widely used in power system.But as sine parameter is surveyed The development of amount technology, the problem of Fourier transformation is present more aobvious protrusion, it is difficult to further meet power system to sine parameter The requirement that high accuracy calculates.
In terms of power system sine parameter measurement, there is many measuring methods, such as zero hand over method, the mensuration based on filtering, Mensuration based on Wavelet Transform, based on neutral net, mensuration based on DFT transform etc..But the specified work of operation of power networks Frequency belongs to the relatively low sinusoidal frequency of frequency near 50Hz.Formed due to the limitation and signal of actual signal treatment technology Complexity, above-mentioned measuring method are low to sinusoidal parameter measurement precision, and anti-harmonic wave and noise jamming are poor.
The content of the invention
Based on the above situation, the present invention proposes a kind of method and system that orthogonal modulation is carried out to electric power signal, obtains Orthogonal modulation sequences, improve the degree of accuracy and anti-harmonic wave, noise jamming that sine parameter calculates.
To achieve these goals, the embodiment of technical solution of the present invention is：
A kind of method that orthogonal modulation is carried out to electric power signal, comprises the following steps：
According to the lower limit of frequency power signal scope, preset sample frequency and preliminary design integer signal period number, described in acquisition The preliminary sequence length of electric power signal；
The electric power signal is tentatively sampled according to the preliminary sequence length, obtains the preliminary of the electric power signal Sequence；
Line frequency preliminary survey is entered to the preliminary sequence, generates the preliminary frequency of the electric power signal, and according to described preliminary The reference frequency of electric power signal described in frequency setting；
According to the preset sample frequency and the reference frequency, the unit period sequence for obtaining the electric power signal is grown Degree；
The default integer signal period number of the electric power signal is obtained according to the preliminary design integer signal period number, and according to The default integer signal period number and the unit period sequence length, obtain the predetermined sequence length of the electric power signal；
According to predetermined sequence starting point and the predetermined sequence length, the electric power signal is obtained from the preliminary sequence The first positive sequence；
Described first positive sequence is reversely exported, obtains the first anti-pleat sequence corresponding with the described first positive sequence；
Respectively by the cosine function of the reference frequency and the SIN function of the reference frequency and the described first positive sequence Row are multiplied, and generate the first real frequency sequence vector and the first empty frequency sequence vector；
Respectively by the cosine function of the reference frequency and the SIN function of the reference frequency and the described first anti-pleat sequence Row are multiplied, and generate the second real frequency sequence vector and the second empty frequency sequence vector；
Integral and calculating, generation first are carried out to the described first real frequency sequence vector and the first empty frequency sequence vector respectively Real frequently vectorial integrated value and the first empty vectorial integrated value of frequency；
Integral and calculating, generation second are carried out to the described second real frequency sequence vector and the second empty frequency sequence vector respectively 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, 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 the first positive sequence；
Compared with the phase of the average initial phase and ± π/4, the phase bit comparison of the phase bit comparison of acquisition and the ± π/4 Difference；
Difference, the predetermined sequence starting point and the unit period sequence length are compared according to the phase, obtained new Predetermined sequence starting point；
According to the new predetermined sequence starting point and the predetermined sequence length, obtained from the preliminary sequence described Second positive sequence of electric power signal；
Described second positive sequence is reversely exported, obtains the second anti-pleat sequence corresponding with the described second positive sequence；
Described second positive sequence is added with the described second anti-pleat sequence, the cosine function for obtaining the electric power signal is adjusted Sequence processed, the described second positive sequence and the second anti-pleat sequence are subtracted each other, the SIN function for obtaining the electric power signal is adjusted Sequence processed.
A kind of system that orthogonal modulation is carried out to electric power signal, including：
Preliminary sequence length modules, it is whole for the lower limit according to frequency power signal scope, preset sample frequency and preliminary design Number signal period number, obtain the preliminary sequence length of the electric power signal；
Preliminary sequence module, for tentatively being sampled to the electric power signal according to the preliminary sequence length, obtain The preliminary sequence of the electric power signal；
Frequency preliminary survey module, for entering line frequency preliminary survey to the preliminary sequence, generate the first cadence of the electric power signal Rate, and according to the reference frequency of the preliminary frequency setting electric power signal；
Unit period sequence length module, for according to the preset sample frequency and the reference frequency, described in acquisition The unit period sequence length of electric power signal；
Predetermined sequence length modules, for obtaining the default of the electric power signal according to the preliminary design integer signal period number Integer signal period number, and according to the default integer signal period number and the unit period sequence length, obtain the electricity The predetermined sequence length of force signal；
First positive block, for according to predetermined sequence starting point and the predetermined sequence length, from described preliminary The first positive sequence of the electric power signal is obtained in sequence；
First anti-pleat block, for the described first positive sequence reversely to be exported, obtain and the described first positive sequence First anti-pleat sequence corresponding to row；
Primary vector sequence generating module, for respectively by the cosine function of the reference frequency and the reference frequency SIN function is multiplied with the described first 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 reference frequency and the reference frequency SIN function is multiplied with the described first anti-pleat sequence, generates the second real frequency sequence vector and the second empty frequency sequence vector；
Primary vector integrated value generation module, for respectively to the described first real frequency sequence vector and the first empty frequency to Measure sequence and carry out integral and calculating, the vectorial integrated value of the real frequency of generation first and the first empty vectorial integrated value of frequency；
Secondary vector integrated value generation module, for respectively to the described second real frequency sequence vector and the second empty frequency to Measure sequence and carry out integral and calculating, 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；
Phase detecting module, for according to default sequence average initial phase transformation rule, by the first phase and institute State the average initial phase that second phase is converted to the described first positive sequence；
Phase compares difference block, compared with by the phase of the average initial phase and ± π/4, obtain and described ± The phase of the phase bit comparison of π/4 compares difference；
Predetermined sequence originates point module, for comparing difference, the predetermined sequence starting point and described according to the phase Unit period sequence length, obtain new predetermined sequence starting point；
Second positive block, for according to the new predetermined sequence starting point and the predetermined sequence length, from The second positive sequence of the electric power signal is obtained in the preliminary sequence；
Second anti-pleat block, for the described second positive sequence reversely to be exported, obtain and the described second positive sequence Second anti-pleat sequence corresponding to row；
Orthogonal modulation module, for the described second positive sequence to be added with the described second anti-pleat sequence, obtain the electricity The cosine function modulation sequence of force signal, the described second positive sequence and the second anti-pleat sequence are subtracted each other, obtain the electricity The SIN function modulation sequence of force signal.
Compared with prior art, beneficial effects of the present invention are：The method that the present invention carries out orthogonal modulation to electric power signal And system, according to the lower limit of frequency power signal scope, preset sample frequency and preliminary design integer signal period number, obtain preliminary sequence Row length, and electric power signal is tentatively sampled, obtain preliminary sequence；Enter line frequency preliminary survey to preliminary sequence, generation is preliminary Frequency, set reference frequency；According to preset sample frequency and reference frequency, unit period sequence length is obtained；According to default whole Number signal period number and unit period sequence length, obtain predetermined sequence length；The first positive sequence is obtained from preliminary sequence, Further obtain the first anti-pleat sequence；The cosine function of reference frequency and SIN function are multiplied with the first positive sequence respectively, The first real frequency sequence vector and the first empty frequency sequence vector are generated, then is multiplied respectively with the first anti-pleat sequence, the real frequency of generation second Sequence vector and the second empty frequency sequence vector；By carrying out integral and calculating to the sequence vector of above-mentioned generation respectively, and then obtain The first phase and second phase, phase-detection is carried out, obtains new predetermined sequence starting point, the is obtained from preliminary sequence Two positive sequences, the second anti-pleat sequence is further obtained, carry out electric power signal orthogonal modulation.Obtain basically identical orthogonal of amplitude Modulation sequence, improve the degree of accuracy and anti-harmonic wave, noise jamming that sine parameter calculates.
Brief description of the drawings
Fig. 1 is the method flow diagram for carrying out orthogonal modulation in one embodiment to electric power signal；
Fig. 2 is the first positive sequence and the first anti-pleat sequence diagram in one embodiment；
Fig. 3 is the system structure diagram for carrying out orthogonal modulation in one embodiment to 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 carrying out orthogonal modulation in one embodiment to electric power signal, as shown in figure 1, comprising the following steps：
Step S101：According to the lower limit of frequency power signal scope, preset sample frequency and preliminary design integer signal period number, Obtain the preliminary sequence length of the electric power signal；
Step S102：The electric power signal is tentatively sampled according to the preliminary sequence length, obtains the electric power The preliminary sequence of signal；
Step S103：Line frequency preliminary survey is entered to the preliminary sequence, generates the preliminary frequency of the electric power signal, and according to The preliminary frequency sets the reference frequency of the electric power signal；
Step S104：According to the preset sample frequency and the reference frequency, the unit for obtaining the electric power signal is all Phase sequence length；
Step S105：The default integer signal period of the electric power signal is obtained according to the preliminary design integer signal period number Number, and according to the default integer signal period number and the unit period sequence length, obtain the default of the electric power signal Sequence length；
Step S106：According to predetermined sequence starting point and the predetermined sequence length, institute is obtained from the preliminary sequence State the first positive sequence of electric power signal；
Step S107：Described first positive sequence is reversely exported, obtained and the described first positive sequence corresponding first Anti- pleat sequence；
Step S108：Respectively by the cosine function of the reference frequency and the SIN function of the reference frequency and described the One positive sequence is multiplied, and generates the first real frequency sequence vector and the first empty frequency sequence vector；
Step S109：Respectively by the cosine function of the reference frequency and the SIN function of the reference frequency and described the One anti-pleat sequence is multiplied, and generates the second real frequency sequence vector and the second empty frequency sequence vector；
Step S110：Integrating meter is carried out to the described first real frequency sequence vector and the first empty frequency sequence vector respectively Calculate, the vectorial integrated value of the real frequency of generation first and the first empty vectorial integrated value of frequency；
Step S111：Integrating meter is carried out to the described second real frequency sequence vector and the second empty frequency sequence vector respectively Calculate, the vectorial integrated value of the real frequency of generation second and the second empty vectorial integrated value of frequency；
Step S112：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 S113：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 described first positive sequence；
Step S114：Compared with the phase of the average initial phase and ± π/4, the phase bit comparison of acquisition and the ± π/4 Phase compare difference；
Step S115：Difference, the predetermined sequence starting point and unit period sequence length are compared according to the phase Degree, obtains new predetermined sequence starting point；
Step S116：According to the new predetermined sequence starting point and the predetermined sequence length, from the preliminary sequence Middle the second positive sequence for obtaining the electric power signal；
Step S117：Described second positive sequence is reversely exported, obtained and the described second positive sequence corresponding second Anti- pleat sequence；
Step S118：Described second positive sequence is added with the described second anti-pleat sequence, obtains the electric power signal Cosine function modulation sequence, the described second positive sequence and the second anti-pleat sequence are subtracted each other, obtain the electric power signal SIN function modulation sequence.
It is evidenced from the above discussion that the present invention obtains the basically identical orthogonal modulation sequences of amplitude, have excellent calculating special Property, improve the degree of accuracy and anti-harmonic wave, noise jamming that sine parameter calculates.
Wherein, for step S101, believed according to the lower limit of frequency power signal scope, preset sample frequency and preliminary design integer Number periodicity, obtain the preliminary sequence length of the electric power signal；
The electric power signal is the sinusoidal signal based on a kind of fundamental wave component.Sinusoidal signal extensively criticizes string function signal and remaining String function signal.
In one embodiment, power system frequency scope is in 45Hz-55Hz, power taking force signal lower-frequency limit fminFor 45Hz；According to being actually needed, the preliminary design integer signal period number C is set2π-start, it is preferable that take C2π-startFor integer 13.
In one embodiment, it is formula (1) to obtain the preliminary sequence length：
In formula, NstartFor the preliminary sequence length, unit dimensionless；(int) round numbers is represented；C2π-startTo be described Preliminary design integer signal period number, unit dimensionless；fminFor the lower limit of the signal frequency range, unit Hz；fnPreset to be described Sample frequency, unit Hz.
For step S102, the electric power signal is tentatively sampled according to the preliminary sequence length, described in acquisition The preliminary sequence of electric power signal.In one embodiment, the electric power signal is the cosine function signal of single fundamental frequency, is obtained The preliminary sequence of the electric power signal is formula (2)：
Wherein, Xstart(n) it is preliminary sequence；A be the electric power signal amplitude, unit v；ω is the electric power signal Frequency, unit rad/s；TnFor sampling interval, unit s；fnFor the preset sample frequency, unit Hz；N is the electric power signal Series of discrete number, unit dimensionless；For signal (preliminary sequence) initial phase, unit rad；NstartFor preliminary sequence length, Unit dimensionless.
For step S103, can by zero hand over method, the algorithm based on filtering, based on Wavelet Transformation Algorithm, based on nerve net The algorithm of network, the frequency algorithm based on DFT transform or at the beginning of entering line frequency to the preliminary sequence based on the frequency algorithm of phase difference Survey, generate the preliminary frequency of the electric power signal.
In one embodiment, it is formula (3) to generate the preliminary frequency：
ωo(3)；
Wherein, ωoFor the preliminary frequency, unit rad/s；
Preferably, the reference frequency is formula (4) equal to the preliminary frequency：
ωs=ωo(4)；
Wherein, ωsFor reference frequency, unit rad/s；ωoFor preliminary frequency, unit rad/s.
For step S104, according to the preset sample frequency and the reference frequency, the list of the electric power signal is obtained Bit period sequence length；
In one embodiment, the unit period sequence length for obtaining the electric power signal is formula (5)：
In formula, N2πFor the unit period sequence length, unit dimensionless；(int) round numbers is represented；fnPreset to be described Sample frequency, unit Hz；fsFor the reference frequency of Hz units；ωsFor the reference frequency of rad/s units.
The error in 1 sampling interval be present in the unit period sequence length integer.
For step S105, the default integer signal of the electric power signal is obtained according to the preliminary design integer signal period number Periodicity, and according to the default integer signal period number and the unit period sequence length, obtain the electric power signal Predetermined sequence length；
In one embodiment, integer signal period number is preset into the preliminary design integer signal period number acquisition that subtracts 1, according to Default integer signal period number and the unit period sequence length, predetermined sequence length is obtained, requires predetermined sequence in principle Length corresponds to integer signal period number；
In one embodiment, it is formula (6) to obtain the predetermined sequence length：
N=(int) (C2πN2π)
(6)；
C2π=C2π-start-1
Wherein, N is the predetermined sequence length, unit dimensionless；(int) round numbers is represented；N2πFor the unit period Sequence length, unit dimensionless；C2πFor the default integer signal period number, unit dimensionless；C2π-startIt is whole for the preliminary design Number signal period number, unit dimensionless；Due to error be present, the predetermined sequence length includes the unit period sequence length Error existing for integer.
For step S106, according to predetermined sequence starting point and the predetermined sequence length, obtained from the preliminary sequence The first positive sequence of the electric power signal is obtained, in one embodiment, predetermined sequence starting point value is the unit period sequence 0.5 times of row length；
In one embodiment, it is formula (7) to obtain the first positive sequence：
Wherein, X1(n) it is the first positive sequence；Xstart(n) it is preliminary sequence；NsetFor the predetermined sequence of dimensionless unit Starting point, initial value 0.5N2π；(int) it is round numbers；A is signal amplitude, unit v；ω is signal frequency, unit rad/s；Tn For sampling interval, unit s；N is series of discrete number, unit dimensionless；For signal initial phase, unit rad；For first just To sequence initial phase, unit rad；Predetermined sequence length described in N；NstartFor preliminary sequence length, unit dimensionless.
The avatars of the first positive sequence, as shown in Figure 2.
For step S107, the described first positive sequence is reversely exported, obtained corresponding with the described first positive sequence First anti-pleat sequence；
In one embodiment, it is formula (8) to obtain the first anti-pleat sequence：
X-1(- n)=X1(N-n)=Acos (- ω Tnn+β1)
(8)；
N=0,1,2,3 ..., N-1
In formula, X-1(- n) is the first anti-pleat sequence；β 1 is the first anti-pleat sequence initial phase, unit rad.In relation, first Anti- pleat sequence initial phase is the cut-off phase of the first positive sequence；N is the predetermined sequence length.
The avatars of the first anti-pleat sequence, as shown in Figure 2.
For step S108, respectively by the SIN function and institute of the cosine function of the reference frequency and the reference frequency State the first positive sequence to be multiplied, generate 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 (9) 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；Ω is the letter of rad/s units The frequency difference of number frequency and reference frequency；N is the predetermined sequence length.
For step S109, respectively by the SIN function and institute of the cosine function of the reference frequency and the reference frequency State the first anti-pleat sequence to be multiplied, generate 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 (10) 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；Ω is the letter of rad/s units The frequency difference of number frequency and reference frequency；Predetermined sequence length described in N.
For step S110, the described first real frequency sequence vector and the first empty frequency sequence vector are integrated respectively Calculate, 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 (11)：
Wherein, R1For the first real vectorial integrated value of frequency；I1For the first empty vectorial integrated value of frequency.
For step S111, the described second real frequency sequence vector and the second empty frequency sequence vector are integrated respectively Calculate, 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.
, will the vectorial integrated value of the described first empty frequency and described the according to default phase transition rule for step S112 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 empty vectorial integrated value of frequency and the described first real vectorial integrated value of frequency Second ratio of the two empty vectorial integrated values of frequency and the described second real vectorial integrated value of frequency；
The opposite number of the arctan function value of first ratio is obtained, generates the first phase, obtains described second The opposite number of the arctan function value of ratio, generates the second phase.
For step S113, 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 described first positive sequence；
In one embodiment, the first phase and the second phase can be converted to by institute by below equation (15) State the average initial phase of the first positive sequence：
Wherein, PH01For the average initial phase of the first positive sequence, unit rad；PH1For the first phase of rad units；PH2 For the second phase of rad units.
The average initial phase is the described first positive sequence initial phase and the described first anti-pleat sequence initial phase is averaged Value.
For step S114, by the phase of the average initial phase and ± π/4 compared with, obtain and ± the π/4 phases The phase compared compares difference；
In one embodiment, phase bit comparison can be carried out by below equation (16)：
Δ PH in formula±π/4Compare difference, unit rad for phase；
For step S115, difference, the predetermined sequence starting point and the unit period sequence are compared according to the phase Row length, obtain new predetermined sequence starting point；
In one embodiment, new predetermined sequence starting point can be calculated by below equation (17)：
N in formulanewFor new predetermined sequence starting point, unit dimensionless；NsetFor the predetermined sequence starting point of dimensionless unit； N2πFor the unit period sequence length of dimensionless unit；(int) it is round numbers.
New predetermined sequence starting point can be by the described second positive sequence initial phase and the described second anti-pleat sequence first phase The average value of positionControl is in ± π/4 or so, it is therefore an objective to ensures the cosine function modulation sequence and the sine The absolute amplitude of FUNCTION MODULATION sequence is essentially identical.
For step S116, according to the new predetermined sequence starting point and the predetermined sequence length, from described preliminary The second positive sequence of the electric power signal is obtained in sequence；
In one embodiment, it is formula (18) to obtain the second positive sequence：
Wherein, X2(n) it is the second positive sequence；Xstart(n) it is preliminary sequence；A is signal amplitude, unit v；ω is signal Frequency, unit rad/s；TnFor sampling interval, unit s；N is the series of discrete number of dimensionless unit；NnewFor (dimensionless unit) New predetermined sequence starting point；For preliminary sequence initial phase, unit rad；For the second positive sequence initial phase, unit rad； N is the predetermined sequence length of dimensionless unit；NstartFor the preliminary sequence length of dimensionless unit.
For step S117, the described second positive sequence is reversely exported, obtained corresponding with the described second positive sequence Second anti-pleat sequence；
In one embodiment, obtain the second anti-pleat to sequence be formula (19)：
X-2(- n)=X2(N-n)=Acos (- ω Tnn+β2)
(19)；
N=0,1,2,3 ..., N-1
In formula, X-2(- n) is the second anti-pleat sequence；β 2 is the second anti-pleat sequence initial phase, unit rad；N is dimensionless list The series of discrete number of position.In relation, the second anti-pleat sequence initial phase is the cut-off phase of the second positive sequence；N is described default Sequence length.
For step S118, the described second positive sequence is added with the described second anti-pleat sequence, obtains the electric power letter Number cosine function modulation sequence；
In one embodiment, the described second positive sequence is added with the described second anti-pleat sequence and obtains cosine function tune Sequence processed, it is formula (20)：
In formula, Xcos(m) it is cosine function modulation sequence；AcosCosine function modulation sequence amplitude (can be negative value), unit v；For cosine function modulation sequence initial phase, unit rad；N is the series of discrete number of dimensionless unit；Described in N Predetermined sequence length.
Described second positive sequence and the second anti-pleat sequence are subtracted each other, the SIN function for obtaining the electric power signal is adjusted Sequence processed；
In one embodiment, the described second positive sequence and the second anti-pleat sequence are subtracted each other and obtains cosine function tune Sequence processed, it is formula (21)：
In formula, Xsin(n) it is SIN function modulation sequence；AsinIt is single for SIN function modulation sequence amplitude (can be negative value) Position v；For SIN function modulation sequence initial phase, unit rad；N is the series of discrete number of dimensionless unit；N institutes State predetermined sequence length.
The SIN function modulation sequence and the cosine function modulation sequence are orthogonal each other, SIN function modulation sequence with The absolute amplitude of cosine function modulation sequence is basically identical, has excellent estimated performance, improves the accurate of sine parameter calculating Degree and anti-harmonic wave, noise jamming.
The system for carrying out orthogonal modulation in one embodiment to electric power signal, as shown in figure 3, including：
Preliminary sequence length modules 301, for the lower limit according to frequency power signal scope, preset sample frequency and preliminary design Integer signal period number, obtain the preliminary sequence length of the electric power signal；
Preliminary sequence module 302, for tentatively being sampled to the electric power signal according to the preliminary sequence length, obtain Obtain the preliminary sequence of the electric power signal；
Frequency preliminary survey module 303, for entering line frequency preliminary survey to the preliminary sequence, generate the preliminary of the electric power signal Frequency, and according to the reference frequency of the preliminary frequency setting electric power signal；
Unit period sequence length module 304, for according to the preset sample frequency and the reference frequency, obtaining institute State the unit period sequence length of electric power signal；
Predetermined sequence length modules 305, for obtaining the electric power signal according to the preliminary design integer signal period number Default integer signal period number, and according to the default integer signal period number and the unit period sequence length, obtain institute State the predetermined sequence length of electric power signal；
First positive block 306, for according to predetermined sequence starting point and the predetermined sequence length, from it is described just The first positive sequence of the electric power signal is obtained in step sequence；
First anti-pleat block 307, for the described first positive sequence reversely to be exported, obtain and the described first forward direction First anti-pleat sequence corresponding to sequence；
Primary vector sequence generating module 308, for respectively by the cosine function of the reference frequency and the reference frequency The SIN function of rate is multiplied with the described first positive sequence, generates the first real frequency sequence vector and the first empty frequency sequence vector；
Secondary vector sequence generating module 309, for respectively by the cosine function of the reference frequency and the reference frequency The SIN function of rate is multiplied with the described first anti-pleat sequence, generates the second real frequency sequence vector and the second empty frequency sequence vector；
Primary vector integrated value generation module 310, for respectively to the described first real frequency sequence vector and first void Frequency sequence vector carries out integral and calculating, the vectorial integrated value of the real frequency of generation first and the first empty vectorial integrated value of frequency；
Secondary vector integrated value generation module 311, for respectively to the described second real frequency sequence vector and second void Frequency sequence vector carries out integral and calculating, the vectorial integrated value of the real frequency of generation second and the second empty vectorial integrated value of frequency；
Phase module 312, 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；
Phase detecting module 313, 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 described first positive sequence；
Phase compares difference block 314, for by compared with the phase of the average initial phase and ± π/4, acquisition and institute The phase for stating the phase bit comparison of ± π/4 compares difference；
Predetermined sequence originates point module 315, for comparing difference, the predetermined sequence starting point and institute according to the phase Unit period sequence length is stated, obtains new predetermined sequence starting point；
Second positive block 316, for according to the new predetermined sequence starting point and the predetermined sequence length, The second positive sequence of the electric power signal is obtained from the preliminary sequence；
Second anti-pleat block 317, for the described second positive sequence reversely to be exported, obtain and the described second forward direction Second anti-pleat sequence corresponding to sequence；
Orthogonal modulation module 318, for the described second positive sequence to be added with the described second anti-pleat sequence, described in acquisition The cosine function modulation sequence of electric power signal, the described second positive sequence and the second anti-pleat sequence are subtracted each other, described in acquisition The SIN function modulation sequence of electric power signal.
Based on the system of the present embodiment shown in Fig. 3, a specific course of work can be discussed further below：
Preliminary sequence length modules 301 are according to the lower limit of frequency power signal scope, preset sample frequency and preliminary design integer Signal period number, obtain preliminary sequence length；Preliminary sequence module 302 is carried out according to preliminary sequence length to the electric power signal Preliminary sampling, obtains preliminary sequence；Frequency preliminary survey module 303 enters line frequency preliminary survey to preliminary sequence, generates preliminary frequency, and root Reference frequency is set according to preliminary frequency；Unit period sequence length module 304 obtains according to preset sample frequency and reference frequency Unit period sequence length；Predetermined sequence length modules 305 obtain default integer signal week according to preliminary design integer signal period number Issue, and according to default integer signal period number and unit period sequence length, obtain predetermined sequence length；First positive sequence Module 306 obtains the first positive sequence according to predetermined sequence starting point and predetermined sequence length from preliminary sequence；First anti-pleat Block 307 reversely exports the first positive sequence, obtains the first anti-pleat sequence；Primary vector sequence generating module 308 divides The cosine function of reference frequency and SIN function are not multiplied with the first positive sequence, generate the first real frequency sequence vector and first Empty frequency sequence vector；Secondary vector sequence generating module 309 is respectively by the cosine function of reference frequency and SIN function and first Anti- pleat sequence is multiplied, and generates the second real frequency sequence vector and the second empty frequency sequence vector；Primary vector integrated value generation module 310 Carry out integral and calculating to the first real frequency sequence vector and the first empty frequency sequence vector respectively, the vectorial integrated value of the real frequency of generation first and The first empty vectorial integrated value of frequency；Secondary vector integrated value generation module 311 is respectively to the second real frequency sequence vector and the second empty frequency Sequence vector carries out integral and calculating, the vectorial integrated value of the real frequency of generation second and the second empty vectorial integrated value of frequency；Phase module 312 According to default phase transition rule, the first empty vectorial integrated value of frequency and the first real vectorial integrated value of frequency are converted into first phase, The second empty vectorial integrated value of frequency and the second real vectorial integrated value of frequency are converted into second phase；Phase detecting module 313 is according to pre- If sequence average initial phase transformation rule, first phase and second phase are converted to the average first phase of the first positive sequence Position；Phase compares difference block 314 will be compared with the phase of average initial phase and ± π/4, the phase of the phase bit comparison of acquisition and ± π/4 Bit comparison difference；Predetermined sequence starting point module 315 compares difference, predetermined sequence starting point and unit period sequence according to phase Length, obtain new predetermined sequence starting point；Second positive block 316 is according to new predetermined sequence starting point and default sequence Row length, the second positive sequence is obtained from preliminary sequence；Second anti-pleat block 317 is reversely defeated by the second positive sequence Go out, obtain the second anti-pleat sequence；Second positive sequence is added by orthogonal modulation module 318 with the second anti-pleat sequence, obtains cosine FUNCTION MODULATION sequence, the second positive sequence and the second anti-pleat sequence are subtracted each other, obtain SIN function modulation sequence.SIN function is adjusted Sequence processed is orthogonal each other with cosine function modulation sequence, the absolute amplitude of SIN function modulation sequence and cosine function modulation sequence It is basically identical.
It is evidenced from the above discussion that the present invention obtains the basically identical orthogonal modulation sequences of amplitude, have excellent calculating special Property, improve the degree of accuracy and anti-harmonic wave, noise jamming that sine parameter calculates.
Wherein, the preliminary sequence length modules 301 according to the lower limit of frequency power signal scope, preset sample frequency and Preliminary design integer signal period number, obtain the preliminary sequence length of the electric power signal；
The electric power signal is the sinusoidal signal based on a kind of fundamental wave component.Sinusoidal signal extensively criticizes string function signal and remaining String function signal.
In one embodiment, power system frequency scope is in 45Hz-55Hz, power taking force signal lower-frequency limit fminFor 45Hz；According to being actually needed, the preliminary design integer signal period number C is set2π-start, it is preferable that take C2π-startFor integer 13.
In one embodiment, it is formula (1) to obtain the preliminary sequence length：
In formula, NstartFor the preliminary sequence length, unit dimensionless；(int) round numbers is represented；C2π-startTo be described Preliminary design integer signal period number, unit dimensionless；fminFor the lower limit of the signal frequency range, unit Hz；fnPreset to be described Sample frequency, unit Hz.
The preliminary sequence module 302 is tentatively sampled according to the preliminary sequence length to the electric power signal, is obtained Obtain the preliminary sequence of the electric power signal.In one embodiment, the electric power signal is believed for the cosine function of single fundamental frequency Number, the preliminary sequence for obtaining the electric power signal is formula (2)：
Wherein, Xstart(n) it is preliminary sequence；A be the electric power signal amplitude, unit v；ω is the electric power signal Frequency, unit rad/s；TnFor sampling interval, unit s；fnFor the preset sample frequency, unit Hz；N is the electric power signal Series of discrete number, unit dimensionless；For signal (preliminary sequence) initial phase, unit rad；NstartFor preliminary sequence length, Unit dimensionless.
The frequency preliminary survey module 303 can by zero friendship method, the algorithm based on filtering, based on Wavelet Transformation Algorithm, be based on The algorithm of neutral net, the frequency algorithm based on DFT transform are carried out based on the frequency algorithm of phase difference to the preliminary sequence Frequency preliminary survey, generate the preliminary frequency of the electric power signal.
In one embodiment, it is formula (3) to generate the preliminary frequency：
ωo(3)；
Wherein, ωoFor the preliminary frequency, unit rad/s；
Preferably, the reference frequency is formula (4) equal to the preliminary frequency：
ωs=ωo(4)；
Wherein, ωsFor reference frequency, unit rad/s；ωoFor preliminary frequency, unit rad/s.
The unit period sequence length module 304 obtains institute according to the preset sample frequency and the reference frequency State the unit period sequence length of electric power signal；
In one embodiment, the unit period sequence length for obtaining the electric power signal is formula (5)：
In formula, N2πFor the unit period sequence length, unit dimensionless；(int) round numbers is represented；fnPreset to be described Sample frequency, unit Hz；fsFor the reference frequency of Hz units；ωsFor the reference frequency of rad/s units.
The error in 1 sampling interval be present in the unit period sequence length integer.
The predetermined sequence length modules 305 obtain the pre- of the electric power signal according to the preliminary design integer signal period number If integer signal period number, and according to the default integer signal period number and the unit period sequence length, described in acquisition The predetermined sequence length of electric power signal；
In one embodiment, integer signal period number is preset into the preliminary design integer signal period number acquisition that subtracts 1, according to Default integer signal period number and the unit period sequence length, predetermined sequence length is obtained, requires predetermined sequence in principle Length corresponds to integer signal period number；
In one embodiment, it is formula (6) to obtain the predetermined sequence length：
N=(int) (C2πN2π)
(6)；
C2π=C2π-start-1
Wherein, N is the predetermined sequence length, unit dimensionless；(int) round numbers is represented；N2πFor the unit period Sequence length, unit dimensionless；C2πFor the default integer signal period number, unit dimensionless；C2π-startIt is whole for the preliminary design Number signal period number, unit dimensionless；Due to error be present, the predetermined sequence length includes the unit period sequence length Error existing for integer.
The first positive block 306 is according to predetermined sequence starting point and the predetermined sequence length, at the beginning of described The first positive sequence of the electric power signal is obtained in step sequence, in one embodiment, predetermined sequence starting point value is described 0.5 times of unit period sequence length；
In one embodiment, it is formula (7) to obtain the first positive sequence：
Wherein, X1(n) it is the first positive sequence；Xstart(n) it is preliminary sequence；NsetFor the predetermined sequence of dimensionless unit Starting point, initial value 0.5N2π；(int) it is round numbers；A is signal amplitude, unit v；ω is signal frequency, unit rad/s；Tn For sampling interval, unit s；N is series of discrete number, unit dimensionless；For signal initial phase, unit rad；It is positive for first Sequence initial phase, unit rad；Predetermined sequence length described in N；NstartFor preliminary sequence length, unit dimensionless.
The avatars of the first positive sequence, as shown in Figure 2.
The first anti-pleat block 307 reversely exports the described first positive sequence, obtains and the described first forward direction First anti-pleat sequence corresponding to sequence；
In one embodiment, it is formula (8) to obtain the first anti-pleat sequence：
X-1(- n)=X1(N-n)=Acos (- ω Tnn+β1)
(8)；
N=0,1,2,3 ..., N-1
In formula, X-1(- n) is the first anti-pleat sequence；β 1 is the first anti-pleat sequence initial phase, unit rad.In relation, first Anti- pleat sequence initial phase is the cut-off phase of the first positive sequence；N is the predetermined sequence length.
The avatars of the first anti-pleat sequence, as shown in Figure 2.
The primary vector sequence generating module 308 is respectively by the cosine function of the reference frequency and the reference frequency SIN function be multiplied with the described first positive sequence, generate 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 (9) 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；Ω is the letter of rad/s units The frequency difference of number frequency and reference frequency；N is the predetermined sequence length.
The secondary vector sequence generating module 309 is respectively by the cosine function of the reference frequency and the reference frequency SIN function be multiplied with the described first anti-pleat sequence, generate 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 (10) 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；Ω is the letter of rad/s units The frequency difference of number frequency and reference frequency；Predetermined sequence length described in N.
The primary vector integrated value generation module 310 is respectively to the described first real frequency sequence vector and the first empty frequency Sequence vector carries out integral and calculating, 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 (11)：
Wherein, R1For the first real vectorial integrated value of frequency；I1For the first empty vectorial integrated value of frequency.
The secondary vector integrated value generation module 311 is respectively to the described second real frequency sequence vector and the second empty frequency Sequence vector carries out integral and calculating, 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.
The phase module 312 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 empty vectorial integrated value of frequency and the described first real vectorial integrated value of frequency Second ratio of the two empty vectorial integrated values of frequency and the described second real vectorial integrated value of frequency；
The opposite number of the arctan function value of first ratio is obtained, generates the first phase, obtains described second The opposite number of the arctan function value of ratio, generates the second phase.
The phase detecting module 313 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 described first positive sequence；
In one embodiment, the first phase and the second phase can be converted to by institute by below equation (15) State the average initial phase of the first positive sequence：
Wherein, PH01For the average initial phase of the first positive sequence, unit rad；PH1For the first phase of rad units；PH2 For the second phase of rad units.
The average initial phase is the described first positive sequence initial phase and the described first anti-pleat sequence initial phase is averaged Value.
The phase compares difference block 314 by compared with the phase of the average initial phase and ± π/4, acquisition with it is described The phase of the phase bit comparison of ± π/4 compares difference；
In one embodiment, phase bit comparison can be carried out by below equation (16)：
Δ PH in formula±π/4Compare difference, unit rad for phase；
The predetermined sequence starting point module 315 compares difference, the predetermined sequence starting point and institute according to the phase Unit period sequence length is stated, obtains new predetermined sequence starting point；
In one embodiment, new predetermined sequence starting point can be calculated by below equation (17)：
N in formulanewFor new predetermined sequence starting point, unit dimensionless；NsetFor the predetermined sequence starting point of dimensionless unit； N2πFor the unit period sequence length of dimensionless unit；(int) it is round numbers.
New predetermined sequence starting point can be by the described second positive sequence initial phase and the described second anti-pleat sequence first phase The average value 0.5 (φ 2+ β 2) of position is controlled in ± π/4 or so, it is therefore an objective to ensures the cosine function modulation sequence and the sine The absolute amplitude of FUNCTION MODULATION sequence is essentially identical.
The second positive block 316 according to the new predetermined sequence starting point and the predetermined sequence length, The second positive sequence of the electric power signal is obtained from the preliminary sequence；
In one embodiment, it is formula (18) to obtain the second positive sequence：
Wherein, X2(n) it is the second positive sequence；Xstart(n) it is preliminary sequence；A is signal amplitude, unit v；ω is signal Frequency, unit rad/s；TnFor sampling interval, unit s；N is the series of discrete number of dimensionless unit；NnewFor (dimensionless unit) New predetermined sequence starting point；For preliminary sequence initial phase, unit rad；For the second positive sequence initial phase, unit rad； N is the predetermined sequence length of dimensionless unit；NstartFor the preliminary sequence length of dimensionless unit.
The second anti-pleat block 317 reversely exports the described second positive sequence, obtains and the described second forward direction Second anti-pleat sequence corresponding to sequence；
In one embodiment, obtain the second anti-pleat to sequence be formula (19)：
X-2(- n)=X2(N-n)=Acos (- ω Tnn+β2)
(19)；
N=0,1,2,3 ..., N-1
In formula, X-2(- n) is the second anti-pleat sequence；β 2 is the second anti-pleat sequence initial phase, unit rad；N is dimensionless list The series of discrete number of position.In relation, the second anti-pleat sequence initial phase is the cut-off phase of the second positive sequence；N is described default Sequence length.
Described second positive sequence is added by the orthogonal modulation module 318 with the described second anti-pleat sequence, described in acquisition The cosine function modulation sequence of electric power signal；
In one embodiment, the described second positive sequence is added with the described second anti-pleat sequence and obtains cosine function tune Sequence processed, it is formula (20)：
In formula, Xcos(m) it is cosine function modulation sequence；AcosCosine function modulation sequence amplitude (can be negative value), unit v；For cosine function modulation sequence initial phase, unit rad；N is the series of discrete number of dimensionless unit；Described in N Predetermined sequence length.
Described second positive sequence and the second anti-pleat sequence are subtracted each other, the SIN function for obtaining the electric power signal is adjusted Sequence processed；
In one embodiment, the described second positive sequence and the second anti-pleat sequence are subtracted each other and obtains cosine function tune Sequence processed, it is formula (21)：
In formula, Xsin(n) it is SIN function modulation sequence；AsinIt is single for SIN function modulation sequence amplitude (can be negative value) Position v；For SIN function modulation sequence initial phase, unit rad；N is the series of discrete number of dimensionless unit；N institutes State predetermined sequence length.
The SIN function modulation sequence and the cosine function modulation sequence are orthogonal each other, SIN function modulation sequence with The absolute amplitude of cosine function modulation sequence is basically identical, has excellent estimated performance, improves the accurate of sine parameter calculating Degree and anti-harmonic wave, noise jamming.
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)

  1. A kind of 1. method that orthogonal modulation is carried out to electric power signal, it is characterised in that comprise the following steps：
    According to the lower limit of frequency power signal scope, preset sample frequency and preliminary design integer signal period number, the electric power is obtained The preliminary sequence length of signal；
    The electric power signal is tentatively sampled according to the preliminary sequence length, obtains the preliminary sequence of the electric power signal Row；
    Line frequency preliminary survey is entered to the preliminary sequence, generates the preliminary frequency of the electric power signal, and according to the preliminary frequency Set the reference frequency of the electric power signal；
    According to the preset sample frequency and the reference frequency, the unit period sequence length of the electric power signal is obtained；
    The default integer signal period number of the electric power signal is obtained according to the preliminary design integer signal period number, and according to described Default integer signal period number and the unit period sequence length, obtain the predetermined sequence length of the electric power signal；
    According to predetermined sequence starting point and the predetermined sequence length, the of the electric power signal is obtained from the preliminary sequence One positive sequence；
    Described first positive sequence is reversely exported, obtains the first anti-pleat sequence corresponding with the described first positive sequence；
    Respectively by the cosine function of the reference frequency and the SIN function of the reference frequency and the described first positive sequence phase Multiply, generate the first real frequency sequence vector and the first empty frequency sequence vector；
    Respectively by the cosine function of the reference frequency and the SIN function of the reference frequency and the described first anti-pleat sequence phase Multiply, generate the second real frequency sequence vector and the second empty frequency sequence vector；
    Integral and calculating, 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 integrated value and the first empty vectorial integrated value of frequency；
    Integral and calculating, 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 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；
    It is described regular according to default phase transition, the described first empty vectorial integrated value of frequency and the described first real vector frequently 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 The step of position, including：Obtain the described first empty vectorial integrated value of frequency and the described first real vectorial integrated value of frequency the first ratio and Second ratio of the described second empty vectorial integrated value of frequency and the described second real vectorial integrated value of frequency；Obtain the anti-of first ratio The opposite number of tan value, the first phase is generated, obtain the opposite number of the arctan function value of second ratio, it is raw Into the second phase；
    According to default sequence average initial phase transformation rule, the first phase and the second phase are converted to described The average initial phase of one positive sequence；
    It is described according to default sequence average initial phase transformation rule, the first phase and the second phase are converted into institute The step of stating the average initial phase of the first positive sequence, including：The average value of the first phase and the second phase is calculated, The average initial phase of the described first positive sequence is averagely worth to according to this；
    Compared with the phase of the average initial phase and ± π/4, the phase of the phase bit comparison of acquisition and the ± π/4 is poor Value；
    Difference, the predetermined sequence starting point and the unit period sequence length are compared according to the phase, obtained new pre- If sequence starting point；
    According to the new predetermined sequence starting point and the predetermined sequence length, the electric power is obtained from the preliminary sequence Second positive sequence of signal；
    Described second positive sequence is reversely exported, obtains the second anti-pleat sequence corresponding with the described second positive sequence；
    Described second positive sequence is added with the described second anti-pleat sequence, obtains the cosine function modulation sequence of the electric power signal Row, the described second positive sequence and the second anti-pleat sequence are subtracted each other, and obtain the SIN function modulation sequence of the electric power signal Row.
  2. 2. the method according to claim 1 that orthogonal modulation is carried out to electric power signal, it is characterised in that the electric power signal is single The cosine function signal of fundamental frequency, as N≤NstartWhen, according to expression formula Obtain the described first positive sequence X1(n), wherein N is the predetermined sequence length, NstartFor the preliminary sequence length, Xstart (n) it is the preliminary sequence, NsetFor the predetermined sequence starting point, n is the series of discrete number of the electric power signal, n=0,1, 2,3 ..., N-1, A be the electric power signal amplitude, ω be the electric power signal frequency, For the electric power The initial phase of signal, TnFor sampling interval, fnFor the preset sample frequency.
  3. 3. the method according to claim 1 that orthogonal modulation is carried out to electric power signal, it is characterised in that when OrWhen, according to expression formulaObtain the phase and compare difference DELTA PH±π/4, whenOrWhen, according to expression formulaIt is poor to obtain the phase It is worth Δ PH±π/4, whenObtain the phase and compare difference DELTA PH±π/4=0, wherein PH01For the average initial phase.
  4. 4. the method according to claim 1 that orthogonal modulation is carried out to electric power signal, it is characterised in that according to expression formulaObtain the new predetermined sequence starting point Nnew, wherein NsetFor the predetermined sequence Starting point, (int) represent round numbers, Δ PH±π/4Compare difference, N for the phase2πFor the unit period sequence length.
  5. A kind of 5. system that orthogonal modulation is carried out to electric power signal, it is characterised in that including：
    Preliminary sequence length modules, believe for the lower limit according to frequency power signal scope, preset sample frequency and preliminary design integer Number periodicity, obtain the preliminary sequence length of the electric power signal；
    Preliminary sequence module, for tentatively being sampled to the electric power signal according to the preliminary sequence length, described in acquisition The preliminary sequence of electric power signal；
    Frequency preliminary survey module, for entering line frequency preliminary survey to the preliminary sequence, the preliminary frequency of the electric power signal is generated, and The reference frequency of the electric power signal is set according to the preliminary frequency；
    Unit period sequence length module, for according to the preset sample frequency and the reference frequency, obtaining the electric power The unit period sequence length of signal；
    Predetermined sequence length modules, for obtaining the default integer of the electric power signal according to the preliminary design integer signal period number Signal period number, and according to the default integer signal period number and the unit period sequence length, obtain the electric power letter Number predetermined sequence length；
    First positive block, for according to predetermined sequence starting point and the predetermined sequence length, from the preliminary sequence Middle the first positive sequence for obtaining the electric power signal；
    First anti-pleat block, for the described first positive sequence reversely to be exported, obtain and the described first positive sequence pair The the first anti-pleat sequence answered；
    Primary vector sequence generating module, for respectively by the cosine function of the reference frequency and the sine of the reference frequency Function is multiplied with the described first 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 reference frequency and the sine of the reference frequency Function is multiplied with the described first anti-pleat sequence, generates the second real frequency sequence vector and the second empty frequency sequence vector；
    Primary vector integrated value generation module, for respectively to the described first real frequency sequence vector and the first empty vectorial sequence of frequency Row carry out integral and calculating, the vectorial integrated value of the real frequency of generation first and the first empty vectorial integrated value of frequency；
    Secondary vector integrated value generation module, for respectively to the described second real frequency sequence vector and the second empty vectorial sequence of frequency Row carry out integral and calculating, 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；
    The phase module, specifically for obtaining the described first empty vectorial integrated value of frequency and the described first real vectorial integrated value of frequency Second ratio of the first ratio and the second empty vectorial integrated value of frequency and the described second real vectorial integrated value of frequency；Obtain described The opposite number of the arctan function value of one ratio, the first phase is generated, obtain the arctan function value of second ratio Opposite number, generate the second phase；
    Phase detecting module, for 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 described first positive sequence；
    The phase detecting module, the average value specifically for calculating the first phase and the second phase are flat according to this It is worth to the average initial phase of the described first positive sequence；
    Phase compares difference block, compared with by the phase of the average initial phase and ± π/4, obtains and ± π/4 The phase of phase bit comparison compares difference；
    Predetermined sequence originates point module, for comparing difference, the predetermined sequence starting point and the unit according to the phase Periodic sequence length, obtain new predetermined sequence starting point；
    Second positive block, for according to the new predetermined sequence starting point and the predetermined sequence length, from described The second positive sequence of the electric power signal is obtained in preliminary sequence；
    Second anti-pleat block, for the described second positive sequence reversely to be exported, obtain and the described second positive sequence pair The the second anti-pleat sequence answered；
    Orthogonal modulation module, for the described second positive sequence to be added with the described second anti-pleat sequence, obtain the electric power letter Number cosine function modulation sequence, the described second positive sequence and the second anti-pleat sequence are subtracted each other, obtain electric power letter Number SIN function modulation sequence.
  6. 6. the system according to claim 5 that orthogonal modulation is carried out to electric power signal, it is characterised in that the electric power signal For the cosine function signal of single fundamental frequency, as N≤NstartWhen, the first positive block is according to expression formulaObtain the described first positive sequence X1(n), wherein N is described Predetermined sequence length, NstartFor the preliminary sequence length, Xstart(n) it is the preliminary sequence, NsetFor the predetermined sequence Starting point, n be the electric power signal series of discrete number, n=0,1,2,3 ..., N-1, A be the electric power signal amplitude, ω is the frequency of the electric power signal, For the initial phase of the electric power signal, TnFor sampling interval, fnTo be described pre- If sample frequency.
  7. 7. the system according to claim 5 that orthogonal modulation is carried out to electric power signal, it is characterised in that when OrWhen, the phase compares difference block according to expression formulaObtain the phase ratio Compared with difference DELTA PH±π/4, whenOrWhen, the phase compares difference block according to expression formulaObtain the phase and compare difference DELTA PH±π/4, whenThe phase compares difference block Obtain the phase and compare difference DELTA PH±π/4=0, wherein PH01For the average initial phase.
  8. 8. the system according to claim 5 that orthogonal modulation is carried out to electric power signal, it is characterised in that the predetermined sequence Point module is originated according to expression formulaObtain the new predetermined sequence starting point Nnew, Wherein NsetFor the predetermined sequence starting point, (int) represents round numbers, Δ PH±π/4Compare difference, N for the phase2πFor institute State unit period sequence length.
CN201510600888.7A 2015-09-18 2015-09-18 The method and system of orthogonal modulation is carried out to electric power signal Active CN105067880B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510600888.7A CN105067880B (en) 2015-09-18 2015-09-18 The method and system of orthogonal modulation is carried out to electric power signal

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510600888.7A CN105067880B (en) 2015-09-18 2015-09-18 The method and system of orthogonal modulation is carried out to electric power signal

Publications (2)

Publication Number Publication Date
CN105067880A CN105067880A (en) 2015-11-18
CN105067880B true CN105067880B (en) 2018-02-16

Family

ID=54497299

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510600888.7A Active CN105067880B (en) 2015-09-18 2015-09-18 The method and system of orthogonal modulation is carried out to electric power signal

Country Status (1)

Country Link
CN (1) CN105067880B (en)

Families Citing this family (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105548694B (en) * 2015-12-04 2018-06-01 广东电网有限责任公司电力科学研究院 Frequency measurement method and system are carried out according to the zero initial phase reference sinusoidal sequence of function
CN105548690B (en) * 2015-12-04 2018-05-04 广东电网有限责任公司电力科学研究院 Based on zero initial phase reference sinusoidal function frequency multiplication sequence frequency measuring method and system
CN105548702B (en) * 2015-12-04 2018-01-30 广东电网有限责任公司电力科学研究院 Obtain the method and system for dropping low-frequency zero initial phase benchmark cosine function sequence
CN105548687B (en) * 2015-12-04 2018-01-30 广东电网有限责任公司电力科学研究院 The method and system of any initial phase orthogonal sequence is obtained from electric power signal
CN105403768B (en) * 2015-12-04 2018-05-15 广东电网有限责任公司电力科学研究院 The method and system of cosine function times frequency sequence is obtained from electric power signal
CN105372489B (en) * 2015-12-04 2018-01-30 广东电网有限责任公司电力科学研究院 The method and system of any initial phase cosine function sequence is obtained from electric power signal
CN105425036B (en) * 2015-12-04 2018-09-18 广东电网有限责任公司电力科学研究院 Obtain the method and system of the reduction frequency cosine function sequence of electric power signal
CN105548703B (en) * 2015-12-04 2018-01-30 广东电网有限责任公司电力科学研究院 The method and system of any initial phase SIN function sequence is obtained from electric power signal
CN105301356B (en) * 2015-12-04 2018-04-03 广州粤能电力科技开发有限公司 Obtain the method and system of zero initial phase benchmark cosine function times frequency sequence
CN105548701B (en) * 2015-12-04 2018-05-04 广东电网有限责任公司电力科学研究院 Obtain the method and system of the reduction frequency orthogonal sequence of electric power signal
CN105388359B (en) * 2015-12-04 2018-07-03 广东电网有限责任公司电力科学研究院 The method and system of orthogonal times of frequency sequence is obtained from electric power signal
CN105467212B (en) * 2015-12-04 2018-05-04 广东电网有限责任公司电力科学研究院 Obtain the method and system of the raising frequency cosine function sequence of electric power signal
CN105425031B (en) * 2015-12-04 2018-01-30 广东电网有限责任公司电力科学研究院 The orthogonal frequency multiplication sequence method of zero initial phase benchmark and system are obtained from electric power signal
CN105548692B (en) * 2015-12-04 2018-03-20 广东电网有限责任公司电力科学研究院 The method and system of frequency measurement is carried out according to zero initial phase benchmark orthogonal sequence
CN105548688B (en) * 2015-12-04 2018-01-30 广东电网有限责任公司电力科学研究院 Frequency measurement method and system are carried out according to zero initial phase benchmark cosine function sequence
CN105548698B (en) * 2015-12-04 2018-04-03 广东电网有限责任公司电力科学研究院 Frequency power signal method for trimming and system
CN105372490B (en) * 2015-12-04 2018-01-30 广东电网有限责任公司电力科学研究院 The method and system of frequency measurement based on zero orthogonal times of frequency sequence of initial phase benchmark
CN105425035B (en) * 2015-12-04 2018-09-04 广东电网有限责任公司电力科学研究院 Obtain the method and system for dropping low-frequency zero initial phase reference sinusoidal sequence of function
CN105548700B (en) * 2015-12-04 2018-06-01 广东电网有限责任公司电力科学研究院 Obtain the method and system of the raising frequency orthogonal sequence of electric power signal
CN105548705B (en) * 2015-12-04 2018-01-30 广东电网有限责任公司电力科学研究院 Obtain the method and system of the raising frequency reference orthogonal sequence of electric power signal
CN105548693B (en) * 2015-12-04 2018-06-01 广东电网有限责任公司电力科学研究院 Obtain the method and system of zero initial phase reference sinusoidal function times frequency sequence
CN105548691B (en) * 2015-12-04 2018-10-16 广东电网有限责任公司电力科学研究院 Obtain the method and system of the raising frequency sine sequence of function of electric power signal
CN105467211B (en) * 2015-12-04 2018-09-04 广东电网有限责任公司电力科学研究院 Obtain the method and system of the raising frequency reference cosine function sequence of electric power signal
CN105548706B (en) * 2015-12-04 2018-09-04 广东电网有限责任公司电力科学研究院 Obtain the method and system of the reduction frequency sine sequence of function of electric power signal
CN105425032B (en) * 2015-12-04 2018-03-20 广东电网有限责任公司电力科学研究院 The method and system of SIN function times frequency sequence is obtained from electric power signal
CN105548699B (en) * 2015-12-04 2018-01-30 广东电网有限责任公司电力科学研究院 Obtain the method and system of the raising frequency reference SIN function sequence of electric power signal
CN105425040B (en) * 2015-12-04 2018-03-23 广东电网有限责任公司电力科学研究院 Based on zero initial phase benchmark cosine function frequency multiplication sequence frequency measuring method and system
CN105548689B (en) * 2015-12-04 2018-01-30 广东电网有限责任公司电力科学研究院 Obtain the method and system of the reduction frequency reference orthogonal sequence of electric power signal

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104502706A (en) * 2014-12-29 2015-04-08 广东电网有限责任公司电力科学研究院 Harmonic amplitude measurement method and system of power signal
CN104502700A (en) * 2014-12-29 2015-04-08 广东电网有限责任公司电力科学研究院 Sinusoidal parameter measurement method and system of power signal
CN104535836A (en) * 2014-12-29 2015-04-22 广东电网有限责任公司电力科学研究院 Electric signal fundamental frequency measuring method and system
CN104635045A (en) * 2015-02-05 2015-05-20 广东电网有限责任公司电力科学研究院 Power signal frequency detection method and system based on phase modulation
CN104635044A (en) * 2015-02-05 2015-05-20 广东电网有限责任公司电力科学研究院 Power signal frequency detection method and system based on amplitude modulation

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104502706A (en) * 2014-12-29 2015-04-08 广东电网有限责任公司电力科学研究院 Harmonic amplitude measurement method and system of power signal
CN104502700A (en) * 2014-12-29 2015-04-08 广东电网有限责任公司电力科学研究院 Sinusoidal parameter measurement method and system of power signal
CN104535836A (en) * 2014-12-29 2015-04-22 广东电网有限责任公司电力科学研究院 Electric signal fundamental frequency measuring method and system
CN104635045A (en) * 2015-02-05 2015-05-20 广东电网有限责任公司电力科学研究院 Power signal frequency detection method and system based on phase modulation
CN104635044A (en) * 2015-02-05 2015-05-20 广东电网有限责任公司电力科学研究院 Power signal frequency detection method and system based on amplitude modulation

Also Published As

Publication number Publication date
CN105067880A (en) 2015-11-18

Similar Documents

Publication Publication Date Title
CN105067880B (en) The method and system of orthogonal modulation is carried out to electric power signal
CN104635045B (en) Frequency power signal detection method and system based on phase-modulation
CN104635044B (en) Power signal frequency detection method and system based on amplitude modulation
CN105203843B (en) The average first phase position detecting method and system of electric power signal
CN105092970B (en) Obtain the method and system of the initial phase datum mark of electric power signal sequence SIN function zero
CN105548699B (en) Obtain the method and system of the raising frequency reference SIN function sequence of electric power signal
CN105067885B (en) Electric power signal is converted to the method and system of zero initial phase signal sequence
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
CN105158560B (en) The method and system of orthogonal reference signals sequence is extracted from electric power signal
CN105137172B (en) Electric power signal is converted to the method and system of zero initial phase cosine signal sequence
CN105548689B (en) Obtain the method and system of the reduction frequency reference orthogonal sequence of electric power signal
CN105548705B (en) Obtain the method and system of the raising frequency reference orthogonal sequence of electric power signal
CN105158559B (en) Electric power signal is converted to the method and system of zero initial phase sinusoidal signal sequence
CN105445551B (en) The cut-off method for detecting phases and system of sinusoidal signal
CN105548691B (en) Obtain the method and system of the raising frequency sine sequence of function of electric power signal
CN105548700B (en) Obtain the method and system of the raising frequency orthogonal sequence of electric power signal
CN105548706B (en) Obtain the method and system of the reduction frequency sine sequence of function of electric power signal
CN105548702B (en) Obtain the method and system for dropping low-frequency zero initial phase benchmark cosine function sequence
CN105445552B (en) The first phase position detecting method and system of sinusoidal signal
CN105425036B (en) Obtain the method and system of the reduction frequency cosine function sequence of electric power signal
CN105548693B (en) Obtain the method and system of zero initial phase reference sinusoidal function times frequency sequence
CN105548698B (en) Frequency power signal method for trimming and system
CN105467211B (en) Obtain the method and system of the raising frequency reference cosine function sequence of electric power signal
CN105425035B (en) Obtain the method and system for dropping low-frequency zero initial phase reference sinusoidal sequence of function

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