# An Improved Walsh Function Algorithm for Use in Sinusoidal and Nonsinusoidal Power Components Measurement

Authors: Saifulnizam Bin Abdul Khalid; Garba Aliyu; Mohd Wazir Mustafa; Hussein Shareef
Journal: Journal of Energy
 (2013)

---

## Abstract

This paper presents an improved Walsh function IWF algorithms as an alternative approach for active and reactive powers measurement in linear and nonlinear, balanced and unbalanced sinusoidal three-phase load system. It takes advantage of Walsh function unified approach, simple algorithm and its intrinsic high level of accuracy as a result of coefficient characteristics and energy behaviour representation. The developed algorithm was modeled on the Matlab Simulink software; different types of load, linear and nonlinear, were also modeled based on practical voltage and current waveforms and tested with the proposed improved Walsh algorithms. The IEEE standard 1459–2000 which is based on fast Fourier transform FFT approach was used as benchmark for the linear load system. The data obtained from laboratory experiment to determine power components in harmonic load systems using Fluke 435 power quality analyzer PQA which complies with IEC/EN61010-1-2001 standard was modeled and used to validate the improved algorithm for nonlinear load measurement. The results showed that the algorithm has the potential to effectively measure three-phase power components under different load conditions.


---

## Body

## 1. Introduction

The awareness of the problems associated with nonlinear devices and their effects on power components measurement algorithms has considerably increased. The use of these power electronics equipment and appliances with microprocessors has been shown to contribute to the distortion of the voltage and current waveforms of the supply system at the distribution end [1]. This distortion affects the smooth operation of most algorithms use for power components measurement and invariably has negative effect on the quality of the power supply system. Accurate measurement and evaluation of energy consumption are of utmost importance for effective planning, billing, monitoring, maintenance, and further development of the electricity supply system. The deregulation of the power sector has led electricity operators to seek for ways to make the best use of the supply systems. Proper evaluation of the energy consumption is one of the important tasks in electric power industries, especially in electric bill and the electrical energy quality estimation and control [2, 3]. Therefore, several companies are into design and manufacture of energy meters using different algorithms to meet the increasing demand of smart metering [4]. In the pricing of electric bill based on the value of the integral of the load active power measured using kilowatt-hour meter, the electric utility board incurs losses in revenue for energy delivered to current harmonic generating customers and also those that cause current asymmetry [5]. Most times, customers that do not generate harmonics but are supplied with distorted or asymmetrical voltages are made to pay not only for the useful energy but also for the energy which may have caused malfunction of their equipment [6]. There is no universally accepted definition for the measurement of reactive power component in the presence of nonlinear and nonsinusoidal waveform [7, 8]. The IEEE working group on appropriate determination of billing system based on sinusoidal and nonsinusoidal situation defined distortion power in terms of total fundamental and harmonic constituents with less unwieldy theory [8].The conventional method of determining power consumption is to evaluate using measured values of the voltage (V), current (I), and power factor (pf). The apparent power S, active power P, and reactive power Q are defined as in [9]. This way of power evaluation requires accurate measurement of the root mean square (RMS) values of the voltage and current before performing the multiplication operation so as to determine the apparent power S. Precise determination of RMS values of voltage and current is complex and poses a serious challenge in electrical measurement [10]. For this method of power calculation to be effective in a three phase network the following has to be satisfied.(i)
The RMS voltages and current on all the three phases must be identical.(ii)
Phase angles between voltage and current must be maintained at 120° and are time invariant.These conditions may not be attainable in the real practical world due to the nature of the loads most of which are nonlinear [11].This paper proposed an improved algorithm as an alternative for power components measurement/estimation using the Walsh function which has the advantage of better accuracy and reliability in energy measurement as it takes into consideration among others, elimination of the effect of harmonic on reactive power measurement. Walsh function transform algorithm has a unique and essentially high level of accuracy as a result of coefficient characteristics and energy behaviour representation [12]. This proposed algorithm might be integrated into any type of energy meter for measurement.The remaining of this paper is organized as follows. Section2 presents Walsh function and its analytical expression, while Section 3 highlights the steps involved in the development of the proposed improved algorithms for measurement using the Walsh function, Section 4 has the modeling and simulation of the improved algorithms, and Section 5 is the simulation with laboratory experiment results. Lastly, Section 6 is the conclusions.


## 2. Walsh Function (WF) Analytical Expression

The Walsh function method for the evaluation of energy parameters was presented as a mathematical tool to analyze energy meter output behaviour for indepth error detection [9]. The method is faster when compared with other techniques like FFT. It does not require the phase shift of π/2 between voltage and current and it has less computation requirement. However, when the load becomes nonlinear and nonsinusoidal the voltage and current waveforms become distorted. This leads to harmonic which is a multiple of the fundamental frequencies of the current and/or voltage. The authors exert that harmonic current affects the reactive power measurement so another attempt using modified approach for power components measurement in both sinusoidal and nonsinusoidal condition was made; the approach reduced computational demand, though the influence of harmonics to the measurement results which was not accounted for, was the main setback of the algorithms [13–15].The attractions to Walsh-function-based technique for use in energy parameter evaluation are the following.(a)
The Walsh transform analyzes signals into rectangular waveform rather than sinusoidal ones and is computed more rapidly when compared with Fourier transform FT.(b)
Walsh-function-based algorithm contains addition and subtractions only and hence results in considerably simplified hardware implementation of power evaluation process.(c)
The IEEE/IEC definition of a phase shift ofπ/2 between the voltage and the current signal mainly used for reactive power evaluation is eliminated from signal processing operation when using Walsh function [8].Generalized Walsh functions and transforms were introduced in 1923 by J. L. Walsh, but their application to engineering and other fields did not happen until recently [16] with some basic and enlightening properties of these functions and transforms considered. This function can be applied among other uses to develop an algorithm that would be applicable to nonlinear loads problem analysis. It is a full orthogonal system with unique properties, which include that it has only two values +1 and −1 over specified normalized period T. This greatly influences the effectiveness of signal processing operation as being related to measurement of power components and characteristics of power distribution system. Analytically the Walsh function is expressed as
(1)Wal(n,β)=(-1)∑k=1m(nm-k+1⊕nm-k)βk,
where n is the order of the function form n=1,2,3,…,nm which is the mth coefficients of the n represented in binary code, that is, n=(n0,n1,n2,n3⋯nm)2,  nm=0,1. With m being the highest order WF serial number in the system, β is the argument of WF that defines the coefficients of βk in binary code β=(βk,βk⋯βk)2, βk=0,1, and k=1,2,3,…,m. From (1) the graphical representation of the first sixteenth order Walsh function is generated as shown in Figure 1.Figure 1
The first 16th order of Walsh function.


## 3. Walsh Function Algorithm

The IEEE standard 1459-2000 for the instantaneous voltages (va,vb,vc) and currents (ia,ib,ic) in a three-phase distribution network with linear balanced or unbalanced load is given as [17]. Instantaneous three-phase sinusoidal voltages are
(2)va=2Vasin(ωt),vb=2Vbsin(ωt-120),vc=2Vcsin(ωt+120).
Instantaneous three-phase currents are
(3)ia=2Iasin(ωt-θa),ib=2Ibsin(ωt-θb-120),ic=2Icsin(ωt-θc+120),
where Va,Vb,Vc and Ia,Ib,Ic are the RMS values of the line to neutral voltages and currents for the phases a, b, and c, respectively. θa, θb and θc are the phase angles between the respective voltage and current. ω=2πf, f=frequencyinHertz, t=1/f.The instantaneous reactive powers for the three phases represented with saypa,pb, and pc are given by
(4)pa=va×ia,pb=vb×ib,pc=vc×ic.
Substituting for (va,vb,vc) and (ia,ib,ic) in the previous expressions and solve trigonometrically the instantaneous powers for phases a,  b, and c of the three-phase linear balanced and unbalanced load network are derived as in (5):
(5)pa=Pa-[Pacos(2ωt)-Qasin(2ωt)],pb=Pb+Pb-3Qbcos(2ωt)+(Pb3+Qb)sin(2ωt),pc=Pc+(Pc+Qc3)cos(2ωt)+(Qc-Pc3)sin(2ωt),
where Pa=VaIacosθa, Pb=VbIbcosθb, Pc=VcIccosθc, Qa=VaIasinθa, Qb=VbIbsinθb, and Qc=VcIcsinθc, respectively.From the instantaneous power of (5) the algorithms for measuring the reactive powers on these three-phase systems say, Qa, Qb, and Qc, using the Walsh function are obtained by multiplying both sides of the equations by the third-order Walsh function, that is, Wal(3,t) and integrate over the period T of power system frequency. The integral of the third order Walsh function with a constant and also a multiplier of cos(2ωt) is equal to zero. Because cos(2ωt) is orthogonal (right angled) with the third Walsh function so the result is obtained as follows:
(6)Qa=-π2T∫0TpaWal(3,t)dt,Qb=π2T∫0TpbWal(3,t)dt-Pb3,Qc=π2T∫0TpcWal(3,t)dt+Pc3.
The algorithms for real or active powers in the three phases a, b, and c are determined by multiplying (5) by the zero-order Walsh function, that is, Wal(0,t) and integrate over the period T. In Walsh algorithm the zero-order function is a constant Wal(0,t)=1, over the period of T as can be seen from Figure 1, so all the integral terms on the right-hand side of (5) that involve product of Wal(0,t) with either cos(2ωt) or sin(2ωt) approach to zero thus giving (7) as
(7)1T∫0TpaWal(0,t)dt=1T∫0TWal(0,t)Padt,1T∫0TpbWal(0,t)dt=1T∫0TWal(0,t)Pbdt,1T∫0TpcWal(0,t)dt=1T∫0TWal(0,t)Pcdt.
Solving the equations further active powers Pa, Pb, and Pc are obtained for three-phase power system as follows:
(8)Pa=1T∫0TWal(0,t)padt,Pb=1T∫0TWal(0,t)pbdt,Pc=1T∫0TWal(0,t)pcdt.
The modeling and simulation of (6) and (8) have been presented earlier in [12] and the authors showed that current harmonic affects the reading of the simulation results. The algorithms have to be further improved as to be able to measure the power components in both linear and nonlinear, sinusoidal nonsinusoidal load conditions. It is worthy to mention, at this juncture, that in AC network source voltages are relatively pure sinusoidal waveforms and that it is mostly at the transmission and distribution end that harmonic and other forms of distortions are introduced. The use of statistical estimators to detect distortions like sags and swells using synthetics analysis that allowed the identification of the variance, skewness, and kurtosis associated with normal sinusoidal waveforms was proposed in [18]. The increase use of nonlinear loads like computers, fluorescent lamps, adjustable speed drive motors, arc furnaces, arc welding machines, electronic control and power converters among others cause harmonics. This lead to overloading of neutral conductors, overheating of transformers, tripping of circuit breakers, power factor correction capacitors over stress and skin effect [19]. Also research has shown that the effect of odd harmonic is more significant in power distribution system while that of even harmonic is negligible [20].


## 4. Proposed Improved Walsh Function Algorithm for Measurement

Harmonic in the power system does not affect the active power measurement but has a great deal of influence on reactive power measurement which invariably affects the power factor and hence the quality of the supply system. To derive the improved Walsh function IWF algorithm for measurement we assume that the load current is contaminated with third-order current harmonic denoted asia3, ib3, and ic3, with θa3, θb3, and θc3 being the phase angles between the fundamental voltages, and the third-order current harmonic waveforms of the phases, Ia3, Ib3, and Ic3 are the RMS values of the third-order current harmonic as given below:
(9)ia3=Ia3sin(3ωt-θa3),ib3=Ib3sin(3ωt-θb3-120),ic3=Ic3sin(3ωt-θc3+120).
The instantaneous powers pa3, pb3, and pc3 for the three phases a, b, and c under this condition are derived as follows:
(10)pa3=Pa+(Pa3-Pa)cos(2ωt)+(Qa3-Qa)sin(2ωt)-Pa3cos(4ωt)-Qa3sin(4ωt),pb3=Pb+(Pb3+Pb-3Qb)cos(2ωt)+(Pb3+Qb3-Qb)sin(2ωt)-(3Qb3-Pb3)cos(4ωt)+(Qb3+Pb3)sin(4ωt),pc3=Pc+(Pc3+Pc-3Qc)cos(2ωt)+(Pc3+Qc3-Qc)sin(2ωt)-(3Qc3-Pc3)cos(4ωt)+(Qc3+Pc3)sin(4ωt),
where Pa3=VaIa3cosθa3, Pb3=VbIb3cosθb3, Pc3=VcIc3cosθc3, Qa3=VaIa3sinθa3, Qb3=VbIb3sinθb3, and Qc3=VcIc3sinθc3, respectively.To obtain the improved algorithm for reactive power under this harmonic condition we apply the Walsh function by multiplying (10) with the third-order WF, that is, Wal(3,t) and integrate over the time T. According to the Walsh functions all the integrals of the right-hand side terms of (10) that involves the multipliers of cos(2ωt), cos(4ωt), and sin(4ωt) and the constant P are all equal to zero [14]. Hence (11) become
(11)1T∫0TpaWal(3,t)dt=1T∫0T(Qa3-Qa)sin(2ωt)×Wal(3,t)dt,1T∫0TpbWal(3,t)dt=1T∫0T(Pb+3Qb3-Qb)sin(2ωt)×Wal(3,t)dt,1T∫0TpcWal(3,t)dt=1T∫0T(Pc+3Qc3-Qc)sin(2ωt)×Wal(3,t)dt.
The product of the 3rd-order WF the expressions, (Qa3-Qa)sin(2ωt), (Pb3+Qb3-Qb)sin(2ωt), and (Pc3+Qc3-Qc)sin(2ωt), results in the full-wave rectification of the terms:
(12)1T∫0TpaWal(3,t)dt=1T∫0T(Qa3-Qa)sin(2ωt)dt,1T∫0TpbWal(3,t)dt=1T∫0T(Pb3+Qb3-Qb)sin(2ωt)dt,1T∫0TpcWal(3,t)dt=1T∫0T(Pc3+Qc3-Qc)sin(2ωt)dt.
Solving for Qa, Qb, and Qc(13)Qa=-π2T∫0TpaWal(3,t)dt+Qa3,Qb=-π2T∫0TpbWal(3,t)dt+Pb3+Qb3,Qc=-π2T∫0TpcWal(3,t)dt+Pc3+Qc3.Qa3, Qb3, and Qc3 in (13) are the reactive components of the distortion power in the phases; these indicate the influence of the third-order current harmonics ia3, ib3, and ic3 on the reactive power measurement algorithms. The third-order Walsh function eliminates the effect of the third-order harmonics of the nonlinear load. Harmonics affect the reactive power measurement of a network. The final terms of (10) are the distortion power terms, that is, Qa3sin(4ωt), (Qb3+Pb3)sin(4ωt), and (Qc3+Pc3)sin(4ωt). They are oscillating with the frequency of 4ω which is similar to the oscillating frequency of the 7th-order Walsh function, Wal(7,t), as can be seen in Figure 2.(a) 7th order WF, (b)sin(4ωt) waveforms.
(a)
(b)To estimate these distortions power terms both sides of (10) are multiplied by the 7th-order WF and then integrated over the period T and simplified to obtain the following:
(14)1T∫0TpaWal(7,t)dt=-1T∫0TQa3sin(4ωt)×Wal(7,t)dt,1T∫0TpbWal(7,t)dt=1T∫0T(Qb3+Pb3)sin(4ωt)×Wal(7,t)dt,1T∫0TpcWal(7,t)dt=1T∫0T(Qc3+Pc3)sin(4ωt)×Wal(7,t)dt.
The 7th-order WF is the odd function with the frequency similar to the frequency of the distortion terms. The product of the 7th-order WF with the distortion terms results in their rectification. So taking cognizance of these rectifying effects, (14) is rewritten as follows:
(15)1T∫0TpaWal(7,t)dt=-1T∫0TQa3sin(4ωt)dt,1T∫0TpbWal(7,t)dt=1T∫0T(Qb3+Pb3)sin(4ωt)dt,1T∫0TpcWal(7,t)dt=1T∫0T(Qc3+Pc3)sin(4ωt)dt.
Solving for Qa3, Qb3, and Qc3 produces the following:
(16)Qa3=-π2T∫0TpaWal(7,t)dt,Qb3=π2T∫0TpbWal(7,t)dt-Pb3,Qc3=π2T∫0TpcWal(7,t)dt-Pc3.
Equations (16) are the improved Walsh function algorithms for measuring the distortion power in a three-phase system. Substituting in (13) algorithms for reactive powers are obtained as follows:
(17)Qa=-π2T[∫0TpaWal(3,t)dt+∫0TpaWal(7,t)dt],Qb=[∫0TpbWal(3,t)dt+∫0TpbWal(7,t)dt],Qc=[∫0TpcWal(3,t)dt+∫0TpcWal(7,t)dt].
These algorithms eliminate the effect of the 3rd- to 7th-order harmonics on the reactive power measurement and also essentially reduced the effect of the higher order current harmonics. In other to reduce the computation involved, the 3rd- and 7th-order Walsh are added together which gives the new improved algorithms. The analytical addition of the 3rd and 7th is represented in Figure 3.(a)Wal(3,t), (b) Wal(7,t), and (c) Wal(3;7,t).
(a)
(b)
(c)From Figure3 it can be observed that the IWF as a result of the addition of standard 3rd- and the 7th-order WF is defined as follows.Wal
(
3
;
7
,
t
) is+1, ift is in the interval [0;T/8], [T/2;5T/8]0, ift is in the interval [T/8;3T/8], [5T/8;7T/8]−1, ift is in the interval [3T/8;T/2], [7T/8;T].The new improved WF algorithms are derived by multiplying (10), with Wal(3;7,t) and then integrate over the time T. Now considering the integrals after the equal sign, the pa, pb, and pc are constants so are equal to zero because IWF is a periodic function; 2nd and 4th integrals are also equal to zero as they include cosine functions that are orthogonal with the IWF. The 3rd and 5th integrals comprise the rectification of the sin functions waveform, so these integrals are not equal to zero. Thus equations are written as follows:
(18)1T∫0TpaWal(3;7,t)dt=1T∫0TWal(3;7,t)(Qa3-Qa)sin(2ωt)dt-1T∫0TWal(3;7,t)Qa3sin(2ωt)dt,1T∫0TpbWal(3;7,t)dt=1T∫0TWal(3;7,t)(Pb3+Qb3-Qb)sin(2ωt)dt+1T∫0TWal(3;7,t)×(Qb3+Pb3)sin(4ωt)dt,1T∫0TpcWal(3;7,t)dt=1T∫0TWal(3;7,t)×(Pc3+Qc3-Qc)sin(2ωt)dt+1T∫0TWal(3;7,t)×(Qc3+Pc3)sin(4ωt)dt.
Solving the right-hand side integrals of (18) yields the new improved algorithms in the following:
(19)Qa,3;7=-πT∫0TpaWal(3;7,t)dt,Qb,3;7=-πT∫0TpbWal(3;7,t)dt+(Pb3+Qb3),Qc,3;7=-πT∫0TpcWal(3;7,t)dt+(Pc3+Qc3).
Equations (7), (16), and (19) are the proposed improved Walsh function IWF algorithms that would be used to measure the active, distortion, and reactive powers, respectively, of a network and also eliminate the effect of higher order harmonic in the three-phase reactive power measurement system. Suffice it to say that in actual cases only lower order harmonics are present in power system signal.


## 5. Modeling the Proposed IWF Algorithm

Equations (7) and (19) are used to create the model for the active and reactive powers measurement based on the proposed improved algorithms using the Matlab Simulink software tool. Then some of the commonly used nonlinear domestic and industrial loads were modeled along and used in the simulation of the improved Walsh function algorithms for active and reactive powers measurement. Examples are compact fluorescent lamps (CFL) and computers. Loads can be resistive R, inductive L, capacitive C, or a combination, for example, RL, RC, LC, RLC depending on the type of load being modeled. RL was used to model the loads used in this simulation (see Figures 4 and 5).Figure 4
The model of the proposed improved Walsh function algorithm.Figure 5
Model of nonlinear load.TheRL nonlinear loads were modeled using universal bridge rectifier; the values for the Rs and Cs were determined using the expressions
(20)Rs>2TsCs,Cs<Pn1000(2πf)Vn2,
where Pn = nominal power of three phase converter VA, Vn = nominal line AC voltage V rms, f = fundamental frequency Hertz, and Ts = sample time.TheR and L of an RL load on dc side of the universal bridge rectifier are calculated using data recorded on the standard 435 fluke power quality analyzer PQA during the laboratory experiment involving measurement of power components of harmonic generating loads.


### 5.1. Simulation and Discussion of Results

For the simulation of the linear loads a synthetic line to neutral voltages (V phases) is chosen as follows:
(21)Va=220∠0°,Vb=220∠-120°,Vc=220∠120°.
Load impedances for two different cases of unbalanced three-phase systems are chosen as shown in cases A and B below for the simulation to verify the algorithms. Knowing the supply voltage and the impedance of the loads (Z) the resistance (R) and inductive reactance (XL) of each phase load were calculated and the obtained results for R and XL were used to configure standard RLC load taken from the Matlab simulink power blocks.Case A:Za=36+j20Ω,Zb=55+j15Ω,Zc=15+j11Ω.Case B:Za=25+j20Ω,Zb=17+j60Ω,Zc=18+j38Ω.Figure 6 shows the voltage and current waveform for the three-phase sinusoidal linear unbalanced load system obtained from simulation with the new improved algorithms.(a) 3-phase voltage waveform and (b) 3-phase unbalance current waveform.
(a)
(b)


### 5.2. Linear Load System

For the linear sinusoidal load system, IEEE standard 1459-2000 which is based on FFT approach was used as benchmark for measurement of active and reactive powers of case A for the three phases and compared with the results using the proposed improved algorithms (7) and (19). The results are displayed in Table 1. Similarly for case B results are shown in Table 2. From the results using the IEEE standard 1459-2000 which is based on FFT as reference it can be observed that the proposed Walsh function algorithm has the potential to measure the active and reactive powers with a high degree of accuracy. In terms of computational complexity the proposed Walsh function algorithm has less computational requirement than the FFT.Table 1
The result of case A.
FFT approach
Active powerP (watts)
Reactive powerQ (var)
Phasea
1181.1
944.8
Phaseb
214.1
746.5
Phasec
492.6
1039.1
Total
1887.8
2730.4
Proposed IWF method
Phasea
1180.7
944.5
Phaseb
214.8
746.7
Phasec
493.5
1040.2
Total
1889
2731.4Table 2
The result of case B.
FFT approach
Active powerP (watts)
Reactive powerQ (var)
Phasea
1027.5
569.5
Phaseb
819.3
223.3
Phasec
2098.6
1538.9
Total
3945.4
2331.7
Proposed IWF method
Phasea
1028
569.6
Phaseb
818.6
223.8
Phasec
2099
1539
Total
3945.6
2332.4


### 5.3. Nonlinear Harmonic Load

Experimental results using harmonic load bank and other harmonic generating loads in the laboratory as recorded by fluke PQA meters were used in modeling of the nonlinear loads. The computer and CFL lamps were modeled asRL loads and simulated using both FFT approach and the proposed improved Walsh algorithms for the measurement of active and reactive powers. The laboratory experimental results were used as benchmark. The nonlinear load currents waveform for the experiment using the Fluke PQA meter and that of the model simulation is recorded as shown in Figures 7 and 8. The results of the simulations and experimental measurement of nonlinear loads are presented in Table 3. It can be observed that the standard FFT approach recorded significant error which could be due to the spectral leakage and picket fence effect phenomenon of the FFT approach when used for measurement in nonlinear harmonic load system. On the other hand, the proposed improved Walsh function algorithm for measurement gave a near accurate result as the influence of the harmonics of the nonlinear load has been address in the improved Walsh algorithms. The Fluke power quality analyzer and FFT graphical user interface analyze signals of one phase at a time so the figure shown is for one of the phases and the same is applicable to the other two phases but has been withheld for convenience, clarity, and space. Walsh function transform is a special form of the Fourier transforms [21, 22]. So it is sometimes called Walsh Fourier.Table 3
Simulation and experimental results.
S/N
Load type
Fluke 435 PQA meter readings
FFT approach
Proposed IWF algorithm
P (watts)
Q (var)
P (watts)
Q (var)
P (watts)
Q (var)
1
CFL lamps  (15.8 watts × 45)
711
364.5
710.8
363.4
711
364
2
Desktop computers (60.3 watts × 70)
422.1
20.3
421.9
19.8
422
20
3
Laptop computers (54.6 watts × 8)
436.8
18.5
435.7
18.3
437
18
4
2 ft. fluorescent lamps (27.8 watt × 15)
417
934.5
416.3
933.8
417
934(a) Nonlinear current waveform and (b) harmonic spectrum (experiment).
(a)
(b)(a) Nonlinear current waveform and (b) harmonic spectrum (Simulation).
(a)
(b)


## 5.1. Simulation and Discussion of Results

For the simulation of the linear loads a synthetic line to neutral voltages (V phases) is chosen as follows:
(21)Va=220∠0°,Vb=220∠-120°,Vc=220∠120°.
Load impedances for two different cases of unbalanced three-phase systems are chosen as shown in cases A and B below for the simulation to verify the algorithms. Knowing the supply voltage and the impedance of the loads (Z) the resistance (R) and inductive reactance (XL) of each phase load were calculated and the obtained results for R and XL were used to configure standard RLC load taken from the Matlab simulink power blocks.Case A:Za=36+j20Ω,Zb=55+j15Ω,Zc=15+j11Ω.Case B:Za=25+j20Ω,Zb=17+j60Ω,Zc=18+j38Ω.Figure 6 shows the voltage and current waveform for the three-phase sinusoidal linear unbalanced load system obtained from simulation with the new improved algorithms.(a) 3-phase voltage waveform and (b) 3-phase unbalance current waveform.
(a)
(b)


## 5.2. Linear Load System

For the linear sinusoidal load system, IEEE standard 1459-2000 which is based on FFT approach was used as benchmark for measurement of active and reactive powers of case A for the three phases and compared with the results using the proposed improved algorithms (7) and (19). The results are displayed in Table 1. Similarly for case B results are shown in Table 2. From the results using the IEEE standard 1459-2000 which is based on FFT as reference it can be observed that the proposed Walsh function algorithm has the potential to measure the active and reactive powers with a high degree of accuracy. In terms of computational complexity the proposed Walsh function algorithm has less computational requirement than the FFT.Table 1
The result of case A.
FFT approach
Active powerP (watts)
Reactive powerQ (var)
Phasea
1181.1
944.8
Phaseb
214.1
746.5
Phasec
492.6
1039.1
Total
1887.8
2730.4
Proposed IWF method
Phasea
1180.7
944.5
Phaseb
214.8
746.7
Phasec
493.5
1040.2
Total
1889
2731.4Table 2
The result of case B.
FFT approach
Active powerP (watts)
Reactive powerQ (var)
Phasea
1027.5
569.5
Phaseb
819.3
223.3
Phasec
2098.6
1538.9
Total
3945.4
2331.7
Proposed IWF method
Phasea
1028
569.6
Phaseb
818.6
223.8
Phasec
2099
1539
Total
3945.6
2332.4


## 5.3. Nonlinear Harmonic Load

Experimental results using harmonic load bank and other harmonic generating loads in the laboratory as recorded by fluke PQA meters were used in modeling of the nonlinear loads. The computer and CFL lamps were modeled asRL loads and simulated using both FFT approach and the proposed improved Walsh algorithms for the measurement of active and reactive powers. The laboratory experimental results were used as benchmark. The nonlinear load currents waveform for the experiment using the Fluke PQA meter and that of the model simulation is recorded as shown in Figures 7 and 8. The results of the simulations and experimental measurement of nonlinear loads are presented in Table 3. It can be observed that the standard FFT approach recorded significant error which could be due to the spectral leakage and picket fence effect phenomenon of the FFT approach when used for measurement in nonlinear harmonic load system. On the other hand, the proposed improved Walsh function algorithm for measurement gave a near accurate result as the influence of the harmonics of the nonlinear load has been address in the improved Walsh algorithms. The Fluke power quality analyzer and FFT graphical user interface analyze signals of one phase at a time so the figure shown is for one of the phases and the same is applicable to the other two phases but has been withheld for convenience, clarity, and space. Walsh function transform is a special form of the Fourier transforms [21, 22]. So it is sometimes called Walsh Fourier.Table 3
Simulation and experimental results.
S/N
Load type
Fluke 435 PQA meter readings
FFT approach
Proposed IWF algorithm
P (watts)
Q (var)
P (watts)
Q (var)
P (watts)
Q (var)
1
CFL lamps  (15.8 watts × 45)
711
364.5
710.8
363.4
711
364
2
Desktop computers (60.3 watts × 70)
422.1
20.3
421.9
19.8
422
20
3
Laptop computers (54.6 watts × 8)
436.8
18.5
435.7
18.3
437
18
4
2 ft. fluorescent lamps (27.8 watt × 15)
417
934.5
416.3
933.8
417
934(a) Nonlinear current waveform and (b) harmonic spectrum (experiment).
(a)
(b)(a) Nonlinear current waveform and (b) harmonic spectrum (Simulation).
(a)
(b)


## 6. Conclusion

A new improved Walsh function algorithm for active and reactive powers measurement was presented. The IEEE standard 1459-2000 fast Fourier transform approach was used as benchmark to validate the algorithm when the load system is linear and sinusoidal. The Fluke 435 power quality analyzer PQA meter was used as reference for nonlinear harmonic load system simulation and measurement. The results show that the proposed improved Walsh function has the potentials of accurately measuring power components under different load conditions. If the algorithm is integrated into measuring instrument it would be able to properly record the active and reactive powers components under sinusoidal and nonsinusoidal load system for proper monitoring, management, and maintenance planning. The approach is simpler with less computation when compared with Fourier transform. The research is continuing in the development of a model for measuring the distortion component of the power system.



---

*Source: 807639-2013-04-10.xml*