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 [15].

Wavelet is an effective tool for nonstational signal processing and has been used in the measurement of active power [68]. 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 [911].

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.

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.

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

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.

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.


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