CN110531149B - Power signal filtering method and system based on waveform regularization - Google Patents
Power signal filtering method and system based on waveform regularization Download PDFInfo
- Publication number
- CN110531149B CN110531149B CN201910818165.2A CN201910818165A CN110531149B CN 110531149 B CN110531149 B CN 110531149B CN 201910818165 A CN201910818165 A CN 201910818165A CN 110531149 B CN110531149 B CN 110531149B
- Authority
- CN
- China
- Prior art keywords
- power signal
- opt
- following
- steps
- signal sequence
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000001914 filtration Methods 0.000 title claims abstract description 28
- 239000011159 matrix material Substances 0.000 claims abstract description 34
- 108010076504 Protein Sorting Signals Proteins 0.000 claims abstract description 26
- 239000013598 vector Substances 0.000 claims abstract description 23
- 238000012937 correction Methods 0.000 claims abstract description 11
- 238000012545 processing Methods 0.000 claims abstract description 6
- 238000004364 calculation method Methods 0.000 description 9
- 238000001514 detection method Methods 0.000 description 9
- 238000000354 decomposition reaction Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 3
- 230000017105 transposition Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000005236 sound signal Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R21/00—Arrangements for measuring electric power or power factor
- G01R21/001—Measuring real or reactive component; Measuring apparent energy
- G01R21/002—Measuring real component
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/01—Arrangements for reducing harmonics or ripples
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/40—Arrangements for reducing harmonics
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Indication And Recording Devices For Special Purposes And Tariff Metering Devices (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
The embodiment of the invention discloses a power signal filtering method and a system based on waveform regularization, wherein the method comprises the following steps: step 1, inputting an actually measured power signal sequence S; step 2, carrying out noise filtering processing on the power signal sequence S, wherein the power signal sequence after noise filtering is SNEW. The method specifically comprises the following steps: sNEW=mOPT+B(S‑LmOPT) (ii) a Wherein m isOPTIs the best prediction vector; b is a correction matrix; l denotes a system matrix.
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.
Therefore, in the switching event detection process, the actually measured power signal used is often affected by noise, and the switching event detection cannot be performed correctly by using the imperfect power signal. Therefore, how to effectively reconstruct the incomplete power signal and filter the influence of noise is the key to the success of the 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 filtering method and a system based on waveform regularization. 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 filtering a power signal based on waveform regularization, comprising:
step 1, inputting an actually measured power signal sequence S;
step 2, carrying out noise filtering processing on the power signal sequence S, wherein the power signal sequence after noise filtering is SNEW. The method specifically comprises the following steps: sNEW=mOPT+B(S-LmOPT) (ii) a Wherein m isOPTIs the best prediction vector; b is a correction matrix; l denotes a system matrix.
A power signal filtering system based on waveform regularization, comprising:
the acquisition module inputs an actually measured power signal sequence S;
the filtering module is used for carrying out noise filtering processing on the power signal sequence S, and the power signal sequence after noise filtering is SNEW. The method specifically comprises the following steps: sNEW=mOPT+B(S-LmOPT) (ii) a Wherein m isOPTIs the best prediction vector; b is a correction matrix; l denotes a system matrix.
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 and relatively mature technology in non-invasive load decomposition, the power signal is often submerged in the pulse noise with strong amplitude during the acquisition and transmission process, and the switching event detection cannot be correctly performed by using the imperfect power signal. Therefore, how to effectively reconstruct the incomplete power signal and filter the influence of noise is the key to the success of the 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 filtering method and a system based on waveform regularization. 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 filtering method based on waveform regularization
Fig. 1 is a schematic flow chart of a power signal filtering method based on waveform regularization according to the present invention. As shown in fig. 1, the method for filtering a power signal based on waveform regularization specifically includes the following steps:
step 1, inputting an actually measured power signal sequence S;
step 2, carrying out noise filtering processing on the power signal sequence S, wherein the power signal sequence after noise filtering is SNEW. The method specifically comprises the following steps: sNEW=mOPT+B(S-LmOPT) (ii) a Wherein m isOPTIs the best prediction vector; b is a correction matrix; l denotes a system matrix.
Before the step 2, the method further comprises:
step 3, obtaining the optimal prediction vector mOPTA correction matrix B and a system matrix L.
The step 3 comprises the following steps:
step 301, generating a delay vector SDThe method specifically comprises the following steps:
SD=[sK+1,sK+2,…,sN,s1,s2,…,sK]
wherein
N: length of the signal sequence S
SNR: signal-to-noise ratio of the signal sequence S
Step 302, obtaining a signal correlation matrix C, specifically:
C=STSD
*T: transposing representation pairs
Step (ii) of303, obtaining a compressed vector SQThe method specifically comprises the following steps:
step 304, obtaining a compressed incidence matrix CQThe method specifically comprises the following steps:
CQ=STSQ
*T: transposing representation pairs
Step 305, obtaining the correction matrix B, specifically:
B=C-TCQ
*-T: inverse matrix transposition of pairs
Step 306, obtaining the system matrix L, specifically:
L=CQC
step 307, iteratively calculating the optimal prediction vector mOPTThe method specifically comprises the following steps:
the first step is as follows: initializing specifically as follows: m is1And k is 1 and is an iteration control parameter.
The second step is that: updating, specifically:
mk+1=mk+B[S-Lmk]。
the third step: and (5) adding 1 to the iteration control parameter k, and repeating the second step. Until the difference between the two iteration results is less than 0.001, obtaining the optimal prediction vector mOPT=mK,mKIs the result of the last update.
FIG. 2 is a structural intent of a power signal filtering system based on waveform regularization
Fig. 2 is a schematic structural diagram of a power signal filtering system based on waveform regularization according to the present invention. As shown in fig. 2, the power signal filtering system based on waveform regularization includes the following structures:
the acquisition module 401 inputs an actually measured power signal sequence S;
a reconstruction module 402 for performing noise filtering processing on the power signal sequence SThe power signal sequence after noise filtering is SNEW. The method specifically comprises the following steps: sNEW=mOPT+B(S-LmOPT) (ii) a Wherein m isOPTIs the best prediction vector; b is a correction matrix; l denotes a system matrix.
The system further comprises:
a calculating module 403 for obtaining the optimal prediction vector mOPTA correction matrix B and a system matrix L.
The calculation module 403 includes the following units:
delay section 4301 generates delay vector SDThe method specifically comprises the following steps:
SD=[sK+1,sK+2,…,sN,s1,s2,…,sK]
wherein
N: length of the signal sequence S
SNR: signal-to-noise ratio of the signal sequence S
The first calculation unit 4302 obtains a signal correlation matrix C, which specifically is:
C=STSD
*T: transposing representation pairs
Second calculation section 4303 obtains compressed vector SQThe method specifically comprises the following steps:
third calculation unit 4304, find compressed correlation matrix CQThe method specifically comprises the following steps:
CQ=STSQ
*T: transposing representation pairs
The fourth calculation unit 4305 obtains the correction matrix B, specifically:
B=C-TCQ
*-T: inverse matrix transposition of pairs
The fifth calculation unit 4306 obtains the system matrix L, which specifically is:
L=CQC
an iteration unit 4307 for iteratively calculating the optimal prediction vector mOPTThe method specifically comprises the following steps:
the first step is as follows: initializing specifically as follows: m is1And (5) setting k to 1 as an iteration control parameter.
The second step is that: updating, specifically:
mk+1=mk+B[S-Lmk]。
the third step: and (5) adding 1 to the iteration control parameter k, and repeating the second step. Until the difference between the two iteration results is less than 0.001, obtaining the optimal prediction vector mOPT=mK,mKIs the result of the last update.
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. Generating a delay vector SD
SD=[sK+1,sK+2,…,sN,s1,s2,…,sK]
Wherein
N: length of the signal sequence S
SNR: signal-to-noise ratio of the signal sequence S
3. Determining a signal correlation matrix
C=STSD
*T: the representation transposes.
4. Finding compressed vectors
5. Determining a compressed correlation matrix
CQ=STSQ
*T: transposing representation pairs
6. Obtaining the correction matrix
B=C-TCQ
*-T: inverse matrix transposition of pairs
7. Determining a system matrix
L=CQC
8. Iterative calculation of optimal prediction vector
The first step is as follows: initializing specifically as follows: m is1And (5) setting k to 1 as an iteration control parameter.
The second step is that: updating, specifically:
mk+1=mk+B[S-Lmk]。
the third step: and (5) adding 1 to the iteration control parameter k, and repeating the second step. Until the difference between the two iteration results is less than 0.001, obtaining the optimal prediction vector mOPT=mK,mKIs the result of the last update.
9. Filtering
SNEW=mOPT+B(S-LmOPT)
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 (1)
1. A power signal filtering method based on waveform regularization is characterized by comprising the following steps:
step 1, inputting an actually measured power signal sequence S;
step 2, generating a delay vector SDThe method specifically comprises the following steps:
SD=[sK+1,sK+2,…,sN,s1,s2,…,sK];
wherein:
n: the length of the signal sequence S;
SNR: the signal-to-noise ratio of the signal sequence S;
step 3, solving a signal correlation matrix C, specifically:
C=STSD;
*T: representing transposing of pairs;
step 4, obtaining a compressed vector SQThe method specifically comprises the following steps:
step 5, solving a compressed incidence matrix CQThe method specifically comprises the following steps:
CQ=STSQ;
*T: representing transposing of pairs;
step 6, obtaining a correction matrix B, specifically:
B=C-TCQ;
*-T: the inversion of the inverse matrix performed on the x is expressed;
step 7, solving a system matrix L, specifically:
L=CQC;
step 8, iteratively solving the optimal prediction vector mOPTThe method specifically comprises the following steps:
the first step is as follows: initializing specifically as follows: m is1S, k 1, which is an iteration control parameter;
the second step is that: updating, specifically:
mk+1=mk+B[S-Lmk];
the third step: adding 1 to the iteration control parameter k, repeating the second step until the difference between the two iteration results is less than 0.001 to obtain the optimal prediction vector mOPT=mK,mKIs the result of the last update;
step 9, carrying out noise filtering processing on the power signal sequence S, wherein the power signal sequence after noise filtering is SNEW(ii) a The method specifically comprises the following steps: sNEW=mOPT+B(S-LmOPT)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910818165.2A CN110531149B (en) | 2019-08-31 | 2019-08-31 | Power signal filtering method and system based on waveform regularization |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910818165.2A CN110531149B (en) | 2019-08-31 | 2019-08-31 | Power signal filtering method and system based on waveform regularization |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110531149A CN110531149A (en) | 2019-12-03 |
CN110531149B true CN110531149B (en) | 2021-06-18 |
Family
ID=68666003
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910818165.2A Expired - Fee Related CN110531149B (en) | 2019-08-31 | 2019-08-31 | Power signal filtering method and system based on waveform regularization |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110531149B (en) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1969459A (en) * | 2004-03-01 | 2007-05-23 | 电力波技术公司 | Digital predistortion system and method for linearizing an RF power amplifier with nonlinear gain characteristics and memory effects |
CN107749304A (en) * | 2017-09-07 | 2018-03-02 | 电信科学技术研究院 | The sustainable renewal method and device of finite impulse response filter coefficient vector |
CN108918932A (en) * | 2018-09-11 | 2018-11-30 | 广东石油化工学院 | Power signal adaptive filter method in load decomposition |
CN109120242A (en) * | 2018-08-24 | 2019-01-01 | 广东石油化工学院 | A kind of coherent noise adaptive filter method and system |
CN109307798A (en) * | 2018-08-29 | 2019-02-05 | 广东石油化工学院 | A kind of power signal filtering method for switch events detection |
CN109391195A (en) * | 2017-08-08 | 2019-02-26 | 西门子股份公司 | System mode prediction |
CN109727272A (en) * | 2018-11-20 | 2019-05-07 | 南京邮电大学 | A kind of method for tracking target based on double branch's space-time regularization correlation filters |
-
2019
- 2019-08-31 CN CN201910818165.2A patent/CN110531149B/en not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1969459A (en) * | 2004-03-01 | 2007-05-23 | 电力波技术公司 | Digital predistortion system and method for linearizing an RF power amplifier with nonlinear gain characteristics and memory effects |
CN109391195A (en) * | 2017-08-08 | 2019-02-26 | 西门子股份公司 | System mode prediction |
CN107749304A (en) * | 2017-09-07 | 2018-03-02 | 电信科学技术研究院 | The sustainable renewal method and device of finite impulse response filter coefficient vector |
CN109120242A (en) * | 2018-08-24 | 2019-01-01 | 广东石油化工学院 | A kind of coherent noise adaptive filter method and system |
CN109307798A (en) * | 2018-08-29 | 2019-02-05 | 广东石油化工学院 | A kind of power signal filtering method for switch events detection |
CN108918932A (en) * | 2018-09-11 | 2018-11-30 | 广东石油化工学院 | Power signal adaptive filter method in load decomposition |
CN109727272A (en) * | 2018-11-20 | 2019-05-07 | 南京邮电大学 | A kind of method for tracking target based on double branch's space-time regularization correlation filters |
Also Published As
Publication number | Publication date |
---|---|
CN110531149A (en) | 2019-12-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108918931B (en) | Power signal self-adaptive filtering method in load decomposition | |
CN108918930B (en) | Power signal self-adaptive reconstruction method in load decomposition | |
CN108918932A (en) | Power signal adaptive filter method in load decomposition | |
CN108918928B (en) | Power signal self-adaptive reconstruction method in load decomposition | |
CN111680590A (en) | Power signal filtering method and system by using contraction gradient | |
CN109063676B (en) | Self-adaptive time-frequency peak value filtering method and system for power signal | |
CN110531149B (en) | Power signal filtering method and system based on waveform regularization | |
CN111666870A (en) | Power signal reconstruction method and system by utilizing quadratic constraint | |
CN108918929B (en) | Power signal self-adaptive filtering method in load decomposition | |
CN110704800A (en) | Power signal reconstruction method and system based on expansion coefficient | |
CN110542855B (en) | Load switch event detection method and system based on discrete cosine transform | |
CN112434567B (en) | Power signal filtering method and system by using noise jitter property | |
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 | |
CN109241874A (en) | Power signal filtering method in Energy Decomposition | |
CN110221119B (en) | Load switch event detection method and system based on power and akie fusion information | |
CN111737645A (en) | Power signal reconstruction method and system by using prediction matrix | |
CN111639606A (en) | Power signal filtering method and system utilizing Dantzig total gradient minimization | |
CN110705426B (en) | Power signal filtering method and system by using deblurring operator | |
CN112307997B (en) | Power signal reconstruction method and system by using main mode decomposition | |
CN110704795A (en) | Power signal reconstruction method and system by using shrinkage operator | |
CN110611310A (en) | Power signal reconstruction method and system by using periodic boundary conditions | |
CN110504678A (en) | A kind of power signal reconstructing method and system based on prediction matrix | |
CN110542790A (en) | power signal filtering method and system based on sparse model | |
CN110783921A (en) | Power signal reconstruction method and system by utilizing dictionary prediction matrix |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210618 |
|
CF01 | Termination of patent right due to non-payment of annual fee |