Research Article  Open Access
An Iterative Decomposed PiecewiseLinear Model Description
Abstract
A model description for the representation of onedimensional piecewiselinear characteristics is presented. The model can be denoted as a decomposed one, because the independent and dependent variables of the PWL characteristic are treated separately. It is also called iterative, because the particular representation of each segment of the PWL characteristic depends on the value of only one parameter included in the mathematical formulation, it gives the possibility of modeling both, univalued and multivalued PWL characteristics.
1. Introduction
In order to be able to simulate nonlinear networks, the nonlinear behavior of the components must be modeled first. The nonlinear behavior cannot exactly be described by familiar analytical nonlinear functions, so it must be approximated by a known function. Some approaches to overcome this problems are the polynomial approximation, tablelookup technique, and PWL approximation.
A lot of research has been done in the area of approximating nonlinear functions using polynomials as well as in the field of calculating transcendental functions. Applying such methods leads to simulators like SPICE [1].
The technique alá tablelookup consisting in storing several sample points of a nonlinear function in a table. If a function value is needed, the entry in the table is used or interpolated.
In PWL approximation, the nonlinear function is repeatedly approximated locally by a linear mapping. This results in a large collection of linear mappings which can be stored in a compact closed form model [2, 3]. These models have been applied in the socalled PWL simulators [4–6].
There are many important proposals in the area of PWL modeling, but no best model description can be found because the advantages or disadvantages of any PWL model depends on its application [7].
The main contributions of the PWL model, which in Section 2 will be exposed, are the possibility of describing all kinds of onedimensional PWL characteristics (univalued and multivalued) having mathematical formulation which is not so heavy to be computed because it is based on an iterative system of linear equations, the advantage of obtaining the model parameters directly from the graphic coordinates and the simulation capability of the model, because it is compatible with iterative simulation techniques [8].
The model is denoted as decomposed, because the independent variable and the dependent variable are included separately in a system of linear equations. The model is also denoted as iterative, because there is a parameter , which selects a specific segment from the set of segments that constitutes the complete PWL characteristic.
2. Model Description
The PWL characteristic, depicted in Figure 1, represents a general onedimensional function characterized by linear segments and coordinates: .
An iterative decomposed formulation for this characteristic can be obtained as follows.
Firstly, the particular problem of obtaining the line equation for any of the kth segments included in the PWL characteristic is analyzed. The graphic exposition of this problem is depicted in Figure 2.
The line equation for the kth segment is given as Then, (1) is factored into the form where and are variables defined by the determinant relations: The relation between (2) and (3) is expressed as the system of linear equations given in The above linear system represents the kth line equation included in the PWL characteristic to be modeled. Finally, a complete mathematical expression which assures a description for the PWL characteristic depicted in Figure 1 is obtained if (4) is generalized for all segments. Thus, the parametric model description can be written: with the condition where is denoted as activation vector and it is defined as follows: for , and for .
The variables and in (5) have an important function because they control the interval where a kth line equation exists. This is expressed in the condition (6) and it is demonstrated as follows.
Firsty, the notations and are defined as Then, mathematical expressions, for the variables and , can be obtained if the following decomposed linear equations, which are included in (5), are algebraically manipulated:
Isolating from (12) and substituting it into (10), (11), can be expressed as
A similar analysis where is isolated from (12) and then susbstituted into (10), (11) yields
If (13) and (14) are subjected to the exposed condition in (6) results The inequations (15) indicate the intervals where a linear segment belonging to a PWL characteristic exists. The shadow area depicted in Figure 2 results by the intersection of these intervals for the kth PWL segment.
3. Electrical Model
The linear system in (5) can be electrically modeled by the circuit shown in Figure 3. In this circuit, the variables and have a physical meaning, and represent the current flow through the conductances whose values are and , respectively. The electrical diagram includes two switches which are controlled simultaneously by . Figure 4 shows the symbol used to represent them. On the one hand, it can be seen that when the switches are open, these variables are equal to zero, it occurs when , on the other hand, when the switches are closed and the variables (currents and ) are different to zero.
4. Example
Consider the theoretical nonlinear element and its PWL characteristic depicted in Figure 5. Determine their iterative decomposed model description.
(a)
(b)
An analysis from Figure 5(b) shows the existence of three linear equations labeled as (1), (2), and (3). The four graphic coordinates are designated as B1, B2, B3, and B4, and summarized in Table 1.

There are three segments, therefore and . Substituting and the coordinates of Table 1 into (5), results: Substituting in (7) and (8), the terms of the activation vector are given by Figure 6 shows the electrical circuit which models the PWL characteristic analyzed in this example. Each th Block is a simplified representation of the circuit shown in Figure 3.
When , , , , and the resulting linear system is The line equation is bounded by the intervals and .
When , , , , and the resulting linear system is The line equation is bounded by the intervals and .
When , , , , and the resulting linear system is
The line equation is bounded by the intervals and .
The above results are summarized in Table 2 and Figures 7, 8, and 9.
5. Conclusions
The proposed iterative and decomposed model can describe any arbitrary PWL function (univalued or multivalued) by considering as model parameters its graphic coordinates. The interval in the model bounds the range where a linear segment is valid and it is determined by the variables and . These variables have a physical meaning (currents) in the electrical model. The equivalent circuit, of the proposed PWL description, has simulation capability because it is compatible with iterative simulation techniques.
Acknowledgments
Victor JimenezFernandez is holder of a scholarship from CONACyTMéxico by the Contract 183249. This work has been partially supported by a CONACyTMéxico research project under Grant 60261.
References
 A. Vladimirescu and S. Liu, “The simulation of MOS integrated circuits SPICE2,” Memorandum no. UCB/ERL M80/7, Electronics Research Laboratory, College of Engineering, University of Berkeley, Berkeley, Calif, USA, February 1980. View at: Google Scholar
 L. O. Chua and S. M. Kang, “Sectionwise piecewiselinear functions: canonical representation, properties, and applications,” Proceedings of the IEEE, vol. 65, no. 6, pp. 915–929, 1977. View at: Google Scholar
 W. M. G. van Bokhoven, Piecewise Linear Modeling and Analysis, Kluwer Technische Boeken, Deventer, The Netherlands, 1981.
 J. T. J. van Eindhoven, A piecewise linear simulator for large scale integrated circuits, Ph.D. dissertation, Eindhoven, The Netherlands, 1984.
 T. A. M. Kevenaar and D. M. W. Leenaerts, “A flexible hierarchical piecewise linear simulator,” Integration, the VLSI Journal, vol. 12, no. 2, pp. 211–235, 1991. View at: Publisher Site  Google Scholar
 M. T. van Stiphout, PLATO: a piecewiselinear analysis tool for mixed level circuit simulation, Ph.D. dissertation, Technische Universiteit Eindhoven, Eindhoven, The Netherlands, 1990.
 T. A. M. Kevenaar and D. M. W. Leenaerts, “A comparison of piecewiselinear model descriptions,” IEEE Transactions, CAS1: Fundamental Theory and Applications, vol. 39, no. 12, pp. 996–1004, 1992. View at: Google Scholar
 L. O. Chua, “Analysis and synthesis of multivalued memoryless nonlinear networks,” IEEE Transactions on Circuit Theory, vol. 14, no. 2, pp. 192–209, 1967. View at: Google Scholar
Copyright
Copyright © 2009 Victor JimenezFernandez et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.