Abstract

The study of fuzzy variational problems has received significant attention over the past decade due to its successful applications in numerous fields, such as image segmentation and optimal control theory. The fuzzy Euler-Lagrange equations provide the necessary optimality conditions for solving fuzzy variational problems explicitly and have been studied under several differentiability conditions. In this paper, we provide a systematic review to recap the history of variational principle in the calculus of variations and compare it with the existing techniques in the fuzzy setting. We begin with the preliminary concepts and definitions of fuzzy theory and scrutinize the Euler-Lagrange’s strategy via systematically searched studies concerning fuzzy variational problems to highlight the importance of improving the existing methods. Finally, we set up the main open problems regarding the limitations of the current approaches, shedding light on future directions.

1. Introduction

Over the past decade, there has been significant development in the study of fuzzy variational problems due to their practicality in implementation. One of its crucial roles is in the field of image processing, which in turn is used in medical imaging, geographic imaging, and forensic science [13]. The fuzzy variational problems are often solved numerically via Euler-Lagrange equations, Jacobi iterative techniques, or energy difference-based algorithms. One of such examples is the utilization of fuzzy Euler-Lagrange equations by Roul et al. [4] to provide solution for cost minimization while optimizing the production rate. In the study, the considered problem is the minimization of a fuzzy cost functional, where the integrand is in terms of production and stock functions which are fuzzy mappings. The model allows the decision-makers to decide how to optimize the production rate given different stock levels while keeping the cost at its optimal value. On the other hand, multicriteria decision-making (MCDM) has been well studied under fuzzy set-theoretic models [510] to analyze and rank alternatives from the finest to the poorest criteria concerning decision-maker preferences [11]. Since most practical optimization problems can be modelled as fuzzy energy functionals, it is therefore advantageous to study fuzzy variational problems for decision-making. However, the existing studies in fuzzy variational problems are mostly aimed at obtaining the explicit solution by deriving the Euler-Lagrange equations in the fuzzy setting [1214], which are essentially based on the intuition that every extremal problem will have a solution. This is in contrast to the standard variational principle in the calculus of variations that focuses on the existence of minimizers or more generally the existence of solutions.

Known as the classical method for solving variational problems in the calculus of variations, the history of the Euler-Lagrange approach can be traced back to one of the ancient problems in mathematics, namely, the isoperimetric problem or Dido’s problem, that is to determine the curve of a given length which encloses a maximal area. Though a circle appeared to be an obvious solution to the isoperimetric problem, there was no rigorous mathematical proof until 1744 when Euler [15] revisited the isoperimetric problem. The problem was eventually solved by means of the standard methods in calculus. However, Euler noted the complicated process of his geometrical approach, which may require a more straightforward strategy in acquiring the variational conditions. In 1755, Euler received a letter from a nineteen-year-old Lagrange who had found a more general approach to the problem, which is purely analytical. Dropping the geometrical approach, Euler and Lagrange developed a systematic way to solve such variational problems which includes Dido’s problem, Galileo’s brachistochrone problem that was formulated in 1638, Fermat’s optical problem in 1662, and Newton’s problem on floating bodies in 1685. This leads to what we know today as the Euler-Lagrange equation and hence the birth of the calculus of variations.

In view of modern calculus of variations, the Euler-Lagrange equation is usually regarded as the classical method, and one of its disadvantages is that it was established based on the intuition that every extremal problem will have a solution. In particular, while being a solution to the Euler-Lagrange equation is necessary for the minimization problem, it is far from being a sufficient condition. This approach seems to be inadequate since there is no guarantee that a minimizer will exist. Moreover, the solutions are assumed to be regular (usually or ). These limitations led to another class of methods, namely the direct method in the calculus of variations, which is a powerful abstract method for proving the existence of minimizers. This approach grew out of many people’s work; in particular, its early development is often attributed to the works of Riemann, Hilbert, Weierstrass, and Tonelli. The essence of the direct method is to show the existence of a minimizer of a given integral functional and subsequently to prove its regularity. Interested readers on the history of calculus of variations may refer to the books of Giaquinta and Hildebrant [16], Goldstine [17], and Monna [18].

Instead of finding the minimizers of variational problems directly, the direct method focuses on finding these minimizers implicitly by considering minimizing sequences and the concept of lower semicontinuity of the given integral functional. In 1920, Tonelli [19] proved that the functional is lower semicontinuous if and only if its integrand is twice continuously differentiable and convex in the last variable for the one-dimensional case. The work of Tonelli was further extended to the -dimensional case by McShane [20]. However, though the regularity condition of the integrand assumed by McShane [20] and Tonelli [19] was sufficient, it was too strong to be true. Therefore, Serrin [21] tried to relax the requirement of differentiability under the assumptions of continuity and nonnegativity. Later, instead of differentiability and continuity, Marcellini and Sbordone [22] assumed the integrand to be a Carathéodory function (measurable in the first variable and jointly continuous in other variables) and gave the necessary and sufficient conditions for lower semicontinuity of the functional. This led to the general existence theorem with minimal regularity conditions. The variational principle in the calculus of variations has since been applied in studies that vary from the optimal design for thin films [23] to the construction of an ideal column [24] and from quantum field theory [25] to softer spacecraft landings [26]. Apart from that, the development in the calculus of variations has contributed to the advancement of multiple fields of mathematics, including topology, functional analysis, and partial differential equations [27, 28].

One of the fundamental drawbacks of crisp functionals in the calculus of variations is the assumption that all the variables defining a particular energy functional would be accessible accurately. This is a significant restriction since, in most practical cases, not all factors are quantifiable. Recent research studies in image segmentation varying from the detection of the objects in an image to contouring vessel structures highlight similar drawbacks of the conventional crisp functional models, where they are highly dependent on the initial images and are unable to handle objects with ill-defined boundaries [29, 30]. Since not all results can be obtained in the crisp sense of Aristotelian logic due to the complexity of problems, the development of fuzzy functionals seems to be promising to handle such problems. In response to this, our study is aimed at discussing fuzzy functionals that consider both fuzzy logic and calculus of variations.

Many practical problems such as the propagation of wave, heat, fluid flow, and optimization theory [31] are formulated using partial differential equations that can be solved by various means of numerical and analytical methods. One of such methods is the variational method using the Lagrangian density function to solve the corresponding partial differential equations. Though there have been studies focusing on solving the fuzzy partial differential equations to model uncertainties and to study the existence of their solutions (see, for instance, [3236]), there has been little attention on the variational principle for solving fuzzy variational problems. In recent years, some studies have been dedicated to establishing the necessary optimality conditions of fuzzy variational problems for finding the extremals [12, 13, 37]. However, the question remains to be considered is the existence of those extremals without resorting to the Euler-Lagrange equations.

Moreover, real-world problems like image segmentation and optimal control theory generally focus on the minimization of energy integrals that implicity correspond to fuzzy variational problems [4, 29, 38]. They are mainly addressed by iteration techniques under the assumption that a minimizer exists. Despite its importance, there is a lack of comprehensive literature that analyzes the study of the existence of minimizers and the sufficient conditions in the fuzzy setting. Thus, it is crucial to establish fuzzy variational principle such that the uncertainties in real-world problems can be addressed more efficiently. Table 1 summarizes the literature and research gaps that motivate this study as a standardized and systematic review to highlight the importance of the study of the existence of minimizer for fuzzy variational problems.

The outline of this paper goes as follows. In Section 2, we present the preliminary concepts and definitions of fuzzy theory that will be used throughout this paper. In Section 3, we discuss the Euler-Lagrange approach, followed by scrutinizing its methodology and presenting systematically searched studies concerning fuzzy variational problems to reveal the necessity for the improvement of such techniques. Finally, in Section 4, we set up the open problems with respect to the limitations of the current methodology to give an overview of possible directions for future research.

2. Standard Definitions in Fuzzy Theory

In this section, we present some standard definitions that follow from Dubo and Prade [39], Klir and Yuan [40], and Zimmermann [41]. These definitions would be used throughout the paper.

Definition 1 (crisp set). A set is called as a crisp set if there is no ambiguity regarding the inclusion or exclusion of its elements.

Example 1. Let be the set of all natural numbers. Any real number will be either an element of the set or not.

Definition 2 (fuzzy set and membership function). Let be the universe of discourse and represents its elements. Then, a fuzzy set in , denoted by , is defined using a set of ordered pairs as follows: where is called the membership function that gives a degree of membership to every element of ranging from to .

A fuzzy set is a mathematical way of representing ambiguous or subjective data obtained through human language. In a fuzzy set, an element can be partially contained in the set whenever its inclusion is ambiguous, and we represent it by the degree of membership that is between and . For instance, consider a collection of tall students of grade from a school. Now, the term “tall” seems subjective. A person with a height of feet might be tall to someone who is feet, but a feet tall person may not consider feet to be tall. The word “tall” is ambiguous and subjective, and such a collection cannot be defined mathematically using the crisp set. Instead, it can be defined via fuzzy sets using linguistic language.

Similarly, let us consider the case of an automatic air-conditioning system along with a temperature sensor working over a fuzzy algorithm to monitor the speed of a motor that maintains an ideal temperature [46]. A fuzzy inference system generally consists of fuzzification for giving fuzzy inputs, knowledge-based if-then rule evaluation followed by fuzzy outputs, and its defuzzification to obtain a single output of the aggregated fuzzy set. Fuzzification refers to the process of storing linguistic variables as fuzzy sets using membership functions. Figure 1 illustrates the linguistic representation of the air temperature using five fuzzy sets for temperature, namely, Very cool, Cool, Ideal, Warm, and Very warm.

With this fuzzy system, the transitions of temperature from one state to another would be easier as with the temperature change. Subsequently, that can be expressed via a corresponding change in the values of the membership functions. The microcontroller then makes decisions using if-then rules based on human expertise to speed up or slow down the motor and adjust the temperature. In case, if it is warm, which is between 24 and 27 degrees, the speed of the motor is very low, and the temp is very high. In response to this, fuzzy inference system will lower the temperature and increase the speed of the motor. Such rule evaluations go inside the microcontroller of the system that results in a fuzzy output. The defuzzification fuzzy output is also done by the microcontroller to give a crisp solution. In this way, a fuzzy system based on knowledge-based if-then rule mimics human behavior that makes it more practical and applicable [47]. Some of the common membership functions used above can be stated mathematically as follows: (1)Triangular Membership Function. The triangular membership function is given by where , , and are the points with membership values 0, 1, and 0, respectively. The graphical illustration for the fuzzy sets defined using is given in Figure 2.(2)Trapezoidal Membership Function. The trapezoidal membership function is given by where , , , and are the points with membership values , , , and , respectively. The graphical illustration for the fuzzy sets defined using is given in Figure 3.(3)Gaussian Membership Function. The Gaussian membership function is given by where denotes the membership function center and denotes its width. The graphical illustration for the fuzzy sets defined using is given in Figure 4.(4)Generalized Bell Membership Function. The generalized bell membership function is given by where variable is usually positive. In the event when is negative, the shape of this turns to an upside-down. Moreover, is a straightforward generalization of the Cauchy distribution. Here, controls the width of the bell shaped curve, is a positive integer responsible for the arch of the curve, and determines the center of the curve. For instance, consider a fuzzy set We know that the numbers in the neighborhood of on the real line can be considered close to . Therefore, one may define its membership function as , considering the center for some and , depending on the choice of the size for the neighborhood of . For instance, if according to a person numbers like and are equivalently close to , then one may choose to represent it via . On the other hand, if according to another person, 0.4 and 0.45 are also close enough to , then, one may increase the arch to , keeping the width as appropriately. In this way, all such fuzzy sets can be defined using by varying the values of and . The graphical illustration for the fuzzy set defined using for several values of can be represented by Figure 5.

Though there are several types of membership functions, the choice highly depends on the distribution of the given data as discussed in the description of generalized bell membership function above. Therefore, the choice of membership function may not be unique [48] and is in fact problem-specific. In terms of practicality, there are some ways in the field of neural networks like genetic algorithms to determine its choice [49].

Remark 3. While fuzzy sets are written as ordered pairs according to Definition 2, the representation for the fuzzy sets can be in various forms depending on the different choices of membership function provided their range lies within and . For special cases such as the triangular membership function (2) and trapezoidal membership function (3), the corresponding fuzzy set is often denoted by and , respectively, as standard notations. Here, and in the trapezoidal membership function are the values of where the membership function peaks at 1, and the membership value is at points and . We shall henceforth denote these fuzzy sets as such.

Although fuzzy sets were introduced in 1965 by Zadeh [50], it was Mizumoto and Tanaka [51] who proved in 1979 that fuzzy sets do not possess addition and multiplication inverses. Subsequently, a collection of fuzzy sets with such membership functions cannot form a field under the ordinary operations of addition and multiplication. This lack of inverses makes it impossible to solve fuzzy equations. In 1980, Yager [52] presented a methodology to obtain approximate solutions based on the highest degree of truth to which the given solution satisfies the given equation. Similarly, numerous studies were conducted to study the properties of membership functions under the crucial notions of convexity and upper semicontinuity so that basic arithmetic operations can be carried out [39, 5254]. To proceed, we first define fuzzy number and its -level set as follows:

Definition 4 (fuzzy number). Let be the universe of discourse. We call a fuzzy set a fuzzy number if its membership function satisfies the following conditions: (i) is normal, i.e., there exists an such that ;(ii) is fuzzy convex, i.e., for all and ;(iii) is upper semicontinuous;(iv)The closure of the support of is compact.

In the above definition, upper semicontinuity refers to the concept of upper semicontinous function in a crisp sense stated as follows:

Definition 5 (upper semicontinuity). Let be a real valued function. It is said to be upper semicontinuous at if for each there exists such that whenever , where denotes delta neighborhood of . The function is said to be upper semicontinuous on if it is upper semicontinuous at every point of .

We now present the terminology of -level sets that aid in decomposing a fuzzy set into its crisp level sets. In particular, -level sets explicitly express the relation between fuzzy sets and crisp sets [55] that permits us to extend arithmetic operations defined on crisp sets to the case of fuzzy sets, which leads to the reduction in computation cost [54].

Definition 6 (-level set). Let be a fuzzy number and be its membership function. Then, the -level set of is defined as for every .

It can be seen from above that the -level set is a crisp set. We shall henceforth drop tilde from and denote it by . The -level set [] = [,] is a closed interval in for each , where and denote the left hand and right hand endpoints of [], respectively. Properties of and (defined for each ) are as follows: (i) is a bounded increasing function;(ii) is a bounded decreasing function;(iii) for every .

In order to understand the -level set, we now consider the following example.

Example 2. Given a triangular membership function with , and and using the convention in Remark 3, we write the fuzzy set as Note that satisfies all the properties of a fuzzy number stated in Definition 4, and hence, is a fuzzy number. We can write the -level set of as Here, and are bounded increasing and decreasing functions of , respectively. The graphical illustration for the -level set of is given in Figure 6.

Since the -level set represents a closed interval in for each , the operation on -level sets follows that of interval arithmetic. Pioneered by Moore and Yang [56], interval analysis is one of the major tools in developing algorithms for machine computing dealing with rigorous error bounds. One can refer to [57, 58] for operations on interval arithmetic and their properties.

In the following section, we trace out the path of the fuzzy Euler-Lagrange equations, which would turn out to be a fuzzified version of the Euler-Lagrange equation with fuzzy terminologies.

3. Euler-Lagrange Equation for Variational Problems

In this section, we first recall the classical approach of variational principle by means of the Euler-Lagrange equation. We scrutinize the strategy and working process of the Euler-Lagrange theorem and compare it with the existing studies concerning fuzzy variational principle to reveal the necessity for its improvement in addressing the current limitations. To begin, we consider the variational problem in the following form: where is an open bounded set, , is a continuously differentiable function, , and is a twice continuously differentiable function. In particular, we seek to minimize the above integral among all functions from the admissible space . In the one-dimensional case and under ideal conditions, it is known that we can find the critical points of a function by letting its derivative to be equal to zero. Analogously, this can be done for solutions of smooth variational problems with vanishing functional derivative. To observe it, consider the above (11) in one dimensional case, that is when :

Assume that is a curve that minimizes the corresponding in (12). Let be a continuously differentiable function equal to zero on the boundary of , i.e., , such that takes the same values of on the boundary for every . Set where has a minimum when , and hence it must satisfy . In other words, the required extremal curve must satisfy

Taking expressions (13) and (14) into (12) and differentiating under the integral sign lead to

Since is assumed to be a function, we integrate by parts to arrive at where we have used the convention . By the fundamental lemma of calculus of variations, we obtain the Euler-Lagrange equation:

By means of the above Euler-Lagrange equation (17) alone, the existence and uniqueness of the solution are not guaranteed. However, if the integrand is convex, then the existence of minimizer is ensured by Theorem 7 below:

Theorem 7 ([59]). Suppose is a twice continuously differentiable function, and where . (i)If (18) has a minimizer , then necessarily Equation (17) holds(ii)Conversely, if satisfies Equation (17) and if is convex for each , then is a minimizer of (18)(iii)Furthermore, if the function is strictly convex for every , then the minimizer of (18), if it exists, is unique

Theorem 7 can also be generalized to higher dimensions. If , then the Euler-Lagrange equation turns into a system of ordinary differential equations and if where is an open bounded set in , then the Euler-Lagrange equation is a single partial differential equation. If , then we have a system of partial differential equations.

However, one of the disadvantages of the classical method above is that it was established based on the intuition that every extremal problem will have a solution. Moreover, it is not necessary that the minimizer belongs to or even , and if is not convex, then the solution of Equation (17) is not essentially an absolute minimizer of in Equation (18). It may be a local minimizer or a local maximizer (stationary point) of . Solving the Euler-Lagrange equation in the vectorial case may even lead to complex calculations.

Its limitation on the existence of minimizer is addressed via the direct method that begins by considering a minimizing sequence such that as . Under the assumptions of the convexity of in the last variable and the coercivity condition of for some , , , and , we prove that

In particular, the functional in (11) is sequentially weakly lower semicontinuous and hence attains its minimizer.

On the other hand, the functions around us may not always be crisp, hence the need for fuzzy functions to deal with ambiguities and uncertainties. While working well with crisp functions, the variational problems in the calculus of variations cannot be applied to solve ill-posed variational problems with uncertain information. This is because the parameters involved in some functions like time and score are generally uncertain, so their representation using real or crisp numbers leads to an impractical presentation of the real world [60]. To give a practical model, fuzzy numbers can be used as they are adequate for modelling the ambiguity present in the parameters, hence the need for an analytical program to deal with uncertainties that come with fuzzy logic.

In 1987, fuzzy differential equations were first considered by Kaleva [61] to model uncertainties. He gave some existence and uniqueness results for fuzzy differential equations using -level sets. Simultaneously, Seikkala [62] solved a first-order fuzzy initial value problem using Zadeh’s principle along with a generalized Lipschitz condition. Using fuzzy numbers, the approach turns out to be effective in solving problems where the initial values are not known precisely, or the nature of errors is not probabilistic due to subjective choices. Similarly, in 1997, Hüllermeier [63] realized that the available knowledge about functional relationships in a dynamical system can be partially available. Therefore, to establish a mathematical modelling for such systems, he utilized vague linguistic information by human experts to design a fuzzy inference system (FIS). The FIS is defined using if-then rules that approximate partially known functions to the desired degree of accuracy. While studying the corresponding initial value problem with fuzzy numbers using the FIS, Hüllermeier observed that the approach gave better results regarding semantical basis and in terms of predictions when compared with other approaches.

Additionally, fuzzy numbers are also studied in reliability science, where subjective judgment and vague data are utilized to ensure that the system performs more effectively. It was observed that the fuzzy logic-based approach overcomes the limitations of traditional reliability approaches such as underestimating the criticality number when failure mode has various effects and neglecting the relative importance of the factors that give the risk priority number more flexibly using linguistic judgments [64, 65].

Fuzzy theory also plays a crucial role in the field of image segmentation. Image segmentation considers the partitions of an image into segments known as pixels to study its characteristics. It is generally used in medical imaging, geographic imaging, and forensic science [13]. In 1989, Mumford and Shah [38] formulated the image segmentation problem as an energy functional minimization problem to obtain piecewise-smooth contour which segments the given gray scale image using the approximate image . Let and . The proposed energy functional is as follows: where denotes the length of the contour , and and are nonnegative fixed parameters. Later, Chan and Vese [66] considered its simplified model by considering the image as a piecewise constant function as follows: where are fixed parameters, and are constants that define the image intensity of , and and denote the regions outside and inside the contour , respectively. Following a similar strategy, Krinidis and Chatzis [29] introduced a fuzzy energy-based active contour to deal with the image segmentation of the objects whose boundaries are blur or discontinuous as follows: where is the degree of membership of to the inside of the contour and is a weighting exponent on each fuzzy membership.

Employing fuzzy theory in the active contouring for image segmentation allows the pixels to belong to multiclusters based on degrees of the membership of each point to the inside or outside the contour. This aids in extracting more information from the original image to enhance the image segmentation abilities [30, 67, 68]. In the above image segmentation problem, the energy functional can be seen as a fuzzy variational problem. Although such problems are often solved numerically by Jacobi iterative techniques or energy difference-based algorithms, the existence of their minimizers is yet to be established. That motivates us to study and advance the theory of fuzzy variational problems.

3.1. Fuzzy Variational Problems

The Euler-Lagrange method discussed at the beginning of this section is relevant only for well-defined crisp functions. It cannot be applied to solve ill-posed problems with uncertain information. The uncertainties that come with fuzzy logic need to be dealt with, requiring the development of an analytical program by considering fuzzy functions instead of crisp functions. Nevertheless, the differences in concepts and definitions do not allow one to substitute the fuzzy function in the Euler-Lagrange equation in place of the crisp function. For instance, suppose and are two fuzzy numbers, and we want to know that which number out of these two is larger. It is impossible to answer it using the concept of the crisp number as we are dealing with fuzzy numbers that do not stand for an individual value but rather a set of possible values. In view of this, we need the concept of partial ordering using -level sets to compare fuzzy numbers, which one can find in Definition 8 below. Similarly, to study fuzzy variational problems, we first define continuity, differentiability, and integrability of a fuzzy function in terms of -level sets [33, 69].

Definition 8 (partial order on ). Let be the set of all fuzzy numbers, and let . Then, if and for each . Similarly, if and for each Moreover, if and hold for each .

Example 3. Let and be two fuzzy numbers and their -level set representations are given in Figure 7. Since and for each , therefore we conclude that

We would use Definition 6 on -level sets in order to further define the continuity, differentiability, and integrability of a fuzzy-valued function. Therefore, we first recall the definition and the notation for fuzzy-valued function as follows:

Definition 9. Let . The function is said to be a fuzzy-valued function (or fuzzy function) if for every , represents a fuzzy number.
Using Definition 6 of -level set, can be represented as where are bounded increasing and decreasing functions of , respectively. Here, for each , represents an -level set.

Definition 10. A fuzzy function : is said to be (i)continuous at if and are continuous functions of for all .(ii)differentiable at if the derivatives of both (,) and (,) with respect to exist for each fixed , and the interval defines the -level set of a fuzzy number for all and . It is denoted by (iii)integrable with respect to if both and are Riemann integrable functions of , and the interval defines the -level set of a fuzzy number for all . It is denoted by

We now consider the following example to observe the above definitions.

Example 4. Let be a fuzzy function defined by where is a fuzzy number such that = for every . For each , we have the -level set . Here, and are continuous. Moreover, and are differentiable; therefore, is differentiable and

Now, we let be a fuzzy function of belonging to the set of fuzzy functions whose first derivatives are continuous with respect to . The fuzzy variational problem (FVP) can be expressed as:

Here, the integrand assigns a fuzzy number corresponding to every point where and are fuzzy functions of . Further, we assume that have continuous first and second partial derivatives with respect to all of its variables. The function is said to be an admissible curve that satisfies the end-point conditions as well as twice continuously differentiable with respect to . Using all the above concepts and under the assumption that is taken in the left hand terms (similarly for the right hand terms) such as and , we have the following theorem by Farhadinia [13].

Theorem 11 (fuzzy Euler-Lagrange equations). Let = be a twice continuously differentiable fuzzy curve satisfying the given endpoints. If corresponds to a relative minimum of the fuzzy functional in (30), then for all , the following equations necessarily hold:

A direct application of Theorem 11 above can be seen in a fuzzy optimal control model by Roul et al. [4], designed to optimize the production rate so that the total cost is minimized. In contrast to the conventional optimal control models, where all the data is taken in the form of crisp numbers in deterministic imperfect dynamic production models, Roul et al. considered an optimal problem of the minimization of a fuzzy cost functional. The fuzzy variational problem is derived from a dynamic production model with imprecise production due to the variation in raw material quality, workers’ efficiencies, machine failure, and imprecise stock level, resulting in the integrand being in terms of fuzzy mappings and the outputs being represented by the corresponding membership functions. When the membership values are set at 1, the results of the crisp model can also be obtained from the fuzzy model.

In terms of practicality, some of the involved parameters can only be estimated up to an extent due to the presence of uncertainties in real-world optimal control models. Uncertainty and inadequate knowledge are inherent characteristics of modelling mathematical systems. To deal with those uncertainties, real-world problems are often modelled with fuzzy numbers and fuzzy mappings, where the corresponding fuzzy numbers are chosen using prior information, expert knowledge, or genetic algorithm [49]. Such models play an important role in stock trading and map matching algorithm for transport. Chang and Liu [70], for example, considered Takagi–Sugeno–Kang fuzzy-based system using a linear combination of the significant technical index to predict the stock price. While testing on the Taiwan Stock Exchange (TSE) and MediaTek, the proposed method is shown to be an effective trading system for stock price prediction compared to other traditional approaches with an accuracy of % to % approximately. Similar outcomes were observed in the study of fuzzy map matching algorithms by Quddus et al. [71]. Map matching algorithms are generally based on a probabilistic or topological approach defined using a global positioning system (GPS). While it provides good results, they are not always accurate in high road density regions or in confusing traffic circles. To address this concern, a comprehensive field test was carried out on different road environments using the fuzzy logic-based algorithm. It has been shown that it is more efficient to deal with qualitative terms and linguistic vagueness in comparison to the well-recognized existing map-matching algorithms. The outcome of the study shows that the fuzzy logic-based algorithm determined % of the links of the vehicle trajectory correctly.

In mathematical modelling, the concept of differentiability also plays a crucial role in its formulations. It aids in describing the dynamic aspects of the considered system and expresses the functional relationships. In this regard, the concept of fuzzy derivative was first discussed by Pur and Ralescu [72] based on the concept of Hukuhara difference [73]. Since then, Hukuhara differentiability has been widely adopted in studies for proving the existence and uniqueness theorems concerning fuzzy differential equations [61, 62]. A similar notion of the fuzzy derivative has been studied in metric space (see, for instance, [72, 7476]). However, there is a disadvantage of the approach based on the Hukuhara difference, namely it may lead to nondecreasing solutions to the corresponding fuzzy differential equation that can result in increased uncertainties [63]. It was first addressed by Hüllermeier [63] where he used the family of multivalued differential inclusions that considers the fuzzy functions on both sides of the fuzzy differential equations. Although the approach addressed the nondecreasing solution problem, it is worth noting that such derivatives might not provide the result as a fuzzy number. Therefore, to overcome this shortcoming, Buckley and Feuring [77] defined the fuzzy derivative using -level sets and evaluated it using classical differentiation of the endpoints of the -level sets. However, this approach suffers the disadvantage of differential inclusions. In particular, the definition of the derivative is used in the classical sense on -level sets without being defined in the fuzzy sense.

A large number of alternative approaches have been developed over the last few decades to define fuzzy derivatives [7880]. One of the most recent approaches is the notion of generalized Hukuhara (gH) differentiability [78]. Theorem 11 has also been studied under the gH-differentiability by Fard et al. [12]. Subsequently, it was further extended to derive the fuzzy fractional Euler-Lagrange equations for fuzzy fractional variational problems using the Caputo-type fuzzy fractional and Atangana-Baleanu fractional derivatives by Fard and Salehi [42] and Zhang et al. [14], respectively. All such studies employed different differentiability concepts to obtain the necessary optimal conditions for fuzzy variational or fractional variational problems, namely the derivation of the Euler-Lagrange in the fuzzy setting. This motivated Heidari et al. [37] to discuss an approach for the sufficiency conditions regarding the existence of a local minimizer. However, the discussed approach is based on the differential inclusions, studied in the classical sense, by converting the fuzzy variational problem into a crisp variational problem and solving it using the Euler-Lagrange equation. No notions like convexity or strict convexity have been studied in a fuzzy sense without resorting to the Euler-Lagrange equations. Moreover, the uniqueness result has yet been studied. A question that remains to be considered is the existence of those extremals without resorting to the Euler-Lagrange equations since these studies are based on the intuition that every extremal problem will have a solution. Despite that, there has been a lack of extensive literature that analyzes the study of the existence and the uniqueness of minimizers in the fuzzy setting, or more generally, the fuzzy variational principle.

4. Conclusions and Future Directions

In the study of fuzzy variational principle, we observed that the necessary optimality conditions had been obtained by following the fundamental pathway established by Euler and Lagrange for problems in the calculus of variations. Considering twice differentiability, the variational principle via the Euler-Lagrange equations provides the necessary optimality conditions under various differentiability notions such as Buckley-Fuering and generalized Hukuhara (gH) derivatives. This was further extended and studied for fuzzy fractional variational principle under Caputo-type and gH Atangana-Baleanu derivatives. Most real-world problems such as image segmentation in medical and geographic imaging and optimal control theory generally focus on solving the corresponding energy integral, which corresponds to a fuzzy variational problem. This reflects the importance of studying and developing the field of fuzzy variational principle for better decision-making. Although decision-making is a well-established branch under fuzzy set-theoretic modeling, there is room for development in the fuzzy functional setting which can be applied to complex fuzzy mathematical problems when dealing with functionals. Some of the crucial questions are as follows: (1)What would the sufficient and necessary conditions be to ensure the existence of minimizers of the fuzzy variational problems?(2)How would the analogy of the direct method work for fuzzy variational problems?(3)Are fuzzy Euler-Lagrange equations (31) and (32) applicable to higher dimensions, or do they instead lead to complicated calculations?

Given the questions above, it is crucial to work on the sufficient conditions to ensure the existence and uniqueness of the minimizer. The classical approach of the fuzzy variational problem to find solutions via the Euler-Lagrange equations under twice differentiability seems restrictive. We can consider more relaxed regularity conditions to enlarge the space of admissible functions.

This systematic study also shows that the current approach for seeking a minimizer among functions satisfying the fuzzy Euler-Lagrange equations may lead to partial conclusions and complicated calculations in higher dimensions since the existence of a minimizer is assumed and not established beforehand. This limitation motivates us to study the variational principle to ensure the existence of minimizers for variational problems in the fuzzy setting. A similar analogy of the Euler-Lagrange theorem can be established in the fuzzy setting to obtain sufficient conditions for the existence of minimizer.

In terms of applications to real problems, it is necessary to establish the theory for functionals in fuzzy calculus of variations to find an optimized curve for optimal decision-making. For this purpose, one may utilize the calculus of variations’ approach to work with higher dimensions without resulting in complex calculations.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Acknowledgments

This work was supported by the Malaysian Ministry of Higher Education under Fundamental Research Grant Scheme (FRGS/1/2020/STG06/UPM/02/3 01-01-20-2296FR) and by Universiti Putra Malaysia under Postgraduate Initiative Grant (GP-IPS/2021/9695100). The authors would like to thank the editor and the reviewers for their helpful comments and suggestions which have led to the improvement of the earlier version of this paper.