Energy and Power Engineering
Vol.11 No.03(2019), Article ID:91468,18 pages
10.4236/epe.2019.113009

Comparative Study of Linear and Non-Linear Controls of Three-Phase Shunt Active Filters for Improving the Quality of Electrical Energy

Mamane Adamou1*, Nabil Alkassoum1, André Foulani1, Amadou Seidou Maiga2

1Department of Physics (FAST), Abdou MOUMOUNI University, Niamey, Niger

2Departments of Applied Physics, Gaston Berger University, Saint-Louis, Senegal

Copyright © 2019 by author(s) and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: December 1, 2018; Accepted: March 25, 2019; Published: March 28, 2019

ABSTRACT

The aim of this work is to comparatively study two types of control of a three-phase shunt active filter (TSAF) in order to realize a system of depollution of electrical networks so as to improve the quality of electrical energy. We used two TSA control models which differ in the method of disturbed currents detection, one linear and the other non-linear. The results show that the non-linear control method, although with high calculation blocks, gives more promising results than the linear control method. When connecting a non-linear load (televisions, lamps, variable speed drives, etc.) directly to the public distribution network, the non-linear nature of the load causes the source current to be deformed due to the presence of harmonic currents and voltages. These harmonics thus generate a high reactive power, and therefore considerable electrical losses in the network. One of the solutions to reduce losses on the electricity grid is to incorporate an electronic circuit called a filter into the electricity grid. This work carried out a comparative study of two types of control of a three-phase shunt active filter (TSAF) in order to carry out a system for cleaning up electricity networks in order to improve the quality of electrical energy. We used two TSA control models that differ in the method of detecting disturbed currents, one linear and the other non-linear. The results show that the non-linear control method, although with high calculation blocks, gives better results than the linear control method.

Keywords:

Power Quality, Active Three-Phase Shunt Filter, Non-Linear Control, Linear Control

1. Introduction

Electrical energy is an energy that can be converted into almost all other existing forms of energy and is managed through an electrical network. The quality of electrical energy received by consumers is very essential for the proper functioning of appliances connected to an electrical network. The increasing use of power electronic components on electrical networks is not without major consequences for the latter. Indeed, an alternating voltage source supplying a non-linear load will undergo harmonic current and/or voltage pollution [1]. This harmonic pollution can have consequences on other devices sensitive to the latter which will be supplied by the same source [2]. To overcome this problem, the techniques initially used were the oversizing of electrical equipment in the network and the use of passive power filters [3]. The limits set by the previous solutions to harmonic pollution control will be addressed by active power filters, which constitute one of the modern solutions to power system pollution control [3] [4].

The performance of an active filter depends on the quality of its control. Several control methods have been proposed in the literature [3] [4] [5] [6] , but these control methods are linear. As the clean-up devices for electricity networks are non-linear, it was necessary to propose non-linear control methods capable of adapting to them.

For a contribution to the improvement of power quality, in this paper, we propose a comparative study by simulation of the linear and non-linear controls of active power filters case of a three-phase active shunt filter (TSAF).

2. Electrical Power Quality and TSAF

The quality of electrical energy is related to any variation in the electrical power supply, resulting in malfunction or damage to users' equipment [5]. These variations in electrical quantities are manifested by voltage dips, overvoltage, harmonic disturbances, imbalances, etc. Figure 1 shows some voltage variations on

Figure 1. Changes in electrical quantities affecting the quality of electrical energy.

an electrical network that can affect the quality of electrical power.

In this work we are interested in improving the quality of electrical power in a balanced three-phase network, by compensating harmonic currents using an active three-phase shunt filter (TSAF). Indeed, the disturbances caused by harmonic currents on a network can also affect the voltage through the impedance of the network because the latter is not zero [3].

A three-phase alternating voltage source, supplying a non-linear load, undergoes harmonic pollution that can have consequences on the devices supplied by the same source. The grid voltage and current will no longer be sinusoidal. The role of the active shunt filter is to inject harmonic currents of the same amplitude as the polluting harmonics into the network, but in phase opposition to the latter [7]. By this principle, TSAF is used to improve the quality of electrical energy. This principle is illustrated by the diagram in Figure 2.

3. Active Three-Phase Shunt Filter (TSAF)

3.1. General Structure of TSAF

The active shunt filter is a modern device for harmonic and often reactive power compensation. It consists of two main parts: power and control. The principle of the active shunt filter lies in the injection of harmonic currents of the same amplitude as those pollutants but in opposition of phase with them. Compensation being in real time, it is therefore essential to define a good control method [8].

Several models of single-phase or three-phase active filters have been proposed in the literature. Their performance is directly linked either to the nature of the storage element (capacitive or inductive) of the active filter or to the control strategy of the latter [3] [4] [5]. In this document we will study an active shunt filter with capacitive storage; we will then apply a linear then non-linear control.

Figure 1 shows the structure of the study system consisting of a network polluted by a non-linear load and the active shunt compensating filter. The source

Figure 2. Improvement of the current waveform of a polluted source of the active power shunt filter.

supplying the system is a balanced three-phase source that can be modeled by the Equations (1).

{ V 1 N = U 2 sin ( ω t ) V 2 N = U 2 sin ( ω t 2 π / 3 ) V 3 N = U 2 sin ( ω t + 2 π 3 ) (1)

The non-linear load is a three-phase diode bridge that supplies a resistive load. The active shunt filter is a voltage inverter i.e. with capacitive storage element, it consists of three arms, each of the arms consists of two bidirectional switches (MOSFET, GTO, IGBT). This structure does not allow the simultaneous closing of the semiconductors of the same arm under penalty of short-circuiting the capacitor [3]. However, in some cases, they may both be open (dead time). The continuity of the currents is then ensured by the conduction of one of the diodes of the same arm. In order to couple the active shunt filter to the mains, LF coupling inductors are used; these inductors are also used to attenuate the harmonics due to the switching of the TSAF switches.

3.2. Modeling the Three-Phase Shunt Active Filter

To set up the control laws of the active shunt filter, it is necessary to model it according to the system in which this latter is implemented. This modeling will allow, thanks to a regulation or a servo system, to optimize the performances of the system [8]. Thus, it was explained in paragraph 2 that the active shunt filter consists of two main parts i.e. power and control.

The TSAF is controlled in two steps, the first to determine the reference currents and the second using the first to determine the control pulses of the TSAF inverter switches. The opening and closing of the TSAF inverter switches depends on the status of the control signals (S1, S2, S3,), as defined by Equation (4). ( [4] [9] ):

S 1 { 1 : T 1 o p e n e d a n d T 4 c l o s e d 0 : T 1 c l o s e d a n d T 4 o p e n d e d

S 2 { 1 : T 2 o p e n e d a n d T 2 c l o s e d 0 : T 5 c l o s e d a n d T 5 o p e n e d

S 3 { 1 : T 3 o p e n e d a n d T 6 c l o s e d 0 : T 3 c l o s e d a n d T 6 o p e n e d

Either S i { 1 : T i o p e n e d a n d T ( i + 3 ) c l o s e d 0 : T i c l o s e d a n d T ( i + 3 ) o p e n e d i = 1, 2, 3 (2)

If V 1 M , V 2 M and V 3 M are the output voltages of each TSAF arm referred to the common point M of the arms and V d c is the voltage across the TSAF storage capacitor, then we can write:

V i M = S i V d c (3)

The modeling of the elements constituting the power part will facilitate the design of the control:

・ Variable side of the TSAF

Consider the diagram in Figure 3; V 1 N , V 2 N et V 3 N are the voltages at the connection points of the TSAF, I F 1 , I F 2 et I F 3 are the currents generated by each arm of the TSAF, V 1 M , V 2 M et V 3 M are the output voltages of each TSAF arm referred to the common point M and V M N is the voltage between M and the neutral (N) of the source. By applying the law of meshes on phases 1, 2 and 3, we obtain the system of Equation (4).

{ V 1 N = L f d I F 1 d t + R f I F 1 + V 1 M + V M N V 2 N = L f d I F 2 d t + R f I F 1 + V 2 M + V M N V 3 N = L f d I F 3 d t + R f I F 1 + V 3 M + V M N (4)

Either V i N = L f d I F i d t + R f I F i + V i M + V M N (5)

With i =1, 2, 3.

The system studied is a balanced three-phase system without zero sequence current, then:

V 1 N + V 2 N + V 3 N = 0 and I F 1 + I F 2 + I F 3 = 0 ,

By summing the three equations of the system of Equation (4), we will obtain Equation (6) giving the relation between V M N and V i M ' s (i = 1, 2, 4).

V M N = 1 3 i = 1 3 V i M (6)

Now according to the Equation (2) we have V i M = S i V d c , considering Equation (5) we can then write Equation (6):

V M N = 1 3 i = 1 3 S i V d c (7)

Figure 3. General structure of an active shunt filter.

Taking into account of Equation (3) and Equation (7) we obtain Equation (8):

V i N = L f d I F i d t + R f I F i + S i V d c 1 3 i = 1 3 S i V d c

V i N = L f d I F i d t + R f I F i + ( S i 1 3 i = 1 3 S i ) V d c

d I F i d t = V i N L f R f L f I F i 1 L f ( S i 1 3 i = 1 3 S i ) V d c (8)

Let us define by V f i the output voltages of the filter referred to the neutral, we have:

V f i = ( S i 1 3 i = 1 3 S i ) V d c = K n i V d c (9)

With K n i = ( S i 1 3 i = 1 3 S i ) n

It can be noticed that i = 1 3 K n i = 0

The mathematical system governing the currents of the active shunt filter is given by the relation (10)

d d t ( I F 1 I F 2 I F 3 ) = R f L f [ 1 0 0 0 1 0 0 0 1 ] ( I F 1 I F 2 I F 3 ) V d c L f ( K n 1 K n 2 K n 3 ) + 1 L f ( V N 1 V N 2 V N 3 ) (10)

With two switches for each of the three inverter arms of the TSAF there are then eight possible sequences for the voltage V f i , the sequential function K n i depends not only on the sequence number n = 0 , 1 , 2 , , 7 but also on the phase number i = 1, 2, 3. Thus, eight possible cases of output voltage of the active filter V f i referred to neutral can be expressed in Table 1 [3] [9].

Let us consider d Q c as the amount of charge stored in the capacitor during a time interval dt, i c the current flowing in the capacitor:

d Q c = i c d t = C d V d c from which

d V d c d t = i c C (11)

Table 1. Possible output voltages of the active filter.

By applying the law of nodes, we will have:

i c = S 1 I F 1 + S 2 I F 2 + S 3 I F 3 , we can then write i c = i = 1 3 S i I F i , it is also verifiable that i = 1 3 S i I F i = i = 1 3 K n i I F i .

The Equation (11) then becomes:

d V d c d t = 1 C i = 1 3 S i I F i = 1 C i = 1 3 K n i I F i (12)

Taking into account the absence of the homo-polar sequence in the current system ( i = 1 3 I F i = 0 ) and the fact that i = 1 3 K n i = 0 , we can then write:

I F 3 = I F 1 I F 2 and K n 3 = K n 1 K n 2 ; we can write the differential equation of the continuous side:

d V d c d t = 1 C ( K n 1 I F 1 + K n 2 I F 2 + K n 3 I F 3 )

d V d c d t = 1 C ( 2 K n 1 I F 1 + 2 K n 2 I F 2 + K n 1 I F 2 + K n 2 I F 1 )

d V d c d t = 1 C ( 2 K n 1 + K n 2 ) I F 1 + 1 C ( 2 K n 2 + K n 1 ) I F 2 (13)

{ d I F 1 d t = V 1 N L f R f L f I F 1 K n 1 L f V d c d I F 2 d t = V 2 N L f R f L f I F 2 K n 2 L f V d c d V d c d t = 1 C ( 2 K n 1 + K n 2 ) I F 1 + 1 C ( 2 K n 2 + K n 1 ) I F 2 (14)

Either in matrix form:

d dt ( I F1 I F2 V dc )=[ R f L f 0 K n1 L f 0 R f L f K n2 L f 1 C ( 2 K n1 + K n2 ) 1 C ( 2 K n2 + K n1 ) 0 ]( I F1 I F2 V dc )+ 1 L f [ V 1N V 2N 0 ] (15)

Notice that this model is a minimal representation in the state space. They are variable in time and non-linear [8] [10].

In the following paragraphs, the two methods, linear and non-linear, will be proposed for the determination of the reference currents of the first step of the TSAF control then the pulse width modulation control will be used for the control of the TSAF inverter switches.

3.3. Linear Three-Phase Shunt Active Filter Control

For the linear control of the three-phase active shunt filter it will be applied d-q transformation method, without implication of the TSAF current for the search of the disturbing currents. This makes the method linear and approximates [8].

The method of detection of disturbance currents by synchronous reference frame or d-q transformation method is a method applicable as well to balanced three-phase systems as to neutral and/or unbalanced three-phase systems. It is based on instantaneous values, hence a good dynamic response time. It consists in transforming the 1-2-3 coordinates of the non-linear load current into d-q coordinates using the Park transform P 32 (18) by setting the frequency of the latter in synchronism with that of the network [3] [4].

The Park transform synchronized with the grid frequency has the effect of transforming the fundamental current component into a DC component while the harmonic current components undergo a shift in the frequency spectrum. The elimination of the DC component is achieved by adding a high-pass filter (HPF). The reference current is obtained by Park reverse transformation P 23 (19).

The algorithm for calculating reference currents using the synchronous reference method is given in Figure 4.

Let [ I c ( 123 ) ] , [ I c ( d q ) ] , [ P 32 ] and [ P 23 ] be the non-linear load current vector in reference frame 1-2-3, in reference frame d-q, the Park transform and the Park inverse transform, respectively.

[ I c ( 123 ) ] = [ I c 1 I c 2 I c 3 ] (16)

[ I c ( d q ) ] = [ I c d I c q ] (17)

[ P 32 ] = 2 3 [ sin ( ω t ) sin ( ω t 2 π / 3 ) sin ( ω t + 2 π / 3 ) cos ( ω t ) cos ( ω t 2 π / 3 ) cos ( ω t + 2 π / 3 ) 1 / 2 1 / 2 1 / 2 ] (18)

[ P 23 ] = 2 3 [ sin ( ω t ) cos ( ω t ) 1 sin ( ω t 2 π / 3 ) cos ( ω t 2 π / 3 ) 1 sin ( ω t + 2 π / 3 ) cos ( ω t + 2 π / 3 ) 1 ] (19)

It should be noticed that [ P 23 ] = [ P 23 ] T = [ P 23 ] 1

The passage of the current vector from reference frame 1-2-3 to reference frame d-q is given by the Equation (20) and the calculation of the reference currents by the Equation (21).

[ I c ( d q ) ] = [ P 32 ] [ I c ( 123 ) ] (20)

Figure 4. Linear control detection algorithm.

[ I c ( 123 ) * ] = [ P 23 ] [ I c ( d q ) ] (21)

[ I c ( 123 ) * ] being the reference current vector.

Linear control models are generally approximate models. In the following paragraph we will study the non-linear control of the TSAF with a view to a comparison with the studies carried out upstream.

3.4. Non-Linear Three-Phase Shunt Active Filter Control

Linear control models are generally approximate models. For a model closer to the real system, it is more interesting to apply a non-linear command. There are several non-linear control approaches [8] [10]. The one used in this paragraph is based on the decoupling and linearization of the equations governing the system modeled above. Thus, it will be carried out a transformation of the base 3 to the base 2 by the Park transformation of the currents of the active shunt filter.

Let us find again the system of Equations (10) resulting from the modeling of the active shunt filter in the d-q reference frame. Let us write the relations allowing the passage from base 2 (d, q) to base 3 (1, 2, 3) of the current vector of the active shunt filter, the sequential function vector K n i , the voltage vector at the connection point V i N and the relation governing the dynamics of the currents of the active shunt filter:

[ I F ( 123 ) ] = [ P 23 ] [ I F ( d q ) ] ,

[ K n ( 123 ) ] = [ P 23 ] [ K n ( d q ) ] ,

[ V N ( 123 ) ] = [ P 23 ] [ V N ( d q ) ]

d d t [ I F ( 123 ) ] = R f L f [ 1 0 0 0 1 0 0 0 1 ] [ I F ( 123 ) ] V d c L f [ K n ( 123 ) ] + 1 L f [ V N ( 123 ) ] (22)

You can then write:

d d t [ P 23 ] [ I F ( d q ) ] = R f L f [ 1 0 0 0 1 0 0 0 1 ] [ I F ( d q ) ] V d c L f [ P 23 ] [ K n ( d q ) ] + 1 L f [ P 23 ] [ V N ( d q ) ] (23)

But

d d t [ P 23 ] [ I F ( d q ) ] = [ I F ( d q ) ] d d t [ P 23 ] + [ P 23 ] d d t [ I F ( d q ) ] (24)

d d t [ I F ( d q ) ] = ( [ P 32 ] d d t [ P 23 ] R f L f [ 1 0 0 0 1 0 0 0 1 ] ) [ I F ( d q ) ] V d c L f [ K n ( d q ) ] + 1 L f [ V N ( d q ) ] (25)

and

[ P 32 ] d d t [ P 23 ] = ω [ P 32 ] [ P 23 ] [ 0 1 0 1 1 0 0 0 1 ]

Hence the Equations (26) and (27) which govern the currents of the active shunt filter in the two-phase plane d-q:

d d t [ I F ( d q ) ] = [ R f L f ω 0 ω R f L f 0 0 0 R f L f ] [ I F ( d q ) ] V d c L f [ K n ( d q ) ] + 1 L f [ V N ( d q ) ] (26)

{ L f d I F d d t + R f I F d = L f ω I F q V d c K n d + V d N L f d I F q d t + R f I F q = ω I F d V d c K n q + V q N (27)

The Equation (28) of the continuous bus can also be rewritten in the two-phase plane d-q:

d V d c d t = 1 C [ K n 1 K n 2 K n 3 ] [ I F 1 I F 2 I F 3 ] (28)

d V d c d t = 1 C [ K n d K n q 0 ] [ P 32 ] [ P 23 ] [ I F d I F q 0 ]

d V d c d t = 1 C K n d I F d + 1 C K n d I F d (29)

The minimum representation in the state space can then be rewritten in the d-q plane:

{ d I F d d t = R f L f I F d + ω I F q V d c L f K n d + V d N L f d I F q d t = R f L f I F q ω I F d V d c L f K n q + V q N L f d V d c d t = 1 C K n d I F d + 1 C K n d I F d (30)

d d t ( I F d I F q V d c ) = [ R f L f ω 0 ω R f L f 0 0 0 R f L f ] ( I F d I F q V d c ) + [ V d c L f 0 0 V d c L f I F d C I F q C ] ( K N d K N q 0 ) + 1 L f [ V d N V q N 0 ] (31)

If we consider the system of Equations (30) we will notice in the first two equations the presence of the current I F d in the differential equation containing the first derivative of the current I F q and vice versa. The third equation is characterized by the multiplication of the components of the sequential function vector by the DC bus voltage V d c . The system is then coupled and non-linear.

To control the filter, the system of Equations (30) would have to be decoupled and linearized [8].

By introducing the new variables E d and E q ,

{ E d = L f ω I F q V d c K n d + V d N E q = ω I F d V d c K n q + V q N

The control law will be defined by K n d and K n q which after transformation from base 2 to base 3, are directly linked to the generation of the control pulses of the active shunt filter by the relation

K n i = ( S i 1 3 i = 1 3 S i ) n .

{ K n d = V d N + L f ω I F q + E d V d c K n q = V q N L f ω I F d + E q V d c (32)

With

E d = L f d I F d d t + R f E q = L f d I F q d t + R f I F q

Thus, the algorithm of non-linear control by the synchronous repository method is given in Figure 5.

4. Simulation Results and Comparison

4.1. Simulation Results

The general structure of the studied system shows in Figure 6 consists of a polluted

Figure 5. Algorithm of detection by non-linear control.

three-phase power source supplying a non-linear load. The active three-phase shunt filter is connected in parallel to the mains in order to decontaminate it.

The power part is identical for both studies. Figure 7, Figure 8 and Figure 9 shows in the Matlab/Simulink environment the contents of the three-phase power supply blocks, the three-phase shunt active filter and the non-linear load, respectively. The common simulation parameters are given in Table 2.

Since the control parts are different, the contents of the disturbance current detection blocks in linear and non-linear control can be seen in Figure 10 and Figure 11.

Figure 6. General structure of the system.

Figure 7. Structure of the three-phase power source.

Figure 8. Three-phase shunt active filter.

Figure 9. Structure of the three-phase non-linear load.

Table 2. Simulation parameters.

Figure 10. Detection of harmonic currents by the linear method.

Figure 11. Detection of harmonic currents by the non-linear method.

The simulation results are given respectively for linear and non-linear control through Figure 12, Figure 13 and Figure 14, which first give the current waveforms and then the harmonic spectra before and after application of the TSAF.

4.2. Comparison of Linear and Non-Linear Controls

The comparison of the two control methods of the active shunt filters will be based on the simulation results: by comparison of the transient and permanent regimes, the waveforms of the output currents, the harmonic distortion rates, the stability of the system. For the same parameters, the harmonic distortion rate after application of the FAS is 1.46% for non-linear control against 1.44% for linear control. The variation in TDH is more significant in non-linear control 28.24% than in linear control 28.05%.

Figure 12 shows that the waveforms of the output voltages I s ( t ) for both types of control became sinusoidal again after application of the active shunt filter. For both types of controls, the transitional regimes are short and stable as shown by the Figure 15.

Table 3 provides a summary of the comparison of linear and non-linear FAST commands.

Figure 12. Waveform of load, harmonic and source currents. (a) Linear control; (b) non linear.

Figure 13. Order of source current harmonics and TDH before FAS action. (a) Linear control; (b) non linear.

Figure 14. Order of source current harmonics and TDH after FAS action. (a) Linear control; (b) non linear.

Figure 15. Transient regime. (a) Linear control; (b) non linear.

Table 3. Comparison of linear and non-linear controls.

5. Conclusion

Linear control models are generally approximate models. Hence, there is the need to explore non-linear models. Most real physics systems are modelled by non-linear equations or systems of equations. Non-linear control is best adopted to control the active filter which is a non-linear system. The compromises will then be the number calculated to raise to manage and above all a stability to maintain. The comparative study, by simulation, of the linear and non-linear controls of three-phase shunt active filters shows that although the non-linear control seems complex, for an improvement of the quality of electrical energy by active filtering, the latter is the best adopted last.

Acknowledgements

We acknowledge support from the African Center of Excellence MITIC of Gaston Berger University in Senegal.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Adamou, M., Alkassoum, N., Foulani, A. and Maiga, A.S. (2019) Comparative Study of Linear and Non-Linear Controls of Three-Phase Shunt Active Filters for Improving the Quality of Electrical Energy. Energy and Power Engineering, 11, 149-166. https://doi.org/10.4236/epe.2019.113009

References

  1. 1. Gupta, S., Das, R.F. and Ashraful, B.C. (2016) A Comparative Study on Harmonics of Different Electric Bulbs. American Journal of Engineering Research (AJER), 5, 156-166.

  2. 2. Abdoulkarimoun, D. (2018) Etudes des perturbations Harmonique et leurs traitements dans le réseau électrique: Cas du réseau BT de la ville de Niamey. PhD Thesis, Abdou Moumouni University, Niger.

  3. 3. Mohamad, A.E.A. (2002) Contribution à l’Etude des Compensateurs Actifs des Réseaux Electriques Basse Tension. PhD Thesis, Louis Pasteur University, Strasbourg.

  4. 4. Hind, D. (2007) Filtrage actif de puissance, PhD Thesis, Mentouri University, Constantine.

  5. 5. Abdelhamid, H. (2010) Contribution à l’étude des filtres hybrides de puissance utilises pour améliorer la qualité de l’énergie dans le réseau électrique de distribution. PhD Thesis, Québec University, Canada.

  6. 6. Shraddha, D., Jalaja, G.P.K., Manashri, Y. and Tuljapurkar, M.D. (2015) Detection of Harmonics in Distribution System Using FFT with Case Studies. International Journal for Scientific Research and Development, 3, 481-484.

  7. 7. Antoni, H.N. (2017) Commande de filtres actifs parallèles sur un réseau fortement perturbé. PhD Thesis, Toulouse University, France.

  8. 8. Sina, S. (2002) Etude d’une commande nonlinéaire sur un filtre actif. PhD Thesis, Québec University, Canada.

  9. 9. Mohamed, M. (2007) Validation des techniques de commande d’un filtre actif parallèle. PhD Thesis, Québec University, Canada.

  10. 10. Slotine, J.-J. and Li, W.P. (1991) Applied Nonlinear Control. Prentice Hall, Upper Saddle River.