Improved Discrete Techniques of Time-Delay and Order Estimation for Large-Scale Interconnected Nonlinear Systems

The selection of a suitable model structure is an essential step in system modeling. Model structure is defined by determining the class, the time delay, and the model order. This paper proposes improved structural estimation procedures for large-scale interconnected nonlinear systems which are composed of a set of interconnected Single-Input Single-Output (SISO) Hammerstein structures and described by discrete-time stochastic models with unknown time-varying parameters. An extensive Determinant Report (DR) algorithm is developed to determine the order of the process. An improved discrete-time technique based onRecursive Extended Least Squares with Varying Time (RELSVT) delay method is proposed to estimate the time delays of the considered system. The developed theoretical analysis and simulation results prove the validity and performance of the proposed algorithms.

Note that the study of the problems related to the description of a dynamic system by a mathematical model has been always a main objective to be met in different fields of science and automatic engineering, more particularly in the identification process.In this case, several automation engineers are interested in the study of this problem and various works dealing with this topic are developed and published in the literature.It is important to note that, before applying the identification methods, it is necessary to have a priori information about the structure of the process to be estimated.However, the order and the time delay of the process may be unknown in a practical estimation problem.Therefore, the estimation of these structural parameters is a fundamental problem.Add to that, in automatic control and signal processing applications, the time-delay and the order estimation of a certain process can be a means to achieve a good model for the control design.In this instance, various structural estimation methods have been proposed in the literature [15][16][17][18][19][20][21][22][23][24][25][26][27][28].
The purpose of this framework is to develop two procedures which can be applied to the class of large-scale interconnected systems with unknown time-varying parameters, in order to formulate the structure estimation problem.We focused on the dynamic large-scale stochastic systems comprised of several SISO interconnected nonlinear subsystems and represented by an interconnection of  Hammerstein structures with unknown parameters.
The rest of this article was organized as follows.The second section describes the dynamics of large-scale interconnected systems by stochastic Hammerstein mathematical models with unknown time-varying parameters.Section 3 2 Mathematical Problems in Engineering details the formulation of the structural identification problem of interconnected nonlinear systems where two iterative procedures are proposed.The main purpose of these developed algorithms is to estimate the order model and the time delay for this class of dynamical systems.Numerical simulation examples are treated to validate our developed analytical results.Conclusions are delineated in Section 4.

Preliminary Study
A dynamic large-scale interconnected system is characterized not only by a large number of dynamics, but also by input and output signals.It can be presented by a mathematical model composed of a large number of equations, which are connected between themselves and contain several signals and parameters.In the literature, several presentations are used to describe the dynamics of this class of systems.We can cite input-output discrete models [11], continuous statespace models [1,2], and block-oriented mathematical models [10,[12][13][14].In this framework, we consider the class of largescale systems which consist of several interconnected Hammerstein subsystems operating in a stochastic environment with unknown time-varying parameters and described by the following dynamic equation: in which the nonlinear functions ℎ where  −1 is the backward shift operator;   () and   () represent, respectively, the output and the input of each interconnected nonlinear system   ;   ,   , and   correspond to the time delays of the system   and the other interaction systems   , ,  = 1, . . ., ,  ̸ = ;  , 1 ,  , 2 , and  , 3 are unknown constant parameters of the nonlinear functions; and   ( −1 , ),   ( −1 , ),   ( −1 , ),   ( −1 , ), and   ( −1 ) are described by with ,  = 1, . . ., ,  ̸ = , and    ,    ,    ,    , and    are the orders of the polynomials   ( −1 , ),   ( −1 , ),   ( −1 , ),   ( −1 , ), and   ( −1 ), respectively.For the sake of simplicity, in what follows, we retain that   ( −1 , ),   ( −1 , ),   ( −1 , ),   ( −1 , ), and   ( −1 ) have an even order   .We assume also that the noise {  ()} is a sequence of independent random variables with zero mean and constant variance  2  .In addition, this sequence is uncorrelated with all input and output signals.Furthermore, we can rewrite the system output in the following developed form: We note that some implementation difficulties are presented in the formulation of the parametric estimation problem in spite of the presence of some redundant parameters in model (4).From [12], we propose that  ,1 ,  ,1 (), and  ,1 () are constant and known, in such a way that  ,1 =  ,1 () =  ,1 () = 1, ∀, in order to avoid the posed problem.In this case, each output   () of the considered interconnected system becomes   (5) The formulation of the parametric estimation problem for this class of large-scale systems, which is described by the discrete model (5), was proved by Elloumi and Kamoun [12].In this instance, the developed parameter estimation procedure is based on the adjustable model and the prediction error method, starting from the knowledge of several measured input-output signals resulting from the considered process.In this case, we supposed that the time delays and the order model of the considered system are known.For this reason, serious errors can be obtained in the synthesis of a control design if the used order model and/or time-delay parameters are wrong.
Note that the dynamic process models are always desirable in order to provide the required design of regulators.In addition, the choice of the model structure (order and time delay) and the parameter estimation are two basic elements in the identification problem.Therefore, lots of works that used various methods to estimate the time delay and the order model have been published in the literature [15][16][17][18][19][20][21][22][23][24][25][26][27][28].

Structural Estimation
In this section, we are proposed to formulate the structural estimation problem for the class of large-scale interconnected nonlinear stochastic systems described by the mathematical model (1).The formulation of this problem can be classified into two subsections.
The first subsection seeks to determine the order of the considered process, where its dynamics are modeled by INARMAX mathematical model (5).The order estimate will be performed on the basis of an extension of DR algorithm that can be applied to this class of complex systems, whereas the formulation of the time-delay estimation problem is considered in the second subsection.In this context, an improved procedure based on RELSVT algorithm is proposed.

Order Estimation.
The adequate value of the order process is one of the major preoccupations in the identification problem.Lots of testing methods have been proposed for determining the orders of discrete dynamic processes, like the Determinant Report (DR) method, the behavior of error function, polynomial approach, test of normality, and so forth [15][16][17].In [15], the developed DR technique is based on the product-moment matrix containing several input-output signals measurements from the system.
In this part, the purpose is to extend this approach in order to estimate the order of large-scale interconnected nonlinear systems, which is described by the discrete-time mathematical model (1).
The organizational chart of the order estimate, represented in Figure 1, seeks to verify certain criteria after describing the system by a mathematical model with an order selected a priori   (  = 1, 2, . . .,  max , where  max is the maximum order).The exact model order of the system corresponds to the value   for which the calculated criterion presents a particular deviation.This procedure seeks to limit the number of possible orders and determine the accurate one, and this is based on the computation of the product-moment matrix with dimension (( 1 +  2 +  3 + 1)  , ( 1 +  2 +  3 + 1)  ), which is defined by with where   represents the measurements number and   depicts the model order of the process, ,  = 1, . . ., ,  ̸ = .The matrix   (  ) possesses the following properties: This matrix is positive definite and it is singular, if   >   .For   <   , DR  is expressed by with Report ( 9) is computed for different   values (  = 1, 2, . . .,  max ).When the Determinant Report value DR  (  ) presents an important increase over the previous value DR  (  − 1), the value   corresponds to the specific order of the considered system.As a result, the DR  procedure can be described as In order to clarify the developed approach and to test the performance of algorithm (11), we treat the following numerical example, which corresponds to a large-scale system constituted of two interconnected Hammerstein subsystems.Each system output is described by In this simulation example, the input   (),  = 1, 2, that applied to the interconnected Hammerstein system   is a high level pseudorandom binary sequence [−1.5, +1.5].
We remark that the evolution curves of the different determinant reports, as delineated in Figures 2 and 3, admit a specific jumping in   = 2.As a result, we can conclude that the model order is equal to 2 for each interconnected subsystem   .Therefore, these results indicate the performance and the efficiency of the developed algorithm (11) for the order estimate of a large-scale interconnected nonlinear system described by the Hammerstein model (1).
It should be mentioned that this approach is used to limit the number of possible model orders before starting the parameters estimation.This method gives the correct model order independently of the number of measurements because from all the obtained values of DR  (  ) only DR  (  ) of the exact order model presents a considerable step.For a system with small noise level, this test still indicates the correct order.So far, this method seems sufficient to be applied for a process with uncorrelated and high disturbances.Nevertheless, in real applications, when the noise is corrupted with all input and/or output signals, this technique becomes not suitable to determine the correct order model.In this sense, other tests methods can be used, which are developed in the literature to estimate the model order for simple systems, like the condition number method, test of signals errors, statistical -test, polynomial approach, test for normality, and so forth.These techniques could be applied in combination with the parameter estimation methods.Let us note that the determinant ratio approach proves to be very robust under all working conditions because, with this test, it is possible to reduce the possibilities of the model order, as well as save a lot of computing time.
Note that the practical value of these different testing order approaches depends on the purpose and further application of the estimated model.However, it must be considered that the use of one testing approach alone may produce a wrong order model; therefore, the use of various methods together can elicit the correct order model with important accuracy.

Time-Delay Estimation.
In reality, there are no physical systems without delay.The time delay as an active research in automatic control and signal processing is widely used in several industrial applications like chemical processes, energy processes, communication processes, and so on.However, it is assumed to be negligible in several researches, in order to simplify the study.To overcome this assumption, a variety of algorithms are introduced into the time-delay estimation, in order to improve the precision and the convergence in the modeling, identification, and control processes [23][24][25][26][27][28].The time-delay identification is a greatly studied problem with several works in the literature.In this instance, different discrete-time and continuous-time techniques for time-delay estimation were proposed in control process and signal processing.Some typical approaches are described: Zhang and Li proposed a time-delay identification approach based on cross-correlation technology in signal processes [25], Zheng and Feng developed a time-delay estimation algorithm

Criterion computation
Interconnected nonlinear system described by a discrete-time mathematical model with time delays d i , based on correlation analysis [26], Gao et al. presented an iterative approach for time-delay estimation based on the output error between the process output and the predictive model [27,28], and so forth.
In what follows, we propose an iterative approach incorporated with a Recursive Extended Least Squares estimator and based on the prediction error method, in order to estimate simultaneously the time delays and the parameters of a large-scale interconnected nonlinear process, which is modeled by (1).
The general procedure of time-delay estimation, as shown in Figure 4, is based on the computation of the criteria  ,  ,  ,  , and  ,  for each of the time delays   ,   , and   .The minimum value of each criterion corresponds to the exact estimated delay.
To apply this procedure, let us consider a large-scale system decomposed into N interconnected nonlinear subsystems and described by the mathematical model (1).Each system output can be written in the following manner: where with in which the delays  max ,  max , and  max correspond to the superior ones which are selected a priori from the real delays (  = 1, . . .,   ,  = 1, 2, 3).Thus, the output   () can be expressed in the following developed form: Equivalently, its matrix form is as follows: where the parameters vector  *   () and the information vector  *  () are defined as The previous vectors, given by ( 19) and (20), are used in a recursive algorithm of parameter estimation, named Recursive Extended Least Squares with Varying Time (RELSVT) delay estimator.As a result, fewer steps are required to estimate the time delays.

(i) Determine the maximum values of 𝑏
(ii) Verify the following relationships: Thereby, the real delay satisfies the relation given as where Ĝ,  (), Ĥ,  (), and K,  () are given as follows.
The minimum values of  ,  ,  ,  , and  ,  indicate, respectively, the exact time delays   ,   , and   of the process.
To validate the proposed approach of time-delay estimation, let us consider the previous numerical example with  1 = 1 and  2 = 2. Figures 5 and 6 represent the evolution curves of  ,  ,  = 1, 2, for each interconnected nonlinear subsystem.
According to the evolution curves of  1, 1 and  2, 2 , it can be remarked that d1 = 1 and d2 = 2.These simulation results prove the validity and feasibility of the proposed algorithm for time-delay estimation of the large-scale interconnected nonlinear processes.
It can be remarked that the estimated delays depend on the estimated parameters of the model.Therefore, (31) is added to an Extended Least Squares algorithm to estimate simultaneously the parameters and the time delays of the model.Note that computer simulations show that discrete solutions properties of the delays equations make them robust to estimation errors in the estimate of parameters vector (19) and the presence of Gaussian-distributed noise.If the parameters are correctly estimated and no disturbances exist, (31) will have a zero minimum.If there is a small estimation error and small noise level, (31) is still enough to identify the delays correctly.As the estimation error and the level noise increase, the delays estimation can be conducted using the maximum likelihood algorithm in the estimation of the process parameters and based on (31) (the maximum likelihood estimator is more robust than the Recursive Extended Least Squares algorithm in the convergence sense).
It should be mentioned that (31) can be added easily to any recursive estimator.The implementation of this equation requires minimum storage and computation time as they contain multiplications and additions.However, in realistic applications, the corrupting noise is correlated with input and/or output signals, so this makes the estimates of the parameters asymptotically biased, as well as a wrong estimated time delay.In this case, we can use other techniques to estimate the time delays and the parameters of the model, such as the correlation analysis, the crosscorrelation method, and generalized  norm estimation.
In conclusion, a priori information of noise distribution is necessary in order to develop an efficient estimation.
In the identification systems, two main problems must be considered.The first one is the identification of the time delay, and the second one is the identification of the parameters.In this case, the time-delay estimation problem is a great research axis with several references dealing with several approaches for time-delay estimation.These approaches can be classified into four classes.The first class denotes the time-delay approximation approach (time domain, frequency domain, and Laguerre domain).In this case, the time delay is Mathematical Problems in Engineering estimated from an approximation of the model.The second class represents the time-delay parameter method (one step, two steps, and sampling method).In this case, the time delay and the parameters of the model are identified alternatively, simultaneously, or using the sampling system.The third class is the area and moments approaches.In this event, we estimate firstly the step or the impulse response.After that, we identify the time delay.The last class represents the high-order statistics approach.The main disadvantage of these techniques is that parameters estimation is required on each iteration step of time-delay estimation, and this makes the implementation difficult.Furthermore, it can be indicated that there is not an efficient solution to the timedelay problem.The general agreement is on which method is the best.Thereby, a comparative study of these techniques for time-delay estimation of large-scale interconnected nonlinear systems will be investigated in the future research.

Conclusions
Two extended discrete-time structural estimation methods, which are applied for estimating the time delay and the order of a dynamical process, are developed.We have particularly focused on the class of large-scale time-varying stochastic systems constituted by various interconnected Hammerstein structures.These systems are described by discrete-time stochastic Hammerstein mathematical models with unknown time-varying parameters.
The main research work developed in this article is divided into two important parts.The first part formulated the problem of the order estimation for the class of large-scale interconnected nonlinear systems, using the determinant ratio test DR  .The second part proposed an extended iterative procedure, which permits estimating the time delays of the considered system.
These developed approaches are validated by treating two numerical simulation examples of a large-scale system composed of two interconnected Hammerstein structures.The obtained results prove the performance and the feasibility of the proposed methods in spite of the presence of interactions signals, Gaussian distribution noise, and parameters variations.
a model of order m i m i = m i + 1

Figure 4 :
Figure 4: General procedure of the delay estimation.