Mathematical Problems in Engineering

Volume 2014 (2014), Article ID 905108, 6 pages

http://dx.doi.org/10.1155/2014/905108

Research Article

Active Power Measurement Based on Multiwavelet Transforms

Xiao-Bing Zhang, Yun-Hui Li, and Xiao-Meng Cui

School of Measurement and Communication, Harbin University of Science and Technology, Harbin 150080, China

Received 29 June 2013; Revised 31 December 2013; Accepted 2 January 2014; Published 11 February 2014

Academic Editor: Claudio R. Fuerte-Esquivel

Copyright © 2014 Xiao-Bing Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This paper discusses a new method for calculating active power in the multiwavelet domain. When the voltage and current waveforms are analyzed using multiwavelet, the active power can be calculated by simply adding the products of the multiwavelet coefficients without having to reconstruct the signals back to the time domain first and then using the traditional integration. From the simulation result, we can see that the results using multiwavelet are better than the ones using wavelet and Fourier Transforms no matter which prefilter is used.

1. Introduction

Active power is important for many purposes such as designing power system equipment, setting tariffs, developing measurement meters, and designing compensation devices for improving the quality of the electric energy. Under sinusoidal conditions, the definitions of active power work well. However, due to the recent widespread use of nonlinear loads, the voltage and current waveforms become nonsinusoidal and therefore their traditional definitions become unsuitable. As a result, many attempts have been made to define active power under this new situation [1–5].

Wavelet is an effective tool for nonstational signal processing and has been used in the measurement of active power [6–8]. However, scalar wavelets cannot contain orthogonality, symmetry, compact support, and higher order of vanishing moments simultaneously. Multiwavelet transform is a new concept in the framework of wavelet transform but has some important differences. It simultaneously possesses orthogonality, compact support, higher order of vanishing moments, and symmetry. It has become a tool for power quality study recently [9–11].

In this paper, a new approach to measure active power based on multiwavelet transforms is studied. In Section 2, we introduce the multiwavelet transform. In Section 3, a new approach to measure active power based on multiwavelet transforms is discussed. In Section 4, an example is given to illustrate validity of our method. At last a conclusion is given in Section 5.

2. Multiwavelet Transform

Unlike a scalar wavelet, a multiwavelet has several mother wavelet functions (scaling functions) which are used to expand a given function [12]. Let denote a compactly supported orthonormal scaling vector , where is the number of scalar scaling functions. Then satisfies a two-scaling dilation equation of the form

for some finite sequence of matrices. Furthermore, the integer shifts of the components of form an orthonormal system; that is,

Let denote the closed span of and define . Then is a multiresolution analysis of [13]. Note that we choose the increasing convention .

Let denote the orthogonal complement of in . Then there exists an orthogonal multiwavelet such that forms an orthogonal basis of . Since , there exists a sequence of matrices such that

Let ; then can be written as a linear combination of the basis in . Consider

for some sequence of vector coefficients. The superscript stands for the complex conjugate transposition. Here denotes the space of finite energy vector sequences with the norm

Because [14], where the symbol denotes an orthogonal direct sum can be written as a linear combination of the basis functions of and

The coefficients and are related to via the following decomposition and reconstruction algorithm:

Define and ; then and are stable bases of and , respectively. And (6) can be written as

Because , where is the levels for decomposition, can be written as

In general, for a signal , it can be expressed as a linear combination of the basis in . Consider

Similar to (9),

where (levels for decomposition) and , are row vectors.

Multiwavelets are a new addition to the body of wavelet theory. Realizable as matrix-valued filter banks leading to wavelet bases, multiwavelets offer simultaneous orthogonality, symmetry, and short support, which is not possible with scalar two-channel wavelet systems. The DGHM multiwavelets which are developed by Donovan et al. [15] are a very popular multiwavelet. They are shown in Figure 1.

Figure 1: GHM multiple scaling functions (top) and multiwavelets (bottom).

Remark 1. It is important to note that the orthogonality conditions for the DGHM scaling function and the associated multiwavelet remain valid even when restricted to a compact interval. This is not true in general. It is actually another favorable property of this multiwavelet system.

The DGHM scaling functions satisfy (1) with four coefficients:

And the multiwavelets satisfy (1) with four coefficients as

GHM multiple scaling functions and multiwavelets are very popular with many advantages over scalars. They are either symmetric or antisymmetric, orthogonal, supported on and , respectively, with approximation order 2. However, it is impossible for scalar wavelets [16].

Let and . By dilations (1) and (3), we have the following recursive relationship between the coefficients and :

Moreover,

Due to their matrix-valued filter banks, multiwavelets differ from scalar wavelets in requiring two or more input streams. That means the multiwavelet decomposition and reconstruction are a multi-input and multioutput (MIMO) system. It is necessary to preprocess the input signal before decomposition and postprocess the output signal after reconstruction [17]. Postfilter is just the inverse matrix of prefilter.

Thus, the above decomposition and reconstruction can be represented in Figure 2. is the prefilter and is the postfilter.

(a)

(b)

Figure 2: Multiwavelets (a) decomposition and (b) reconstruction ().

3. Active Power Representation Using Multiwavelet Transform

It is known that the active power is calculated by averaging the voltage and current product with a running window over one cycle of the fundamental frequency

where and are instantaneous voltage and current, and is the period of the specified fundamental frequency.

Let , , and . Then (16) becomes

Suppose that the instantaneous voltage and current have been discretized and expressed using scaling functions and multiwavelets over a running unit interval . Consider

where and are the coarsest resolution level of the decomposition; and are the approximation coefficients that represent the smoothed part of the signal; and are the detail coefficients that represent the oscillatory part of the same signal.

Substituting (18) into (17), we have

Define , , , and ; then (19) can be rewritten as

We further assume the following.

(1)The signal data is zero outside the decomposition window so that the integration range can be extended to the whole real space .(2)The voltage and current signals are already in the same approximation space.(3)The voltage and current waveforms have been decomposed to the same coarsest level; that is, in (20).(4)The multiwavelet analysis of voltage and current signals employs the same type of scaling functions and orthogonal multiwavelets.

Remark 2. Assumption can be satisfied very easy. For example, the periodic signals such as can satisfy assumption .

Assumptions are decided by us. Let the voltage signal and the current signal be and , respectively. They are in the space and have been decomposed to the same coarsest level using the same multiwavelet GHM.

Expand (20) and take into account the following orthogonal properties:

We get

Remark 3. The orthogonal properties (21) can effectively reduce or eliminate the aliasing phenomenon. They ensure the orthogonality of the frequency band.

Using (14), we have

where and with .

From (23), we can see that, in the multiwavelet domain, the power delivered is the addition of the power calculated at each decomposition level, and the power at a decomposition level is the addition of the products of the corresponding multiwavelet coefficients.

It is interesting to note that (23) resembles to a certain extent the formula for calculating active power by wavelet [7] and by Fourier series source [18].

4. Numerical Example

A numerical example is considered in this section. All computations in this section are carried out by Matlab 7.8.0.347.

The example considers the case of sinusoidal voltage ( Hz) and a nonsinusoidal current containing the fundamental, third, and fifth harmonic components

In the field of power systems, it is generally considered up to the 50th harmonic signals; that is,  Hz. According to the Nyquist sampling theorem, the sampling frequency is greater than or equal to 2 times the highest frequency; that is,  Hz. So  Hz is selected as the sampling frequency. You can also choose the other frequencies to be greater than  Hz.

After sampling the voltage and current waveforms with a sampling frequency ( Hz), a prefilter is used to produce the initial coefficients. Then the multiwavelet transform is applied in order to get the approximations and the details.

Since we want to extract the fundamental signal in original signal, it needs the fundamental signal to be in the center of the frequency bank. The lowest decomposition approximation frequency range is 0–100 Hz, which means that the width of approximation frequency bank is  Hz. The width of approximation frequency bank is , where is the decomposition level. So we have . If it is further decomposed to level 6, the fundamental signal is then at the boundary of frequency range 0–50 Hz and 50–100 Hz. If the decomposition level is 4, then the approximation frequency range of the signal is 0–200 Hz, resulting in 3rd harmonic signal that also falls in this frequency range, and the fundamental signal cannot be separated. Thus five wavelet levels are chosen, and this represents the best number of levels for decomposition; therefore, each frequency component of the analyzed waveform can be extracted in a single band.

The voltage and current signals are analyzed using multiwavelet GHM, Daubechies (Db4, Db10), and Fourier Transform (FT). When using wavelet and multiwavelet, the boundary should be handled first. There are many handle methods such as zero padding, symmetric extension, and periodic extension. Zero padding does not preserve orthogonality. It is shown in [19] that the finite multiwavelet transform based on symmetric extension will then preserve the symmetry across scales as in the scalar case. The periodic extension approach will work for multiwavelets and will preserve orthogonality unless the data are truly periodic. Because our data are periodic, here we are using the periodic extension. When using multiwavelet, it also needs a prefilter to transfer the scalar signal into a vector. We are using the two popular prefilters, one is the oversampling [20] and other is GHM.init [17].

The calculation results of the active power using GHM with oversampling prefilter, GHM with GHM.init prefilter, Db4, Db10, and FT are listed in Table 1. We can see that an active power of about 500 W is obtained, and the difference is too small. Furthermore, we can see that the active power measured using multiwavelet no matter which prefilter is used is better than using wavelet (Db4 and Db10) and Fourier Transforms.

Table 1: Active power calculation results.

Method name
Prefilter
Theory values
Calculated values
Error

GHM
Oversampling
500
500.3804

GHM
GHM.init
500
500.3755

db4

500
501.8237

db10

500
499.4937

FT

500
501.2168

Remark 4. From the results, we can see that the active power measured using multiwavelet no matter which prefilter is used is better than using wavelet (Db4 and Db10) and Fourier Transforms. However, multiwavelet transform needs prefiltering and postfiltering for one-dimensional signals. Calculation of multiwavelet transform is greater than wavelet transform. Fast algorithm for multiwavelet transform will be our next work.

5. Conclusions

A new method of active power measurement was mathematically examined in the paper. It is shown that the active power can be calculated by simply adding the products of multiwavelet coefficients without having to spend time on synthesizing the processed coefficients back to time-domain signals first and then using the traditional integration. The new method is convenient and beneficial when multiwavelet transform has to be performed in the true real-time environment or to compress electric data of large volume. The new method is better than wavelet and Fourier Transforms.

Conflict of Interests

The authors declare that there is no conflict of interests concerning the publication of this paper.

Acknowledgment

This work was partially supported by National Natural Science Foundation of China (51277043).

References

W.-K. Yoon and M. J. Devaney, “Power measurement using the wavelet transform,” IEEE Transactions on Instrumentation and Measurement, vol. 47, no. 5, pp. 1205–1210, 1998.

A. E. Emanuel, “Powers in non-sinusoidal situations, A review of definitions and physical meaning,” IEEE Transactions on Power Delivery, vol. 5, no. 3, pp. 1377–1383, 1990.

X. Zhang, Y. Yan, W. Cao, and H. Wu, “A novel rational measurement method for electric energy of distorted signals,” in Proceedings of the International Conference on New Trends in Information and Service Science (NISS '09), pp. 1046–1051, Beijing, China, July 2009.

A. Ortiz, M. Mañana, C. J. Renedo, S. Pérez, and F. Delgado, “A new approach to frequency domain power measurement based on distortion responsibility,” Electric Power Systems Research, vol. 81, no. 1, pp. 202–208, 2011.

F. Vatansever and A. Ozdemir, “Power parameters calculations based on wavelet packet transform,” International Journal of Electrical Power and Energy Systems, vol. 31, no. 10, pp. 596–603, 2009.

X.-B. Zhang, Y.-H. Liang, and J.-R. Cui, “A novel active power measuring method for power systems with distorted signals,” in Proceedings of the IEEE International Conference on Information and Automation (ICIA '08), pp. 5–9, Zhangjiajie, China, June 2008.

T. X. Zhu, “A new approach to active power calculation using wavelet coefficients,” IEEE Transactions on Power Systems, vol. 21, no. 1, pp. 435–437, 2006.

W. G. Morsi and M. E. El-Hawary, “Reformulating power components definitions contained in the IEEE standard 1459-2000 using discrete wavelet transform,” IEEE Transactions on Power Delivery, vol. 22, no. 3, pp. 1910–1916, 2007.

S. Kaewarsa, K. Attakitmongcol, and T. Kulworawanichpong, “Recognition of power quality events by using multiwavelet-based neural networks,” International Journal of Electrical Power and Energy Systems, vol. 30, no. 4, pp. 254–260, 2008.

Q. Yinglin, T. Lijun, and W. Chunbao, “Power quality disturbances detection and location based on multi-wavelet and neighboring coefficient de-noising,” in Proceedings of the 4th International Conference on Electric Utility Deregulation and Restructuring and Power Technologies (DRPT '11), pp. 424–427, WeiHai, China, July 2011.

R. Kapoor and M. K. Saini, “Multiwavelet transform based classification of PQ events,” European Transactions on Electrical Power, vol. 22, no. 4, pp. 518–532, 2012.

J. S. Geronimo, D. P. Hardin, and P. R. Massopust, “Fractal functions and wavelet expansions based on several scaling functions,” Journal of Approximation Theory, vol. 78, no. 3, pp. 373–401, 1994.

S. G. Mallat, “Theory for multiresolution signal decomposition: the wavelet representation,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 11, no. 7, pp. 674–693, 1989.

F. Keinert, Wavelets and Multiwavelets, Chapman & Hall/CRC, Boca Raton, Fla, USA, 2004.

G. C. Donovan, J. S. Geronimo, D. P. Hardin, and P. R. Massopust, “Construction of orthogonal wavelets using fractal interpolation functions,” SIAM Journal on Mathematical Analysis, vol. 27, no. 4, pp. 1158–1192, 1996.

I. Daubechies, Ten Lectures on Wavelets, SIAM, Philadelphia, Pa, USA, 1992.

X.-G. Xia, J. S. Geronimo, D. P. Hardin, and B. W. Suter, “Why and how prefiltering for discrete multiwavelet transforms,” in Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP '96), vol. 3, pp. 1578–1581, May 1996.

D. E. Scott, An Introduction to Circuit Analysis: A System Approach, McGraw-Hill, New York, NY, USA, 1987.

T. Xia and Q. Jiang, “Optimal multifilter banks: design, related symmetric extension transform, and application to image compression,” IEEE Transactions on Signal Processing, vol. 47, no. 7, pp. 1878–1889, 1999.

A. Vasily Strela and A. T. Walden, “Signal and image denoising via wavelet thresholding: orthogonal and biorthogonal, scalar and multiple wavelet transforms,” Tech. Rep. TR-98-01, Imperial College, Statistics Section, 1998.