# Fusion Control of Flexible Logic Control and Neural Network

Authors: Lihua Fu; Dan Wang
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

Based on the basic physical meaning of errorE and error variety EC, this paper analyzes the logical relationship between them and uses Universal Combinatorial Operation Model in Universal Logic to describe it. Accordingly, a flexible logic control method is put forward to realize effective control on multivariable nonlinear system. In order to implement fusion control with artificial neural network, this paper proposes a new neuron model of Zero-level Universal Combinatorial Operation in Universal Logic. And the artificial neural network of flexible logic control model is implemented based on the proposed neuron model. Finally, stability control, anti-interference control of double inverted-pendulum system, and free walking of cart pendulum system on a level track are realized, showing experimentally the feasibility and validity of this method.


---

## Body

## 1. Introduction

In recent years, fuzzy control has made a rapid development, and it has found a considerable number of successful industrial applications [1–3]. But fuzzy control has two shortcomings in the process of controlling some practical complex systems. One is that the number of control rules increases exponentially with the increase of the number of inputs, and the other one is that the precision of control system is low [4].To reduce the dimension of control model, hierarchical fuzzy logic control divides the collection of control rules into several collections based on different functions [5, 6]. Compound control combines fuzzy control and other relatively mature control methods to realize the effective control [7], such as Fuzzy-PID Compound Control [8], fuzzy predication control [9], adaptive fuzzy H
∞ control [10], and so forth. The basic idea of adaptive fuzzy control based on variable universe [11, 12] is to keep the form of rules and varies universe of discourse according to the control error. Though a great deal of research has been done to improve the performance of fuzzy control, most of these methods are based on the basic idea that fuzzy controller is a piecewise approximator. However, to date, there has been relatively little research conducted on the internal relations among input variables of fuzzy controllers.Based on analysis of the logical relationship between the system’s errorE and error variety E
C, this paper indicates that the relationship is just universal combinatorial relation in Universal Logic [13], and the simple Universal Combinatorial Operation can be used instead of complex fuzzy rule-based reasoning process. As a result, a flexible logic control method is proposed to realize effective control on multivariable nonlinear system.Artificial neural network is widely used in modelling and controlling thanks to its properties of self-learning, self-organizing, and self-adapting [14]. In order to realize fusion control with artificial neural network, this paper attempts to study the neuron model of Zero-level Universal Combinatorial Operation in Universal Logic and propose a new neural model. Based on this neuron model, the artificial neural network of flexible logic control model is implemented. Finally, stability control, anti-interference control of double inverted-pendulum system, and free walking of cart pendulum system on a level track are realized to prove the feasibility and validity of this method.The rest of the paper is organized as follows. Section2 introduces necessary background on Universal Combinatorial Operation Model and flexible logic control method and gives and proves some important theorems of Universal Combinatorial Operation Model in the interval [
a
,
b
]. Section 3 puts forward a new neuron model of Zero-level Universal Combinatorial Operation. Based on this neuron model, the artificial neural network of flexible logic control model is implemented in Section 4. The designed flexible logic control model is applied to treat the double inverted-pendulum system in Section 5. Finally, concluding remarks are given in Section 6.


## 2. Universal Combinatorial Operation Model

### 2.1. Universal Combinatorial Operation Model

In order to deal reasonably with the complex relation between factors in complex system,T-norm, S-norm or Mean operators are taken as Aggregation Operators.However,T-norm result is not bigger than the minimum value, and S-norm result is not less than the maximum value. As a result, T-norm, or S-norm can only handle mutually conflictive relation. In contrast, Mean operators can vary only between the minimum and maximum values based on its “tradeoff” concept, so it can only handle mutually consistent relation [13].Universal Logic [13], proposed by Professor He et al., is a kind of flexible logic. It considers the continuous change of not only the truth value of propositions, which is called truth value flexibility, but also the relation between propositions, which is called relation flexibility. Based on fuzzy logic, it puts forward two important coefficients: generalized correlation coefficient “h” and generalized self-correlation coefficient “k”. The flexible change of universal logic operations is based on “h” and “k”. So Universal Logic provides a new theoretical foundation to realize more effective control for complex systems.Universal Combinatorial Operation Model is the combinatorial connective of Universal Logic. In this paper, we will only consider generalized correlation coefficient
h. So Zero-level Universal Combinatorial Operation Model is defined as follows.Definition 1 (see [13]).
Set mappingC
e
:
[
0,1
]
×
[
0,1
]
→
[
0,1
] and
(1)
C
e
(
x
,
y
,
h
)
=
ite
{
Γ
e
[
(
x
m
+
y
m
-
e
m
)
1
/
m
]
∣
x
+
y
<
2
e
;
1
-
(
Γ
1
-
e
[
(
1
-
x
)
m
+
(
1
-
y
)
m
-
(
1
-
e
)
m
]
)
1
/
m
∣
x
+
y
>
2
e
;
e
Γ
e
[
(
x
m
+
y
m
-
e
m
)
1
/
m
]
∣
x
+
y
<
2
e
;
Γ
e
[
(
x
m
+
y
m
-
e
m
)
1
/
m
]
}
.
SoC
e is Zero-level Universal Combinatorial Operation Model, denoted by C
h
e, where m
=
(
3
-
4
h
)
/
(
4
h
(
1
-
h
)
),  h
∈
[
0,1
],  m
∈
R,  e
∈
[
0,1
].Note 1.
Conditional expressionite
{
β
∣
α
;
γ
} means if α is true, then β; otherwise γ
.
ite
{
β
1
∣
α
1
;
β
2
∣
α
2
;
γ
}
=
ite
{
β
1
∣
α
1
;
ite
{
β
2
∣
α
2
;
γ
}
}. Amplitude limiting function Γ
1
[
x
]
=
ite
{
1
∣
x
>
1
;
0
∣
x
<
0
or
imaginary
number
;
x
}.Universal Combinatorial Operation Model is a cluster of combinatorial operators, which is determined by general correlation coefficient
h between propositions. In practical application, according to general correlation between propositions, we can take the corresponding one from the cluster. As generalized correlation coefficient
h is equal to some special values, the corresponding combinatorial operators are given as follows.(1)
Whenh
=
1, it means two propositions attract each other to the maximum extent. And C
e
(
x
,
y
,
1
)
=
ite
{
min
⁡
(
x
,
y
)
∣
x
+
y
<
2
e
;
max
⁡
(
x
,
y
)
∣
x
+
y
>
2
e
;
e
} is Zadeh combination C
3
e.
(2)
Whenh
=
0.75, it means two propositions are independently correlated. And C
e
(
x
,
y
,
0.75
)
=
ite
{
x
y
/
e
∣
x
+
y
<
2
e
;
(
x
+
y
-
x
y
-
e
)
/
(
1
-
e
)
∣
x
+
y
>
2
e
;
e
} is probability combination C
2
e.
(3)
Whenh
=
0.5, it means two propositions reject each other to the maximum extent or restrain each other to the minimum extent. And C
e
(
x
,
y
,
0.5
)
=
Γ
1
[
x
+
y
-
e
] is bounded combination C
1
e.
(4)
Whenh
=
0, it means two propositions restrain each other to the maximum extent. And C
e
(
x
,
y
,
0
)
=
ite
{
0
∣
x
,
y
<
e
;
1
∣
x
,
y
>
e
;
e
} is drastic combination C
0
e.


### 2.2. Universal Combinatorial Operation Model in Any Interval[
a
,
b
]

In practical control application, fuzzy domain of fuzzy variables,E and E
C, is mostly symmetrical, such as [
-
6
,
6
]. However, the conventional Universal Combinatorial Operation Model has been limited in the interval [
0
,
1
]. To this end, Chen, based on the basic idea of Universal Logic, sets up Fractal Logic in his doctoral dissertation [15], which can make inference in any interval [
a
,
b
].The combinatorial operation model inFractal Logic is described below.Definition 2 (see [15]).
Set mappingG
N
:
[
a
,
b
]
→
[
a
,
b
] and
(2)
G
N
(
x
)
=
b
+
a
-
x
.
Then G
N is normal universal Not operation model in any interval [
a
,
b
], denoted by G
N.For the above definition, normal universalNot operation model has the following characters.(1)
Closure:
(3)
G
N
(
x
)
∈
[
a
,
b
]
.
(2)
Two polar law:
(4)
G
N
(
a
)
=
b
,
G
N
(
b
)
=
a
.
(3)
Symmetric involution:
(5)
G
N
(
G
N
(
x
)
)
=
x
.Definition 3 (see [15]).
Set mappingG
C
e
:
[
a
,
b
]
×
[
a
,
b
]
→
[
a
,
b
] and
(6)
G
C
e
(
x
,
y
,
h
)
=
ite
{
-
(
e
-
a
)
m
)
×
(
(
b
-
a
)
m
)
-
1
]
1
/
m
min
⁡
(
-
(
e
-
a
)
m
)
×
(
(
b
-
a
)
m
)
-
1
]
1
/
m
e
,
(
b
-
a
)
×
[
max
⁡
(
0
,
(
(
x
-
a
)
m
+
(
y
-
a
)
m
-
(
e
-
a
)
m
(
x
-
a
)
m
)
×
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m
+
a
)
∣
x
+
y
<
2
e
;
b
+
a
-
min
⁡
(
-
(
e
-
a
)
m
)
×
(
(
b
-
a
)
m
)
-
1
]
1
/
m
e
′
,
(
b
-
a
)
×
[
(
(
b
-
a
)
m
)
-
1
max
⁡
(
0
,
(
(
b
-
x
)
m
+
(
b
-
y
)
m
-
(
b
-
e
)
m
(
b
-
x
)
m
)
×
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m
+
a
)
∣
x
+
y
>
2
e
;
e
-
(
e
-
a
)
m
)
×
(
(
b
-
a
)
m
)
-
1
]
1
/
m
}
.
ThenG
C
e is Zero-level Universal Combinatorial Operation Model in any interval [
a
,
b
], denoted by G
C
h
e, where m
=
(
3
-
4
h
)
/
(
4
h
(
1
-
h
)
), h
∈
[
0,1
], m
∈
R, e
,
e
′
∈
[
a
,
b
], and e
′
=
G
N
(
e
).


### 2.3. Demonstration of Attributes of Universal Combinatorial Operation Model in Any Interval

According to the definition ofUniversal Combinatorial Operation Model in any interval, the following characters [15] are attained.(1)
G
C
e
(
x
,
y
,
h
) conforms to the combination axiom:
(i)
Boundary condition
G
C
1:
Ifx
,
y
<
e, then G
C
e
(
x
,
y
,
h
)
≤
min
⁡
(
x
,
y
).
Ifx
,
y
>
e, then G
C
e
(
x
,
y
,
h
)
≥
max
⁡
(
x
,
y
).
Ifx
+
y
=
2
e, then G
C
e
(
x
,
y
,
h
)
=
e.
Otherwise,min
⁡
(
x
,
y
)
≤
G
C
e
(
x
,
y
,
h
)
≤
max
⁡
(
x
,
y
).
(ii)
Monotonicity
G
C
2:
G
C
e
(
x
,
y
,
h
) increases monotonously along with x and y.
(iii)
Continuity
G
C
3:
Whenh
∈
(
0,1
), G
C
e
(
x
,
y
,
h
) is continuous for all x and y.
(iv)
Commutative law
G
C
4:
(7)
G
C
e
(
x
,
y
,
h
)
=
G
C
e
(
y
,
x
,
h
)
.
(v)
Law of identical element
G
C
5:
(8)
G
C
e
(
x
,
e
,
h
)
=
x
.
(2)
Closure:
(9)
G
C
e
(
x
,
y
,
h
)
∈
[
a
,
b
]
.
(3)
Inverse law:
(10)
G
C
e
(
x
,
2
e
-
x
,
h
)
=
e
.
(4)
Renunciation law
(11)
G
C
e
(
e
,
e
,
h
)
=
e
.Theorem 4.
G
N
(
G
C
G
N
(
e
)
(
G
N
(
x
)
,
G
N
(
y
)
,
h
)
)
=
G
C
e
(
x
,
y
,
h
)
,
x
,
y
∈
[
a
,
b
]
,
e
∈
[
a
,
b
].Proof.
x
,
y
∈
[
a
,
b
]
,
e
∈
[
a
,
b
], according to the closure of normal universal Not operation and universal combinatorial operation: G
N
(
x
)
,
G
N
(
y
)
∈
[
a
,
b
],  G
N
(
e
)
∈
[
a
,
b
] and G
C
G
N
(
e
)
(
G
N
(
x
)
,
G
N
(
y
)
,
h
)
∈
[
a
,
b
] and according to the definition of universal combinatorial operation:(1)
whenx
+
y
<
2
e
(12)
G
N
(
x
)
+
G
N
(
y
)
=
(
b
+
a
-
x
)
+
(
b
+
a
-
y
)
=
2
(
b
+
a
)
-
(
x
+
y
)
>
2
(
b
+
a
-
e
)
=
2
G
N
(
e
)
.
Then, according to the definition ofG
C
e
(
x
,
y
,
h
):
(13)
G
C
G
N
(
e
)
(
G
N
(
x
)
,
G
N
(
y
)
,
h
)
=
b
+
a
-
min
⁡
(
×
(
(
b
-
a
)
m
)
-
1
)
]
1
m
a
b
+
a
-
G
N
(
e
)
,
(
b
-
a
)
×
[
(
(
b
-
a
)
m
)
-
1
max
⁡
(
0
,
(
(
b
-
G
N
(
x
)
)
m
+
(
b
-
G
N
(
y
)
)
m
-
(
b
-
G
N
(
e
)
)
m
)
×
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m
+
a
×
(
(
b
-
a
)
m
)
-
1
)
]
1
m
)
.
According to the definition of central generalized negation operation:(14)
G
N
(
x
)
=
b
+
a
-
x
,
(15)
G
N
(
y
)
=
b
+
a
-
y
,
(16)
G
N
(
e
)
=
b
+
a
-
e
.
Substituting (14), (15), and (16) separately into (13):
(17)
G
C
G
N
(
e
)
(
G
N
(
x
)
,
G
N
(
y
)
,
h
)
=
b
+
a
-
min
⁡
(
×
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m
e
,
(
b
-
a
)
×
[
×
(
(
b
-
a
)
m
)
-
1
)
max
⁡
(
0
,
(
(
x
-
a
)
m
+
(
y
-
a
)
m
-
(
e
-
a
)
m
(
x
-
a
)
m
)
×
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m
+
a
)
.
And then(18)
G
N
(
G
C
G
N
(
e
)
(
G
N
(
x
)
,
G
N
(
y
)
,
h
)
)
=
b
+
a
-
(
b
+
a
-
min
⁡
(
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m
e
,
(
b
-
a
)
×
[
max
⁡
(
0
,
(
(
x
-
a
)
m
+
(
y
-
a
)
m
-
(
e
-
a
)
m
)
×
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m

+
a
)
)
=
min
⁡
(
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m
e
,
(
b
-
a
)
×
[
max
⁡
(
0
,
(
(
x
-
a
)
m
+
(
y
-
a
)
m
-
(
e
-
a
)
m
)
×
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m
+
a
)
=
G
C
e
(
x
,
y
,
h
)
(2)
whenx
+
y
>
2
e
(19)
G
N
(
x
)
+
G
N
(
y
)
=
(
b
+
a
-
x
)
+
(
b
+
a
-
y
)
=
2
(
b
+
a
)
-
(
x
+
y
)
<
2
(
b
+
a
-
e
)
=
2
G
N
(
e
)
.
So, according to the definition ofG
C
e
(
x
,
y
,
h
):
(20)
G
C
G
N
(
e
)
(
G
N
(
x
)
,
G
N
(
y
)
,
h
)
=
min
⁡
(
×
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m
G
N
(
e
)
,
(
b
-
a
)
×
[
max
⁡
(
0
,
(
(
G
N
(
x
)
-
a
)
m
+
(
G
N
(
y
)
-
a
)
m
-
(
G
N
(
e
)
-
a
)
m
)
×
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m
+
a
)
.
Substituting (14), (15), and (16) separately into (20):
(21)
G
C
G
N
(
e
)
(
G
N
(
x
)
,
G
N
(
y
)
,
h
)
=
min
⁡
(
-
(
b
-
e
)
m
)
×
(
(
b
-
a
)
m
)
-
1
]
1
/
m
b
+
a
-
e
,
(
b
-
a
)
×
[
max
⁡
(
0
,
(
(
b
-
x
)
m
+
(
b
-
y
)
m
-
(
b
-
e
)
m
)
×
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m
+
a
)
.
And then(22)
G
N
(
G
C
G
N
(
e
)
(
G
N
(
x
)
,
G
N
(
y
)
,
h
)
)
=
b
+
a
-
min
⁡
(
-
(
b
-
e
)
m
)
×
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m
b
+
a
-
e
,
(
b
-
a
)
×
[
max
⁡
(
0
,
(
(
b
-
x
)
m
+
(
b
-
y
)
m
-
(
b
-
e
)
m
)
×
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m
+
a
)
=
G
C
e
(
x
,
y
,
h
)
(3)
whenx
+
y
=
2
e
(23)
G
N
(
x
)
+
G
N
(
y
)
=
(
b
+
a
-
x
)
+
(
b
+
a
-
y
)
=
2
(
b
+
a
)
-
(
x
+
y
)
=
2
(
b
+
a
-
e
)
=
2
G
N
(
e
)
.According to the definition of G
C
e
(
x
,
y
,
h
):
(24)
G
C
G
N
(
e
)
(
G
N
(
x
)
,
G
N
(
y
)
,
h
)
=
G
N
(
e
)
.
Then:
(25)
G
N
(
G
C
G
N
(
e
)
(
G
N
(
x
)
,
G
N
(
y
)
,
h
)
)
=
G
N
(
G
N
(
e
)
)
=
e
.
From the above, the theorem is true.Lemma 5.
G
C
G
N
(
e
)
(
G
N
(
x
)
,
G
N
(
y
)
,
h
)
=
G
N
(
G
C
e
(
x
,
y
,
h
)
).Proof.
According to Theorem4 and involution law of normal universal Not operation in any interval [
a
,
b
], the theorem can be proved simply.Lemma 6.
C
e
(
x
,
y
,
h
)
=
1
-
C
1
-
e
(
1
-
x
,
1
-
y
,
h
).Proof.
Setting the interval[
a
,
b
] of x
,
y as [
0,1
], the lemma can be proved simply.Lemma 7.
If the interval[
a
,
b
] is symmetrical about e, then G
C
e
(
x
*
,
y
*
,
h
)
=
(
G
C
e
(
x
,
y
,
h
)
)
*, where x
* and x are symmetrical about e, namely, x
*
=
2
e
-
x, y
*
,
e
*, and (
G
C
e
(
x
,
y
,
h
)
)
* are similar, e
∈
[
a
,
b
], h
∈
[
0,1
].Proof.
Since the interval[
a
,
b
] is symmetrical about e, then a
+
b
=
2
e and G
N
(
x
)
=
a
+
b
-
x
=
2
e
-
x; thus
(26)
G
N
(
x
)
=
x
*
.
Similarly(27)
G
N
(
y
)
=
y
*
,
G
N
(
e
)
=
a
+
b
-
e
=
2
e
-
e
=
e
.
From (26) and (27) the following could be obtained:
(28)
G
C
G
N
(
e
)
(
G
N
(
x
)
,
G
N
(
y
)
,
h
)
=
G
C
e
(
x
*
,
y
*
,
h
)
,
G
N
(
G
C
e
(
x
,
y
,
h
)
)
=
(
G
C
e
(
x
,
y
,
h
)
)
*
.
And then from Lemma5
(29)
G
C
e
(
x
*
,
y
*
,
h
)
=
(
G
C
e
(
x
,
y
,
h
)
)
*
.
So the theorem is true.Lemma 8.
C
0.5
(
1
-
x
,
1
-
y
,
h
)
=
1
-
C
0.5
(
x
,
y
,
h
), where h
∈
[
0,1
].Lemma 9.
If the interval[
a
,
b
] is symmetrical about the original point, then G
C
0
(
-
x
,
-
y
,
h
)
=
-
G
C
0
(
x
,
y
,
h
), h
∈
[
0,1
].This lemma indicates that, when the interval[
a
,
b
] is symmetrical about the origin and identity element e is 0, Universal Combinatorial Operation
G
C
e
(
x
,
y
,
h
) is also symmetrical about the origin.As pointed out in the literature [4], the logical relationship between input variables, error E and error variety E
C of normal two-dimensional fuzzy controller, is universal combinatorial relation in Universal Logic. Consequently, the complex reasoning process based on fuzzy rule can be replaced by the simple universal combinatorial operation, and a flexible logic control model is presented accordingly. In fuzzy control, the domains of input variables and output variable are generally symmetric to the origin, such as [
-
5
,
5
]. Obviously, it is the prerequisite of control model that the operation model be symmetric to the origin. Therefore, Lemma 9 provides a basis for the Universal Combinatorial Operation’s application in control.


### 2.4. Flexible Logic Control Method

Xiao et al. have put forward a concept of fuzzy composed variable to reduce effectively fuzzy control rules in multivariable nonlinear system [16]. According to the characteristics of controlled system and the internal relationship between input variables, the core is to construct a fuzzy composed variable by the fuzzy logic system to synthetically reflect the deviation between reference and the process output.Four output variables in single inverted pendulum are considered, which are the displacement and speed of cart,x and x
′, and the angle and angle speed between pendulum bar and vertical line, θ and θ
′. Four input variables are involved for the fuzzy controller. In the input variables of control system, the angle and angle speed, θ and θ
′, directly reflected the motion of pendulum. Therefore, according to the angle and angle speed of pendulum and the language rules shown in Table 1, a fuzzy composed variable, the error G
E
θ of pendulum, can be defined to synthetically describe the motion of pendulum with the fuzzy logic system. Similarly, a fuzzy composed variable, the error G
E
x of cart, can be defined according to the displacement and speed of cart, x and x
′, so as to synthetically describe the motion of cart. For multivariable system, it need not define, respectively, fuzzy logic system for every fuzzy composed variable. We can use a uniform fuzzy rule table, such as Table 1, and just select different quantization factors to obtain different fuzzy composed variables.Table 1
Fuzzy rules defining the fuzzy composed variableG
E
θ.
Error of pendulum rod
Angle speed of pendulum rod
NB
NM
NS
ZE
PS
PM
PB
Angle of pendulum rod
NB
NB
NB
NB
NM
NM
NS
ZE
NM
NB
NB
NM
NM
NS
ZE
PS
NS
NB
NM
NM
NS
ZE
PS
PM
ZE
NM
NM
NS
ZE
PS
PM
PM
PS
NM
NS
ZE
PS
PM
PM
PB
PM
NS
ZE
PS
PM
PM
PB
PB
PB
ZE
PS
PM
PM
PB
PB
PBNote 2.
In the paper, the fuzzy controller discussed has inputs, such as errore and error variety e
c, and output variable u is the control signal. The variables, e,  e
c, and u, are crisp values from the practical process. The fuzzy language variables, E, E
C, and U, are the corresponding fuzzy ones, and the fuzzy domains are unified as [
-
1
,
1
] with fuzzy subsets, such as negative big (NB), negative middle (NM), negative small (NS), zero (ZE), positive small (PS), positive middle (PM), and positive big (PB).Obviously, the language rules shown in Table1 actually reflect the essential relationship between error E and error variety E
C. The contents in Table 1 can be approximately divided into four parts, which separately give the language rules used to define composed variables under four conditions, such as error E and error variety E
C are both negative, error E is negative but error variety E
C is positive, error E is positive but error variety E
C is negative, and error E and error variety E
C are both positive.Both errorE and error variety E
C reflect the deviation between the reference and the output; and then a composed variable E
′ can be defined to synthetically describe the control deviation of system, based on the essential relationship between them.According to the physical meanings of errorE and error variety E
C, we can get the following conclusions.(1)
When errorE and error variety E
C are both positive, the control deviation of system is positive and tends to further increase positively. So the value of composed variable E
′ should not be smaller than their maximum value and the direction is positive. The combination rules are shown by the lower right bold corner of Table 1.
(2)
When errorE is positive but error variety E
C is negative, the control deviation of system is positive but tends to decrease. So the value of composed variable E
′ should be between E and E
C in this case. The combination rules are shown by the lower left bold corner of Table 1.
(3)
Similarly, the value of composed variableE
′ under two other conditions can also be obtained.Through above analysis, it is easy to know that the essential relationship among errorE, error variety E
C and the composed variable E
′ is just a kind of universal combinatorial one in Universal Logic. As a result, we have
(30)
E
′
=
G
C
e
(
E
,
E
C
,
h
)
.Obviously, now the identity elemente is zero.Moreover, according to the concept of negative feedback control, there is only a difference of a single sign between composed variableE
′ describing control deviation and the output variable U of controller, namely;
(31)
U
=
-
E
′
.Therefore, the relationship among errorE, error variety E
C, and the output variable U is
(32)
U
=
-
G
C
e
(
E
,
E
C
,
h
)
,
where E
,
E
C
,
U
∈
[
-
1,1
], e
=
0, and h
∈
[
0,1
]. The control method is called Flexible Logic Control Method [4].At the same time, a weighted factorα
,
α
∈
[
0,1
] is introduced to meet the requirements of different controlled objects. By adjusting the value of α, it is possible to change the weighting degrees for error E and error variety E
C. When the general correlation coefficient
h is 0.5, there is
(33)
U
=
-
G
C
e
(
α
E
,
(
1
-
α
)
E
C
,
0.5
)
=
Γ
-
1
1
[
α
E
+
(
1
-
α
)
E
C
-
e
]
=
α
E
+
(
1
-
α
)
E
C
.And the formula (33) is just the fuzzy control method proposed by Long and Wang [17]. He used a linear equation, such as (33), to describe fuzzy control rules. But the relationship among E, E
C, and U is not only linear. So (32) is a cluster of operators determined by general correlation coefficient
h, and (33) is only a special operator in the cluster as h is equal to 0.5. As a result, flexible logic control method can realize the effective control for complex system.


## 2.1. Universal Combinatorial Operation Model

In order to deal reasonably with the complex relation between factors in complex system,T-norm, S-norm or Mean operators are taken as Aggregation Operators.However,T-norm result is not bigger than the minimum value, and S-norm result is not less than the maximum value. As a result, T-norm, or S-norm can only handle mutually conflictive relation. In contrast, Mean operators can vary only between the minimum and maximum values based on its “tradeoff” concept, so it can only handle mutually consistent relation [13].Universal Logic [13], proposed by Professor He et al., is a kind of flexible logic. It considers the continuous change of not only the truth value of propositions, which is called truth value flexibility, but also the relation between propositions, which is called relation flexibility. Based on fuzzy logic, it puts forward two important coefficients: generalized correlation coefficient “h” and generalized self-correlation coefficient “k”. The flexible change of universal logic operations is based on “h” and “k”. So Universal Logic provides a new theoretical foundation to realize more effective control for complex systems.Universal Combinatorial Operation Model is the combinatorial connective of Universal Logic. In this paper, we will only consider generalized correlation coefficient
h. So Zero-level Universal Combinatorial Operation Model is defined as follows.Definition 1 (see [13]).
Set mappingC
e
:
[
0,1
]
×
[
0,1
]
→
[
0,1
] and
(1)
C
e
(
x
,
y
,
h
)
=
ite
{
Γ
e
[
(
x
m
+
y
m
-
e
m
)
1
/
m
]
∣
x
+
y
<
2
e
;
1
-
(
Γ
1
-
e
[
(
1
-
x
)
m
+
(
1
-
y
)
m
-
(
1
-
e
)
m
]
)
1
/
m
∣
x
+
y
>
2
e
;
e
Γ
e
[
(
x
m
+
y
m
-
e
m
)
1
/
m
]
∣
x
+
y
<
2
e
;
Γ
e
[
(
x
m
+
y
m
-
e
m
)
1
/
m
]
}
.
SoC
e is Zero-level Universal Combinatorial Operation Model, denoted by C
h
e, where m
=
(
3
-
4
h
)
/
(
4
h
(
1
-
h
)
),  h
∈
[
0,1
],  m
∈
R,  e
∈
[
0,1
].Note 1.
Conditional expressionite
{
β
∣
α
;
γ
} means if α is true, then β; otherwise γ
.
ite
{
β
1
∣
α
1
;
β
2
∣
α
2
;
γ
}
=
ite
{
β
1
∣
α
1
;
ite
{
β
2
∣
α
2
;
γ
}
}. Amplitude limiting function Γ
1
[
x
]
=
ite
{
1
∣
x
>
1
;
0
∣
x
<
0
or
imaginary
number
;
x
}.Universal Combinatorial Operation Model is a cluster of combinatorial operators, which is determined by general correlation coefficient
h between propositions. In practical application, according to general correlation between propositions, we can take the corresponding one from the cluster. As generalized correlation coefficient
h is equal to some special values, the corresponding combinatorial operators are given as follows.(1)
Whenh
=
1, it means two propositions attract each other to the maximum extent. And C
e
(
x
,
y
,
1
)
=
ite
{
min
⁡
(
x
,
y
)
∣
x
+
y
<
2
e
;
max
⁡
(
x
,
y
)
∣
x
+
y
>
2
e
;
e
} is Zadeh combination C
3
e.
(2)
Whenh
=
0.75, it means two propositions are independently correlated. And C
e
(
x
,
y
,
0.75
)
=
ite
{
x
y
/
e
∣
x
+
y
<
2
e
;
(
x
+
y
-
x
y
-
e
)
/
(
1
-
e
)
∣
x
+
y
>
2
e
;
e
} is probability combination C
2
e.
(3)
Whenh
=
0.5, it means two propositions reject each other to the maximum extent or restrain each other to the minimum extent. And C
e
(
x
,
y
,
0.5
)
=
Γ
1
[
x
+
y
-
e
] is bounded combination C
1
e.
(4)
Whenh
=
0, it means two propositions restrain each other to the maximum extent. And C
e
(
x
,
y
,
0
)
=
ite
{
0
∣
x
,
y
<
e
;
1
∣
x
,
y
>
e
;
e
} is drastic combination C
0
e.


## 2.2. Universal Combinatorial Operation Model in Any Interval[
a
,
b
]

In practical control application, fuzzy domain of fuzzy variables,E and E
C, is mostly symmetrical, such as [
-
6
,
6
]. However, the conventional Universal Combinatorial Operation Model has been limited in the interval [
0
,
1
]. To this end, Chen, based on the basic idea of Universal Logic, sets up Fractal Logic in his doctoral dissertation [15], which can make inference in any interval [
a
,
b
].The combinatorial operation model inFractal Logic is described below.Definition 2 (see [15]).
Set mappingG
N
:
[
a
,
b
]
→
[
a
,
b
] and
(2)
G
N
(
x
)
=
b
+
a
-
x
.
Then G
N is normal universal Not operation model in any interval [
a
,
b
], denoted by G
N.For the above definition, normal universalNot operation model has the following characters.(1)
Closure:
(3)
G
N
(
x
)
∈
[
a
,
b
]
.
(2)
Two polar law:
(4)
G
N
(
a
)
=
b
,
G
N
(
b
)
=
a
.
(3)
Symmetric involution:
(5)
G
N
(
G
N
(
x
)
)
=
x
.Definition 3 (see [15]).
Set mappingG
C
e
:
[
a
,
b
]
×
[
a
,
b
]
→
[
a
,
b
] and
(6)
G
C
e
(
x
,
y
,
h
)
=
ite
{
-
(
e
-
a
)
m
)
×
(
(
b
-
a
)
m
)
-
1
]
1
/
m
min
⁡
(
-
(
e
-
a
)
m
)
×
(
(
b
-
a
)
m
)
-
1
]
1
/
m
e
,
(
b
-
a
)
×
[
max
⁡
(
0
,
(
(
x
-
a
)
m
+
(
y
-
a
)
m
-
(
e
-
a
)
m
(
x
-
a
)
m
)
×
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m
+
a
)
∣
x
+
y
<
2
e
;
b
+
a
-
min
⁡
(
-
(
e
-
a
)
m
)
×
(
(
b
-
a
)
m
)
-
1
]
1
/
m
e
′
,
(
b
-
a
)
×
[
(
(
b
-
a
)
m
)
-
1
max
⁡
(
0
,
(
(
b
-
x
)
m
+
(
b
-
y
)
m
-
(
b
-
e
)
m
(
b
-
x
)
m
)
×
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m
+
a
)
∣
x
+
y
>
2
e
;
e
-
(
e
-
a
)
m
)
×
(
(
b
-
a
)
m
)
-
1
]
1
/
m
}
.
ThenG
C
e is Zero-level Universal Combinatorial Operation Model in any interval [
a
,
b
], denoted by G
C
h
e, where m
=
(
3
-
4
h
)
/
(
4
h
(
1
-
h
)
), h
∈
[
0,1
], m
∈
R, e
,
e
′
∈
[
a
,
b
], and e
′
=
G
N
(
e
).


## 2.3. Demonstration of Attributes of Universal Combinatorial Operation Model in Any Interval

According to the definition ofUniversal Combinatorial Operation Model in any interval, the following characters [15] are attained.(1)
G
C
e
(
x
,
y
,
h
) conforms to the combination axiom:
(i)
Boundary condition
G
C
1:
Ifx
,
y
<
e, then G
C
e
(
x
,
y
,
h
)
≤
min
⁡
(
x
,
y
).
Ifx
,
y
>
e, then G
C
e
(
x
,
y
,
h
)
≥
max
⁡
(
x
,
y
).
Ifx
+
y
=
2
e, then G
C
e
(
x
,
y
,
h
)
=
e.
Otherwise,min
⁡
(
x
,
y
)
≤
G
C
e
(
x
,
y
,
h
)
≤
max
⁡
(
x
,
y
).
(ii)
Monotonicity
G
C
2:
G
C
e
(
x
,
y
,
h
) increases monotonously along with x and y.
(iii)
Continuity
G
C
3:
Whenh
∈
(
0,1
), G
C
e
(
x
,
y
,
h
) is continuous for all x and y.
(iv)
Commutative law
G
C
4:
(7)
G
C
e
(
x
,
y
,
h
)
=
G
C
e
(
y
,
x
,
h
)
.
(v)
Law of identical element
G
C
5:
(8)
G
C
e
(
x
,
e
,
h
)
=
x
.
(2)
Closure:
(9)
G
C
e
(
x
,
y
,
h
)
∈
[
a
,
b
]
.
(3)
Inverse law:
(10)
G
C
e
(
x
,
2
e
-
x
,
h
)
=
e
.
(4)
Renunciation law
(11)
G
C
e
(
e
,
e
,
h
)
=
e
.Theorem 4.
G
N
(
G
C
G
N
(
e
)
(
G
N
(
x
)
,
G
N
(
y
)
,
h
)
)
=
G
C
e
(
x
,
y
,
h
)
,
x
,
y
∈
[
a
,
b
]
,
e
∈
[
a
,
b
].Proof.
x
,
y
∈
[
a
,
b
]
,
e
∈
[
a
,
b
], according to the closure of normal universal Not operation and universal combinatorial operation: G
N
(
x
)
,
G
N
(
y
)
∈
[
a
,
b
],  G
N
(
e
)
∈
[
a
,
b
] and G
C
G
N
(
e
)
(
G
N
(
x
)
,
G
N
(
y
)
,
h
)
∈
[
a
,
b
] and according to the definition of universal combinatorial operation:(1)
whenx
+
y
<
2
e
(12)
G
N
(
x
)
+
G
N
(
y
)
=
(
b
+
a
-
x
)
+
(
b
+
a
-
y
)
=
2
(
b
+
a
)
-
(
x
+
y
)
>
2
(
b
+
a
-
e
)
=
2
G
N
(
e
)
.
Then, according to the definition ofG
C
e
(
x
,
y
,
h
):
(13)
G
C
G
N
(
e
)
(
G
N
(
x
)
,
G
N
(
y
)
,
h
)
=
b
+
a
-
min
⁡
(
×
(
(
b
-
a
)
m
)
-
1
)
]
1
m
a
b
+
a
-
G
N
(
e
)
,
(
b
-
a
)
×
[
(
(
b
-
a
)
m
)
-
1
max
⁡
(
0
,
(
(
b
-
G
N
(
x
)
)
m
+
(
b
-
G
N
(
y
)
)
m
-
(
b
-
G
N
(
e
)
)
m
)
×
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m
+
a
×
(
(
b
-
a
)
m
)
-
1
)
]
1
m
)
.
According to the definition of central generalized negation operation:(14)
G
N
(
x
)
=
b
+
a
-
x
,
(15)
G
N
(
y
)
=
b
+
a
-
y
,
(16)
G
N
(
e
)
=
b
+
a
-
e
.
Substituting (14), (15), and (16) separately into (13):
(17)
G
C
G
N
(
e
)
(
G
N
(
x
)
,
G
N
(
y
)
,
h
)
=
b
+
a
-
min
⁡
(
×
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m
e
,
(
b
-
a
)
×
[
×
(
(
b
-
a
)
m
)
-
1
)
max
⁡
(
0
,
(
(
x
-
a
)
m
+
(
y
-
a
)
m
-
(
e
-
a
)
m
(
x
-
a
)
m
)
×
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m
+
a
)
.
And then(18)
G
N
(
G
C
G
N
(
e
)
(
G
N
(
x
)
,
G
N
(
y
)
,
h
)
)
=
b
+
a
-
(
b
+
a
-
min
⁡
(
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m
e
,
(
b
-
a
)
×
[
max
⁡
(
0
,
(
(
x
-
a
)
m
+
(
y
-
a
)
m
-
(
e
-
a
)
m
)
×
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m

+
a
)
)
=
min
⁡
(
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m
e
,
(
b
-
a
)
×
[
max
⁡
(
0
,
(
(
x
-
a
)
m
+
(
y
-
a
)
m
-
(
e
-
a
)
m
)
×
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m
+
a
)
=
G
C
e
(
x
,
y
,
h
)
(2)
whenx
+
y
>
2
e
(19)
G
N
(
x
)
+
G
N
(
y
)
=
(
b
+
a
-
x
)
+
(
b
+
a
-
y
)
=
2
(
b
+
a
)
-
(
x
+
y
)
<
2
(
b
+
a
-
e
)
=
2
G
N
(
e
)
.
So, according to the definition ofG
C
e
(
x
,
y
,
h
):
(20)
G
C
G
N
(
e
)
(
G
N
(
x
)
,
G
N
(
y
)
,
h
)
=
min
⁡
(
×
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m
G
N
(
e
)
,
(
b
-
a
)
×
[
max
⁡
(
0
,
(
(
G
N
(
x
)
-
a
)
m
+
(
G
N
(
y
)
-
a
)
m
-
(
G
N
(
e
)
-
a
)
m
)
×
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m
+
a
)
.
Substituting (14), (15), and (16) separately into (20):
(21)
G
C
G
N
(
e
)
(
G
N
(
x
)
,
G
N
(
y
)
,
h
)
=
min
⁡
(
-
(
b
-
e
)
m
)
×
(
(
b
-
a
)
m
)
-
1
]
1
/
m
b
+
a
-
e
,
(
b
-
a
)
×
[
max
⁡
(
0
,
(
(
b
-
x
)
m
+
(
b
-
y
)
m
-
(
b
-
e
)
m
)
×
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m
+
a
)
.
And then(22)
G
N
(
G
C
G
N
(
e
)
(
G
N
(
x
)
,
G
N
(
y
)
,
h
)
)
=
b
+
a
-
min
⁡
(
-
(
b
-
e
)
m
)
×
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m
b
+
a
-
e
,
(
b
-
a
)
×
[
max
⁡
(
0
,
(
(
b
-
x
)
m
+
(
b
-
y
)
m
-
(
b
-
e
)
m
)
×
(
(
b
-
a
)
m
)
-
1
)
]
1
/
m
+
a
)
=
G
C
e
(
x
,
y
,
h
)
(3)
whenx
+
y
=
2
e
(23)
G
N
(
x
)
+
G
N
(
y
)
=
(
b
+
a
-
x
)
+
(
b
+
a
-
y
)
=
2
(
b
+
a
)
-
(
x
+
y
)
=
2
(
b
+
a
-
e
)
=
2
G
N
(
e
)
.According to the definition of G
C
e
(
x
,
y
,
h
):
(24)
G
C
G
N
(
e
)
(
G
N
(
x
)
,
G
N
(
y
)
,
h
)
=
G
N
(
e
)
.
Then:
(25)
G
N
(
G
C
G
N
(
e
)
(
G
N
(
x
)
,
G
N
(
y
)
,
h
)
)
=
G
N
(
G
N
(
e
)
)
=
e
.
From the above, the theorem is true.Lemma 5.
G
C
G
N
(
e
)
(
G
N
(
x
)
,
G
N
(
y
)
,
h
)
=
G
N
(
G
C
e
(
x
,
y
,
h
)
).Proof.
According to Theorem4 and involution law of normal universal Not operation in any interval [
a
,
b
], the theorem can be proved simply.Lemma 6.
C
e
(
x
,
y
,
h
)
=
1
-
C
1
-
e
(
1
-
x
,
1
-
y
,
h
).Proof.
Setting the interval[
a
,
b
] of x
,
y as [
0,1
], the lemma can be proved simply.Lemma 7.
If the interval[
a
,
b
] is symmetrical about e, then G
C
e
(
x
*
,
y
*
,
h
)
=
(
G
C
e
(
x
,
y
,
h
)
)
*, where x
* and x are symmetrical about e, namely, x
*
=
2
e
-
x, y
*
,
e
*, and (
G
C
e
(
x
,
y
,
h
)
)
* are similar, e
∈
[
a
,
b
], h
∈
[
0,1
].Proof.
Since the interval[
a
,
b
] is symmetrical about e, then a
+
b
=
2
e and G
N
(
x
)
=
a
+
b
-
x
=
2
e
-
x; thus
(26)
G
N
(
x
)
=
x
*
.
Similarly(27)
G
N
(
y
)
=
y
*
,
G
N
(
e
)
=
a
+
b
-
e
=
2
e
-
e
=
e
.
From (26) and (27) the following could be obtained:
(28)
G
C
G
N
(
e
)
(
G
N
(
x
)
,
G
N
(
y
)
,
h
)
=
G
C
e
(
x
*
,
y
*
,
h
)
,
G
N
(
G
C
e
(
x
,
y
,
h
)
)
=
(
G
C
e
(
x
,
y
,
h
)
)
*
.
And then from Lemma5
(29)
G
C
e
(
x
*
,
y
*
,
h
)
=
(
G
C
e
(
x
,
y
,
h
)
)
*
.
So the theorem is true.Lemma 8.
C
0.5
(
1
-
x
,
1
-
y
,
h
)
=
1
-
C
0.5
(
x
,
y
,
h
), where h
∈
[
0,1
].Lemma 9.
If the interval[
a
,
b
] is symmetrical about the original point, then G
C
0
(
-
x
,
-
y
,
h
)
=
-
G
C
0
(
x
,
y
,
h
), h
∈
[
0,1
].This lemma indicates that, when the interval[
a
,
b
] is symmetrical about the origin and identity element e is 0, Universal Combinatorial Operation
G
C
e
(
x
,
y
,
h
) is also symmetrical about the origin.As pointed out in the literature [4], the logical relationship between input variables, error E and error variety E
C of normal two-dimensional fuzzy controller, is universal combinatorial relation in Universal Logic. Consequently, the complex reasoning process based on fuzzy rule can be replaced by the simple universal combinatorial operation, and a flexible logic control model is presented accordingly. In fuzzy control, the domains of input variables and output variable are generally symmetric to the origin, such as [
-
5
,
5
]. Obviously, it is the prerequisite of control model that the operation model be symmetric to the origin. Therefore, Lemma 9 provides a basis for the Universal Combinatorial Operation’s application in control.


## 2.4. Flexible Logic Control Method

Xiao et al. have put forward a concept of fuzzy composed variable to reduce effectively fuzzy control rules in multivariable nonlinear system [16]. According to the characteristics of controlled system and the internal relationship between input variables, the core is to construct a fuzzy composed variable by the fuzzy logic system to synthetically reflect the deviation between reference and the process output.Four output variables in single inverted pendulum are considered, which are the displacement and speed of cart,x and x
′, and the angle and angle speed between pendulum bar and vertical line, θ and θ
′. Four input variables are involved for the fuzzy controller. In the input variables of control system, the angle and angle speed, θ and θ
′, directly reflected the motion of pendulum. Therefore, according to the angle and angle speed of pendulum and the language rules shown in Table 1, a fuzzy composed variable, the error G
E
θ of pendulum, can be defined to synthetically describe the motion of pendulum with the fuzzy logic system. Similarly, a fuzzy composed variable, the error G
E
x of cart, can be defined according to the displacement and speed of cart, x and x
′, so as to synthetically describe the motion of cart. For multivariable system, it need not define, respectively, fuzzy logic system for every fuzzy composed variable. We can use a uniform fuzzy rule table, such as Table 1, and just select different quantization factors to obtain different fuzzy composed variables.Table 1
Fuzzy rules defining the fuzzy composed variableG
E
θ.
Error of pendulum rod
Angle speed of pendulum rod
NB
NM
NS
ZE
PS
PM
PB
Angle of pendulum rod
NB
NB
NB
NB
NM
NM
NS
ZE
NM
NB
NB
NM
NM
NS
ZE
PS
NS
NB
NM
NM
NS
ZE
PS
PM
ZE
NM
NM
NS
ZE
PS
PM
PM
PS
NM
NS
ZE
PS
PM
PM
PB
PM
NS
ZE
PS
PM
PM
PB
PB
PB
ZE
PS
PM
PM
PB
PB
PBNote 2.
In the paper, the fuzzy controller discussed has inputs, such as errore and error variety e
c, and output variable u is the control signal. The variables, e,  e
c, and u, are crisp values from the practical process. The fuzzy language variables, E, E
C, and U, are the corresponding fuzzy ones, and the fuzzy domains are unified as [
-
1
,
1
] with fuzzy subsets, such as negative big (NB), negative middle (NM), negative small (NS), zero (ZE), positive small (PS), positive middle (PM), and positive big (PB).Obviously, the language rules shown in Table1 actually reflect the essential relationship between error E and error variety E
C. The contents in Table 1 can be approximately divided into four parts, which separately give the language rules used to define composed variables under four conditions, such as error E and error variety E
C are both negative, error E is negative but error variety E
C is positive, error E is positive but error variety E
C is negative, and error E and error variety E
C are both positive.Both errorE and error variety E
C reflect the deviation between the reference and the output; and then a composed variable E
′ can be defined to synthetically describe the control deviation of system, based on the essential relationship between them.According to the physical meanings of errorE and error variety E
C, we can get the following conclusions.(1)
When errorE and error variety E
C are both positive, the control deviation of system is positive and tends to further increase positively. So the value of composed variable E
′ should not be smaller than their maximum value and the direction is positive. The combination rules are shown by the lower right bold corner of Table 1.
(2)
When errorE is positive but error variety E
C is negative, the control deviation of system is positive but tends to decrease. So the value of composed variable E
′ should be between E and E
C in this case. The combination rules are shown by the lower left bold corner of Table 1.
(3)
Similarly, the value of composed variableE
′ under two other conditions can also be obtained.Through above analysis, it is easy to know that the essential relationship among errorE, error variety E
C and the composed variable E
′ is just a kind of universal combinatorial one in Universal Logic. As a result, we have
(30)
E
′
=
G
C
e
(
E
,
E
C
,
h
)
.Obviously, now the identity elemente is zero.Moreover, according to the concept of negative feedback control, there is only a difference of a single sign between composed variableE
′ describing control deviation and the output variable U of controller, namely;
(31)
U
=
-
E
′
.Therefore, the relationship among errorE, error variety E
C, and the output variable U is
(32)
U
=
-
G
C
e
(
E
,
E
C
,
h
)
,
where E
,
E
C
,
U
∈
[
-
1,1
], e
=
0, and h
∈
[
0,1
]. The control method is called Flexible Logic Control Method [4].At the same time, a weighted factorα
,
α
∈
[
0,1
] is introduced to meet the requirements of different controlled objects. By adjusting the value of α, it is possible to change the weighting degrees for error E and error variety E
C. When the general correlation coefficient
h is 0.5, there is
(33)
U
=
-
G
C
e
(
α
E
,
(
1
-
α
)
E
C
,
0.5
)
=
Γ
-
1
1
[
α
E
+
(
1
-
α
)
E
C
-
e
]
=
α
E
+
(
1
-
α
)
E
C
.And the formula (33) is just the fuzzy control method proposed by Long and Wang [17]. He used a linear equation, such as (33), to describe fuzzy control rules. But the relationship among E, E
C, and U is not only linear. So (32) is a cluster of operators determined by general correlation coefficient
h, and (33) is only a special operator in the cluster as h is equal to 0.5. As a result, flexible logic control method can realize the effective control for complex system.


## 3. Neuron Model of Zero-Level Universal Combinatorial Operation

The uniform neuron model of generalized logic operators in any intervals [a
,
b] is established in the literature [18]. Seven neuron models of logic operation are given separately, such as Not, And, Or, Implication, Equation, Average, and Combination.But the neuron model of combination logic operation is too complicated for practical applications.In this section, a new neuron model of zero-level universal combinatorial operation is constructed based on Lemma6 in Section 2.3.It is learnt from Lemma6 in Section 2.3 that
(34)
C
e
(
x
,
y
,
h
)
=
1
-
C
1
-
e
(
1
-
x
,
1
-
y
,
h
)
.Therefore, the model of zero-level universal combinatorial operation can be represented by the neuron shown in Figure1.Figure 1
The neuron model of zero-level universal combinatorial operation.The neuron model is composed of several subneurons which are interconnected. In the model, there are three input parameters:x, y, and e. And the output z, net denotes the weighted sum of all inputs, and all unmarked weights are 1. The transfer functions of the subneurons marked ④, ⑤, ⑥ are f
(
x
)
=
x
m with the parameter m and m
=
(
3
-
4
h
)
/
(
4
h
(
1
-
h
)
). And the transfer function of the sub-neuron marked ⑦ is f
-
1
(
x
). The sub-neurons marked ①, ②, ③, ⑧, and ⑨ have the following transfer functions, respectively:
(35)
f
1
(
u
,
v
,
e
)
=
{
u
u
+
v
<
2
e
1
-
u
u
+
v
>
2
e
e
otherwise
,
f
2
(
u
,
v
,
e
)
=
{
v
u
+
v
<
2
e
1
-
v
u
+
v
>
2
e
e
otherwise
,
f
3
(
u
,
v
,
e
)
=
{
e
u
+
v
<
2
e
1
-
e
u
+
v
>
2
e
e
otherwise
,
f
4
(
u
,
v
)
=
Γ
v
[
u
]
,
f
5
(
u
,
v
,
w
,
e
)
=
{
w
u
+
v
<
2
e
1
-
w
u
+
v
>
2
e
w
otherwise
.We will discuss if the neuron model shown in Figure1 has realized zero-level universal combinatorial operation.(
1)
x
+
y
<
2
e. According to the definition of transfer functions f
1, f
2, and f
3, it is easy to find that the outputs of the sub-neurons marked ①, ②, ③ separately are x, y, and e. According to the definition of transfer function f, it is found that the outputs of the sub-neuron marked ④, ⑤, ⑥ separately are x
m, y
m, and e
m. Therefore, the output of the sub-neuron marked ⑦ is (
x
m
+
y
m
-
e
m
)
1
/
m. According to the definition of transfer function f
4, the output of the sub-neuron marked ⑧ is Γ
e
[
(
x
m
+
y
m
-
e
m
)
1
/
m
]. According to the definition of transfer function f
5, the output of the sub-neuron marked ⑨ is Γ
e
[
(
x
m
+
y
m
-
e
m
)
1
/
m
]. When x
+
y
<
2
e, the output of the neuron is
(36)
Γ
e
[
(
x
m
+
y
m
-
e
m
)
1
/
m
]
.(
2)
x
+
y
>
2
e. It is similarly learnt that the output of the sub-neuron marked ⑦ is (
(
1
-
x
)
m
+
(
1
-
y
)
m
-
(
1
-
e
)
m
)
1
/
m, the output of the sub-neuron marked ⑧ is Γ
1
-
e
[
(
(
1
-
x
)
m
+
(
1
-
y
)
m
-
(
1
-
e
)
m
)
1
/
m
], and the output of the sub-neuron marked ⑨ is 1
-
Γ
1
-
e
[
(
(
1
-
x
)
m
+
(
1
-
y
)
m
-
(
1
-
e
)
m
)
1
/
m
]. When x
+
y
>
2
e, the output of the neuron is
(37)
1
-
Γ
1
-
e
[
(
(
1
-
x
)
m
+
(
1
-
y
)
m
-
(
1
-
e
)
m
)
1
/
m
]
.(
3)
x
+
y
=
2
e. It is similarly learnt that the output of the sub-neuron marked ⑦ is e, the output of the sub-neuron marked ⑧ is e, and the output of the sub-neuron marked ⑨ is e. When x
+
y
=
2
e, the output of neuron is e.In summary, the neuron model shown in Figure1 has fully realized the model of zero-level universal combinatorial operation.In practical control application, the fuzzy domains of system variables are generally expressed as[
-
n
,
n
], for example, [
-
6
,
6
]. To simplify the neuron model, we can unify the fuzzy domains of system variables as [
-
1
,
1
]. Obviously, the identity element e is 0. Therefore, the model of zero-level universal combinatorial operation in the interval [
-
1
,
1
] can be obtained by the definition of universal combinatorial operation model in any intervals [
a
,
b
] in Section 2.2,
(38)
G
C
0
(
x
,
y
,
h
)
=
ite
{
min
⁡
(
0,2
[
max
⁡
(
0
,
(
x
+
1
)
m
+
(
y
+
1
)
m
-
1
2
m
)
]
1
/
m
-
1
[
max
⁡
(
0
,
(
x
+
1
)
m
+
(
y
+
1
)
m
-
1
2
m
)
]
1
/
m
)
∣
x
+
y
<
0
;
-
min
⁡
(
0,2
[
max
⁡
(
0
,
(
1
-
x
)
m
+
(
1
-
y
)
m
-
1
2
m
)
]
1
/
m
-
1
[
max
⁡
(
0
,
(
1
-
x
)
m
+
(
1
-
y
)
m
-
1
2
m
)
]
1
/
m
)
∣
x
+
y
>
0
;
0
}
.Therefore, it is easy to obtain the neuron model of zero-level universal combinatorial operation in the interval[
-
1
,
1
] according to the neuron model given in Figure 1, as shown in Figure 2.Figure 2
The neuron model of zero-level universal combinatorial operation in the interval [-
1
,
1].The neuron model is composed of interconnected multiple sub-neurons, and its internal structure is a small artificial neural network. In the model,x and y are the inputs, z is the output, net
1, net
2, net
3, and net
4 are the weighted sum of all inputs, the values of connecting weights without specific marking are 1, and ①
,
②
,
…
,
⑧ are marks of some sub-neurons for convenient discussion. The sub-neurons marked ③ and ④ have the transfer function f
(
x
)
=
x
m with the parameter m. And the sub-neuron marked ⑤ also has the transfer function f
3
(
x
)
=
max
⁡
(
0
,
x
/
(
2
m
)
) with the parameter m, where m
=
(
3
-
4
h
)
/
(
4
h
(
1
-
h
)
). The sub-neuron marked ⑥ also has the transfer function f
-
1
(
x
) with the parameter m, and the sub-neurons marked ①, ②, ⑦, and ⑧ have the following transfer functions, respectively:
(39)
f
1
(
u
,
v
)
=
{
u
u
+
v
<
0
-
u
u
+
v
>
0
0
otherwise
,
f
2
(
u
,
v
)
=
{
v
u
+
v
<
0
-
v
u
+
v
>
0
0
otherwise
,
f
4
(
u
)
=
min
⁡
(
0
,
u
)
,
f
5
(
u
,
v
,
w
)
=
{
w
u
+
v
<
0
-
w
u
+
v
>
0
w
otherwise
.


## 4. Realization of Artificial Neural Network of Flexible Logic Control Model

In the practical complex system, the control objects are always not less than one, so multiple system variables are selected for feedback control. The basic idea of flexible logic control model presented in Section2.4 is to design sub-goal controller based on the control object, and each subgoal controller is designed with the flexible logic control method. Then the output variable of control system is the weighted sum of the outputs of these sub-goal controllers.Obviously, most of the operators involved in this model are universal combinatorial operation. Therefore, it is easy to implement the artificial neural network of the flexible logic control model by using the neuron model of zero-level universal combinatorial operation in Section3. In this section, the artificial neural network structure of flexible logic control model is presented using the double inverted-pendulum system as the control object.The objective is to maintain the rods in an upright position and the cart in an appointed position in the rail. There are six output variables and one input variable in double inverted-pendulum, which arex, x
′, θ
1, θ
1
′, θ
2, θ
2
′, and u. The variables, E
x, E
C
x, E
θ
1, E
C
θ
1, E
θ
2, E
C
θ
2, and U, are the corresponding fuzzy ones, and the fuzzy domains are unified as [
-
1
,
1
].We can design three subcontrollers with the flexible logic control method. One is to maintain the cart in an appointed position with two input variables,E
x and E
C
x. The other ones are to maintain, respectively, the rods in an upright position with two input variables E
θ
i and E
C
θ
i
(
i
=
1
,
2
). And we lead into weighted factors, such as α
x, α
θ
1, and α
θ
2. The three sub-controllers are designed as follows:
(40)
U
x
=
-
G
C
0
(
α
x
E
x
,
(
1
-
α
x
)
E
C
x
,
h
x
)
,
(41)
α
x
=
(
α
s
_
x
-
α
0
_
x
)
|
E
x
|
+
α
0
_
x
,
(42)
U
θ
i
=
-
G
C
0
(
α
θ
i
E
θ
i
,
(
1
-
α
θ
i
)
E
C
θ
i
,
h
θ
i
)
,
(43)
α
θ
i
=
(
α
s
_
θ
i
-
α
0
_
θ
i
)
|
E
θ
i
|
+
α
0
_
θ
i
.When the control signalU is combined, we lead into three weighted factors, such as k
x, k
θ
1, and k
θ
2, for the three sub-controllers. According to (40) and (42), we can get the output of the controller as follows:
(44)
U
=
k
x
U
x
+
k
θ
1
U
θ
1
+
k
θ
2
U
θ
2
=
k
x
(
-
G
C
0
(
α
x
E
x
,
(
1
-
α
x
)
E
C
x
,
h
x
)
)
+
∑
i
=
1
2
k
θ
i
(
-
G
C
0
(
α
θ
i
E
θ
i
,
(
1
-
α
θ
i
)
E
C
θ
i
,
h
θ
i
)
)
,
where E
x, E
C
x, E
θ
i, E
C
θ
i, and U
∈
[
-
1,1
], G
C
0
(
x
,
y
,
h
) is (38), h
x, h
θ
1, and h
θ
2 are general correlation coefficients, h
x, h
θ
1, and h
θ
2
∈
[
0,1
], α
x, α
θ
1, and α
θ
2
∈
[
0,1
], k
x, k
θ
1, and k
θ
2
∈
[
-
1,1
], 0
≤
α
0
_
x
≤
α
s
_
x
≤
1, 0
≤
α
0
_
θ
1
≤
α
s
_
θ
1
≤
1, and 0
≤
α
0
_
θ
2
≤
α
s
_
θ
2
≤
1.By the above analysis, we can obtain the control model of double inverted-pendulum, as shown in Figure3. Based on the neuron model of zero-level universal combinatorial operation, the artificial neural network structure of the control model could be obtained, as shown in Figure 4.Figure 3
The flexible logic control model of double inverted-pendulum.Figure 4
The neural network structure of flexible logic control model.The displacement and speed of cart,x and x
′, the angle and angle speed between the lower pendulum bar and vertical line, θ
1 and θ
1
′, and the angle and angle speed between the upper pendulum bar and vertical line, θ
2 and θ
2
′, are the input variables of the artificial neural network. The output variable is the control signal. ①, ②,…, ⑥ are marks of some sub-neurons for convenient discussion. Net is the weighted sum of these inputs of the corresponding sub-neurons. The connection weight values of the unmarked sub-neurons are 1. The sub-neurons marked ①, ②, and ③ are the neuron model of zero-level universal combinatorial operation in the interval [
-
1
,
1
], as given in Section 3. The other sub-neurons have the following transfer functions, respectively:
(45)
f
1
(
u
)
=
(
α
s
-
α
0
)
|
u
|
+
α
0
,
f
2
(
u
)
=
1
-
u
,
f
3
(
u
,
v
)
=
u
v
,
where the transfer function f
1 has the parameters such as α
0 and α
s, and the transfer functions of the sub-neurons marked ④, ⑤, and ⑥ separately have the parameters α
0
_
x,  α
s
_
x,  α
0
_
θ
1,  α
s
_
θ
1,  α
0
_
θ
2, and α
s
_
θ
2.From the Figure4, we can see that the artificial neural network structure of the control model is a feed forward neural network. Genetic algorithms (GAs) are a robust and efficient optimization technique based on the mechanism of natural selection and natural genetics [19]. One of the important features of GAs is that they are a population-based search technique. Instead of moving from one single point to another like traditional mathematical programming techniques, GAs always maintain and manipulate a solution set. Therefore, GAs are used to train the artificial neural network.


## 5. Results of Experiments

This section takes a double inverted-pendulum physical system, for example, to show the feasibility and validity of the flexible logic control method based on artificial neural network.Owing to the rapidity and the absolute instability of the inverted-pendulum system, it requires a high real-time processing frequency. Therefore, the sampling interval of this system is set as 5 ms. When the system is running for 20 s, the control effect comparison could be conducted. The physical parameters of the system are given in Table2.Table 2
The physical parameters of the double inverted-pendulum.
Symbol
Value
m
0
0.924 kg
m
1
0.185 kg
m
2
0.2 kg
l
1
0.283 m
l
2
0.245 m
L
1
0.483 m
f
0
0.1 N·s/m
J
1
0.00547 kg·m2
J
2
0.00549 kg·m2The flexible logic control method based on artificial neural network is applied into the above double inverted-pendulum physical system. The network parameters are optimized by genetic algorithms. The definition of the fitness function is shown as (46) and the control parameters of the system are shown in Table 3. Consider
(46)
Dis
=
∑
i
=
0
N
(
(
(
x
2
(
i
)
/
25
)
+
(
x
′
2
(
i
)
/
50
)
+
(
θ
1
2
(
i
)
/
5
)
+
(
θ
1
′
2
(
i
)
/
10
)
+
θ
2
2
(
i
)
+
θ
2
′
2
(
i
)
)
×
(
N
)
-
1
)
,
fitness
=
1
(
Dis
/
10
)
.Table 3
The control parameters of the control model.
Symbol
Value
Description
K
e
_
x
29.4118
Quantification factor forE
x
K
c
_
x
16.3529
Quantification factor forE
C
x
h
x
0.2902
General correlation coefficient betweenE
x and E
C
x
a
0
_
x
0.2157
Minimum value ofa
x
a
s
_
x
0.5258
Maximum value ofa
x
K
e
_
θ
1
73.3333
Quantification factor forE
θ
1
K
c
_
θ
1
8.0000
Quantification factor forE
C
θ
1
h
θ
1
0.8275
General correlation coefficient betweenE
θ
1 and E
C
θ
1
a
0
_
θ
1
0.7569
Minimum value ofa
θ
1
a
s
_
θ
1
0.7701
Maximum value ofa
θ
1
K
e
_
θ
2
47.4510
Quantification factor forE
θ
1
K
c
_
θ
2
14.9412
Quantification factor forE
C
θ
1
h
θ
2
0.5674
General correlation coefficient betweenE
θ
2 and E
C
θ
2
a
0
_
θ
2
0.6745
Minimum value ofa
θ
2
a
s
_
θ
2
0.7530
Maximum value ofa
θ
2
K
x
0.4902
Weighted factor of the subcontroller for cart
K
θ
1
−0.3645
Weighted factor of the subcontroller for rod1
K
θ
2
0.9686
Weighted factor of the subcontroller for rod2
K
u
8.7843
Proportion factor forUThree experiments as the stability control, anti-interference control, and the free movement in a level track of cart pendulum system for the above double inverted-pendulum physical system are carried out by using the controller made of these control parameters.


### 5.1. Stability Control Experiment

Under the same initial state such asx
(
0
)
=
0 m, θ
1
(
0
)
=
0.05 rad and θ
2
(
0
)
=
0.05 rad, the stability control of the double inverted-pendulum physical system has been realized. The experimental results are shown in Figures 5, 6, 7, and 8 (running time: 20 s).The car displacement and speed curve when the system is steady.
(a)
The cart displacement curve
(b)
The cart speed curveThe lower rod angle curve and angular velocity curve when the system is steady.
(a)
The lower rod angle curve
(b)
The lower rod angular velocity curveThe upper rod angle curve and angular velocity curve when the system is steady.
(a)
The upper rod angle curve
(b)
The upper rod angular velocity curveFigure 8
The controlled quantity curve when the system is steady.From the experimental results, we can find when the system is in a stable state, the deviation of the cart displacement could be about 0.01 m with a good stability control effect.


### 5.2. Anti-Interference Control Experiment

After knocking the upper pendulum twice when the system reaches a steady state, it shows strong anti-interference ability. The experimental results are shown in Figures9, 10, 11, and 12.The car displacement and speed curve when the system is interfered.
(a)
The cart displacement curve
(b)
The cart speed curveThe lower rod angle curve and angular velocity curve when the system is interfered.
(a)
The lower rod angle curve
(b)
The lower rod angular velocity curveThe upper rod angle curve and angular velocity curve when the system is interfered.
(a)
The upper rod angle curve
(b)
The upper rod angular velocity curveFigure 12
The controlled quantity curve when the system is interfered.From the experimental results, we can find the system becomes stable again after 1.9 s.


### 5.3. Free Movement Control Experiment of Cart Pendulum System

When the system reaches a steady state, the target positionx
d is changed on-line, such as from 0 m to −0.2 m. So the free movement control of cart pendulum on a level track could be realized. The experimental results are shown in Figures 13, 14, 15, and 16.The car displacement and speed curve when the cart is moving in a level track.
(a)
The cart displacement curve
(b)
The cart speed curveThe lower rod angle curve and angular velocity curve when the cart is moving in a level track.
(a)
The lower rod angle curve
(b)
The lower rod angular velocity curveThe upper rod angle curve and angular velocity curve when the cart is moving in a level track.
(a)
The upper rod angle curve
(b)
The upper rod angular velocity curveFigure 16
The controlled quantity curve when the cart is moving in a level track.From the experimental results, we can see the cart moves to the new target location after 2.7 s and the system remains stable.


### 5.4. Self-Adaptive Experiment of Control Model

Some physical parameters have been changed in another double inverted-pendulum physical system. The physical parameters of the system are shown in Table4. The control model also can realize the stability control of the new double inverted-pendulum after changing some parameters. The experimental results are shown in Figures 17, 18, 19, and 20. The experimental results show that the control model based on artificial neural network has excellent self-adaptability and portability.Table 4
The physical parameters of another double inverted-pendulum.
Symbol
Value
m
0
0.595 kg
m
1
0.158 kg
m
2
0.142 kg
l
1
0.1237 m
l
2
0.2265 m
L
1
0.151 m
f
0
0.1 N·s/m
J
1
0.0028394 kg·m2
J
2
0.0024329 kg·m2The car displacement and speed curve when the system is steady.
(a)
The cart displacement curve
(b)
The cart speed curveThe lower rod angle curve and angular velocity curve when the system is steady.
(a)
The lower rod angle curve
(b)
The lower rod angular velocity curveThe upper rod angle curve and angular velocity curve when the system is steady.
(a)
The upper rod angle curve
(b)
The upper rod angular velocity curveFigure 20
The controlled quantity curve when the system is steady.


### 5.5. Comparison of the Experimental Results

Cheng et al. [20] put forward a parameter fuzzy control method in 1996. The core idea of this method is to find the synthetical relationship among state variables by using modern control theory and form composed error and composed variety of error so as to construct a fuzzy controller.The parameter fuzzy control method is used to stabilize the above double inverted-pendulum physical system. The control parameters of this system are shown in Table5.Table 5
The control parameters of the parameter fuzzy controller.
Symbol
Value
k
1
8.3089
k
2
33.4311
k
3
−114.1740
k
4
0.0978
k
5
26.8817
k
6
−22.1896
a
0
0.3333
a
s
0.7518Real-time control is carried out in the above double inverted-pendulum physical system. When the initial state of the system isx
(
0
)
=
0 m, θ
1
(
0
)
=
0.05 rad, and θ
2
(
0
)
=
0.05 rad, the running time is 20 s. The comparison of the experimental results for two control methods is shown in Table 6.Table 6
The comparison of the control effects for two control methods.
Control method
Control effects
Evaluation indexes
Steady state deviation
Adjusting time (s)
The cart displacement (m)
The angle of the lower rod (rod)
The angle of the upper rod (rod)
Fuzzy control parameter method
0.125
0.01
0.004
2.37
Flexible logic control method
0.01
0.004
0.002
1.7From Table6, the adjusting time and the steady-state errors of the displacement and the angle are compared with parameter fuzzy control method about the stability control in double inverted-pendulum physical system.In summary, the flexible logic control model based on artificial neural networks has good control effects. And the controlling system has good stability and anti-interference ability.


## 5.1. Stability Control Experiment

Under the same initial state such asx
(
0
)
=
0 m, θ
1
(
0
)
=
0.05 rad and θ
2
(
0
)
=
0.05 rad, the stability control of the double inverted-pendulum physical system has been realized. The experimental results are shown in Figures 5, 6, 7, and 8 (running time: 20 s).The car displacement and speed curve when the system is steady.
(a)
The cart displacement curve
(b)
The cart speed curveThe lower rod angle curve and angular velocity curve when the system is steady.
(a)
The lower rod angle curve
(b)
The lower rod angular velocity curveThe upper rod angle curve and angular velocity curve when the system is steady.
(a)
The upper rod angle curve
(b)
The upper rod angular velocity curveFigure 8
The controlled quantity curve when the system is steady.From the experimental results, we can find when the system is in a stable state, the deviation of the cart displacement could be about 0.01 m with a good stability control effect.


## 5.2. Anti-Interference Control Experiment

After knocking the upper pendulum twice when the system reaches a steady state, it shows strong anti-interference ability. The experimental results are shown in Figures9, 10, 11, and 12.The car displacement and speed curve when the system is interfered.
(a)
The cart displacement curve
(b)
The cart speed curveThe lower rod angle curve and angular velocity curve when the system is interfered.
(a)
The lower rod angle curve
(b)
The lower rod angular velocity curveThe upper rod angle curve and angular velocity curve when the system is interfered.
(a)
The upper rod angle curve
(b)
The upper rod angular velocity curveFigure 12
The controlled quantity curve when the system is interfered.From the experimental results, we can find the system becomes stable again after 1.9 s.


## 5.3. Free Movement Control Experiment of Cart Pendulum System

When the system reaches a steady state, the target positionx
d is changed on-line, such as from 0 m to −0.2 m. So the free movement control of cart pendulum on a level track could be realized. The experimental results are shown in Figures 13, 14, 15, and 16.The car displacement and speed curve when the cart is moving in a level track.
(a)
The cart displacement curve
(b)
The cart speed curveThe lower rod angle curve and angular velocity curve when the cart is moving in a level track.
(a)
The lower rod angle curve
(b)
The lower rod angular velocity curveThe upper rod angle curve and angular velocity curve when the cart is moving in a level track.
(a)
The upper rod angle curve
(b)
The upper rod angular velocity curveFigure 16
The controlled quantity curve when the cart is moving in a level track.From the experimental results, we can see the cart moves to the new target location after 2.7 s and the system remains stable.


## 5.4. Self-Adaptive Experiment of Control Model

Some physical parameters have been changed in another double inverted-pendulum physical system. The physical parameters of the system are shown in Table4. The control model also can realize the stability control of the new double inverted-pendulum after changing some parameters. The experimental results are shown in Figures 17, 18, 19, and 20. The experimental results show that the control model based on artificial neural network has excellent self-adaptability and portability.Table 4
The physical parameters of another double inverted-pendulum.
Symbol
Value
m
0
0.595 kg
m
1
0.158 kg
m
2
0.142 kg
l
1
0.1237 m
l
2
0.2265 m
L
1
0.151 m
f
0
0.1 N·s/m
J
1
0.0028394 kg·m2
J
2
0.0024329 kg·m2The car displacement and speed curve when the system is steady.
(a)
The cart displacement curve
(b)
The cart speed curveThe lower rod angle curve and angular velocity curve when the system is steady.
(a)
The lower rod angle curve
(b)
The lower rod angular velocity curveThe upper rod angle curve and angular velocity curve when the system is steady.
(a)
The upper rod angle curve
(b)
The upper rod angular velocity curveFigure 20
The controlled quantity curve when the system is steady.


## 5.5. Comparison of the Experimental Results

Cheng et al. [20] put forward a parameter fuzzy control method in 1996. The core idea of this method is to find the synthetical relationship among state variables by using modern control theory and form composed error and composed variety of error so as to construct a fuzzy controller.The parameter fuzzy control method is used to stabilize the above double inverted-pendulum physical system. The control parameters of this system are shown in Table5.Table 5
The control parameters of the parameter fuzzy controller.
Symbol
Value
k
1
8.3089
k
2
33.4311
k
3
−114.1740
k
4
0.0978
k
5
26.8817
k
6
−22.1896
a
0
0.3333
a
s
0.7518Real-time control is carried out in the above double inverted-pendulum physical system. When the initial state of the system isx
(
0
)
=
0 m, θ
1
(
0
)
=
0.05 rad, and θ
2
(
0
)
=
0.05 rad, the running time is 20 s. The comparison of the experimental results for two control methods is shown in Table 6.Table 6
The comparison of the control effects for two control methods.
Control method
Control effects
Evaluation indexes
Steady state deviation
Adjusting time (s)
The cart displacement (m)
The angle of the lower rod (rod)
The angle of the upper rod (rod)
Fuzzy control parameter method
0.125
0.01
0.004
2.37
Flexible logic control method
0.01
0.004
0.002
1.7From Table6, the adjusting time and the steady-state errors of the displacement and the angle are compared with parameter fuzzy control method about the stability control in double inverted-pendulum physical system.In summary, the flexible logic control model based on artificial neural networks has good control effects. And the controlling system has good stability and anti-interference ability.


## 6. Conclusion

Based on the basic physical meaning of errorE and error variety E
C, this paper analyzes the logical relationship between them and uses universal combinational operation model to describe it. And a flexible logic control method is put forward to realize effective control on multivariable nonlinear system.In order to implement the fusion control ofUniversal Logic and artificial neural networks, this paper puts forward a new neuron model of zero-level universal combinational operation. Based on this neuron model, artificial neural network is designed for flexible logic control model. Meanwhile, the stability control of the double inverted-pendulum physical system, anti-interference control, and free movement of the cart pendulum have been realized in a level track by using the proposed control model.



---

*Source: 913549-2014-01-20.xml*