# Reference Function Based Spatiotemporal Fuzzy Logic Control Design Using Support Vector Regression Learning

Authors: Xian-Xia Zhang; Ye Jiang; Shiwei Ma; Bing Wang
Journal: Journal of Applied Mathematics
 (2013)

---

## Abstract

This paper presents a reference function based 3D FLC design methodology using support vector regression (SVR) learning. The concept of reference function is introduced to 3D FLC for the generation of 3D membership functions (MF), which enhance the capability of the 3D FLC to cope with more kinds of MFs. The nonlinear mathematical expression of the reference function based 3D FLC is derived, and spatial fuzzy basis functions are defined. Via relating spatial fuzzy basis functions of a 3D FLC to kernel functions of an SVR, an equivalence relationship between a 3D FLC and an SVR is established. Therefore, a 3D FLC can be constructed using the learned results of an SVR. Furthermore, the universal approximation capability of the proposed 3D fuzzy system is proven in terms of the finite covering theorem. Finally, the proposed method is applied to a catalytic packed-bed reactor and simulation results have verified its effectiveness.


---

## Body

## 1. Introduction

Many industrial processes and systems are “distributed’’ in space [1] and are usually called spatially distributed systems. Recently, a novel three-dimensional fuzzy logic controller (3D FLC) [2] has been developed for the control of such systems. The 3D FLC employs a three-dimensional (3D) fuzzy set [2], which is composed of the traditional fuzzy set plus a third dimension for the spatial information, and carries out a 3D rule inference engine; thus, it has the inherent capability to process spatiotemporal dynamic systems. The control strategy of the 3D FLC is similar to how human operators or experts control the temperature in a space domain. Actually, it is a kind of spatiotemporal fuzzy control system with the traditional model-free advantage.Currently, most 3D FLC designs are based on expert knowledge [2–5], which requires that the human knowledge to the control solution must exist and be structured [6]. However, in many real-world applications, experts may have problems structuring the knowledge. 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 [6]. On the contrast, data sets hidden with effective control rules are usually available. The motivation of this study is to design a 3D FLC using spatiotemporal data information.To date, few literatures are found to be focused on data-based 3D FLC design methods. In [7], a table look-up scheme was employed to design 3D FLC in terms of input-output pairs. In [8], a fuzzy c-means algorithm (FCM) and gradient-descent approach were used to design a data-based 3D FLC, where FCM was used to learn the initial 3D fuzzy rule base and then the gradient-descent approach was used to optimize the parameters of MFs. In [9, 10], a clustering and linear support vector regression based 3D FLC design method was proposed, where the nearest neighborhood clustering was used to construct the antecedent part of 3D fuzzy rules and a linear support vector regression (SVR) was used to learn the consequent parameters. These methods either yield lots of fuzzy control rules (e.g. in [7]) or require additional algorithm to reduce redundant 3D fuzzy sets or 3D fuzzy rules [8–10]. As a complementary, Zhang et al. proposed a data-based 3D FLC design method using SVR learning [11], where the learned support vectors and associated learning parameters are directly used to design antecedent part and consequent part of 3D fuzzy rules. The best advantage of this method is that reasonable 3D fuzzy control rules are directly extracted and constructed by SVR learning. The limitation of the design is that Gaussian shape membership function (MF) is the only choice for MF design.In this study, we focus on a reference function based 3D FLC design using SVR learning, which integrates the merits of SVR learning and flexible MF choice. Utilizing the concept of reference function, the 3D FLC can cope with more kinds of MFs, for example, Symmetric triangle, Gaussian, Cauchy, Laplace, Hyperbolic Secant, and Squared Sinc. A nonlinear mathematical description of a reference function based 3D FLC can be derived, and spatial fuzzy basis functions are defined. Via relating spatial fuzzy basis functions of a 3D FLC to kernel functions (KFs) of an SVR, an equivalence relationship between a 3D FLC and an SVR is established. Therefore, a 3D FLC can be constructed using the learned results of an SVR. In addition, the universal approximation capability of the proposed 3D fuzzy system is proven in terms of the finite covering theorem.The paper is organized as follows. Preliminaries about the reference function, 3D MF generated by reference function, and the nonlinear mapping of reference function based 3D FLC are addressed in Section2. In Section 3, the methodology and design scheme of the reference function based 3D FLC design using SVR learning are presented. Then, the finite covering theorem is used to prove that the 3D FLC is a universal approximator in Section 4. In Section 5, a catalytic packed-bed reactor is presented as an example to illustrate the proposed 3D FLC and validate its effectiveness. In Section 6, conclusions are given.


## 2. Preliminaries

### 2.1. Reference Function

Definition of Reference Function (see [12, 13]). A function υ:R→[0,1] is a reference function if and only if the following two conditions hold:(1)
υ(x)=υ(-x),(2)
υ(0)=1.Many functions may be reference functions. For instance, Symmetric triangle, Gaussian, Cauchy, Laplace, Hyperbolic Secant, and Squared Sinc as listed in Table 1 are reference functions. The reference functions can be used to generate 3D MFs, which provide a way for 3D FLC to access more kinds of 3D MFs.Table 1
Reference functions.
Classification
Mathematical expression
Symmetric triangle
υ
(
x
)
=
max
(
1
-
d
|
x
|
,
0
), d>0
Gaussian
υ
(
x
)
=
e
-
d
x
2, d>0
Cauchy
υ
(
x
)
=
1
1
+
d
x
2, d>0
Laplace
υ
(
x
)
=
e
-
d
|
x
|, d>0
Hyperbolic Secant
υ
(
x
)
=
2
e
d
x
+
e
-
d
x, d>0
Squared Sinc
υ
(
x
)
=
si
n
2
(
d
x
)
d
2
x
2, d>0


### 2.2. Reference Function Based 3D FLC

#### 2.2.1. 3D MF Generated by Reference Function

A 3D MF is an extension of a traditional MF by adding a third coordinate for the spatial information. In detail, the 3D MF has three coordinates: one is for the universe of discourse of the variable, another one is for the spatial information, and the third one is for the membership degree. If finite sensors are used, the 3D MF can be considered as the assembly of the traditional 2D MFs at each sensing location. In this way, we can generate a 3D MF by location transformation of a reference function at each sensing location.For example, we have a spatial input variablex1(z) defined in a discrete spatial domain Z={z1,z2,…,zp}. A 3D MF of x1(z) can be an assembly of the traditional 2D MFs at each sensing location z=zj  (j=1,2,…,p). The MF distribution of x1(z) at sensing location z=z1 can be shown in Figure 1, where β1(x1(z1)) is a Gaussian type reference function; β11(x1(z1)),  β12(x1(z1)), and β13(x1(z1)) are generated by location transformation of β1(x1(z1)).Figure 1
Gaussian MF distribution of spatial input variablex1(z) at sensing location z=z1. β11(x1(z1)),  β12(x1(z1)), and β13(x1(z1)) are generated by the Gaussian type reference function β1(x1(z1)).


#### 2.2.2. Reference Function Based 3D FLC as a Nonlinear Mapping

The basic structure of a 3D FLC is composed of 3D fuzzifier, 3D rule inference, and defuzzifier. Due to its unique 3D nature, some detailed operations of a 3D FLC are different from a traditional one for spatial information expression, processing, and compression. For their detailed operations, one can refer to [2]. Once each component of a reference function based 3D FLC is set, the nonlinear mathematical description of the 3D FLC can be derived (see Appendix A for a brief derivation). Assuming that we employ 3D singleton fuzzifier, 3D fuzzy rules as shown in (A.4) of Appendix A, “product’’ t-norm and “weighted aggregation’’ dimension reduction [3] in the 3D rule inference, singleton fuzzy sets for the output variable, and “linear’’ defuzzifier [14], the reference function based 3D FLC can be mathematically expressed as
(1)u(xz)=b0+∑l=1Nbl∑j=1paj∏i=1sβil(xi(zj))=b0+∑l=1Nbl∑j=1paj∏i=1sβi(xi(zj)-τijl),
where xi(zj) denotes the input of the ith spatial input variable xi(z) from the sensing location z=zj; βil(xi(zj))=βi(xi(zj)-τijl) denotes MF generated by the reference function βi(xi(zj)); τijl denotes a location parameter, that is, the location transformation of the MF of xi(zj) with respect to βi(xi(zj)); aj denotes the spatial weight from thejth sensing location; b0 and bl are constants; p denotes sensor number.In (1), let
(2)Ψl(xz)=∑j=1paj∏i=1sβil(xi(zj));
then (1) can be rewritten as
(3)u(xz)=b0+∑l=1NblΨl(xz).
We define Ψl(xz) as a Spatial Fuzzy Basis Function (SFBF) [11]. Each SFBF corresponds to a 3D fuzzy rule, and all the SFBFs correspond to a 3D rule base. Mathematically, a 3D FLC is a linear combination of all the SFBFs. Furthermore, we rewrite (2) into (4)
(4)Ψl(xz)=∑j=1pajφl(x(zj)),
where φl(x(zj))=∏i=1sβil(xi(zj)).From (4), we can find that, at each sensing location, there exists a traditional FBF [15] φl(x(zj)); in the whole space domain, multiple traditional FBFs are assembled by the spatial weights a1,…,ap into a SFBF Ψl(xz). All the spatial information expression and processing as well as the fuzzy linguistic expression and rule inference are integrated into SFBFs.Equation (1) (or (3)) shows that the reference function based 3D FLC is a nonlinear mapping from the input space xz∈Ω⊂Rp×s to the output space u(xz)∈U⊂R. In particular, using (3) a reference function based 3D FLC can be represented by a three-layer network structure as show in Figure 2.Figure 2
Three-layer network structure of a 3D FLC.


## 2.1. Reference Function

Definition of Reference Function (see [12, 13]). A function υ:R→[0,1] is a reference function if and only if the following two conditions hold:(1)
υ(x)=υ(-x),(2)
υ(0)=1.Many functions may be reference functions. For instance, Symmetric triangle, Gaussian, Cauchy, Laplace, Hyperbolic Secant, and Squared Sinc as listed in Table 1 are reference functions. The reference functions can be used to generate 3D MFs, which provide a way for 3D FLC to access more kinds of 3D MFs.Table 1
Reference functions.
Classification
Mathematical expression
Symmetric triangle
υ
(
x
)
=
max
(
1
-
d
|
x
|
,
0
), d>0
Gaussian
υ
(
x
)
=
e
-
d
x
2, d>0
Cauchy
υ
(
x
)
=
1
1
+
d
x
2, d>0
Laplace
υ
(
x
)
=
e
-
d
|
x
|, d>0
Hyperbolic Secant
υ
(
x
)
=
2
e
d
x
+
e
-
d
x, d>0
Squared Sinc
υ
(
x
)
=
si
n
2
(
d
x
)
d
2
x
2, d>0


## 2.2. Reference Function Based 3D FLC

### 2.2.1. 3D MF Generated by Reference Function

A 3D MF is an extension of a traditional MF by adding a third coordinate for the spatial information. In detail, the 3D MF has three coordinates: one is for the universe of discourse of the variable, another one is for the spatial information, and the third one is for the membership degree. If finite sensors are used, the 3D MF can be considered as the assembly of the traditional 2D MFs at each sensing location. In this way, we can generate a 3D MF by location transformation of a reference function at each sensing location.For example, we have a spatial input variablex1(z) defined in a discrete spatial domain Z={z1,z2,…,zp}. A 3D MF of x1(z) can be an assembly of the traditional 2D MFs at each sensing location z=zj  (j=1,2,…,p). The MF distribution of x1(z) at sensing location z=z1 can be shown in Figure 1, where β1(x1(z1)) is a Gaussian type reference function; β11(x1(z1)),  β12(x1(z1)), and β13(x1(z1)) are generated by location transformation of β1(x1(z1)).Figure 1
Gaussian MF distribution of spatial input variablex1(z) at sensing location z=z1. β11(x1(z1)),  β12(x1(z1)), and β13(x1(z1)) are generated by the Gaussian type reference function β1(x1(z1)).


### 2.2.2. Reference Function Based 3D FLC as a Nonlinear Mapping

The basic structure of a 3D FLC is composed of 3D fuzzifier, 3D rule inference, and defuzzifier. Due to its unique 3D nature, some detailed operations of a 3D FLC are different from a traditional one for spatial information expression, processing, and compression. For their detailed operations, one can refer to [2]. Once each component of a reference function based 3D FLC is set, the nonlinear mathematical description of the 3D FLC can be derived (see Appendix A for a brief derivation). Assuming that we employ 3D singleton fuzzifier, 3D fuzzy rules as shown in (A.4) of Appendix A, “product’’ t-norm and “weighted aggregation’’ dimension reduction [3] in the 3D rule inference, singleton fuzzy sets for the output variable, and “linear’’ defuzzifier [14], the reference function based 3D FLC can be mathematically expressed as
(1)u(xz)=b0+∑l=1Nbl∑j=1paj∏i=1sβil(xi(zj))=b0+∑l=1Nbl∑j=1paj∏i=1sβi(xi(zj)-τijl),
where xi(zj) denotes the input of the ith spatial input variable xi(z) from the sensing location z=zj; βil(xi(zj))=βi(xi(zj)-τijl) denotes MF generated by the reference function βi(xi(zj)); τijl denotes a location parameter, that is, the location transformation of the MF of xi(zj) with respect to βi(xi(zj)); aj denotes the spatial weight from thejth sensing location; b0 and bl are constants; p denotes sensor number.In (1), let
(2)Ψl(xz)=∑j=1paj∏i=1sβil(xi(zj));
then (1) can be rewritten as
(3)u(xz)=b0+∑l=1NblΨl(xz).
We define Ψl(xz) as a Spatial Fuzzy Basis Function (SFBF) [11]. Each SFBF corresponds to a 3D fuzzy rule, and all the SFBFs correspond to a 3D rule base. Mathematically, a 3D FLC is a linear combination of all the SFBFs. Furthermore, we rewrite (2) into (4)
(4)Ψl(xz)=∑j=1pajφl(x(zj)),
where φl(x(zj))=∏i=1sβil(xi(zj)).From (4), we can find that, at each sensing location, there exists a traditional FBF [15] φl(x(zj)); in the whole space domain, multiple traditional FBFs are assembled by the spatial weights a1,…,ap into a SFBF Ψl(xz). All the spatial information expression and processing as well as the fuzzy linguistic expression and rule inference are integrated into SFBFs.Equation (1) (or (3)) shows that the reference function based 3D FLC is a nonlinear mapping from the input space xz∈Ω⊂Rp×s to the output space u(xz)∈U⊂R. In particular, using (3) a reference function based 3D FLC can be represented by a three-layer network structure as show in Figure 2.Figure 2
Three-layer network structure of a 3D FLC.


## 2.2.1. 3D MF Generated by Reference Function

A 3D MF is an extension of a traditional MF by adding a third coordinate for the spatial information. In detail, the 3D MF has three coordinates: one is for the universe of discourse of the variable, another one is for the spatial information, and the third one is for the membership degree. If finite sensors are used, the 3D MF can be considered as the assembly of the traditional 2D MFs at each sensing location. In this way, we can generate a 3D MF by location transformation of a reference function at each sensing location.For example, we have a spatial input variablex1(z) defined in a discrete spatial domain Z={z1,z2,…,zp}. A 3D MF of x1(z) can be an assembly of the traditional 2D MFs at each sensing location z=zj  (j=1,2,…,p). The MF distribution of x1(z) at sensing location z=z1 can be shown in Figure 1, where β1(x1(z1)) is a Gaussian type reference function; β11(x1(z1)),  β12(x1(z1)), and β13(x1(z1)) are generated by location transformation of β1(x1(z1)).Figure 1
Gaussian MF distribution of spatial input variablex1(z) at sensing location z=z1. β11(x1(z1)),  β12(x1(z1)), and β13(x1(z1)) are generated by the Gaussian type reference function β1(x1(z1)).


## 2.2.2. Reference Function Based 3D FLC as a Nonlinear Mapping

The basic structure of a 3D FLC is composed of 3D fuzzifier, 3D rule inference, and defuzzifier. Due to its unique 3D nature, some detailed operations of a 3D FLC are different from a traditional one for spatial information expression, processing, and compression. For their detailed operations, one can refer to [2]. Once each component of a reference function based 3D FLC is set, the nonlinear mathematical description of the 3D FLC can be derived (see Appendix A for a brief derivation). Assuming that we employ 3D singleton fuzzifier, 3D fuzzy rules as shown in (A.4) of Appendix A, “product’’ t-norm and “weighted aggregation’’ dimension reduction [3] in the 3D rule inference, singleton fuzzy sets for the output variable, and “linear’’ defuzzifier [14], the reference function based 3D FLC can be mathematically expressed as
(1)u(xz)=b0+∑l=1Nbl∑j=1paj∏i=1sβil(xi(zj))=b0+∑l=1Nbl∑j=1paj∏i=1sβi(xi(zj)-τijl),
where xi(zj) denotes the input of the ith spatial input variable xi(z) from the sensing location z=zj; βil(xi(zj))=βi(xi(zj)-τijl) denotes MF generated by the reference function βi(xi(zj)); τijl denotes a location parameter, that is, the location transformation of the MF of xi(zj) with respect to βi(xi(zj)); aj denotes the spatial weight from thejth sensing location; b0 and bl are constants; p denotes sensor number.In (1), let
(2)Ψl(xz)=∑j=1paj∏i=1sβil(xi(zj));
then (1) can be rewritten as
(3)u(xz)=b0+∑l=1NblΨl(xz).
We define Ψl(xz) as a Spatial Fuzzy Basis Function (SFBF) [11]. Each SFBF corresponds to a 3D fuzzy rule, and all the SFBFs correspond to a 3D rule base. Mathematically, a 3D FLC is a linear combination of all the SFBFs. Furthermore, we rewrite (2) into (4)
(4)Ψl(xz)=∑j=1pajφl(x(zj)),
where φl(x(zj))=∏i=1sβil(xi(zj)).From (4), we can find that, at each sensing location, there exists a traditional FBF [15] φl(x(zj)); in the whole space domain, multiple traditional FBFs are assembled by the spatial weights a1,…,ap into a SFBF Ψl(xz). All the spatial information expression and processing as well as the fuzzy linguistic expression and rule inference are integrated into SFBFs.Equation (1) (or (3)) shows that the reference function based 3D FLC is a nonlinear mapping from the input space xz∈Ω⊂Rp×s to the output space u(xz)∈U⊂R. In particular, using (3) a reference function based 3D FLC can be represented by a three-layer network structure as show in Figure 2.Figure 2
Three-layer network structure of a 3D FLC.


## 3. Reference Function Based 3D FLC Design Using SVR Learning

### 3.1. Design Methodology

The design methodology can be depicted by Figure3. The SFBFs from a 3D FLC are input to an SVR as the KFs, and the learned spatial support vectors as leading spatiotemporal data points from the SVR are imported for the design of a 3D fuzzy rule base. The design theory will involve two crucial issues. The first one is whether SFBFs from a 3D FLC can be used as KFs in an SVR. The second one is what the relationship between a 3D FLC and an SVR would be like on the basis of the first issue.Figure 3
Design methodology of reference function based 3D FLC design using SVR learning.


#### 3.1.1. Spatial Fuzzy Basis Function as Mercer KF

When relating the SFBFs with the KFs in an SVR, for instance, SFBFs are regarded as KFs, the SVR and the 3D FLC will have the same network structures and then have the same mathematical expressions, which will be discussed in Section3.1.2. Generally speaking, a function satisfying Mercer theorem can be used as a KF for an SVR [16]. In this study, we will prove that an SFBF is a Mercer KF.In (3), we rewrite Ψl(xz) into K(xz,τl), which can be further expressed as follows:
(5)K(xz,τl)=K(xz-τl)=∑j=1pajKj(x(zj),τjl),Kj(x(zj),τjl)=Kj(x(zj)-τjl)=∏i=1sβi(xi(zj)-τijl),
where xz∈Rp×s is a spatial input, τl∈Rp×s is the location transformation parameter of 3D MF in the lth rule, and Kj(x(zj),τjl) and K(xz,τl) are translation invariant KFs [13].In terms of [13], if the reference functions are positive definite functions, then we do get a Mercer kernel. The reference functions as listed in Table 1 are positive definite functions. Using these reference functions to generate MF, from [13], we can conclude that Kj(x(zj),τjl) is a Mercer kernel. Since the linear combination of KFs is still a KF [16], we can derive that K(xz,τl) is still a Mercer KF. Therefore, SFBFs are Mercer KFs, which can be used as KFs for SVR learning.


#### 3.1.2. Mathematical Equivalence of a Spatial SVR and a 3D FLC

Once the SFBFs from the 3D FLC are employed as the KFs for an SVR, an inherent equivalence relationship will be built between the SVR and the 3D FLC (comparing Figure2 and Figure 10). By combining (1) and (B.7), we have the following mathematical expressions:
(6)u(xz)=b0+∑l=1NblΨl(xz)=b0+∑l=1Nbl∑j=1pajφl(x(zj))=b+∑l=1N(αl*-αl)K(xz,xzl),
where b0=b, bl=(αl*-αl), and Ψl(xz)=K(xz,τl)=K(xz,xzl).From (6), we can find that each spatial support vector xzl and its associated learning parameter (αl*-αl) correspond to one 3D fuzzy rule, where xzl is applied to set the center of the 3D MF of the 3D fuzzy set C-il  (i=1,…,s) in the lth rule, that is, the location transformation of the 3D fuzzy set with respect to reference function β(xz,τl), and (αl*-αl) is used to set bl (the constant for the consequent set of the lth rule in 3D FLC).


### 3.2. Design Scheme

The design of a reference function based 3D FLC consists of five parts: data collection, KF generation, SVR learning, 3D fuzzy rule construction, and 3D fuzzy controller integration, as shown in Figure4.Figure 4
Design scheme of reference function based 3D FLC using SVR learning.(
1
) Data Collection. A set of spatiotemporal data will be collected. The data should contain effective control laws. Essentially, the reference function based 3D FLC design is a fuzzy modeling [17] that extracts fuzzy control rules from the spatiotemporal data.(
2
) KF Generation. Before SVR learning, KFs should be properly designed. In this step, via properly selecting reference function, SFBFs from a 3D FLC will be formulated (as in (4)) to set KFs for SVR learning.(
3) SVR Learning. With proper KFs, the SVM learning algorithm directly executes the spatiotemporal data set and yields spatial support vectors xz1,…,xzN and their associated learning parameters α1*-α1,…,αN*-αN.(
4) 3D Fuzzy Rule Construction. The spatial support vectors and their associated learning parameters, as leading control laws, are used to construct 3D fuzzy control rules. In detail, the spatial support vector xzl is employed to construct the antecedent part of the lth rule; αl*-αl is employed to construct the consequent part of the lth rule. The form of each 3D fuzzy rule is shown as below
(7)R-l:ifxzisclosetoxzlthenuiscloseto(αl*-αl).
It is shown that the result of the SVM learning can be easily interpreted using structured linguistic knowledge. Finally, a 3D rule base with N rules is established.(
5) 3D Fuzzy Controller Integration. Once the 3D rule base is established, a 3D FLC can be achieved by integrating other components including 3D fuzzifier, 3D rule inference, and defuzzifier. The detailed settings are given in Section 2.2.2. Finally, we obtain a complete 3D FLC, which can be used as a controller for a spatially distributed dynamic system.


## 3.1. Design Methodology

The design methodology can be depicted by Figure3. The SFBFs from a 3D FLC are input to an SVR as the KFs, and the learned spatial support vectors as leading spatiotemporal data points from the SVR are imported for the design of a 3D fuzzy rule base. The design theory will involve two crucial issues. The first one is whether SFBFs from a 3D FLC can be used as KFs in an SVR. The second one is what the relationship between a 3D FLC and an SVR would be like on the basis of the first issue.Figure 3
Design methodology of reference function based 3D FLC design using SVR learning.


### 3.1.1. Spatial Fuzzy Basis Function as Mercer KF

When relating the SFBFs with the KFs in an SVR, for instance, SFBFs are regarded as KFs, the SVR and the 3D FLC will have the same network structures and then have the same mathematical expressions, which will be discussed in Section3.1.2. Generally speaking, a function satisfying Mercer theorem can be used as a KF for an SVR [16]. In this study, we will prove that an SFBF is a Mercer KF.In (3), we rewrite Ψl(xz) into K(xz,τl), which can be further expressed as follows:
(5)K(xz,τl)=K(xz-τl)=∑j=1pajKj(x(zj),τjl),Kj(x(zj),τjl)=Kj(x(zj)-τjl)=∏i=1sβi(xi(zj)-τijl),
where xz∈Rp×s is a spatial input, τl∈Rp×s is the location transformation parameter of 3D MF in the lth rule, and Kj(x(zj),τjl) and K(xz,τl) are translation invariant KFs [13].In terms of [13], if the reference functions are positive definite functions, then we do get a Mercer kernel. The reference functions as listed in Table 1 are positive definite functions. Using these reference functions to generate MF, from [13], we can conclude that Kj(x(zj),τjl) is a Mercer kernel. Since the linear combination of KFs is still a KF [16], we can derive that K(xz,τl) is still a Mercer KF. Therefore, SFBFs are Mercer KFs, which can be used as KFs for SVR learning.


### 3.1.2. Mathematical Equivalence of a Spatial SVR and a 3D FLC

Once the SFBFs from the 3D FLC are employed as the KFs for an SVR, an inherent equivalence relationship will be built between the SVR and the 3D FLC (comparing Figure2 and Figure 10). By combining (1) and (B.7), we have the following mathematical expressions:
(6)u(xz)=b0+∑l=1NblΨl(xz)=b0+∑l=1Nbl∑j=1pajφl(x(zj))=b+∑l=1N(αl*-αl)K(xz,xzl),
where b0=b, bl=(αl*-αl), and Ψl(xz)=K(xz,τl)=K(xz,xzl).From (6), we can find that each spatial support vector xzl and its associated learning parameter (αl*-αl) correspond to one 3D fuzzy rule, where xzl is applied to set the center of the 3D MF of the 3D fuzzy set C-il  (i=1,…,s) in the lth rule, that is, the location transformation of the 3D fuzzy set with respect to reference function β(xz,τl), and (αl*-αl) is used to set bl (the constant for the consequent set of the lth rule in 3D FLC).


## 3.1.1. Spatial Fuzzy Basis Function as Mercer KF

When relating the SFBFs with the KFs in an SVR, for instance, SFBFs are regarded as KFs, the SVR and the 3D FLC will have the same network structures and then have the same mathematical expressions, which will be discussed in Section3.1.2. Generally speaking, a function satisfying Mercer theorem can be used as a KF for an SVR [16]. In this study, we will prove that an SFBF is a Mercer KF.In (3), we rewrite Ψl(xz) into K(xz,τl), which can be further expressed as follows:
(5)K(xz,τl)=K(xz-τl)=∑j=1pajKj(x(zj),τjl),Kj(x(zj),τjl)=Kj(x(zj)-τjl)=∏i=1sβi(xi(zj)-τijl),
where xz∈Rp×s is a spatial input, τl∈Rp×s is the location transformation parameter of 3D MF in the lth rule, and Kj(x(zj),τjl) and K(xz,τl) are translation invariant KFs [13].In terms of [13], if the reference functions are positive definite functions, then we do get a Mercer kernel. The reference functions as listed in Table 1 are positive definite functions. Using these reference functions to generate MF, from [13], we can conclude that Kj(x(zj),τjl) is a Mercer kernel. Since the linear combination of KFs is still a KF [16], we can derive that K(xz,τl) is still a Mercer KF. Therefore, SFBFs are Mercer KFs, which can be used as KFs for SVR learning.


## 3.1.2. Mathematical Equivalence of a Spatial SVR and a 3D FLC

Once the SFBFs from the 3D FLC are employed as the KFs for an SVR, an inherent equivalence relationship will be built between the SVR and the 3D FLC (comparing Figure2 and Figure 10). By combining (1) and (B.7), we have the following mathematical expressions:
(6)u(xz)=b0+∑l=1NblΨl(xz)=b0+∑l=1Nbl∑j=1pajφl(x(zj))=b+∑l=1N(αl*-αl)K(xz,xzl),
where b0=b, bl=(αl*-αl), and Ψl(xz)=K(xz,τl)=K(xz,xzl).From (6), we can find that each spatial support vector xzl and its associated learning parameter (αl*-αl) correspond to one 3D fuzzy rule, where xzl is applied to set the center of the 3D MF of the 3D fuzzy set C-il  (i=1,…,s) in the lth rule, that is, the location transformation of the 3D fuzzy set with respect to reference function β(xz,τl), and (αl*-αl) is used to set bl (the constant for the consequent set of the lth rule in 3D FLC).


## 3.2. Design Scheme

The design of a reference function based 3D FLC consists of five parts: data collection, KF generation, SVR learning, 3D fuzzy rule construction, and 3D fuzzy controller integration, as shown in Figure4.Figure 4
Design scheme of reference function based 3D FLC using SVR learning.(
1
) Data Collection. A set of spatiotemporal data will be collected. The data should contain effective control laws. Essentially, the reference function based 3D FLC design is a fuzzy modeling [17] that extracts fuzzy control rules from the spatiotemporal data.(
2
) KF Generation. Before SVR learning, KFs should be properly designed. In this step, via properly selecting reference function, SFBFs from a 3D FLC will be formulated (as in (4)) to set KFs for SVR learning.(
3) SVR Learning. With proper KFs, the SVM learning algorithm directly executes the spatiotemporal data set and yields spatial support vectors xz1,…,xzN and their associated learning parameters α1*-α1,…,αN*-αN.(
4) 3D Fuzzy Rule Construction. The spatial support vectors and their associated learning parameters, as leading control laws, are used to construct 3D fuzzy control rules. In detail, the spatial support vector xzl is employed to construct the antecedent part of the lth rule; αl*-αl is employed to construct the consequent part of the lth rule. The form of each 3D fuzzy rule is shown as below
(7)R-l:ifxzisclosetoxzlthenuiscloseto(αl*-αl).
It is shown that the result of the SVM learning can be easily interpreted using structured linguistic knowledge. Finally, a 3D rule base with N rules is established.(
5) 3D Fuzzy Controller Integration. Once the 3D rule base is established, a 3D FLC can be achieved by integrating other components including 3D fuzzifier, 3D rule inference, and defuzzifier. The detailed settings are given in Section 2.2.2. Finally, we obtain a complete 3D FLC, which can be used as a controller for a spatially distributed dynamic system.


## 4. Universal Approximation of Reference Function Based 3D FLC

The reference function based 3D FLC design method is used to construct a 3D FLC from spatiotemporal data hidden with effective control laws. In other words, the 3D FLC aims at approximating an unknown nonlinear control function. In this study, we use the finite covering theorem to prove that the 3D FLC is a universal approximator; that is, it can approximate continuous control functions to arbitrary accuracy.The universal approximation capability of the reference function based 3D FLC can be described by Theorem1.Theorem 1.
Letg~(xz):Rp×s→R be a continuous function defined on a compact Ω. For each ε>0, there exists a reference function based 3D FLC u(xz) such that
(8)supxz∈Ω(|u(xz)-g~(xz)|)<ε.From(A.4), it is shown that R-0 is an universal rule, namely, for any spatial input xz,  R-0 will be fired. In (A.4), the fired rule R-0 will produce the constant b0. Let
(9)F(xz)=u(xz)-b0=∑l=1Nbl∑j=1paj∏i=1sβi(xi(zj)-τijl)=∑l=1NblK(xz,τl).F(xz) can be regarded as a 3D FLC generated by rule base {R-1,R-2,…,R-N}. Then, Theorem 1 can be restated as Theorem 2 as follows.Theorem 2.
Under the condition of Theorem1, let g(xz):Rp×s→R be a continuous function defined on a compact Ω. For any constant b0, one has g~(xz)=g(xz)+b0. For each ε>0, there exists a reference function based 3D FLC F(xz) such that
(10)supxz∈Ω(|F(xz)-g(xz)|)<ε.Before the proof of Theorem2, we first present some preparation work.Whenxz=x-z∈Ω, the firing level of the fired rule R-l(l=1,…,N) is
(11)μφl(x-z)=∑j=1paj∏i=1sβi(x-i(zj)-τijl).
The inference result of R-l is given by
(12)Dl(u)={0ifμφl(x-z)=0oru≠bl,blμφl(x-z)inothercase,
where u is the output variable of the 3D FLC, which is corresponding to the “u” of the consequent part of the fired rule. The composition result of all the fired rules is given as follows:
(13)D(u)=⋃l=1NDl(u).Based on the above preparation, Lemma3 is presented as follows.Lemma 3.
Under the condition of Theorem1 there exists a reference function based 3D FLC F(xz) such that
(14)D(u)|u-g(xz)|≤D(u)*εforeachu∈R.Proof.
Letaz∈Ω. As g(·) is continuous at az, for each i=1,…,s there exists a δazi>0 such that
(15)|xzi-azi|≤δazii=(1,…,s)⟺|g(xz)-g(az)|≤ε.
For each az∈Ω, set
(16)Oaz={xz∣|xzi-azi|≤δazi(i=1,…,s)}.
Then, Oazis open on Rp×s and Ω⊆⋃az∈ΩOaz. As Ω is compact, there exists a finite subfamily Oaz1,Oaz2,…,Oazk such that
(17)Ω⊆Oaz1⋃Oaz2⋃⋯⋃Oazk.
We can construct a 3D FLC F(xz), defined by
(18)μφl(xz)={∑j=1pajl∏i=1sβi(xi(zj)-τijl)xz∈Oazl0xz∉Oazl∑j=1pajl∏i=1sβi(xi(zj)-τijl)vvvvvvv1(l=1,…,N),(19)bl={g(xz)xz∈Oazl0xz∉Oazl(l=1,…,N).
When xz=x-z∈Ω, we have the following.(1)
IfD(u)=0, the lemma is trivial.(2)
IfD(u)>0, then D(u)=⋃l=1NDl(u)>0; hence, there exists a k∈[1,N] such that Dk(u)>0. Therefore, we further have that the following.
(a)
FromDk(u)≠0, it follows that μφk(x-z)≠0; that is, x-z∈Oazk. In terms of the continuity of g(·), we have |g(x-z)-g(azk)|≤ε.(b)
FromDk(u)≠0, it follows that Dk(u)=bkμφk(x-z)≠0, and then bk≠0. We have u=bk=g(azk).Hence, we have
(20)|u-g(x-z)|≤|u-g(azk)|+|g(azk)-g(x-z)|<0+ε=ε.In terms of Lemma3, the proof of Theorem 2 can be given as follows.Proof of Theorem2.
Consider(21)|F(x-z)-g(x-z)|≤∑D(u)|u-g(x-z)|∑D(u)≤ε*∑D(u)∑D(u)≤ε.


## 5. Application

### 5.1. A Catalytic Packed-Bed Reactor

This designed 3D FLC is applied to a catalytic packed-bed reactor [1, 4, 18] shown in Figure 5, where a reaction of the form A→B takes place 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:
(22)∂Tr∂t=-∂Tr∂z+1PeT∂2Tr∂z2-BTBCCAexp(γrTr1+Tr)+βT(u-Tr),∂CA∂t=-∂CA∂z+1PeC∂2CA∂z2-BCCAexp(γrTr1+Tr)
subject to the boundary conditions
(23)z=0,PeTTr=∂Tr∂z,PeC(CA-1)=∂CA∂z,z=1,∂Tr∂z=0,∂CA∂z=0,
where Tr, CA, and u denote the dimensionless temperature, the concentration of reactant A, and jacket temperature, respectively; t and z denote the dimensionless time and space; PeT and PeC are the heat and mass Peclet numbers, BT is a dimensionless heat of reaction, BC is a dimensionless preexponential factor, γr is a dimensionless activation energy, and βT is a dimensionless heat transfer coefficient. The values of the process parameters are given as follows:
(24)PeT=5.0,PeC=5.0,BC=0.00001,BT=1.0,βT=15.62,γr=22.14.
The control problem is to maintain a desired reaction rate via tuning the jacket temperature to control catalyst temperature. In this application, the reactor began to work at one steady state; because of the requirement of operation conditions, the reference value of temperature is increased by 8%. Thus, the control objective is to make the temperature of reactor well track the new reference value along the space domain.Figure 5
Sketch of a catalytic packed-bed reactor.


### 5.2. Design of Reference Function Based 3D FLC

(
1) Data Collection. The spatiotemporal data is collected from the catalytic packed-bed reactor, which is controlled by an expert based 3D FLC [4]. Five point sensors are located along the length of the reactor with Z=[00.40.60.81] for collecting the spatial distribution of the temperature Tr. Two spatial inputs are error and error in change; that is, e*(Z,k)=Tsd(Z)-Ts(Z,k) and Δe*(Z)=e*(Z,k)-e*(Z,k-1). The detailed design of the expert based 3D FLC, including 3D fuzzifier, 3D rule inference, and defuzzifier, can refer to [4]. The scaling factors for the spatial error, the spatial error in change, and the incremental output are set as 0.5, 0.1, and 0.3, respectively. The sample period is 0.1 s, and the sampling duration is 6 s. Thus, we have 60 input-output data pairs (shown in Figure 6), each of which is represented by ([e(z)k,Δe(z)k],Δuk),  k=1,…,60.Spatial-temporal data set: (a) spatial errore(z); (b) spatial error in change Δe(z); (c) incremental output Δu (z: spatial dimension; k: serial number of input-output data).
(a)
(b)
(c)(
2) KF Generation. The reference functions (Symmetric triangle, Gaussian, Cauchy, Laplace, Hyperbolic Secant, and Squared Sinc) as listed in Table 1 are employed, respectively, to generate 3D MFs and then formulate SFBFs. SFBFs are used for KFs in an SVR learning.(
3) SVR Learning. With the spatiotemporal data set as above, the SVM learning algorithm is used for the support vector learning. It should be noted that the SFBFs in (3) are taken as the KFs. As a result, spatial support vectors are extracted and their associated learning parameters are obtained. For instance, when the Gaussian type reference function is used for KF generation, eight support vectors were learned from 60 spatiotemporal data pairs when C=1000,  ε=0.00005, and d=0.1, as shown in Figure 7, where the spatiotemporal input data is decomposed into multiple two-dimensional graphical representations over the space domain.Figure 7
Input data and support vectors from each sensing location. Circle (o): support vector; dot (+): input data.(
4) 3D Fuzzy Rule Construction. In terms of the learned results of the SVR in the previous step, we establish 3D fuzzy rules. For instance, with the Gaussian type reference function, eight 3D fuzzy rules are constructed. The first four 3D fuzzy rules are presented as follows.R-1: if e(z) is close to [-0.0580-0.0963-0.0988-0.1000-0.1005]′ and Δe(z) is close to [-0.0580-0.0963-0.0988-0.1000-0.1005]′, then Δu is close to -9.5172.R-2: if e(z) is close to [-0.0479-0.0817-0.0839-0.0848-0.0851]′ and Δe(z) is close to [0.01010.01460.01490.01520.0154]′, then Δu is close to -298.9862.R-3: if e(z) is close to [-0.0414-0.0714-0.0731-0.0735-0.0737]′ and Δe(z) is close to [0.00650.01040.01090.01130.0114]′, then Δu is close to 249.2288.R-4: if e(z) is close to [-0.0359-0.0626-0.0638-0.0639-0.0639]′ and Δe(z) is close to [0.00550.00870.00920.00960.0098]′, then Δu is close to 148.1335.The first four 3D fuzzy rules and their associated 3D MF distributions can be depicted in Figure 8, which show the inherent spatial nature of the 3D fuzzy control system.Figure 8
The first four 3D fuzzy rules and their associated MFs (b=-0.0752).(
5) 3D Fuzzy Controller Integration. Based on the 3D fuzzy rules established in step (4), we obtain a complete 3D FLC by combining 3D fuzzifier, 3D rule inference, and defuzzifier. The resultant 3D FLC will be used as a controller for the catalytic packed-bed reactor.


### 5.3. Control Performance Validation

The designed reference function based 3D FLC using SVR learning is validated on the nonlinear catalytic packed-bed reactor. We employed six kinds of reference functions, that is, Symmetric triangle, Gaussian, Cauchy, Laplace, Hyperbolic Secant, and Squared Sinc, and finally produced six 3D FLCs. With the reference function based 3D FLC as the controller, the catalyst temperature varying with time and space, manipulated input, and the catalyst temperature at steady state are presented in Figures9(a)–9(f). The control performance is given in Table 2, where steady-state error (SSE), integral of the absolute error (IAE), and integral of time multiplied by absolute error (ITAE) [2] are used as the performance criteria. In terms of Figures 9(a)–9(f) and Table 2, we can find that different reference functions will yield different control performance. In this application, Gaussian, Cauchy, Hyperbolic Secant, and Squared Sinc reference functions result in good control performance, while Symmetric triangle and Laplace reference functions lead to poor control performance. The results illustrate that KF will influence the control performance; thus, in the actual application, we should choose proper KF to design a 3D FLC.Table 2
Performance comparisons.
Reference function based 3D FLC
3D fuzzy rules
ISS (×10−3)
IAE (×10−1)
ITAE (×10−1)
Symmetric triangle reference function
7
103.4
16.163
65.457
Gaussian reference function
8
4.7
1.864
3.526
Cauchy reference function
8
4.7
1.881
3.537
Laplace reference function
7
108
16.525
67.408
Hyperbolic Secant reference function
28
4.8
1.918
3.567
Squared Sinc reference function
5
4.7
1.859
3.515
Expert-based 3D FLC
25
4.8
1.931
3.676Control performance comparisons. From left to right in (a)–(g): catalyst temperature varying with time and space, manipulated input, and catalyst temperature at steady state.
(a)
Controlled by a reference function based 3D FLC with Symmetric triangle reference function
(b)
Controlled by a reference function based 3D FLC with Gaussian reference function
(c)
Controlled by a reference function based 3D FLC with Cauchy reference function
(d)
Controlled by a reference function based 3D FLC with Laplace reference function
(e)
Controlled by a reference function based 3D FLC with Hyperbolic Secant reference function
(f)
Controlled by a reference function based 3D FLC with Squared Sinc reference function
(g)
Controlled by an expert based 3D FLCFigure 10
Three-layer network structure of an SVR.As a comparison, we do another control experiment; that is, the expert knowledge-based 3D FLC from [4] is taken as the controller. As for its detailed design including 3D MF, 3D rule base, 3D inference, fuzzification, and defuzzification, one can refer to [4]. The scaling factors for the spatial error, the spatial error in change, and the incremental output are set as 0.5, 0.1, and 0.3, respectively. The controlled profiles and control performance are given in Figure 9(g) and Table 2, respectively.From Figure9 and Table 2, we can find that with a proper reference function, the reference function based 3D FLC has nearly the same control performance as the expert knowledge-based 3D FLC. It means that the proposed spatial SVR learning method can well extract the control laws hidden in a spatiotemporal input–output dataset and formulate them in the form of 3D fuzzy rules.


## 5.1. A Catalytic Packed-Bed Reactor

This designed 3D FLC is applied to a catalytic packed-bed reactor [1, 4, 18] shown in Figure 5, where a reaction of the form A→B takes place 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:
(22)∂Tr∂t=-∂Tr∂z+1PeT∂2Tr∂z2-BTBCCAexp(γrTr1+Tr)+βT(u-Tr),∂CA∂t=-∂CA∂z+1PeC∂2CA∂z2-BCCAexp(γrTr1+Tr)
subject to the boundary conditions
(23)z=0,PeTTr=∂Tr∂z,PeC(CA-1)=∂CA∂z,z=1,∂Tr∂z=0,∂CA∂z=0,
where Tr, CA, and u denote the dimensionless temperature, the concentration of reactant A, and jacket temperature, respectively; t and z denote the dimensionless time and space; PeT and PeC are the heat and mass Peclet numbers, BT is a dimensionless heat of reaction, BC is a dimensionless preexponential factor, γr is a dimensionless activation energy, and βT is a dimensionless heat transfer coefficient. The values of the process parameters are given as follows:
(24)PeT=5.0,PeC=5.0,BC=0.00001,BT=1.0,βT=15.62,γr=22.14.
The control problem is to maintain a desired reaction rate via tuning the jacket temperature to control catalyst temperature. In this application, the reactor began to work at one steady state; because of the requirement of operation conditions, the reference value of temperature is increased by 8%. Thus, the control objective is to make the temperature of reactor well track the new reference value along the space domain.Figure 5
Sketch of a catalytic packed-bed reactor.


## 5.2. Design of Reference Function Based 3D FLC

(
1) Data Collection. The spatiotemporal data is collected from the catalytic packed-bed reactor, which is controlled by an expert based 3D FLC [4]. Five point sensors are located along the length of the reactor with Z=[00.40.60.81] for collecting the spatial distribution of the temperature Tr. Two spatial inputs are error and error in change; that is, e*(Z,k)=Tsd(Z)-Ts(Z,k) and Δe*(Z)=e*(Z,k)-e*(Z,k-1). The detailed design of the expert based 3D FLC, including 3D fuzzifier, 3D rule inference, and defuzzifier, can refer to [4]. The scaling factors for the spatial error, the spatial error in change, and the incremental output are set as 0.5, 0.1, and 0.3, respectively. The sample period is 0.1 s, and the sampling duration is 6 s. Thus, we have 60 input-output data pairs (shown in Figure 6), each of which is represented by ([e(z)k,Δe(z)k],Δuk),  k=1,…,60.Spatial-temporal data set: (a) spatial errore(z); (b) spatial error in change Δe(z); (c) incremental output Δu (z: spatial dimension; k: serial number of input-output data).
(a)
(b)
(c)(
2) KF Generation. The reference functions (Symmetric triangle, Gaussian, Cauchy, Laplace, Hyperbolic Secant, and Squared Sinc) as listed in Table 1 are employed, respectively, to generate 3D MFs and then formulate SFBFs. SFBFs are used for KFs in an SVR learning.(
3) SVR Learning. With the spatiotemporal data set as above, the SVM learning algorithm is used for the support vector learning. It should be noted that the SFBFs in (3) are taken as the KFs. As a result, spatial support vectors are extracted and their associated learning parameters are obtained. For instance, when the Gaussian type reference function is used for KF generation, eight support vectors were learned from 60 spatiotemporal data pairs when C=1000,  ε=0.00005, and d=0.1, as shown in Figure 7, where the spatiotemporal input data is decomposed into multiple two-dimensional graphical representations over the space domain.Figure 7
Input data and support vectors from each sensing location. Circle (o): support vector; dot (+): input data.(
4) 3D Fuzzy Rule Construction. In terms of the learned results of the SVR in the previous step, we establish 3D fuzzy rules. For instance, with the Gaussian type reference function, eight 3D fuzzy rules are constructed. The first four 3D fuzzy rules are presented as follows.R-1: if e(z) is close to [-0.0580-0.0963-0.0988-0.1000-0.1005]′ and Δe(z) is close to [-0.0580-0.0963-0.0988-0.1000-0.1005]′, then Δu is close to -9.5172.R-2: if e(z) is close to [-0.0479-0.0817-0.0839-0.0848-0.0851]′ and Δe(z) is close to [0.01010.01460.01490.01520.0154]′, then Δu is close to -298.9862.R-3: if e(z) is close to [-0.0414-0.0714-0.0731-0.0735-0.0737]′ and Δe(z) is close to [0.00650.01040.01090.01130.0114]′, then Δu is close to 249.2288.R-4: if e(z) is close to [-0.0359-0.0626-0.0638-0.0639-0.0639]′ and Δe(z) is close to [0.00550.00870.00920.00960.0098]′, then Δu is close to 148.1335.The first four 3D fuzzy rules and their associated 3D MF distributions can be depicted in Figure 8, which show the inherent spatial nature of the 3D fuzzy control system.Figure 8
The first four 3D fuzzy rules and their associated MFs (b=-0.0752).(
5) 3D Fuzzy Controller Integration. Based on the 3D fuzzy rules established in step (4), we obtain a complete 3D FLC by combining 3D fuzzifier, 3D rule inference, and defuzzifier. The resultant 3D FLC will be used as a controller for the catalytic packed-bed reactor.


## 5.3. Control Performance Validation

The designed reference function based 3D FLC using SVR learning is validated on the nonlinear catalytic packed-bed reactor. We employed six kinds of reference functions, that is, Symmetric triangle, Gaussian, Cauchy, Laplace, Hyperbolic Secant, and Squared Sinc, and finally produced six 3D FLCs. With the reference function based 3D FLC as the controller, the catalyst temperature varying with time and space, manipulated input, and the catalyst temperature at steady state are presented in Figures9(a)–9(f). The control performance is given in Table 2, where steady-state error (SSE), integral of the absolute error (IAE), and integral of time multiplied by absolute error (ITAE) [2] are used as the performance criteria. In terms of Figures 9(a)–9(f) and Table 2, we can find that different reference functions will yield different control performance. In this application, Gaussian, Cauchy, Hyperbolic Secant, and Squared Sinc reference functions result in good control performance, while Symmetric triangle and Laplace reference functions lead to poor control performance. The results illustrate that KF will influence the control performance; thus, in the actual application, we should choose proper KF to design a 3D FLC.Table 2
Performance comparisons.
Reference function based 3D FLC
3D fuzzy rules
ISS (×10−3)
IAE (×10−1)
ITAE (×10−1)
Symmetric triangle reference function
7
103.4
16.163
65.457
Gaussian reference function
8
4.7
1.864
3.526
Cauchy reference function
8
4.7
1.881
3.537
Laplace reference function
7
108
16.525
67.408
Hyperbolic Secant reference function
28
4.8
1.918
3.567
Squared Sinc reference function
5
4.7
1.859
3.515
Expert-based 3D FLC
25
4.8
1.931
3.676Control performance comparisons. From left to right in (a)–(g): catalyst temperature varying with time and space, manipulated input, and catalyst temperature at steady state.
(a)
Controlled by a reference function based 3D FLC with Symmetric triangle reference function
(b)
Controlled by a reference function based 3D FLC with Gaussian reference function
(c)
Controlled by a reference function based 3D FLC with Cauchy reference function
(d)
Controlled by a reference function based 3D FLC with Laplace reference function
(e)
Controlled by a reference function based 3D FLC with Hyperbolic Secant reference function
(f)
Controlled by a reference function based 3D FLC with Squared Sinc reference function
(g)
Controlled by an expert based 3D FLCFigure 10
Three-layer network structure of an SVR.As a comparison, we do another control experiment; that is, the expert knowledge-based 3D FLC from [4] is taken as the controller. As for its detailed design including 3D MF, 3D rule base, 3D inference, fuzzification, and defuzzification, one can refer to [4]. The scaling factors for the spatial error, the spatial error in change, and the incremental output are set as 0.5, 0.1, and 0.3, respectively. The controlled profiles and control performance are given in Figure 9(g) and Table 2, respectively.From Figure9 and Table 2, we can find that with a proper reference function, the reference function based 3D FLC has nearly the same control performance as the expert knowledge-based 3D FLC. It means that the proposed spatial SVR learning method can well extract the control laws hidden in a spatiotemporal input–output dataset and formulate them in the form of 3D fuzzy rules.


## 6. Conclusions

A reference function based 3D FLC design methodology using SVR learning is proposed for spatially distributed dynamic systems. Utilizing the concept of reference function, the 3D FLC can access more kinds of 3D MFs, such as Symmetric triangle, Gaussian, Cauchy, Laplace, Hyperbolic Secant, and Squared Sinc. Based on the mathematical expressions of reference function based 3D FLC, we define spatial fuzzy basis functions and then find an equivalence relationship between a 3D FLC and an SVR by connecting spatial fuzzy basis functions in the 3D FLC to KFs in the SVR. On the basis of the equivalence relationship, a 3D FLC can be designed using the SVR learning; that is, the learned spatial support vectors as the optimal leading data points can be directly used for 3D fuzzy control rule generation. The proposed reference function based 3D FLC design can be carried out in five steps: data collection, KF generation, SVR learning, 3D fuzzy rule construction, and 3D fuzzy controller integration. Besides, the universal approximation capability of the proposed 3D fuzzy system is discussed. Finally, effectiveness of the proposed 3D FLC design methodology is validated on a catalytic packed-bed reactor.



---

*Source: 410279-2013-12-12.xml*