Summary of the invention
Based on this, be necessary at the problems referred to above, a kind of digital filtering method of three-phase alternating current systematic survey data is provided.
A kind of digital filtering method of three-phase alternating current systematic survey data comprises the steps:
With three-phase voltage signal V
A, V
B, V
CBe converted to the voltage signal V under the two-phase α β coordinate system
α, V
β
According to the transfer function of setting respectively to described voltage signal V
α, V
βCarry out low-pass filtering and obtain voltage signal V'
α, V'
β
Calculate described voltage signal V'
α, V'
βFundamental frequency phase place δ;
Calculate the fundamental frequency phase delay angle δ of described low-pass filtering according to described transfer function
0
According to described fundamental frequency phase place δ with described voltage signal V'
α, V'
βBe converted to the direct voltage V under rotation dq two phase coordinate systems
d, V
q
With δ-δ
0For angle of transformation with described direct voltage V
d, V
qBe converted to the voltage signal V' under three phase coordinates
A, V'
B, V'
C
The digital filtering method of above-mentioned three-phase alternating current systematic survey data, be α β two phase coordinates with three-phase voltage signal from the coordinate transform of ABC three-phase at first, under α β coordinate system, carry out low-pass filtering then, filtered fundamental frequency value is DC quantity under the dq rotating coordinate system by α β coordinate transform, and the fundamental frequency phase delay δ that low-pass filtering produces
0Compensate returning the ABC coordinate time by the dq coordinate transform, thereby solve the fundamental frequency phase delay problem of traditional filtering mode effectively, improve the phase control precision of electric power system.
Embodiment
Be described in detail below in conjunction with the embodiment of accompanying drawing to the digital filtering method of three-phase alternating current systematic survey data of the present invention.
Fig. 1 shows the digital filtering method flow chart of the three-phase alternating current systematic survey data of an embodiment, comprises the steps:
Step S101: with three-phase voltage signal V
A, V
B, V
CBe converted to the voltage signal V under the two-phase α β coordinate system
α, V
β
In one embodiment, also comprise before step S101: the measurement data of three-phase alternating current system is sampled obtains voltage signal V under three phase coordinates
A, V
B, V
C
Particularly, in step S101, will measure the voltage data V that sampling obtains
A, V
B, V
C, be static α β two phase coordinates from the coordinate transform of static ABC three-phase, the formula of conversion is as follows:
V
α=V
A
;
V
α=V
B-V
C
Wherein, ABC → α β coordinate transform is orthogonal transform, and α and A same axis are set.
Step S102: according to the transfer function of setting respectively to described voltage signal V
α, V
βCarry out low-pass filtering and obtain voltage signal V'
α, V'
β
Particularly, computational process comprises following formula:
V'
α=H(jω)·V
α
;
V'
β=H(jω)·V
β
In the formula, the transfer function of H (s) for setting,
Step S103: calculate described voltage signal V'
α, V'
βFundamental frequency phase place δ.
Particularly, computational process comprises following formula:
δ=arctg(V'
β/V'
α);
In the formula, δ is the fundamental frequency phase place, V'
α, V'
βBe respectively the voltage signal of α β coordinate after the filtering.
Step S104: the fundamental frequency phase delay angle δ that calculates described low-pass filtering according to described transfer function
0
Particularly, the angle of transfer function H (j ω) is the phase delay angle δ of low-pass filtering link
0, namely pass through low-pass filtering after, the fundamental frequency phase place δ of voltage signal is with respect to A, B, C voltage in three phases signal delay δ under three phase coordinates under the two-phase α β coordinate system
0, δ
0Computational process comprises following formula:
The transfer function of described low-pass filtering is H (j ω), and its general formula can be expressed as:
Wherein, a
0, a
1, b
0, b
1, b
2Be default constant, the mould value of H (j ω) | H (j ω) | equaling 1, ω is the fundamental frequency angular frequency.
Step S105: according to described fundamental frequency phase place δ with described voltage signal V'
α, V'
βBe converted to the direct voltage V under rotation dq two phase coordinate systems
d, V
q
Particularly, computational process comprises following formula:
V
d+jV
q=(V
α+jV
β)e
-jδ;
In the formula, V'
α, V'
βBe respectively the voltage signal of α β coordinate after the filtering, δ is the fundamental frequency phase place, V
d, V
qBe respectively the voltage signal of dq rotating coordinate system.
Step S106: with δ-δ
0For angle of transformation with described direct voltage V
d, V
qBe converted to the ac voltage signal V' under three phase coordinates
A, V'
B, V'
C
Particularly, variable is returned static ABC three phase coordinates from the coordinate transform of rotation dq two-phase, the leading d axle of q axle is set, angle of transformation is δ-δ
0, its computational process comprises as follows:
In the formula, θ=δ-δ
0, by above-mentioned processing procedure, the fundamental frequency phase delay δ that produces in the low-pass filtering process
0Be compensated in dq → ABC coordinate transform link.
The digital filtering method of three-phase alternating current systematic survey data of the present invention, Measurement and Data Processing at the three-phase alternating current system, at first the three-phase voltage signal (containing high order harmonic component) with sampling is α β two phase coordinates from the coordinate transform of ABC three-phase, under α β coordinate system, carry out low-pass filtering then, filtered fundamental frequency value is DC quantity under the dq rotating coordinate system by α β coordinate transform, and the fundamental frequency phase delay δ that low-pass filtering produces
0Compensate returning the ABC coordinate time by the dq coordinate transform, thereby effectively solve the fundamental frequency phase delay problem of traditional filtering mode.
For clear more technical scheme of the present invention, set forth a simulation example below in conjunction with accompanying drawing.
Referring to shown in Figure 2, Fig. 2 is simulation example laboratory model construction schematic diagram, and wherein the transfer function of low-pass filtering adopts the simplest first-order filtering H (j ω)=G/ (1+Ts), parameter G=1.606, T=0.004 second.
Then Shu Ru ABC three-phase voltage signal is V
A, V
B, V
C, fundamental frequency (50Hz) phase voltage effective value is that 1.0,7 subharmonic content are that 10%, 19 subharmonic content are 180%, and total percent harmonic distortion (THD) is 180.278%, and the voltage signal of exporting after the filtering is V'
A, V'
B, V'
C
As shown in Figure 3, Fig. 3 is the The simulation experiment result schematic diagram, V among the figure
ABe the A phase voltage signal waveform of input, wherein, V'
ABe the A phase voltage signal waveform that carries out through filtering method of the present invention exporting after the filtering, V'
AThe fundamental frequency effective value be 1.0, and phase place and V
AIdentical, 7 subharmonic content are that 1.812%, 19 subharmonic content are 12.007%, and overall percent harmonic distortion (THD) is 12.143%.
Experimental result by above-mentioned simulation example as can be seen, even transfer function H (s) adopts without optimizing and the simplest single order low-pass filter function, filtering method of the present invention also can make total percent harmonic distortion (THD) reduce more than 90%, and the amplitude of fundamental frequency and phase place remain unchanged.
The above embodiment has only expressed several execution mode of the present invention, and it describes comparatively concrete and detailed, but can not therefore be interpreted as the restriction to claim of the present invention.Should be pointed out that for the person of ordinary skill of the art without departing from the inventive concept of the premise, can also make some distortion and improvement, these all belong to protection scope of the present invention.Therefore, the protection range of patent of the present invention should be as the criterion with claims.