Research Article | Open Access
Dragutin Kostić, Vladimir Šinik, "Approximate Algorithm for Determining Pulse Edges of a PWM Inverter Based on Natural Sampling", Mathematical Problems in Engineering, vol. 2009, Article ID 495360, 23 pages, 2009. https://doi.org/10.1155/2009/495360
Approximate Algorithm for Determining Pulse Edges of a PWM Inverter Based on Natural Sampling
The paper presents a new algorithm for determination of pulse edges of a modulated wave of a PWM voltage inverter which offers a possibility that natural sampling is realized with an arbitrary accuracy without applying an iterative procedure. The basic idea is to express the angles which determine pulse edges of the modulated signal as polynomials of amplitude modulation index. Geometric interpretation of sampling of the polynomial algorithm is identical with the geometric interpretation of natural algorithm, but the transcendental equation whose solution defines pulse edges of the modulated signal is replaced by a simple procedure of finding values of a polynomial whose coefficients are determined in advance by an exact procedure. This approach gives the possibility of digital implementation of polynomial sampling method using the low-cost microprocessor platforms.
Pulse width modulation is an important method of the technique of control of power converters. At the current technological level, three areas within the pulse width modulation possess certain autonomy as regards the accomplished solutions . These are switching strategy of pulse width modulation, concerned with generation of edges of switching pulses [2–10], range of regulation of amplitude of the fundamental harmonic of inverter output voltage [11–15], and degree of harmonic distortions of inverter output voltage and current [16–21]. The first PWM techniques are based on the method of the modulating function.
The modulating function contains information regarding the desired waveform, whereas the signal carrier contains information concerning the switching frequency. Pulse edges of a modulated signal are determined by crossings of the modulating function and signal carrier, which makes the basis of all analog implementations of the modulating function methods involving natural sampling. Signal at the output of a pulse width modulator defines the switching function of the inverter branch which controls the on-state of the switch within the branch. Synthesis of the ac voltage at inverter's output, regulation of the amplitude, and frequency are accomplished by the switching actions in inverter branches.
The basic PWM method is related to the mode of determination of pulse edges. Determination of pulse edges based on crossing of the modulating and carrying signals belongs to the group of natural sampling techniques.
The basic shortcoming of the natural sampling technique is tied to the transcendental relation which connects the angles determining pulse edges and modulating function. This aggravated the application of natural sampling technique in digital and microprocessor systems applying PWM methods.
Regular sampling method is the basic solution for digital realization of the modulating function method. The regular sampling method is considerably more flexible since the modulating function is specified at discrete points; thus it is possible to calculate its values in advance and interpret them by means of digital words [5, 6, 22]. The natural sampling method has the advantage over the regular method in so far as the harmonic content linearity, and range of regulation of the amplitude of the fundamental harmonic are concerned.
Middle of the eighties marks the beginning of the application of the space vector modulation as the vector approach to PWM for three-phase inverters. The basic advantages of SVM are related to the following.(i)Expanded linear range of modulation without injection of the third harmonics into the sine modulating function [23–25].(ii)Lower harmonics content in relation to the regular methods based on the sine modulating function . (iii)Lower switching losses are conditioned by only one change of state [24, 27] (iv)Simple digital implementation[28, 29]
On the other side in the papers [30, 31], it is shown that the conventional method of SVM is equivalent to the regular method with modified sine-triangle uniform sampling PWM with triplen harmonic injected. In the papers [32, 33] a detailed analysis of carrier based PWM methods is given and their relation with the space vector approach as well. Using the simple transformations, modulating functions translate from the time domain to the complex domain, where the reference vector rotates with the angular speed equal to angular frequency of time modulating functions.
Finally, in the paper  it was demonstrated that the calssical method SVM is identical with the method of modulation sine function injected with triplen harmonics and regular sampling technique.
Taking into account the relevant literature[34–36], we propose an aproximate algorithm for PWM sampling according to the modulating function method. Unlike some other works, in which the natural sampling method was rejected [2, 3, 15], we have accepted the elements of this technique but have substituted the transcendental equation by an algebraic expression derived, applying an exact procedure.
2. Modulating Function Method
General characteristic of this group of methods is that the synthesis of the pulse at the output of the modulator is carried out on the basis of the reference signal and signal carrier.
The reference signal is a periodic function having the frequency equal to the fundamental frequency of the inverter output voltage and amplitude which is proportional to the amplitude of the fundamental harmonic. The modulating function is analytical representation of the reference signal. In general, this is a function of the following general form :
Parameter stands for the modulation amplitude index, and parameter for the fundamental frequency of inverter ac voltage Function can be continuous or discontinuous having first-order interruptions in terms of . However, a great number of modulating functions is of the following simpler form:
In classical PWM methods the carrier signal is an ac signal of triangular half-periods, of frequency , and amplitude equal to one. The ratio of the carrier signal frequency and the reference signal frequency determines the frequency index of modulation. In cases when this is an integer, this index is pulse number which is characteristic for the synchronized methods of modulating function. Onother class of these methods is characterized by a continuous variation of the fundamental frequency for a constant frequency of the carrier signal. In this case, owing to nonperiodicity of pulses, appearance of subharmonics is possible having particularly harmful influence on the load at lower fundamental frequencies.
The basic characteristics of the signal at modulator output, irrespective of the method of edge determination, are the binary form and modulated width. They are transferred to the states of switching elements by means of the switching function.
Pulse edges at modulator output are most frequently determined by the natural or regular sampling method on the basis of the reference and carrier signals.
In case of a three-phase inverter modulation function may contain, in addition to the basic component, components of the order 3th which neutralize across the load. Application of the vector modulating function  extends the amplitude range.
2.1. Regular Sampling Method
In the regular sampling method the reference signal is discretized, and the pulse edges are determined by comparing the carrier signal with the reference signal modified in this way, which implies a digital realization of the method. Therefore, the regular sampling method has the advantage that microprocessor control is possible compared to the natural sampling method where the pulse edges should be determined by iterative solution of the transcendental equation obtained from the condition for crossing of the reference and carrier (triangular) signals. Figure 1 presents the synthesis of the output voltage of a half-bridge inverter by applying the modulating function method and regular sampling.
The interval of radians corresponding to the period of the modulating function will be divided in intervals of length , as shown in Figure 1. Denote the center of th interval by . At these points modulating function is discretized to a set of values which are compared with a triangular signal carrier of the frequency equal to the switching frequency of the inverter. The average value of the switching function of the inverter branch is
The pulse edges centered at are
Switching function of the half-bridge is fully determined if the angles of the pulse edges are known. From (2.4) it can be concluded that for complete determination of the switching function, it is necessary to determine pulse widths which correspond to the cycle of the switching frequency. Alternative voltage at the output is determined by means of switching function by the following expression:
The voltage waveform represents a copy of the periodic sequence of modulator output pulses. This voltage may be expressed using the Fourier series:
From expressions (3.12a)–(3.12d) it is possible to obtain the approximate solutions of transcendental equations (3.9). In view of the phase synchronization between the reference and carrier signal, it is sufficient to apply some of the expression times to obtain the th successive pulse edge angle. The remaining angles are determined by simply using the symmetry property from the natural sampling method, which has been preserved in expressions (3.12a)–(3.12d).Using the property of the symmetry of angles about the quarter-period of the reference signal, for an odd frequency index value, it is sufficient to apply the stated expressions and calculate the successive pulse edge angle.
3.2. Influence of Polynomial Degree to the Level of Approximation of Natural Sampling
The deviations of angles obtained using formulae (2.13), (2.14), (2.15), and (2.16) from the values found by applying the iterative procedure to expression (3.9) represent the measure of the remainder term . The maximum values of absolute remainder terms for for the amplitude index value varying from zero to one are shown in Table 1. As can be seen, polynomial approximations yield the same accuracy in determining the pulse edge angles as that achieved by applying some of the iterative numerical procedures (e.g., the Newton-Raphson or Gauss-Saidel methods).
However, unlike the mentioned numerical method, formulae (2.13), (2.14), (2.15), and (2.16) are directly applicable in the form of a microprocessor algorithm. The symmetry properties of the PWM sequence are preserved even for the low values of frequency index , and this permits the algorithm to be applied to the cases where the regular method shows a tendency toward degrading the harmonic and amplitude performances of inverter output voltage regulation. It can be seen from Table 1 that the maximum of remainder term is inversely proportional to the frequency index and directly proportional to the amplitude index.
3.3. Modified Polynomial Approximation by Using Chebyshev Polynomials
To achieve a higher accuracy of pulse edge with polynomial of lower order, we will apply the economization procedure by using Chebyshev polynomials . In this way we obtain the following polynomials:
When used as a PWM sampling algorithm, the modified polynomial permits computer time savings and yields a maximum error of 0.1297 degrees, which is comparable to the use of a nonmodified third-degree polynomial.
4. Analysis of the Results
Because of the odd symmetry of the ac component of the pulse sequence generated by the method proposed, the Fourier expansion of inverter ac voltage contains only sinusoidal components. The amplitude of the nth harmonic of voltage is given by expression:
After substituting the expressions for pulse edge angles, into (4.1) and then into (2.11), we find the current distortion factor as a function of the amplitude index, with the index figuring as a parameter. The regulation amplitude factor is obtained directly by applying formula (3.5) for , depending on the amplitude index . As the second-order polynomial algorithm gives relatively small deviations of pulse edge angles from those obtained using the natural method, the same conclusion applies to the amplitude of the fundamental harmonic generated by the second-order polynomial algorithm. The use of polynomial algorithm (3.14) produces a considerable decrease in the distortion factor in the upper part of amplitude modulation. The results obtained by numerical analysis and simulations clearly indicate the advantages of the proposed sampling method, at low-frequency index values, with respect to both the distortion factor and the maximum achievable modulation amplitude factor.
When the second-degree polynomial algorithms are applied to the sinusoidal modulating function, the results obtained for are already close to the ideal natural sampling method (Tables 3 and 4). At higher frequency index values, all algorithms converge to the same distortion factor values and to a common amplitude characteristic of voltage regulation, as confirmed by experiment.
All that has been mentioned above can be applied to the results for a three-phase inverter, since its bridge structure consists of three half-bridge structures. The proposed polynomial algorithms exhibit their major advantages when applied to the vector modulating function as well. At frequency index values , a suitably chosen algorithm may give the maximum values of amplitude characteristic, which are considerably better compared with the regular sampling method. In addition, the distortion factor is decreased. To illustrate this, let us say that the polynomial algorithm (3.13) for and yields the reduced amplitude of the fundamental harmonic of 0.5672 and a distortion factor DIS = 5.6761%. In case the regular sampling method is applied under identical conditions, the reduced amplitude is 0.5521 with a distortion factor DIS = 8.4087%. Algorithm (3.14) gives the reduced amplitude of 0.5773 and a distortion factor DIS = 10.0245%.
The odd values of frequency index imply the use of algorithm (3.14). As far as the amplitude characteristic is concerned, the best results for are achieved with synchronized half-periods of opposite signs (), whereas for half-periods of the same sign () should be synchronized. These results are presented in Figures 3 and 4.
5. Implementation of Polynomial Algorithm and Experimental Results
It is known that the basic characteristics of the digital implementation of the natural sampling are the iterative procedure and the work with the data in the floating point numbers format, which is a consequence of the transcendental nature of the expressions by which the angles of impulses edges are determined. Algebraic nature of the formula for the angles (2.13), (2.14), (2.15), (2.16), (3.13) and (3.14) give the possibility to replace the iterative procedure faster method of direct computing.
By transition to the time domain we find the formulas which determine the position of the th impulse inside the appropriate periods of switching frequency. The software implementation of a polynomial algorithm is achieved by simulation of the process which is shown in Figure 5 with the characteristic times determined by the following formulas which represents the times of intersection of carrier signal and reference signal, given in the explicit algebraic form of polynome with as the argument
In the formula (5.1) we suppose that where is the angle of the th sample of the modulation function determined in such a way that the polynomial algorithm becomes a regular algorithm with the characteristic times, given in the following formula:
The synthesis of impulses on the modulator output is performed by simulation of the process which is shown in the Figure 6. Detailed microprocessor implementation of regular symmetrical algorithm is given in [43–45].
Modern microcontrollers are equiped with integrated peripheral units designed for the synthesis of the modulated signals which control the state of switching of inverter branch. Intel’s microcontroller 8XC196MC/MD is optimized for implementation of PWM using the method of modulation function . It contains the generator of the modulated signal which significantly simplifies the programe implementation of selection algorithm and reduces the need for external hardware devices, enabling the synthesis of three independent complementary sets of modulated signals, with common switching frequency and conservatively determined “dead’’ time. Digital word which determines the period of switching frequency is defined by
whereis the value of the register of maximum frequency of output inverter's voltage; is the maximum value of frequency in (Hz); pulse number; is the clock frequency on the pin XTAL1 in (MHz).
The times of leading impulses edge are determined in the form of digital words which are calculated using
Similarly, we find the formula for digital words of the trailing edges of impulses, the value of which is set in the appropriate registers WG_COMPx when the counter WG_COUNT reaches the value 0001 H:
The values of the coefficient are normalized and are represented as double-complementary quantities in fixed point formats in the following way:
In this kind of representation all previously calculated coefficients are shown with the accuracy greater than in the range of changes from to within which are also all polynomial coefficients.
The coefficients of the leading and trailing edges are placed in the form of pairs of memory tables so that one table corresponds to the leading impulse edges and the other to the trailing impulse edges. The number of tables depends on the polynomial degree and the number of different values of the pulse number. The coefficients which correspond to the different phases are determined from the same table by applying the suitable index addressing of the data from the tables. The basic element which determines the conditions under which a polynomial algorithm can be applied is that the time necessary for realization of the program must be shorter than half of the carrier signal period. Concerning the memory table it is possible to perform certain optimizations using the properties of the symmetry of pulse angles about the half-period. The amplitude modulation index can be obtained also from the memory table which shows the output amplitude dependency on the frequency. In the experiment performed this dependency is linear, and it maps the ratio .
The experimental model (Figure 7) is made for the practical realization of the polynomial algorithm, and it contains three basic units.
Experimentally, the possibility of application of the polynomial algorithm as the modification of the regular selection algorithm for the switching frequency of 2,4 kHz is confirmed. For higher switching frequencies, the times needed for online computation of the positions of the impulse edges are critical regarding their asymmetry about the switching segment. The applications of the polynomial algorithm of the first degree for the output frequency of 40 Hz with the pulse number are achieved. Oscilloscope recordings of the line voltage and current, obtained by the synthesis according to polynomial algorithm with the vector modulation function, are shown in the Figure 8.
The suggested new sampling PWM algorithm by using the method of modulating function offers a possibility that natural sampling can be realized without any iterative procedure. It is shown that the angle of pulse edges can be approximated by a polynomial, the highest degree of the polynomial determining the measure of approximation of the obtained solution compared to natural sampling. The economizing procedure of using Chebyshev polinomials offered the possibility of using polinomials of a lower order, usually of the first- or second-order, still acomplishing a good approximation of the natural sampling. In terms of simplicity, this algorithm is comparable to a version of regular sampling method. As regards the performance concerning the quality of inverter’s output voltage, improvement by application of the suggested algorithm, which retains all good properties of the natural technique of PWM, is confirmed by simulation. In addition, possibility of implementation of polynomial algorithm in the real application is experimentally confirmed.
- A. H. Bonnett and D. J. Rhodes, “PWM inverter and motor applications: a quick reference of bibliographies and abstracts,” IEEE Transactions on Industry Applications, vol. 34, no. 1, pp. 217–221, 1998.
- Y. H. Kim and M. Ehsani, “An algebric algorithm for micro computer-based (direct) inverter pulse with modulation,” IEEE Transactions on Industry Applications, vol. 23, no. 4, pp. 654–659, 1987.
- S. Legowski and A. M. Trzynadlowski, “Incremental method of pulse width modulation for three-phase inverters,” in Proceedings of the 21st IEEE Industry Applications Society Annual Meeting (IAS '86), pp. 593–600, Denver, Colo, USA, 1986.
- A. M. Trzynadlowski, S. Legowski, and R. Kirilin, “Random pulse width modulation technique for voltage-controlled power inverters,” International Journal of Electronics, vol. 68, no. 6, pp. 1027–1037, 1990.
- J. Holtz, “Pulsewidth modulation for electronic power conversion,” Proceedings of the IEEE, vol. 82, no. 8, pp. 1194–1214, 1994.
- J. Holtz, “Pulsewidth modulation—a survey,” IEEE Transactions on Industrial Electronics, vol. 39, no. 5, pp. 410–420, 1992.
- L. M. Tolbert and T. G. Habetier, “Novel multilevel inverter carrier-based PWM method,” IEEE Transactions on Industry Applications, vol. 35, no. 5, pp. 1098–1107, 1999.
- J. M. Burdio, A. Martinez, and J. R. Garcia, “A synthesis method for generating switched electronic converters,” IEEE Transactions on Power Electronics, vol. 13, no. 6, pp. 1056–1068, 1998.
- B.-H. Kwon, T.-W. Kim, and J.-H. Youm, “A novel SVM-based hysteresis current controller,” IEEE Transactions on Power Electronics, vol. 13, no. 2, pp. 297–307, 1998.
- B. Mwinyiwiwa, Z. Wolanski, and B.-T. Ooi, “Microprocessor-implemented SPWM for multiconverters with phase-shifted triangle carriers,” IEEE Transactions on Industry Applications, vol. 34, no. 3, pp. 487–494, 1998.
- J. A. Houldsworth and D. A. Grant, “The use harmonic distortion to increase the output voltage of a three-phase PWM inverters,” IEEE Transactions on Industry Applications, vol. 20, no. 5, pp. 1224–1228, 1984.
- M. A. Boost and P. D. Ziogas, “State-of-the-art carrier PWM techniques: a critical evaluation,” IEEE Transactions on Industry Applications, vol. 24, no. 2, pp. 271–280, 1988.
- A. M. Trzynadlowski, “Nonsinusoidal modulating functions for three-phase inverters,” IEEE Transactions on Power Electronics, vol. 4, no. 3, pp. 331–338, 1989.
- D.-C. Lee and G.-M. Lee, “A novel overmodulation technique for space-vector PWM inverters,” IEEE Transactions on Power Electronics, vol. 13, no. 6, pp. 1144–1151, 1998.
- A. M. Hava, R. J. Kerkman, and T. A. Lipo, “Carrier-based PWM-VSI overmodulation strategies: analysis, comparison, and design,” IEEE Transactions on Power Electronics, vol. 13, no. 4, pp. 674–689, 1998.
- H. S. Patel and R. G. Hoft, “Generalized techniques of harmonic elimination and voltage control in thyristor inverters —I: harmonic elimination,” IEEE Transactions on Industry Applications, vol. 9, no. 3, pp. 310–317, 1973.
- G. S. Buja and G. B. Indri, “Optimal pulse width modulation for feeding ac motors,” IEEE Transactions on Industry Applications, vol. 13, no. 1, 1977.
- G. Franzo, M. Mazzucchelli, L. Puglisi, and G. Sciutto, “Analysis of PWM techniques using uniform sampling in variable-speed electrical drives with large speed range,” IEEE Transactions on Industry Applications, vol. 21, no. 4, pp. 966–974, 1985.
- T. Kato, “Precise PWM waveform analysis of inverter for selected harmonic elimination,” in Proceedings of the 21st IEEE Industry Applications Society Annual Meeting (IAS '86), pp. 611–616, Denver, Colo, USA, 1986.
- I. Takahashi and H. Mochikawa, “A new control of PWM invereter waveform for minimum loss operation of an induction motor drive,” IEEE Transactions on Industry Applications, vol. 21, no. 3, pp. 580–587, 1985.
- M. Grötzbach and R. Redmann, “Line current harmonics of VSI-fed adjustable-speed drives,” IEEE Transactions on Industry Applications, vol. 36, no. 2, pp. 683–690, 2000.
- S. R. Bowes, “Developments in PWM switching strategies for microprocessor-controlled inverter drives,” in Proceedings of the 22nd IEEE Industry Applications Society Annual Meeting (IAS '87), pp. 323–329, Atlanta, Ga, USA, 1987.
- Z. Yu, A. Mohammed, and I. Panahi, “A review of three PWM techniques,” in Proceedings of the American Control Conference (ACC '97), pp. 257–261, Alberquerque, NM, USA, June 1997.
- V. R. Stefanovic and S. N. Vukosavic, “Space-vector PWM voltage control with optimized switching strategy,” in Proceedings of IEEE Industry Applications Society Annual Meeting (IAS '92), pp. 1025–1033, 1992.
- A. M. Massoud, S. J. Finney, and B. W. Williams, “Control techniques for multilevel voltage source inverters,” in Proceedings of the 34th IEEE Annual Power Electronics Specialists Conference (PESC '03), vol. 1, pp. 171–176, Acapulco, NM, USA, June 2003.
- B. K. Bose, Modern Power Electronics and AC Drives Englewood Cliffs, Prentice-Hall, Englewood Cliffs, NJ, USA, 2002.
- D. G. Holmes, “The general relationship between regular-sampled pulse-width-modulation and space vector modulation for hard switched converters,” in Proceedings of IEEE Industry Applications Society Annual Meeting (IAS '92), pp. 1002–1009, 1992.
- S. R. Bowes and M. J. Munt, “Microprocessor control of PWM inverters,” IEE Proceedings B: Electric Power Applications, vol. 128, no. 6, pp. 293–305.
- G. Narayanan and V. T. Ranganathan, “Extension of operation of space vector PWM strategies with low switching frequencies using different overmodulation algorithms,” IEEE Transactions on Power Electronics, vol. 17, no. 5, pp. 788–798, 2002.
- J. T. Boys and P. G. Handley, “Harmonic analysis of space vector modulated PWM waveforms,” IEE Proceedings B, vol. 137, no. 4, pp. 197–204, 1990.
- S. R. Bowes and Y.-S. Lai, “The relationship between space-vector modulation and regular-sampled PWM,” IEEE Transactions on Industrial Electronics, vol. 44, no. 5, pp. 670–679, 1997.
- A. M. Trzynadlowski, R. L. Kirilin, and S. F. Legowski, “Space vector PWM technique with minimum switching losses and a variable pulse rate,” IEEE Transactions on Industrial Electronics, vol. 44, no. 2, pp. 173–181, 1997.
- K. Zhou and D. Wang, “Relationship between space-vector modulation and three-phase carrier-based PWM: a comprehensive analysis,” IEEE Transactions on Industrial Electronics, vol. 49, no. 1, pp. 186–196, 2002.
- A. Kwasinski, P. T. Krein, and P. L. Chapman, “Time domain comparison of pulse-width modulation schemes,” IEEE Power Electronics Letters, vol. 1, no. 3, pp. 64–68, 2003.
- D. G. Holmes and T. A. Lipo, Pulse Width Modulation for Power Converters—Principle and Practice, IEEE Series on Power Engineering, IEEE Press/Wiley InterScience, New York, NY, USA, 2003.
- I. Deslauriers, N. Avdiu, and B. T. Ooi, “Naturally sampled triangle carrier PWM bandwidth limit and output spectrum,” IEEE Transactions on Power Electronics, vol. 20, no. 1, pp. 100–106, 2005.
- A. M. Trzynadlowski and S. Legowski, “Minimum-loss vector PWM strategy for three-phase inverters,” IEEE Transactions on Power Electronics, vol. 9, no. 1, pp. 26–34, 1994.
- F. Blaabjerg, J. K. Pedersen, and P. Thoegersen, “Improved modulation techniques for PWM-VSI drives,” IEEE Transactions on Industrial Electronics, vol. 44, no. 1, pp. 87–95, 1997.
- A. M. Hava, Carrier based PWM voltage source inverter in overmodulation range, Ph.D. thesis, University of Wisconsin, Madison, Wis, USA, 1998.
- D. Kostic, New method for the optimisation PWM voltage source inverters, Ph.D. thesis, University of Belgrade, Belgrade, Serbia, 2001.
- G. M. Fihtengolc, Differential and Integral Calculs, vol. 2, Nauka, Moscow, Russia, 1966.
- L. Fox and J. B. Parker, Chebyshev Polynomials in Numerical Analysis, Oxford University Press, London, UK, 1972.
- S. R. Bowes and A. Midoun, “Suboptimal switching strategies for microprocessor-controlled PWM inverter drives,” IEE Proceedings B, vol. 132, no. 3, pp. 133–148, 1985.
- B. K. Bose and H. A. Sutherland, “A high-performance pulse-width modulator for an inverter-fed drive system using a microcomputer,” IEEE Transactions on Industrial Electronics, vol. 19, no. 2, pp. 235–243, 1983.
- D. Vincenti, P. D. Ziogas, and R. V. Patel, “A PC-based pulse-width modulator for static converters,” IEEE Transactions on Industrial Electronics, vol. 37, no. 1, pp. 57–69, 1990.
- 8XC196∗∗ Microcontroler-User manual, INTEL 1996.
Copyright © 2009 Dragutin Kostić and Vladimir Šinik. 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.