Discrete Dynamics in Nature and Society

Volume 2010, Article ID 408346, 12 pages

http://dx.doi.org/10.1155/2010/408346

## Parameter Identification of a Class of Economical Models

Instituto de Matemática Multidisciplinar, Universidad Politécnica de Valencia, 46022 Valencia, Spain

Received 8 October 2009; Revised 13 May 2010; Accepted 19 May 2010

Academic Editor: Aura Reggiani

Copyright © 2010 Begoña Cantó 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.

#### Abstract

In mathematical modeling of economic problems, it is common to assume that economic laws influence the process and analyze the type of solution obtained. Sometimes these laws have unknown parameters and they worth studying if these parameters can be determined. In this paper, a structured system is considered and the identifiability is analyzed. In case the parameters can be uniquely identified, an algorithm for obtaining the model parameters is presented and, finally, the existence of a balanced growth solution is studied.

#### 1. Introduction

Linear ordinary differential or difference equations are frequently used in modeling many engineering, economic, and social processes. These equations describe the system behavior with a relative degree of accuracy and yet are simple enough to be solved. It is well known that any differential or difference equation of order can be written as a system of first-order differential or difference equations. This process works for any higher-order equation, linear or not, and it consists of expressing the top derivative or difference as a function of the lower ones. An order equation gives rise to a first order system in variables whose coefficient matrix is called companion matrix. The characteristic polynomial of an order linear equation is the same as the characteristic polynomial of the companion matrix. Then, a companion matrix is associated with the monomial characteristic polynomial

In the last years, the use of dynamical models in economic theory is quite usual, in stark contrast, with the classical static economic models. Many of these models are mathematically formulated by higher-order difference equations (see [1]), which derive to a structured representation, where the companion matrix can be used.

As already mentioned, economic problems are usually modeled by structured systems, that is, systems where the matrices have a number of fixed zero entries while the rest of the entries are not known. This kind of system is widely considered in the literature, than in control theory (see [2, 3]), as in computing methods to solving it. When the size of the system is large, for solving these structured systems, fast algorithms are used. In [4] an exhaustive analysis on numerical properties of several fast algorithms was shown. In [5] an algorithm for solving dynamic economic models developing first-order and second-order iterations was studied.

In structured systems, the analysis of the unknown entries plays an important role. For example, the economic models shown in [1] have a parameterized structure, where the parameters are given information on the economic process. A first step for the validation of the model is to identify the parameters that appear in the system. Before identifying the parameters, it should be examined whether the parameters can be determined from the data. The identification of parameters allows solving the model. And finally, it is necessary to analyze the obtained solution. In particular, in economic models, the existence of a balanced growth solution is interesting. In the economical sense, balanced growth refers to simultaneous coordinated expansion of several sectors.

The aim of this paper is to study the identifiability problem and the balanced growth solution of an economic problem. The structure of the paper is as follows. In Section 2 we treat the identification problem of a companion structured system. In Section 3, we show a structured Leontief model associated to a determinate economical process and give some results about the identifiability of its parameters. This analysis leads us to construct an algorithm to obtain the parameters of the model. Finally, in the last section we study the balanced growth solution of an economic model.

#### 2. Identification Problem of a Structured System

In the first part of this paper, we discuss a problem of identifiability for a linear differential or difference system or where the coefficient matrices have a fixed structure, is a companion matrix and is a monomial vector given by with the parameter vector belonging to a subset Note that, the matrix is a companion matrix where the last row is the vector .

In general, an interpretation of a structured linear system is a family of linear systems, with the same structure, together a set of parameters. In this work, this kind of systems is denoted by

A formal definition of structured linear dynamic systems can be found in [6]. Furthermore, denotes the input-output behavior of the system . Note that, of the system can completely be characterized by the Markov parameters, defined as

Usually, in the analysis of parametric models, studied two problems can be: identifiability and estimability. The difference between estimability and identifiability is that in the identifiability all possible options of the input-output signal are evaluated, while estimability is the ability to estimate accurately the parameters of a given data set.

A model is identifiable if and only if, for every parameter set there is a unique input-output behavior. Several papers deal with identification of models describing behavior of biomedical, chemical, or other systems, [7, 8].

The problem of the structural identifiability of the model consists of the determination of all parameter sets which give the same input-output structure. The structural identifiability of depends on the number of solutions of . If the parameters can be determined uniquely from the data, that is, the equation has only the solution , then the system is *structurally globally identifiable*. On the other hand, if there exists a finite or contable number of solutions for the parameters, the system is *structurally locally identifiable*. Otherwise, the system is *structurally unidentifiable* ( see [9, 10]).

Now, we consider the system , with and given by (1.1) and (2.3). In order to analyze the identification problem of this system, we use the Markov parameters to know its input-output behavior . By the structure of the companion matrix and by a technical process, we have the following result.

Proposition 2.1. *Consider the structured system The Markov parameters , are given by
*

Directly from the above result, we obtain the following proposition which gives a condition to solve the identification problem for the system .

Proposition 2.2. *The structured system is globally identifiable.*

In the next section, we analyze the identification problem for economic models.

#### 3. Economic Dynamic Model

First, we consider the nonlinear discrete-time dynamic model given in [1] where , , is the income at time . Parameters used in this model have the following economic interpretation: the parameter represents the autonomous expenditures, the parameter , and represents the consumer's reaction against the increase or decrease of his income. Finally, the investment is described using parameters and , where represents the difference between the control and the investment function.

This model is interesting because for different values of parameters we obtain well-known linear dynamic models as Hick-Samuelson, Keynes, and the nonlinear dynamic Puu model (see [11]).

In particular, for , , with and , , we obtain the Hick-Samuelson model which matrix expression is given by the companion matrix If we do not have investment, that is, , we obtain the Keynes model which has the following coefficient matrix

Note that these models are structured systems of the kind where the coefficient matrix depends on a vector of parameters which could be identified to analyze its solution. Then, by Proposition 2.2 these models are globally identifiable.

Now, we consider the Leontief model. In the Leontief economic process is assumed that there is a market for different goods in which each industry has only one production process. Assuming that is the production level vector and is the demand level vector, the control process (see, e.g., [12, 13]), the system has the structured representation given by where is the capital coefficient matrix, is the technological coefficient matrix, and is the demand coefficient matrix (excluding investment). The -entry of the matrices and represent the amounts of material input and capital of the good necessary for the production of one unit of the industry, respectively. Since the Leontief model is an economic model, and This kind of systems is a dynamic generalized system because the capital coefficient matrix can be or not be singular. The singularity of this matrix arises because no output from one sector is used in the production of some products.

In this paper, the system (3.6) is a structured linear dynamic system where the capital matrix is a nonsingular diagonal matrix, , and the matrix has a companion matrix structure. By the nature of the technological matrix and where the unknown demand with almost one column of type (2.3)exists.

Since capital matrix is invertible, the structured generalized system can be transformed into a structured standard system where is such that it has almost one monomial column of kind (2.3), with .

An important point to note here is the structure of matrix . In the literature, this kind of matrix is known as the modified Leslie (Lefkovitch) matrix, which appears in different areas such as the study of recruitment, survival and population growth rate (see, e.g., [14, 15]). For example, in [16] it is used to describe a discrete age-structured population model where no emigration during harvesting is considered.

Some properties of the spectrum of are given in the following result, which are useful to study the stability of the system .

Proposition 3.1. *Consider the matrix given in (3.9). *(i)*If or , with , then *(ii)*If or , then *

*Proof. *If we structure as , with , , , and , then
(i)If or , with , then we simplify
where is a polynomial of degree .Hence, .(ii)By (i) it is straightforward.

In order to solve the identifiability problem of a family of structured systems of kind (3.6), we show some results on structured standard system .

The system (3.9)-(3.10) belongs to a kind of structured systems. In particular, the identifiability of the parameters of a standard system is concerned with the determination of them from the external behavior of the system. That is, to determine the input-output behavior of a model , we can use the Markov parameters associated to the system .

The parameter identification process for the structured system is followed from the structure of the vectors obtained in the following proposition.

Proposition 3.2. *Consider the structured system The Markov parameters , are given by
*

We solve the identification problem for the system in the next result.

Proposition 3.3. *The structured system where and are defined by (3.9) and (3.10) with is globally identifiable.*

*Proof. *We consider two structured systems and defined by (3.9)-(3.10) with input-output behaviorand we shall prove that By the structure of the Markov parameters obtained in Proposition 3.2, for it is show that
Moreover, the rest of unknown entries of are identified from the first equality . Hence, .

Consider the associated standard Leontief model (3.8), in this case given by and has the column By Proposition 2.2 the entries of these matrices are identified and hence we can also identify the initial parameters of the vector . Thus, the initial Leontief model is globally identifiable.

Proposition 3.4. *The Leontief economic model given by (3.6) where is a non-singular diagonal matrix, and is a companion matrix, is globally identifiable.*

The following algorithm allow us determine the parameters of the Leontief model.

*Algorithm**Step 1. *Introduce the size of the state vector: . Introduce the capital matrix . Introduce the matrices that determine the known external behavior of the system with the structure (3.9)-(3.10). And introduce the position of the monomial column of matrix : , for . *Step 2. *Choose the th column of and denote it as .*Step 3. *Introduce .*Step 4. *Construct the following system:
*Step 5. *Solve the above system and obtain the parameters for and and , for , and .*Step 6. *Use the parameters obtained in Step 5 and construct the matrices (3.17)-(3.18) and from them obtain and , , and .

With this process, we have determinate all parameters of the system.*Step 7. *Finally, construct the matrices of the system (3.6).Note that, this algorithm can also be applied when the Hick-Samuelson and the Keynes models are considered, that is, when and , .

Next we present an academic example to clarify the above algorithm. The model considers a number of consumers such that at any given time each consumer holds quantities of different services or activities. We assume that this kind of market can be modeled by a structured Leontief model. Using the data, the problem is to assure the uniqueness in the parameter of the dynamical model. In this example to determine these parameters we use the above algorithm.*Example 3.5. *Tourism demand in a country involves different activities. In this example the external behavior of the foreign demand for tourist services during a period of years is known and we want to apply the algorithm in order to identify the structured model.

The external behavior of the foreign demand for tourist services during a period of four years is given by
and consider that the capital coefficient matrix (given in thousands)

Then we apply the above algorithm.*Step 1. * Introduce , the matrices , and , for .*Steps 2-3*

Denote as and introduce .*Steps 4-5*

Solving the system
obtain
*Steps 6-7*

From (3.17)-(3.18) the system is identifiable and its coefficient matrices are

An interpretation of these results can be, for example the following. It seems reasonable that the fact that capital matrix is diagonal implies that to obtain one unit of production in each industry, an investor is only necessary. On the other hand, the technological matrix is determined as each service depends on one material input except the last service, which needs a part of all material input involved in the process. Finally, the structure of the demand matrix indicates that only the last service has consumers, given in thousands. This service includes all others.

#### 4. Balanced Growth Solution

In economic problems it is important to assure the existence of a solution which output of each sector increases by a constant percentage per unit of time. This kind of solution is called balanced growth solution, and it satisfies with the balanced growth rate .

The balanced growth problem is solved in [17] for autonomous systems. The interest in obtaining balanced growth solutions for some dynamic model has to do the following question: is it possible to obtain a feedback such that the closed-loop system has a balanced growth solution? A first approximation of the solution to this problem has been given in [18]. In this work a characterization to the existence of a balanced growth solution has been obtained. This characterization is based on the existence of a feedback such that the closed-loop system has a positive eigenvalue with a positive eigenvector.

To analyze the balanced growth solution of the system (3.6), we consider the associated system given by

The balanced growth problem for Leontief structured model is treated in the following result.

Proposition 4.1. *The balanced growth problem for Leontief structured model (3.6) is solvable.*

*Proof. *Note that technological matrix is an irreducible matrix, that is, its associated directed graph is strongly connected, which follows by the companion matrix structure of . Moreover, as and , we can assure that . These conditions on lead us to the fact that matrix is positive. If we construct a feedback with , then . Therefore, the coefficient matrix of the closed-loop system is positive. Hence, there exists a positive eigenvalue with a positive eigenvector , that is,
and taking we have

Considering the state-feedback , we show that is a solution of the system (3.6). Thus, the balanced growth problem for Leontief structured model (3.6) is solvable.

The economic meaning of the existence of a positive eigenvector can be given, for instance, when there is one economy where each sector of this economy depends on all others directly or indirectly for either its intermediate product or its capital. We refer to [17] for more details for it and in the case when the economy can be divided into several sub-conomies.

#### Acknowledgments

The authors are grateful to the anonymous referees for their constructive remarks. This paper is supported in part by Grant MMT2007-64477.

#### References

- L. I. Dobrescu, M. Neamtu, A. L. Ciurdariu, and D. Opriş, “A dynamic economic model with discrete time and consumer sentiment,”
*Discrete Dynamics in Nature and Society*, vol. 2009, Article ID 509561, 18 pages, 2009. View at Publisher · View at Google Scholar - J.-M. Dion, C. Commault, and J. van der Woude, “Generic properties and control of linear structured systems: a survey,”
*Automatica*, vol. 39, no. 7, pp. 1125–1144, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - T. Boukhobza, F. Hamelin, and D. Sauter, “Observability of structured linear systems in descriptor form: a graph-theoretic approach,”
*Automatica*, vol. 42, no. 4, pp. 629–635, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. P. Brent, “Stability of fast algorithms for structured linear systems,” in
*Fast Reliable Algorithms for Matrices with Structure*, T. Kailath, Ed., pp. 103–116, SIAM, Philadelphia, Pa, USA, 1999. View at Google Scholar · View at MathSciNet - A. Ludwig, “The Gauss-Seidel-quasi-Newton method: a hybrid algorithm for solving dynamic economic models,”
*Journal of Economic Dynamics & Control*, vol. 31, no. 5, pp. 1610–1632, 2007. View at Publisher · View at Google Scholar · View at MathSciNet - J. M. van den Hof, “Structural identifiability of linear compartmental systems,”
*IEEE Transactions on Automatic Control*, vol. 43, no. 6, pp. 800–818, 1998. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - L. G. Hanin, “Identification problem for stochastic models with application to carcinogenesis, cancer detection and radiation biology,”
*Discrete Dynamics in Nature and Society*, vol. 7, no. 3, pp. 177–189, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - B. R. Jayasankar, A. Ben-Zvi, and B. Huang, “Identifiability and estimability study for a dynamic solid oxide fuel cell model,”
*Computers and Chemical Engineering*, vol. 33, no. 2, pp. 484–492, 2009. View at Publisher · View at Google Scholar · View at Scopus - É. Walter and L. Pronzato,
*Identification of Parametric Models from Experimental Data*, Communications and Control Engineering Series, Springer, Berlin, Germany, 1997. View at MathSciNet - A. Ben-Zvi, P. J. McLellan, and K. B. McAuley, “Identifiability of linear time-invariant differential-algebraic systems. I. The generalized Markov parameter approach,”
*Industrial and Engineering Chemistry Research*, vol. 42, no. 25, pp. 6607–6618, 2003. View at Google Scholar · View at Scopus - T. Puu and I. Sushko, “A business cycle model with cubic nonlinearity,”
*Chaos, Solitons & Fractals*, vol. 19, no. 3, pp. 597–612, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - B. Cantó, C. Coll, and E. Sánchez, “Positive $N$-periodic descriptor control systems,”
*Systems & Control Letters*, vol. 53, no. 5, pp. 407–414, 2004. View at Publisher · View at Google Scholar · View at MathSciNet - M. S. Silva and T. P. de Lima, “Looking for nonnegative solutions of a Leontief dynamic model,”
*Linear Algebra and Its Applications*, vol. 364, pp. 281–316, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. R. Skalski, K. E. Ryding, and J. J. Millspaugh,
*Wildlife Demography*, Elsevier Academic Press, London, UK, 2005. - L. Rogers-Bennett and D. W. Rogers, “A semi-empirical growth estimation method for matrix models of endangered species,”
*Ecological Modelling*, vol. 195, no. 3-4, pp. 237–246, 2006. View at Publisher · View at Google Scholar · View at Scopus - D. Sadhukhan, B. Mondal, and M. Maiti, “Discrete age-structured population model with age dependent harvesting and its stability analysis,”
*Applied Mathematics and Computation*, vol. 201, no. 1-2, pp. 631–639, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - L. Zeng, “Some applications of spectral theory of nonnegative matrices to input-output models,”
*Linear Algebra and Its Applications*, vol. 336, pp. 205–218, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - B. Cantó, C. Coll, and E. Sánchez, “Positive solutions of a discrete-time descriptor system,”
*International Journal of Systems Science*, vol. 39, no. 1, pp. 81–88, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet