Abstract

This paper presents the use of an iteration method to solve the identifiability problem for a class of discretized linear partial differential algebraic equations. This technique consists in replacing the partial derivatives in the PDAE by differences and analyzing the difference algebraic equations obtained. For that, the theory of discrete singular systems, which involves Drazin inverse matrix, is used. This technique can also be applied to other differential equations in mathematical physics.

1. Introduction

Mathematical models are used extensively to simulate the behavior of experimental and practical situations. These models can be constructed in different ways and differ greatly in complexity but they have in common that they can be used to predict the behavior of the initial system. Usually, the model is constructed using some of the properties of a system by tinkering with adjustable parameters and sometimes the equations of the models involved unknown parameters. In this case, it is interesting to obtain these parameters to make the model accurate. This property is known as identifiability problem. Given a parameterized state-space model, structural identifiability is concerned with whether the unknown parameters within the model can be identified uniquely from the experiment considered. Thus structural identifiability analysis is an important step in the modeling process, and it is necessary for theoretical prerequisites to experiment design and system identification, that is, to estimate the unknown parameters of the model using experimental data. A wide variety of these models is often defined by differential algebraic equations (DAEs) or partial differential algebraic equations (PDAEs). These equations relate internal variables and their derivatives with inputs to the system and appear in several scientific fields. Recently, several researches have been devoted to PDAEs, and some theoretical and numerical results are obtained [14].

In this paper, we limit ourselves to solve the identifiability problem for second order partial differential equations which has added additional algebraic equations. The relations can be described by the general explicit formulation (see [1, 5]) using a matrix equation in the form where is a smooth function, are matrices in , all of which can be singular, and is a unknown parameter vector. The solution of (1.1), denoted by , depends on the parameter vector , the spatial variable and the time value .

For system (1.1) initial values and boundary conditions may depend on the parameter vector to be estimated. In general, if an arbitrary initial condition is imposed, it may be inconsistent with the PDAE; for that, we assume that initial values may be decomposed into the data which can be prescribed arbitrarily and the consistent data. The boundary values are in similar form.

The outline of this paper is as follows. Problem description and general formulations are discussed in Section 2, also it contains the numerical procedure to discretize the partial differential algebraic equation. Section 3 contains the solution to the identifiability problem of a structural discrete system and the application of the obtained results to solve the initial identifiability problem. For that, an explicit expression of the Drazin inverse of a structured matrix with parameters is given and the structure of the Markov parameters of the singular system is obtained. Finally, in last the section, an algorithm and an illustrative example are shown.

2. Problem Description and Discretization

In this section we consider the second-order partial differential algebraic equations (PDAEs) given by where , , and , , , and , . The differential algebraic aspect of the system is due to singularity of .

The basic idea to solve the identifiability problem is the use of an iteration method based on the numerical procedure that is useful to discretize partial differential algebraic equations. This technique consists in replacing the partial derivatives in the PDAE by differences. For that, we discretize our domain, both time and space can be discretized by an uniform grid with gridparameter , where is the distance between two neighboured nodes of the grid, and we discretize the interval by a one-dimensional grid with stepsize .

In this case we use the forward-looking difference operator in time and the approximation is and the second-order derivative is approximated by the second-order central difference operator in space given by

To obtain the discretized system we take . The discretized problem can by written as where and .

Moreover, we need to add some equations to show the transition among the variables. These new equations are linearly independent of the other equations and have the form

To obtain a structure which simplifies the solution we write the full system of equations as follows: where the state vector with and , the control vector , and the coefficient matrices are given by with . These matrices have a fixed structure with the parameter vector belonging to a subset , which suggests the relation among variables in the process. It is well-known (see [6]) that the system has solution if there exists an scalar such that (regularity condition). The matrix is a symmetric, tridiagonal, and regular matrix. This follows from the eigenvalues of given by and using the expressions of and , since it is easy to check that This establishes the regularity of and this fact guaranties the regularity of the system (2.6).

3. The Identification Problem

The system (2.6) is a singular discrete-time system which, from the regularity of , we can rewrite as follows: where and . This system is denoted by . Since there exists a scalar such that , the explicit solution of the system is given by where denotes the Drazin inverse of matrix , is the index of , and is an admissible initial condition, , with where , .

The input-output behavior of the system is obtained from the Markov parameters associated which the system. These parameters are defined as

3.1. Identifiability of Structural Singular Systems

To get a good formulation of the model we need to know if the unknown parameters can be determined uniquely from the obtained experiment data. This is known as identifiability problem.

Our system (2.6) has a fixed structure, and we want to solve the identifiability problem when this structure holds on. The identifiability helps us test the unique relationship between parameter sets and model response and guarantees that the parameters can be estimated under ideal conditions.

A model is globally structurally identifiable if and only if there exists a unique input-output behavior for every parameter set [3]. That is, the system is globally structurally identifiable if and only if, for almost any two candidates parameter vector values , imply , where denotes the input-output behavior of the system . From now on, we use the concept globally identifiable as globally structurally identifiable.

For linear models there are many well established techniques to analyze structural identifiability; see, for example, [4, 7, 8] and the references therein. The main contribution of this part consists in using the special structure of the Markov parameters associated with to the identifiability analysis. For that, we need to know the structure of .

By definition of to obtain we use the inverse of the tridiagonal matrix. A concise expression of is obtained as follows: where and , for , with and . Note that , where is the principal minor of order of the matrix and the symmetry property of matrix is transferred to . For more information on the inverse of a tridiagonal matrix see [9, 10].

Using the above inverse, an explicit expression of the Drazin inverse of is given in the following result.

Proposition 3.1. Consider the matrices given in (2.7) and construct . Then, the Drazin inverse matrix is given by

Proof. Write the matrix in a block form where and are tridiagonal matrices and is given by
Matrix can be obtained from the Schur complement of matrix , . Then, we obtain From the above expression and by definition of , Using the expression to obtain the Drazin inverse given in [11] we obtain Using the expression given in (3.5) to obtain and by simple mathematical calculations expression (3.6) is obtained.
Now, the next step is to check that this matrix verifies the properties of the Drazin inverse matrix. First, we need to know the index of . From the symmetry of the inverse matrix obtained in (3.5) we obtain , then .
By technical calculations we check that this matrix verifies the Drazin properties (i) , (ii) , and (iii) , to prove the proposition.

It is important to focus our attention on the structure obtained in (3.6) since it is useful to analyze the input-output behaviour in terms of the vector parameter . This fact allows to study directly the identification problem.

Proposition 3.2. Consider the system . Then, the first Markov parameter is given by

Proof. First we prove for all that Using the expression given in (3.5) we have that and since we have Expression (3.12) is easily seen using the above relation and the Drazin inverse of matrix given in (3.6).

Using Proposition 3.2 we give the structure of the rest of Markov parameters.

Proposition 3.3. Consider the system . Then, the Markov parameters , , are given by

Next, we solve the identifiability problem in the next result.

Proposition 3.4. The system is identifiable.

Proof. We consider two structured systems and with , , and such that they have the same input-output behavior (io), that is, , .
From expression (3.12), since , and taking , we have that On the other hand, from expression (3.16), since , Hence, and the system is identifiable.

4. Algorithm and Example

In this section we introduce an algorithm to obtain the parameters of a PDAE given by (2.1), using the Markov parameters technique.

Step 1. Introduce the size of the state vector, . Introduce the matrices that determine the known input-output behavior of the process.

Step 2. For , if , then rename .

Step 3. Construct .

Step 4. Construct .

Step 5. Introduce and .

Step 6. For , define .

Step 7. For , define .

For , apply the following.

Step 8. Prove that verifies If the result is false, then stop the process. This means that the given input-output behavior is not suitable to a system of type (2.1).

Step 9. Prove that , , verifies If the result is false, then stop the process. This means that the given input-output behavior is not suitable to a system of type (2.1).

Step 10. Obtain and .

Step 11. Obtain the matrices of the system (2.6).

Step 12. Obtain the matrices of the system (2.1).

This algorithm has been implemented in MATLAB. To illustrate this algorithm, we present the following academic example.

Example 4. Consider the PDAE given in (2.1), with . The discrete system associated with this PDAE is given by (2.6) with . The input-output behavior of a process is given by
Then we apply the above algorithm. First, we introduce the size of the state vector, . Introduce the matrices that determine the known input-output behavior of the process, that is, Applying MATLAB algorithm we obtain , , and . The comparison is ok; ; . The discrete model (2.6) is given by

5. Conclusions

In this paper we show some results to identify and determine parameters in a mathematical model with a particular structure. To obtain these results we study a discrete model associated with a partial differential algebraic equation. This model has been obtained using finite difference method from the partial differential equation. Some results about the Markov parameters of the discrete system are given, and the identifiability problem is solved in the discrete case. By a similar way, the given results can be applied to different models with a similar structure to different areas.

In this case, the results are used to know if a partial differential algebraic equation with unknown parameters is identifiable and to obtain the parameter values. Finally, an algorithm, implemented in MATLAB, is given to obtain the parameters.

Acknowledgments

The authors are very grateful to the referees for their comments and suggestions. The paper is supported by Grant PAID-05-10-003-295 and Grant MTM2010-18228.