Next Article in Journal
Accessing Information Asymmetry in Peer-to-Peer Lending by Default Prediction from Investors’ Perspective
Previous Article in Journal
Preparation of a Zirconia-Based Ceramic Membrane and Its Application for Drinking Water Treatment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Goal Programming Models with Linear and Exponential Fuzzy Preference Relations

1
College of Science & Theoretical studies, Saudi Electronic University, Riyadh 11673, Saudi Arabia
2
Department of Statistics and Operations Research, Aligarh Muslim University, Aligarh 202002, India
3
Department of Mathematics & Statistics, Integral University, Lucknow 226026, India
4
Department of Engineering Mathematics and Numerics of Optimization, Brandenburg University of Technology Cottbus-Senftenberg, Platz der Deutschen Einheit 1, 03046 Cottbus, Germany
*
Author to whom correspondence should be addressed.
Current Address: Department of Mathematics, Madanapalle Institute of Technology & Science, Chittoor 517325, India.
Symmetry 2020, 12(6), 934; https://doi.org/10.3390/sym12060934
Submission received: 30 April 2020 / Revised: 25 May 2020 / Accepted: 27 May 2020 / Published: 3 June 2020

Abstract

:
Goal programming (GP) is a powerful method to solve multi-objective programming problems. In GP the preferential weights are incorporated in different ways into the achievement function. The problem becomes more complicated if the preferences are imprecise in nature, for example ‘Goal A is slightly or moderately or significantly important than Goal B’. Considering such type of problems, this paper proposes standard goal programming models for multi-objective decision-making, where fuzzy linguistic preference relations are incorporated to model the relative importance of the goals. In the existing literature, only methods with linear preference relations are available. As per our knowledge, nonlinearity was not considered previously in preference relations. We formulated fuzzy preference relations as exponential membership functions. The grades or achievement function is described as an exponential membership function and is used for grading levels of preference toward uncertainty. A nonlinear membership function may lead to a better representation of the achievement level than a linear one. Our proposed models can be a useful tool for different types of real life applications, where exponential nonlinearity in goal preferences exists. Finally, a numerical example is presented and analyzed through multiple cases to validate and compare the proposed models. A distance measure function is also developed and used to compare proposed models. It is found that, for the numerical example, models with exponential membership functions perform better than models with linear membership functions. The proposed models will help decision makers analyze and plan real life problems more realistically.

1. Introduction

Multi-objective mathematical programming is a powerful mathematical procedure and applicable in decision making to a wide range of real life problems [1,2,3,4,5]. A multi-objective mathematical programming model consists of multiple objective functions, constraints, and decision variables and can be stated as follows:
Max ( Min ) : f 1 ( x ) , f 2 ( x ) , , f K ( x ) ,
subject to A x b
x 0 , x S .
where f k ( x ) for all k = 1 , 2 , , K is a function of a decision variable x and S is the set of feasible solutions. Here, Equation (2) represents the constraints where A and b are the m × n and n × 1 matrices, respectively. The ideal solution for a multiple objective programming problem would be to find a solution x, which optimizes each individual objective function of the problem, simultaneously. However, with the conflicting objectives in the model, a feasible solution that optimizes one objective may not optimize any of the other remaining objective functions. Considering this difficulty, several methodologies have been developed to handle the problem of multiple objectives.
Goal Programming (GP) is one of the widely used technique to solve multi-objective optimization problems. It was initially proposed by Charnes and Cooper [6], Charnes et al. [7] and further extended by Caballero et al. [8], Chang [9], Ignizio [10,11], Romero [12,13] and others. GP deals with the over (positive) and under (negative) achievement of the goals in such a manner that each objective is satisfied with the least deviation from its respective goal. In classical GP, the aspiration level of goals are considered as precise and deterministic. However, in real life, when the aspiration levels are not precise in nature, the concept of fuzzy uncertainty can be applied. The fuzzy set theory was first introduced by Zadeh [14]. Bellman and Zadeh [15] further developed a decision theory in a fuzzy environment. Moreover, the uncertainty involved due to vagueness and ambiguity is being dealt within the fuzzy domain.
When fuzzy set theory is introduced in goal programming, it became fuzzy goal programming. Using fuzzy goal programming, the degree of satisfaction of each objective function is maximized which helps the decision maker to decide the best efficient solution. Many researchers proposed and applied fuzzy goal programming models, for example Aköz and Petrovic [16], Anukokila et al. [17], Chen and Tsai [18], Cheng [19], Dalman and Bayram [20], Díaz-Madronero et al. [21], Hocine et al. [22], Jadidi et al. [23], Jamalnia and Soukhakian [24], Jiménez et al. [25], Pramanik and Roy [26], and Jana et al. [27]. It is not always possible that the goals priority is symmetric in nature. The decisions become more complicated, when the problem has imprecise preferences for the priorities of the goals and their aspiration levels. In 2007, Aköz and Petrovic [16] proposed a fuzzy goal programming model with an imprecise linear goal hierarchy. They used a linear membership function for three type of linguistic fuzzy preference relations, i.e., Goal A is “slightly” or “moderately” or “significantly more important” than Goal B. The preference relations are of asymmetric nature in importance of goals. After Aköz and Petrovic [16], the fuzzy preference relation was applied and extended by many authors, e.g., Cheng [19], Petrovic and Aköz [28], Torabi and Moghaddam [29], Khalili-Damghani and Sadi-Nezhad [30], Díaz-Madronero et al. [31], Sheikhalishahi and Torabi [32], Bilbao-Terol et al. [33,34], Arenas-Parra et al. [35] Hasan et al. [36], and Hashmi et al. [37]. Khalili-Damghani and Sadi-Nezhad [30], Khalili-Damghani et al. [38] extended the work of Aköz and Petrovic [16]; they extended the linear membership function for three linguistic preference terms to ten statements. Further, it was applied by Khalili-Damghani et al. [39] in a supply chain network model.
Fuzzy membership functions are used to define the aspiration levels of a decision maker. It describes the nature of fuzzy values, data, preferences, etc. Since the behavior of different fuzzy quantities may vary, the corresponding fuzzy membership functions are of different types, for example, linear, exponential, logistic, hyperbolic, parabolic, S-curve, and tangent. Decision makers may choose the membership function which fits best to their problem. Many authors applied and analyzed their problem models using linear and nonlinear membership functions. A linear membership function is perhaps the simplest and most common one, as it bounds the upper and lower acceptance levels. In nonlinear membership functions, the marginal rate of increase or decrease of membership values as a function of model parameters is not constant. Hence, it reflects the practical situations better than the linear case. Some of the authors who applied nonlinear membership functions are Anoop Kumar Dhingra and Moskowitz [40], Peidro and Vasant [41], Zangiabadi and Maleki [42], Gupta et al. [43], Ehsani et al. [44], Dhodiya and Tailor [45], and Singh and Yadav [46].
In this paper, Standard goal programming models with fuzzy goal hierarchies are proposed. The problem considered is a multi-objective optimization model having linguistic preferences among the objectives, for example “significantly more important”, “moderately more important”, and “slightly more important”, and can be expressed as follows:
FuzzyHierMax : f 1 ( x ) , f 2 ( x ) , , f K ( x ) ,
subject to A x b
x 0 , x S .
Here, “FuzzyHierMax” denotes the fuzzy hierarchical maximization for achieving the satisfactory solution considering the preferences among the objectives. The proposed GP models are the extensions to existing models of Aköz and Petrovic [16], Arenas-Parra et al. [35]. The preferences among the goals f 1 ( x ) , f 2 ( x ) , , f K ( x ) are assumed as fuzzy, having linear and exponential membership functions. The results for different models due to linear and exponential membership functions are further presented and compared. In real life problems, preferences among goals may be of exponential nature. To deal with nonlinearity, we formulated the exponential membership function for ten types of linguistic preference relations, which is the extension to the linear membership function which was earlier proposed by Aköz and Petrovic [16], Khalili-Damghani and Sadi-Nezhad [30], and Khalili-Damghani et al. [38]. Four goal programming models are formulated along with the solution approach. As per our knowledge, in existing literature, there is no GP model with proposed exponential membership function for relative preferences between the goals. In the proposed GP model, we take account of maximization of the minimum and the sum of achievements of goals using approaches from Arenas-Parra et al. [35], Torabi and Hassini [47]. Hence, this work address the research gap, and helps the decision maker to plan more realistically by introducing exponential type preference relations. To analyze the performance of proposed models, we developed a distance measure function. It helps to determine the degree of closeness of the result to the ideal solution.
In this paper, Section 2 presents the proposed goal programming formulations. In Section 3, the solution approach along with flow diagram is presented. In Section 4, an experimental study is performed along with efficiency measure analysis to validate the proposed models. Finally, Section 5 presents the conclusion.

2. Goal Programming Formulation

This section presents the goal programming formulation of proposed models with fuzzy hierarchy among the goals. Before moving to the GP, we first go through the formulation of normalized function and membership functions of preference relations for linear as well as exponential case. To avoid biasedness towards goals with high and low aspiration levels, Arenas-Parra et al. [35] used a normalized function corresponding to each objective. Let us assume that f k ̲ and f k , k = 1 , , K , are the anti-ideal and aspiration values of objective f k , respectively. An anti-ideal value is the worst value of an objective function, whereas an aspiration value is the desired value of an objective function by the decision maker. Now, the normalized function ( n k ) for objective function f k can be modeled as:
n k = 1 ( d k + d k + ) ( f k f k ̲ ) .
In Equation (7), d k and d k + represent the negative and positive deviation of the objectives f k from f k . The deviations must satisfy the condition, d k · d k + = 0 , k 1 , , K . The normalized value satisfies the following property
0 n k 1 .
If the kth objective is more important than the lth objective, then the kth normalized deviation is greater than or equal to that of the lth objective. Hence, we have
n l n k .
However, the normalization constraint may lead to an infeasible solution, as it is a strict relation. To avoid this infeasibility, Arenas-Parra et al. [35] applied the approach of Aköz and Petrovic [16], which expresses preferences among goals in terms of fuzzy relations.

2.1. Representation of Fuzzy Preference Relations as Membership Functions

In 2007, Aköz and Petrovic [16] proposed linear membership functions for three type of fuzzy imprecise preferences namely “slightly more important than”, “moderately more important than”, and “significantly more important than”, to describe the imprecise relative importance relations among the goals. Later in 2013, Khalili-Damghani and Sadi-Nezhad [30] and Khalili-Damghani et al. [38] extended these linear membership functions to ten different linear membership functions. They introduced ten type of linguistic imprecise preference statements alongwith their formulations. The hierarchy structure of ten different linguistic preference relations for linear as well as exponential type is presented in Figure 1. A brief presentation of the figures and linear formulations of membership functions is given in Figure 2 and Equations (10a)–(10j), respectively.
In the present work, the formulation of a type of nonlinear membership function, i.e., exponential membership function, for fuzzy goal hierarchies is presented. More formally, we develop a procedure using a wide range of linguistic terms, different fuzzy priority relations for the decision makers preferences on the normalized values of the goals, and their associated membership functions. The details of the linguistic terms and their associated fuzzy nonlinear membership functions ( μ R ˜ q , q = 1 , , 10 ) are shown in Table 1. The table has four columns with type of linguistic term, notation for the fuzzy relation, membership function, and finally the transform function to formulate the preference membership functions.
μ R ˜ 1 ( k , l ) = 0 , if 1 n k n l < 0 ; 1 , if n k n l = 0 ; 0 , if 0 < n k n l 1 .
μ R ˜ 2 ( k , l ) = 0 , if 1 n k n l 0.5 ; 2 ( n k n l + 0.5 ) , if 0.5 n k n l 0 ; 2 ( n k n l + 0.5 ) , if 0 n k n l 0.5 ; 0 , if 0.5 n k n l 1 .
μ R ˜ 3 ( k , l ) = 2 ( n k n l + 1 ) , if 1 n k n l 0.5 ; 1 , if 0.5 n k n l 1 .
μ R ˜ 4 ( k , l ) = ( n k n l + 1 ) , if 1 n k n l 0 ; 1 , if 0 n k n l 1 .
μ R ˜ 5 ( k , l ) = 2 3 ( n k n l + 1 ) , if 1 n k n l 0.5 ; 1 , if 0.5 n k n l 1 .
μ R ˜ 6 ( k , l ) = ( n k n l + 1 2 ) , if 1 n k n l 1 .
μ R ˜ 7 ( k , l ) = 1 , if 1 n k n l 0.5 ; 2 3 ( n k n l + 0.5 ) , if 0.5 n k n l 1 .
μ R ˜ 8 ( k , l ) = 0 , if 1 n k n l 0 ; ( n k n l ) , if 0 n k n l 1 .
μ R ˜ 9 ( k , l ) = 0 , if 1 n k n l 0.5 ; 2 ( n k n l 0.5 ) , if 0.5 n k n l 1 .
μ R ˜ 10 ( k , l ) = 0 , if 1 n k n l < 1 ; 1 , if n k n l = 1 .
The need for introducing the grade or achievement function as a nonlinear membership function and used for grading levels of preference towards uncertainty is that it may lead to a better representation of the achievement level than a linear one. It is also good enough to represent practical situations because sometimes the marginal rate of increase or decrease of membership values is not constant always, as in the case of linear membership functions. Furthermore, the nature of the behavior of the nonlinear membership function allows itself to change its achievement values according to the parameter values, which also allow the decision maker to execute his/her strategy in a very convenient manner with flexibility. Keeping all these facts in mind, we use the exponential membership function for the achievement of the fuzzy imprecise preferences amongst the goal relations. The proposed exponential fuzzy relation membership functions are shown in Figure 3. The exponential membership function of the fuzzy relations between the goals can be formulated as given in Equations (11a)–(11j).
μ R ˜ 1 ( k , l ) = 0 , if 1 n k n l < 0 ; 1 , if n k n l = 0 ; 0 , if 0 < n k n l 1 .
μ R ˜ 2 ( k , l ) = 0 , if 1 n k n l 0.5 ; 1 e s · 2 ( n k n l + 0.5 ) e s 1 e s , if 0.5 n k n l 0 ; 1 e s · ( 2 ( n k n l + 0.5 ) ) e s 1 e s , if 0 n k n l 0.5 ; 0 , if 0.5 n k n l 1 .
μ R ˜ 3 ( k , l ) = 1 e s · 2 ( n k n l + 1 ) e s 1 e s , if 1 n k n l 0.5 ; 1 , if 0.5 n k n l 1 .
μ R ˜ 4 ( k , l ) = 1 e s · ( n k n l + 1 ) e s 1 e s , if 1 n k n l 0 ; 1 , if 0 n k n l 1 .
μ R ˜ 5 ( k , l ) = 1 e s · ( 2 3 ( n k n l + 1 ) ) e s 1 e s , if 1 n k n l 0.5 ; 1 , if 0.5 n k n l 1 .
μ R ˜ 6 ( k , l ) = 1 e s · ( n k n l + 1 2 ) e s 1 e s , if 1 n k n l 1 ; .
μ R ˜ 7 ( k , l ) = 1 , if 1 n k n l 0.5 ; 1 e s · ( 2 3 ( n k n l + 0.5 ) ) e s 1 e s , if 0.5 n k n l 1 .
μ R ˜ 8 ( k , l ) = 0 , if 1 n k n l 0 ; 1 e s · ( n k n l ) e s 1 e s , if 0 n k n l 1 .
μ R ˜ 9 ( k , l ) = 0 , if 1 n k n l 0.5 ; 1 e s · 2 ( n k n l 0.5 ) e s 1 e s , if 0.5 n k n l 1 .
μ R ˜ 10 ( k , l ) = 0 , if 1 n k n l < 1 ; 1 , if n k n l = 1 .
In an exponential membership function, s is the measurement of the degree of fuzziness. It can be determined by the decision maker. The parameter s is the preference of the decision maker in a fuzzy environment and can be set flexibly in different conditions. Practically, s can be selected according to the level of impreciseness of the problem. A high imprecision in decision making is reflected by a high value of s. Therefore, if the decision environment is associated with a high level of uncertainty, the decision maker should assign a higher value to s. It becomes counter-wise for the cases in which the uncertainty for decision environment is low.

2.2. Proposed Approach

Now, we explain and formulate four standard goal programming models with fuzzy preference relations (SGP-PR). Model 1 is from Arenas-Parra et al. [35] and Models 2–4 are extended goal programming models with imprecise linear and exponential goal hierarchies.
Model 1:
This model is derived from Arenas-Parra et al. [35] and considers a linear membership function for fuzzy preference relations.
Model 2:
This model is an extension of Model 1 with an integration of the proposed formulated exponential membership functions for linguistic preferences among the objectives. This model uses nonlinearities in preferences among the objectives.
Model 3:
It is another extension of Model 1, with a new objective function considering a linear membership function for linguistic preferences between the objectives.
Model 4:
It is an extension of Model 2, incorporating a new objective function and an exponential membership function for linguistic preferences.
The following are the formulations of standard goal programming models with fuzzy hierarchies:
Model 1: This model was initially proposed by Arenas-Parra et al. [35] and later utilized by Bilbao-Terol et al. [33]. Bilbao-Terol et al. [33] applied the extended model in forest planning. Here, Model 1 is the same model proposed by Arenas-Parra et al. [35] with a linear membership function for linguistic preferences among the objectives as explained. The objective functions are defined as:
Objective - 1 : Maximize k = 1 K n k ,
Objective - 2 : Maximize k = 1 K l = 1 K b k l μ R ˜ ( k , l ) .
Here, b k l , for k , l = 1 , 2 , , K with k l ; are binary parameters, taking the value 1 if there exists an importance relation between the goal f k and f l , and 0 otherwise. A convex combination of the two objectives by applying the respective weights α and ( 1 α ) to obtain a compromise solution between the two objectives leads to the following formulation:
Maximize Z = α k = 1 K n k + ( 1 α ) k = 1 K l = 1 K b k l μ R ˜ ( k , l ) ,
subject to n k = 1 ( d k + d k + ) ( f k f k ̲ ) ,
f k + d k d k + = f k ,
( n k n l + 1 ) μ R ˜ 4 ( k , l ) , b k l = 1 and R ˜ ( k , l ) = R ˜ 4 ,
( n k n l + 1 ) 2 μ R ˜ 6 ( k , l ) , b k l = 1 and R ˜ ( k , l ) = R 6 ˜ ,
( n k n l ) μ R ˜ 8 ( k , l ) , b k l = 1 and R ˜ ( k , l ) = R 8 ˜ ,
μ R ˜ ( k , l ) 1 , b k l = 1 ,
A x b ,
x 0 , n k 0 , k = 1 , , K ,
n k 1 , k = 1 , , K ,
d k · d k + = 0 , k = 1 , , K ,
μ R ˜ ( k , l ) 0 , b k l = 1 .
Model 2: This model is the extension of Model 1. It utilizes an exponential membership function for the fuzzy hierarchies between the goals. The formulation is presented below:
Maximize Z = α k = 1 K n k + ( 1 α ) k = 1 K l = 1 K b k l μ R ˜ ( k , l ) ,
subject to n k = 1 ( d k + d k + ) ( f k f k ̲ ) ,
f k + d k d k + = f k ,
1 e s · ( n k n l + 1 ) e s 1 e s μ R ˜ 4 ( k , l ) , b k l = 1 and R ˜ ( k , l ) = R ˜ 4 ,
1 e s · ( n k n l + 1 2 ) e s 1 e s μ R ˜ 6 ( k , l ) , b k l = 1 and R ˜ ( k , l ) = R 6 ˜ ,
1 e s · ( n k n l ) e s 1 e s μ R ˜ 8 ( k , l ) , b k l = 1 and R ˜ ( k , l ) = R 8 ˜ ,
μ R ˜ ( k , l ) 1 , b k l = 1
A x b ,
x 0 , n k 0 , k = 1 , , K ,
n k 1 , k = 1 , , K ,
d k · d k + = 0 , k = 1 , , K ,
μ R ˜ ( k , l ) 0 , b k l = 1 .
Model 3: It is derived from a hybridization of the Torabi and Hassini [47] approach with the goal programming model of Arenas-Parra et al. [35]. Objective 1 in Model 1 only maximizes the sum of satisfaction levels of normalized values; it does not considers a minimum satisfaction of all the objectives. By integrating the Torabi and Hassini [47] approach, we get more flexibility as it considers both the minimum satisfaction and the sum of satisfaction of all the objectives.
Maximize Z = γ 1 ( λ ) + γ 2 k = 1 K n k + γ 3 k = 1 K l = 1 K b k l μ R ˜ ( k , l ) ,
subject to n k = 1 ( d k + d k + ) ( f k f k ̲ ) ,
λ n k ,
f k + d k d k + = f k ,
( n k n l + 1 ) μ R ˜ 4 ( k , l ) , b k l = 1 and R ˜ ( k , l ) = R ˜ 4 ,
( n k n l + 1 ) 2 μ R ˜ 6 ( k , l ) , b k l = 1 and R ˜ ( k , l ) = R 6 ˜ ,
( n k n l ) μ R ˜ 8 ( k , l ) , b k l = 1 and R ˜ ( k , l ) = R 8 ˜ ,
μ R ˜ ( k , l ) 1 , b k l = 1 ,
A x b ,
x 0 , n k 0 , k = 1 , , K ,
n k 1 , k = 1 , , K ,
d k · d k + = 0 , k = 1 , , K ,
μ R ˜ ( k , l ) 0 , b k l = 1 .
Model 4: This model is similar to Model 3, the only difference being in the linguistic preferences, which are represented here by an exponential membership function. Hence, Model 4 is formulated as follows:
Maximize Z = γ 1 ( λ ) + γ 2 k = 1 K n k + γ 3 k = 1 K l = 1 K b k l μ R ˜ ( k , l ) ,
subject to n k = 1 ( d k + d k + ) ( f k f k ̲ ) ,
λ n k ,
f k + d k d k + = f k ,
1 e s · ( n k n l + 1 ) e s 1 e s μ R ˜ 4 ( k , l ) , b k l = 1 and R ˜ ( k , l ) = R ˜ 4 ,
1 e s · ( n k n l + 1 2 ) e s 1 e s μ R ˜ 6 ( k , l ) , b k l = 1 and R ˜ ( k , l ) = R 6 ˜ ,
1 e s · ( n k n l ) e s 1 e s μ R ˜ 8 ( k , l ) , b k l = 1 and R ˜ ( k , l ) = R 8 ˜ ,
μ R ˜ ( k , l ) 1 , b k l = 1 ,
A x b ,
x 0 , n k 0 , k = 1 , , K ,
n k 1 , k = 1 , , K ,
d k · d k + = 0 , k = 1 , , K ,
R ˜ ( k , l ) 0 b k l = 1 .
Here, λ = min{ n k } denotes the minimum normalized function value among the objectives. This formulation has a new achievement function defined as a convex combination of the lower bound ( λ ), the sum of normalized values of objectives n k , and the sum of membership functions for linguistic preferences to ensure an adjustable balanced compromise solution. In Model 3 and 4, γ 1 , γ 2 , γ 3 are weights such that γ 1 + γ 2 + γ 3 = 1 . If we put γ 1 = 0 , then Models 3 and 4 reduce to Models 1 and 2, respectively. Hence, it gives wider applicability to the decision maker by considering extended objective functions.

2.3. Comparison with Other Approaches

A review related to a comparison of the proposed GP models with models from the literature is presented, so as to focus on our present work. The following comparisons are presented based on three aspects.
  • Comparison with respect to the objective function
The objective function defined by Arenas-Parra et al. [35] is a convex combination of two objectives, the sum of the achievement of a normalized function and the sum of the satisfaction degree of fuzzy membership functions representing imprecise goal relations. In our work, Models 1 and 2 have the same objective function as in Arenas-Parra et al. [35], but Models 3, and 4 have an extended objective function. The new objective is the integration of Torabi and Hassini [47] with the Arenas-Parra et al. [35] approach. It considers an addition of the maximization of the minimum achievement of normalized function values of the objectives. Models 3 and 4 can be reduced to Models 1 and 2 by considering weights associated to the minimum normalized functions as zero.
2.
Comparison with respect to Fuzzy preference membership function
Aköz and Petrovic [16] proposed a linear membership function for three linguistic terms. Later, these linear membership functions were extended to ten linguistic terms by Khalili-Damghani and Sadi-Nezhad [30], Khalili-Damghani et al. [38]. To our knowledge, no work in the literature considers nonlinear membership functions for three or ten linguistic preference relations. Here, in our work, we propose the formulation of ten different linguistic terms as exponential membership functions, which is more practical for real life problems.
3.
Comparison with respect to GP models
In our work, we formulate four GP models. Here, Model 1 derived from Arenas-Parra et al. [35] is the base for the other three models. These models are obtained by incorporating different combinations of linear and nonlinear membership functions and objective functions. These proposed models are the extension of the Aköz and Petrovic [16], Arenas-Parra et al. [35] model.

3. Solution Approach

The solution approach to the SGP-PR is consists of multiple steps. Initially, the crisp model of the problem is formulated with all the objectives. Then, to know the best (aspiration) and worst values, either the decision maker is asked or the crisp model is solved individually for each objective function. To proceed further, the solution obtained should be feasible. Take the optimum value as the aspiration value of the objective function. Then, determine all the objective function values for the obtained solution. This will provide us with a payoff table (table of extreme solutions), and the values of aspiration and tolerance level. Formulate all the normalized functions for the objective functions. Construct the transform function for linguistic preference relations by using normalized functions of the objective functions. Define the linguistic preference relations among the different goals and formulate the membership grades of the preference relations. It may be linear or exponential; choose the desired function for fuzzy preference relation.
After, formulating all the functions, now, formulate the problem as SGP-PR as Models 1–4, as detailed in the previous section. Solve the formulation using any optimization software package such as LINGO (Schrage and LINDO Systems [48], Tan et al. [49]), NEOS (Server [50]), AMPL (Fourer et al. [51]), etc. For this paper, the obtained mathematical programming problem was modeled in AMPL language (Fourer et al. [51]) and solved using solvers available on NEOS server online facility provided by Wisconsin Institutes for Discovery at the University of Wisconsin in Madison for solving optimization problems (see Server [50], Gropp, W. Moré [52], Czyzyk et al. [53], Dolan [54]).
As the solution is obtained, check the feasibility and optimality. If it is infeasible, stop and repeat the steps by changing the aspiration and tolerance values, If not, present the solution to the DM. The flow chart of the procedure is presented in Figure 4.

4. Experimental Study

To validate and compare the proposed goal programming models, we applied the GP models on a numerical example. The solution approach used is already detailed in previous section. The main focus is on the validity and the performance of the proposed models. We considered a multi-objective mathematical programming problem, which was already used by Aköz and Petrovic [16], Chen and Tsai [18], Cheng [19], and Arenas-Parra et al. [35]:
Goal 1 : f 1 ( x ) : = 4 x 1 + 2 x 2 + 8 x 3 + x 4 35 ,
Goal 2 : f 2 ( x ) : = 4 x 1 + 7 x 2 + 6 x 3 + 2 x 4 100 ,
Goal 3 : f 3 ( x ) : = x 1 6 x 2 + 5 x 3 + 10 x 4 120 ,
Goal 4 : f 4 ( x ) : = 5 x 1 + 3 x 2 + 2 x 4 70 ,
Goal 5 : f 5 ( x ) : = 4 x 1 + 4 x 2 + 4 x 3 40 ,
subject to 7 x 1 + 5 x 2 + 3 x 3 + 2 x 4 98 ,
7 x 1 + x 2 + 2 x 3 + 6 x 4 117 ,
x 1 + x 2 + 2 x 3 + 6 x 4 130 ,
9 x 1 + x 2 + 6 x 4 105 ,
x i 0 , i = 1 , , 4 .
The type of preference relations between the goals are assumed as:
  • Goal 1 is significantly more important than Goal 2.
  • Goal 2 is significantly more important than Goal 4.
  • Goal 2 is significantly more important than Goal 5.
  • Goal 3 is fully more important than Goal 2.
The solution steps are as follows:
Step:1 Start
Step:2 The optimization model is formulated as single objective problem and result values are obtained for all the objectives. It will provide us the best and worst values for the objectives.
Solving the numerical example using the software packages for each objective function independently and finding the best and worst value of five objective functions, we get f 1 ̲ = 487.66 , f 1 = 261.33 , f 2 = 100 , f 2 ̲ = 0 , f 3 = 120.16 , f 3 ̲ = 117.6 , f 4 = 70 , f 4 ̲ = 0 and f 5 = 40 , f 5 ̲ = 0 . The solution obtained is feasible and optimal.
Step:3 Formulate the normalized function and membership grades for linear and exponential preference relations.
n 1 = 1 d 1 + 226.33 , n 2 = 1 d 2 100 , n 3 = 1 d 3 237.76 , n 4 = 1 d 4 70 , n 5 = 1 d 5 40 ,
Membership functions for linear preference relations:
μ R ˜ 6 ( 1 , 2 ) = ( n 1 n 2 + 1 ) 2 , μ R ˜ 6 ( 2 , 4 ) = ( n 2 n 4 + 1 ) 2 ,
μ R ˜ 6 ( 2 , 5 ) = ( n 2 n 5 + 1 ) 2 , μ R ˜ 8 ( 3 , 2 ) = ( n 3 n 2 )
Membership functions for exponential preference relations:
μ R ˜ 6 ( 1 , 2 ) = 1 e s · ( n 1 n 2 + 1 2 ) e s 1 e s , μ R ˜ 6 ( 2 , 4 ) = 1 e s · ( n 2 n 4 + 1 2 ) e s 1 e s ,
μ R ˜ 6 ( 2 , 5 ) = 1 e s · ( n 2 n 5 + 1 2 ) e s 1 e s , μ R ˜ 8 ( 3 , 2 ) = 1 e s · ( n 3 n 2 ) e s 1 e s
Step:4 Now, formulate the problem as Models 1–4, as detailed in Section 2.2. The formulation of Model 1 becomes:
Maximize Z = α ( n 1 + n 2 + n 3 + n 4 + n 5 ) + ( 1 α ) ( μ R ˜ 6 ( 1 , 2 ) + μ R ˜ 6 ( 2 , 4 ) + μ R ˜ 6 ( 2 , 5 ) + μ R ˜ 8 ( 3 , 2 ) ) ,
subject to 4 x 1 + 2 x 2 + 8 x 3 + x 4 + d 1 d 1 + = 35 ,
4 x 1 + 7 x 2 + 6 x 3 + 2 x 4 + d 2 d 2 + = 100 ,
x 1 6 x 2 + 5 x 3 + 10 x 4 + d 3 d 3 + = 120 ,
5 x 1 + 3 x 2 + 2 x 4 + d 4 d 4 + = 70 ,
4 x 1 + 4 x 2 + 4 x 3 + d 5 d 5 + = 40 ,
n 1 = 1 d 1 + 226.33 , n 2 = 1 d 2 100 , n 3 = 1 d 3 237.76 ,
n 4 = 1 d 4 70 , n 5 = 1 d 5 40 ,
( n 1 n 2 + 1 ) 2 μ R ˜ 6 ( 1 , 2 ) ,
( n 2 n 4 + 1 ) 2 μ R ˜ 6 ( 2 , 4 ) ,
( n 2 n 5 + 1 ) 2 μ R ˜ 6 ( 2 , 5 ) ,
( n 3 n 2 ) μ R ˜ 8 ( 3 , 2 ) ,
7 x 1 + 5 x 2 + 3 x 3 + 2 x 4 98 ,
7 x 1 + x 2 + 2 x 3 + 6 x 4 117 ,
x 1 + x 2 + 2 x 3 + 6 x 4 130 ,
9 x 1 + x 2 + 6 x 4 105 ,
0 μ R ˜ 6 ( 1 , 2 ) , μ R ˜ 6 ( 2 , 4 ) , μ R ˜ 6 ( 2 , 5 ) , μ R ˜ 8 ( 3 , 2 ) 1 ,
x i 0 , n k 0 , i = 1 , , 4 , k = 1 , , 5 ,
0 n k 1 , k = 1 , , 5 .
d k · d k + = 0 ,
μ R ˜ ( k , l ) 0 for all b k l = 1 .
The formulation for Model 2 becomes the same model as Model 1 with replacement of Equations (18i)–(18l) by Equations (19a)–(19d).
1 e s · ( n 1 n 2 + 1 2 ) e s 1 e s μ R ˜ 6 ( 1 , 2 ) ,
1 e s · ( n 2 n 4 + 1 2 ) e s 1 e s μ R ˜ 6 ( 2 , 4 ) ,
1 e s · ( n 2 n 5 + 1 2 ) e s 1 e s μ R ˜ 6 ( 2 , 5 ) ,
1 e s · ( n 3 n 2 ) e s 1 e s μ R ˜ 8 ( 3 , 2 ) ,
The formulation of Model 3 is
Maximize Z = γ 1 ( λ ) + γ 2 ( n 1 + n 2 + n 3 + n 4 + n 5 ) + γ 3 ( μ R ˜ 6 ( 1 , 2 ) + μ R ˜ 6 ( 2 , 4 ) + μ R ˜ 6 ( 2 , 5 ) + μ R ˜ 8 ( 3 , 2 ) ) ,
subject to λ n 1 , n 2 , n 3 , n 4 , n 5 ,
( 18 b ) ( 18 u )
The formulation for Model 4 is as follows:
Maximize Z = γ 1 ( λ ) + γ 2 ( n 1 + n 2 + n 3 + n 4 + n 5 ) + γ 3 ( μ R ˜ 6 ( 1 , 2 ) + μ R ˜ 6 ( 2 , 4 ) + μ R ˜ 6 ( 2 , 5 ) + μ R ˜ 8 ( 3 , 2 ) ) ,
subject to λ n 1 , n 2 , n 3 , n 4 , n 5 ,
( 18 b ) ( 18 h )
( 19 a ) ( 19 d )
( 18 m ) ( 18 u )
Step:5 As the formulation is done, then all the models are modeled in AMPL language (Fourer et al. [51]) and solved by the CONOPT solver (Drud [55]) using the NEOS server online facility provided by Wisconsin Institutes for Discovery at the University of Wisconsin in Madison for solving Optimization problems (see Gropp, W. Moré [52], Czyzyk et al. [53], Dolan [54], and Server [50]). The solution obtained is feasible and optimal.
Step:6 The method is stopped and the solution is presented to the DM.

4.1. Results and Discussion

The results obtained by formulating and solving Models 1–4 on the numerical example are presented in Table 2, Table 3, Table 4 and Table 5, respectively. To analyze the result in more detail, multiple values of weight parameter α and ( γ 1 , γ 2 , γ 3 ) are considered. The decision maker is free to choose these weight parameters α and ( γ 1 , γ 2 , γ 3 ), such that γ 1 + γ 2 + γ 3 = 1 in between 0 snf 1. Ten cases of α and six cases of weight combinations for ( γ 1 , γ 2 , γ 3 ) as c1 = ( 0.1,0.1,0.8), c2 = ( 0.1,0.3,0.6), c3 = ( 0.1,0.8,0.1), c4 = ( 0.3,0.3,0.3), c5 = ( 0.3,0.5,0.2), and c6 = ( 0.6,0.3,0.1) on his/her preferences are considered.
On solving models, many result values for normalized function, decision variables, individual objective function, and membership grade of preference relations between the multiple objectives are obtained. The detailed result values are presented in Table 2, Table 3, Table 4 and Table 5.
In Table 2 and Table 3, the value of n k is increasing and for μ R q ( k , l ) is decreasing from α = 0 to 1. In Table 2, the values for n k are the same for α = 0–0.3 as 2.5828 and α = 0.5–1 as 4.7655. In Table 3, the values of n k for α = 0–0.2 and α = 0.5–1 are same as Model 1. For all the values of μ R q ( k , l ) , Model 2 is higher and thus better than Model 1 for all values of α .
In addition, all decision variables of Table 2, Table 3, Table 4 and Table 5 for Models 3 and 4 present one more decision variable λ , which is the minimum value of the normalized function ( n k ) among all the objective functions f k . Table 4 and Table 5 present the solution values of Models 3 and 4, respectively. On comparing Models 3 and 4, for n k and μ R q ( k , l ) , we find that Model 4 performs better than Model 3 for all cases of ( γ 1 , γ 2 , γ 3 ). The value of n k is the same for c1, c3, c4, c5, and c6 for Models 3 and 4, whereas μ R q ( k , l ) M o d e l 4 > μ R q ( k , l ) M o d e l 3 for all combinations of ( γ 1 , γ 2 , γ 3 ). In the case of c2, n k of Model 4 has a higher value than Model 3.
To get more insight into the analysis of Models 1–4 with linear and exponential preference relations, we further solved the numerical example for five types of preference relations combinations. The considered preference relation types are as follows:
Preference Relation Type-1 (PRT-1):
The preference relations between the goals are:
  • Goal 1 is significantly more important than Goal 2.
  • Goal 2 is significantly more important than Goal 4.
  • Goal 2 is significantly more important than Goal 5.
  • Goal 3 is fully more important than Goal 2.
Preference Relation Type-2 (PRT-2):
The preference relations between the goals are:
  • Goal 1 is significantly more important than Goal 2.
  • Goal 2 is significantly more important than Goal 4.
  • Goal 2 is significantly more important than Goal 5.
  • Goal 3 is significantly more important than Goal 2.
Preference Relation Type-3 (PRT-3):
The preference relations between the goals are:
  • Goal 1 is significantly more important than Goal 2.
  • Goal 2 is partially equal to Goal 4.
  • Goal 2 is significantly more important than Goal 5.
  • Goal 3 is fully more important than Goal 2.
Preference Relation Type-4 (PRT-4):
The preference relations between the goals are:
  • Goal 1 is significantly more important than Goal 2.
  • Goal 2 is completely more important than Goal 4.
  • Goal 2 is moderately more important than Goal 5.
  • Goal 3 is fully more important than Goal 2.
Preference Relation Type-5 (PRT-5):
The preference relations between the goals are:
  • Goal 1 is significantly more important than Goal 2.
  • Goal 2 is significantly more important than Goal 4.
  • Goal 2 is significantly more important than Goal 5.
  • Goal 3 is partially more important than Goal 2.
A total of twenty five different models are formulated. After modeling and solving the problem models, the results are presented in Table 6, Table 7, Table 8, and Table 9. The values of n k and μ R q ( k , l ) are determined for each preference relation combination on Models 1–4. The values of n k is increasing while μ R q ( k , l ) is decreasing from α = 0 to 1 for all PR type in Models 1 and 2. In both models, at α = 0, the value of n k is the same for PRT-1, -2, and -4 as 2.5828, whereas, at α = 1, it is teh same for PRT-1, -3, and -4 as 4.7655 and PRT-2 and -5 as 4.786. The value of μ R q ( k , l ) is 0 for all PR Type. On comparing the preference relation membership grades ( μ R q ( k , l ) ) for Models 1 and 2 for all PR types, Model 2 performs best.
For Model 3, in c1 the value of n k is the same for PRT1, -2, and -4. The highest value is for PRT5 at 3.2. In Model 4, for c1 all the values are greater than the values of Model 3. For Model 3, c2 has the value of 4.4577 for PRT1 and 04 but in Model 4 PRT1 and 04 have lesser values compared to Model 3. For c3, all the values in Models 3 and 4 are the same from PRT1–5. The highest one is for c5 as 4.7616. For c4, all the PRTs have same values for Models 3 and 4, except for PRT3 the value is different. For c4, all the values of Models 3 and 4 are same, but for c5 only PRT4 and -5 values are different.
If we compare μ R q ( k , l ) for Models 3 and 4, most of the values for Model 4 for all weights combinations are greater than the values of Model 3. Although the decision can be taken based on the requirement by the decision maker.

4.2. Efficiency Analysis

The proposed models were applied to solve the numerical problem with different types of preference relations combinations. The problem models are linear and nonlinear programming models. While solving the models, we found that the solution time is in seconds for both linear and nonlinear models. The computational time has been found acceptable.
Here, Models 1 and 3 and Models 2 and 4 consist of linear and exponential type of preference relations, respectively. We have previously analyzed the models based on values of n k and μ R q ( k , l ) which is not a good basis for analyzing the efficiency and performance of proposed models. Thus, to determine the efficiency of different models, we determine the degree of closeness of the solution value to the ideal solution. For this the concept of distance function is utilized, see the works of Pramanik and Roy [56], Zheng et al. [57], and Zhao et al. [58]. In this paper, a new distance function is developed to select the priority solution for the proposed models.
D ( x ) = ( k ( 1 n k ) 2 + q ( 1 μ R q ( k , l ) ) 2 ) 1 / 2
The smaller value of the distance function gives the most satisfying solution as it represents the better solution among all other solutions. The results of numerical example for the distance function is presented in Table 10 and Table 11 and Figure 5 and Figure 6.
Furthermore, the proposed approach in Models 1 and 2 considers a utility function divided into two criteria: (1) sum of individual normalized function; and (2) sum of membership function for the achievement level of preference relations. Based on these two criteria with a different value of α , the performance of different models with five different combinations of preference relations is analyzed. The proposed approach of Models 3 and 4 considers a different and new utility function with addition of maximization of λ , i.e., minimum value among all the normalized functions. Models 1 and 3 consider membership function for linear preference relations, whereas Models 2 and 4 consider membership function for Exponential preference relation.
In Table 10 and Figure 5, the results for distance function values are shown for different values of α and for type of preference relation model. In Figure 5, it can be seen that for all models of PRT-1–5, Model 2 performs best for all values of α . In PRT-2 and -5, the difference of values for Models 1 and 2 is slightly greater than for PRT-1, -3, and -5. For PRT-1 and -4, at α = 0.3 , the difference is more in between Model 1 and 2 for distance function values. It can be seen that, for PRT-5, at α = 0.1–−0.3, the performance of model is best among all other models. It gives the lowest value, 0.838 and 0.6339 for linear and exponential cases, respectively.
In Table 11 and Figure 6, the results for distance function values are shown for different values of ( γ 1 , γ 2 , γ 3 ) and for type of preference relation model. For PRT-1, -2, -4, and -5, Model 4 with exponential preference relation performs best. In the case of PRT-3 for c2, Model 3 is better than Model 4; for the rest of the weight combinations Model 4 is better than Model 3. For PRT-5, the distance function values at c2 is the lowest among all distance function values for different types of preference relations. Overall, the performance ranking is M o d e l 2 > M o d e l 1 and M o d e l 4 > M o d e l 3 . Thus, We can conclude that, for the given experimental study, the performance pattern for the type of membership function is E x p o n e n t i a l > L i n e a r . For both examples, the results show that the proposed models perform best for the exponential case in comparison to the linear one.

5. Conclusions, Limitations, and Future Directions

This work proposes standard goal programming models with fuzzy hierarchies between the goals. This work is based on several existing methods (see, [16,30,35,38]) and provides an extension to the existing methods by introducing exponential fuzzy preference relations and incorporating new objective function. The exponential membership function for fuzzy preference relations of ten types of linguistic terms is formulated. Four models, namely Models 1–4, based on linear, exponential preference relations, and a new objective are proposed. A numerical example is solved and results are analyzed. A comparison is also made using a newly developed distance measure function for the proposed models to analyze the performance of the proposed models.
It is concluded that, for the experimental numerical problem, solutions are better in case of models having exponential membership function for preference relations in comparison to linear. The proposed models perform better than Model 1, which was proposed by Arenas-Parra et al. [35]. Thus, the efficiency of models with respect to the satisfaction level of the decision maker can be ordered as M o d e l 2 > M o d e l 1 and M o d e l 4 > M o d e l 3 . Another flexibility in the proposed models is that the solution can be chosen according to the priority of the minimum normalized value, the sum of normalized values, and the sum of membership functions of the preference relations. That is, the decision maker can adopt the best solution by having a sensitivity analysis of weight parameters. On choosing the weight γ 1 = 0 , Models 3 and 4 are reduced to Models 1 and 2, respectively.
This study also hase some limitations. As, the results are presented for a single numerical example under different scenarios of preference relations. Further investigation on multiple numerical examples and real life problems are needed to verify the superiority of proposed linear and nonlinear models. The results may change as per the consideration of type of linear and nonlinear preference relations. This research can be extended by the implementation in real life problems from industry such as in scheduling, planning, transportation, manufacturing, supply chain, etc., which can be a challenging work. The inclusion of exponential preferences provides wider applicability in real life scenario. The present work can be further extended to advanced goal programming models using fuzzy, stochastic uncertainties, and different nonlinear membership functions such parabolic, hyperbolic, S-curve, etc. for preference relations.

Author Contributions

Conceptualization, Formal analysis, Investigation, Methodology, Software, Validation, Visualization, Writing—original draft by M.G.H.; Software, Visualization, Writing—review & editing by A.Q.; Writing—review & editing by A.F., M.F.K. and S.S.H. Funding acquisition by M.F.K. All authors have read and agreed to the published version of the manuscript.

Funding

The first author is thankful to Saudi Electronic University, Riyadh Saudi Arabia, for providing a financial assistance to carry out this research work.

Acknowledgments

The authors are thankful to the anonymous referees and the journal editor for their fruitful constructive comments and suggestions for the improvement of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hwang, C.L.; Masud, A.S.M. Multiple Objective Decision Making—Methods and Applications: A State-of-the-Art Survey; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2012; Volume 164. [Google Scholar]
  2. Lai, Y.J.; Hwang, C.L. Fuzzy multiple objective decision making. In Fuzzy Multiple Objective Decision Making; Springer: Berlin/Heidelberg, Germany, 1994; pp. 139–262. [Google Scholar]
  3. Sakawa, M.; Nishizaki, I.; Katagiri, H. Fuzzy Stochastic Multiobjective Programming; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2011; Volume 159. [Google Scholar]
  4. Sakawa, M.; Yano, H.; Nishizaki, I.; Nishizaki, I. Linear and Multiobjective Programming with Fuzzy Stochastic Extensions; Springer: Berlin/Heidelberg, Germany, 2016. [Google Scholar]
  5. Tzeng, G.H.; Huang, J.J. Fuzzy Multiple Objective Decision Making; Chapman and Hall/CRC: London, UK, 2016. [Google Scholar]
  6. Charnes, A.; Cooper, W.W. Management models and industrial applications of linear programming. Manag. Sci. 1957, 4, 38–91. [Google Scholar] [CrossRef]
  7. Charnes, A.; Cooper, W.W.; Ferguson, R.O. Optimal estimation of executive compensation by linear programming. Manag. Sci. 1955, 1, 138–151. [Google Scholar] [CrossRef]
  8. Caballero, R.; Gómez, T.; Ruiz, F. Goal programming: Realistic targets for the near future. J. Multi-Criteria Decis. Anal. 2009, 16, 79–110. [Google Scholar] [CrossRef]
  9. Chang, C.T. Efficient structures of achievement functions for goal programming models. Asia-Pac. J. Oper. Res. 2007, 24, 755–764. [Google Scholar] [CrossRef]
  10. Ignizio, J.P. A review of goal programming: A tool for multiobjective analysis. J. Oper. Res. Soc. 1978, 29, 1109–1119. [Google Scholar] [CrossRef]
  11. Ignizio, J.P. Generalized goal programming An overview. Comput. Oper. Res. 1983, 10, 277–289. [Google Scholar] [CrossRef]
  12. Romero, C. Handbook of Critical Issues in Goal Programming; Elsevier Science: Amsterdam, The Netherlands, 1991; pp. 1–12. [Google Scholar] [CrossRef]
  13. Romero, C. A general structure of achievement function for a goal programming model. Eur. J. Oper. Res. 2003, 153, 675–686. [Google Scholar] [CrossRef]
  14. Zadeh, L. Fuzzy sets. Inf. Control 1965, 8, 338–353. [Google Scholar] [CrossRef] [Green Version]
  15. Bellman, R.; Zadeh, L. Decision Making in a Fuzzy Environment. Manag. Sci. 1970, 17, 141–164. [Google Scholar] [CrossRef]
  16. Aköz, O.; Petrovic, D. A fuzzy goal programming method with imprecise goal hierarchy. Eur. J. Oper. Res. 2007, 181, 1427–1433. [Google Scholar] [CrossRef]
  17. Anukokila, P.; Radhakrishnan, B.; Anju, A. Goal programming approach for solving multi-objective fractional transportation problem with fuzzy parameters. RAIRO-Oper. Res. 2019, 53, 157–178. [Google Scholar] [CrossRef] [Green Version]
  18. Chen, L.H.; Tsai, F.C. Fuzzy goal programming with different importance and priorities. Eur. J. Oper. Res. 2001, 133, 548–556. [Google Scholar] [CrossRef]
  19. Cheng, H.W. A satisficing method for fuzzy goal programming problems with different importance and priorities. Qual. Quant. 2013, 47, 485–498. [Google Scholar] [CrossRef]
  20. Dalman, H.; Bayram, M. Interactive goal programming algorithm with Taylor series and interval type 2 fuzzy numbers. Int. J. Mach. Learn. Cybern. 2019, 10, 1563–1579. [Google Scholar] [CrossRef]
  21. Díaz-Madroñero, M.; Pérez-Sánchez, M.; Satorre-Aznar, J.R.; Mula, J.; López-Jiménez, P.A. Analysis of a wastewater treatment plant using fuzzy goal programming as a management tool: A case study. J. Clean. Prod. 2018, 180, 20–33. [Google Scholar] [CrossRef]
  22. Hocine, A.; Kouaissah, N.; Bettahar, S.; Benbouziane, M. Optimizing renewable energy portfolios under uncertainty: A multi-segment fuzzy goal programming approach. Renew. Energy 2018, 129, 540–552. [Google Scholar] [CrossRef]
  23. Jadidi, O.; Zolfaghari, S.; Cavalieri, S. A new normalized goal programming model for multi-objective problems: A case of supplier selection and order allocation. Int. J. Prod. Econ. 2014, 148, 158–165. [Google Scholar] [CrossRef]
  24. Jamalnia, A.; Soukhakian, M.A. A hybrid fuzzy goal programming approach with different goal priorities to aggregate production planning. Comput. Ind. Eng. 2009, 56, 1474–1486. [Google Scholar] [CrossRef]
  25. Jiménez, M.; Bilbao-Terol, A.; Arenas-Parra, M. A model for solving incompatible fuzzy goal programming: An application to portfolio selection. Int. Trans. Oper. Res. 2018, 25, 887–912. [Google Scholar] [CrossRef]
  26. Pramanik, S.; Roy, T.K. Multiobjective Transportation Model with Fuzzy Parameters: Priority based Fuzzy Goal Programming Approach. J. Transp. Syst. Eng. Inf. Technol. 2008, 8, 40–48. [Google Scholar] [CrossRef]
  27. Jana, R.K.; Sharma, D.K.; Chakraborty, B. A hybrid probabilistic fuzzy goal programming approach for agricultural decision-making. Int. J. Prod. Econ. 2016, 173, 134–141. [Google Scholar] [CrossRef]
  28. Petrovic, D.; Aköz, O. A fuzzy goal programming approach to integrated loading and scheduling of a batch processing machine. J. Oper. Res. Soc. 2008, 59, 1211–1219. [Google Scholar] [CrossRef]
  29. Torabi, S.A.; Moghaddam, M. Multi-site integrated production-distribution planning with trans-shipment: A fuzzy goal programming approach. Int. J. Prod. Res. 2012, 50, 1726–1748. [Google Scholar] [CrossRef]
  30. Khalili-Damghani, K.; Sadi-Nezhad, S. A decision support system for fuzzy multi-objective multi-period sustainable project selection. Comput. Ind. Eng. 2013, 64, 1045–1060. [Google Scholar] [CrossRef]
  31. Díaz-Madroñero, M.; Mula, J.; Jiménez, M. Fuzzy goal programming for material requirements planning under uncertainty and integrity conditions. Int. J. Prod. Res. 2014, 52, 6971–6988. [Google Scholar] [CrossRef]
  32. Sheikhalishahi, M.; Torabi, S. Maintenance supplier selection considering life cycle costs and risks: A fuzzy goal programming approach. Int. J. Prod. Res. 2014, 52, 7084–7099. [Google Scholar] [CrossRef]
  33. Bilbao-Terol, A.; Jiménez, M.; Arenas-Parra, M. A group decision making model based on goal programming with fuzzy hierarchy: An application to regional forest planning. Ann. Oper. Res. 2016, 245, 137–162. [Google Scholar] [CrossRef]
  34. Bilbao-Terol, A.; Arenas-Parra, M.; Cañal-Fernández, V.; Jiménez, M. A sequential goal programming model with fuzzy hierarchies to Sustainable and responsible portfolio selection problem. J. Oper. Res. Soc. 2016, 67, 1259–1273. [Google Scholar] [CrossRef]
  35. Arenas-Parra, M.; Bilbao-Terol, A.; Jiménez, M. Standard goal programming with fuzzy hierarchies: A sequential approach. Soft Comput. 2016, 20, 2341–2352. [Google Scholar] [CrossRef]
  36. Hasan, M.G.; Qayyum, Z.; Hasan, S.S. Multi-Objective Annualized Hours Manpower Planning Model: A Modified Fuzzy Goal Programming Approach. Ind. Eng. Manag. Syst. 2019, 18, 52–66. [Google Scholar] [CrossRef]
  37. Hashmi, N.; Jalil, S.A.; Javaid, S. A model for two-stage fixed charge transportation problem with multiple objectives and fuzzy linguistic preferences. Soft Comput. 2019, 1–15. [Google Scholar] [CrossRef]
  38. Khalili-Damghani, K.; Sadi-Nezhad, S.; Tavana, M. Solving multi-period project selection problems with fuzzy goal programming based on TOPSIS and a fuzzy preference relation. Inf. Sci. 2013, 252, 42–61. [Google Scholar] [CrossRef]
  39. Khalili-Damghani, K.; Tavana, M.; Amirkhan, M. A fuzzy bi-objective mixed-integer programming method for solving supply chain network design problems under ambiguous and vague conditions. Int. J. Adv. Manuf. Technol. 2014, 73, 1567–1595. [Google Scholar] [CrossRef]
  40. Anoop, K.D.; Moskowitz, H. Application of fuzzy theories to multiple objective decision making in system design. Eur. J. Oper. Res. 1991, 53, 348–361. [Google Scholar] [CrossRef]
  41. Peidro, D.; Vasant, P. Transportation planning with modified S-curve membership functions using an interactive fuzzy multi-objective approach. Appl. Soft Comput. 2011, 11, 2656–2663. [Google Scholar] [CrossRef] [Green Version]
  42. Zangiabadi, M.; Maleki, H. Fuzzy goal programming technique to solve multiobjective transportation problems with Some non-linear membership functions. Iran. J. Fuzzy Syst. 2013, 10, 61–74. [Google Scholar]
  43. Gupta, P.; Mehlawat, M.; Mittal, G. A fuzzy approach to multicriteria assignment problem using exponential membership functions. Int. J. Mach. Learn. Cybern. 2013, 4. [Google Scholar] [CrossRef]
  44. Ehsani, E.; Kazemi, N.; Olugu, E.U.; Grosse, E.E.H.; Schwindl, K. Applying fuzzy multi-objective linear programming to a project management decision with nonlinear fuzzy membership functions. Neural Comput. Appl. 2016, 28, 2193–2206. [Google Scholar] [CrossRef]
  45. Dhodiya, J.M.; Tailor, A.R. Genetic algorithm based hybrid approach to solve fuzzy multi-objective assignment problem using exponential membership function. SpringerPlus 2016, 5. [Google Scholar] [CrossRef] [Green Version]
  46. Singh, S.; Yadav, S. Intuitionistic fuzzy multi-objective linear programming problem with various membership functions. Ann. Oper. Res. 2017. [Google Scholar] [CrossRef]
  47. Torabi, S.A.; Hassini, E. An interactive possibilistic programming approach for multiple objective supply chain master planning. Fuzzy Sets Syst. 2008, 159, 193–214. [Google Scholar] [CrossRef]
  48. Schrage, L.E.; LINDO Systems, I. Optimization Modeling with LINGO; Duxbury Press: Belmont, CA, USA, 1997. [Google Scholar]
  49. Tan, R.R.; Aviso, K.B.; Promentilla, M.A.B.; Yu, K.D.S.; Santos, J.R. Programming in LINGO. In Input-Output Models for Sustainable Industrial Systems; Springer: Berlin/Heidelberg, Germany, 2019; pp. 29–46. [Google Scholar]
  50. Server, N. State-of-the-Art Solvers for Numerical Optimization. 2016. Available online: https://neos-server.org/neos/index.html (accessed on 30 May 2020).
  51. Fourer, R.; Gay, D.M.; Kernighan, B.W. AMPL-A Modeling Language for Mathematical Programming, 2nd ed.; Cengage Learning; Cengage Learning: Stamford, CT, USA, 2003. [Google Scholar]
  52. Gropp, W.; Moré, J. Optimization environments and the NEOS server. In Approximation Theory and Optimization; Buhman, M.D., Iserles, A., Eds.; Cambridge University Press: Cambridge, UK, 1997; pp. 167–182. [Google Scholar]
  53. Czyzyk, J.; Mesnier, M.P.; More, J.J. NEOS server. IEEE Comput. Sci. Eng. 1998, 5, 68–75. [Google Scholar] [CrossRef]
  54. Dolan, E. The NEOS Server 4.0 Administrative Guide; Technical Report, Memorandum ANL/MCS-TM-250; Mathematics and Computer Science Division, Argonne National Laboratory: Argonne, IL, USA, 2001. [Google Scholar]
  55. Drud, A.S. CONOPT—Alarge-scale GRG code. ORSA J. Comput. 1994, 6, 207–216. [Google Scholar] [CrossRef]
  56. Pramanik, S.; Roy, T.K. Fuzzy goal programming approach to multilevel programming problems. Eur. J. Oper. Res. 2007, 176, 1151–1166. [Google Scholar] [CrossRef]
  57. Zheng, Y.; Liu, J.; Wan, Z. Interactive fuzzy decision making method for solving bilevel programming problem. Appl. Math. Model. 2014, 38, 3136–3141. [Google Scholar] [CrossRef]
  58. Zhao, X.; Zheng, Y.; Wan, Z. Interactive intuitionistic fuzzy methods for multilevel programming problems. Expert Syst. Appl. 2017, 72, 258–268. [Google Scholar] [CrossRef]
Figure 1. Hierarchy structure of linguistic preference relations.
Figure 1. Hierarchy structure of linguistic preference relations.
Symmetry 12 00934 g001
Figure 2. Linear membership functions for linguistic preferences among objectives.
Figure 2. Linear membership functions for linguistic preferences among objectives.
Symmetry 12 00934 g002
Figure 3. Exponential membership functions for linguistic preferences among objectives.
Figure 3. Exponential membership functions for linguistic preferences among objectives.
Symmetry 12 00934 g003
Figure 4. Flow chart of solution approach.
Figure 4. Flow chart of solution approach.
Symmetry 12 00934 g004
Figure 5. Graph of distance function values of Models 1 and 2 for PRT 1–5.
Figure 5. Graph of distance function values of Models 1 and 2 for PRT 1–5.
Symmetry 12 00934 g005aSymmetry 12 00934 g005b
Figure 6. Graph of distance function values of Model 3 and 4 for PRT 1–5.
Figure 6. Graph of distance function values of Model 3 and 4 for PRT 1–5.
Symmetry 12 00934 g006aSymmetry 12 00934 g006b
Table 1. Linguistic relative preferences of objective k over l.
Table 1. Linguistic relative preferences of objective k over l.
Linguistic TermFuzzy RelationMembership ValuesTransform Function
R ˜ q , q = 1 , , 10 μ R ˜ q
Exactly equal R ˜ 1 μ R ˜ 1
Partially equal R ˜ 2 μ R ˜ 2
Partially more important than R ˜ 3 μ R ˜ 3
Slightly more important than R ˜ 4 μ R ˜ 4
Moderately more important than R ˜ 5 μ R ˜ 5 μ k ( X ) μ l ( X )
Significantly more important than R ˜ 6 μ R ˜ 6 k , l 1 , , K
Completely more important than R ˜ 7 μ R ˜ 7
Fully more important than R ˜ 8 μ R ˜ 8
Extremely more important than R ˜ 9 μ R ˜ 9
Incomparable R ˜ 10 μ R ˜ 10
Table 2. Solution results for Model 1 (Linear).
Table 2. Solution results for Model 1 (Linear).
VariablesModel 1 (Linear)
α 00.10.20.30.40.50.60.70.80.91
Z2.70852.6962.68342.670862.84683.16623.4863.80594.12574.44564.7655
f 1 1212121245.913946.381646.381646.381646.381646.381646.3816
f 2 24242424100100.526100.526100.526100.526100.526100.526
f 3 120120120120120120120120120120120
f 4 2424242457.016857.105357.105357.105357.105357.105357.1053
f 5 000039.6639404040404040
x 1 00000000000
x 2 00008.25638.28948.28948.28948.28948.28948.2894
x 3 00001.65961.71051.71051.71051.71051.71051.7105
x 4 1212121216.123916.118416.118416.118416.118416.118416.1184
n 1 11110.95170.94970.94970.94970.94970.94970.9497
n 2 0.240.240.240.241111111
n 3 11111111111
n 4 0.34280.34280.34280.34280.81450.81570.81570.81570.81570.81570.8157
n 5 00000.9915111111
μ 12 0.880.880.880.880.47580.47480.47480.47480.47480.47480.4748
μ 24 0.44850.44850.44850.44850.59270.59210.59210.59210.5921050.59210.5921
μ 32 0.760.760.760.760000000
μ 25 0.620.620.620.620.50420.50.50.50.50.50.5
n k 2.58282.58282.58282.58284.75794.76554.76554.76554.76554.76554.7655
μ R q ( k , l ) 2.70852.70852.70852.70851.57281.56691.56691.56691.56691.56690
Table 3. Solution results for Model 2 (Exponential).
Table 3. Solution results for Model 2 (Exponential).
VariablesModel 2 (Exponential)
α 00.10.20.30.40.50.60.70.80.91
Z3.07073.02192.97312.94743.09693.34643.63023.9144.19784.48164.7655
f 1 121212353546.381646.381646.381646.381646.381646.3816
f 2 24242445.746384.8031100.526100.526100.526100.526100.526100.526
f 3 120120120134.416120120120120120120120
f 4 24242427.166254.179857.105357.105357.105357.105357.105357.1053
f 5 00013.312731.0515404040404040
x 1 00000.2852000000
x 2 0000.69446.97848.28948.28948.28948.28948.28948.2894
x 3 0002.63370.49911.71051.71051.71051.71051.71051.7105
x 4 12121212.541415.90916.118416.118416.118416.118416.118416.1184
n 1 111110.94970.94970.94970.94970.94970.9497
n 2 0.240.240.240.45740.848111111
n 3 11111111111
n 4 0.34280.34280.34280.38800.77390.81570.81570.81570.81570.81570.8157
n 5 0000.33280.7762111111
μ 12 0.92580.92580.92580.85040.69260.5980.59800.59800.59800.59800
μ 24 0.57180.57180.57180.65510.65730.70680.70680.70680.70680.70680
μ 32 0.84210.84210.84210.66240.223000000
μ 25 0.73090.73090.73090.68040.65620.62240.62240.62240.62240.62240
n k 2.58282.58282.58283.17834.39834.76554.76554.76554.76554.76554.7655
μ R q ( k , l ) 3.07073.07073.07072.84842.22931.92731.92731.92731.92731.92730
Table 4. Solution results for Model 3 (Linear).
Table 4. Solution results for Model 3 (Linear).
VariablesModel 3 (Linear)
( γ 1 , γ 2 , γ 3 )c1c2c3c4c5c6
Z2.42512.47854.05062.14442.94082.0899
f 1 123546.381646.381646.381665.1907
f 2 2487.7187100.526100.526100.526121.692
f 3 120120120120120120
f 4 2454.952757.105357.105357.105360.6625
f 5 031.820340404053.5177
x 1 000000
x 2 07.48228.28948.28948.28949.62345
x 3 00.47281.71051.71051.71053.7559
x 4 1216.25316.118416.118416.118415.8961
λ 00.78500.81570.81570.81570.8666
n 1 110.94970.94970.94970.8666
n 2 0.240.87711111
n 3 111111
n 4 0.34280.78500.81570.81570.81570.8666
n 5 00.79551111
μ 12 0.880.56140.47480.47480.47480.4333
μ 24 0.44850.54600.59210.59210.59210.5666
μ 32 0.760.12280000
μ 25 0.620.54080.50.50.50.5
n k 2.58284.45774.76554.76554.76554.7332
μ R q ( k , l ) 2.70851.77111.56691.56691.56691.5
Table 5. Solution results for Model 4 (Exponential).
Table 5. Solution results for Model 4 (Exponential).
VariablesModel 4 (Exponential)
( γ 1 , γ 2 , γ 3 )c1c2c3c4c5c6
Z2.71482.74444.08672.25263.01292.1262
f 1 123546.381646.381646.381665.1907
f 2 2470.4306100.526100.526100.526121.692
f 3 120120120120120120
f 4 2450.369857.105357.105357.105360.6625
f 5 027.261440404053.5177
x 1 01.69160000
x 2 04.49498.28948.28948.28949.6234
x 3 00.62871.71051.71051.71053.7559
x 4 1214.213416.118416.118416.118415.8961
λ 00.68150.81570.81570.81570.8666
n 1 110.94970.94970.94970.8666
n 2 0.240.70431111
n 3 111111
n 4 0.34280.71950.81570.81570.81570.8666
n 5 00.68151111
μ 12 0.92580.75430.59800.59800.59800.5562
μ 24 0.57180.61510.70680.70680.70680.6843
μ 32 0.84210.4049005.38E-170
μ 25 0.73090.63330.62240.62240.62240.6224
n k 2.58284.10544.76554.76554.76554.7332
μ R q ( k , l ) 3.07072.40771.92731.92731.92731.8631
Table 6. Solution results of Model 1 (Linear) for all cases of preference relations.
Table 6. Solution results of Model 1 (Linear) for all cases of preference relations.
PR TypeModel 1 (Linear)
α 00.10.20.30.40.50.60.70.80.91
PRT1 n k 2.58282.58282.58282.58284.75794.76554.76554.76554.76554.76554.7655
μ R q ( k , l ) 2.70852.70852.70852.708571.57281.56691.566961.566961.566961.566960
PRT2 n k 2.58284.45774.45774.45774.75794.76554.76554.76554.7864.7864.786
μ R q ( k , l ) 2.82852.20972.20972.20972.07282.06692.06692.06692.0122.0120
PRT3 n k 2.82.82.82.84.15544.15544.63644.76554.76554.76554.7655
μ R q ( k , l ) 3.13.13.13.12.42222.42221.82491.60641.60641.60640
PRT4 n k 2.58282.58282.58282.58284.75794.76554.76554.76554.76554.76554.7655
μ R q ( k , l ) 2.73142.73142.73142.73141.60511.59761.59761.59761.59761.59760
PRT5 n k 1.49492.854.76164.76164.76164.7864.7864.7864.7864.7864.786
μ R q ( k , l ) 32.9252.57592.57592.57592.55762.55762.55762.55762.55760
Table 7. Solution results of Model 2 (Exponential) for all cases of preference relations.
Table 7. Solution results of Model 2 (Exponential) for all cases of preference relations.
PR TypeModel 2 (Exponential)
α 00.10.20.30.40.50.60.70.80.91
PRT1 n k 2.58282.58282.58283.178374.39834.76554.76554.76554.76554.76554.7655
μ R q ( k , l ) 3.07073.07073.07072.848452.22931.92731.927371.927371.927371.927370
PRT2 n k 2.58282.58283.24.57164.75794.76554.76554.76554.7864.7864.786
μ R q ( k , l ) 3.15433.15433.03722.63782.55542.54982.54982.54982.49872.49870
PRT3 n k 2.75262.82.83.14774.15544.15544.76554.76554.76554.76554.7655
μ R q ( k , l ) 3.37653.37573.37573.23722.7572.7571.96121.96121.96121.96120
PRT4 n k 2.58282.58282.82613.18764.45774.76554.76554.76554.76554.76554.7655
μ R q ( k , l ) 3.02573.02572.97342.84712.19121.94721.94721.94721.94721.94720
PRT5 n k 1.49493.24.76164.76164.76164.7864.7864.7864.7864.7864.786
μ R q ( k , l ) 3.24493.19552.93822.93822.93822.9212.9212.9212.9212.9210
Table 8. Solution results of Model 3 (Linear) for all cases of preference relations.
Table 8. Solution results of Model 3 (Linear) for all cases of preference relations.
VariablesModel 3 (Linear)
PR Type( γ 1 , γ 2 , γ 3 )c1c2c3c4c5c6
PRT1 n k 2.58284.45774.76554.76554.76554.7332
μ R q ( k , l ) 2.70851.771131.56691.56691.56691.5
PRT2 n k 2.58284.75794.7864.7864.78514.7556
μ R q ( k , l ) 2.82852.07282.0122.0121.99651.9592
PRT3 n k 2.84.15544.76554.49874.73324.7332
μ R q ( k , l ) 3.12.42221.60642.00851.66651.6665
PRT4 n k 2.58284.45774.76554.76554.76554.7655
μ R q ( k , l ) 2.73141.80011.59761.59761.59761.5976
PRT5 n k 3.24.76164.7864.78514.78514.7851
μ R q ( k , l ) 2.88122.57592.55762.54852.54852.5485
Table 9. Solution results for Model 4 (Exponential) for all cases of preference relations.
Table 9. Solution results for Model 4 (Exponential) for all cases of preference relations.
VariablesModel 4 (Exponential)
PR Type( γ 1 , γ 2 , γ 3 )c1c2c3c4c5c6
PRT1 n k 2.58284.10544.76554.76554.76554.7332
μ R q ( k , l ) 3.07072.40771.92731.92731.92731.8631
PRT2 n k 2.994.75794.7864.7864.78514.7556
μ R q ( k , l ) 3.07872.55542.49872.49872.48392.4483
PRT3 n k 4.71544.71744.76554.73324.73324.7332
μ R q ( k , l ) 2.01842.01651.96122.00082.00082.0008
PRT4 n k 2.71214.0974.76554.76554.76554.7332
μ R q ( k , l ) 3.00012.40241.94721.94721.94721.8709
PRT5 n k 4.41264.76164.7864.78514.78514.7835
μ R q ( k , l ) 3.01362.93822.9212.91262.91262.91
Table 10. Solution results for distance measure for Models 1 and 2.
Table 10. Solution results for distance measure for Models 1 and 2.
PR TypeDistance Measure for Model 1 and 2
α 00.10.20.30.40.50.60.70.80.91
PRT1Model 11.59051.59051.59051.59051.31261.31471.31471.31471.31471.31472.009
Model 21.51511.51511.51511.21291.02851.19431.19431.19431.19431.19432.009
PRT2Model 11.57691.57691.57690.95110.98640.98920.98920.98921.00791.00792.0056
Model 21.50871.50871.18160.72740.75180.75430.75430.75430.76870.76872.0056
PRT3Model 11.39091.39091.39091.39091.06631.06631.24241.3031.3031.3032.009
Model 21.35911.33451.33451.17890.90840.90841.18651.18651.18651.18652.009
PRT4Model 11.62851.62851.62851.62851.30911.31111.31111.31111.31111.31112.009
Model 21.56571.56571.40321.22581.0691.19491.19491.19491.19491.19492.009
PRT5Model 11.93781.46260.8380.8380.8380.84750.84750.84750.84750.84752.0056
Model 21.88151.1710.63390.63390.63390.64230.64230.64230.64230.64232.0056
Table 11. Solution results for distance measure for Models 3 and 4.
Table 11. Solution results for distance measure for Models 3 and 4.
PR TypeDistance Measure for Model 3 and 4
( γ 1 , γ 2 , γ 3 )c1c2c3c4c5c6
PRT1Model 31.59051.21721.31471.31471.31471.3395
Model 41.51510.9821.19431.19431.21431.2143
PRT2Model 31.57690.98641.00791.00791.01511.0322
Model 41.26820.75180.76870.76870.77540.7912
PRT3Model 31.39091.06631.3031.25451.29531.2953
Model 41.19281.19181.18651.1861.1861.186
PRT4Model 31.62851.22841.31111.31111.31111.3422
Model 41.47631.02581.19491.19491.19491.2217
PRT5Model 31.2580.8380.84750.85190.85190.8519
Model 40.66690.63390.64230.64610.64610.6472

Share and Cite

MDPI and ACS Style

Khan, M.F.; Hasan, M.G.; Quddoos, A.; Fügenschuh, A.; Hasan, S.S. Goal Programming Models with Linear and Exponential Fuzzy Preference Relations. Symmetry 2020, 12, 934. https://doi.org/10.3390/sym12060934

AMA Style

Khan MF, Hasan MG, Quddoos A, Fügenschuh A, Hasan SS. Goal Programming Models with Linear and Exponential Fuzzy Preference Relations. Symmetry. 2020; 12(6):934. https://doi.org/10.3390/sym12060934

Chicago/Turabian Style

Khan, Mohammad Faisal, Md. Gulzarul Hasan, Abdul Quddoos, Armin Fügenschuh, and Syed Suhaib Hasan. 2020. "Goal Programming Models with Linear and Exponential Fuzzy Preference Relations" Symmetry 12, no. 6: 934. https://doi.org/10.3390/sym12060934

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop