Identifiability for a Class of Discretized Linear Partial Differential Algebraic Equations

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.


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 1-4 .

Mathematical Problems in Engineering
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 f t, x ∈ R n is a smooth function, E / 0, F p / 0, G p are matrices in R n×n , all of which can be singular, and p α, β ∈ R 2 is a unknown parameter vector.The solution of 1.1 , denoted by z z t, x , depends on the parameter vector p, the spatial variable x ∈ Ω and the time value t ∈ 0, T .
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.

Problem Description and Discretization
In this section we consider the second-order partial differential algebraic equations PDAEs given by where and G p βI E, β > 4α.The differential algebraic aspect of the system is due to singularity of E.
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 Δx, where Δx is the distance between two neighboured nodes of the grid, and we discretize the interval 0, T by a one-dimensional grid with stepsize Δt.
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 Δx Δt 1.The discretized problem can by written as where p 2α − β and q −α.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

2.5
To obtain a structure which simplifies the solution we write the full system of equations as follows: where the state vector x k z 1 k T z 2 k T T with z 1 k z 1 k,j j 0,...,N and z 2 k z 2 k,N−j j 0,...,N , the control vector u k u k ∈ R, and the coefficient matrices are given by with l N 1 n 1 .These matrices have a fixed structure with the parameter vector p p, q belonging to a subset P ⊆ R 2 , 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 λ ∈ C such that det λE − A p / 0 regularity condition .The matrix A p is a symmetric, tridiagonal, and regular matrix.This follows from the eigenvalues of A p given by and using the expressions of p and q, since it is easy to check that This establishes the regularity of A p and this fact guaranties the regularity of the system 2.6 .

The Identification Problem
The system 2.6 is a singular discrete-time system which, from the regularity of A p , we can rewrite as follows: where E p A −1 p E and B p A −1 p B. This system is denoted by S p .Since there exists a scalar λ ∈ C such that det λE − I / 0, the explicit solution of the system is given by where E D p denotes the Drazin inverse of matrix E p , ν is the index of E p , and x 0 is an admissible initial condition, x 0 ∈ X 0 , with where The input-output behavior of the system S p is obtained from the Markov parameters associated which the system.These parameters are defined as 3.4

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 inputoutput behavior for every parameter set 3 .That is, the system S p is globally structurally identifiable if and only if, for almost any two candidates parameter vector values p, q ∈ P, io p io q imply p q, where io • denotes the input-output behavior of the system S p .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 S p to the identifiability analysis.For that, we need to know the structure of E D p .
By definition of E p to obtain E D p we use the inverse of the tridiagonal matrix.A concise expression of A −1 p a ij i,j is obtained as follows: where α |j − i| and d i pd i−1 − q 2 d i−2 , for i 2, . . ., n, with d 0 1 and d 1 p.Note that d i M i p , where M i p is the principal minor of order i of the matrix A p and the symmetry property of matrix A p is transferred to A −1 p .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 E p is given in the following result.

6 Mathematical Problems in Engineering
Proof.Write the matrix A p in a block form

3.7
where A 1 and A 4 are tridiagonal matrices and A 2 is given by

3.8
Matrix A −1 p can be obtained from the Schur complement of matrix Then, we obtain From the above expression and by definition of

3.10
Using the expression to obtain the Drazin inverse given in 11 we obtain

3.11
Using the expression given in 3.5 to obtain A −1 4 and by simple mathematical calculations expression 3.6 is obtained.Now, the next step is to check that this matrix E D p verifies the properties of the Drazin inverse matrix.First, we need to know the index of E p .From the symmetry of the inverse matrix obtained in 3.5 we obtain rank E p rank E p 2 , then ind E p 1.By technical calculations we check that this matrix verifies the Drazin properties i 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 p.This fact allows to study directly the identification problem.

3.12
Proof.First we prove for all i 1, . . ., n that q a i−1,n a i 1,n pa i,n 0.

3.13
Using the expression given in 3.5 we have that Expression 3.12 is easily seen using the above relation and the Drazin inverse of matrix E D p given in 3.6 .
Using Proposition 3.2 we give the structure of the rest of Markov parameters.
Proposition 3.3.Consider the system S p .Then, the Markov parameters V j, p v j, p i i , j > 1, are given by

3.16
Next, we solve the identifiability problem in the next result.
Proposition 3.4.The system S p is identifiable.
Proof.We consider two structured systems S p 1 and S p 2 with p 1 , p 2 ∈ P, p 1 p 1 , q 1 , and p 2 p 2 , q 2 such that they have the same input-output behavior io , that is, V j, p 1 V j, p 2 , j ≥ 1.
From expression 3.12 , since v 1, p n − q/p v 1, p n−1 , and taking i n, we have that 3.17 On the other hand, from expression 3.16 , since v 2,

3.18
Hence, p 1 p 2 and the system is identifiable.

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, n.Introduce the matrices {V j v j i i 1,...,n , j 1, . . ., n} that determine the known input-output behavior of the process.
Step 8. Prove that V 1 verifies

4.1
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 .

4.2
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 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.Consider the PDAE given in 2.1 , with n 1 n 2 1.The discrete system associated with this PDAE is given by 2.6 with A p ∈ R 6 .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, n 6. Introduce the matrices {V j v j i i 1,...,n , j 1, . . ., n} that determine the known input-output behavior of the process, that is, v :, 1 0; 0; 0; −

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.

Proposition 3 . 1 .
Consider the matrices given in 2.7 and construct E p A −1 p E. Then, the Drazin inverse matrix E D p e ij i,j is given by