Abstract

The uncertain nonlinear systems can be modeled with fuzzy equations by incorporating the fuzzy set theory. In this paper, the fuzzy equations are applied as the models for the uncertain nonlinear systems. The nonlinear modeling process is to find the coefficients of the fuzzy equations. We use the neural networks to approximate the coefficients of the fuzzy equations. The approximation theory for crisp models is extended into the fuzzy equation model. The upper bounds of the modeling errors are estimated. Numerical experiments along with comparisons demonstrate the excellent behavior of the proposed method.

1. Introduction

Since the uncertainties in the nonlinear systems can be transformed into fuzzy set theory [1], the fuzzy systems are good models for uncertainty systems. The fuzzy models are based on fuzzy rules. These fuzzy rules give information of the uncertain nonlinear systems. Any nonlinear system can be approximated by several piecewise linear systems (Takagi-Sugeno fuzzy model [2]) or known nonlinear systems (Mamdani fuzzy model [3]). The uncertain nonlinear systems can be modeled by the fuzzy models with simple linear or nonlinear models.

The nonlinear systems can be also modeled with difference equations and algebraic systems. Interpolation methodology has been broadly utilized for function approximation as well as system identification [4, 5]. In [6], the fuzzy polynomial interpolation is applied for system modeling. The coefficients of the polynomials are the fuzzy numbers [7], such that the uncertainties are interpolated with the fuzzy set theory. The theory problem associated with polynomial interpolation is discussed in [8]. It concludes that the interpolation of the function includes time complexity at data points.

In [9], two-dimensional polynomial interpolation is demonstrated. The constraint associated with multivariable interpolation has been investigated in [6], where the Newton-form interpolation is employed. In [10], the multivariate Vandermonde matrix is utilized. Smooth function approximation has been broadly implemented currently [11, 12]. It yields a model by utilizing Lagrange interpolating polynomials at the points of product grids [7, 13]. However if it involves uncertainties in the interpolation points, the above suggested techniques will not work appropriately.

The fuzzy equation can be regarded as a generalized form of the fuzzy polynomial. Compared with the normal fuzzy systems, the fuzzy equations are more easy to be applied, because the uncertainties are direct fuzzy parameters of the fuzzy equations. However, these parameters are not easy to be obtained. There are several approaches to construct the fuzzy equations. Reference [14] utilized the parametric form of fuzzy numbers and restored the original fuzzy equations using crisp linear systems. In [15], the extension principle is implemented and it suggests that the coefficients can be either real or complex fuzzy numbers. Whatsoever, the validation of the solution is not assured. Reference [16] prescribed the homeotypic analysis methodology. Reference [17] inducted Newton’s technique. In [18], the solution of fuzzy equations is extracted using the fixed point methodology. One of the well-known methods is termed as -level [19]. By using the method of superimposition of sets, fuzzy numbers can be solved. Recently, fuzzy fractional differential and integral equations have been extensively studied in [2022]. However, the above methods are very complex.

The numerical solution associated with fuzzy equation can be fetched using the iterative technique [23], interpolation technique [24], and Runge-Kutta technique [25]. It can also be implemented to fuzzy differential equations [26]. These methods are also difficult to be applied. Both neural networks as well as fuzzy logic are considered to be the universal estimators which can estimate any nonlinear function to any notified precision [27]. Current outcomes demonstrate that the fusion methodology of these two different techniques appears to be highly efficient for nonlinear systems identification [28]. Neural networks can also be implemented for resolving the fuzzy equation. In [29], the simple fuzzy quadratic equation is resolved by the neural network method. References [30, 31] elaborated the outcomes of [29] into fuzzy polynomial equation. In [32, 33], the solution of dual fuzzy equations is obtained by neural networks. A matrix pattern associated with the neural learning has been quoted in [34]. However, these techniques are not general; they cannot give the fuzzy coefficients directly with neural networks [35, 36].

We use the neural network method to approximate the coefficients of the fuzzy equations as in our previous paper [37]. In this paper, the standard backpropagation method is modified, such that the fuzzy numbers in the fuzzy ideations can be trained. The approximation theory of the crisp models is extended into the fuzzy equations. The upper bounds of the modeling errors with fuzzy equations are estimated. Finally, we use some real examples to show the effectiveness of our approximation method.

2. Nonlinear System Modeling with Fuzzy Equations

A general discrete-time nonlinear system can be described asHere we consider as the input vector, is regarded as an internal state vector, and is the output vector. and are noted as generalized nonlinear smooth functions . Define and Supposing is nonsingular at the instance , , this will create a path towards the following model:where is a nonlinear difference equation exhibiting the plant dynamics and and are computable scalar input and output, respectively. The nonlinear system which is represented by (2) is implied as a NARMA model. The input of the system with incorporated nonlinearity is considered to be

Taking into consideration the nonlinear systems as mentioned in (2), it can be simplified as the following linear-in-parameter model:where is considered to be the linear parameter and is nonlinear function. The variables related to this function are quantifying input and output.

Many nonlinear systems can be expressed by linear-in-parameter models such as Lagrangian mechanical systems. The parameters of these models are uncertain and the uncertainties satisfy the fuzzy set theory [1]. In this way, the inconvenience problems in nonlinear modeling such as complexity and uncertainty are solved by the fuzzy logic theory and linear-in-parameter structure. The modeling process with the fuzzy equation is to find the fuzzy coefficients of the linear-in-parameter model such that the fuzzy equation can represent the uncertain nonlinear system.

We assume that the model of the nonlinear system (4) has uncertainties in the parameter . The following definitions will be used in this paper.

Definition 1 (fuzzy number). A fuzzy number is a function ; in such a way,    is normal, (there prevail in such a way that );    is convex, , , ;    is upper semicontinuous on ; that is, , , , , is a neighborhood; the set is compact.

In order to demonstrate the fuzzy numbers, the membership functions are utilized. The most widely discussed membership functions are noted to be the triangular functionand trapezoidal function

On a par with crisp variable, the fuzzy variable possesses three essential operations: , , and . They signify these operations: sum, subtract, and multiply.

The fuzzy variable which contains the dimension of is dependent on the membership function, (5) has three variables, and (6) has four variables. In order to define consistency operations, we first apply -level operation to the fuzzy number.

Definition 2 (-level). The -level associated with a fuzzy number is stated aswhere , .

Therefore Since , is bounded mentioned as . The -level of between and is explained as

Definition 3 (fuzzy operations). If , , the fuzzy operations are as follows:
Sum: Subtract:Multiply:Scalar multiplication: .

Definition 4 (absolute value). Absolute value of a triangular fuzzy number is

Now we utilize the following fuzzy equation to model the uncertain nonlinear system (1):Because is fuzzy number, we apply the fuzzy operation .

Taking into consideration a particular case, has polynomial form,Equation (15) is termed as fuzzy polynomial.

Modeling with fuzzy equation (or fuzzy polynomial) can be regarded as fuzzy interpolation. In this paper, we utilize the fuzzy equation (14) to model the uncertain nonlinear system (1), in such a manner that the output related to the plant can approach the desired output ,

This modeling object can be regarded as a way to detect for the following fuzzy equation:where

3. Fuzzy Parameter Estimation with Neural Networks

We design a neural network to represent the fuzzy equation (14); see Figure 1. The input to the neural network is and the output is the fuzzy number The weights are The main idea is to detect appropriate weight of neural network in such a manner that the output of the neural network converges to the desired output .

In order to simplify the operation of the neural network, we use the triangular fuzzy number (5) in this paper. The input fuzzy number is first applied to -level as in (7):Then we haveThe neural network output iswhere , , , and .

In order to train the weights, we need to define a cost function for the fuzzy numbers. The error of the training iswhere , , and The cost function is defined as is considered to be a scalar function. It is quite obvious that means

The vital positiveness that lies within the least mean square (22) is that it has a self-correcting feature that makes it suitable to function for arbitrarily vast duration without shifting from its constraints. The mentioned gradient algorithm is subjected to cumulative series of errors and is convenient for long runs in absence of an additional error rectification procedure. It is more robust in statistics, identification, and signal processing [38].

Now we use gradient method to train the weight defined in (5), or , . We compute asAccording to the chain rulewhere , andwhere  , andalso we havewhere , and where  , and

The coefficient is updated aswhere and is the training rate . For the requirement of increasing the training process, the adding of the momentum term is mentioned aswhere .

Learning Algorithm(1)Step : choose the training rates , , and the stop criterion The initial fuzzy vector is selected randomly. The initial learning iteration is , and the initial learning error .(2)Repeat the following steps for , until all training data are applied.(a)Forward calculation: calculate the -level of the fuzzy output with the -level of the fuzzy input , and the fuzzy connection weight .(b)Backpropagation: adjust fuzzy parameters , , by using the cost function for the -level of the fuzzy output , and the fuzzy target output .(c)Stop criterion: calculate the cycle error , . . If , let ; a new training cycle is initiated. Go to (a).

4. Upper Bounds of the Modeling Errors

In this section, we extend some well-known approximation theories into fuzzy equation modeling. We first define the modeling error in the sense of fuzzy number.

Definition 5. The distance between two fuzzy numbers, , is defined as the Hausdorff metric :

Lemma 6. If is a compact set, then is uniformly support-bounded; that is, there is a compact set , such that, ,

Lemma 7. Let , and , ; then one has (i) if is continuous, holds; (ii) if is continuous, then .

Proof. We need to only prove (ii) since (i) comes from [39]. At first, we demonstrate for . In fact, since , and is closed by the continuity of , hence . On the other hand, for arbitrarily given , there is a sequence and a , such that , . The continuity of implies . But , so . Hence . Thus .
Consideringwe obtain the fact thatholds. Since it may be easily proved that , therefore,which implies the lemma.

Lemma 8. Let be a compact set and , be continuous on ; moreover Then, for each compact set , we have .

Proof. Because of the facts that is a compact set and , are continuous on , then there are , , such thatSupposing , we haveIn the first case of (39), because ,holds, which contradicts (37). In the second case of (39), since , we obtainwhich also contradicts (37). Therefore, (39) is not true; hence , so ; that is, . The proof is completed.

Theorem 9. Let be a continuous function; then for each compact set (the set of all the bounded fuzzy set) and , there are and , , such thatwhere is a finite number.

Proof. The proof of the theorem can be followed from the results below.

If the function , we can extend by the extension principle to the fuzzy function which is also written as as follows:where is called the extended function. Moreover, stands for the set of bounded closed intervals of . ObviouslyMoreoverSo from now on, we suppose

Theorem 10. Let be a continuous function; then for each compact set , , and arbitrary , there are and , , such thatwhere is a finite number. The lower and the upper limits of the -level set of fuzzy function diminish to , but the center goes to .

Proof. Because is a compact set, hence, by Lemma 6, we consider to be the compact set corresponding to . , by the conclusions in [27], there are and , , such thatholds. Let , ; then By Theorem 11, we imply that (47) holds.

Theorem 11. Supposing is compact, the corresponding compact set of and are the continuous functions which satisfy the condition that, for given ,holds. Then, , .

Proof. Let and . Because are continuous, hence , holds by Lemma 7. Therefore, we obtain the following facts by the conclusions from [40]:Because, for : , we have thatholds. And considering the conditions in the theorem, we obtainTherefore, by (51), (52), and Lemma 8, the followingholds, which proves the theorem.

5. Numerical Examples

In this section, we use four examples to show how to apply fuzzy equations to model uncertain nonlinear systems. Here the learning rate in (30) is ; the momentum constant in (31) is .

Example 1. There are three tanks connected to a pipeline at a constant ; see Figure 2. We want to pump water from one tank to the other two tanks. This system satisfies the following fuzzy relationwhere , , and are the quantity of flow and is the elapsed time. is the high of the pipe, , , , and are the characteristic coefficients of the pump; they satisfy triangular uncertainty (5):We use four types of input to train the neural network, whose weights are the fuzzy parameters of (55). The input data arewhere and satisfy the triangle function (5), is the trapezoidal function (6), and is a crisp number. These inputs are applied to (55); the corresponding outputs data areNow these input/output pairs are used repeatedly to train the neural network. The weights are , , , and The training results are shown in Table 1. We can see that the fuzzy characteristic coefficients converge to their real values.

Example 2. The heat source is on the left of the insulating materials; see Figure 3. The conductivity coefficients of these materials are , , and ; here is the elapsed time. The thermal balance is as follows [41]:The real parameters areWe have the following three types of input:where and satisfy the triangle function (5) and is the trapezoidal function (6). The corresponding outputs data areWe repeat them for neural network training. The approximation errors are shown in Figure 4. After times, the parameters , , and converge.

Example 3. The water in the pipe is divided into three pipes , , and ; see Figure 5. The areas of the pipes are uncertain; they satisfy the triangle function (5). , , and . The water velocities in the pipes are controlled by the valves parameters , , , and [42]. The discharge of the water cross in pipe isWe have the following three types of input:where , , and satisfy the triangle function (5). The corresponding outputs data areWe repeat the above input/output times. The neural approximation results of , , and are shown in Table 2.

Example 4. The deformation of a solid cylindrical rod made of special foam depends on stiffness , three different forces on it , , and , the positions of the forces , , and , and the diameter of the rod [36]; see Figure 6. The area of the rod is The external forces are the function of , , , and [43]. According to the tension relations [43],where , , andWe have the following three types of input:where , , and satisfy the triangle function (5). The corresponding outputs data areWe use a neural network to approximate , , and . Convergence is reached after times training. The results are shown in Table 3. To compare our method with the other numerical methods, we use the fuzzy cubic spline method [44] to approximate the solutions of the fuzzy equation (66). The comparison results are shown in Figure 7. We can see that both the neural networks based algorithm (our method) and the fuzzy cubic spline method can approximate the solutions of the fuzzy equations. The approximation errors of the neural networks based algorithm are much smaller than the fuzzy cubic spline method. Also the convergence speed of the neural networks based algorithm is faster. The fuzzy cubic spline method is not robust at the beginning stage.

6. Conclusion

In this paper, the uncertainties in nonlinear systems are modeled by the fuzzy equations and fuzzy numbers. However, the parameters of the fuzzy numbers and fuzzy equations cannot be obtained directly. We construct neural models with the similar structure as the fuzzy equations. By modifying the backpropagation method, the neural networks are trained. The coefficients of the fuzzy equations are approximated by the neural networks. The upper bounds of the fuzzy modeling errors are proven. We successfully extend the approximation theory of crisp models to fuzzy equation model. The novel modeling method is validated with four benchmark examples.

Competing Interests

The authors declare that they have no competing interests.