Approximation Properties and Continuous Differentiability of a Class of SISO Fuzzy Systems Barnabas Bede1,2
Matthew P. Peterson
Department of Mathematics DigiPen Institute of Technology 9931 Willows Rd. NE Redmond, WA, 98052, USA 2 Γbuda Univeristy BΓ©csi ΓΊt 96/B, H-134 Budapest, Hungary Email:1
[email protected] 2
[email protected] DigiPen Institute of Technology 9931 Willows Rd. NE Redmond, WA, 98052, USA Email:
[email protected] 1
AbstractβThe goal of the present paper is to investigate approximation and smoothness properties of Larsen type single input single output (SISO) fuzzy systems, that is, fuzzy logic systems using the maximum as aggregation for the individual rule outputs, product (Goguen) t-norm as the conjunctive operator and center of gravity defuzzification. We prove that the function providing the output of the above considered Larsen type fuzzy system is capable of approximating any continuous function. Also, it is continuously differentiable under very relaxed conditions as e.g. continuous differentiability of the antecendents except at their core, and continuous differentiability of the consequences of fuzzy rules except at their core and the enpoints of their support. We show practical examples regarding the approximation and smoothness of the Larsen type operators, showing also by an example that the conditions on the antecedent part cannot be in general weakened further without loosing continuous differentiability. The present paper provides a theoretical background for claims in literature stating that the output of a fuzzy control system of this type is smooth.
are smooth, so will be the output of the TS fuzzy system. The case of a fuzzy system of Larsen type is not immediate and it is interesting, because of the usage of the maximum operator which is known to destroy differentiability.
I. I NTRODUCTION
The fuzzy set π΄π can be seen as a model for the expression π₯ is about π₯π , π = 1, ..., π. The consequence part π΅π will be considered integrable and it is constructed such that π΅π (π¦π ) = 1 and it has the support
The study of approximation capability of fuzzy systems started from B. Koskoβs paper [7]. Most of the papers available in literature use the Takagi-Sugeno approach [13], [11], [2] or they are using sum as the aggregation method for the fuzzy rules [9], [14], [6]. In this situation smoothness properties have been investigated in e.g., [3]. Also, minβ aggregation was considered [4], [5], however, smoothness properties were not investigated for this case. Our goal in the present paper is to investigate the approximation problem of an unknown bounded function π : [π, π] β [πΌ, π½] by a SISO Fuzzy System of Larsen type, i.e., having maximum as aggregation operator of the individual rule outputs, product playing the role of conjunction, and the defuzzification being the most widespread method, i.e., the center of gravity. This problem of approximation capability was also considered in detail in the recent monograph [1]. In the present paper we also consider the problem of continuous differentiability of Larsen type SISO fuzzy systems. In many works it is mentioned that Mamdani type fuzzy controllers, [8] have a smooth output. Also, smoothness for Takagi Sugeno systems is relatively easily addressed [3]. If the antecedents
II. A PPROXIMATION PROPERTIES OF L ARSEN FUZZY SYSTEMS
Let us construct now the framework of our problem. Let π = π₯0 < π₯1 < ... < π₯π < π₯π+1 = π be a partition of the input domain. Let π (π₯π ) = π¦π , π = 0, ..., π + 1. We consider fuzzy systems with the support of each fuzzy set in the antecedent part, satisfying π΄π (π₯π ) = 1 and let us suppose that there exist π > 0, π β β, π < π, such that the π-level sets satisfy (π΄π )π β [π₯πβπ , π₯π+π ], π = 1, ..., π.
(π΅π )π β [min{π¦πβπ , ..., π¦π+π }, max{π¦πβπ , ..., π¦π+π }] for π = 1, ..., π. The fuzzy set π΅π can be seen as a model for the expression π¦ is about π¦π , π = 1, ..., π. The fuzzy rule base If π₯ is π΄π then π¦ is π΅π , π = 1, ..., π, helps us construct the Larsen SISO fuzzy system Β΄ π½ βπ [ π=1 π΄π (π₯) β
π΅π (π¦)] β
π¦ β
ππ¦ . πΉ (π, π₯) = πΌΒ΄ π½ β π [ π=1 π΄π (π₯) β
π΅π (π¦)] β
ππ¦ πΌ The following theorem is a generalization of Theorem 7.16 that was proposed recently in [1]. In the above cited Theorem the conditions imposed on the antecedents and consequences are stronger, namely, they restrict the supports of the antecedents and consequences to be small.
We will provide error estimates in terms of the modulus of continuity defined as
where π§πβπ = min{π¦πβπ , ..., π¦π+π },
π(π, πΏ) = sup{β£π (π₯) β π (π¦)β£ : β£π₯ β π¦β£ β€ πΏ}. Theorem 1: Any continuous function π : [π, π] β [πΌ, π½] can be approximated by the Larsen fuzzy system Β΄ π½ βπ [ π=1 π΄π (π₯) β
π΅π (π¦)] β
π¦ β
ππ¦ πΉ (π, π₯) = πΌΒ΄ π½ β π [ π=1 π΄π (π₯) β
π΅π (π¦)] β
ππ¦ πΌ
π§π+π = max{π¦πβπ , ..., π¦π+π }. Taking now into account that (π΅π )π β [min{π¦πβπ , ..., π¦π+π }, max{π¦πβπ , ..., π¦π+π }] we obtain
with any membership functions for the antecedents and consequences π΄π , π΅π , π = 1, ..., π such that there exist π > 0, π β β, π < π, such that (i) π΄π continuous, π΄π (π₯π ) = 1
β€
+
β€
π¦π = π (π₯π ), π = 1, ..., π. Moreover the following error estimate holds true β₯πΉ (π, π₯) β π (π₯)β₯ β€ 2ππ (π, πΏ) + π2 (π½ β πΌ)2 π
πΌ
.
π=1
β£πΉ (π, π₯) β π (π₯)β£ Β΄ π½ β [ π (π΄ (π₯) β
π΅ (π¦))] β
(π¦ β π (π₯))ππ¦ π πΌ π=1 π = Β΄ π½ βπ (π΄π (π₯) β
π΅π (π¦)) β
ππ¦ π=1
Β΄ π½ βπ β£[ π=1 (π΄π (π₯) β
π΅π (π¦))] β
(π¦ β π (π₯))β£ ππ¦ β€ πΌ . Β΄ π½ βπ π=1 (π΄π (π₯) β
π΅π (π¦)) β
ππ¦ πΌ By Lemma 7.10 in [1] we have β£πΉ (π, π₯) β π (π₯)β£ Β΄ π½ βπ [ (π΄π (π₯) β
π΅π (π¦))] β
β£π¦ β π (π₯)β£ ππ¦ β€ πΌ π=1 . Β΄ π½ βπ (π΄ (π₯) β
π΅ (π¦)) β
ππ¦ π π π=1 πΌ Also, any fixed value of π₯ β [π, π] belongs to the π-level set of at most 2π antecedents. Suppose that π₯ β (π΄πβπ )π βͺ ... βͺ (π΄π+π )π . Then we have Β΄ π§π+π βπ
β£πΉ (π, π₯) β π (π₯)β£
π=1 (π΄π (π₯) β
π΅π (π¦)) β
β£π¦ β π (π₯)β£ ππ¦ Β΄ π½ βπ π=1 (π΄π (π₯) β
π΅π (π¦)) β
ππ¦ πΌ Β΄ βπ π=1 (π΄π (π₯) β
π΅π (π¦)) β
β£π¦ β π (π₯)β£ ππ¦ [πΌ,π§πβπ ]βͺ[π§π+π ,πΌ] + , Β΄ π½ βπ π=1 (π΄π (π₯) β
π΅π (π¦)) β
ππ¦ πΌ
β€
π§πβπ
π2 (π½ β πΌ)2
π=1 (π΄π (π₯)
β
π΅π (π¦)) β
ππ¦
Finally we obtain
with
)β1
Proof: We have
πΌ
β
π΅π (π¦)) β
ππ¦
β£πΉ (π, π₯) β π (π₯)β£ β€ 2ππ(π, πΏ) + π2 (π½ β πΌ)2 π
π=1,...,π
π΄π (π₯) β
π΅π (π¦)] β
ππ¦
π=1 (π΄π (π₯)
β£π¦ β π (π₯)β£ ππ¦
π=1 (π΄π (π₯) β
π΅π (π¦)) β
π(π, 2ππΏ)ππ¦ Β΄ π½ βπ π=1 (π΄π (π₯) β
π΅π (π¦)) β
ππ¦ πΌ
πΌ
πΌ
πΏ = max {π₯π β π₯πβ1 }
[
[πΌ,π§πβπ ]βͺ[π§π+π ,π½]
Β΄ π½ βπ
+ Β΄ π½ βπ
with
π=
Β΄
πΌ
(π΅π )π β [min{π¦πβπ , ..., π¦π+π }, max{π¦πβπ , ..., π¦π+π }],
π π½ β
π2
Β΄ π½ βπ
(ii) π΅π integrable, π΅π (π¦π ) = 1,
(Λ
π=1 (π΄π (π₯) β
π΅π (π¦)) β
π(π, π§π+π β π§πβπ )ππ¦ Β΄ π½ βπ π=1 (π΄π (π₯) β
π΅π (π¦)) β
ππ¦ πΌ
π§πβπ
(π΄π )π β [π₯πβπ , π₯π+π ], π = 1, ..., π;
and
β£πΉ (π, π₯) β π (π₯)β£
Β΄ π§π+π βπ
(Λ π=
π π½ β
πΌ π=1
)β1 (π΄π (π₯) β
π΅π (π¦)) β
ππ¦
.
Corollary 2: A Larsen SISO fuzzy system πΉ (π, π₯) given as above, can approximate any continuous function π (π₯) with arbitrary accuracy. Proof: If πΏ β 0 and π β 0 then it is easy to see that πΉ (π, π₯) β π (π₯), and the result is immediate. Let us remark that results in [1] do not cover the important case of Gaussian antecedents and consequences. This case can be covered now easily. Corollary 3: A Larsen SISO fuzzy system πΉ (π, π₯) with Gaussian antecedents and consequences can approximate any continuous function π (π₯) with arbitrary accuracy. Proof: The proof is immediate from the previous Theorem, since for Gaussian membership functions we can find a set of parameter values π, π, π such that we have β£πΉ (π, π₯) β π (π₯)β£ arbitrarily small. III. C ONTINUOUS DIFFERENTIABILITY OF L ARSEN FUZZY SYSTEMS
The fact that a fuzzy system provides a smooth output is very intuitive and this was mentioned in several works as e.g. in [10]. In fact this is widely known as being one of the main advantages of fuzzy controllers of Mamdani types or Larsen types. Yet a proof of this property can not be found in the literature as to the best of the authors knowledge. In the present paper we address this problem and we prove that the Larsen type SISO fuzzy system described above, under very
relaxed conditions is continuously differentiable. The proof is elementary but non-trivial. Theorem 4: Let π : [π, π] β β be a monotone function and let π¦π = π (π₯π ), π = 1, ..., π.. We consider the Larsen type SISO fuzzy system Β΄ π½ βπ [ π=1 π΄π (π₯) β
π΅π (π¦)] β
π¦ β
ππ¦ πΉ (π, π₯) = πΌΒ΄ π½ β π [ π=1 π΄π (π₯) β
π΅π (π¦)] β
ππ¦ πΌ with any membership functions for the antecedents and consequences π΄π , π΅π , π = 1, ..., π satisfying (i) π΄π monotone increasing and differentiable on (ββ, π₯π ) and monotone decreasing and differentiable on (π₯π , β), with the closure of its support being (π΄π )0 = [π₯πβ1 , π₯π+1 ], π = 1, ..., π; (ii) π΅π strictly increasing and differentiable on [min{π¦πβ1 , π¦π , π¦π+1 }, π¦π ) and strictly decreasing and differentiable on (π¦π , max{π¦πβ1 , π¦π π¦π+1 }],
Λ +
π½
π(π₯)
Λ =
π(π₯)
π΄β²1 (π₯) β
π΅1 (π¦) β
π¦ππ¦ +
πΌ
and π β² (π₯) = Λ +
π΄β²2 (π₯) β
π΅2 (π¦) β
π¦ππ¦ β πβ² (π₯) β
π΄2 (π₯) β
π΅2 (π(π₯)) β
π(π₯)
Λ
πΌ
π½
π(π₯)
Λ
=
π(π₯)
Λ
π½
π(π₯)
π΄β²2 (π₯) β
π΅2 (π¦) β
π¦ππ¦
π΄β²1 (π₯)β
π΅1 (π¦)ππ¦ +πβ² (π₯)β
π΄1 (π₯)β
π΅1 (π(π₯))β
π(π₯)
π΄β²2 (π₯) β
π΅2 (π¦)ππ¦ β πβ² (π₯) β
π΄2 (π₯) β
π΅2 (π(π₯)) β
π(π₯) π(π₯)
πΌ
π΄β²1 (π₯) β
π΅1 (π¦)ππ¦ +
Λ
π½
π(π₯)
π΄β²2 (π₯) β
π΅2 (π¦)ππ¦.
As a conclusion the function πΉ (π, π₯) s differentiable and its derivative is πΉ β² (π, π₯) (π΅π )0 = [min{π¦πβ1 , π¦π , π¦π+1 }, max{π¦πβ1 , π¦π , π¦π+1 }], π = 1, ..., π; ) (Λ Λ π½ π(π₯) Then the Larsen type system given above is continuous and β² β² π΄1 (π₯) β
π΅1 (π¦) β
π¦ππ¦ + π΄2 (π₯) β
π΅2 (π¦) β
π¦ππ¦ = continuously differentiable function (class πΆ 1 ) on [π, π]. πΌ π(π₯) Proof: Let us remark that for the SISO fuzzy system ) (Β΄ Β΄π½ π(π₯) considered in this theorem, i.e., π΄ (π₯) β
π΅ (π¦)ππ¦ + π΄ (π₯) β
π΅ (π¦)ππ¦ 1 1 2 2 πΌ π(π₯) β
(Β΄ Β΄ π½ βπ )2 Β΄ π(π₯) π½ [ π=1 π΄π (π₯) β
π΅π (π¦)] β
π¦ β
ππ¦ π΄1 (π₯) β
π΅1 (π¦)ππ¦ + π(π₯) π΄2 (π₯) β
π΅2 (π¦)ππ¦ πΉ (π, π₯) = πΌΒ΄ π½ β πΌ π [ π=1 π΄π (π₯) β
π΅π (π¦)] β
ππ¦ (Λ ) πΌ Λ π½ π(π₯) β² β² only two rules are fired at a time, i.e., for any input value β π΄1 (π₯) β
π΅1 (π¦)ππ¦ + π΄2 (π₯) β
π΅2 (π¦)ππ¦ πΌ π(π₯) π₯ β [π, π] there exists π β {1, 2, ..., π β 1} such that π₯ β (Β΄ ) Β΄π½ (π΄π )0 β©(π΄π+1 )0 . Without restricting the generality we suppose π(π₯) π΄ (π₯) β
π΅ (π¦) β
π¦ππ¦ + π΄ (π₯) β
π΅ (π¦) β
π¦ππ¦ 1 1 2 2 πΌ π(π₯) that π = 1. β
(Β΄ )2 . Β΄π½ First let us assume that π₯ β (π₯1 , π₯2 ). Then we have π(π₯) π΄1 (π₯) β
π΅1 (π¦)ππ¦ + π(π₯) π΄2 (π₯) β
π΅2 (π¦)ππ¦ πΌ Β΄π½ (π΄ (π₯) β
π΅ (π¦) β¨ π΄ (π₯) β
π΅ (π¦)) β
π¦ππ¦ 1 1 2 2 πΉ (π, π₯) = πΌΒ΄ π½ For the case π₯ β (π₯π , π₯π+1 ) the proof is exactly the same (π΄1 (π₯) β
π΅1 (π¦) β¨ π΄2 (π₯) β
π΅2 (π¦))ππ¦ as above. πΌ Now let us consider the extreme case π₯ = π₯1 . Then since We have to prove that this function is differentiable with its derivative being continuous. Let us remark first that the π΄2 (π₯) is differentiable at π₯1 and since (π΄2 )0 = [π₯1 , π₯3 ] we have π΄β²2 (π₯1 ) = 0. Also we observe that π(π₯) is right consequences π΅1 and π΅2 overlap on (π₯1 , π₯2 ). We define π(π₯) to be the unique solution in π¦ of the equation continuous at π₯1 . We obtain the derivative to the right of π₯1 as (Β΄ ) (Β΄ ) π΄1 (π₯) β
π΅1 (π¦) = π΄2 (π₯) β
π΅2 (π¦). π½ β² π½ π΄ (π₯) β
π΅ (π¦) β
π¦ππ¦ π΄ (π₯) β
π΅ (π¦)ππ¦ 1 1 1 1 πΌ πΌ We can easily see that πΉπβ² (π, π₯) = (Β΄ )2 π½ ( )β1 ( ) π΄1 (π₯) β
π΅1 (π¦)ππ¦ πΌ π΅2 π΄1 (π₯) π(π₯) = π΅1 π΄2 (π₯) (Β΄ ) (Β΄ ) π½ π½ β² π΄ (π₯) β
π΅ (π¦) β
π¦ππ¦ π΄ (π₯) β
π΅ (π¦)ππ¦ 1 1 1 1 πΌ πΌ is differentiable and we can rewrite β = 0. (Β΄ )2 Β΄π½ Β΄ π(π₯) π½ π΄ (π₯) β
π΅ (π¦)π¦ππ¦ + π΄ (π₯) β
π΅ (π¦)π¦ππ¦) π΄ (π₯) β
π΅ (π¦)ππ¦ 1 1 2 1 1 πΌ π(π₯) 2 πΌ . πΉ (π, π₯) = Β΄ π(π₯) Β΄π½ π΄1 (π₯) β
π΅1 (π¦)ππ¦ + π(π₯) π΄2 (π₯) β
π΅2 (π¦)ππ¦) Similar reasoning can be performed at the left of π₯1 and πΌ we get πΉπβ² (π, π₯) = 0. The choice of π₯1 did not restrict the We denote that numerator as π (π₯) and the denominator as generality of our proof, and we obtain similarly continuous π(π₯). Then we observe that derivative at each π₯π , π = 1, ..., π and one-sided differentiabilΛ π(π₯) ity can be obtained at the endpoints π₯0 and π₯π+1 . Finally we π΄β²1 (π₯)β
π΅1 (π¦)β
π¦ππ¦+πβ² (π₯)β
π΄1 (π₯)β
π΅1 (π(π₯))β
π(π₯) π β² (π₯) = obtain that πΉ (π, π₯) is continuously differentiable. πΌ
Figure 1.
Figure 2.
Gaussian antecedents of a SISO fuzzy system
Gaussian consequences of a SISO fuzzy system
Figure 3. Larsen SISO fuzzy system associated to π (π₯) = π₯2 , Gaussian membership functions
Figure 4.
Smooth splines as antecedents of a SISO fuzzy system
IV. E XAMPLES Remark 5: It is an interesting remark that the fuzzy system constructed above does not have have the consequences π΅π differentiable, nor the target function π does not need to be differentiable and moreover it does not even need it to be continuous, yet the output will be differentiable. This indeed proves the hypothesys that we have provided at the begining of the paper that fuzzy systems and as a consequence fuzzy controllers have a smoothing effect in practical problems. Further let us remark that the requirement that π is monotone is restrictive in the proof of the previous theorem mention and it is necessary to have only two overlapping antecedents each time.
We exemplify our theoretical findings through a few examples. Surely a visual confirmation of smoothness properties is not necessarry, however it is useful to understand properties and limitations of Larsen type fuzzy systems. First we consider π : [0, 1] β [0, 1], π (π₯) = π₯2 and a Larsen fuzzy system with Gaussian Antecedents Fig. 1, Gaussian Consequences Fig. 2. This exemplifies Theorem 1 (see Fig. 3) To exemplify smooth monotonic spline antecedents and consequences (we use here the monotonic splines described in [12]). The antecedents used for our construction are shown
Figure 5.
Smooth splines as consequences of a SISO fuzzy system
Figure 6. Larsen SISO fuzzy system associated to π (π₯) = π₯2 smooth spline membership
in Fig. 4, the consequences are shown in Fig. 5 and the results are to be found in Fig. 7. Let us observe that differentiability of the antecedents π΄π , π = 1, ..., π is required at the endpoints of the support π₯πβ1 , π₯π+1 but differentiability at the core π₯π is not required. To examplify that we consider monotone splines that are not differentiable at the core, both as antecedents and consequences for a Larsen fuzzy system. We can also see in Fig. 10 a counterexample. If we consider triangular Antecedents and Consequences we do not obtain continuously differentiable output. Let us observe that in Theorem 4 besides monotony and boundedness, nothing is required on the target function π (π₯).
Figure 7.
Figure 8.
Antecedents that are not differentiable at the core
Consequences that are not differentiable at the core
As a conclusion we obtain continuously differentiable output even if the target function is discontinuous. We consider { 2 π₯ if π₯ < 0.5 π (π₯) = π₯ otherwise together with the spline antecedents in Fig. 7. The Larsen SISO fuzzy system in this case is illustrated in Fig. 11 V. C ONCLUSION AND FURTHER RESEARCH Approximation and Smoothness of Larsen fuzzy systems is investigated and our findings are that the output of a Larsen fuzzy system is able to approximate any continuous function with arbitrary accuracy and it is continuously differentiable, under very relaxed conditions. For future research we plan to
Figure 9. Larsen SISO fuzzy system associated to π (π₯) = π₯2 sharp splines
Figure 10. Larsen SISO fuzzy system associated to π (π₯) = π₯2 with triangular Antecedents and Consequences.
release the requirements in the Theorem such as monotonicity of th etarget function, etc. R EFERENCES [1] B. Bede, Mathematics of Fuzzy Sets and Fuzzy Logic, Springer, 2013. [2] B.S. Chen, H.C.H. Lee ; Y.C. Chang, π» β tracking design of uncertain nonlinear SISO systems: adaptive fuzzy approach, IEEE Transactions on Fuzzy Systems, 4(1996), 32 - 43. [3] M. Landajo, M. J. Rio, R. Perez, A note on smooth approximation capabilities of fuzzy systems, IEEE Transactions on Fuzzy Systems, 9(2001) 229-237. [4] Yong-Ming Li, Zhong-Ke Shi, Zhi-Hui Li, Approximation theory of fuzzy systems based upon genuine many-valued implications-SISO cases, Fuzzy Sets and Systems 130 (2002) 147β157.
Figure 11. Larsen SISO fuzzy system associated to a discontinuous function.
[5] Yong-Ming Li, Zhong-Ke Shi, Zhi-Hui Li, Approximation theory of fuzzy systems based upon genuine many-valued implications-MIMO cases, Fuzzy Sets and Systems 130 (2002) 159-174. [6] L. T. KΓ³czy and A. Zorat. Fuzzy systems and approximation. Fuzzy Sets and Systems, 85(1997) 203β222. [7] B. Kosko, Fuzzy Systems as Universal Approximators, IEEE Transactions on Computers, 43(1994), 1329-1333. [8] E.H. Mamdani, S. Assilian, An experiment in linguistic synthesis with a fuzzy logic controller, J. Man Machine Stud., 7(1975) 1-13. [9] S. Mitaim, B. Kosko, The shape of fuzzy sets in adaptive function approximation, IEEE Transactions on Fuzzy Systems, 9(2001) 637-656. [10] C Moewes, R Kruse, Fuzzy Control for Knowledge-Based Interpolation, Combining Experimentation and Theory, Springer, 2012 [11] A.H. Sonbol, M.S. Fadali, TSK Fuzzy Function Approximators: Design and Accuracy Analysis, IEEE Transactions on Systems Man and Cybernetics, 42(2012) 702-712. [12] L. Stefanini, L. Sorini, M.L. Guerra, Parametric representation of fuzzy numbers and application to fuzzy calculus, Fuzzy Sets and Systems 157 (2006) 2423 β 2455. [13] M. Sugeno, An introductory survey of fuzzy control Information Sciences, 36 (1985), 59β83. [14] D. Tikk, L.T. KΓ³czy, T.D. Gedeon, A survey on universal approximation and its limits in soft computing techniques, International Journal of Approximate Reasoning, 33(2003) 185202.