# A Novel Feature Extraction Method for Nonintrusive Appliance Load Monitoring

Authors: Khaled Chahine; Khalil El Khamlichi Drissi
Journal: Applied Computational Intelligence and Soft Computing
 (2013)

---

## Abstract

Improving energy efficiency by monitoring household electrical consumption is of significant importance with the climate change concerns of the present time. A solution for the electrical consumption management problem is the use of a nonintrusive appliance load monitoring (NIALM) system. This system captures the signals from the aggregate consumption, extracts the features from these signals and classifies the extracted features in order to identify the switched-on appliances. This paper focuses solely on feature extraction through applying the matrix pencil method, a well-known parametric estimation technique, to the drawn electric current. The result is a compact representation of the current signal in terms of complex numbers referred to as poles and residues. These complex numbers are shown to be characteristic of the considered load and can thus serve as features in any subsequent classification module. In the absence of noise, simulations indicate an almost perfect agreement between theoretical and estimated values of poles and residues. For real data, poles and residues are used to determine a feature vector consisting of the contribution of the fundamental, the third, and the fifth harmonic currents to the maximum of the total load current. The result is a three-dimensional feature space with reduced intercluster overlap.


---

## Body

## 1. Introduction

The reason behind the drive for the installation of smart meters in homes and businesses is that they facilitate for consumers to monitor their energy consumption, thereby making it easier for them to save energy, carbon emissions, and money. To help customers as well as utilities in the monitoring process, researchers have been studying load disaggregation schemes for more than two decades.One method of load disaggregation is distributed direct sensing which requires a sensor at each device or appliance in order to measure consumption. The one-sensor-per-device requirement is both the blessing and the curse of this method, for it is highly accurate but expensive. To overcome the limitations associated with the direct sensing approach, researchers have explored methods to infer disaggregated energy usage via a single sensor. Pioneering work in this area is non intrusive appliance load monitoring (NIALM), first introduced by Hart in the late 1980s [1]. In contrast to the direct sensing methods, NIALM relies solely on single-point measurements of voltage and current on the power feed entering the household. NIALM consists of four steps: data acquisition, event detection, feature extraction, and event classification. The raw current and voltage waveforms are transformed into a feature vector, that is, a more compact and meaningful representation that may include real power, reactive power, current-voltage phase difference, and harmonics (e.g., [2]). These extracted features are monitored for changes, identified as events (e.g., an appliance turning “on” or “off”), and classified down to the appliance or device category level using a classification algorithm, which usually compares the features to a preexisting database of signatures. Several reviews of feature extraction methods for electric loads in residential and commercial buildings can be found in the literature [3, 4].Based on the degree of nonintrusiveness, the literature draws a distinction between manual-setup NIALM (MS-NIALM) and automatic-setup NIALM (AS-NIALM) systems. While the former requires switching individual appliances on and off manually to learn their signatures, the latter sets itself up using prior information about potential appliances. AS-NIALM hence extracts the signatures and labels them without any sort of manual intervention which would greatly facilitate mass installation of smart meters. To the authors’ knowledge, no AS-NIALM system has hitherto been implemented. It is hence the main goal of this work to pave the way for such a solution.In this paper, the matrix pencil method, a well-known parametric estimation technique, is applied to the electric current drawn by some elementary linear and nonlinear electric loads driven by a sinusoidal voltage source as well as real loads. The result is a compact representation of the current in terms of complex numbers referred to as poles and residues [5, 6]. These complex numbers are shown to be characteristic of the considered load and thus can serve as features for the subsequent classification phase [7]. For both synthetic and real data, results indicate that poles and residues extracted by the MPM allow an almost perfect reconstruction of drawn electric currents. Results obtained from a database of a household indicate that the extracted features succeed in reducing the intercluster overlap of different appliances.The objectives of this paper are summarized in the following two points:(1)
show that the reduced number of poles and residues estimated by MPM enable an accurate reconstruction of synthetic and real signals;(2)
show that the fundamental and higher harmonic currents determined from poles and residues yield a feature space with reduced intercluster overlap.The rest of the paper is organized as follows. Section2 presents the signal model and the principle of the MPM. Sections 3 and 4 show the validation on simulated and real data, respectively. Finally, Section 5 provides the summary and conclusion.


## 2. Feature Extraction

### 2.1. Signal Model

For a sinusoidal driving voltage of the formv(t)=V2sin(ωt), the drawn electric current can be modeled as a linear combination of d cisoids (complex-valued sinusoidal signals) weighted by complex residues according to the following signal model:
(1)i(t)≈∑m=1d‍rmexp{(αm+j2πfm)t}+b(t),
where rm is the residue of the mth cisoid, αm is its attenuation factor, fm is its frequency, and b(t) is additive white Gaussian noise. After sampling, the time variable, t, is replaced by tk=kts, where ts=6.25×10-4 is the chosen sampling period. The discrete current signal becomes
(2)i(k)≈∑m=1d‍rmzmk+b(k)k=1,2,…,N,
where
(3)zm=exp{(αm+j2πfm)ts}m=1,2,…,d
is the mth complex pole. Under matrix form, the signal model is expressed by
(4)i=Ar+b
with the following notational definitions:
(5)i=[i(1)i(2)…i(N)]T,A=[a1a2…ad],am=[zmzm2…zmN]T,r=[r1r2…rd]T,b=[b(1)b(2)…b(N)]T.
The superscript T denotes the transpose operator.The feature extraction problem can now be stated as follows. Given the electric current data sequence{i(k)}k=1N, use a feature extraction method to extract the complex poles {zm}m=1d and residues {rm}m=1d of the load.


### 2.2. Matrix Pencil Method (MPM)

This section briefly recalls the principle of MPM which is a linear prediction method tailored to the parameter estimation of the damped/undamped exponential model. Starting from the signal model given in (1), MPM chooses a free parameter, L, known as the pencil parameter such as d≤L≤N-d. The proper choice of L results in significant robustness against noise. The next step is to construct a Hankel data matrix:
(6)H=[i(1)i(2)⋯i(L+1)i(2)i(3)⋯i(L+2)⋮⋮⋱⋮i(N-L)i(N-L+1)⋯i(N)].
Two matrices are then obtained by removing the last and first columns of H. In MATLAB notation, they are given as follows:
(7)H→=H(:,1:L),H←=H(:,2:L+1).
The matrix pencil for the two matrices H→ and H← is defined as their linear combination H←-λH→, with λ a scalar parameter. In the absence of noise and owing to the assumed signal model, it is easily verified that H→ and H← admit the following Vandermonde decomposition:
(8)H→=Z1RZ2,H←=Z1RZ0Z2,
where
(9)Z1=[z1z2⋯zdz12z22⋯zd2⋮⋮⋱⋮z1N-Lz2N-L⋯zdN-L],Z2=[1z1⋯z1L-11z2⋯z2L-1⋮⋮⋱⋮1zd⋯zdL-1],Z0=diag{z1,z2,…,zd},R=diag{r1,r2,…,rd},
revealing the fundamental shift-invariance property in the column and row spaces. The matrix pencil can then be written as
(10)H←-λH→=Z1R[Z0-λI]Z2,
where I is the identity matrix. Hence, each value of λ=zm is a rank-reducing number of the pencil. The estimates of zm are, therefore, the generalized eigenvalues (GEVs) of the matrix pair [H←,H⟶].Once the complex poles{zm}m=1d are determined, the complex residues can be estimated using a least squares fit having the following solution:
(11)r=(AHA)-1AHi.For noisy data, total least squares matrix pencil method (TLSMPM) is usually preferred in which the singular value decomposition is used to prefilter the complex signals, and then conventional procedures follow. For more details, the reader can refer to [8].


## 2.1. Signal Model

For a sinusoidal driving voltage of the formv(t)=V2sin(ωt), the drawn electric current can be modeled as a linear combination of d cisoids (complex-valued sinusoidal signals) weighted by complex residues according to the following signal model:
(1)i(t)≈∑m=1d‍rmexp{(αm+j2πfm)t}+b(t),
where rm is the residue of the mth cisoid, αm is its attenuation factor, fm is its frequency, and b(t) is additive white Gaussian noise. After sampling, the time variable, t, is replaced by tk=kts, where ts=6.25×10-4 is the chosen sampling period. The discrete current signal becomes
(2)i(k)≈∑m=1d‍rmzmk+b(k)k=1,2,…,N,
where
(3)zm=exp{(αm+j2πfm)ts}m=1,2,…,d
is the mth complex pole. Under matrix form, the signal model is expressed by
(4)i=Ar+b
with the following notational definitions:
(5)i=[i(1)i(2)…i(N)]T,A=[a1a2…ad],am=[zmzm2…zmN]T,r=[r1r2…rd]T,b=[b(1)b(2)…b(N)]T.
The superscript T denotes the transpose operator.The feature extraction problem can now be stated as follows. Given the electric current data sequence{i(k)}k=1N, use a feature extraction method to extract the complex poles {zm}m=1d and residues {rm}m=1d of the load.


## 2.2. Matrix Pencil Method (MPM)

This section briefly recalls the principle of MPM which is a linear prediction method tailored to the parameter estimation of the damped/undamped exponential model. Starting from the signal model given in (1), MPM chooses a free parameter, L, known as the pencil parameter such as d≤L≤N-d. The proper choice of L results in significant robustness against noise. The next step is to construct a Hankel data matrix:
(6)H=[i(1)i(2)⋯i(L+1)i(2)i(3)⋯i(L+2)⋮⋮⋱⋮i(N-L)i(N-L+1)⋯i(N)].
Two matrices are then obtained by removing the last and first columns of H. In MATLAB notation, they are given as follows:
(7)H→=H(:,1:L),H←=H(:,2:L+1).
The matrix pencil for the two matrices H→ and H← is defined as their linear combination H←-λH→, with λ a scalar parameter. In the absence of noise and owing to the assumed signal model, it is easily verified that H→ and H← admit the following Vandermonde decomposition:
(8)H→=Z1RZ2,H←=Z1RZ0Z2,
where
(9)Z1=[z1z2⋯zdz12z22⋯zd2⋮⋮⋱⋮z1N-Lz2N-L⋯zdN-L],Z2=[1z1⋯z1L-11z2⋯z2L-1⋮⋮⋱⋮1zd⋯zdL-1],Z0=diag{z1,z2,…,zd},R=diag{r1,r2,…,rd},
revealing the fundamental shift-invariance property in the column and row spaces. The matrix pencil can then be written as
(10)H←-λH→=Z1R[Z0-λI]Z2,
where I is the identity matrix. Hence, each value of λ=zm is a rank-reducing number of the pencil. The estimates of zm are, therefore, the generalized eigenvalues (GEVs) of the matrix pair [H←,H⟶].Once the complex poles{zm}m=1d are determined, the complex residues can be estimated using a least squares fit having the following solution:
(11)r=(AHA)-1AHi.For noisy data, total least squares matrix pencil method (TLSMPM) is usually preferred in which the singular value decomposition is used to prefilter the complex signals, and then conventional procedures follow. For more details, the reader can refer to [8].


## 3. Validation on Synthetic Data

### 3.1. Linear Loads

To validate MPM as a feature extraction method, we shall first compare its poles and residues with those obtained from the theoretical expressions of the following linear elementary loads: series RC, series RL, parallel RL, and series RLC. The RC and RL circuits lead to first order differential equations in time whereas the RLC circuit leads to a second-order differential equation. Using Euler’s formula and rearranging allow rewriting the current expression obtained from the solution of the differential equation characterizing the load in the form of (1). The poles and residues of each elementary load can then be readily identified. Tables 1, 2, 3, and 4 give the residues, attenuation factors, and frequencies of the four studied elementary loads. As can be seen from these tables, first-order circuits (RL and RC) are characterized by two pure imaginary conjugate poles representing their forced response and one real pole representing their natural response, whereas the second-order circuit (RLC) has, besides the two pure imaginary conjugate poles of its forced response, two conjugate complex poles related to its natural response. The expressions of the dependent parameters are given in the appendix.Table 1
The residues, attenuation factors, and frequencies of the series RC load.
m
r
m
α
m
f
m (Hz)
1
V
2
2
j
R
cos
(
ϕ
)
e
-
j
ϕ
0
+
50
2
-
V
2
2
j
R
cos
(
ϕ
)
e
j
ϕ
0
-
50
3
-
1
R
(
v
C
0
-
V
2
sin
(
ϕ
)
cos
(
ω
t
0
-
ϕ
)
)
e
t
0
/
τ
-
1
τ
0Table 2
The residues, attenuation factors, and frequencies of the series RL load.
m
r
m
α
m
f
m (Hz)
1
V
2
2
j
R
cos
(
ϕ
)
e
-
j
ϕ
0
+
50
2
-
V
2
2
j
R
cos
(
ϕ
)
e
j
ϕ
0
-
50
3
i
L
0
-
V
2
R
cos
(
ϕ
)
sin
(
ω
t
0
-
ϕ
)
e
t
0
/
τ
-
1
τ
0Table 3
The residues, attenuation factors, and frequencies of the parallel RL load.
m
r
m
α
m
f
m (Hz)
1
V
2
2
j
R
cos
(
ϕ
)
e
-
j
ϕ
0
+
50
2
-
V
2
2
j
R
cos
(
ϕ
)
e
j
ϕ
0
-
50
3
i
L
0
-
V
2
R
cos
(
ϕ
)
sin
(
ω
t
0
-
ϕ
)
0
0Table 4
The residues, attenuation factors, and frequencies of the series RLC load.
m
r
m
α
m
f
m (Hz)
1
V
2
2
j
R
cos
(
ϕ
)
e
-
j
ϕ
0
+
50
2
-
V
2
2
j
R
cos
(
ϕ
)
e
j
ϕ
0
-
50
3
A
e
-
k
1
ω
0
t
0
k
1
ω
0
0
4
B
e
-
k
2
ω
0
t
0
k
2
ω
0
0


### 3.2. Nonlinear Loads

A nonlinear load is one for which the relationship between the current through the load and the voltage across the load is a nonlinear function. A simple view of the nature of nonlinear loads can be presented using Ohm’s Law, which states that the voltage is the product of the load resistance and the current (V=RI). For a linear load, the resistance (R) is a constant; for a nonlinear load, the resistance varies. When AC power is supplied to a nonlinear load, the result is the creation of currents that do not oscillate at the supply frequency. These currents are called harmonics. Harmonics occur at multiples of the supply (fundamental) frequency. For instance, if the fundamental frequency is 50 Hz, the so-called second harmonic is 100 Hz, the third harmonic is 150 Hz, and so on. Any number of harmonics can be created by a particular piece of equipment depending on that equipment’s electrical characteristics. Therefore, the current drawn by nonlinear loads can still be represented by (1) where harmonics appear in the form of pole-residue couples at frequency multiples of 50 Hz.


### 3.3. Results

Assuming zero initial conditions (iL0=0 and/or vC0=0), the following numerical values were used to determine the electric current data sequence from which MPM extracted poles and residues: {R=100 Ω, C=0.1 mF} for the series RC circuit, {R=10 Ω, L=100mH} for both the series and parallel RL circuits, and {R=1 Ω, L=20mH, C=60 mF} for the series RLC circuit. A duration of ten periods or 0.2 seconds was chosen for the current which at ts=6.25×10-4 is equivalent to 320 samples, and MPM was applied at each period. Figures 1, 2, 3, and 4 show the current obtained from the analytic expression of poles and residues in the tables above and its reconstruction obtained from the poles and residues extracted by MPM. An almost perfect agreement can be seen between the two curves indicating the accuracy of the characteristic complex numbers extracted by MPM. In addition, the figures show the forced and natural responses of each of the four elementary circuits.Figure 1
The analytic and reconstructed currents of the series RC circuit along with its forced and natural responses.Figure 2
The analytic and reconstructed currents of the series RL circuit along with its forced and natural responses.Figure 3
The analytic and reconstructed currents of the parallel RL circuit along with its forced and natural responses.Figure 4
The analytic and reconstructed currents of the series RLC circuit along with its forced and natural responses.To evaluate the performance of MPM on nonlinear loads, we considered the current shown in Table5. It consists of a fundamental and four harmonics and hence can be represented by ten pairwise complex conjugate pole-residue couples. We then used MPM to extract these ten couples which served to reconstruct the current as shown in Figure 5. As can be seen, MPM is successful in estimating the pole-residue couples of the load.Table 5
Current composition of a nonlinear load.
I
1
I
5
I
7
I
11
I
13
100%
18.9%
11%
5.9%
4.8%Figure 5
The analytic and reconstructed currents of a nonlinear load. The inset zooms in on a half period of the drawn current in order to show the accuracy of reconstruction.


## 3.1. Linear Loads

To validate MPM as a feature extraction method, we shall first compare its poles and residues with those obtained from the theoretical expressions of the following linear elementary loads: series RC, series RL, parallel RL, and series RLC. The RC and RL circuits lead to first order differential equations in time whereas the RLC circuit leads to a second-order differential equation. Using Euler’s formula and rearranging allow rewriting the current expression obtained from the solution of the differential equation characterizing the load in the form of (1). The poles and residues of each elementary load can then be readily identified. Tables 1, 2, 3, and 4 give the residues, attenuation factors, and frequencies of the four studied elementary loads. As can be seen from these tables, first-order circuits (RL and RC) are characterized by two pure imaginary conjugate poles representing their forced response and one real pole representing their natural response, whereas the second-order circuit (RLC) has, besides the two pure imaginary conjugate poles of its forced response, two conjugate complex poles related to its natural response. The expressions of the dependent parameters are given in the appendix.Table 1
The residues, attenuation factors, and frequencies of the series RC load.
m
r
m
α
m
f
m (Hz)
1
V
2
2
j
R
cos
(
ϕ
)
e
-
j
ϕ
0
+
50
2
-
V
2
2
j
R
cos
(
ϕ
)
e
j
ϕ
0
-
50
3
-
1
R
(
v
C
0
-
V
2
sin
(
ϕ
)
cos
(
ω
t
0
-
ϕ
)
)
e
t
0
/
τ
-
1
τ
0Table 2
The residues, attenuation factors, and frequencies of the series RL load.
m
r
m
α
m
f
m (Hz)
1
V
2
2
j
R
cos
(
ϕ
)
e
-
j
ϕ
0
+
50
2
-
V
2
2
j
R
cos
(
ϕ
)
e
j
ϕ
0
-
50
3
i
L
0
-
V
2
R
cos
(
ϕ
)
sin
(
ω
t
0
-
ϕ
)
e
t
0
/
τ
-
1
τ
0Table 3
The residues, attenuation factors, and frequencies of the parallel RL load.
m
r
m
α
m
f
m (Hz)
1
V
2
2
j
R
cos
(
ϕ
)
e
-
j
ϕ
0
+
50
2
-
V
2
2
j
R
cos
(
ϕ
)
e
j
ϕ
0
-
50
3
i
L
0
-
V
2
R
cos
(
ϕ
)
sin
(
ω
t
0
-
ϕ
)
0
0Table 4
The residues, attenuation factors, and frequencies of the series RLC load.
m
r
m
α
m
f
m (Hz)
1
V
2
2
j
R
cos
(
ϕ
)
e
-
j
ϕ
0
+
50
2
-
V
2
2
j
R
cos
(
ϕ
)
e
j
ϕ
0
-
50
3
A
e
-
k
1
ω
0
t
0
k
1
ω
0
0
4
B
e
-
k
2
ω
0
t
0
k
2
ω
0
0


## 3.2. Nonlinear Loads

A nonlinear load is one for which the relationship between the current through the load and the voltage across the load is a nonlinear function. A simple view of the nature of nonlinear loads can be presented using Ohm’s Law, which states that the voltage is the product of the load resistance and the current (V=RI). For a linear load, the resistance (R) is a constant; for a nonlinear load, the resistance varies. When AC power is supplied to a nonlinear load, the result is the creation of currents that do not oscillate at the supply frequency. These currents are called harmonics. Harmonics occur at multiples of the supply (fundamental) frequency. For instance, if the fundamental frequency is 50 Hz, the so-called second harmonic is 100 Hz, the third harmonic is 150 Hz, and so on. Any number of harmonics can be created by a particular piece of equipment depending on that equipment’s electrical characteristics. Therefore, the current drawn by nonlinear loads can still be represented by (1) where harmonics appear in the form of pole-residue couples at frequency multiples of 50 Hz.


## 3.3. Results

Assuming zero initial conditions (iL0=0 and/or vC0=0), the following numerical values were used to determine the electric current data sequence from which MPM extracted poles and residues: {R=100 Ω, C=0.1 mF} for the series RC circuit, {R=10 Ω, L=100mH} for both the series and parallel RL circuits, and {R=1 Ω, L=20mH, C=60 mF} for the series RLC circuit. A duration of ten periods or 0.2 seconds was chosen for the current which at ts=6.25×10-4 is equivalent to 320 samples, and MPM was applied at each period. Figures 1, 2, 3, and 4 show the current obtained from the analytic expression of poles and residues in the tables above and its reconstruction obtained from the poles and residues extracted by MPM. An almost perfect agreement can be seen between the two curves indicating the accuracy of the characteristic complex numbers extracted by MPM. In addition, the figures show the forced and natural responses of each of the four elementary circuits.Figure 1
The analytic and reconstructed currents of the series RC circuit along with its forced and natural responses.Figure 2
The analytic and reconstructed currents of the series RL circuit along with its forced and natural responses.Figure 3
The analytic and reconstructed currents of the parallel RL circuit along with its forced and natural responses.Figure 4
The analytic and reconstructed currents of the series RLC circuit along with its forced and natural responses.To evaluate the performance of MPM on nonlinear loads, we considered the current shown in Table5. It consists of a fundamental and four harmonics and hence can be represented by ten pairwise complex conjugate pole-residue couples. We then used MPM to extract these ten couples which served to reconstruct the current as shown in Figure 5. As can be seen, MPM is successful in estimating the pole-residue couples of the load.Table 5
Current composition of a nonlinear load.
I
1
I
5
I
7
I
11
I
13
100%
18.9%
11%
5.9%
4.8%Figure 5
The analytic and reconstructed currents of a nonlinear load. The inset zooms in on a half period of the drawn current in order to show the accuracy of reconstruction.


## 4. Validation on Real Data

### 4.1. Reconstruction Results

In this section, the validation of MPM is carried out on currents of three representative loads: a television set, a vacuum cleaner, and an economy lamp. As for the case of synthetic data, MPM was applied at each period. Figures6, 7, and 8 show the current drawn by the appliances and its reconstruction based on the pole-residue estimates of MPM. The close agreement shown in the figures indicate that the exponential model of (1) and its parameters estimated by the MPM accurately predict the response of the actual loads. It is worth mentioning that the number of pole-residue couples d increases with the nonlinearity of the load. For instance, the current of the vacuum cleaner could be accurately reconstructed from four pole-residue couples, whereas that of the economy lamp needed up to twelve couples.Figure 6
The measured and reconstructed currents of the television set.Figure 7
The measured and reconstructed currents of the vacuum cleaner.Figure 8
The measured and reconstructed currents of the economy lamp.


### 4.2. Feature Space

The feature space contains 900 signatures uniformly distributed among the following nine appliances: incandescent lamp, halogen lamp, economy lamp, water heater, electric convector, oven, two-burner hot plate, television set, and computer. As shown in Figure9, each signature (represented by a point in the the three-dimensional feature space) is characterized by three pole-residue products corresponding to the maxima of the fundamental, third, and fifth harmonic currents. The restriction to three frequencies has the sole aim of representing the feature space graphically. From the feature space, ten clusters representing the nine appliances can be clearly distinguished. The additional cluster is due to the two-burner hot plate which is represented by two clusters, one for each burner. It can hence be concluded that the studied appliances can be fairly distinguished using the fundamental and higher harmonics.Figure 9
The feature space showing the disaggregated contribution of the fundamental and harmonic currents to the maximum of the total current for several appliances.


## 4.1. Reconstruction Results

In this section, the validation of MPM is carried out on currents of three representative loads: a television set, a vacuum cleaner, and an economy lamp. As for the case of synthetic data, MPM was applied at each period. Figures6, 7, and 8 show the current drawn by the appliances and its reconstruction based on the pole-residue estimates of MPM. The close agreement shown in the figures indicate that the exponential model of (1) and its parameters estimated by the MPM accurately predict the response of the actual loads. It is worth mentioning that the number of pole-residue couples d increases with the nonlinearity of the load. For instance, the current of the vacuum cleaner could be accurately reconstructed from four pole-residue couples, whereas that of the economy lamp needed up to twelve couples.Figure 6
The measured and reconstructed currents of the television set.Figure 7
The measured and reconstructed currents of the vacuum cleaner.Figure 8
The measured and reconstructed currents of the economy lamp.


## 4.2. Feature Space

The feature space contains 900 signatures uniformly distributed among the following nine appliances: incandescent lamp, halogen lamp, economy lamp, water heater, electric convector, oven, two-burner hot plate, television set, and computer. As shown in Figure9, each signature (represented by a point in the the three-dimensional feature space) is characterized by three pole-residue products corresponding to the maxima of the fundamental, third, and fifth harmonic currents. The restriction to three frequencies has the sole aim of representing the feature space graphically. From the feature space, ten clusters representing the nine appliances can be clearly distinguished. The additional cluster is due to the two-burner hot plate which is represented by two clusters, one for each burner. It can hence be concluded that the studied appliances can be fairly distinguished using the fundamental and higher harmonics.Figure 9
The feature space showing the disaggregated contribution of the fundamental and harmonic currents to the maximum of the total current for several appliances.


## 5. Conclusion

This paper presented a novel feature extraction method for non intrusive appliance load monitoring. First, the poles and residues estimated by the matrix pencil method were shown to enable accurate reconstruction of synthetic and real current signals. Second, these complex numbers were used to determine a three-dimensional feature space with reduced intercluster overlap. Future research will make use of the extracted features for the classification phase.



---

*Source: 686345-2013-05-02.xml*