Abstract
Many industrial processes and physical systems are spatially distributed systems. Recently, a novel 3-D FLC was developed for such systems. The previous study on the 3-D FLC was concentrated on an expert knowledge-based approach. However, in most of situations, we may lack the expert knowledge, while input-output data sets hidden with effective control laws are usually available. Under such circumstance, a data-driven approach could be a very effective way to design the 3-D FLC. In this study, we aim at developing a new 3-D FLC design methodology based on clustering and support vector machine (SVM) regression. The design consists of three parts: initial rule generation, rule-base simplification, and parameter learning. Firstly, the initial rules are extracted by a nearest neighborhood clustering algorithm with Frobenius norm as a distance. Secondly, the initial rule-base is simplified by merging similar 3-D fuzzy sets and similar 3-D fuzzy rules based on similarity measure technique. Thirdly, the consequent parameters are learned by a linear SVM regression algorithm. Additionally, the universal approximation capability of the proposed 3-D fuzzy system is discussed. Finally, the control of a catalytic packed-bed reactor is taken as an application to demonstrate the effectiveness of the proposed 3-D FLC design.
1. Introduction
Many industrial processes and physical systems such as industrial chemical reactor [1, 2], semiconductor manufacturing [3], and thermal processing [4] are “distributed” in space. They are usually called spatially distributed systems, or distributed parameter systems [1]. The states, controls, and outputs of such systems depend on the space position as well as on the time [2]. Traditionally, model-based methods are used to control such systems, where a good mathematical model is definitely required. However, the process model may not be easily obtained in many complex situations, and then, a model-free control method has to be used. This leads to the recent development of the novel three-dimensional fuzzy-logic control (3-D FLC) [5–8], which has the inherent capability to process spatiotemporal dynamic systems. The 3-D FLC uses one kind of three-dimensional (3-D) fuzzy set (shown in Figure 1), which is composed of the traditional fuzzy set and a third dimension for the spatial information, and executes a 3-D rule inference engine. It is actually a kind of spatiotemporal fuzzy-control system with the traditional model-free advantage.
To date, the 3-D FLC design has been focused on an expert-knowledge-based approach [5], that is, the fuzzy-rule design is from human experts’ knowledge. In this approach, human knowledge to the control solution must exist, and be structured. Practically, experts may have problems structuring the knowledge [9]. Sometimes, although experts have the structured knowledge, they may sway between extreme cases: offering too much knowledge in the field of expertise, or tending to hide their knowledge [9]. Thus, we often lack expert knowledge for control that is usually hidden in an input-output data set. Under this circumstance, a data-driven design becomes a good choice for the 3-D FLC, that is, extraction of fuzzy rules from a spatiotemporal input-output data set. Since the research on the 3-D FLC is just at the beginning stage, extracting 3-D fuzzy control rules from a spatiotemporal data set is still a challenging and open problem for spatially distributed systems.
Traditional data-driven FLC design methods have been developed in the past three decades. They are usually composed of three parts: rule generation, structure optimization, and parameter optimization [10]. For instance, grid partitioning of multidimensional space [11] and clustering technique [12] can be used to generate rules automatically; reducing redundancy variable [12], fusing similar clusters [13], and fusing similar fuzzy set [14] can be applied to reduce the rule number and realize the structure optimization; genetic algorithm [15] and gradient decent approach [16] can be adopted for fine tuning of membership function and realize the parameter optimization. For a complete review of data-driven fuzzy system design, one can further refer to [10]. These methods provide useful solutions to a traditional FLC design.
In this study, we aim at developing a new data-driven 3-D FLC design method based on clustering and SVM-regression learning. The initial 3-D rule base is first generated by a nearest-neighborhood-clustering method from a spatiotemporal data set via defining Frobenius norm as a distance. Then, the initial 3-D rule base is simplified based on similarity measure technique defined for 3-D fuzzy sets and 3-D fuzzy rules. Subsequently, an SVM-regression learning algorithm is used to learn the parameters of the rule consequent parts. In addition, the universal approximation capability of the proposed 3-D fuzzy system is discussed.
The paper is organized as follows. Preliminaries about 3-D FLC and SVM regression are addressed in Section 2. In Section 3, a clustering and SVM-regression learning-based 3-D fuzzy control design methodology is presented in detail. In Section 4, the universal approximation capability of the proposed 3-D fuzzy system is presented. In Section 5, a catalytic packed-bed reactor is presented as an example to illustrate the proposed design scheme of a 3-D FLC and validate its effectiveness. Finally, conclusions are given in Section 6.
2. Preliminaries
2.1. 3-D FLC
The 3-D FLC is designed to have the inherent capability to deal with spatial information and its basic structure is shown in Figure 2. It has a similar functional structure similar to the traditional FLC, which consists of three basic blocks: fuzzification, rule inference, and defuzzification. However, it will differ in the detailed operations because of the spatial processing requirement. Generally, the 3-D FLC will be involved with the following basic designs: 3-D membership function (MF), 3-D fuzzification, 3-D rule base, 3-D rule inference, and defuzzification. One can refer to [5] for detailed description. Once each component of a 3-D FLC is set, a precise mathematical formula of the 3-D FLC can be derived.
Assumed that we have 3-D fuzzy rules represented by the following expression: where denotes the th spatial input variable (), is the input of from the sensing location (); denotes one-dimensional space in a discrete space domain ; denotes a 3-D fuzzy set, ; denotes output variable (the control action); denotes a traditional fuzzy set.
If Gaussian type 3-D membership functions (MF) are used to describe 3-D fuzzy sets in (2.1), then we have where denotes the Gaussian type 3-D MF of the spatial input in the lth rule; and are the center and width of , respectively; and denote center and width of the Gaussian type 2D MF of the spatial input at the sensing location . The Gaussian type 3-D MF can be regarded as an assembly of multiple Gaussian type 2-D MFs over the space domain . Then, the Gaussian type 2-D MF of the spatial input at the sensing location is given as where and .
Furthermore, if we employ singleton fuzzification, “product” -norm and “weighted aggregation” dimension reduction [6] in the 3-D rule inference, singleton fuzzy sets for the output variable, and “center of sets” defuzzification [17], the 3-D FLC can be mathematically expressed as where is a spatial input vector with as the input domain, as the number of sensors, and as the number of spatial inputs; is the number of rules; is the spatial weight from the spatial point [6]; is the nonzero value in the singleton fuzzy set of the output variable for the rule.
In (2.4), let then (2.4) can be rewritten as Similar to a traditional FLC [16], we define as a spatial fuzzy basis function (SFBF). Each SFBF corresponds to a 3-D fuzzy rule, and all the SFBFs correspond to a 3-D rule base. Mathematically, a 3-D FLC is a linear combination of all the SFBFs.
Equation (2.6) shows that the 3-D FLC is a nonlinear mapping from the input space to the output space . It provides us a way to understand and analyze the 3-D FLC from the point of view of function approximation. In Section 4, we will prove that the 3-D FLC has a universal approximation property based on the nonlinear mapping in (2.6).
2.2. Linear SVM Regression
An SVM is a learning algorithm that originated from theoretical foundations of the statistical learning theory [18] and has been widely used in many practical applications, such as bioinformatics, machine vision, text categorization, handwritten character recognition, time series analysis, and so on. The distinct advantage of the SVM over other machine learning algorithms is that it has a good generalization ability and can simultaneously minimize the empirical risk and the expected risk [19]. The SVM algorithms can be categorized into two categories: SVM classification and SVM regression. In this study, we are concerned with the SVM regression with -insensitive loss function [20].
Suppose we have a training set consisting of pairs , where the inputs are -dimensional vectors, and the labels are continuous values. In -SVM regression, the goal is to find a function so that for all training patterns has a maximum deviation from the target values and has a maximum margin. The -insensitive loss function is defined as follows: The -insensitive loss function defines an tube [9].
The regression problem can be formulated as a convex optimization problem as follows: subject to where and are slack variables, and the constant is a design parameter chosen by the user, which determines the trade off between the complexity of and the approximate error.
The above optimization problem can be solved in a dual space. By introducing the Lagrange multipliers, the primal optimization problem can be formulated in its dual form as follows: subject to Solving the dual quadratic programming problem, we can find an optimal weight vector and an optimal bias of the regression hypersurface given as follows: Then, the best regression hypersurface is given by The training pattern with nonzero is called support vector (SV).
3. Clustering and SVM-Regression Learning-Based 3-D FLC Design
Clustering and SVM-regression learning-based 3-D FLC design is a novel design of a 3-D FLC by integrating a nearest-neighborhood-clustering and an SVM-regression. The design methodology can be depicted by Figure 3. Firstly, a nearest-neighborhood-clustering method with Frobenius norm defined as a distance is employed to mine the underlying knowledge of the spatiotemporal data set S and yield the initial structure, that is, antecedent part of 3-D fuzzy rules. Because the obtained input space partition may have redundancy in terms of highly overlapping MFs, it is necessary to optimize the obtained initial fuzzy partition. Then, a similarity measure technique is utilized to merge similar 3-D fuzzy sets and to merge similar 3-D fuzzy rules, and then to simplify the initial rule structure. Finally, a linear SVM-regression algorithm is used to learn the parameters of the consequent parts based on an equivalence relationship between a linear SVM regression and a 3-D FLC.
The spatiotemporal data set S from a spatially distributed system is composed of n spatiotemporal input-output data pairs given as follows: where denotes the value of s spatial input variables at the kth sampling time, denotes the value of th spatial input variable at the kth sampling time , denotes the output value at the kth sampling time, n denotes the number of sampling time, and p denotes the number of sensors. Since infinite sensors are used, is a matrix with p rows and s columns.
3.1. Initial Structure Learning
3.1.1. Nearest Neighborhood Clustering Method
Clustering method is one of the data-driven learning tools for unlabeled data. It can mine underlying knowledge (or data structure) from a dataset that is difficult for humans to manually identify. One of the simplest clustering algorithms is the nearest-neighborhood clustering algorithm [16]. However, the existing nearest-neighborhood clustering algorithm has not the capability to deal with spatiotemporal data. In this study, we expand its capability to deal with spatiotemporal data set , which is of matrix form. The key point is that the Frobenius norm given in (3.2) is used for defining a distance in a nearest neighborhood clustering algorithm.
The nearest neighborhood clustering algorithm is summarized as follows.(i)Step 1: Begin from the first spatiotemporal data . Let the first cluster center be , the number of data pairs be 1, and the threshold be for generating new fuzzy rules.(ii)Step 2: Suppose that the kth spatiotemporal data is considered, when N clusters have been generated and their centers are , and respectively. Firstly, compute the distance between and each center of N clusters using . Then, compute the threshold using Hence, the corresponding cluster center is taken as the nearest neighborhood cluster of .(iii)Step 3: (a) If , then is taken as a new cluster center, and let , and . (b) If belongs to the cluster with the center . The center of cluster is tuned by introducing a learning rate as follows: and let .(iv)Step 4: Let . If , then quit. Otherwise, back to Step 2.
3.1.2. Rule Extraction and 3-D MF Construction
After clustering learning, we obtain an input space partition with cluster centers . Then, we will produce antecedent part of rule base and construct 3-D MFs in terms of the partition. Each cluster corresponds to a 3-D fuzzy rule. Assumed that we employ Gaussian type 3-D MF. Then, the cluster center corresponds to the center of Gaussian type 3-D MFs in the antecedent part. Thus, the number of fuzzy rules is equal to the number of clusters . In addition, we determine the width of the Gaussian MFs in terms of the domain of variables. For instance, the width of the Gaussian type 3-D MFs from the same sensing location are defined as where and are the maximum and the minimum bound values of the th spatial input variable, respectively.
3.2. Structure Simplication
After the initial structure learning, the obtained fuzzy partition of the input space and fuzzy rules may have redundancy in terms of highly overlapping MFs. In this step, we will simplify the fuzzy partition and fuzzy rules. The crucial technique for simpification is similarity measure. The previous similarity measure techniques [14, 21, 22] developed for traditional fuzzy sets and traditional fuzzy rules are not suitable to 3-D fuzzy sets and 3-D fuzzy rules. In this study, we will define a new similarity measure technique.
3.2.1. Similarity Measure
Firstly, we define the similarity of two 3-D fuzzy sets and as below. where is a distance between and . Since Gaussian type 3-D MFs are chosen, the following simple expression can be used to approximate the distance: where and are center and width of the Gaussian type 3-D MF () at sensing location , respectively.
Based on the similarity measure, we can merge similar 3-D fuzzy sets, or merge similar 3-D fuzzy rules.
(i) Merge of Two Similar 3-D Fuzzy Sets and
Firstly, the similarity between and is computed according to (3.6). If is higher than a threshold, we can conclude that and are similar, and then merge them into a new 3-D fuzzy set . The center and width of are viewed as the average values of and , and are given as the following:
(ii) Merge of Two Similar 3-D Fuzzy Rules and
The similarity and is inferred by measuring their similarity in the antecedent part. For instance, the similarity computation between and is given by
where and have the same rule form as in (2.1), () denotes the 3-D fuzzy set for the th spatial input variable in the th (th) rule. If is higher than a threshold, we can conclude that and are similar, and then merge them into a new 3-D fuzzy rule . The merging of two 3-D fuzzy rules is realized by merging the two fuzzy sets of each spatial input variable in the two 3-D fuzzy rules, respectively.
3.2.2. Similarity Measure-Based Structure Simplification
Based on the similarity measure, the simplification task includes removing 3-D fuzzy sets similar to the universal set, merging similar 3-D fuzzy sets, and merging similar rules. The detailed procedure of structure simplification is summarized as follows.(i)Step 1: Given a 3-D fuzzy rule base . Firstly, set proper thresholds: for removing 3-D fuzzy sets that are similar to the universal set, for merging similar 3-D fuzzy sets, and for merging 3-D fuzzy rules with similar antecedents.(ii)Step 2: Calculate with . Let and select and .(iii)Step 3: If , merge and into a new 3-D fuzzy set , set and , and back to step 2. If no more two 3-D fuzzy sets have the similarity with (), then go to step 4.(iv)Step 4: Remove the 3-D fuzzy set similar to the universal set and the rule with membership function that is always near zero over the space domain.(v)Step 5: Calculate the similarity of two rules with . Let .(vi)Step 6: If , merge the th and the th rules into a new rule and substitute them. Let , and back to step 5. If no more rules have similarity with , then quit.
Generally speaking, the threshold is higher than the threshold , while the choice of a suitable threshold depends on the application. The lower is set, the less fuzzy sets and less fuzzy rules are yielded in the resulting rule base. In this study, we set , and .
3.3. Parameter Learning
After the structure simplification, we obtain a rule base with optimized antecedent parts. For a complete rule base, the rest task is to determine the consequent part parameters. In this study, we employ an SVM regression algorithm to learn the consequent part parameter in the 3-D FLC.
Firstly, the original input samples are transformed into new samples. Utilizing the spatial fuzzy basis functions in (2.5), we can transform each spatial input sample in into a new input sample . Then, the original data set in (3.1) can be transformed into a new data set as follows:
Secondly, an equivalence relationship of an SVM regression and a 3-D FLC can be derived based on the new data set . From (2.13), the final decision function of an SVM can be described with the following form: where and are associated learning parameters in a SVM, The training pattern with nonzero is called support vector (SV). Furthermore, (3.11) can further be expressed by In (3.12), the bias term in a 3-D FLC can be realized by adding a fuzzy rule as follows: where is a universal 3-D fuzzy set, whose fuzzy degree is 1 over the space domain for any spatial input , . From (3.12), we can see that an SVM will be equivalent to a 3-D FLC if (3.14) holds.
Finally, a linear SVM regression is employed to learn the consequent part parameters. Using (3.14), the parameters in consequent parts are obtained in terms of the SVM learning, that is,
4. Universal Approximation of Clustering and SVM-Regression Learning-Based 3-D FLC
In essence, the clustering and SVM-regression learning-based 3-D FLC design is a fuzzy modeling that extracts fuzzy control rules and constructs a 3-D FLC from spatiotemporal data hidden with effective control laws. In other words, the proposed 3-D FLC aims at approximating an unknown nonlinear control function. Thus, in this subsection, we are concerned with its universal approximation capability. The universal approximation capability of the SVM learning-based 3-D FLC can be described by the following theorem.
Theorem 4.1. Suppose that the input universe of discourse is a compact set in . Then, for any given real continuous function on and arbitrary , there exists a 3-D FLC as described in (2.4) satisfying the following inequality:
The proof of the theorem is given in the appendix by using Stone-Weierstrass theorem [23]. Theorem 4.1 indicates that the clustering and SVM-regression learning-based 3-D FLC is a universal approximator, that is, it can approximate continuous control functions to arbitrary accuracy.
5. Application
5.1. A Catalytic Packed-Bed Reactor
We take a catalytic packed-bed reactor [1, 5] as an example. The reactor is long and thin as shown in Figure 4. It is fed with gaseous reactant from the right side, and the zero-order gas phase reaction is carried out on the catalyst. The reaction is endothermic, and a jacket is used to heat the reactor. A dimensionless model that describes this nonlinear tubular chemical reactor is provided as follows: subject to the boundary conditions where , and u denote the dimensionless temperature of the gas, the catalyst, and jacket, respectively. The values of the process parameters are given as follows:
The concerned control problem is to control the catalyst temperature throughout the reactor to track a spatial reference profile () in order to maintain a desired degree of reaction rate using the measurements of catalyst temperature from five sensing locations and manipulating one spatially distributed heating source (). The mathematical model (5.1)-(5.2) is only for the process simulation for evaluation of the control scheme. The method of lines [24] is used to simulate the model.
In this application, we aim at extracting 3-D fuzzy rules from a spatiotemporal data set using clustering and SVM regression learning algorithm and constructing a complete 3-D FLC without any prior knowledge.
5.2. Design of a Clustering and SVM-Regression Learning-Based 3-D FLC
5.2.1. Spatiotemporal Data Collection
The spatiotemporal input-output data set is collected from the catalytic packed-bed reactor controlled by expert-knowledge-based 3-D FLC [5], where pseudorandom quinary signal (PRQS) [25] with maximum length of 124 as perturbed signal is added to the control input. Each spatiotemporal input-output data pair consists of a spatial error input , a spatial error in change input , and an incremental output , where ; and denotes the qth and q-1th sampling time, respectively. The detailed design of the expert-knowledge-based 3-D FLC, including fuzzification, 3-D rule inference, and defuzzification, can refer to [5]. The scaling factors for the spatial error, the spatial error in change, and the incremental output are set as 2.0, 0.001, and 0.8716, respectively. The parameters of PRQS are chosen with the following settings: the number of the levels is 5, the length of the period is 124, the sampling time is 0.2 s, and the minimum switching time (i.e., clock period) is 0.2 s.
Two groups of data sets are obtained by adding PRQS signal with different scaling factor (i.e., 0.447 and 0.1) to the control input. The first group with 150 data pairs is generated for training by adding PRQS perturbation signal with a scaling factor 0.447, and the other group with 150 data pairs is generated for test by adding PRQS perturbation signal with a scaling factor 0.1. To evaluate the performance, we employ the following root-mean-squared error (RMSE) as the criteria: where denotes the number of samples, denotes actual output, and denotes expected output.
5.2.2. Design of a Clustering and SVM-Regression Learning-Based 3-D FLC
The design procedure of the proposed 3-D FLC is given as follows: (i)Employ the nearest neighborhood clustering algorithm to deal with the spatiotemporal data set for the input space partition with and , and then generate 16 3-D fuzzy rules with 32 3-D fuzzy sets, where the width of Gaussian type 3-D fuzzy sets is from (3.5).(ii)Simplify the 3-D fuzzy sets and 3-D fuzzy rules based on similarity measure (as described in Section 3.2) with , , and , and then obtain 15 3-D fuzzy rules with 15 3-D fuzzy sets. The distributions of Gaussian type 3-D fuzzy sets are shown in Figure 5.(iii)SVM algorithm described in Section 3.2 is used to learn the consequent part parameters with and . The RMSE in (5.5) for training and test are listed in Table 1. From Table 1, we can find that: smaller yielded more support vectors and led to reasonable training and test performance; while larger yielded less support vectors and led to worse training and test performance. almost had no influence on the training and test performance, once was fixed. In this study, we choose and . Finally, a complete 3-D FLC is constructed with 15 3-D fuzzy rules and 15 3-D fuzzy sets as shown in Figure 6. Using the linguistic hedges approach [14, 21], we can interpret these 3-D fuzzy rules using linguistic words. For instance, the first 3 fuzzy rules are interpreted as follows.(a) IF is less than POSITIVE SMALL and is more than POSITIVE SMALL, THEN is sort of POSITIVE MEDIUM. (b) IF is very ZERO and is very NEGATIVE SMALL, THEN is very ZERO. (c) IF is sort of POSITIVE SMALL and is more than POSITIVE MEDIUM, THEN is more than POSITIVE MEDIUM.
(a)
(b)
(a)
(b)
5.2.3. Control-Performance Validation
The designed clustering and SVM regression learning-based 3-D FLC is applied to the control of the catalytic packed-bed reactor, where simulation time is 10 s. We select the same quantitative performance criteria as in [5]: steady-state error (SSE), integral of the absolute error (IAE), and integral of time multiplied by absolute error (ITAE). The control performance is given in Table 2, and the control profile is given in Figures 7 and 8, where (a), (b), and (c) represent catalyst temperature evolution profile, manipulated input, and catalyst temperature profiles in steady state, respectively. We can find that the proposed 3-D FLC has comparable control performance to the expert-knowledge-based 3-D FLC in [5] both in ideal condition and in disturbed condition.
(a)
(b)
(c)
(a)
(b)
(c)
In addition, we do more control experiments when the SVM-learning algorithm adopts different and . According to the experimental results (see the last three columns in Table 1), we can find that the proposed 3-D FLC shows good control performance when a smaller is chosen.
The above simulation results demonstrate that the proposed design method of a clustering and SVM-regression learning-based 3-D FLC is effective. It provides a beneficial complementary design method to 3-D FLCs.
6. Conclusions
In this paper, we have proposed a new 3-D FLC design methodology based on clustering and SVM regression learning from a spatiotemporal data set. The 3-D FLC design is divided into three steps. Firstly, an initial rule structure is extracted by a nearest neighborhood clustering method, which is modified to be suitable for spatio-temporal data. Secondly, the initial structure is simplified via using similarity measure technique, which is defined for 3-D fuzzy sets and 3-D fuzzy rules. Thirdly, the parameters of the rule consequent parts are learned by a spatial fuzzy basis function-based SVM regression learning algorithm. Besides, the universal approximation capability of the proposed 3-D fuzzy system is discussed. Finally, effectiveness of the proposed 3-D FLC design methodology is validated on a catalytic packed-bed reactor.
Appendix
A. Proof of the Clustering and SVM-Regression Learning-Based 3-D FLC as a Universal Approximator
Let be a set of 3-D FLCs defined in , which is a compact set in . Then, Preliminary 1 is given as follows.
Preliminary 1
Let be a semimetric [26] with the following definition
Therefore, is a metric space. Since there is at least one fuzzy rule in the rule base of a 3-D FLC, is non-empty. Thus, is strictly defined.
Subsequently, we will prove that is dense in using Stone-Weierstrass theorem, where is a set of real continuous functions defined in a compact set . The Stone-Weierstrass theorem is first stated here as follows.
Stone-Weierstrass Theorem (see [16, 23])
Let be a set of real continuous functions on a compact set . If is an algebra, that is, the set is closed under addition, multiplication, and scalar multiplication; separates points on , that is, for every , there exists such that ; and vanishes at no point of , that is, for each there exists such that ; then, the uniform closure of consists of all real continuous functions on , that is, is dense in .
Proof. Firstly, we prove is an algebra. Let , then we can write them as Subsequently, we have three derivation procedures.
(i) Addition
where denotes . Equation (A.3) has the same form as (3.12), then .
(ii) Multiplication
where denotes . In terms of algebraic operation, the product of functions in Gaussian form is also a function in Gaussian form. Thus, (A.4) has the same form as (3.12), and .
(iii) Scalar Multiplication
For arbitrary , we have
Equation (A.5) has the same form as (3.12), and .
Finally, by combining (A.3)~(A.5) together, we can conclude that is an algebra.
Secondly, we will prove that separates point on by constructing a simple 3-D FLC as in (3.12), namely, holds for arbitrarily given with .
We choose two fuzzy rules, that is, .
Let
We have
Since , there must be some and such that . Thus, we have . For arbitrary , holds, therefore, we have. If we choose and , then
Therefore, separates point on .
Finally, we prove vanishes at no point of .
For any 3-D FLC expressed as in (3.12), if we choose and , then for any , we have .
Therefore, vanishes at no point of .
By combining the results from to together, Theorem 4.1 is proven.
Acknowledgments
This work was supported partly by the National Specialised Research Fund for the Doctoral Programme of Higher Education under Grant 20113705120003 and by Shandong Natural Science Foundation under Grant ZR2010FM018 and under Grant ZR2010FM022.