# Extended Fully Fuzzy Linear Regression to Analyze a Solid Cantilever Beam Moment

Authors: Seyedehnegar Seyedmonir; Mostafa Bayrami; Saeid Jafarzadeh Ghoushchi; Amir Alipour Yengejeh; Hakimeh Morabbi Heravi
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

There are several procedures such as possibilistic and least-square methods to estimate regression models. In this study, first, a fully fuzzy regression equation is converted into a fully fuzzy linear framework. By considering a least-square approach, a model is suggested based on matrix equations for solving fully fuzzy regression models. The main advantage of this method over existing ones is that this method considered values based on their specification, and all linear problems can be easily solved. Moreover, a case study for solid mechanics about the quantity of beam momentum is considered. In this example, the inner data are force values, and the output is momentum values.


---

## Body

## 1. Introduction

Regression assessment is concerned with statistical methodologies’ collection for simulation, solving, and researching the correlation between response variables and regressor or predictor ones. The examples of regression analysis applications are various and occur in almost all applied fields consisting of engineering, chemical, physics, biology, social science, management, and economics. Among various regression models, linear one is so popular. Additionally, the main nature of regression models is the estimation of the parameters employing possibility and least square approaches. However, the classic regression may sometimes fail to cover all observation datasets, particularly those which are imprecise or vague. In other words, some observations are considered in the fuzzy form [1–4]. Therefore, how to estimate regression coefficients and make the subsequent prediction under a classic environment is a dominant challenge to the classical regression analysis.Tanaka [5] was probably a pioneer who proposed a solution for regression problems with the fuzzy dependent variable and crisp independent one in 1982. Afterward, some investigations have been published by different authors concentrated on various variants of fuzzy regression, their characteristics, and applications, namely, Savic and Pedrycz [6]; Kacprzyk [7]; Chang and Lee [8]; Peters [9]; Nather [10]; Kim and Bishu [11]; Diamond and Korner [12]; Ming et al [13]; Tanaka et al. [14]; Ezzati et al. [15]; Wang and Tsaur [16]; D’Urso and Gastaldi [17]; Hong et al. [18]; Modarres et al. [19];Wang et al. [20]; Mosleh et al. [21]; Hassanpour et al. [22]; Nong [23]; Chung [24]; Mosleh et al. [25]; Ezzati et al. [26]; Ebrahimnejad and Tavana [27]; and Nazari-Shirkouhi and Keramati [28]. A few studies have focused on the situations that consider both input and output variables in the fuzzy form. For instance, Sakawa and Yano [29] discussed it for the first time; then, Diamond [30] obtained the regression coefficients in the fuzzy environment using a method called shape-preserving operations, while Kao and Chyu [31] did it by introducing another method named two-stage methodology. Furthermore, D’Urso [32] proposed an iterative least square approach to study all possible combinations of fuzzy/crisp input-fuzzy/crisp output in detail. Finally, Mosleh et al. [25] applied the neural network method to evaluate fully fuzzy regression models.The estimation of regression parameters in linear models leads to a series of equations in that solving is somewhat difficult in fuzzy forms and especially fully fuzzy states. Therefore, several methods have been proposed by some researchers such as Ma et al. [33]; Wang et al. [34]; and Dehghan and Hashemi [35] for treating fuzzy linear systems. Moreover, a convenient way has recently been suggested by Kumar et al. [36], and its updated version in matrix equations was proposed by Otadi and Mosleh [37] for solving the fully fuzzy matrix equations (FFMEs) with triangular fuzzy numbers/arrays. A fuzzy semiparametric regression model was introduced by Hesamian and Akbari [38]. They employed fuzzy predictors, a smooth function, responses, and also exact coefficients. The suggested methodology was based on novel kernel-based signed-distance measures in the space of fuzzy numbers. They demonstrated that the defined signed-distance measure was more reasonable and productive than a common signed-distance which was extensively used in fuzzy set theory. In another study [39], these researchers proposed a new fuzzy additive regression model with nonfuzzy predictors and fuzzy responses. They defined a backfitting stepwise regression approach with kernel smoothing to predict a fuzzy smooth function corresponding to each predictor. Hose and Hanss [40] proposed a fuzzy linear least square for the identification of possibilistic regression models. They reported that their technique applied to a wide range of problems that could be formulated as a linear least-squares problem.In this study, fully fuzzy linear regression matrix equations obtained by the least-square approach through the latter one are solved. It should be noted that all computations of this article are performed by RStudio and LINGO software. This study consists of some sections as follows. In Section2, preparatives of fuzzy backgrounds have been presented. The FFME method is explained in Section 3. The general form of fully fuzzy linear regression (FFLR) is considered, and a new simple method for estimating its parameters is argued in Section 4. Finally, a numerical example is presented to reveal the application ideas of this article in Section 5.


## 2. Preparatives

In this section, some initial definitions utilized in fuzzy operations are reviewed. The fuzzy set theory was first introduced by Prof. Zadeh [41]. Then, many scholars from around the world used it in different fields of study. A triangular fuzzy number (TFN) is known as a type of trapezoidal fuzzy number wide-spread employed in the literature. The mathematical membership function μv˜:U⟶0,1 that is defined for the TFN v˜=α,β,γ is given as follows [42] (refer to Figure 1).Figure 1
Triangular fuzzy numbers.Definition 1.
v˜=α,β,γ is a triangular number where β is the center point, α is the left width, and β is the right one. The membership function of the above triangular fuzzy number is as follows:(1)μv˜x=x−αβ−α,α≤x≤βx−γβ−γ,β≤x≤γ0,otherwise.Definition 2.
The triangular fuzzy numberα,β,γ is nonnegative if and only if α≥0.Definition 3.
The equality of two triangular fuzzy numbers ofu˜=α,β,γ and v˜=α′,β′,γ′ is realized if and only if α=α′, β=β′, and γ=γ′.Definition 4.
In the case that all the components of the matrixM˜=m˜ij are fuzzy numbers, M˜ is known as a fuzzy number matrix. Moreover, according to the above definitions, if each fuzzy element M˜ is positive (negative), then M˜ is regarded as a fuzzy positive (negative) matrix M˜≻0, (M˜≺0). The same definitions can be expressed for nonnegative and nonpositive fuzzy matrices.Definition 5.
Supposem˜=α,β,γ and n˜=α′,β′,γ′ as two triangular fuzzy numbers. Therefore,(i)
m˜⊕n˜=α+α′,β+β′,γ+γ′.(ii)
−m˜=−γ,−β,−α.(iii)
m˜Θn˜=α−γ′,β−β′,γ−β′.As to fuzzy multiplication, it is obvious that fuzzy multiplication will result in a fuzzy number. Here, by delivering some special privileges because of the augmentation guideline, it is a bit different from the classical fuzzy one [37]. Therefore, all fuzzy multiplications and divisions in this study are denoted via term(2)t˜⊗s˜=ν,θ,λ,with(3)θ=β.β′,ν≔minα.α′,α.γ′,γ.α′,γ.γ′,λ≔maxα.α′,α.γ′,γ.α′,γ.γ′.Ift˜ and s˜ are triangular fuzzy numbers, by considering s˜ as a nonnegative fuzzy number,(iv)
t˜⊗s˜=α.α′,ββ′,γ.γ′,α≥0,α.γ′,ββ′,γ.γ′,α<0,γ≥0,α.γ′,ββ′,γ.α′,α<0,γ<0.(v)
m˜⊘n˜=α/γ′,β/β′,γ/α′.


## 3. General Fully Fuzzy Matrix Equation

The following equation is called anm×n fully fuzzy linear system (FFLS) in which r˜ij,1≤i,j≤n, and y˜i,1≤i≤n, are arbitrary triangular fuzzy numbers.(4)r˜11⊗x˜1⊕r˜12⊗x˜2⊕⋯⊕r˜1n⊗x˜n=y˜1r˜21⊗x˜1⊕r˜22⊗x˜2⊕⋯⊕r˜2n⊗x˜n=y˜2⋮⋮⋮⋮r˜m1⊗x˜1⊕r˜m2⊗x˜2⊕⋯⊕r˜mn⊗x˜n=y˜n.Also, a matrix notation of the above system can be written as(5)r˜11r˜12…r˜1nr˜21r˜22…r˜2n⋮⋮⋮⋮r˜m1r˜n2…r˜mnu˜1u˜2⋮u˜n=y˜1y˜2⋮y˜n,where the unknown elements u˜ given by u˜i=vi,ui,wi are called a solution of the fuzzy matrix system if(6)R˜⊗u˜=y˜.For solving (6), a new notation R˜=F,A,E has been suggested (6) where F, A, and E are three crisp matrices with the equivalent size of R˜ as well as(7)b˜=q1,b1,p1,q2,b2,p2,…,qn,bn,pn'.So, (6) can be rewritten as(8)Fv=qAu=bEw=p.u˜ can be both a coherent and constant solution if v≥0,u−v≥0, and w−u≥0.Even based on (8), we have(9)F,A,Ev,u,w=q,b,p.Moreover, FFME (8) can be written in the following terms if it is supposed that fik,aik,eik⊗vk,uk,wk=mikj,tikj,nikj,1≤i,j,k≤n, concerning(10)∑i=1nmikj,tikj,sikj=qk,bk,pk,1≤k≤n.Now, by adding artificial variables, namely,zl,l=1,2,…,3n2, the following linear programming can result.Minimizez1+z2+⋯+z3n2,(11)subject to∑k=1nm1k1+z1=q1,∑k=1nm2k1+z2=q2,⋮∑k=1nmnk1+zn=qn,∑k=1nt1k1+zn+1=b1,⋮∑k=1nsnkn+z3n2=pn,vk,wk,uk≥0,uk−vk≥0,zl≥0,1≤i,j,k≤n,1≤l≤3n2.


## 4. Fully Fuzzy Linear Regression

A fuzzy linear regression usual configuration is described as follows:(12)y˜=a˜0⊕a˜1⊗x˜1⊕a˜2⊗x˜2⊕,…,⊕a˜n⊗x˜n,where x˜i,1≤i≤n, are fuzzy independent variables (input), a˜i,1≤i≤n, are regression coefficients (unknown parameters), and y˜i,1≤i≤n, are fuzzy dependent variables (output).The fully fuzzy regression technique’s key objective is to achieve a model including the least prediction error which has an acceptable agreement with observations. In the current fuzzy least-square approach, the FFLR is normally gained by computing the regression parametersa˜i,1≤i≤n, based on an observed dataset X˜i,y˜i, where X˜i=1˜,x˜1,x˜2,…,x˜nt; likewise, X˜i⊗a˜=y˜i.(13)y˜i=a˜0⊕a˜i⊗x˜i1⊕a˜2⊗x˜i2⊕,…,⊕a˜n⊗x˜in,i=1,2,…,m,where a˜=a˜1,a˜2,…,a˜nt. In this work, the least-square regression configuration is considered to minimize errors of the predicted values based on their features. This methodology of fuzzy regression is an augmentation of the traditional least squares for limiting the distance of the perceptions of the info yield informational index. By gathering the above m equations, in the matrix form, the accompanying FFLS is acquired:(14)y˜=X˜⊗a˜,where X˜=X˜1,X˜2,…,X˜mt, y˜=y˜1,y˜2,…,y˜mt, and a˜=a˜1,a˜2,…,a˜nt.The following theorem states how a meaningful prediction for the solution of (15) can be computed.Theorem 1.
LetX be a nonnegative fuzzy matrix; a˜ is the solution of (15) if(15)X˜T⊗X˜⊗a˜=X˜T⊗y˜.Owing to this theorem, a fully fuzzy regression analysis is converted into a fully fuzzy linear system for computinga˜. Therefore, a system of FFMEs is solved based on LP instead of solving directly. It should be noted that we consider X˜T⊗X˜=A˜ and X˜T⊗y˜=b˜.


## 5. Numerical Examples

To illustrate the proposed approach, two examples are introduced as follows.


### 5.1. Example 1

As the first example, the estimation fory=Lnx is represented. The variables are gathered in Table 1.Table 1
Observations for variables of the example.
No.X˜≡v,x,z independent variableY˜≡g,y,h response variable11.00,1.10,1.400.00,0.10,0.3421.10,1.30,1.700.10,0.26,0.5331.20,1.60,2.000.18,0.47,0.6941.30,1.90,2.300.26,0.64,0.8351.40,2.20,2.600.34,0.79,0.9661.50,2.50,2.900.41,0.92,1.06As shown in Table1, all observations for two variables are considered triangular fuzzy members. Therefore, an FFLR model will be appropriate for showing the relationship between them:(16)y˜=a˜0⊕a˜1⊗x˜,with unknown fuzzy parameters a˜0 and a˜1. Hence, X˜T⊗X˜⊗a˜=X˜T⊗y˜, and we have(17)6.00,6.00,6.007.50,10.60,12.907.50,10.60,12.909.55,20.16,29.31a˜0a˜1=1.28,3.17,4.411.74,6.44,10.24.Leta˜0=c0,a0,b0anda˜1=c1,a1,b1 be nonnegative triangular fuzzy numbers. Thus, the above FMME is rewritten as(18)6.00,6.00,6.00⊗c0,a0,b0⊕7.50,10.60,12.90⊗c1,a1,b1=1.28,3.17,4.41,7.50,10.60,12.90⊗c0,a0,b0⊕9.55,20.16,29.31⊗c1,a1,b1=1.74,6.44,10.24.By applying fuzzy addition and multiplication operations introduced in Section2,(19)6.00c0+7.50c1,6.00a0+10.60a1,6.00b0+12.90b1=1.28,3.17,4.41,7.50c0+9.55c1,10.60a0+20.16a1,12.90b0+29.31b1=1.74,6.44,10.24.According to the suggested method, instead of FMME, we have the following crisp system:(20)6.00c0+7.50c1=1.28,6.00a0+10.60a1=3.17,6.00b0+12.90b1=4.41,7.50c0+9.55c1=1.74,10.60a0+20.90a1=6.44,12.90b0+29.31b1=10.24,Based on LP, we have the following.Minimizez1+z2+,⋯,+z6,(21)subject to6.00c0+7.50c1+z1=1.28,6.00a0+10.60a1+z2=3.17,6.00b0+12.90b1+z3=4.41,7.50c0+9.55c1+z4=1.74,10.60a0+20.90a1+z5=6.44,12.90b0+29.31b1+z6=10.24,a0−c0≥0,a1−c1≥0,b0−a0≥0,b1−a1≥0,ci≥0,ai≥0,bi≥0,zj≥0,i=0,1,j=1,2,…,6.The optimal solution to the above problem is(22)c0=0.00,a0=0.00,c1=0.00,b0=0.17,a1=0.30,b1=0.34.Also, fuzzy solution problem (15) is the regression equation:(23)a˜0=0.00,0.00,0.00,a˜1=0.17,0.30,0.34,y˜=0.00,0.00,0.00⊕0.17,0.30,0.34⊗x˜.Figure2 represents the response of this example in a range that includes all uncertainties.Figure 2
Chart of response variables vs. independent variable for example 1.The figure consists of three components including lower, middle, and upper linear triple functionsY (lower), Y (middle), and Y (upper), respectively. The figure shows that Y (lower) has a considerable deviation from the main function, while Y (middle) and Y (upper) are suitable approximations for the Ln function, especially after x=1.25.


### 5.2. Example 2

This section includes a case study from the solid mechanics field to show an application of this method. All computations have been performed by RStudio and LINGO software. Figure3 represents the schematic form of the considered cantilever beam. The main output of this system is the value of the moment. Based on statics [34], it can be written as follows:(24)M⇀=r⇀×F⇀,where M⇀ is the moment vector, r⇀ is the distance vector from the origin in which moment acts around to the effect point of force, and finally, F⇀ is the force vector that causes the moment. As it is obvious in Figure 3, two types of forces can be considered including F that is an external force and W which is the weight of the beam that is applied to the center of beam mass. The moment caused by weight can be considered as a constant value, but the moment caused by an external force can vary since its value and effect point can be different. Also, if the distance vector is considered perpendicular to the force vector, the equation can be shown in the scalar form. Therefore, equation (24) can be rewritten as follows:Figure 3
Schematic form of the considered cantilever beam.(25)M=MF+Mw.By substituting the moment components in equation (25), it will be(26)M=RF+rw⟶M=RF+Mo.If we considerF as X˜ and M as Y˜, R and Mo will be a˜1 and a˜0, respectively. Table 2 shows the variables of an example in this area. In this example, eight components are considered as the independent variables, and consequently, eight response variables will be achieved.Table 2
Observations for variables of the example.
No.X˜≡v,x,z independent variableY˜≡g,y,h response variable1(1.5, 2, 2.5)(3.5, 4, 4.5)2(3, 3.5, 4)(5, 5.5, 6)3(4.5, 5.5, 6.5)(6.5, 7.5, 8.5)4(6.5, 7, 7.5)(6, 6.5, 7)5(8, 8.5, 9)(8, 8.5, 9)6(9.5, 10.5, 11.5)(7, 8, 9)7(10.5, 11, 11.5)(10, 10.5, 11)8(12, 12.5, 13)(9, 9.5, 10)As shown in Table2, all observations for two variables are considered triangular fuzzy members. Therefore, an FFLR model will be appropriate for showing the relationship between them:(27)y˜=a˜0⊕a˜1⊗x˜,with unknown fuzzy parameters a˜0 and a˜1. Hence, X˜T⊗X˜⊗a˜=X˜T⊗y˜, and we have(28)8,8,855.5,60.5,65.555.5,60.5,65.5482.25,555.25,635.25a˜0a˜1=55,60,65432,504.5,584.Leta˜0=c0,a0,b0anda˜1=c1,a1,b1 be nonnegative triangular fuzzy numbers. Thus, the above FMME is rewritten as(29)8,8,8⊗c0,a0,b0⊕55.5,60.5,65.5⊗c1,a1,b1=55,60,6555.5,60.5,65.5⊗c0,a0,b0⊕482.25,555.25,635.25⊗c1,a1,b1=432,504.5,584.By applying fuzzy addition and multiplication operations introduced in Section2,(30)8c0+55.5c1,8a0+60.5a1,8b0+65.5b1=55,60,6555.5c0+482.25c1,60.5a0+555.25a1,65.5b0+635.25b1=432,504.5,584.According to the suggested method, instead of FMME, we have the following crisp system:(31)8c0+55.5c1=55,8a0+60.5a1=60,8b0+65.5b1=65,55.5c0+482.25c1=432,60.5a0+555.25a1=504.5,65.5b0+635.25b1=584,and by using LP.Minimizez1+z2+,⋯,+z6,(32)subject  to8c0+55.5c1+z=551,8a0+60.5a1+z2=60,8b0+65.5b1+z3=655,55.5c0+482.25c1+z4=432,60.5a0+555.25a1+z5=504.5,65.5b0+635.25b1+z6=584,a0−c0≥0,a1−c1≥0,ci≥0,ai≥0,bi≥0,zj≥0,i=0,1,j=1,2,…,6.The optimal solution to the above problem is(33)c0=0.876,a0=0.876,b0=1.030,c1=0.795,a1=0.813,b1=0.813.Also, fuzzy solution problem (15) is the regression equation:(34)a˜0=0.876,0.876,1.030,a˜1=0.795,0.813,0.813,y˜=0.876,0.876,1.030⊕0.795,0.813,0.813⊗x˜.These results are shown in Figure4. This figure represents the response of the case study in a range that includes all uncertainties.Figure 4
Chart of response variables vs. independent variable for the cantilever beam.As a sample, ifX˜1=7.75,8.25,9.25 and X˜2=14.5,15,15.5, based on regression equation achievement, we have Y˜1=7.037,7.583,8.550 and Y˜2=12.404,13.071,13.632, respectively. Figure 4 consists of three components including lower, middle, and upper linear triple functions Y (lower), Y (middle), and Y (upper), respectively.


## 5.1. Example 1

As the first example, the estimation fory=Lnx is represented. The variables are gathered in Table 1.Table 1
Observations for variables of the example.
No.X˜≡v,x,z independent variableY˜≡g,y,h response variable11.00,1.10,1.400.00,0.10,0.3421.10,1.30,1.700.10,0.26,0.5331.20,1.60,2.000.18,0.47,0.6941.30,1.90,2.300.26,0.64,0.8351.40,2.20,2.600.34,0.79,0.9661.50,2.50,2.900.41,0.92,1.06As shown in Table1, all observations for two variables are considered triangular fuzzy members. Therefore, an FFLR model will be appropriate for showing the relationship between them:(16)y˜=a˜0⊕a˜1⊗x˜,with unknown fuzzy parameters a˜0 and a˜1. Hence, X˜T⊗X˜⊗a˜=X˜T⊗y˜, and we have(17)6.00,6.00,6.007.50,10.60,12.907.50,10.60,12.909.55,20.16,29.31a˜0a˜1=1.28,3.17,4.411.74,6.44,10.24.Leta˜0=c0,a0,b0anda˜1=c1,a1,b1 be nonnegative triangular fuzzy numbers. Thus, the above FMME is rewritten as(18)6.00,6.00,6.00⊗c0,a0,b0⊕7.50,10.60,12.90⊗c1,a1,b1=1.28,3.17,4.41,7.50,10.60,12.90⊗c0,a0,b0⊕9.55,20.16,29.31⊗c1,a1,b1=1.74,6.44,10.24.By applying fuzzy addition and multiplication operations introduced in Section2,(19)6.00c0+7.50c1,6.00a0+10.60a1,6.00b0+12.90b1=1.28,3.17,4.41,7.50c0+9.55c1,10.60a0+20.16a1,12.90b0+29.31b1=1.74,6.44,10.24.According to the suggested method, instead of FMME, we have the following crisp system:(20)6.00c0+7.50c1=1.28,6.00a0+10.60a1=3.17,6.00b0+12.90b1=4.41,7.50c0+9.55c1=1.74,10.60a0+20.90a1=6.44,12.90b0+29.31b1=10.24,Based on LP, we have the following.Minimizez1+z2+,⋯,+z6,(21)subject to6.00c0+7.50c1+z1=1.28,6.00a0+10.60a1+z2=3.17,6.00b0+12.90b1+z3=4.41,7.50c0+9.55c1+z4=1.74,10.60a0+20.90a1+z5=6.44,12.90b0+29.31b1+z6=10.24,a0−c0≥0,a1−c1≥0,b0−a0≥0,b1−a1≥0,ci≥0,ai≥0,bi≥0,zj≥0,i=0,1,j=1,2,…,6.The optimal solution to the above problem is(22)c0=0.00,a0=0.00,c1=0.00,b0=0.17,a1=0.30,b1=0.34.Also, fuzzy solution problem (15) is the regression equation:(23)a˜0=0.00,0.00,0.00,a˜1=0.17,0.30,0.34,y˜=0.00,0.00,0.00⊕0.17,0.30,0.34⊗x˜.Figure2 represents the response of this example in a range that includes all uncertainties.Figure 2
Chart of response variables vs. independent variable for example 1.The figure consists of three components including lower, middle, and upper linear triple functionsY (lower), Y (middle), and Y (upper), respectively. The figure shows that Y (lower) has a considerable deviation from the main function, while Y (middle) and Y (upper) are suitable approximations for the Ln function, especially after x=1.25.


## 5.2. Example 2

This section includes a case study from the solid mechanics field to show an application of this method. All computations have been performed by RStudio and LINGO software. Figure3 represents the schematic form of the considered cantilever beam. The main output of this system is the value of the moment. Based on statics [34], it can be written as follows:(24)M⇀=r⇀×F⇀,where M⇀ is the moment vector, r⇀ is the distance vector from the origin in which moment acts around to the effect point of force, and finally, F⇀ is the force vector that causes the moment. As it is obvious in Figure 3, two types of forces can be considered including F that is an external force and W which is the weight of the beam that is applied to the center of beam mass. The moment caused by weight can be considered as a constant value, but the moment caused by an external force can vary since its value and effect point can be different. Also, if the distance vector is considered perpendicular to the force vector, the equation can be shown in the scalar form. Therefore, equation (24) can be rewritten as follows:Figure 3
Schematic form of the considered cantilever beam.(25)M=MF+Mw.By substituting the moment components in equation (25), it will be(26)M=RF+rw⟶M=RF+Mo.If we considerF as X˜ and M as Y˜, R and Mo will be a˜1 and a˜0, respectively. Table 2 shows the variables of an example in this area. In this example, eight components are considered as the independent variables, and consequently, eight response variables will be achieved.Table 2
Observations for variables of the example.
No.X˜≡v,x,z independent variableY˜≡g,y,h response variable1(1.5, 2, 2.5)(3.5, 4, 4.5)2(3, 3.5, 4)(5, 5.5, 6)3(4.5, 5.5, 6.5)(6.5, 7.5, 8.5)4(6.5, 7, 7.5)(6, 6.5, 7)5(8, 8.5, 9)(8, 8.5, 9)6(9.5, 10.5, 11.5)(7, 8, 9)7(10.5, 11, 11.5)(10, 10.5, 11)8(12, 12.5, 13)(9, 9.5, 10)As shown in Table2, all observations for two variables are considered triangular fuzzy members. Therefore, an FFLR model will be appropriate for showing the relationship between them:(27)y˜=a˜0⊕a˜1⊗x˜,with unknown fuzzy parameters a˜0 and a˜1. Hence, X˜T⊗X˜⊗a˜=X˜T⊗y˜, and we have(28)8,8,855.5,60.5,65.555.5,60.5,65.5482.25,555.25,635.25a˜0a˜1=55,60,65432,504.5,584.Leta˜0=c0,a0,b0anda˜1=c1,a1,b1 be nonnegative triangular fuzzy numbers. Thus, the above FMME is rewritten as(29)8,8,8⊗c0,a0,b0⊕55.5,60.5,65.5⊗c1,a1,b1=55,60,6555.5,60.5,65.5⊗c0,a0,b0⊕482.25,555.25,635.25⊗c1,a1,b1=432,504.5,584.By applying fuzzy addition and multiplication operations introduced in Section2,(30)8c0+55.5c1,8a0+60.5a1,8b0+65.5b1=55,60,6555.5c0+482.25c1,60.5a0+555.25a1,65.5b0+635.25b1=432,504.5,584.According to the suggested method, instead of FMME, we have the following crisp system:(31)8c0+55.5c1=55,8a0+60.5a1=60,8b0+65.5b1=65,55.5c0+482.25c1=432,60.5a0+555.25a1=504.5,65.5b0+635.25b1=584,and by using LP.Minimizez1+z2+,⋯,+z6,(32)subject  to8c0+55.5c1+z=551,8a0+60.5a1+z2=60,8b0+65.5b1+z3=655,55.5c0+482.25c1+z4=432,60.5a0+555.25a1+z5=504.5,65.5b0+635.25b1+z6=584,a0−c0≥0,a1−c1≥0,ci≥0,ai≥0,bi≥0,zj≥0,i=0,1,j=1,2,…,6.The optimal solution to the above problem is(33)c0=0.876,a0=0.876,b0=1.030,c1=0.795,a1=0.813,b1=0.813.Also, fuzzy solution problem (15) is the regression equation:(34)a˜0=0.876,0.876,1.030,a˜1=0.795,0.813,0.813,y˜=0.876,0.876,1.030⊕0.795,0.813,0.813⊗x˜.These results are shown in Figure4. This figure represents the response of the case study in a range that includes all uncertainties.Figure 4
Chart of response variables vs. independent variable for the cantilever beam.As a sample, ifX˜1=7.75,8.25,9.25 and X˜2=14.5,15,15.5, based on regression equation achievement, we have Y˜1=7.037,7.583,8.550 and Y˜2=12.404,13.071,13.632, respectively. Figure 4 consists of three components including lower, middle, and upper linear triple functions Y (lower), Y (middle), and Y (upper), respectively.


## 6. Conclusion

In this paper, based on the least-square technique in fuzzy regression models, the existing method for estimating coefficients was reviewed in case of having full fuzzy input and output variables. Due to the computation complexity of FFLR, most of the methods have always been criticized. Therefore, to decrease a considerable amount of difficulties, linear programming has been introduced in this study. Besides its simplicity, it can be used when the equation matrix is nonnegative as well. To help understand the FFME application, two numerical examples are provided.



---

*Source: 2684816-2021-10-11.xml*