CN105548687A - Method and system acquiring random initial phase orthogonal sequence from electric power signal - Google Patents

Abstract

The invention relates to a method and a system acquiring a random initial phase orthogonal sequence from an electric power signal. A zero initial phase reference cosine function modulation sequence and a zero initial phase reference sine function modulation sequence are acquired, an initial phase is set, a first multiplication sequence, a second multiplication sequence, a third multiplication sequence and a fourth multiplication sequence are acquired according to a cosine function and a sine function of the initial phase, the zero initial phase reference cosine function modulation sequence and the zero initial phase reference sine function modulation sequence, a cosine function sequence of the initial phase is acquired by subtracting the second multiplication sequence from the first multiplication sequence, and a sine function sequence of the initial phase is acquired by adding the fourth multiplication sequence with the third multiplication sequence. According to the method, the initial phase of the orthogonal sequence can be set according actual demands, influence of input sequence randomness and initial phase indetermination can be avoided, and positive significance on improving sine parameter calculation accuracy is realized.

Description

Method and system for acquiring any initial phase orthogonal sequence from power signal

Technical Field

The invention relates to the technical field of power systems, in particular to a method and a system for acquiring an arbitrary initial phase orthogonal sequence from a power signal.

Background

The frequency measurement, the phase measurement, the amplitude measurement and the like of the power system are all the measurement of sine parameters in nature. Fast Fourier Transform (FFT) and Discrete Fourier Transform (DFT) algorithms are basic mathematical methods for sinusoidal parameter calculation, and have wide application in power systems. However, with the development of sinusoidal parameter measurement technology, the problems of the fast fourier transform algorithm and the discrete fourier transform algorithm are more prominent, and it is difficult to further meet the requirement of the power system on high accuracy calculation of sinusoidal parameters.

In the aspect of measuring sinusoidal parameters of the power system, improved parameter measuring methods such as a zero-crossing method, a filtering-based measuring method, a wavelet transformation-based measuring method, a neural network-based measuring method, a DFT transformation-based measuring method and the like exist. The rated power frequency of the power grid is near 50Hz (hertz), and the frequency is low sinusoidal. Due to the limitations of the actual signal processing technology and the complexity of signal composition, such as the influence of data quantization background noise generated by signal discrete sampling, the problem of spectrum leakage caused by signal sequence truncation is objectively difficult to avoid, the influence of the problem of arbitrary and uncertain initial phases of signals, the influence of the problems of direct current, subharmonic and subharmonic in the signals and the like, the measurement precision of the algorithms is low, and the harmonic and noise interference resistance is poor.

Disclosure of Invention

In view of the above, it is necessary to provide a method and system for acquiring an arbitrary initial phase orthogonal sequence from a power signal, which can ensure high accuracy of sinusoidal parameter calculation.

In order to solve the technical problems, the invention adopts the following technical scheme:

a method for obtaining an arbitrary initial phase quadrature sequence from a power signal, comprising the steps of:

obtaining a preliminary sequence length according to the lower limit of the power signal frequency range, a preset sampling frequency and a preset integer signal period number;

sampling the power signal according to the length of the preliminary sequence to obtain a preliminary sequence of the power signal;

carrying out frequency initial measurement on the initial sequence to obtain an initial frequency of the power signal, and obtaining a reference frequency according to the initial frequency;

obtaining the unit cycle sequence length of the power signal according to the preset sampling frequency and the reference frequency;

obtaining a preset sequence length according to the preset integer signal cycle number and the unit cycle sequence length, wherein the preset sequence length is an odd number;

obtaining a first forward sequence from the preliminary sequence according to the preset sequence length and a preset starting point, and obtaining a first reverse pleat sequence according to the first forward sequence;

obtaining a first positive phase from the first forward sequence and a first anti-phase from the first anti-aliasing sequence;

obtaining a first average initial phase according to the first positive phase and the first negative phase;

obtaining a phase comparison value according to the first average initial phase and a preset phase value, and obtaining a new initial point according to the phase comparison value, the preset initial point and the unit cycle sequence length;

according to the preset sequence length and the new starting point, obtaining a second forward sequence from the preliminary sequence, and obtaining a second reverse pleat sequence according to the second forward sequence;

obtaining a second positive phase from the second forward sequence and a second anti-phase from the second anti-aliased sequence;

obtaining a second average initial phase according to the second positive phase and the second negative phase;

adding the second forward sequence and the second inverse pleat sequence to obtain a sum sequence, and obtaining a cosine function modulation sequence according to the sum sequence and a cosine function value of the second average initial phase;

subtracting the second reverse-folding sequence from the second forward sequence to obtain a difference sequence, and obtaining a sine function modulation sequence according to the difference sequence and the sine function value of the second average initial phase;

outputting from the center point of the cosine function modulation sequence to obtain a zero initial phase reference cosine function modulation sequence, and outputting from the center point of the sine function modulation sequence to obtain a zero initial phase reference sine function modulation sequence;

setting an initial phase, multiplying the cosine function of the initial phase by the zero initial phase reference cosine function modulation sequence to obtain a first multiplication sequence, multiplying the sine function of the initial phase by the zero initial phase reference sine function modulation sequence to obtain a second multiplication sequence, multiplying the sine function of the initial phase by the zero initial phase reference cosine function modulation sequence to obtain a third multiplication sequence, and multiplying the cosine function of the initial phase by the zero initial phase reference sine function modulation sequence to obtain a fourth multiplication sequence;

and subtracting the second multiplication sequence from the first multiplication sequence to obtain a cosine function sequence of the initial phase, and adding the fourth multiplication sequence to the third multiplication sequence to obtain a sine function sequence of the initial phase.

A system for obtaining an arbitrary initial phase quadrature sequence from a power signal, comprising:

the preliminary sequence length determining module is used for obtaining a preliminary sequence length according to the lower limit of the power signal frequency range, a preset sampling frequency and a preset integer signal period number;

the preliminary sequence acquisition module is used for sampling the electric power signal according to the length of the preliminary sequence to obtain a preliminary sequence of the electric power signal;

a reference frequency determining module, configured to perform frequency initial measurement on the preliminary sequence to obtain a preliminary frequency of the power signal, and obtain a reference frequency according to the preliminary frequency;

a unit cycle sequence length determining module, configured to obtain a unit cycle sequence length of the power signal according to the preset sampling frequency and the reference frequency;

a preset sequence length determining module, configured to obtain a preset sequence length according to the preset integer signal cycle number and the unit cycle sequence length, where the preset sequence length is an odd number;

a first sequence obtaining module, configured to obtain a first forward sequence from the preliminary sequence according to the preset sequence length and a preset starting point, and obtain a first anti-aliasing sequence according to the first forward sequence;

a first positive and negative phase determination module for obtaining a first positive phase from the first forward sequence and a first negative phase from the first deconvolution sequence;

a first average initial phase determining module, configured to obtain a first average initial phase according to the first positive phase and the first negative phase;

a new starting point determining module, configured to obtain a phase comparison value according to the first average initial phase and a preset phase value, and obtain a new starting point according to the phase comparison value, the preset starting point, and the unit cycle sequence length;

a second sequence obtaining module, configured to obtain a second forward sequence from the preliminary sequence according to the preset sequence length and the new starting point, and obtain a second inverse-folding sequence according to the second forward sequence;

a second forward-reverse phase determination module, configured to obtain a second positive phase according to the second forward sequence and obtain a second reverse phase according to the second deconvolution sequence;

a second average initial phase determining module, configured to obtain a second average initial phase according to the second positive phase and the second negative phase;

a cosine function modulation sequence determining module, configured to add the second forward sequence and the second inverse-pleated sequence to obtain a sum sequence, and obtain a cosine function modulation sequence according to the sum sequence and a cosine function value of the second average initial phase;

a sinusoidal function modulation sequence determining module, configured to subtract the second forward sequence from the second inverse-folding sequence to obtain a difference sequence, and obtain a sinusoidal function modulation sequence according to the difference sequence and a sinusoidal function value of the second average initial phase;

the zero initial phase modulation sequence acquisition module is used for outputting from the center point of the cosine function modulation sequence to acquire a zero initial phase reference cosine function modulation sequence and outputting from the center point of the sine function modulation sequence to acquire a zero initial phase reference sine function modulation sequence;

a multiplication sequence obtaining module, configured to set an initial phase, multiply the cosine function of the initial phase by the zero initial phase reference cosine function modulation sequence to obtain a first multiplication sequence, multiply the sine function of the initial phase by the zero initial phase reference sine function modulation sequence to obtain a second multiplication sequence, multiply the sine function of the initial phase by the zero initial phase reference cosine function modulation sequence to obtain a third multiplication sequence, and multiply the cosine function of the initial phase by the zero initial phase reference sine function modulation sequence to obtain a fourth multiplication sequence;

and the initial phase orthogonal sequence determining module is used for subtracting the second multiplication sequence from the first multiplication sequence to obtain a cosine function sequence of the initial phase, and adding the fourth multiplication sequence to the third multiplication sequence to obtain a sine function sequence of the initial phase.

According to the method and the system for acquiring the orthogonal sequence of the arbitrary initial phase from the power signal, the cosine function sequence of the initial phase and the sine function sequence of the initial phase are orthogonal to each other to form the orthogonal sequence. The method can set the initial phase of the orthogonal sequence according to actual needs, avoids the influence of the problem that the input sequence is arbitrary and uncertain initial phase, and has positive significance for improving the accuracy of sine parameter calculation.

Drawings

FIG. 1 is a schematic flow chart illustrating an embodiment of a method for obtaining an arbitrary initial phase orthogonal sequence from a power signal according to the present invention;

FIG. 2 is a schematic representation of a preliminary sequence, a first forward sequence, and a first reverse pleat sequence of the present invention;

FIG. 3 is a schematic diagram of a zero initial phase reference point according to the present invention;

fig. 4 is a schematic structural diagram of an embodiment of a system for acquiring an arbitrary initial phase orthogonal sequence from a power signal according to the present invention.

Detailed Description

In order to further explain the technical means and effects of the present invention, the following description of the present invention with reference to the accompanying drawings and preferred embodiments will be made for clarity and completeness.

As shown in fig. 1, a method for obtaining an arbitrary initial phase orthogonal sequence from a power signal includes the steps of:

s101, obtaining a preliminary sequence length according to the lower limit of the power signal frequency range, a preset sampling frequency and a preset integer signal period number;

s102, sampling the electric power signal according to the length of the preliminary sequence to obtain the preliminary sequence of the electric power signal;

s103, carrying out frequency initial measurement on the initial sequence to obtain an initial frequency of the power signal, and obtaining a reference frequency according to the initial frequency;

s104, obtaining the unit cycle sequence length of the power signal according to the preset sampling frequency and the reference frequency;

s105, obtaining a preset sequence length according to the preset integer signal cycle number and the unit cycle sequence length, wherein the preset sequence length is an odd number;

s106, obtaining a first forward sequence from the preliminary sequence according to the preset sequence length and a preset starting point, and obtaining a first reverse pleat sequence according to the first forward sequence;

s107, obtaining a first positive phase according to the first forward sequence, and obtaining a first reverse phase according to the first reverse pleat sequence;

s108, obtaining a first average initial phase according to the first positive phase and the first negative phase;

s109, obtaining a phase comparison value according to the first average initial phase and a preset phase value, and obtaining a new initial point according to the phase comparison value, the preset initial point and the unit cycle sequence length;

s110, obtaining a second forward sequence from the preliminary sequence according to the preset sequence length and the new starting point, and obtaining a second reverse-pleat sequence according to the second forward sequence;

s111, obtaining a second positive phase according to the second positive sequence, and obtaining a second reverse phase according to the second reverse-folding sequence;

s112, obtaining a second average initial phase according to the second positive phase and the second negative phase;

s113, adding the second forward sequence and the second inverse pleat sequence to obtain a sum sequence, and obtaining a cosine function modulation sequence according to the sum sequence and a cosine function value of the second average initial phase;

s114, subtracting the second forward sequence from the second inverse pleat sequence to obtain a difference sequence, and obtaining a sine function modulation sequence according to the difference sequence and the sine function value of the second average initial phase;

s115, outputting from the center point of the cosine function modulation sequence to obtain a zero initial phase reference cosine function modulation sequence, and outputting from the center point of the sine function modulation sequence to obtain a zero initial phase reference sine function modulation sequence;

s116, setting an initial phase, multiplying the cosine function of the initial phase by the zero initial phase reference cosine function modulation sequence to obtain a first multiplication sequence, multiplying the sine function of the initial phase by the zero initial phase reference sine function modulation sequence to obtain a second multiplication sequence, multiplying the sine function of the initial phase by the zero initial phase reference cosine function modulation sequence to obtain a third multiplication sequence, and multiplying the cosine function of the initial phase by the zero initial phase reference sine function modulation sequence to obtain a fourth multiplication sequence;

and S117, subtracting the second multiplication sequence from the first multiplication sequence to obtain a cosine function sequence of the initial phase, and adding the fourth multiplication sequence to the third multiplication sequence to obtain a sine function sequence of the initial phase.

The actual power signal is a sinusoidal signal mainly composed of fundamental wave components, and if not specifically stated, the power signal refers to the fundamental wave signal, and the power signal frequency refers to the fundamental wave frequency. For step S101, the power system frequency range is generally 45Hz (Hertz) -55Hz, so the lower limit f of the power signal frequency range_minMay be taken to be 45 Hz. Presetting integer signal period number C_2πCan be set according to actual needs, for example, C_2πIs taken as13. The preliminary sequence length may be calculated according to equation (1):

$N_{start}=\left(\mathrm{int}\right)\frac{C_{2\mathrm{\π}}f}{f_{min}}---\left(1\right)$

wherein N is_startIs the preliminary sequence length; (int) represents rounding; c_2πIs a preset integer signal period number; f. of_minThe lower limit of the frequency range of the power signal, in Hz; f is the preset sampling frequency in Hz.

For step S102, the power signal can be expressed by a cosine function signal of a single fundamental frequency, and then the initial sequence is formula (2):

wherein, X_start(n) is a preliminary sequence; a is the signal amplitude in v; omega_iIs the signal frequency, T is the sampling interval time, f is the preset sampling frequency, the unit Hz, n is the sequence discrete number,for the initial phase of the preliminary sequence, N_startIs the preliminary sequence length.

For step S103, a frequency preliminary measurement may be performed on the preliminary sequence by a zero-crossing method, a filtering-based algorithm, a wavelet transform algorithm, a neural network-based algorithm, a DFT transform-based frequency algorithm, or a phase difference-based frequency algorithm to obtain a preliminary frequency ω_o. In one embodiment, the reference frequency ω_s＝ω_o。

For step S104, in one embodiment, the length of the unit cycle sequence of the power signal is calculated as formula (3):

$N_{2\mathrm{\π}}=\left(\mathrm{int}\right)\frac{2\mathrm{\π}f}{{\mathrm{\ω}}_s}---\left(3\right)$

wherein N is_2πIs the length of the unit period sequence; (int) is an integer; f is a preset sampling frequency in Hz; omega_sIs the reference frequency. The unit period sequence length integer has an error within 1 sampling interval.

For step S105, the preset sequence length is odd, and in one embodiment, the preset sequence length is calculated by equation (4):

wherein, N is the length of the preset sequence, (int) is an integer, N_2πIs the length of the unit period sequence, C_2πFor a predetermined integer signal periodAnd (4) counting.

For step S106, in an embodiment, the preset starting point may be 0.5 times the length of the unit cycle sequence, and the first forward sequence is represented by equation (5):

wherein, X_start(n) is a preliminary sequence, X_+start(n) is a first forward sequence, P_startTo a predetermined starting point, N_2πIs the unit period sequence length, (int) is an integer, A is the signal amplitude, and the unit v, omega_iIs the signal frequency, T is the sampling interval time, n is the sequence discrete number,the sequence length is preset for the first forward sequence initial phase, N.

The first deconvolution sequence is of formula (6):

$\begin{array}{c}{X}_{-\mathrm{start}}(-n)=X_{+\mathrm{start}}(N-n)=A\cos (-{\mathrm{\ω}}_i\mathrm{Tn}+\mathrm{\β}1)\\ n=\mathrm{0,1,2,3},.....,N-1\end{array}---\left(6\right)$

wherein, X_-start(-n) is the first deconvolution sequence, X_+start(n) is the first forward sequence, A is the signal amplitude in v, ω_iFor signal frequency, T is the sampling interval time, N is the sequence discrete number, β 1 is the first unwrapped sequence initial phase, and N is the predetermined sequence length, as shown in fig. 2, it is a schematic diagram of the initial sequence, the first forward sequence, and the first unwrapped sequence.

For step S107, in one embodiment, a first positive phase is obtained from the result of the quadrature mixing and integration calculation on the first forward sequence; a first anti-phase is obtained from the results of quadrature mixing and integration calculations performed on the first anti-aliasing sequence. When the mixing interference frequency of the quadrature mixing is not considered, the quadrature mixing is expressed as equation (7), and the integral calculation is expressed as equation (8):

wherein R is_+start(n) is a first positive real-frequency mixing sequence, I_+start(n) is the first positive virtual mixing sequenceColumn, R_-start(-n) is the first inverse real mixing sequence, I_-start(-n) is the first inverse virtual mixing sequence, cos ([ omega ])_sTn) or cos (-omega)_sTn) is a discrete cosine function of the reference frequency sin (ω)_sTn) or sin (-omega)_sTn) is a discrete sine function of the reference frequency and omega is the signal frequency omega_iWith reference frequency omega_sT is the sampling interval time, n is the sequence discrete number,the first forward sequence initial phase is β 1, and N is a predetermined sequence length.

Wherein R is_+startFirst positive real-frequency integral value, unit dimensionless, I_+startIs the first positive imaginary frequency integral value, with dimensionless units, R_-startIs the first inverse real-frequency integral value, unit dimensionless, I_-startIs the first inverse virtual mixing integral value, the unit is dimensionless, omega is the signal frequency omega_iWith reference frequency omega_sT is the sampling interval time, N is the sequence discrete number, N is the preset sequence length,the first forward sequence initial phase, β 1 the first unwrapped sequence initial phase, and N the predetermined sequence length.

In one embodiment, the calculation of the first positive phase and the first negative phase is expressed by equation (9):

wherein the pH is_+startIs the first positive phase, PH_-startIs a first antiphase, R_+startIs the first positive real-frequency integral value, unit dimensionless, I_+startIs the first positive imaginary frequency integral value, with dimensionless units, R_-startIs the first inverse real-frequency integral value, unit dimensionless, I_-startIs the first inverse virtual mixing integral value, the unit is dimensionless, omega is the signal frequency omega_iWith reference frequency omega_sT is the sampling interval time, N is the preset sequence length,the first forward sequence initial phase, β 1 is the first unwrapped sequence initial phase.

For step S108, in one embodiment, the first average initial phase calculation method is expressed by equation (10):

wherein the pH is_start-avgIs the first average initial phase, PH_+startIs the first positive phase, PH_-startIn order to be in a first opposite phase,the first forward sequence initial phase, β 1 is the first unwrapped sequence initial phase.

For step S109, in one embodiment, the preset phase value may be ± pi/4; the step of obtaining a phase comparison value according to the first average initial phase and the preset phase value may include:

if the first average initial phase is greater than or equal to 0 and less than or equal to pi/2, subtracting the first average initial phase according to pi/4 to obtain a phase comparison value;

and if the first average initial phase is greater than or equal to-pi/2 and less than or equal to 0, subtracting the first average initial phase according to-pi/4 to obtain a phase comparison value.

Specifically formula (11):

${\mathrm{\ΔPH}}_{com}=\{\begin{array}{cc}\frac{\mathrm{\π}}{4}-{\mathrm{PH}}_{start-avg}& 0\≤{\mathrm{PH}}_{start-avg}\≤\frac{\mathrm{\π}}{2}\\ -\frac{\mathrm{\π}}{4}-{\mathrm{PH}}_{start-avg}& -\frac{\mathrm{\π}}{2}\≤{\mathrm{PH}}_{start-avg}\≤0\\ 0& {\mathrm{PH}}_{start-avg}=\±\frac{\mathrm{\π}}{4}\end{array}---\left(11\right)$

wherein, △ PH_comIs the phase comparison value, in units of rad, PH_start-avgIs the first average initial phase.

In one embodiment, the new starting point is calculated as equation (12):

$P_{new}=P_{start}+\left(\mathrm{int}\right)\left(\frac{{\mathrm{\ΔPH}}_{com}}{2\mathrm{\π}}{N}_{2\mathrm{\π}}\right)---\left(12\right)$

wherein, P_newAs a new starting point, in dimensionless units, P_start△ PH as a predetermined starting point_comFor phase comparison, in units rad, N_2πIs the length of the unit period sequence, and (int) is an integer.

For step S110, the second forward sequence and the second reverse-pleated sequence are equation (13):

wherein, X_+end(n) is a second forward sequence, X_-end(-n) is the second deconvolution sequence, P_newIs a new starting point and has no unit dimension,is the second forward sequence initial phase, β 2 is the second inverse-pleated sequence initial phase, ω_iIs the signal frequency, T is the sampling interval time, N is the sequence discrete number, and N is the preset sequence length.

For step S111, in one embodiment, a second positive phase is obtained from the results of the quadrature mixing and digital filtering of the second forward sequence; a second anti-phase is obtained based on the results of the quadrature mixing and digital filtering of the second anti-aliasing sequence. I.e., the second positive phase and the second negative phase, are based on the results of the quadrature mixing and digital filtering calculations. In one embodiment, the orthogonally mixed second forward sequence and second inverse-pleated sequence may be digitally filtered by a 6-stage rectangular window arithmetic mean filter of 2 filter parameters.

When the mixing interference frequency of the quadrature mixing is not considered, the quadrature mixing is expressed as equation (14), and the filtering calculation of the 6-level rectangular window arithmetic mean filter of the 2 filtering parameters is expressed as equation (15):

wherein R is_+end(n) is a second positive real-frequency mixing sequence, I_+end(n) is a second positive virtual audio mixing sequence, R_-end(-n) is the second inverse real mixing sequence, I_-end(-n) is the second inverse virtual mixing sequence, cos ([ omega ])_sTn) or cos (-omega)_sTn) is a discrete cosine function of the reference frequency sin (ω)_sTn) or sin (-omega)_sTn) is a discrete sine function of the reference frequency and omega is the signal frequency omega_iWith reference frequency omega_sT is the sampling interval time, n is the sequence discrete number,the first forward sequence initial phase, β 1 the first unwrapped sequence initial phase, and N the predetermined sequence length.

Wherein R is_+endFiltering for second positive real digitalFinal value, unit dimensionless; i is_+endThe unit is a second positive virtual frequency digital filtering final value and is dimensionless; r_-endThe second inverse digital filtering final value is a unit dimensionless; i is_-endThe second inverse virtual frequency digital filtering final value is a unit dimensionless; omega is the signal frequency omega_iWith reference frequency omega_sThe frequency difference of (2); k (omega) is the amplitude gain of the digital filtering at the frequency difference omega, and the unit is dimensionless; t is sampling interval time;is the second forward sequence initial phase, β 2 is the second inverse pleat sequence initial phase, N_D1For filtering parameter 1, i.e. for N_D1Adding the continuous discrete values, and then taking the arithmetic mean value of the continuous discrete values as the current filtering value to be output; n is a radical of_D2For filter parameter 2, i.e. for N_D2Adding the continuous discrete values, and then taking the arithmetic mean value of the continuous discrete values as the current filtering value to be output; n is a radical of_DThe sequence length is used for digital filtering, the number of the sequence length is the sum of filtering parameters of a 6-stage rectangular window arithmetic mean filter, and the sequence length is less than or equal to a preset sequence length N.

In one embodiment, the filter parameter N_D1The value is 1.5 times of the length of the unit period sequence of the reference frequency, so that the frequency mixing interference frequency generated by 1/3 subharmonic waves is deeply inhibited; filter parameter N_D2The value is 2 times of the length of the unit period sequence of the reference frequency, so that the frequency mixing interference frequency generated by direct current, 1/2 frequency division, subharmonic and the like is deeply inhibited. The 6-stage rectangular window arithmetic mean filter filtering calculation of 2 filtering parameters needs to use 10.5 times of the length of the signal period sequence.

Filter parameter N_D1And a filter parameter N_D2Calculated as formula (16):

$\begin{array}{c}{N}_{D1}=\left(\mathrm{int}\right)\left(1.5N_{2\mathrm{\π}}\right)\\ N_{D2}=2N_{2\mathrm{\π}}\end{array}---\left(16\right)$

wherein N is_D1Is a digital filtering parameter 1, the unit is dimensionless, (int) is an integer, N_D2For the digital filter parameter 2, unit dimensionless, N_2πIs the unit period sequence length.

In one embodiment, the calculation method of the second positive phase and the second negative phase is expressed by formula (17):

wherein the pH is_+endIs the second positive phase, PH_-endIs a second opposite phase, R_+endIs the second positive real-frequency integral value, unit dimensionless, I_+endIs the second positive imaginary frequency integral value, with dimensionless units, R_-endIs the second inverse real-frequency integral value, unit dimensionless, I_-endIs the second inverse virtual mixing integral value, the unit is dimensionless, omega is the signal frequency omega_iWith reference frequency omega_sFrequency difference of (1), T is sampling interval time, N_DThe sequence length is used for the digital filtering,the second forward sequence initial phase, and β 2 the second unwrapped sequence initial phase.

For step S112, the second average initial phase calculation method, expressed as equation (18):

wherein the pH is_end-avgIs the second average initial phase, PH_+endIs the second positive phase, PH_-endIn order to be in the second opposite phase,the second forward sequence initial phase β 2 is the second forward sequence initial phase.

For step S113, the cosine function modulation sequence is expressed as equation (19):

wherein, X_cos(n) is a cosine function modulation sequence; x_+end(n) is a second forward sequence; x_-end(-n) is a second deconvolution sequence; PH value_end-avgA second average initial phase; a is the cosine function modulation sequence amplitude, unit v;modulating the initial phase, omega, of the sequence for the cosine function_iIs the signal frequency, T is the sampling interval time, N is the sequence discrete number, N isThe length of the sequence is preset, and the sequence length is preset,the second forward sequence initial phase, and β 2 the second unwrapped sequence initial phase.

For step S114, the sine function modulation sequence is expressed as equation (20):

wherein, X_sin(n) is a sine function modulation sequence, X_+end(n) is a second forward sequence, X_-end(-n) is the second deconvolution sequence, PH_end-avgIs the second average initial phase, a is the sine function modulation sequence amplitude, unit v,modulating the initial phase, omega, of the sequence for the cosine function_iSignal frequency, T is sampling interval time, N is sequence discrete number, N is preset sequence length,the second forward sequence initial phase, and β 2 the second unwrapped sequence initial phase.

For step S115, in one embodiment, the zero initial phase reference cosine function modulation sequence is expressed by equation (21):

$\begin{array}{c}X{0}_{cos}\left(n\right)=X_{cos}(\frac{N-1}{2}+n)=Acos\left({\mathrm{\ω}}_i{T}_{n}n\right)\\ n=0,1,2,3,\mathrm{.....},\frac{N-1}{2}-1\end{array}---\left(21\right)$

wherein, X0_sin(n) is zero initial phase reference cosine function modulation sequence, A is cosine function modulation sequence amplitude, unit v, omega_iSignal frequency, T is sampling interval time, N is a sequence discrete number, and N is a preset sequence length.

In one embodiment, the zero initial phase reference sine function modulation sequence is expressed as equation (22):

$\begin{array}{c}X{0}_{\sin }\left(n\right)=X_{sin}(\frac{N-1}{2}+n)=Asin\left({\mathrm{\ω}}_i{T}_{n}n\right)\\ n=0,1,2,3,\mathrm{.....},\frac{N-1}{2}-1\end{array}---\left(22\right)$

wherein, X0_sin(n) is a zero initial phase reference sine function modulation sequence, A is the amplitude of the sine function modulation sequence and the unit v, omega_iSignal frequency, T is sampling interval time, N is a sequence discrete number, and N is a preset sequence length. The zero initial phase reference point pattern is shown in fig. 3.

For step S116, the initial phaseIn the range of 0 to. + -. π/2. In one embodiment, the expressions of the first multiplication sequence, the second multiplication sequence, the third multiplication sequence and the fourth multiplication sequence are obtained as (23):

wherein X1(n) is a first multiplication sequence, X2(n) is a second multiplication sequence, X3(n) is a third multiplication sequence, X4(n) is a fourth multiplication sequence,as a result of the initial phase being the first phase,as a sine function of the initial phase, X0_cos(n) is the zero initial phase reference cosine function modulation sequence,as a cosine function of the initial phase, X0_sin(n) is the zero initial phase reference sine functionThe number modulates the sequence, N is the discrete number of the sequence, and N is the preset sequence length.

For step S117, the cosine function sequence of the initial phase and the sine function sequence of the initial phase are expressed by equation (24):

wherein,is a sequence of cosine functions of the initial phase,is a sequence of sine functions of the initial phase,andthe sequences are orthogonal to each other.

Based on the same inventive concept, the invention also provides a system for acquiring any initial phase orthogonal sequence from the power signal, and the following describes the specific implementation of the system in detail with reference to the accompanying drawings.

As shown in fig. 4, a system for obtaining an arbitrary initial phase orthogonal sequence from a power signal includes:

a preliminary sequence length determining module 101, configured to obtain a preliminary sequence length according to a lower limit of a power signal frequency range, a preset sampling frequency, and a preset integer signal cycle number;

a preliminary sequence obtaining module 102, configured to sample an electric power signal according to the length of the preliminary sequence, so as to obtain a preliminary sequence of the electric power signal;

a reference frequency determining module 103, configured to perform frequency initial measurement on the preliminary sequence to obtain a preliminary frequency of the power signal, and obtain a reference frequency according to the preliminary frequency;

a unit cycle sequence length determining module 104, configured to obtain a unit cycle sequence length of the power signal according to the preset sampling frequency and the reference frequency;

a preset sequence length determining module 105, configured to obtain a preset sequence length according to the preset integer signal cycle number and the unit cycle sequence length, where the preset sequence length is an odd number;

a first sequence obtaining module 106, configured to obtain a first forward sequence from the preliminary sequence according to the preset sequence length and a preset starting point, and obtain a first anti-aliasing sequence according to the first forward sequence;

a first positive and negative phase determining module 107, configured to obtain a first positive phase according to the first forward sequence and obtain a first negative phase according to the first anti-aliasing sequence;

a first average initial phase determining module 108, configured to obtain a first average initial phase according to the first positive phase and the first negative phase;

a new starting point determining module 109, configured to obtain a phase comparison value according to the first average initial phase and a preset phase value, and obtain a new starting point according to the phase comparison value, the preset starting point, and the unit cycle sequence length;

a second sequence obtaining module 110, configured to obtain a second forward sequence from the preliminary sequence according to the preset sequence length and the new starting point, and obtain a second inverse-convolution sequence according to the second forward sequence;

a second positive and negative phase determining module 111, configured to obtain a second positive phase according to the second forward sequence and obtain a second inverse phase according to the second inverse-folding sequence;

a second average initial phase determining module 112, configured to obtain a second average initial phase according to the second positive phase and the second negative phase;

a cosine function modulation sequence determining module 113, configured to add the second forward sequence and the second inverse-pleated sequence to obtain a sum sequence, and obtain a cosine function modulation sequence according to the sum sequence and a cosine function value of the second average initial phase;

a sinusoidal function modulation sequence determining module 114, configured to subtract the second forward sequence from the second inverse-folding sequence to obtain a difference sequence, and obtain a sinusoidal function modulation sequence according to the difference sequence and the sinusoidal function value of the second average initial phase;

a zero initial phase modulation sequence obtaining module 115, configured to output from the cosine function modulation sequence center point to obtain a zero initial phase reference cosine function modulation sequence, and output from the sine function modulation sequence center point to obtain a zero initial phase reference sine function modulation sequence;

a multiplication sequence obtaining module 116, configured to set an initial phase, multiply the cosine function of the initial phase by the zero initial phase reference cosine function modulation sequence to obtain a first multiplication sequence, multiply the sine function of the initial phase by the zero initial phase reference sine function modulation sequence to obtain a second multiplication sequence, multiply the sine function of the initial phase by the zero initial phase reference cosine function modulation sequence to obtain a third multiplication sequence, and multiply the cosine function of the initial phase by the zero initial phase reference sine function modulation sequence to obtain a fourth multiplication sequence;

an initial phase orthogonal sequence determining module 117, configured to subtract the second multiplication sequence from the first multiplication sequence to obtain a cosine function sequence of the initial phase, and add the fourth multiplication sequence to the third multiplication sequence to obtain a sine function sequence of the initial phase.

In one embodiment, the cosine function modulationThe sequence determination module 113 may be based on an expressionObtaining a cosine function modulation sequence X_cos(n) wherein X_+end(n) is a second forward sequence, X_-end(-n) is the second deconvolution sequence, PH_end-avgThe second average initial phase.

In one embodiment, the sine function modulation sequence determination module 114 may be based on an expressionObtaining a sine function modulation sequence X_sin(n) wherein X_+end(n) is a second forward sequence, X_-end(-n) is the second deconvolution sequence, PH_end-avgThe second average initial phase.

In one embodiment, the zero initial phase modulation sequence obtaining module 115 may obtain the zero initial phase modulation sequence according to an expressionObtaining the zero initial phase reference cosine function modulation sequence X0_cos(n) wherein,n is a discrete number of sequences, N is a predetermined sequence length, X_cosAnd (n) is a cosine function modulation sequence.

In one embodiment, the zero initial phase modulation sequence obtaining module 115 may obtain the zero initial phase modulation sequence according to an expressionObtaining the zero initial phase reference sine function modulation sequence X0_sin(n) wherein,n is a discrete number of sequences, N is a predetermined sequenceLength, X_sinAnd (n) is a sine function modulation sequence.

The multiplication sequence acquisition module 116 may be according to an expressionObtaining a first multiplication sequence X1(n), according to the expressionObtaining a second multiplication sequence X2 (n); according to the expressionObtaining a third multiplication sequence X3(n) according to the expressionObtaining a fourth multiplication sequence X4 (n);

wherein,as a result of the initial phase being the first phase,as a sine function of the initial phase, X0_cos(n) is the zero initial phase reference cosine function modulation sequence,as a cosine function of the initial phase, X0_sin(n) is the zero initial phase reference sine function modulation sequence,n is a discrete number of sequences, and N is a preset sequence length.

Other technical features of the system of the present invention are the same as those of the method of the present invention, and are not described herein again.

The technical features of the embodiments described above may be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the embodiments described above are not described, but should be considered as being within the scope of the present specification as long as there is no contradiction between the combinations of the technical features.

The above-mentioned embodiments only express several embodiments of the present invention, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the inventive concept, which falls within the scope of the present invention. Therefore, the protection scope of the present patent shall be subject to the appended claims.

Claims (10)

1. A method for obtaining an arbitrary initial phase orthogonal sequence from a power signal, comprising the steps of:

obtaining a preliminary sequence length according to the lower limit of the power signal frequency range, a preset sampling frequency and a preset integer signal period number;

sampling the power signal according to the length of the preliminary sequence to obtain a preliminary sequence of the power signal;

carrying out frequency initial measurement on the initial sequence to obtain an initial frequency of the power signal, and obtaining a reference frequency according to the initial frequency;

obtaining the unit cycle sequence length of the power signal according to the preset sampling frequency and the reference frequency;

obtaining a preset sequence length according to the preset integer signal cycle number and the unit cycle sequence length, wherein the preset sequence length is an odd number;

obtaining a first forward sequence from the preliminary sequence according to the preset sequence length and a preset starting point, and obtaining a first reverse pleat sequence according to the first forward sequence;

obtaining a first positive phase from the first forward sequence and a first anti-phase from the first anti-aliasing sequence;

obtaining a first average initial phase according to the first positive phase and the first negative phase;

obtaining a phase comparison value according to the first average initial phase and a preset phase value, and obtaining a new initial point according to the phase comparison value, the preset initial point and the unit cycle sequence length;

according to the preset sequence length and the new starting point, obtaining a second forward sequence from the preliminary sequence, and obtaining a second reverse pleat sequence according to the second forward sequence;

obtaining a second positive phase from the second forward sequence and a second anti-phase from the second anti-aliased sequence;

obtaining a second average initial phase according to the second positive phase and the second negative phase;

adding the second forward sequence and the second inverse pleat sequence to obtain a sum sequence, and obtaining a cosine function modulation sequence according to the sum sequence and a cosine function value of the second average initial phase;

subtracting the second reverse-folding sequence from the second forward sequence to obtain a difference sequence, and obtaining a sine function modulation sequence according to the difference sequence and the sine function value of the second average initial phase;

outputting from the center point of the cosine function modulation sequence to obtain a zero initial phase reference cosine function modulation sequence, and outputting from the center point of the sine function modulation sequence to obtain a zero initial phase reference sine function modulation sequence;

setting an initial phase, multiplying the cosine function of the initial phase by the zero initial phase reference cosine function modulation sequence to obtain a first multiplication sequence, multiplying the sine function of the initial phase by the zero initial phase reference sine function modulation sequence to obtain a second multiplication sequence, multiplying the sine function of the initial phase by the zero initial phase reference cosine function modulation sequence to obtain a third multiplication sequence, and multiplying the cosine function of the initial phase by the zero initial phase reference sine function modulation sequence to obtain a fourth multiplication sequence;

and subtracting the second multiplication sequence from the first multiplication sequence to obtain a cosine function sequence of the initial phase, and adding the fourth multiplication sequence to the third multiplication sequence to obtain a sine function sequence of the initial phase.

2. The method of claim 1, wherein the method comprises obtaining the arbitrary initial phase orthogonal sequence from the power signal according to the expressionObtaining a first multiplication sequence X1(n), according to the expressionObtaining a second multiplication sequence X2 (n); according to the expressionObtaining a third multiplication sequence X3(n) according to the expressionObtaining a fourth multiplication sequence X4 (n);

wherein,as a result of the initial phase being the first phase,as a sine function of the initial phase, X0_cos(n) is the zero initial phase reference cosine function modulation sequence,as a cosine function of the initial phase, X0_sin(n) is the zero initial phase reference sine function modulation sequence, n is 0,1,2, 3.N is a discrete number of sequences, and N is a preset sequence length.

3. The method of claim 2, wherein the method comprises obtaining the arbitrary initial phase orthogonal sequence from the power signal according to the expressionObtaining the zero initial phase reference cosine function modulation sequence X0_cos(n), wherein n is 0,1,2, 3. -,n is a discrete number of sequences, N is a predetermined sequence length, X_cosAnd (n) is a cosine function modulation sequence.

4. The method of claim 2, wherein the method comprises obtaining the arbitrary initial phase orthogonal sequence from the power signal according to the expressionObtaining the zero initial phase reference sine function modulation sequence X0_sin(n), wherein n is 0,1,2, 3. -,n is a discrete number of sequences, NFor a predetermined sequence length, X_sinAnd (n) is a sine function modulation sequence.

5. The method of any one of claims 1 to 4, wherein the method comprises the steps of:

according to the expressionObtaining a cosine function modulation sequence X_cos(n) wherein X_+end(n) is a second forward sequence, X_-end(-n) is the second deconvolution sequence, PH_end-avgA second average initial phase;

according to the expressionObtaining a sine function modulation sequence X_sin(n) wherein X_+end(n) is a second forward sequence, X_-end(-n) is the second deconvolution sequence, PH_end-avgThe second average initial phase.

6. A system for obtaining an arbitrary initial phase quadrature sequence from a power signal, comprising:

the preliminary sequence length determining module is used for obtaining a preliminary sequence length according to the lower limit of the power signal frequency range, a preset sampling frequency and a preset integer signal period number;

the preliminary sequence acquisition module is used for sampling the electric power signal according to the length of the preliminary sequence to obtain a preliminary sequence of the electric power signal;

a reference frequency determining module, configured to perform frequency initial measurement on the preliminary sequence to obtain a preliminary frequency of the power signal, and obtain a reference frequency according to the preliminary frequency;

a unit cycle sequence length determining module, configured to obtain a unit cycle sequence length of the power signal according to the preset sampling frequency and the reference frequency;

a preset sequence length determining module, configured to obtain a preset sequence length according to the preset integer signal cycle number and the unit cycle sequence length, where the preset sequence length is an odd number;

a first sequence obtaining module, configured to obtain a first forward sequence from the preliminary sequence according to the preset sequence length and a preset starting point, and obtain a first anti-aliasing sequence according to the first forward sequence;

a first positive and negative phase determination module for obtaining a first positive phase from the first forward sequence and a first negative phase from the first deconvolution sequence;

a first average initial phase determining module, configured to obtain a first average initial phase according to the first positive phase and the first negative phase;

a new starting point determining module, configured to obtain a phase comparison value according to the first average initial phase and a preset phase value, and obtain a new starting point according to the phase comparison value, the preset starting point, and the unit cycle sequence length;

a second sequence obtaining module, configured to obtain a second forward sequence from the preliminary sequence according to the preset sequence length and the new starting point, and obtain a second inverse-folding sequence according to the second forward sequence;

a second forward-reverse phase determination module, configured to obtain a second positive phase according to the second forward sequence and obtain a second reverse phase according to the second deconvolution sequence;

a second average initial phase determining module, configured to obtain a second average initial phase according to the second positive phase and the second negative phase;

a cosine function modulation sequence determining module, configured to add the second forward sequence and the second inverse-pleated sequence to obtain a sum sequence, and obtain a cosine function modulation sequence according to the sum sequence and a cosine function value of the second average initial phase;

a sinusoidal function modulation sequence determining module, configured to subtract the second forward sequence from the second inverse-folding sequence to obtain a difference sequence, and obtain a sinusoidal function modulation sequence according to the difference sequence and a sinusoidal function value of the second average initial phase;

the zero initial phase modulation sequence acquisition module is used for outputting from the center point of the cosine function modulation sequence to acquire a zero initial phase reference cosine function modulation sequence and outputting from the center point of the sine function modulation sequence to acquire a zero initial phase reference sine function modulation sequence;

a multiplication sequence obtaining module, configured to set an initial phase, multiply the cosine function of the initial phase by the zero initial phase reference cosine function modulation sequence to obtain a first multiplication sequence, multiply the sine function of the initial phase by the zero initial phase reference sine function modulation sequence to obtain a second multiplication sequence, multiply the sine function of the initial phase by the zero initial phase reference cosine function modulation sequence to obtain a third multiplication sequence, and multiply the cosine function of the initial phase by the zero initial phase reference sine function modulation sequence to obtain a fourth multiplication sequence;

and the initial phase orthogonal sequence determining module is used for subtracting the second multiplication sequence from the first multiplication sequence to obtain a cosine function sequence of the initial phase, and adding the fourth multiplication sequence to the third multiplication sequence to obtain a sine function sequence of the initial phase.

7. The system of claim 6, wherein the multiplication sequence obtaining module obtains the arbitrary initial phase quadrature sequence from the power signal according to the expressionObtaining a first multiplication sequence X1(n), according to the expressionObtaining a second multiplication sequence X2 (n); according to the expressionObtaining a third multiplication sequence X3(n) according to the expressionObtaining a fourth multiplication sequence X4 (n);

wherein,as a result of the initial phase being the first phase,as a sine function of the initial phase, X0_cos(n) is the zero initial phase reference cosine function modulation sequence,as a cosine function of the initial phase, X0_sin(n) is the zero initial phase reference sine function modulation sequence, n is 0,1,2, 3.N is a discrete number of sequences, and N is a preset sequence length.

8. The system of claim 7, wherein the zero initial phase modulation sequence obtaining module is according to the expressionObtaining the zero initial phase reference cosine function modulation sequence X0_cos(n), wherein n is 0,1,2, 3. -,n is a discrete number of sequences, N is a predetermined sequence length, X_cosAnd (n) is a cosine function modulation sequence.

9. The system for deriving an arbitrary initial phase quadrature sequence from a power signal as claimed in claim 7, wherein said zero initial phase modulation sequence is derived moduloBlock according to expressionObtaining the zero initial phase reference sine function modulation sequence X0_sin(n), wherein n is 0,1,2, 3. -,n is a discrete number of sequences, N is a predetermined sequence length, X_sinAnd (n) is a sine function modulation sequence.

10. A system for deriving an arbitrary initial phase quadrature sequence from a power signal as claimed in any one of claims 6 to 9, wherein:

the cosine function modulation sequence determining module is according to an expressionObtaining a cosine function modulation sequence X_cos(n) wherein X_+end(n) is a second forward sequence, X_-end(-n) is the second deconvolution sequence, PH_end-avgA second average initial phase;

the sine function modulation sequence determination module is based on an expressionObtaining a sine function modulation sequence X_sin(n) wherein X_+end(n) is a second forward sequence, X_-end(-n) is the second deconvolution sequence, PH_end-avgThe second average initial phase.