Mathematical Problems in Engineering

Volume 2016 (2016), Article ID 6435320, 9 pages

http://dx.doi.org/10.1155/2016/6435320

## Identification of Coupled Map Lattice Based on Compressed Sensing

^{1}Information Security Center, State Key Laboratory of Networking and Switching Technology, Beijing University of Posts and Telecommunications, Beijing 100876, China^{2}National Engineering Laboratory for Disaster Backup and Recovery, Beijing University of Posts and Telecommunications, Beijing 100876, China

Received 25 January 2016; Accepted 29 February 2016

Academic Editor: Yang Tang

Copyright © 2016 Dong Xie 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

A novel approach for the parameter identification of coupled map lattice (CML) based on compressed sensing is presented in this paper. We establish a meaningful connection between these two seemingly unrelated study topics and identify the weighted parameters using the relevant recovery algorithms in compressed sensing. Specifically, we first transform the parameter identification problem of CML into the sparse recovery problem of underdetermined linear system. In fact, compressed sensing provides a feasible method to solve underdetermined linear system if the sensing matrix satisfies some suitable conditions, such as restricted isometry property (RIP) and mutual coherence. Then we give a low bound on the mutual coherence of the coefficient matrix generated by the observed values of CML and also prove that it satisfies the RIP from a theoretical point of view. If the weighted vector of each element is sparse in the CML system, our proposed approach can recover all the weighted parameters using only about samplings, which is far less than the number of the lattice elements . Another important and significant advantage is that if the observed data are contaminated with some types of noises, our approach is still effective. In the simulations, we mainly show the effects of coupling parameter and noise on the recovery rate.

#### 1. Introduction

Coupled map lattice (CML), developed by Kaneko [1–3], is a dynamical system that can be described as an array of smaller finite-dimensional subsystems endowed with local interactions. The evolution of the subsystem at one lattice element depends not only on the states of the element itself, but also on the state of the other elements. If the coupling interaction is a function of all the elements, then the CML is called globally coupled map lattice. So far, a very large number of remarkable results about this topic were obtained by researchers working in different areas of physics, biology, mathematics, and engineering. However, most of the papers mainly concentrated on analysing the evolutionary behaviours and properties of CML [4–6]. And the topic about the parameter identification or reconstruction of CML from measured data has received relatively little attention. Over the past fifteen years, several methods for the identification problem have been presented [7–11]. Billings and Coca [8] introduced a novel approach to the identification of CML and the method exploited the regularity of the CML model. Subsequently, Billings et al. [10] pointed out some drawbacks of the method [7, 9] and presented a multiresolution approximation for the underlying spatiotemporal dynamics using B-spline wavelet and scaling functions. Ning et al. [11] first transformed the nonlinear spatiotemporal system into a class of multi-input-multi-output (MIMO) partially linear systems and then proposed an effective online identification algorithm for this nonlinear spatiotemporal system. Subsequently, Ning et al. [12] proposed a novel kernel-based learning algorithm and it can be used to effectively achieve accurate derivative estimation for nonlinear functions in the time domain. Kohar et al. [13] studied the usefulness of coupled redundancy as a mechanism for reduction in local noise in coupled map lattices and investigated the role of network topology, coupling strength, and iteration number in this mechanism. But sometimes in order to improve the efficiency of parameter identification, we need to reduce the number of observations or samplings under the condition of accurate identification. In addition, from an engineering perspective, we cannot ensure that the received data are truly generated by the CML system because our sampled data may be contaminated with noise when we perform some indispensable tasks (e.g., data sampling and data transmission). To the best of our knowledge, most of papers about this topic (including the papers mentioned above) did not consider the above two fascinating and significant problems. And the parameter identification of CML using a very few number of observations or samplings is still regarded as a very intractable problem. Another challenging problem is that if the observed data are contaminated with noise, how can we accurately identify all the weighted parameters?

The compressed sensing theory, introduced in [14–17], breaks the limitation of Nyquist-Shannon sampling theorem and is widely used in various fields, such as information theory, medical imaging, earth sciences, and astronomy. The standard model of compressed sensing is an underdetermined linear system , where , , , and . It is clear that, given a coefficient matrix and an observation vector , there exist infinite such that . In order to recover the unique vector from the observed vector , we need two conditions; that is, the vector is sparse and the matrix must satisfy some special properties, such as restricted isometry property (RIP) and mutual coherence [18]. In other words, if the matrix satisfies the RIP, then we can recover the -dimensional sparse vector from -dimensional observed vector . Further, the RIP is sufficient for a variety of algorithms to be able to successfully recover a sparse vector from noisy samplings. These two useful theoretical results about compressed sensing provide insights to solve the above stubborn problems about the identification of CML.

Different from the existing methods about parameter identification of CML system, we independently present a new approach using the compressed sensing theory in this paper. Specifically, we first transform the parameter identification problem into the reconstruction problem of underdetermined linear system through simple mathematical transformation. Subsequently, we give a low bound on the mutual coherence of sensing matrix using the inverse formula of Cauchy-Schwarz inequality and then prove that it satisfies the RIP with suitable parameters. At last we exploit the classic Orthogonal Matching Pursuit algorithm in compressed sensing to recover the weighted parameters of CML system. The result shows that if the weighted vector of every lattice element has at most nonzero entries, then we just need observations, which is far less than , to accurately identify the weighted vector. Another critical advantage of our approach is that the reconstruction algorithm can still work if the sampled data are contaminated with some types of noises (e.g., Gaussian noise, uniform noise). We investigate the effects of coupling parameter and noise on the recovery rate in our simulations. We found an interesting fact that when the sampled data are contaminated with uniform noise, the recovery effects are better than that if the observed data are authentic.

The remainder of this paper is organized in the following manner. We recall the standard CML model in Section 2. In Section 3 we introduce some backgrounds about the compressed sensing theory. We present the theoretical results about our proposed identification method based on compressed sensing in Section 4 and show our simulations in Section 5. Finally, we draw our conclusions in Section 6.

#### 2. The CML Model

We give a brief review of standard CML model in this section. Generally speaking, the widely used CMLs can be divided into two categories, that is, the diffusive (DCML) model and the global (GCML) model. The difference between these two models is whether the evolution of the subsystem at any lattice element depends on the states of all the other elements. In a standard GCML model, if we set some weighted parameters to 0, then it can be seen as a DCML model. Without loss of generality, we consider the following GCML model: where defines the state of the lattice element at a discrete time step , is the coupling parameter, and are maps with respect to the local dynamics and nonlocal system, and is the unknown weighted vector of the lattice element . The two widely used forms of , namely, and , are corresponding to synchronous and asynchronous updates of the lattice, respectively. For the sake of simplicity, we just consider the synchronous case in our paper; that is, . In many scenarios, is the standard logistic map, which can be written as where and the range eventually lies in the interval (0,1).

#### 3. Compressed Sensing

In the field of compressed sensing, the researchers mainly concentrate on designing the excellent sensing matrices and the efficient reconstruction algorithms. Compressed sensing provides a method to obtain a unique solution from an underdetermined linear system taking advantage of the prior knowledge that the true solution is sparse. In mathematical terms, we consider the following underdetermined linear system: where and . This equation depicts the process that the original -dimensional data is compressed into an -dimensional data . In this formula, denotes the sensing matrix and we call observed vector or observed values. A vector is said to be -sparse when it has at most nonzero entries. For a -sparse original vector , if the sensing matrix satisfies some suitable condition, then there exist some algorithms to recover the vector from the observed vector . Next, we recall some fundamental results about sensing matrices and reconstruction algorithms in the following two subsections.

##### 3.1. Sensing Matrix

In this subsection, we address two significant concepts frequently used in compressed sensing theory, that is, RIP and mutual coherence. The famous RIP was first presented by [14], which is defined as follows.

*Definition 1. *A matrix satisfies the RIP of order if there exists a such that holds for all -sparse vectors .

In fact, the RIP sufficiently guarantees that the original vector can be exactly recovered from the observed values. Furthermore, if the matrix satisfies the RIP, then it is also sufficient for a variety of algorithms to be able to recover the spare vector from observed data contaminated with noise. The following lemma [19] denotes noiseless recovery if the matrix satisfies the RIP.

Lemma 2. *Let be a sensing matrix in (3) that satisfies the RIP of order with . Then the original sparse vector can be exactly recovered.*

Baraniuk et al. [20] have proved that the random sensing matrix, whose entries are identically and independently sampled from a Gaussian distribution or a sub-Gaussian distribution, satisfies the RIP with overwhelming probability. Similar to a Gaussian random matrix, a Bernoulli matrix or a Random Discrete Fourier Transform (RDFT) matrix can also be used as sensing matrix. However, given a specific matrix , it is very hard to verify whether satisfies the RIP or not. The concept of mutual coherence provides a more simple choice [18, 21].

*Definition 3. *Given a matrix , denotes the th column of . The mutual coherence of the matrix , , is the largest absolute normalized inner product between any two columns of . Namely,

From Cauchy-Schwarz inequality , it is clear that the mutual coherence of any matrix is bounded by 1. In order to complete our proof, we will review the Gershgorin circle theorem [22] here.

Lemma 4. *The eigenvalues of an matrix with entries , , , lie in the union of discs centered at and with radius .*

##### 3.2. Recovery Algorithm

A straightforward approach to obtain the original -sparse vector from the above underdetermined system (3) can be viewed as the optimization problem of which is called optimization problem. This problem is very difficult to solve and it is NP-hard in general case [14]. Since the convex property of norm, a classic method used in compressed sensing is to replace with ; that is,

In fact, Candès [19] has pointed out that the solution to (6) is that of (7) when satisfies the RIP with . Although some convex optimization techniques are powerful for sparse vector reconstruction, there also exist a variety of greedy algorithms [23–25] for solving this problem. We describe the classic Orthogonal Matching Pursuit (OMP) algorithm as follows because our simulations need to use it.

In Algorithm 1, denotes the support set and is the vector obtained by collecting the entries of corresponding to the set . is the submatrix of composed of the column vectors of corresponding to the support set . During each pass through the loop, the algorithm first chooses the optimal index and updates the support set and then computes the vector which minimizes such that .