Abstract
This paper aims to develop a new idea for the classification of power quality disturbances. The method is based on Stockwell’s-Transform (ST) and Type-2 Fuzzy Kernel Support Vector Machine (T2FK-SVM). Through the introduction of ST and its properties, we propose a classification plan for nine types of power quality disturbances. Firstly, features of disturbance signals extracted through the ST. Secondly, features extracted by using the ST are applied as input to T2FK-SVM classifier for automatic classification of the power quality (PQ) disturbances. Design of Kernel is a main part of many kernel based methods such as Support Vector Machine (SVM), so by using of Type-2 Fuzzy sets as a kernel of SVM, the total accuracy of classification enhanced.This method can reduce the features of the disturbance signals significantly, and so less time and memory is required for classification by the T2FK-SVM method. Six single event and two complex event as well normal voltage selected as reference are considered for the classification. The simulation results showed accurate classification, fast learning and execution in the detection and classification of PQ events. Results are compared with other methods and the robustness of proposed method evaluated under noisy conditions. Finally, proposed method is also implemented on real time PQ disturbances to confirm the validity of this method in practical conditions.
Keywords
Introduction
Power supply quality issues and the resulting problems are the consequences of increasing use of solid state switching devices, non-linear and power electronically switched loads, unbalanced power systems, lighting controls, computer and data processing equipment as well as industrial plant rectifiers and inverters. These electronic type loads cause voltage distortion, inrush, pulse type current phenomenon with excessive harmonics and high current distortion [1].
The problem of power quality event classification has three key parts: processing method, feature selection and classification. The goal of processing method is to convert raw data into feature extractable data. The second part of the PQ event classification is optimal feature selection in this segment, a suitable set of feature must be selected from converted signal. These sets of features must have smaller size to reduce classification time. Classification process is the last part of PQ event classification. In this part extracted features must be classified into the accurate clusters.
In literature, many signal processing methods are introduced for analyzing PQ events and extracting features: Short time Fourier transform (STFT) method [2], discrete wavelet transform (DWT) method [3], empirical mode decomposition (EMD) method [4], Hilbert Huang transform (HHT) method [5] and S-transform (ST) method [6]
Artificial intelligence (AI) based methods can be classify PQ events accurately and efficiently. Many AI based methods are applied for PQ events classification: support vector machine (SVM) method [7, 8], artificial neural network (ANN) system [9], fuzzy-expert system [10, 11], neuro-fuzzy system [12], genetic algorithm (GA) method [13] and particle swarm optimization (PSO) method [14].
PQ events are non-stationary signals, features extraction is an important for event classification. STFT, DWT and ST are popular methods for feature extraction that used in many recent works. PQ event detection by ST used in [15] and a new method based on the binary feature matrix is designed for making a decision regarding the disturbance type. Detection and classification of voltage disturbances using a Fuzzy-ARTMAP-wavelet network introduced in [16] in this research wavelet extract features and Fuzzy-ARTMAP neural network classified extracted features. Fuzzy C-means and adaptive particle swarm optimization for classification of PQ events proposed in [17].In this article authors used windowing techniques for non-stationary disturbance signals. Balanced neural tree proposed by Authors of [18] for PQ event classification.The authors used Hilbert transform and empirical mode decomposition for feature extraction.
This article is described as follows: S-transform is described in section 2, and section 3 addresses Type-2 Fuzzy Kernel Based Support Vector Machine method. Problem formulation, simulation results and conclusion are given in section 4-6 to show the accuracy and efficiency of the proposed method.
S-transform
The S-transform [19] based on two-advanced signal processing tools; the short-time Fourier transforms (STFT) and the wavelet transform. It can be viewed as a frequency dependent STFT or a phase corrected wavelet transform [19]. The frequency dependent window used for analysis of a signal so the S-transform has better performance than other time-frequency transforms. Furthermore, it provides superior time-frequency localization property computing both amplitude and phase spectrum of discrete data samples [19].The S-transform of a signal h(t) is defined as:
Where
The inverse S-transform is given by
Another representation of S-transform is an amplitude and phase correction of the CWT (continuous wavelet transform) as follows
Equation (6) shows distribution of the time-frequency resolution in the time frequency plane like wavelet transform but with Fourier transform a direct link exists. So the S-transform is given as follows
A discrete explanation of the S-transform is obtained from (8) as follows
The H[n] is the DFT of h (k) and can be calculated by FFT.
Where
Support vector machines
SVM is a powerful and reliable method for classification of patterns [20, 21]. Given the training data (x1, y1) , . . . , (x l , y l ) , x ∈ R M , y i ∈ {-1, + 1} for two class problem, the decision functions of form sgn ((w T x i ) + w0) constructs by SVM with the maximum margin, where w is the normal vector of the separating hyperplane in the canonical form and w0 is a bias term [21]. The distances of the point closest to the hyperplanes of both -1 and+1 are calculated as 1/∥ w ∥/∥ w ∥. The separating margin is defined to be 2/||w||.
In many real world applications, data is noisy and a linear hyperplane cannot separate. To allow deteriorating from margin, ξ
i
≥ 1 are introduced as follows:
For the training data x1 if 0 ≤ ξ i < 1, the data do not have the maximum margin but are still correctly classified. But if ξ i ≥ 1, the data are misclassified. So by leaving intramargin the noise points occurring to near the boundaries region the separation margin is improved and that generalization performance is increased [21].
SVM builds the constraint primal quadratic optimization problem that minimizes the training and generalization error as follows:
SVM maps the inputs x into some higher dimensional space by means of a nonlinear feature mapping φ (x) for solving the classification problem separated by only highly complex decision boundaries in the input space. Thus the problem changes into linearly separable case at the feature space. If only scalar product
K (x, x
i
) = x
T
x
i
in (8) is a linear kernel. One of the most common kernels is Gaussian kernel defined as:
Multi-class SVM classifiers carries out by combining two class SVMs. The most efficient methods are 1-against-r scheme (with r represent the number of classes) and the 1-against-1 scheme [23]. 1-against-r scheme users (r1) two classifiers: each machine is trained as a classifier for one class against all other classes. 1-against-1 scheme creates a multi-class classifier, i.e. (r (r1)/2) two class machines are created. Each machine is trained as a classifier for one class against other classes. In order to classify test data, pairwise competition between all the machines is performed; each winner competes against another winner until a single winner remains. This final winner determines the class of the test data [23].
Performance of SVM reduced in the uncertain condition and the optimal design of kernel function can enhance accuracy and performance. In this research interval type-2 fuzzy function used as kernel of SVM classifier. This kernel can efficiently classify power quality events in uncertain conditions.
An interval type-2 fuzzy set
J
x
is called primary membership of
The upper membership function (UMF) and lower membership function (LMF) of
Interval type-2 fuzzy set concept applied, to represent the uncertainty of the kernel. The primary membership J
k
m
of a kernel K
m
can be represented as a membership interval with all secondary grades of the primary memberships equaling to one. An interval type-2 fuzzy kernel set can be represented as [24]:
RBF kernel set with different kernel parameters would be extended to interval type-2 fuzzy kernel type as shown in Fig. 1 [24].

Example of a an interval type-2 fuzzy kernel.
In this part, an interval of primary membership of kernels computed with various kernel parameters which show different degrees of membership. A kernel set to interval type-2 fuzzy kernel sets to show uncertainty by using this kernel in SVM. The centroid type-reducer used for the type-reduction of interval type-2 fuzzy kernels. Considering all the M kernels, the fuzzy relation of interval type-2 fuzzy kernels are given as follows:
The upper output and lower output of the kernel defined as follows:
In the defuzzification process final output of the kernel computed for the learning part of SVM, which is performed by defuzzification operation on the interval set
By using interval type-2 fuzzy kernels, the robustness and reliability of SVM improved.
The proposed method for PQ event classification have three key parts: preprocessing, feature extraction and classification. Preprocessing part used the equation of each event to generate event waveform. ST applied for feature extraction and T2FK-SVM classified features. The diagram of proposed method presented in Fig. 2.

Diagram of proposed method.
In this paper nine types of events considered as follows: normal voltage, voltage sag, swell, interruption, oscillatory transient, flicker, sag with harmonic and swell with harmonic.
Feature extraction
The samples of PQ event are huge for classification process so valuable features must be extracted from events. In this part features of each signal extracted by ST. Five features extracted from the ST output of each event as follows:
F1 = max (A) where A is the amplitude versus time graph from the ST output matrix under disturbance.
F2 = min (A) where A is the amplitude versus time graph from the ST output matrix under disturbance.
where B is the amplitude versus time graph of the ST output in normal condition.
F4= total harmonic distortion (THD) =
F5 = estimated frequency at the maximum amplitude obtained from ST output matrix. So a single vector with five elements is generated for each PQ event and fed to classifier stage. The number of features that fed to classifier is low (only 5 inputs). The ST output contour for three PQ events showed in Figs. 3–5.

(a) Voltage sag signal and it’s ST, (b) Magnitude– time and amplitude– frequency diagram from ST.

(a) Voltage oscillatory transient signal and it’s ST, (b) Magnitude– time and amplitude– frequency diagram from ST.

(a) Voltage swell with harmonics signal and it’s ST, (b) Magnitude– time and amplitude– frequency diagram from ST.
In this stage features of events classified by T2FK-SVM. The Kernel function of T2FK-SVM is interval type-2 fuzzy function. Firstly, parameters of T2FK-SVM set by trial and error. Secondly SVM trained by using optimized parameters in the previous step and 50% of events that fed to the inputs of T2FK-SVM. Finally performance of T2FK-SVM evaluated by 50% of generated events.
Simulation and analysis
Data generation
Signal generation by parametric equations for event classification has many advantages. It was possible to change training and testing signal parameters in a wide range and in a controlled manner. The signals simulated with this way were very close to the real world conditions. On the other hand, different signals belonging to the same class gave possibility to estimate the generalization ability of classifiers based on neural networks [25].
Simulation results
As mentioned nine types of events considered in this research.Two hundred of each event generated based on parametric equations of Table 1 by using random parameters. The ranges of equations parameter variation are presented in Table 1. These events sampled by 50Hz. Each event has 16 cycles and 32 samples so the sampling frequency is 1.6 kHz. To assess accuracy of proposed method mean absolute error (MAE) measure introduced as follows:
Overall accuracy and misclassification error
The overall accuracy and misclassification error of each event demonstrated in Table 2.
Classification results under noisy condition
An another important issue related to the proposed method is noisy conditions which are present in an electrical distribution network. So, the performance of the proposed method was tested in different noise Conditions. Gaussian white noise is widely used in the researches on the PQ issues [26]. To evaluate the performance of the proposed method in different noise environments, different levels of noises with the signal to noise ratio (SNR) values from 20 and 40 dB were considered. The value of SNR is defined as follows:

Voltage sag under 20dB noise.
The greatest advantage of the proposed method based on the ST and T2FK-SVM classifier has a smaller feature size for each event. Hence, total simulation time reduced significantly. To verify the efficiency and robustness of the proposed method, overall accuracy compared with [27–30] and error of each method demonstrated in Table 3. As seen from Table 3 the proposed method has significant results and better accuracy than other methods.
Comparison with other methods.
Comparison with other methods.
To evaluate the performance of the proposed method for real-world measurements, 34 waveforms with different types of multiple disturbances taken from the DOE were presented to the classifier. In Figs. 7 and 8 real voltage sag and swell are presented. After applying proposed method, total of 29 waveforms were correctly classified. From these 34 waveforms, 19 are corrupted by sag, 6 are corrupted by swell and 9 are corrupted by harmonics.

Real voltage sag.

Real voltage swell.
The large number of PQ events occurred in the electrical power distribution system. So, there is a requirement for a new method that can be used for automated classification of PQ events. This paper introduced a novel algorithm to classify the PQ events. This technique is based on the S-transform (ST) and type-2 fuzzy kernel support vector machine (T2FK-SVM) classifier. PQ events generated by equations and features are extracted by ST the extracted features classified by SVM. The kernel function of SVM is type-2 fuzzy functions that have a wide range of uncertainty. The parameters of type-2 fuzzy functions are set by trial and error method. The T2FK-SVM is trained by 50% of data and tested by remain data. The main advantage of proposed method is small features size, low execution time and significant accuracy. So, this method requires less memory space at both the training and testing processes.
