1 Introduction

A family of fuzzy numbers constitutes an important subclass of fuzzy sets having countless applications in all cases where imprecise real values are modeled by their fuzzy counterparts. To avoid problems in processing and calculations on fuzzy numbers described by complicated membership functions the suitable approximations are commonly applied. In particular, the interval (see Chanas 2001; Grzegorzewski 2002, 2012), triangular (see Abbasbandy et al. 2010; Ban 2011; Ban and Coroianu 2015a, 2016; Yeh 2017) or trapezoidal approximation (see Abbasbandy and Amirfakhrian 2006; Abbasbandy and Asady 2004; Abbasbandy and Hajjari 2009b; Ban 2008, 2009a, b; Ban et al. 2011a, b; Ban and Coroianu 2011, 2012, 2014; Coroianu 2011, 2012; Grzegorzewski 2008a, b, 2010; Grzegorzewski and Mrówka 2005, 2007; Grzegorzewski and Pasternak-Winiarska 2009, 2014; Yeh 2007, 2008a, b) is very popular, mainly because of simplicity of the output representation (by no more than four points). More recently (see Ban et al. 2016; Yeh and Chu 2011), the studies were extended by employing approximations by L-R fuzzy numbers. Fuzzy number approximation via shadowed sets was discussed in Grzegorzewski (2013). For the review of the fuzzy numbers approximation methods and their applications, we refer the reader to the monograph Ban et al. (2015).

One of the reasons that the trapezoidal fuzzy numbers are so popular in applications is that each such fuzzy number is represented completely by four real numbers only. Of course, the representation simplicity is a strong advantage but one may consider whether the shape reduction going so far is not too impoverish and restrictive in some situations. Therefore, the problem of the piecewise linear approximation of fuzzy numbers by the so-called 1-knot fuzzy numbers was considered in Coroianu et al. (2013). Each such 1-knot fuzzy number is completely characterized by six points on the real line. This way we get fuzzy numbers which are still simple enough, but simultaneously, having more “degrees of freedom”, we may obtain approximations that are more flexible to preserve some important properties of the input fuzzy numbers.

In this paper, a generalization of the results presented in Coroianu et al. (2013) is given. Now, instead of six points on the real line and piecewise linear sides each consisting of two segments that characterize 1-knot fuzzy numbers, we consider n-knot fuzzy numbers (where \(n\ge 2\)) which enables to quantify the uncertainty at n intermediate levels between 0 and 1. Such fuzzy numbers were already introduced in paper Báez-Sánchez et al. (2012) and were called polygonal fuzzy numbers. In this way, we obtain a subfamily of fuzzy numbers with piecewise linear membership functions, where each side consists of \(n+1\) segments. Hence, the output of the approximation is still simple but more flexible for reconstructing a fuzzy input than obtained using other methods discussed above. Moreover, it turns out that basic characteristics of the approximations converge to corresponding characteristics of the input fuzzy numbers. When comparing n-knot approximation with the 1-knot, discussed in Coroianu et al. (2013), its natural disadvantage is the longer processing time caused by more knots but its obvious advantage is better accuracy.

It should be stressed that the n-knot fuzzy number is the most natural and desired fuzzy structure for the computer representation which is by definition discrete, even if the original object it represents is a continuous one. Actually, a typical way a fuzzy object is stored in a database is a set of its \(\alpha \)-cuts. Here a natural question arises about the adequate choice of those \(\alpha \)-cuts, even for a fixed set of \( \alpha \)-levels. The method of the nearest piecewise linear approximation of fuzzy numbers proposed in this paper enables to avoid unjustified subjectivity and to choose the required \(\alpha \)-cuts in an objective way based on the nearness criteria crucial to any reasonable approximation.

The paper is organized as follows. In Sect. 2, we recall basic terminology connected with fuzzy numbers and define the \(\alpha \) -piecewise linear n-knot fuzzy numbers. In Sect. 3, we introduce some auxiliary results and present a convenient reparametrization of piecewise linear fuzzy numbers useful for solving the approximation problem. Then in Sect. 4, we discuss briefly the existence and properties of the nearest piecewise linear fuzzy number for the Euclidean distance and a fixed knot setting.

Next section, i.e., Sect. 5, contains a broad study on the convergence results concerning the approximation operator. In particular, we consider some special cases, like the so-called naïve approximator (i.e., an n-knot piecewise linear fuzzy number that interpolates the sides of a fuzzy number at the knots) or the approximator with equidistant knots. We prove there some theorems on the rate of convergence but we also examine the stability of basic characteristics like the expected interval, expected value, value of fuzzy number and its ambiguity.

Section 6 is both theoretical and strongly user-oriented as well. We give there not only practical approximation algorithms and illustrative examples but we also provide a simulation study on the approximation accuracy of the computer calculations on fuzzy numbers and stability of some fuzzy number characteristics (for the practical implementation of the presented algorithms we refer the reader to FuzzyNumbers package for R by Gagolewski 2015). Therefore, Sect. 6 appears in some sense as the core of the paper by linking theory with applications as the latter ones are so desired in real-world problems.

Finally, Sect. 7 concludes the paper. Some open problems and directions for the further research are also sketched there.

2 Piecewise linear fuzzy numbers

Fuzzy numbers are particular cases of convex fuzzy sets on the real line. The membership function of a fuzzy number A is given by

$$\begin{aligned} A(x)=\left\{ \begin{array}{lll} 0 &{}\quad \text {if} &{} x<a_{1}, \\ l_{A}(x) &{}\quad \text {if} &{} a_{1}\le x<a_{2}, \\ 1 &{}\quad \text {if} &{} a_{2}\le x\le a_{3}, \\ r_{A}(x) &{}\quad \text {if} &{} a_{3}<x\le a_{4}, \\ 0 &{}\quad \text {if} &{} x>a_{4}, \end{array} \right. \end{aligned}$$
(1)

where \(a_{1},a_{2},a_{3},a_{4}\in {\mathbb {R}}\), \(l_{A}:\mathbb {[} a_{1},a_{2}]\longrightarrow [0,1]\) is a nondecreasing upper semicontinuous function, \(l_{A}(a_{1})=0\), \(l_{A}(a_{2})=1\), called the left side of the fuzzy number, and \(r_{A}:[a_{3},a_{4}] \longrightarrow [0,1]\) is a nonincreasing upper semicontinuous function, \(r_{A}(a_{3})=1\), \(r_{A}(a_{4})=0\), called the right side of the fuzzy number A. The \(\alpha \)-cut of A, \(\alpha \in (0,1] \), is a crisp set defined as

$$\begin{aligned} A_{\alpha }=\{x\in {\mathbb {R}}:A(x)\ge \alpha \}. \end{aligned}$$

The case of the 0-cut is defined in slightly different way, and it actually coincides with the support of the fuzzy number. More exactly we have

$$\begin{aligned} \mathrm {supp}(A)=A_{0}=\overline{\{x\in {\mathbb {R}}:A(x)>0\}}. \end{aligned}$$

From the convexity property, it easily results that every \(\alpha \)-cut of a fuzzy number is a compact interval

$$\begin{aligned} A_{\alpha }=[A_{L}(\alpha ),A_{U}(\alpha )], \end{aligned}$$

where \(A_{L}(\alpha )=\inf \{x\in {\mathbb {R}}:A(x)\ge \alpha \}\) and \( A_{U}(\alpha )=\sup \{x\in {\mathbb {R}}:A(x)\ge \alpha \}\).

The 1-cut of A is also called the core of A and the following notation is commonly used

$$\begin{aligned} A_{1}=\mathrm {core}(A)=\{x\in {\mathbb {R}}:A(x)=1\}. \end{aligned}$$

When \(\mathrm {core}(A)\) is reduced to a single element, then A is called a unimodal fuzzy number [in some works like the book Hanss (2005), the authors reduce the term fuzzy numbers to unimodal fuzzy numbers only and call multimodal fuzzy numbers as fuzzy intervals]. Further on we denote by \({\mathbb {F}}(\mathbb {R)}\) the set of all fuzzy numbers.

However, fuzzy numbers with simple membership functions are often preferred in practice. For example, triangular or trapezoidal fuzzy numbers are most often used to rank fuzzy numbers (see, e.g., Abbasbandy and Hajjari 2009a; Ban and Coroianu 2015b; Facchinetti and Ricci 2004). Then, the same classes or even more generally the well-known class of L-R fuzzy numbers are used in fuzzy arithmetic (see, e.g., Carlsson and Fullér 2011; Hanss 2005; Hong and Hwang 1997; Kolesárová 1995). We can also mention here the Bodjanova (2005) fuzzy numbers used in statistical problems or in multicriteria decision making (see Ban and Ban 2012). Finally, we mention the recently introduced so-called parametric fuzzy numbers also known as semi-trapezoidal fuzzy numbers (see Nasibov and Peker 2008; Yeh 2009, 2011). Apart from the aforementioned applications, such fuzzy numbers are very suitable in approximation as we already discussed this in Introduction. Another subclass of \({\mathbb {F}}(\mathbb {R)}\), very useful and convenient especially in computer processing, may be defined by considering fuzzy numbers with piecewise linear sides. We consider a particular type of such fuzzy numbers introduced in Báez-Sánchez et al. (2012) and being referred there as polygonal fuzzy numbers. Consider the following definition.

Definition 1

(see Definition 5 and Example 6 in Báez-Sánchez et al. 2012) Fix \( n\in {\mathbb {N}}_{0}\), and let \({\mathfrak {A}}_{n}=\{(\alpha _{0},\alpha _{1},\dots ,\alpha _{n+1})\in {[0,1]}^{n+2}:0=\alpha _{0}<\alpha _{1}<\dots<\alpha _{n}<\alpha _{n+1}=1)\}\), and \({\mathfrak {S}}_{n}=\{(s_{1},\dots ,s_{2n+4})\in {\mathbb {R}}^{2n+4}:s_{1}\le \dots \le s_{2n+4}\}\). Given \({\varvec{\alpha }}\in {\mathfrak {A}}_{n}\) and \({\mathbf {s}} \in {\mathfrak {S}}_{n}\), an \({\varvec{\alpha }}\)-piecewise linear n-knot fuzzy number (also known as polygonal fuzzy number according to Báez-Sánchez et al. 2012) \(S({\varvec{\alpha }},{\mathbf {s}})\) is defined by:

$$\begin{aligned} S({\varvec{\alpha }},{\mathbf {s}})_{L}(\beta )= & {} s_{i+1}+(s_{i+2}-s_{i+1})\nonumber \\&\quad \frac{\beta -\alpha _{i}}{\alpha _{i+1}-\alpha _{i}}, \end{aligned}$$
(2)
$$\begin{aligned} S({\varvec{\alpha }},{\mathbf {s}})_{U}(\beta )= & {} s_{2n+4-i}+(s_{2n+3-i}-s_{2n+4-i})\nonumber \\&\quad \frac{\beta -\alpha _{i}}{\alpha _{i+1}-\alpha _{i}}, \end{aligned}$$
(3)

for some \(i\in \{0,\dots ,n\}\) such that \(\beta \in \left[ \alpha _{i},\alpha _{i+1}\right] \).

Please note that by the connection between the functions \(l_{S({\varvec{\alpha }}, {\mathbf {s}})}, r_{S({\varvec{\alpha }}, {\mathbf {s}})}\) and \(S( {\varvec{\alpha }}, {\mathbf {s}})_{L}, S({\varvec{\alpha }}, {\mathbf {s}})_{U}\) , the membership function of \(S({\varvec{\alpha }}, {\mathbf {s}})\) is also piecewise linear in the case when \({\mathbf {s}}\) is strictly monotone (as an example see Fig. 1).

Fig. 1
figure 1

Plot of an exemplary 3-knot piecewise linear fuzzy number \(S({\varvec{\alpha }}, {\mathbf {s}})\), with \({\varvec{\alpha }} =(0.3, 0.5, 0.7)\) and \({\mathbf {s}}=(1,1.5,2,2.4,2.5,4,5,5.5,6.5,7)\)

Let \({\mathbb {F}}^{\pi _{n}({\varvec{\alpha }})}(\mathbb {R)}\) denote the set of all \({\varvec{\alpha }}\)-piecewise linear n-knot fuzzy numbers (for fixed n and \({\varvec{\alpha }}\)). It is worth noting that the class introduced in Definition 1 generalizes well-known subfamilies of fuzzy numbers. Actually, for \(n=0\) and \(s_1=s_4\) we get “crisp” real numbers, for \(n=0\) and \(s_1=s_2, s_3=s_4\) we obtain “crisp” real intervals; if \(n=0\) and \(s_2=s_3\), we get triangular fuzzy numbers, assuming only \(n=0\) we obtain trapezoidal fuzzy numbers, while for \(n=1\) we receive 1-knot piecewise linear fuzzy numbers, discussed in Coroianu et al. (2013).

Further on we assume that two fuzzy numbers A and B are equal (and denote it as \(A=B\)) if \(A_{L}(\beta )=B_{L}(\beta )\) and \(A_{U}(\beta )=B_{U}( \beta )\) almost everywhere, \(\beta \in [0,1]\).

Let \(A,B\in {\mathbf {F}}\left( {\mathbb {R}}\right) \), and \(\lambda \in {\mathbb {R}}\). We define the sum \(A+B\) and the scalar multiplication \(\lambda \cdot A\) (see, e.g., Diamond and Kloeden 1994, p. 40) for every \(\beta \in [0,1]\) as:

$$\begin{aligned} \left( A+B\right) _{\beta }= & {} A_{\beta }+B_{\beta } \\= & {} \left[ A_{L}\left( \beta \right) +B_{L}\left( \beta \right) , A_{U}\left( \beta \right) +B_{U}\left( \beta \right) \right] \end{aligned}$$

and

$$\begin{aligned} \left( \lambda \cdot A\right) _{\beta }=\lambda A_{\beta }=\left\{ \begin{array}{ll} \left[ \lambda A_{L}\left( \beta \right) ,\lambda A_{U}\left( \beta \right) \right] &{}\quad \text {if }\lambda \ge 0, \\ \left[ \lambda A_{U}\left( \beta \right) ,\lambda A_{L}\left( \beta \right) \right] &{}\quad \text {if }\lambda <0. \end{array} \right. \end{aligned}$$

In the case of \({\varvec{\alpha }}\)-piecewise linear n-knot fuzzy numbers \(S^{\prime }=S({\varvec{\alpha }}, {\mathbf {s}}^{\prime })\) and \(S^{\prime \prime }=S({\varvec{\alpha }},{\mathbf {s}}^{\prime \prime })\), we obtain

$$\begin{aligned} S^{\prime }+S^{\prime \prime }=S({\varvec{\alpha }}, {\mathbf {s}}^{\prime }+ {\mathbf {s}}^{\prime \prime }). \end{aligned}$$

Moreover, we have

$$\begin{aligned} \lambda S=\left\{ \begin{array}{ll} S({\varvec{\alpha }}, \lambda {\mathbf {s}}) &{}\quad \text {if }\lambda \ge 0, \\ S({\varvec{\alpha }}, \lambda (s_{2n+4}, s_{2n+3},\dots ,s_{1})) &{}\quad \text {if } \lambda <0. \end{array} \right. \end{aligned}$$

Thus, \({\mathbb {F}}^{\pi _{n}({\varvec{\alpha }})}(\mathbb {R)}\) is closed under these two extension principle-based arithmetic operations.

3 Some auxiliary results

Many papers show that the most suitable and commonly used metric for fuzzy numbers approximation is an extension of the Euclidean distance d defined by (see Grzegorzewski 1998)

$$\begin{aligned} d^{2}(A,B)= & {} \int _{0}^{1}(A_{L}(\beta )-B_{L}(\beta ))^{2}\mathrm{{d}}\beta \nonumber \\&+\int _{0}^{1}(A_{U}(\beta )-B_{U}(\beta ))^{2}\mathrm{{d}}\beta . \end{aligned}$$
(4)

Let \(L^{2}[0,1]\) denote the space of square integrable functions on [0, 1]. It is well known (cf. Yeh 2011) that we can embed the space \( \left( {\mathbb {F}}(\mathbb {R)},d,+,\cdot \right) \) into the Hilbert space \( \left( L^{2}[0,1]\times L^{2}[0,1],{\widetilde{d}},\oplus ,\odot \right) \), such that for \({\mathbf {f}}=(f_{1},f_{2})\), \({\mathbf {g}}=(g_{1},g_{2})\in L^{2}[0,1]\times L^{2}[0,1]\) we have

$$\begin{aligned} {\widetilde{d}}^{2}(\mathbf {f,g)}= & {} \int _{0}^{1}\big (f_{1}(\alpha )-g_{1}(\alpha )\big )^{2}\,\mathrm{{d}}\alpha \nonumber \\&+\int _{0}^{1}\big (f_{2}(\alpha )-g_{2}(\alpha )\big )^{2}\,\mathrm{{d}}\alpha , \end{aligned}$$
(5)

and

$$\begin{aligned} A\oplus B= & {} A+B, \\ \lambda \odot A= & {} \lambda \cdot A \end{aligned}$$

for any \(A,B\in {\mathbb {F}}(\mathbb {R)}\) and \(\lambda \in [0,\infty )\), as \((A_{L},A_{U})\), \((B_{L},B_{U})\in L^{2}[0,1]\times L^{2}[0,1]\) (for the proof see, e.g., Yeh 2011). Because there is no risk of confusion, in all what follows, we use \(+\) and \(\cdot \) to define the addition and scalar multiplication in the space \(L^{2}[0,1]\times L^{2}[0,1]\), instead of \(\oplus \) and \(\odot \).

Note that the inner product which generates \({\widetilde{d}}\) is given by (see, e.g., Yeh 2008a)

$$\begin{aligned} \left\langle {\mathbf {f}},{\mathbf {g}}\right\rangle =\int _{0}^{1}f_{1}(\alpha )\,g_{1}(\alpha )\,\mathrm{{d}}\alpha +\int _{0}^{1}f_{2}(\alpha )\,g_{2}(\alpha )\,\mathrm{{d}}\alpha . \end{aligned}$$
(6)

Assuming a fixed n and a fixed family of \(\alpha \)-cuts \({\varvec{\alpha }} =(\alpha _0,\dots ,\alpha _{n+1})\in {\mathfrak {A}}_n\), let us define the following set of vectors \({\mathbf {e}}_1,\dots ,{\mathbf {e}}_{2n+4}\in L^{2}[0,1]\times L^{2}[0,1]\), \({\mathbf {e}}_i = (e_{i,1}, e_{i,2})\), \( i=1,\dots ,2n+4\), such that for \(\beta \in [0,1]\):

  • \(e_{1,1}(\beta ) = e_{1,2}(\beta )=1\),

  • for \(i=2,\ldots ,n+2\):

    $$\begin{aligned} e_{i,1}(\beta )=&\left\{ \begin{array}{lll} 0 &{}\quad \text {for} &{} \beta <\alpha _{i-2}, \\ \frac{\beta -\alpha _{i-2}}{\alpha _{i-1}-\alpha _{i-2}} &{}\quad \text {for} &{} \beta \in [\alpha _{i-2},\alpha _{i-1}], \\ 1 &{}\quad \text {for} &{} \beta >\alpha _{i-1}, \end{array} \right. \\ e_{i,2}(\beta )=&1, \end{aligned}$$
  • \(e_{n+3,1}(\beta ) = 0\), \(e_{n+3,2}(\beta )=1\),

  • for \(i=n+4,\ldots ,2n+4\):

    $$\begin{aligned} e_{i,1}(\beta )=&0, \\ e_{i,2}(\beta )=&\left\{ \begin{array}{lll} 1 &{}\quad \text {for} &{} \beta <\alpha _{2n-i+4}, \\ \frac{\alpha _{2n-i+5}-\beta }{\alpha _{2n-i+5}-\alpha _{2n-i+4}} &{}\quad \text {for} &{} \beta \in [\alpha _{2n-i+4},\\ &{} &{} \alpha _{2n-i+5}], \\ 0 &{}\quad \text {for} &{} \beta >\alpha _{2n-i+5}, \end{array} \right. \end{aligned}$$

The set of vectors defined above is extremely useful for a reparametrization of piecewise linear fuzzy numbers, which appears very convenient in further considerations. Namely, given any piecewise linear n-knot fuzzy number \(S= \mathrm {S}({\varvec{\alpha }},{\mathbf {s}})\), let \({\varvec{\delta }} =(\delta _1,\dots ,\delta _{2n+4})\) with \(\delta _1=s_1\) and \(\delta _i := s_{i}-s_{i-1} \) for \(i=2,\dots ,2n+4\). It can be seen that we have \( S=\sum _{i=1}^{2n+4} \delta _i\,{\mathbf {e}}_i\). From now on we use the notation \(\mathrm {S}_d({\varvec{\alpha }}, {\varvec{\delta }})\) to represent S in its new parametrization.

Lemma 1

The set \(\{{\mathbf {e}}_{1},{\mathbf {e}}_{2},\dots ,{\mathbf {e}}_{2n+4}\}\) is linearly independent in \(L^{2}[0,1]\times L^{2}[0,1]\).

Proof

First of all let us note that we have

$$\begin{aligned} {\mathbb {F}}^{\pi _{n}({\alpha })}(\mathbb {R)}=\left\{ \sum \limits _{i=1}^{2n+4}\delta _{i}{\mathbf {e}}_{i}:\delta _{1}\in {\mathbb {R}} ,\delta _{2},\ldots ,\delta _{2n+4}\in [0,\infty )\right\} .\nonumber \\ \end{aligned}$$
(7)

Therefore, to obtain the desired conclusion it suffices to find \(2n+4\) vectors in \({\mathbb {F}}^{\pi _{n}({\alpha })}(\mathbb {R)}\) which are linearly independent in \(L^{2}[0,1]\times L^{2}[0,1]\). Obviously such vectors exist as for example \(\mathrm {S}_{d}({\alpha },\varvec{ \xi }_{i})\), \(i=1,\ldots ,2n+4\), where \(\{\varvec{\xi }_{1},\ldots ,\varvec{\xi }_{2n+4}\}\) is the canonical base of \({\mathbb {R}}^{2n+4}\). \(\square \)

Now we are in position to present the main result of this section.

Lemma 2

For any \({\varvec{\alpha }}\in {\mathfrak {A}}_n\) the set \({\mathbb {F}}^{\pi _{n}({\varvec{\alpha }})}(\mathbb {R)}\) is a closed convex subset of the space \(L^{2}[0,1]\times L^{2}[0,1]\) endowed with the topology generated by metric \({\widetilde{d}}\).

Proof

If we denote \({\mathbb {V}}=\mathrm {span}\{{\mathbf {e}}_{1},{\mathbf {e}}_{2},\dots , {\mathbf {e}}_{2n+4}\}\) then let us consider the linear bijective transformation \(i:{\mathbb {R}}^{2n+4}\rightarrow {\mathbb {V}}\), \( i(x_{1},x_{2},\ldots ,x_{2n+4}\mathbf {)=}\sum \nolimits _{i=1}^{2n+4}x_{i}{\mathbf {e}} _{i}\). We consider on \({\mathbb {V}}\) the same metric d as on \( L^{2}[0,1]\times L^{2}[0,1]\) (hence by (6) \({\mathbb {V}}\) is an inner product space) and \({\mathbb {R}}^{2n+4}\) is endowed with the Euclidean topology. Since i and \(i^{-1}\) are linear transformations between finite-dimensional normed spaces, it is immediate that both are continuous and hence homeomorphisms. We observe that \({\mathbb {F}}^{\pi _{n}( {\varvec{\alpha }})}(\mathbb {R)=}i(\Omega )\) where \(\Omega =\{ (x_{1},\ldots ,x_{2n+4})\in {\mathbb {R}}^{2n+4}:x_{i}\ge 0,i=2,\ldots ,2n+4\}\). Obviously \(\Omega \) is closed in \({\mathbb {R}} ^{2n+4}\) (actually \(\Omega \) is a polyhedral subset of \({\mathbb {R}} ^{2n+4}\) and all such sets are closed in \({\mathbb {R}}^{2n+4}\)) and since \(\ i^{-1}\) is continuous it results that \({\mathbb {F}}^{\pi _{n}(\varvec{ \alpha })}(\mathbb {R)=}i(\Omega )\) is a closed subset of \({\mathbb {V}}\). On the other hand, \({\mathbb {V}}\) is a finite linear subspace of \(L^{2}[0,1]\times L^{2}[0,1]\) and hence \({\mathbb {V}}\) is closed in \(L^{2}[0,1]\times L^{2}[0,1]\) . By elementary topology, it easily results now that \({\mathbb {F}}^{\pi _{n}( {\varvec{\alpha }})}(\mathbb {R)}\) is closed in \(L^{2}[0,1]\times L^{2}[0,1]\).

As the convexity of \({\mathbb {F}}^{\pi _{n}({\varvec{\alpha }})}(\mathbb {R)}\) is obvious, the proof is complete. \(\square \)

Remark 1

From the above proof, it follows that a sequence of \({\varvec{\alpha }}\)-piecewise linear n-knot fuzzy numbers, \(\left( \mathrm {S}_{d}({\varvec{\alpha }},\varvec{\delta }_{m})\right) _{m=1,2,\ldots }\) with \(\varvec{\delta }_{m}=(\delta _{m,i})_{i=1,\dots ,2n+4}\), converges to \(\mathrm {S}_{d}({\varvec{\alpha }},\varvec{\delta })\), where \(\varvec{\delta }=(\delta _{i})_{i=1,\dots ,2n+4}\) if and only if for any i we have \(\delta _{m,i}\rightarrow \delta _{i}\), which is equivalent with \(\left\| \varvec{\delta }_{m}-\varvec{ \delta }\right\| \rightarrow 0\) (here \(\left\| \cdot \right\| \) denotes the Euclidean norm over the space \({\mathbb {R}}^{2n+4}\)).

Remark 2

Lemma 2 can be deduced from Corollary 12 in Báez-Sánchez et al. (2012). Although in Báez-Sánchez et al. (2012), a different metric is employed one can easily prove that this metric is equivalent with the Euclidean metric on the span of \({\mathbb {F}}^{\pi _{n}({\varvec{\alpha }})}(\mathbb {R) }\). However, we prefer the present proof because it is more suitable to our forthcoming results. Especially, we refer here to the set \(\{{\mathbf {e}}_{1}, {\mathbf {e}}_{2},\dots ,{\mathbf {e}}_{2n+4}\}\) of linearly independent vectors in \(L^{2}[0,1]\times L^{2}[0,1]\) which will be used to approach the main results of the paper.

4 The best approximation for fixed \(\varvec{\alpha }\)

In this section, we generalize the theoretical results from paper Coroianu et al. (2013) concerning the existence, uniqueness and characterization of the approximation as well as the properties of the derived approximation operator by extending them to the case of piecewise linear n-knot fuzzy numbers.

If membership functions of fuzzy numbers under study are too complicated, we usually approximate them by some simpler forms that are more useful for processing and easier to interpret. In this section, we discuss how to approximate an arbitrary fuzzy number by a piecewise linear n-knot fuzzy number described by fixed \(\alpha \)-cuts.

Hence, given a fuzzy number \(A\in {\mathbb {F}}({\mathbb {R}})\), fixed n, and a fixed family of \(\alpha \)-cuts \({\varvec{\alpha }}\in {\mathfrak {A}}_n\), we are interested in finding \(S_{{\varvec{\alpha }}}^{*}(A)\in {\mathbb {F}} ^{\pi _{n}({\varvec{\alpha }})}(\mathbb {R)}\) such that:

$$\begin{aligned} d(A,S_{{\varvec{\alpha }}}^{*}(A)) =\min \limits _{S\in {\mathbb {F}} ^{\pi _{n}({\varvec{\alpha }})}(\mathbb {R)}} d(A,S). \end{aligned}$$
(8)

Using similar reasoning as in paper Coroianu et al. (2013), one can prove that the above problem has a unique solution. By showing this fact, we are able to define the approximation operator \(\Pi _{{{\varvec{\alpha }}} }^{n}:{\mathbb {F}}({\mathbb {R}})\rightarrow {\mathbb {F}}^{\pi _{n}(\varvec{ \alpha })}({\mathbb {R}})\) which assigns to a fuzzy number A its unique best piecewise linear approximator (for given n and \({\varvec{\alpha }}\)).

Theorem 1

If A is an arbitrary fuzzy number, then there exists the unique fuzzy number \(\Pi _{{\varvec{\alpha }}}^{n}(A)\in {\mathbb {F}}^{\pi _{n}({\varvec{\alpha }})}(\mathbb {R)}\) satisfying (8).

Taking into account Lemma 2, we obtain the proof (we omit this proof since it is identical with the proof of Theorem 8 in Coroianu et al. 2013).

Now let us denote with \(\Phi =(\phi _{i,j})_{i,j=1,\dots ,2n+4}\), the Gram matrix associated with the set \(\{{\mathbf {e}}_{1},\dots ,{\mathbf {e}}_{2n+4}\}\), i.e., \(\phi _{i,j}=\left\langle {\mathbf {e}}_{i},{\mathbf {e}}_{j}\right\rangle \). Since these vectors are linearly independent, it follows that \(\Phi \) is invertible. Moreover, for fixed A let \({\mathbf {b}}=(b_1,\dots ,b_{2n+4})\) such that \(b_{i}=\left\langle A,{\mathbf {e}}_{i}\right\rangle \).

We have the following characterization of the best approximation (see, e.g., Yeh 2009, Fact 2.1). If \((X,\left\langle \cdot ,\cdot \right\rangle )\) is a Hilbert space, \(\Omega \) is a closed convex subset of X, and \(x\in X\), then \({x}^*\in \Omega \) is the unique best approximation of x relatively to the set \(\Omega \) if and only if \( \left\langle x-{x}^*,y-{x}^*\right\rangle \)\(\le 0\) for any \(y\in \Omega \). Note that the notation \({x}^*=P_{\Omega }(x)\) is often used to denote that \({x}^*\) is the projection of x onto \(\Omega \).

Theorem 2

Let \({\varvec{\alpha }}\in {\mathfrak {A}}_n\) be fixed and let A denote a fuzzy number. Then \(\Pi _{{\varvec{\alpha }}}^{n}(A)= \mathrm {S}_{d}({\varvec{\alpha }},{\varvec{\delta }}^*)\), \(\varvec{ \delta }^*\in {\mathbb {R}}\times {\mathbb {R}}_{+}^{2n+3}\), is the unique best approximation of A relatively to the set \({\mathbb {F}}^{\pi _{n}({\varvec{\alpha }})}(\mathbb {R)}\) with respect to the Euclidean metric d, if and only if there exists a vector \({\mathbf {z}}^*\in {\mathbb {R}}^{2n+4}\) such that all the following requirements hold

  1. (i)

    \(\Phi \,{{\varvec{\delta }}^*}^T - {{\mathbf {z}}^*}^T={\mathbf {b}}^T\),

  2. (ii)

    \(z^*_{1}=0\) and \(z^*_{i}\ge 0\), \(\forall i>1\),

  3. (iii)

    \(\delta _{i}^{*}=0\) or \(z_{i}^{*}=0\), \(\forall i>1\).

We omit the proof because it uses a similar reasoning as the proof of Theorem 9 in Coroianu et al. (2013). Basically the only difference is that now we have a Gram matrix of dimension \(2n+4\).

It turns out that our approximation operator has some very important properties from the well-known list presented in Grzegorzewski and Mrówka (2005).

Theorem 3

For any \({\varvec{\alpha }}\in {\mathfrak {A}}_n\), the \(\Pi _{{\varvec{\alpha }}}^{n}\) operator fulfills the following properties:

  1. (i)

    Identity, i.e., \(\Pi _{{\varvec{\alpha }}}^{n}(A)=A\), \(\forall A\in {\mathbb {F}}^{\pi _{n}({\varvec{\alpha }})}(\mathbb {R)}\);

  2. (ii)

    Invariance to translation, i.e.,

    \(\Pi _{{\varvec{\alpha }}}^{n}(A+z)=\Pi _{{\varvec{\alpha }}}^{n}(A)+z\), \(\forall A\in {\mathbb {F}}({\mathbb {R}})\), \(\forall z\in {\mathbb {R}}\);

  3. (iii)

    Scale invariance, i.e.,

    \(\Pi _{{\varvec{\alpha }}}^{n}(\lambda \,A) =\lambda \, \Pi _{{\varvec{\alpha }}}^{n}(A)\)f, \(\forall A\in {\mathbb {F}}(\mathbb {R)}\), \(\forall \lambda \in {\mathbb {R}}\);

  4. (iv)

    Lipschitz-continuity, i.e.,

    \(d(\Pi _{{\varvec{\alpha }}}^{n}(A),\Pi _{{\varvec{\alpha }}}^{n}(B))\le d(A,B)\), \(\forall A,B\in {\mathbb {F}}(\mathbb {R)}\).

The proofs are identical with those from the particular case of piecewise linear 1-knot approximation (see Theorem 10 in Coroianu et al. 2013).

5 Some remarks on convergence

For the piecewise linear approximation operator introduced in the previous section, we can prove useful approximation results. First of all we can find a rate of convergence for the sequence of piecewise linear approximations by letting \(n\rightarrow \infty \). Additionally, we can estimate the convergence rate in the approximation of important characteristics associated to fuzzy numbers such as the value, ambiguity and expected interval, respectively. The key element in the process of obtaining the convergence results is the use of the well-known naïve approximator.

5.1 The naïve approximator

In what follows, we propose some convergence results for sequences of piecewise approximations. However, firstly, for given fuzzy number A, \( n\in {\mathbb {N}}\), and \({\varvec{\alpha }}_{n}=(\alpha _{n,0},\alpha _{n,1},\ldots ,\alpha _{n,n+1})\in {\mathfrak {A}}_{n}\), let us consider a fuzzy number \(N_{{\varvec{\alpha }}_{n}}^{n}(A)\in {\mathbb {F}}^{\pi _{n}( {\varvec{\alpha }})}(\mathbb {R)}\) defined by

$$\begin{aligned} \left( N_{{\varvec{\alpha }}_{n}}^{n}(A)\right) _{L}(\beta )&= A_{L}(\alpha _{n,i}) \\&\quad \; +\frac{A_{L}(\alpha _{n,i+1})-A_{L}(\alpha _{n,i})}{\alpha _{n,i+1}-\alpha _{n,i}}(\beta -\alpha _{n,i}), \\ \left( N_{{\varvec{\alpha }}_{n}}^{n}(A)\right) _{U}(\beta )&= A_{U}(\alpha _{n,i}) \\&\quad \; +\frac{A_{U}(\alpha _{n,i+1})-A_{U}(\alpha _{n,i})}{\alpha _{n,i+1}-\alpha _{n,i}}(\beta -\alpha _{n,i}) \end{aligned}$$

for \(i=0,\dots ,n\) and \(\beta \in [\alpha _{i},\alpha _{i+1}]\). This is the so-called naïve approximator of A, i.e., an n-knot piecewise linear fuzzy number that interpolates both \(A_{L}\) and \(A_{U}\) at the knots (it minimizes the \(L_2\) metric exactly at finite number of points, while the approximator suggested in this paper minimizes the metric “globally” , i.e., at all alpha cuts). This operator was also proposed in Báez-Sánchez et al. (2012) where it was used as an approximator with respect to the Hausdorff metric. We will use this operator to obtain approximation results for the sequence \(\left( \Pi _{{\varvec{\alpha }}_{n}}^{n}\right) _{n\ge 1}\) with respect to the Euclidean metric. Please note that the convergence results can be deduced from Báez-Sánchez et al. (2012) using some simple inequalities between the Hausdorff and Euclidean metrics. However, we will use a direct approach which will give us precise estimations for the error of approximation. Clearly, \(N_{{\varvec{\alpha }}}^{n}(A)\) preserves the support and core of A. Since \(N_{{\varvec{\alpha }}_{n}}^{n}(A)\in {\mathbb {F}}^{\pi _{n}( {\varvec{\alpha }})}(\mathbb {R)}\) for any \(n\in {\mathbb {N}}\), let us find an estimation for \(d(N_{{\varvec{\alpha }}_{n}}^{n}(A),A)\) for some \(n\ge 1\). Let us choose arbitrary \(\beta \in [0,1]\). Then there exists \( i\in \{0,\dots ,n\}\) such that \(\beta \in [\alpha _{n,i},\alpha _{n,i+1}]\). If \(\left( N_{{\varvec{\alpha }}_{n}}^{n}(A)\right) _{L}(\beta )\ge A_{L}(\beta )\) then

$$\begin{aligned}&\left| \left( N_{{\varvec{\alpha }}_{n}}^{n}(A)\right) _{L}(\beta ) - A_{L}(\beta )\right| = \left( N_{{\varvec{\alpha }}_{n}}^{n}(A)\right) _{L}(\beta )-A_{L}(\beta ) \\&\quad {\le }\left( N_{{\varvec{\alpha }}_{n}}^{n}(A)\right) _{L}(\alpha _{n,i+1}){-}A_{L}(\alpha _{n,i}) = A_{L}(\alpha _{n,i+1}){-}A_{L}(\alpha _{n,i}){.} \end{aligned}$$

On the other hand, if \(\left( N_{{\varvec{\alpha }}_{n}}^{n}(A)\right) _{L}(\beta )<A_{L}(\beta )\) then

$$\begin{aligned} \left| \left( N_{{\varvec{\alpha }}_{n}}^{n}(A)\right) _{L}(\beta ) - A_{L}(\beta )\right|&= A_{L}(\beta )-\left( N_{{\varvec{\alpha }} _{n}}^{n}(A)\right) _{L}(\beta ) \\&\le A_{L}(\alpha _{n,i+1}){-}\left( N_{{\varvec{\alpha }}_{n}}^{n}(A)\right) _{L}(\alpha _{n,i}) \\&= A_{L}(\alpha _{n,i+1})-A_{L}(\alpha _{n,i}) \end{aligned}$$

and therefore from the two inequalities we get

$$\begin{aligned} \left| \left( N_{{\varvec{\alpha }}_{n}}^{n}(A)\right) _{L}(\beta )-A_{L}(\beta )\right| \le \left| A_{L}(\alpha _{n,i+1})-A_{L}(\alpha _{n,i})\right| . \end{aligned}$$

By similar reasonings, we obtain

$$\begin{aligned} \left| \left( N_{{\varvec{\alpha }}_{n}}^{n}(A)\right) _{U}(\beta )-A_{U}(\beta )\right| \le \left| A_{U}(\alpha _{n,i+1})-A_{U}(\alpha _{n,i})\right| . \end{aligned}$$

This implies that in general, for \(\beta \in [0,1]\), we have

$$\begin{aligned}&\left| \left( N_{{\varvec{\alpha }}_{n}}^{n}(A)\right) _{L}(\beta )-A_{L}(\beta )\right| \nonumber \\&\qquad \qquad \le \max _{i=0,\dots ,n}\left| A_{L}(\alpha _{n,i+1})-A_{L}(\alpha _{n,i})\right| , \end{aligned}$$
(9)
$$\begin{aligned}&\left| \left( N_{{\varvec{\alpha }}_{n}}^{n}(A)\right) _{U}(\beta )-A_{U}(\beta )\right| \nonumber \\&\qquad \qquad \le \max _{i=0,\dots ,n}\left| A_{U}(\alpha _{n,i+1})-A_{U}(\alpha _{n,i})\right| . \end{aligned}$$
(10)

Theorem 4

Let us consider a sequence \(\left( {\varvec{\alpha }}_{n}\right) _{n=1,2,\dots }\), \({\varvec{\alpha }}_{n}\in {\mathfrak {A}}_{n}\), where \({\varvec{\alpha }}_{n}=({\alpha }_{n,i})_{i=0,\dots ,n+1}\) for all n, and such that \(\left\| {\varvec{\alpha }}_{n}\right\| \rightarrow 0\). If A is a fuzzy number with continuous side functions then \(\lim _{n\rightarrow \infty }^{(d)}\Pi _{{\varvec{\alpha }}_{n}}^{n}(A)=A\). Moreover, the rate of convergence is given by the inequalities

$$\begin{aligned} d(A,\Pi _{{\varvec{\alpha }}_{n}}^{n}(A))\le & {} \sqrt{2} \max _{i=0,...,n}(\max \{\left| A_{L}(\alpha _{n,i+1})-A_{L}(\alpha _{n,i})\right| ,\nonumber \\&\left| A_{U}(\alpha _{n,i+1})-A_{U}(\alpha _{n,i})\right| \})\text {.} \end{aligned}$$
(11)

where \(n\ge 1\).

Proof

We will use in our proof the naïve approximator of A, introduced just above. Since functions \(\left( N_{{\varvec{\alpha }}_{n}}^{n}(A)\right) _{L}\), \(\left( N_{{\varvec{\alpha }} _{n}}^{n}(A)\right) _{U}\), \(A_{L}\) and \(A_{U}\) are all continuous, by the mean value theorem for integrals there exists \(\xi _{n}\in [0,1]\) such that

$$\begin{aligned} d^{2}(N_{{\varvec{\alpha }}_{n}}^{n}(A),A) =&\left( \left( N_{{\varvec{\alpha }}_{n}}^{n}(A)\right) _{L}(\xi _{n})-A_{L}(\xi _{n})\right) ^{2} \\&+ \left( \left( N_{{\varvec{\alpha }}_{n}}^{n}(A)\right) _{U}(\xi _{n})-A_{U}(\xi _{n})\right) ^{2}, \end{aligned}$$

which by relations (9)–(10) implies (for any \(n\ge 1\))

$$\begin{aligned} d(N_{{\varvec{\alpha }}_{n}}^{n}(A),A)\le & {} \sqrt{2}\max _{i=0,...,n}(\max \{\left| A_{L}(\alpha _{n,i+1})-A_{L}(\alpha _{n,i})\right| ,\nonumber \\&\left| A_{U}(\alpha _{n,i+1})-A_{U}(\alpha _{n,i})\right| \})\text {.} \end{aligned}$$

Since \(A_{L}\) and \(A_{U}\) are continuous and hence uniformly continuous and \( \left\| {\varvec{\alpha }}_{n}\right\| \rightarrow 0\), from the above inequalities it easily results \(d(N_{{\varvec{\alpha }}_{n}}^{n}(A),A)\rightarrow 0\).

Now let us prove that \(\lim _{n\rightarrow \infty }^{(d)}\Pi _{\varvec{ \alpha }_{n}}^{n}(A)=A\). Since \(N_{{\varvec{\alpha }}_{n}}^{n}(A)\in {\mathbb {F}}^{\pi _{n}({\varvec{\alpha }})}(\mathbb {R)}\) for any \(n\in {\mathbb {N}}\) thus, by the definition of \(\Pi _{{\varvec{\alpha }} _{n}}^{n}(A)\), it follows that

$$\begin{aligned} d(A,\Pi _{{\varvec{\alpha }}_{n}}^{n}(A))\le d(A,N_{{\varvec{\alpha }}_{n}}^{n}(A)) \end{aligned}$$

and as \(d(N_{{\varvec{\alpha }}_{n}}^{n}(A),A)\rightarrow 0\), it is immediate now that \(\lim _{n\rightarrow \infty }^{(d)}\Pi _{\varvec{ \alpha }_{n}}^{n}(A)=A\) with the convergence rate given by (11). \(\square \)

Recall that if \((X,d_{1})\) and \((Y,d_{2})\) are metric spaces then a function \(f:X\rightarrow Y\) is called lipschitzian, if there exists a real positive constant C such that

$$\begin{aligned} d_{2}(f(x^{\prime }),f(x^{\prime \prime }))\le C\,d_{1}(x^{\prime },x^{\prime \prime }), \end{aligned}$$
(12)

for all \(x^{\prime },x^{\prime \prime }\in X\). It is well known that Lipschitz functions are continuous (thus, the property (iv) in Theorem 3 also implies the operator’s continuity).

Proposition 1

Let us consider a sequence \(\left( {\varvec{\alpha }}_{n}\right) _{n=1,2,\dots }\), \({\varvec{\alpha }}_{n}\in {\mathfrak {A}}_n \), where \({\varvec{\alpha }}_{n}=(\varvec{\alpha }_{n,i})_{i=0,...,n+1}\). If A denotes a fuzzy number such that \(A_{L}\) and \(A_{U}\) are Lipschitz functions with Lipschitz constants \(C_{1}\) and \( C_{2}\), respectively, then we have

$$\begin{aligned} d(A,\Pi _{{\varvec{\alpha }}_{n}}^{n}(A))\le C\left\| {\varvec{\alpha }}_{n}\right\| , \end{aligned}$$

where \(C=\sqrt{2}\max \{C_{1},C_{2}\}\) and \(n\in {\mathbb {N}}\).

Proof

The hypothesis implies that

$$\begin{aligned}&|A_{L}(\alpha _{n,i}) -A_{L}(\alpha _{n,i-1})| \le C_{1}\Vert {\varvec{\alpha }}_{n}\Vert \\&|A_{U}(\alpha _{n,i})-A_{U}(\alpha _{n,i-1}) | \le C_{2}\Vert {\varvec{\alpha }}_{n}\Vert . \end{aligned}$$

for \(i\in \{1,\ldots n+1\}\), which leads to the desired conclusion. \(\square \)

5.2 Equidistant knots

From the above theorem, we easily get some important propositions for the case of equidistant knots.

Proposition 2

Let us consider a sequence \(\left( {\varvec{\alpha }}_{n}\right) _{n=1,2,\dots }\), where \({\varvec{\alpha }}_{n}\in {\mathfrak {A}}_n\), \(\alpha _{n,i}=\frac{i}{n+1}\) for all \(i=0,\dots ,n+1\). If A denotes a fuzzy number with continuous sides, then for all \(n\in {\mathbb {N}}\) we have

$$\begin{aligned} d(A,\Pi _{{\varvec{\alpha }}_{n}}^{n}(A))\le & {} \sqrt{2} \max _{i=0,...,n}(\max \{\left| A_L\left( \tfrac{i+1}{n+1}\right) - A_L\left( \tfrac{i}{n+1}\right) \right| ,\nonumber \\&\left| A_U\left( \tfrac{i+1}{n+1}\right) - A_U\left( \tfrac{i}{n+1}\right) \right| \}). \end{aligned}$$

Proof

The proof is immediate by Theorem 4 and noting that \( \alpha _{n,i}=\tfrac{i}{n+1}\) for every \(i\in \{0,1,\ldots n+1\}\). \(\square \)

Example 1

Consider a fuzzy number A with \(\mathrm {supp} \,A=[-5,20]\), \(\mathrm {core}\,A=[3,6]\) and \(\alpha \)-cuts given by

$$\begin{aligned} A_\alpha = \left[ -5+8\,\mathtt {qbeta}(\alpha ; 0.4, 3), 6+14\,(1-\alpha )^4 \right] , \end{aligned}$$

where \(\mathtt {qbeta}(\alpha ; a, b)\) denotes the quantile function of the Beta distribution \(\mathrm {B}(a,b)\), i.e., the inverse of \(F(x)=\int _0^x t^{(a-1)} (1-t)^{(b-1)}\,\mathrm{{d}}t/\beta (a,b)\), where \(\beta (a,b)\) denotes the Euler beta function.

In Fig. 2, we show the distance between A and the naïve approximator, best-Euclidean approximator, and the theoretical upper bound given in Proposition 2, expressed as functions of n. We see that the best-Euclidean one has the fastest convergence rate. It is worth noting that the situation illustrated in Fig. 2 is not unique but it shows a typical behavior observed while studying many numerical examples. \(\boxdot \)

Fig. 2
figure 2

Distance functions from Example 1. Note the logarithmic scale on the y axis

Corollary 1

Let us consider a sequence \(\left( {\varvec{\alpha }}_{n}\right) _{n\ge 1}\), \({\varvec{\alpha }}_{n}\in {\mathfrak {A}}_n\), where \({\varvec{\alpha }}_{n}=\mathbf {(} i/(n+1))_{i=0,\dots , n+1} \). If A denotes a fuzzy number such that \(A_{L}\) and \(A_{U}\) are Lipschitz functions with Lipschitz constants \(C_{1}\) and \( C_{2}\), respectively, then we have

$$\begin{aligned} d(A,\Pi _{{\varvec{\alpha }}_{n}}^{n}(A))\le C\mathbf {/(}n+1),\quad n\in {\mathbb {N}}, \end{aligned}$$

where \(C=\sqrt{2}\max \{C_{1},C_{2}\}\).

Proof

The proof is immediate by applying Proposition 1 and noting that \(\left\| {\varvec{\alpha }}_{n}\right\| =1/(n+1)\) for any \(n\in {\mathbb {N}},n\ge 1\). \(\square \)

Thus, we may conclude that in case of the fuzzy number approximation using equidistant knots the convergence rate is at most linear.

5.3 Convergence w.r.t. some important characteristics

Another important consequence of Theorem 4 is the convergence with respect to some important characteristics of a fuzzy number such as its value, ambiguity, expected interval or expected value. Let us recall briefly their definitions.

The value of a fuzzy number A with respect to a nondecreasing function \(c:[0,1]\rightarrow [0,1]\), called reducing function, where usually \(c(0)=0\) and \(c(1)=1\), is

$$\begin{aligned} \mathrm {Val}_{c}(A)=\int _{0}^{1}c(\beta )\left( A_{L}(\beta )+A_{U}(\beta )\right) \mathrm{{d}}\beta , \end{aligned}$$

and might be seen as a typical value of the magnitude represented by the fuzzy number A. The next index, called the ambiguity, is given by

$$\begin{aligned} \mathrm {Amb}_{c}(A)=\int _{0}^{1}c(\beta )\left( A_{U}(\beta )-A_{L}(\beta )\right) \mathrm{{d}}\beta , \end{aligned}$$

and characterizes the vagueness of A. These two characteristics were proposed for the first time by Delgado et al. (1988) to simplify the representation of fuzzy numbers.

The expected interval of a fuzzy number A is given by (see Dubois and Prade 1987)

$$\begin{aligned} \mathrm {EI}(A)=\left[ \int _{0}^{1}A_{L}(\beta )\mathrm{{d}}\beta ,\int _{0}^{1}A_{U}(\beta )\mathrm{{d}}\beta \right] \end{aligned}$$

while the expected value of A is just the middle of its expected interval, i.e.,

$$\begin{aligned} \mathrm {EV}(A)=\frac{1}{2} \left( \int _{0}^{1}A_{L}(\beta )\mathrm{{d}}\beta +\int _{0}^{1}A_{U}(\beta )\mathrm{{d}}\beta \right) . \end{aligned}$$

It is worth noting that the invariance of the above-mentioned characteristics applies a key role in the fuzzy number approximation (see Grzegorzewski and Mrówka 2005). Hence, the following result is of interest.

Theorem 5

Let us consider a sequence \(\left( {\varvec{\alpha }}_{n}\right) _{n=1,2,\dots }\), \({\varvec{\alpha }}_{n}\in {\mathfrak {A}}_{n}\), where \({\varvec{\alpha }}_{n}=({\alpha }_{n,i})_{i=0,\dots ,n+1}\) for all n, and such that \(\left\| \varvec{ \alpha }_{n}\right\| \rightarrow 0\). If A denotes a fuzzy number with continuous side functions and \(c:[0,1]\rightarrow [0,1]\) is a continuous reducing function, then we have

$$\begin{aligned}&\lim \limits _{n\rightarrow \infty }\mathrm {Val}_{c}(\Pi _{{\varvec{\alpha }} _{n}}^{n}(A)) = \mathrm {Val}_{c}(A), \end{aligned}$$
(13)
$$\begin{aligned}&\lim \limits _{n\rightarrow \infty }\mathrm {Amb}_{c}(\Pi _{{\varvec{\alpha }} _{n}}^{n}(A)) = \mathrm {Amb}_{c}(A), \end{aligned}$$
(14)
$$\begin{aligned}&\lim \limits _{n\rightarrow \infty }\mathrm {EI}(\Pi _{{\varvec{\alpha }} _{n}}^{n}(A)) = \mathrm {EI}(A), \end{aligned}$$
(15)
$$\begin{aligned}&\lim \limits _{n\rightarrow \infty }\mathrm {EV}(\Pi _{{\varvec{\alpha }}_{n}}^{n}(A)) = \mathrm {EV}(A). \end{aligned}$$
(16)

Proof

Noting that \(0\le c\le 1\) and making use of the Schwarz inequality, we get

$$\begin{aligned}&\left| \int _{0}^{1}c(\beta )A_{L}(\beta )\mathrm{{d}}\beta -\int _{0}^{1}c(\beta )\left( \Pi _{{\varvec{\alpha }}_{n}}^{n}(A)\right) _{L}(\beta )\right| \\&\quad \quad \le \int _{0}^{1}\left| A_{L}(\beta )-\left( \Pi _{{\varvec{\alpha }}_{n}}^{n}(A)\right) _{L}(\beta )\right| \mathrm{{d}}\beta \\&\quad \quad \le \left( \int _{0}^{1}\left( A_{L}(\beta )-\left( \Pi _{{\varvec{\alpha }}_{n}}^{n}(A)\right) _{L}(\beta )\right) ^{2}\mathrm{{d}}\beta \right) ^{1/2} \\&\quad \quad \le d(A,\Pi _{{\varvec{\alpha }}_{n}}^{n}(A)). \end{aligned}$$

Repeating the reasoning, we obtain

$$\begin{aligned}&\left| \int _{0}^{1}c(\beta )A_{U}(\beta )\mathrm{{d}}\beta -\int _{0}^{1}c(\beta )\left( \Pi _{{\varvec{\alpha }}_{n}}^{n}(A)\right) _{U}(\beta )\right| \\&\quad \quad \le \left( \int _{0}^{1}\left( A_{U}(\beta )-\left( \Pi _{{\varvec{\alpha }}_{n}}^{n}(A)\right) _{U}(\beta )\right) ^{2}\mathrm{{d}}\beta \right) ^{1/2} \\&\quad \quad \le d(A,\Pi _{{\varvec{\alpha }}_{n}}^{n}(A)). \end{aligned}$$

By the last two inequalities, we easily observe that

$$\begin{aligned}&\max \left\{ \left| \mathrm {Amb}_{c}(A)-\mathrm {Amb}_{c}(\Pi _{{\varvec{\alpha }}_{n}}^{n}(A))\right| , \right. \\&\quad \quad \quad \left. \left| \mathrm {Val}_{c}(A)-\mathrm {Val}_{c}(\Pi _{{\varvec{\alpha }}_{n}}^{n}(A))\right| \right\} \le 2d(A,\Pi _{{\varvec{\alpha }}_{n}}^{n}(A)). \end{aligned}$$

which by Theorem 4 implies that (13) and (14) hold. Moreover, we have

$$\begin{aligned}&\max \left\{ \left| \mathrm {Amb}_{c}(A)-\mathrm {Amb}_{c}(\Pi _{{\varvec{\alpha }}_{n}}^{n}(A))\right| , \right. \\&\quad \quad \quad \left. \left| \mathrm {Val}_{c}(A)-\mathrm {Val}_{c}(\Pi _{{\varvec{\alpha }}_{n}}^{n}(A))\right| \right\} \\&\le 2\sqrt{2}\max _{i=0,...,n}(\max \{\left| A_{L}(\alpha _{n,i+1})-A_{L}(\alpha _{n,i})\right| ,\nonumber \\&\left| A_{U}(\alpha _{n,i+1})-A_{U}(\alpha _{n,i})\right| \}) \end{aligned}$$

for \(n\ge 1\). Hence, we easily obtain the convergence of the expected interval and expected value, i.e., (15) and (16) hold, respectively. \(\square \)

6 Computer implementation and applications

In this section, we propose an algorithm to compute the nearest piecewise linear n-knot approximation. The reasonings are inspired by the particular case studied in Coroianu et al. (2013). Later on, we will show that best-Euclidean piecewise linear approximations are a better alternative than the classical naïve approximation in the implementation of fuzzy arithmetic on a computer.

6.1 Algorithm

Fix \({\varvec{\alpha }}\). To perform the calculations, we should derive formulas for the elements \(\phi _{i,j}=\left\langle {\mathbf {e}}_{i},{\mathbf {e}} _{j}\right\rangle \), \(i,j=1,\dots ,2n+4\), of the Gram matrix \( \Phi =[\phi _{i,j}] \). It might be seen that we have:

$$\begin{aligned} \begin{array}{ll} \phi _{1,1} = 2, \\ \phi _{i,j} = 2-(\alpha _{j-1}+\alpha _{j-2})/2 &{}\quad \text { for }i<j\le n+2, \\ \phi _{j,j} = 2-(2\alpha _{j-1}+\alpha _{j-2})/3 &{}\quad \text { for }j\le n+2, \\ \phi _{i,n+3} = 1 &{}\quad \text { for }i\le n+3, \\ \phi _{i,j} = (\alpha _{2n-j+4}+\alpha _{2n-j+5})/2 &{}\quad \text { for }j\ge n+4, i< j \\ \phi _{j,j} = (2\alpha _{2n-j+4}+\alpha _{2n-j+5})/3 &{}\quad \,\, \text {for}j\ge n+4. \\ \end{array} \end{aligned}$$

The \(\Phi \) matrix is, of course, symmetric and—as already stated—invertible. Note that the elements above the diagonal do not depend on the row index i.

Fix A. Let \({\mathbf {w}}=(w_1,\dots ,w_{2n+3})\), \(\mathbf {w^{\prime }} =(w^{\prime }_0,\dots ,w^{\prime }_{2n+3})\) be such that:

  • \(w^{\prime }_0 = 0\),

  • For \(i=1,\dots ,n+1\) do: \(w_{i} := \int \nolimits _{\alpha _{i-1}}^{\alpha _{i}} A_L(\beta )\,d\beta \),

    $$\begin{aligned} w_{i}^{\prime }:= \frac{\int \nolimits _{\alpha _{i-1}}^{\alpha _{i}} \beta A_L(\beta )\,d\beta - \alpha _{i-1}\,w_{i}}{\alpha _{i}-\alpha _{i-1}}, \end{aligned}$$
  • \(w_{n+2}=w^{\prime }_{n+2}=0\),

  • For \(i=n+3,\dots ,2n+3\) do: \(w_{i} := \int \nolimits _{\alpha _{2n-i+3}}^{\alpha _{2n-i+4}} A_U(\beta )\,d\beta \),

    $$\begin{aligned} w_{i}^{\prime }:= \frac{\alpha _{2n-i+4}\,w_{i} -\int \nolimits _{\alpha _{2n-i+3}}^{\alpha _{2n-i+4}} \beta A_U(\beta )\,d\beta }{\alpha _{2n-i+4}-\alpha _{2n-i+3}}. \end{aligned}$$

It may be shown easily that \({\mathbf {b}}=(b_1,\dots ,b_{2n+4})\) such that \( b_{i}=\left\langle A,{\mathbf {e}}_{i}\right\rangle \), is determined by \(b_{1} := \sum _{i=1}^{2n+3} w_i\) and \(b_{i} := b_{i-1} - w_{i-2}^{\prime }- w_{i-1} + w_{i-1}^{\prime }\) for \(i=2,\dots ,2n+4\).

Moreover, if \(A\ge 0\) (i.e., if \(\inf \mathrm {supp}\,A\ge 0\)), then \(w_i\ge w_i^{\prime }\) and it follows \(b_1\ge b_2\ge \dots \ge b_{2n+4}\ge 0\). Note that \({\mathbf {b}}\) contains the whole information about A needed to solve our approximation problem. Moreover, it may be seen that for \(n=1\) the above derivations for \(\Phi \) and \({\mathbf {b}}\) are equivalent to those presented in our previous paper (Coroianu et al. 2013), and for \(n=0\) with those obtained by Ban (2009a).

Please note that if the solution to the system of \(2n+4\) linear equations \( \Phi \,{\breve{{\varvec{\delta }}}}{}^T={\mathbf {b}}^T\) is such that \(\breve{ \delta }_2,\dots ,\breve{\delta }_{2n+4}\ge 0\), then the problem of determining the nearest piecewise linear approximation is immediate: we have \( {\varvec{\delta }}^*=\breve{{\varvec{\delta }}}\).

Example 2

Consider a fuzzy number A with support [1, 4], core [2, 3], and \(\alpha \)-cuts given by \({A}_\alpha = [1+\sqrt{\alpha }, 4-\alpha ^{2}]\). Let us determine the best approximation of A by a piecewise linear 3-knot fuzzy number with \({\varvec{\alpha }}=(0.25,0.5,0.75)\). We have

$$\begin{aligned} {\mathbf {w}}= & {} \left( \frac{1}{3},\frac{1+\sqrt{2}}{6},\frac{3-2 \sqrt{2}+3\sqrt{3}}{12},\right. \\&\quad \left. \frac{11-3 \sqrt{3}}{12},0,\frac{155}{192}, \frac{173}{192},\frac{185}{192},\frac{191}{192}\right) ,\\ {\mathbf {w}}^{\prime }= & {} \left( 0,\frac{7}{40},\frac{19+4 \sqrt{2}}{120},\frac{15+16 \sqrt{2}-6 \sqrt{3}}{120}, \right. \\&\quad \left. \frac{-11+12 \sqrt{3}}{40}, 0, \frac{317}{768}, \frac{117}{256}, \frac{373}{768},\frac{383}{768}\right) \end{aligned}$$

and thus

$$\begin{aligned} {\mathbf {b}}= & {} \left( \frac{16}{3},\frac{207}{40}, \frac{599-16 \sqrt{2}}{120},\frac{565+16 \sqrt{2}-36 \sqrt{3}}{120},\right. \\&\quad \left. \frac{407+36 \sqrt{3}}{120}, \frac{11}{3}, \frac{2513}{768},\frac{1855}{768},\frac{379}{256},\frac{383}{768})\right) . \end{aligned}$$

Solving \(\Phi \,{{\varvec{\delta }}^*}^T={\mathbf {b}}^T\) we get

$$\begin{aligned} {\varvec{\delta }}^*= & {} \left( \frac{145-84 \sqrt{2}+54 \sqrt{3}}{105},\frac{6}{35} (-2+14 \sqrt{2}-9 \sqrt{3}), \right. \\&\quad -\frac{2}{35} (38+70 \sqrt{2}-81 \sqrt{3}), \frac{80}{7}+4\sqrt{2}-\frac{342 \sqrt{3}}{35},\\&\quad -\frac{2}{35} (328+42 \sqrt{2}-225 \sqrt{3}), \frac{12833+896\sqrt{2}-7488 \sqrt{3}}{1120},\\&\qquad \left. \frac{7}{16}, \frac{5}{16},\frac{3}{16},\frac{1}{16}\right) \ge 0, \end{aligned}$$

and hence we directly obtain the desired approximator.

Here is the R code that calculates this result numerically via the piecewiseLinearApproximation() method from the FuzzyNumbers Gagolewski (2015) package.

figure a

\(\boxdot \)

However, if \(\breve{\delta }_{i}<0\) for some \(i>1\), then we have to find the index set \(K^{*}\subseteq \{2,\dots ,2n+4\}\) corresponding to the optimal solution. Intuitively, this set indicates between which knots (left or right) \(\alpha \)-cut bounds are constant functions. Generally, there are \( 2^{2n+3}\) possible selections of the index sets and we know that at least one of them (as situation with \(\delta _{k}^{*}=z_{k}^{*}=0\) is possible) leads to the solution fulfilling conditions from Theorem 2). Thus, although theoretically correct, in practice we cannot look for each possible K and check whether it gives the desired result. Below we postulate an algorithm that finds the solution in up to \(2n+4\) steps.

First, please note that if \(\mathrm {S}_d({\varvec{\alpha }},\breve{ {\varvec{\delta }}})\not \in {\mathbb {F}}^{\pi _n({\varvec{\alpha }})}({\mathbb {R}})\) then, by definition, we have \(\langle A, {\mathbf {e}}_i\rangle = \langle \mathrm {S}_d({\varvec{\alpha }},\breve{{\varvec{\delta }}}), {\mathbf {e}}_i\rangle = b_i\) for all i. Thus, finding the best linear approximation of A is the same as approximating the object \(\mathrm {S}_d({\varvec{\alpha }} ,\breve{{\varvec{\delta }}})\) (corresponding to a pair of two square integrable piecewise linear functions).

Please note that it may be tempting to assume that as \(n\rightarrow \infty \), then for all \({\varvec{\alpha }}_i\in {\mathfrak {A}}_i\) such that \({\varvec{\alpha }}_i\subset {\varvec{\alpha }}_{i+1}\) we necessarily approach the solution with \(z^*_{j,n}=0\) for all j. This is, unfortunately, not true, as a counterexample may easily be constructed (see Example 3).

It may be noted that if the conditions in Theorem 2 hold for a given index set \(K^*\), then \({\mathbf {z}}^*\) and \({\varvec{\delta }}^*\) may be obtained by solving the following system of \(2n+4\) linear equations for \( {\mathbf {x}}\):

$$\begin{aligned} \left\{ \begin{array}{lll} \sum _{i=1}^{2n+4}\,\left( \phi _{1,i}\,{\mathbf {1}}_{i\not \in K^*} - {\mathbf {1}} _{i=1}\,{\mathbf {1}}_{i\in K^*}\right) \,x_{i} = b_{1} &{}&{}\\ \vdots &{} \vdots &{} \vdots \\ \sum _{i=1}^{2n+4}\,\left( \phi _{2n+4,i}\,{\mathbf {1}}_{i\not \in K^*} - \mathbf { 1}_{i=2n+4}\,{\mathbf {1}}_{i\in K^*}\right) \,x_{i} &{} = &{} b_{2n+4} \\ &{} &{} \end{array} \right. \end{aligned}$$

where \({\mathbf {1}}_p=1\) if p is true and 0 if p is false. This is done by setting \(\delta _i^* = x_i\,{\mathbf {1}}_{i\not \in K^*}\) and \(z_i^* = x_i\, {\mathbf {1}}_{i\in K^*}\). The above system of equations may be written as

$$\begin{aligned} \Phi ^{(K^*)}\,{\mathbf {x}}^T={\mathbf {b}}^T, \end{aligned}$$

where the i-th column of \(\Phi ^{(K^*)}\) is exactly the same as in \(\Phi \) if \(i\not \in K^*\) and set to \((0,\dots ,0,\phi _{i,i}^{(K^*)}=-1,0,\dots ,0)^T\) otherwise. Moreover, \(\Phi ^{(K^*)}\) is always invertible.

The algorithm that finds the solution to the approximation problem of our interest is of “greedy” type. It relies on adding in each step to a temporary K set such index \(i>1\) at which an intermediate solution \(x_i\) has the smallest negative value.

Heuristic Algorithm for determining best piecewise linear approximation of A given n and \({\varvec{\alpha }}\) (cf. Coroianu et al. (2013):

  1. 1.

    Calculate \(\Phi \) and \({\mathbf {b}}\) (according to \(A_L\), \(A_U\), n, and \({\varvec{\alpha }}\));

  2. 2.

    \(K^{(1)}:=\emptyset \);

  3. 3.

    for\(i=1,2,\dots \):

    1. 3.1.

      Solve \(\Phi \,{\mathbf {x}}^T={\mathbf {b}}^T\) for \({\mathbf {x}}\);

    2. 3.2.

      \(m^{(i)} := \min \{ \arg \min _{i=2,3,\dots ,2n+4}\{ x_i \} \}\);

    3. 3.3.

      if\((x_{m^{(i)}} \ge 0)\):

      1. 3.3.1.

        \({\varvec{\delta }}^* := (x_1, x_2\,{\mathbf {1}}_{2\not \in K^{(i)}}, \dots , x_{2n+4}\,{\mathbf {1}}_{2n+4\not \in K^{(i)}})\);

      2. 3.3.2.

        return\(\mathrm {S}_d({\varvec{\alpha }}, \varvec{ \delta }^*)\) as result and stop;

    4. 3.4.

      \(\phi _{i,m}:=0\) for \(i\ne m\);

    5. 3.5.

      \(\phi _{m,m}:=-1\);

    6. 3.6.

      \(K^{(i+1)} := K^{(i)}\cup \{m^{(i)}\}\);

    \(\boxdot \)

Note that \(\mathrm {S}_d({\varvec{\alpha }},{\varvec{\delta }}^*) =\mathrm {S}({\varvec{\alpha }}, \mathtt {cumsum}({\varvec{\delta }}^*))\), where \(\mathtt {cumsum}({\varvec{\delta }}^*) = (\delta _1^*,\delta _1^*+\delta _2^*,\dots , \delta _1^*+\delta _2^*+\dots +\delta _{2n+4}^*)\). Moreover, an example may be constructed easily for which at some step \(m < 0\), \(m > \arg \min _{i=2,3,\dots ,2n+4}\{ x_i \}\) and \(m\not \in K^*\).

Example 3

Consider a fuzzy number A with support [1, 3], core \(\{2\}\), and \(\alpha \)-cuts given by \({A}_\alpha = [1+{\alpha }^{0.2}, 3-\alpha ^{0.2}]\).

Let us determine the best approximation of A by a piecewise linear 3-knot fuzzy number with \({\varvec{\alpha }}=(0.75,0.8,0.9)\). The solution is obtained for \(K=\{3,4,8,9\}\). Here is the output generated by the FuzzyNumbers Gagolewski (2015) package.

figure b

Please, note that adding more knots between 0.75 and 0.9 does not change the resulting piecewise linear fuzzy number (in the sense of \(=\)). \(\boxdot \)

Remark 3

One can easily observe that the finding of the piecewise linear n-knot approximation of a fuzzy number using Theorem 2 directly may be represented as an instance of the Boolean satisfiability problem (SAT) which in general is NP-hard. The proposed heuristic algorithm drastically simplifies the process. Note that if it converges, the obtained solution is guaranteed to be optimal. Whether it always converges is still an open question. However, extensive simulation studies indicate that this is indeed the case.

Also note that from the computer processing perspective, calculating \({\mathbf {b}}\) is the most sensitive part of the procedure in terms of numeric error, due to the fact that \(w_i\) and \(w_i^{\prime }\) have to be obtained by some numerical quadratures, like the adaptive routine integrate() in the R language. This requires the alpha-cut bounds to be well-behaving functions. Interestingly, the trapezoidal rule of integration (the Newton-Cotes formula of degree 1, without subdivisions) will lead us to the naïve piecewise linear approximator (not optimal in general), in which we just probe the alpha-cut bounds at points in \({\varvec{\alpha }}\).

6.2 Computing on piecewise linear FNs

Suppose we have a set of fuzzy numbers \(\{A_1,\ldots ,A_k\}\) and we would like to compute a series of operations on them, obtaining \( B=f(A_1,\ldots ,A_k)\). Let each operation be defined using the extension principle and rely on a proper transformations of their \(\alpha \) -cuts.

Most often such a task is performed numerically and not symbolically. Thus, the side functions of the fuzzy numbers must be discretized at a fixed, possibly large number of \(\alpha \)-cuts (see, e.g., Hanss 2005). Such a process involves nothing else than taking a naïve approximation of \(A_i\) , \(i=1,\dots ,k\). Here, the values of membership functions at given \(\alpha \) -cuts (e.g., equidistant ones) are exact (of course, up to numeric error involved in calculating the operations). Of course, the larger the number of knots, the lower the computation speed but higher the accuracy.

In practice, we are interested in finding a trade-off between these two factors. An important question is whether, for fixed n and \({\varvec{\alpha }}\in {\mathfrak {A}}_n\), by considering the nearest-Euclidean approximation we will obtain results of higher quality than in case of a naïve approximation.

a) Experiment setup:

Let \(N=N_{{\varvec{\alpha }}}^n(f(N_{{\varvec{\alpha }}}^n(A_1),\dots ,N_{{\varvec{\alpha }}}^n(A_k)))\) and \(\Pi = \Pi _{{\varvec{\alpha }}}^n(f(\Pi _{{\varvec{\alpha }}}^n(A_1),\dots ,\Pi _{{\varvec{\alpha }}}^n(A_k)))\) denote piecewise linear n-knot fuzzy numbers determined by taking, respectively, the naïve and nearest-Euclidean approximations of the input fuzzy numbers and then applying f only at given \({\varvec{\alpha }}\)-cuts. Moreover, let \(\Pi \text {-post}=\Pi _{{\varvec{\alpha }}}^n(B)\) denote the nearest-Euclidean piecewise linear approximation of the exact result, B. To verify this, let us perform a series of numerical experiments. First of all, one should be interested in:

  • Comparing d(BN) and \(d(B, \Pi )\) where d denotes the Euclidean metric, i.e., determining whether it is better to use the nearest-Euclidean or the naïve approximation,

  • Calculating \(d(B, \Pi \text {-post})\) and \(d(\Pi , \Pi \text {-post})\), to determine a possible “error” which might appear if we perform computations directly on approximated fuzzy numbers.

To investigate these two issues, we performed a large number of computer simulations involving the following randomly generated types of fuzzy numbers:

  • FNs with sides being power functions, i.e., \(\alpha \mapsto \alpha ^p\), where p is uniformly distributed on the interval (0, 10), i.e., \(p\sim U(0,10)\);

  • FNs with sides defined via quantile functions of a beta distribution, \( \alpha \mapsto \mathtt {qbeta}(\alpha , p, q)\), \(p,q\sim U(0.1,4)\).

Supports \([a_1,a_4]\) and cores \([a_2,a_3]\) were also generated by \((a_1,a_2,a_3,a_4):=\mathrm {sort}(a_1^{\prime },a_2^{\prime },a_3^{\prime },a_4^{\prime })\), with \(a_1^{\prime },a_2^{\prime },\)\(a_3^{\prime },a_4^{\prime }\) following the uniform distribution.

Fig. 3
figure 3

Results of the numerical experiment: \(d(B, \Pi )\), d(BN), \(d(B, \Pi \text {-post})\) (left subfigure) and \(d(\Pi , \Pi \text {-post})\) (right subfigure)

b) Distance between exact results and the results computed on approximated fuzzy numbers:

As all the considered scenarios exhibited the same patterns, below we will only present results for the scenario with \(f(A_1,\dots ,A_9)=(A_1+A_2)A_3-A_4(A_5-A_6)\log (A_7)+\exp (A_8)/2^{A_9}\), with \(M=1000\) Monte Carlo iterations, and \(n=3\) equidistant knots. Figure 3 depicts the boxplots of \(d(B, \Pi )\), d(BN) and \(d(B,\Pi \text {-post})\), as well as \(d(\Pi , \Pi \text {-post})\). It is worth noting that \(d(B,\Pi )<d(B,N)\) in each case. In other words, applying computations on nearest-Euclidean approximated inputs leads to smaller errors (as measured by the Euclidean distance) than for the naïve approximations. Figure 4 depicts the difference between an exemplary exact output and ones that are computed using two types of piecewise linear approximations.

Fig. 4
figure 4

An exemplary result of applying \(f(A_1,\dots ,A_9) =(A_1+A_2)A_3-A_4(A_5-A_6)\log (A_7)+\exp (A_8)/2^{A_9}\): the exact one, and ones computed on best-Euclidean and naïvely approximated inputs

Fig. 5
figure 5

Scatter plots for \(d(B, \Pi )\), d(BN), and \(d(B, \Pi \text {-post} ) \)

Next, Fig. 5 presents scatter plots for each pair of observations from \(\{d(B, \Pi ), d(B, N), d(B, \Pi \text {-post})\}\). Interestingly, we observe high linear association between the errors (Pearson’s \(r>0.99\)). Applying linear regression, we get the following relationships for \(n=3\):

$$\begin{aligned} d(B, \Pi )\simeq & {} 1.00544\, d(B, \Pi \text {-post}) \\ d(B, N)\simeq & {} 3.31601\, d(B, \Pi \text {-post}) \\ d(B, \Pi )\simeq & {} 0.30237\, d(B, N) \end{aligned}$$

which are quite nicely fitted (\(R^2>0.99\)). Moreover, for \(n=10\) equidistant knots we have

$$\begin{aligned} d(B, \Pi )\simeq & {} 1.00474\, d(B, \Pi \text {-post}) \\ d(B, N)\simeq & {} 2.29118\, d(B, \Pi \text {-post}) \\ d(B, \Pi )\simeq & {} 0.43851\, d(B, N). \end{aligned}$$

We see that the error induced by the nearest-Euclidean approximation is more than 2–3 times smaller than that accompanying the naïve approximation.

Moreover, in Table 1 we list basic summary statistics of the \(L_2\) error measure between the exact solution and ones obtained by means on acting on approximated fuzzy numbers. We see that with a large number of knots we get much more precise results and that the best \(L_2\) approximator is far better than the naïve one.

Table 1 Basic summary statistics of the error measure between the solution to \(f(A_1,\dots ,A_9)=(A_1+A_2)A_3-A_4(A_5-A_6)\log (A_7)+\exp (A_8)/2^{A_9}\) and one obtained by the two studied approximators for different n (0 – trapezoidal, 1, 3, and 10)
Table 2 Basic characteristics of the empirical distribution of \(F_1=d_M(\mathrm {supp}\,B,\mathrm {supp}\,\Pi )\), \(F_2=d_M(\mathrm {supp}\,B,\mathrm {supp}\,\Pi \text {-post})\), \(F_3=d_M(\mathrm {core}\,B,\mathrm {core}\,\Pi )\) and \(F_4=d_M(\mathrm {core} \,B,\mathrm {core}\,\Pi \text {-post})\)

c) Preservation of fuzzy numbers’ characteristics:

We might be also interested whether the nearest-Euclidean approximation preserves better some important characteristics of the concerned fuzzy numbers. Of course, the support and core of B and N are the same. To measure the error for \(\Pi \) and \(\Pi \text {-post}\), we will use the Moore’s interval metric Moore (1962), given by \(d_M([a,b],[c,d])=\max \{|a-c|,|b-d|\}\).

Table 2 shows some basic summary statistics for those measures of interest. We observe again the almost perfect correlation (\(r>0.999\)) between \(d_M(\mathrm {supp}\,B\), \(\mathrm {supp}\,\Pi )\) and \(d_M(\mathrm {supp}\,B,\mathrm {supp}\,\Pi \text {-post})\) as well as between \(d_M(\mathrm {core}\,B,\mathrm {core}\,\Pi )\) and \(d_M(\mathrm {core}\,B,\mathrm {core}\,\Pi \text {-post})\). Indeed

$$\begin{aligned} d_M(\mathrm {supp}\,B,\mathrm {supp}\,\Pi ) \simeq&1.000676\, d_M(\mathrm {supp}\,B,\mathrm {supp}\,\Pi \text {-post}) \\ d_M(\mathrm {core}\,B,\mathrm {core}\,\Pi ) \simeq&1.007172\, d_M(\mathrm {core}\,B,\mathrm {core}\,\Pi \text {-post}). \end{aligned}$$

The above-mentioned results are very interesting especially that the stability of the support and core is very important in fuzzy number approximation (for more details we refer the reader to Coroianu et al. 2014a).

Lastly, let us investigate the difference between B’s ambiguity, width, expected value, or value and the corresponding characteristics of \(\Pi , N\), and \(\Pi \text {-post}\). The relevant boxplots are depicted in Figs. 6 and 7. We see that the nearest-Euclidean approximation preserves these fuzzy numbers characteristics much better than the naïve approximation. Although these characteristics are not ideally preserved, they are reconstructed quite accurately.

Fig. 6
figure 6

\(\mathrm {amb}\,B-\mathrm {amb}\,\Pi \), \(\mathrm {amb}\,B-\mathrm {amb} \,N\), and \(\mathrm {amb}\,B-\mathrm {amb}\,\Pi \text {-post}\) (left subfigure) as well as \(\mathrm {width}\,B-\mathrm {width}\,\Pi \), \(\mathrm {width}\,B- \mathrm {width}\,N\), and \(\mathrm {width}\,B-\mathrm {width}\,\Pi \text {-post}\) (right subfigure)

Fig. 7
figure 7

\(\mathrm {val}\,B-\mathrm {val}\,\Pi \), \(\mathrm {val}\,B-\mathrm {val} \,N\), and \(\mathrm {val}\,B-\mathrm {val}\,\Pi \text {-post}\) (left subfigure) as well as \(\mathrm {ev}\,B-\mathrm {ev}\,\Pi \), \(\mathrm {ev}\,B-\mathrm {ev}\,N\) , and \(\mathrm {ev}\,B-\mathrm {ev}\,\Pi \text {-post}\) (right subfigure)

7 Conclusions

The problem of fuzzy number approximation by the piecewise linear fuzzy numbers introduced in Báez-Sánchez et al. (2012) was considered. In this paper, the nearest piecewise linear approximation with respect to the Euclidean metric was concerned. The properties of the approximator, including the asymptotic ones, were investigated. The practical implementation of the approximation algorithm is available in the FuzzyNumbers package for R by Gagolewski (2015) .

The results on convergence indicate some advantages of the piecewise approximation. Let us recall that the general idea is not only to approximate an arbitrary fuzzy number by another fuzzy number with a simpler representation, but to find such an approximation which also possesses some interesting properties. For example, in the papers Ban (2008), Coroianu (2011), Grzegorzewski and Mrówka (2005, 2007), Yeh (2008a), the trapezoidal approximation preserving the expected interval is studied. Although this type of approximation is important in some applications, the expected interval invariance may imply that other important characteristics, such as the value or ambiguity, are not generally preserved there (see, e.g., Ban 2008). Additionally, in some cases as a result of the trapezoidal approximation we obtain a degenerated triangular fuzzy number such that the output support and core is relatively far from the input support and core, respectively. Similarly, the trapezoidal approximation preserving the ambiguity (see Ban and Coroianu 2012), trapezoidal approximation preserving the ambiguity and value (see Ban et al. 2011a) or the trapezoidal approximation preserving the core (see Abbasbandy and Hajjari 2009b), entails the loss of invariance of some other important characteristics of fuzzy numbers.

When discussing about our best-Euclidean piecewise linear approximator, one may ask naturally whether it is really better than the naïve approximator, i.e., an n-knot piecewise linear fuzzy number that interpolates both \(A_{L}\) and \(A_{U}\) at the knots. Although, as usually in such cases, there is no overall champion, it is worth stressing that

  • In the case of the best-Euclidean approximator, we have a guarantee that the solution is the closest possible in terms of the \(L_2\) metric, while there is no such a guarantee in the naïve case.

  • The best-Euclidean approximator has faster rate of convergence than the naïve one. Convergence rate is determined w.r.t. \(L_2\) metric, thus it is obvious that the suggested approximator will be better in this case.

  • One can easily observe that the finding of the piecewise linear n-knot approximation of a fuzzy number using Theorem 2 directly may be represented as an instance of the Boolean satisfiability problem (SAT) which in general is NP-hard. The proposed heuristic algorithm drastically simplifies the process. Note that if it converges, the obtained solution is guaranteed to be optimal. Whether it always converges is still an open question. However, extensive simulation studies indicate that this is indeed the case.

  • The algorithm proposed in Sect. 6 makes the computations as easy as in the naïve case (from the practitioner’s perspective). However, practical computation on the best-Euclidean approximations are more exact (see simulation study), especially if the series of arithmetic operations are involved, while the naïve approximators lead to a great error propagation.

Many interesting problems related to the approximation method suggested in this paper are still open. Firstly, a study of the global uniform convergence of the functions \(\Pi _{{\alpha _{\mathbf{n }}}}(A)_{L}\) and \(\Pi _{{\alpha _{\mathbf{n }}}}(A)_{U}\) which will imply the convergence with respect to the support and core is strongly desired.

Next important problem is a comparison between different methods that one can apply to approximate fuzzy numbers. In particular, it would be interesting to compare our general piecewise linear approximation and the approximation of fuzzy numbers by using the F-transform (see Coroianu and Stefanini 2016; Stefanini and Sorini 2012) or with approximation by the Bernstein operators of max-product kind (see Coroianu et al. 2014b). Here we can discuss many issues as, for instance, which method is easier to implement on the computer. Then, of course, we could compare their convergence rates (with respect to different kind of metrics such as the Chebyshev or the Euclidean), and their errors in approximating the important characteristics.

Finally, as it was mentioned above, the convergence of the algorithm proposed in Sect. 6 is still an open question.