CN110704795A - Power signal reconstruction method and system by using shrinkage operator - Google Patents
Power signal reconstruction method and system by using shrinkage operator Download PDFInfo
- Publication number
- CN110704795A CN110704795A CN201910866873.3A CN201910866873A CN110704795A CN 110704795 A CN110704795 A CN 110704795A CN 201910866873 A CN201910866873 A CN 201910866873A CN 110704795 A CN110704795 A CN 110704795A
- Authority
- CN
- China
- Prior art keywords
- opt
- signal sequence
- power signal
- correction
- vector
- 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.)
- Withdrawn
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/327—Testing of circuit interrupters, switches or circuit-breakers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Operations Research (AREA)
- Computing Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The embodiment of the invention discloses a power signal reconstruction method and a system based on a prediction matrix, wherein the method comprises the following steps: step 1, inputting an actually measured power signal sequence S; step 2, carrying out data reconstruction on the power signal sequence S, wherein the reconstructed power signal sequence is SNEW. The method specifically comprises the following steps:wherein α is a contractile factor; mu is a correction factor; zOPTIs the best prediction vector; xOPTIs the best correction vector.
Description
Technical Field
The present invention relates to the field of power, and in particular, to a method and a system for reconstructing a power signal.
Background
With the development of smart grids, the analysis of household electrical loads becomes more and more important. Through the analysis of the power load, a family user can obtain the power consumption information of each electric appliance and a refined list of the power charge in time; the power department can obtain more detailed user power utilization information, can improve the accuracy of power utilization load prediction, and provides a basis for overall planning for the power department. Meanwhile, the power utilization behavior of the user can be obtained by utilizing the power utilization information of each electric appliance, so that the method has guiding significance for the study of household energy consumption evaluation and energy-saving strategies.
The current electric load decomposition is mainly divided into an invasive load decomposition method and a non-invasive load decomposition method. The non-invasive load decomposition method does not need to install monitoring equipment on internal electric equipment of the load, and can obtain the load information of each electric equipment only according to the total information of the electric load. The non-invasive load decomposition method has the characteristics of less investment, convenience in use and the like, so that the method is suitable for decomposing household load electricity.
In the non-invasive load decomposition algorithm, the detection of the switching event of the electrical equipment is the most important link. The initial switch event detection takes the change value of the active power P as the judgment basis of the switch event detection, and is convenient and intuitive. This is because the power consumed by any one of the electric devices changes, and the change is reflected in the total power consumed by all the electric devices. The method needs to set a reasonable threshold value of the power change value, and also needs to solve the problems existing in the practical application of the event detection method, for example, a large peak appears in the instantaneous power value at the starting time of some electric appliances (the starting current of a motor is far larger than the rated current), which causes the inaccurate steady-state power change value of the electric appliances, thereby influencing the judgment of the detection of the switching event; moreover, the transient process of different household appliances is long or short (the duration and the occurrence frequency of impulse noise are different greatly), so that the determination of the power change value becomes difficult; due to the fact that the active power changes suddenly when the quality of the electric energy changes (such as voltage drop), misjudgment is likely to happen. Meanwhile, in the process of acquiring and transmitting the power signal, the operation state of the related instrument and equipment may be temporarily in an abnormal state, which often causes the loss of the power signal.
Therefore, the actual measurement power signal used in the switching event detection process is often incomplete, and the switching event detection cannot be performed correctly by using the incomplete power signal. Therefore, how to effectively reconstruct the incomplete power signal is the key to the success of this method. The existing common method has insufficient attention to the problem, and no effective measure is taken to solve the problem.
Disclosure of Invention
The invention aims to provide a power signal reconstruction method and a system by using a contraction operator. The method has the advantages of good robustness and simple calculation.
In order to achieve the purpose, the invention provides the following scheme:
a method of power signal reconstruction using a puncturing operator, comprising:
step 1, inputting an actually measured power signal sequence S;
step 2, carrying out data reconstruction on the power signal sequence S, wherein the reconstructed power signal sequence is SNEW. The method specifically comprises the following steps:wherein α is a contractile factor; mu is a correction factor; zOPTIs the best prediction vector; xOPTIs the best correction vector.
A power signal reconstruction system using a puncturing operator, comprising:
the acquisition module inputs an actually measured power signal sequence S;
a reconstruction module for performing data reconstruction on the power signal sequence S, wherein the reconstructed power signal sequence is SNEW. The method specifically comprises the following steps:wherein α is a contractile factor; mu is a correction factor; zOPTIs the best prediction vector; xOPTIs the best correction vector.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
although the switching event detection method has wide application in non-invasive load decomposition and is relatively mature in technology, the power signal is often lost in the acquisition and transmission process and is often submerged in pulse noise with strong amplitude, and the switching event detection cannot be correctly performed by using the incomplete power signal. Therefore, how to effectively reconstruct the incomplete power signal is the key to the success of this method. The existing common method has insufficient attention to the problem, and no effective measure is taken to solve the problem.
The invention aims to provide a power signal reconstruction method and a system by using a contraction operator. The method has the advantages of good robustness and simple calculation.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the embodiments will be briefly described below. It is obvious that the drawings in the following description are only some embodiments of the invention, and that for a person skilled in the art, other drawings can be derived from them without inventive effort.
FIG. 1 is a schematic flow diagram of the process of the present invention;
FIG. 2 is a schematic diagram of the system of the present invention;
FIG. 3 is a flow chart illustrating an embodiment of the present invention.
Detailed Description
The technical solution in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. It is to be understood that the described embodiments are merely exemplary of the invention, and not restrictive of the full scope of the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
FIG. 1 is a schematic flow chart of a power signal reconstruction method using a contraction operator
Fig. 1 is a schematic flow chart of a power signal reconstruction method using a contraction operator according to the present invention. As shown in fig. 1, the method for reconstructing a power signal by using a shrinking operator specifically includes the following steps:
step 1, inputting an actually measured power signal sequence S;
step 2, carrying out data reconstruction on the power signal sequence S, wherein the reconstructed power signal sequence is SNEW. The method specifically comprises the following steps:wherein α is a contractile factor; mu is a correction factor; zOPTIs the best prediction vector; xOPTIs the best correction vector.
Before the step 2, the method further comprises:
step 3, obtaining the shrinkage factor alpha, the correction factor mu and the optimal prediction vector ZOPTAnd the optimum correction vector XOPT。
The step 3 comprises the following steps:
step 301, obtaining the shrinkage factor α, specifically as
Wherein
si: the ith element [ i ═ 1,2, …, N of the signal sequence S]
mS: mean value of the signal sequence S
σS: mean square error of the signal sequence S
N: length of the signal sequence S
Step 302, obtaining the correction factor μ, specifically:
wherein
S: the signal sequence
First correction matrix
Step 303, iteratively calculating the optimal prediction vector ZOPTAnd the optimum correction vector XOPTThe method specifically comprises the following steps:
a first step of initialization, in particular
Z1As S: initialized prediction vectors
X1When the ratio is 0: initialized correction vector
n is 1: iterative control parameter
Wherein
0: all zero vector
Second step, update, specifically
Xn+1=Γ[H(Xn)]+Xn
Wherein
The contraction operator is used to perform a contraction operation,
representing a puncturing operation on each element of the matrix B
B: a parameter matrix for representing a contraction operator expression
Thirdly, iterative judgment, specifically comprising
And (4) turning the iteration control parameter n plus 1 to a second step to continuously perform iteration updating until the difference between the iteration results of two adjacent times is less than one thousandth. At this time, the iterative control parameter n is K, and the optimal prediction vector Z is obtainedOPT=ZK+1And the optimum correction vector XOPT=XK+1。
FIG. 2 structural intent of a power signal reconstruction system using a shrink operator
FIG. 2 is a schematic diagram of a power signal reconstruction system using a contraction operator according to the present invention. As shown in fig. 2, the power signal reconstruction system using the puncturing operator includes the following structure:
the acquisition module 401 inputs an actually measured power signal sequence S;
a reconstruction module 402, configured to perform data reconstruction on the power signal sequence S, where the reconstructed power signal sequence is SNEW. The method specifically comprises the following steps:wherein α is a contractile factor; mu is a correction factor; zOPTIs the best prediction vector; xOPTIs the best correction vector.
The system further comprises:
a calculating module 403 for obtaining the shrinkage factor α, the correction factor μ, and the optimal prediction vector ZOPTAnd the optimum correction vector XOPT。
The calculation module 403 further includes the following units:
a first calculation unit 4031 for calculating the shrinkage factor α, specifically
Wherein
si: the ith element [ i ═ 1,2, …, N of the signal sequence S]
mS: mean value of the signal sequence S
σS: mean square error of the signal sequence S
N: length of the signal sequence S
The second calculation unit 4032 calculates the correction factor μ, specifically:
wherein
S: the signal sequence
An iteration unit 4033 for iteratively calculating the optimal prediction vector ZOPTAnd the optimum correction vector XOPT,
The method specifically comprises the following steps:
a first step of initialization, in particular
Z1As S: initialized prediction vectors
X1When the ratio is 0: initialized correction vector
n is 1: iterative control parameter
Wherein
0: all zero vector
Second step, update, specifically
Xn+1=Γ[H(Xn)]+Xn
Wherein
representing a puncturing operation on each element of the matrix B
B: a parameter matrix for representing a contraction operator expression
Thirdly, iterative judgment, specifically comprising
And (4) turning the iteration control parameter n plus 1 to a second step to continuously perform iteration updating until the difference between the iteration results of two adjacent times is less than one thousandth. At this time, the iterative control parameter n is K, and the optimal prediction vector Z is obtainedOPT=ZK+1And the optimum correction vector XOPT=XK+1。
The following provides an embodiment for further illustrating the invention
FIG. 3 is a flow chart illustrating an embodiment of the present invention. As shown in fig. 3, the method specifically includes the following steps:
1. inputting a sequence of measured power signals
S=[s1,s2,…,sN-1,sN]
Wherein:
s: real vibration and sound signal data sequence with length N
siI is 1,2, …, N is measured vibration sound signal with serial number i
2. Determining a shrinkage factor
Wherein
si: the ith element [ i ═ 1,2, …, N of the signal sequence S]
mS: mean value of the signal sequence S
σS: mean square error of the signal sequence S
N: length of the signal sequence S
3. Calculating a correction factor
Wherein
S: the signal sequence
4. Obtaining an optimal prediction vector and an optimal correction vector
First, initialization
Z1As S: initialized prediction vectors
X1When the ratio is 0: initialized correction vector
n is 1: iterative control parameter
Wherein
0: all zero vector
Second, update
Xn+1=Γ[H(Xn)]+Xn
Wherein
The contraction operator is used to perform a contraction operation,
representing a puncturing operation on each element of the matrix B
B: a parameter matrix for representing a contraction operator expression
Third, iterative judgment
And (4) turning the iteration control parameter n plus 1 to a second step to continuously perform iteration updating until the difference between the iteration results of two adjacent times is less than one thousandth. At this time, the iterative control parameter n is K, and the optimal prediction vector Z is obtainedOPT=ZKAnd the optimum correction vector XOPT=XK。+1
5. Data reconstruction
Carrying out data reconstruction on the power signal sequence S, wherein the reconstructed power signal sequence is SNEW. The method specifically comprises the following steps:wherein α is a contractile factor; mu is a correction factor; zOPTIs the best prediction vector; xOPTIs the best correction vector.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. For the system disclosed by the embodiment, the description is simple because the system corresponds to the method disclosed by the embodiment, and the relevant part can be referred to the method part for description.
The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.
Claims (5)
1. A method for power signal reconstruction using a puncturing operator, comprising:
step 1, inputting an actually measured power signal sequence S;
step 2, carrying out data reconstruction on the power signal sequence S, wherein the reconstructed power signal sequence is SNEW. The method specifically comprises the following steps:wherein α is a contractile factor; mu is a correction factor; zOPTIs the best prediction vector; xOPTIs the best correction vector.
2. The method of claim 1, wherein prior to step 2, the method further comprises:
step 3, obtaining the shrinkage factor alpha, the correction factor mu and the optimal prediction vector ZOPTAnd the optimum correction vector XOPT。
3. The method of claim 2, wherein step 3 comprises:
step 301, obtaining the shrinkage factor α, specifically as
Wherein
si: the ith element [ i ═ 1,2, …, N of the signal sequence S]
mS: mean value of the signal sequence S
σS: mean square error of the signal sequence S
N: length of the signal sequence S
Step 302, obtaining the correction factor μ, specifically:
wherein
S: the signal sequence
Step 303, iteratively calculating the optimal prediction vector ZOPTAnd the optimum correction vector XOPTThe method specifically comprises the following steps:
a first step of initialization, in particular
Z1As S: initialized prediction vectors
X1When the ratio is 0: initialized correction vector
n is 1: iterative control parameter
Wherein
0: all zero vector
Second step, update, specifically
Xn+1=Γ[H(Xn)]+Xn
Wherein
The contraction operator is used to perform a contraction operation,
representing a puncturing operation on each element of the matrix B
B: a parameter matrix for representing a contraction operator expression
And thirdly, performing iterative judgment, specifically, turning the iterative control parameter n plus 1 to the second step for continuous iterative updating until the difference between the two adjacent iterative results is less than one per thousand. At this time, the iterative control parameter n is K, and the optimal prediction vector Z is obtainedOPT=ZK+1And the optimum correction vector XOPT=XK+1。
4. A system for power signal reconstruction using a puncturing operator, comprising:
the acquisition module inputs an actually measured power signal sequence S;
a reconstruction module for performing data reconstruction on the power signal sequence S, wherein the reconstructed power signal sequence is SNEW. The method specifically comprises the following steps:wherein α is a contractile factor; mu is a correction factor; zOPTIs the best prediction vector; xOPTIs the best correction vector.
5. The system of claim 4, further comprising:
a calculation module for obtaining the contraction factor alpha, the correction factor mu and the optimal prediction vector ZOPTAnd the optimum correction vector XOPT。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910866873.3A CN110704795A (en) | 2019-09-13 | 2019-09-13 | Power signal reconstruction method and system by using shrinkage operator |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910866873.3A CN110704795A (en) | 2019-09-13 | 2019-09-13 | Power signal reconstruction method and system by using shrinkage operator |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110704795A true CN110704795A (en) | 2020-01-17 |
Family
ID=69194984
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910866873.3A Withdrawn CN110704795A (en) | 2019-09-13 | 2019-09-13 | Power signal reconstruction method and system by using shrinkage operator |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110704795A (en) |
-
2019
- 2019-09-13 CN CN201910866873.3A patent/CN110704795A/en not_active Withdrawn
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108918930B (en) | Power signal self-adaptive reconstruction method in load decomposition | |
CN111680590A (en) | Power signal filtering method and system by using contraction gradient | |
CN111666870A (en) | Power signal reconstruction method and system by utilizing quadratic constraint | |
CN110704800A (en) | Power signal reconstruction method and system based on expansion coefficient | |
CN110704795A (en) | Power signal reconstruction method and system by using shrinkage operator | |
CN110542855B (en) | Load switch event detection method and system based on discrete cosine transform | |
CN110531149B (en) | Power signal filtering method and system based on waveform regularization | |
CN110221119B (en) | Load switch event detection method and system based on power and akie fusion information | |
CN110244115B (en) | Load switch event detection method and system based on signal connectivity | |
CN110514884B (en) | Power signal filtering method and system based on delay vector | |
CN110705426B (en) | Power signal filtering method and system by using deblurring operator | |
CN110702981A (en) | Load switch event detection method and system using classification tree | |
CN110196354B (en) | Method and device for detecting switching event of load | |
CN110611310A (en) | Power signal reconstruction method and system by using periodic boundary conditions | |
CN111400649A (en) | Harmonic responsibility quantification method and device, computer equipment and storage medium | |
CN110504678A (en) | A kind of power signal reconstructing method and system based on prediction matrix | |
CN111832474A (en) | Power signal filtering method and system by using energy scale | |
CN110749841A (en) | Load switch event detection method and system by utilizing conversion space factor | |
CN111639606A (en) | Power signal filtering method and system utilizing Dantzig total gradient minimization | |
CN111830405A (en) | Load switch event detection method and system by using frequency difference | |
CN110783921A (en) | Power signal reconstruction method and system by utilizing dictionary prediction matrix | |
CN112200061A (en) | Power signal reconstruction method and system by using first-order smooth model | |
CN111679119A (en) | Power signal reconstruction method and system by using near-end matching mapping | |
CN110704802A (en) | Power signal reconstruction method and system by utilizing global regularization | |
CN110784019A (en) | Power signal filtering method and system by using Mercer projection |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WW01 | Invention patent application withdrawn after publication | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20200117 |