Next Article in Journal
Applications of the q-Sălăgean Differential Operator Involving Multivalent Functions
Next Article in Special Issue
An Interactive Consensus Model in Group Decision Making with Heterogeneous Hesitant Preference Relations
Previous Article in Journal
Supersymmetric Polynomials and a Ring of Multisets of a Banach Algebra
Previous Article in Special Issue
Evaluation of Cryptocurrencies for Investment Decisions in the Era of Industry 4.0: A Borda Count-Based Intuitionistic Fuzzy Set Extensions EDAS-MAIRCA-MARCOS Multi-Criteria Methodology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of Integrated Linear Programming Fuzzy-Rough MCDM Model for Production Optimization

1
Higher Colleges of Technology Sharjah, University City, Sharjah P.O. Box 7947, United Arab Emirates
2
Faculty of Economics, PIM University, Despota Stefana Lazarevića bb, 78000 Banja Luka, Bosnia and Herzegovina
3
Environmental Protection and Energy Efficiency Fund of Republic of Srpska, King Alfonso XIII Number 21, 78000 Banja Luka, Bosnia and Herzegovina
4
Ministry of Spatial Planning, Construction and Ecology of the Republic of Srpska, Trg Republike Srpske, 78000 Banja Luka, Bosnia and Herzegovina
5
Sustainable Transportation Research Group, Civil Engineering, School of Engineering, University of Kwazulu Natal, Durban 4041, South Africa
6
Faculty of Transport and Traffic Engineering, University of East Sarajevo, Vojvode Mišića 52, 74000 Doboj, Bosnia and Herzegovina
7
College of Engineering and Technology, American University of the Middle East, Egaila 54200, Kuwait
*
Author to whom correspondence should be addressed.
Axioms 2022, 11(10), 510; https://doi.org/10.3390/axioms11100510
Submission received: 13 September 2022 / Revised: 24 September 2022 / Accepted: 26 September 2022 / Published: 27 September 2022
(This article belongs to the Special Issue Fuzzy Set Theory and Its Applications in Decision Making)

Abstract

:
One of the most common tools for achieving optimization and adequate production process management is linear programming (LP) in various forms. However, there are specific cases of the application of linear programming when production optimization implies several potential solutions instead of one. Exactly such a problem is solved in this paper, which integrates linear programming and a Multi-Criteria Decision-Making (MCDM) model. First, linear programming was applied to optimize production and several potential solutions lying on the line segment AB were obtained. A list of criteria was created and evaluated using the Improved Fuzzy Stepwise Weight Assessment Ratio Analysis (IMF SWARA). To obtain the final solution, a novel Rough compromise ranking of alternatives from distance to ideal solution (R-CRADIS) method was developed and verified through comparative analysis. The results show that the integration of linear programming and a Fuzzy-Rough MCDM model can be an exceptional solution for solving specific optimization problems.
MSC:
03B52; 90C08; 90B50; 91B06; 91B38

1. Introduction

Today, linear programming is one of the most effective mathematical tools for solving economic problems, certainly with the support of computer systems. Along with the development of computer systems, there has been an intensified application of linear programming and other related functions to solve specific economic problems. By applying linear programming, greater efficiency is achieved, and production costs can be significantly reduced, as indicated in [1] in which the authors applied the Mixed-Integer Linear Programming model. This is also confirmed by the authors Chandrawat et al. [2] applying fuzzy linear programming to optimize production costs because, in order to optimize these processes, we must make changes that are allow managers to be well informed [3]. However, taking into account that market conditions change very quickly and that there is uncertainty in every segment of the economic system, integrated models are most frequently applied to optimize various problems, as is the case in this paper, in which we presented a novel Rough CRADIS model in integration with the IMF SWARA method and linear programming. The application of linear programming and multi-criteria models is not rare [4,5,6,7]. Cheng et al. [8] created an integration model of multi-criteria decision analysis (MCDA) and inexact mixed integer linear programming (IMILP) to support the selection of an optimal landfill location and waste-flow-allocation pattern so that the total system cost can be minimized. This model includes qualitative and quantitative indicators, which is an advantage considering that the disadvantage of various multi-objective programming models is that they are basically mathematical and often ignore qualitative and subjective factors [8]. The paper [9] uses an integrated fuzzy multi-criteria model and a multi-objective programming approach for supplier selection and order allocation in a green supply chain. First, the fuzzy analytic hierarchy process and fuzzy TOPSIS were applied in order to analyze the significance of criteria and determine the best green suppliers. Then, multi-objective linear programming was used to consider and formulate various constraints such as quality control, capacity and other objectives. The objective of the mathematical model is simultaneously to maximize the total value of purchasing and to minimize the total cost of purchasing. The subject of research in the paper [10] refers to an integrated approach of multi-attribute utility theory (MAUT) and linear programming (LP) for evaluating and selecting the best suppliers and defining optimal order quantities among the selected ones to maximize total utility.
The following can be identified as special contributions of this paper:
(1) Formation of a novel integrated model consisting of linear programming, IMF SWARA method and Rough CRADIS approach.
(2) Extending the CRADIS approach with Rough Numbers (R-CRADIS) and presenting it for the first time in the literature, which is an enrichment of the entire field that treats multi-criteria problems.
(3) Solving the special case of linear programming obtaining multiple optimal solutions and integration with a MCDM model in order to achieve the desired optimum.
(4) Carrying out a sensitivity analysis based on which decision-makers can make decisions in real time, by looking at the real needs of the company and market requirements at any moment, and by considering the simulated values of criteria and obtaining new optimal solutions.
Motivation and objectives of the paper can be manifested through the following. When applying LP for the optimization and management of production processes in this special case, several potential solutions are obtained instead of one which is usually the case. In order to determine one solution that is optimal under the given circumstances and taking into account various factors, integration with a novel MCDM model was carried out. Moreover, the aim is to develop a model ensuring managers use real-time decision-making.
The rest of the paper consists of the following sections. In Section 2, materials and methods are given. The overall flow of the research is presented graphically and the steps of the developed method and applied methods are explained. The development of a novel Rough CRADIS approach is presented in detail. Section 3 introduces the optimization of production with the aim of maximizing profit. A special case of linear programming and its integration with a fuzzy-rough MCDM model is given. Section 4 is also of great importance because it presents the verification of the proposed model through sensitivity analysis, comparative analysis and calculation of statistical SCC and WS coefficients to determine the correlation of the ranks obtained. Finally, Section 5 provides concluding considerations with guidelines for future actions.

2. Materials and Methods

2.1. Methodology Steps

Figure 1 presents the overall methodology of this paper, which is divided into four different phases (I, II, III and IV). The first phase represents the expressed need for solving special problems by LP optimization. Additionally, in this phase, solving the problem of optimizing the production plan is presented with the aim of determining the maximum profit.
When applying LP for the optimization and management of production processes in this special case, several potential solutions are obtained instead of one which is usually the case. In order to determine one solution that is optimal under the given circumstances and taking into account various factors, the integration with a novel MCDM model was carried out. This represents the second phase of the research. First, four criteria were defined: capacity utilization of the first group of machines M (C1), capacity utilization of the second group of machines N (C2), supplying the market with articles of group A (C3) supplying the market with articles of group B (C4). Alternatives represent all potential solutions for x1 = 2000–2500, and for x2 = 1000–2000. A total of 21 potential solutions were created, and then the data were transformed into an initial decision matrix. The third phase involves the development of the novel Rough CRADIS approach (in which computation was made in Microsoft Excel), which represents a contribution from the methodological aspect. The developed approach is used to evaluate potential solutions, while the IMF SWARA method is used to calculate the weights of criteria. The fourth phase is determining the verification of the proposed model through several steps. Using sensitivity analysis and the formation of 40 scenarios, it is possible for decision-makers to determine the optimal solution in real time, which is also one of the contributions of this research both from the methodological and professional aspect. Then the obtained results were compared with four other Rough MCDM methods: Rough WASPAS [11], Rough SAW [12], Rough ARAS [13] and Rough MARCOS [14]. An additional test of the rank stability of potential solutions was also performed by creating 20 new scenarios in which the worst alternative was removed, and the entire procedure was repeated. Finally, statistical correlation tests were calculated for both sensitivity analysis and comparative analysis. Thus, SCC [15] and WS coefficients were calculated [16].

2.2. Linear Programming

The general task of LP: find a non-negative solution
X = ( x 1 , x 2 , , x n ) , ( x i 0 , i = 1 , 2 , , n )
of the systems of linear inequalities or equations with constraints:
a 11 x 1 + a 12 x 2 + + a 1 n x n { = } b 1 a 21 x 1 + a 22 x 2 + + a 2 n x n { = } b 2 a m 1 x 1 + a m 2 x 2 + + a m n x n { = } b m
so that the linear objective function reaches an extremum, i.e., a minimum or maximum value.
F = F ( x 1 , x 2 , , x n ) = c 1 x 1 + c 2 x 2 + + c n x n
The constraints in the LP problem form a domain of admissible solutions, which contains a set of solutions C ¯ . Each point in C ¯ is a candidate for a solution in the LP problem and is denoted as an admissible solution. A set C ¯ represents a set of admissible solutions. Among all points in the set C ¯ , the point(s) that optimize(s) the objective function of the LP problem is called the optimal solution of the LP problem.

2.3. IMF SWARA Method

The IMF SWARA method has been created by Vrtagić et al. [17] and consists of the following steps [14,18]:
Step 1: Sorting criteria in descending order based on their expected significance.
Step 2: Relatively smaller significance of each criterion (criterion Cj) is determined in relation to the previous one (Cj − 1), and this is repeated for each subsequent criterion. The comparative significance of the average value is denoted by j ¯ . Linguistics and the TFN scale for assessment of criteria is shown in Table 1.
Step 3: Determining the fuzzy coefficient j ¯ from Equation (4):
j ¯ = { 1 ¯ j = 1 j ¯ 1 ¯ j > 1
The comparative significance of the average value— j ¯ can then be found.
Step 4: Determining the calculated weights j ¯ from Equation (5):
j ¯ = { 1 ¯ j = 1 j 1 ¯ j ¯ j > 1
j ¯ is the fuzzy coefficient from the previous step.
Step 5: Calculation of the fuzzy weight coefficients, from Equation (6):
w j ¯ = j ¯ j = 1 n j ¯
where w j ¯ represents the fuzzy relative weight of the criteria j, and n represents the number of criteria.

2.4. Development of a Novel Rough CRADIS Approach

In this section of the paper, the CRADIS approach [19,20,21] was extended with Rough numbers. The developed Rough CRADIS approach is described through the following steps and equations.
Step 1. The first step includes the formation of an initial Rough decision matrix consisting of n criteria and m alternatives as shown in Equation (7),
ξ = [ [ ξ 11 L , ξ 11 U ] [ ξ 21 L , ξ 21 U ] [ ξ m 1 L , ξ m 1 U ] [ ξ 12 L , ξ 12 U ] [ ξ 22 L , ξ 22 U ] [ ξ m 2 L , ξ m 2 U ] [ ξ 1 n L , ξ 1 n U ] [ ξ 2 n L , ξ 2 n U ] [ ξ m n L , ξ m n U ] ]
where ξ i j denotes the performance of alternative m in relation to criterion n.
Step 2. Determination of a Rough normalized decision matrix (8) using Equations (9) and (10):
ς = [ [ ς 11 L , ς 11 U ] [ ς 21 L , ς 21 U ] [ ς m 1 L , ς m 1 U ] [ ς 12 L , ς 12 U ] [ ς 22 L , ς 22 U ] [ ς m 2 L , ς m 2 U ] [ ς 1 n L , ς 1 n U ] [ ς 2 n L , ς 2 n U ] [ ς m n L , ς m n U ] ]
ς i j = ξ i j max ξ i j = [ ξ i j L max ξ i j U , ξ i j U max ξ i j L ] , j B
where B marks a set of benefit criteria, and
ς i j = min ξ i j ξ i j = [ min ξ i j L ξ i j U , min ξ i j U ξ i j L ] , j C
where C marks a set of cost criteria.
Step 3. Calculation of the weighted normalized Rough matrix (11) using Equation (12):
υ = [ [ υ 11 L , υ 11 U ] [ υ 21 L , υ 21 U ] [ υ m 1 L , υ m 1 U ] [ υ 12 L , υ 12 U ] [ υ 22 L , υ 22 U ] [ υ m 2 L , υ m 2 U ] [ υ 1 n L , υ 1 n U ] [ υ 2 n L , υ 2 n U ] [ υ m n L , υ m n U ] ]
υ i j = ς i j w j = [ ς i j L w j L , ς i j U w j U ]
Step 4. Calculation of an ideal τ i d (the largest value of υ i j ) and anti-ideal solution τ a i d (the smallest value υ i j ).
τ i d = [ τ i d L , τ i d U ] 1 × n = max υ i j
τ a i d = [ τ a i d L , τ a i d U ] 1 × n = min υ i j
Step 5. Determination of deviations from ideal and anti-ideal solutions, Equations (15) and (16).
δ i j + = [ δ i j + L , δ i j + U ] m × n = max τ i d υ i j = [ max τ i j L υ i j U , max τ i j U υ i j L ]
δ i j = [ δ i j L , δ i j U ] m × n = υ i j min τ a i d = [ υ i j L min τ i j U , υ i j U min τ i j L ]
Step 6. Calculation of the degree of deviation of each alternative from ideal and anti-ideal solutions.
i + = [ i + L , i + U ] 1 × m = j = 1 n δ i j +
i = [ i L , i U ] 1 × m = j = 1 n δ i j
Step 7. Determination of the utility function of each alternative relative to the optimal alternatives.
i + = [ i + L , i + U ] 1 × m = [ o + i + ] = [ o + L i + U , o + U i + L ] o + = j = 1 n ( min ( δ i j + ) )
i = [ i L , i U ] 1 × m = [ o i ] = [ o L i U , o U i L ] o = j = 1 n ( min ( δ i j ) )
Step 8. Sorting alternatives by calculating the average utility function using Equation (21). The best alternative is the one with the highest value.
χ i = [ χ i L , χ i U ] 1 × m = [ i + + i 2 ] = [ i + L + i L 2 , i + U + i U 2 ]

3. Integration of Linear Programming and a Fuzzy-Rough MCDM Model for Production Optimization

3.1. Special Case of Linear Programming Optimization

The company produces two articles, A and B, on two groups of machines, M and N. In the period considered, the first group of machines has a capacity of 6000 working hours, and the second group has a capacity of 3000 working hours. The processing time of a unit of article A is 1.5 working hours on the first group of machines, and 1 working hour on the second group of machines. The processing time of a unit of article B is 1.5 working hours on the first group of machines, and 0.5 working hour on the second group of machines. For the production realization, the company has enough raw materials and labor, but it can only put 2500 units of article A and 3000 units of article B on the market. By selling, the company makes a profit of 2000 monetary units per unit of article A and 1000 monetary units per unit of article B. It is necessary to determine the optimal production plan in order to achieve the maximum profit for the company.
First of all, it is necessary to set up a mathematical model in accordance with Equations (1)–(3). Constraints of the model according to the parameters of the study are as follows:
3 2 x 1 + 3 2 x 2 6000     x 1 + 1 2 x 2 3000 x 1 2500 x 2 3000 x 1 , x 2 0
while the objective function of maximizing the company’s profit is presented as follows:
F max = 2000 x 1 + 1000 x 2
By applying the simplex method of linear programming (Simplex Method calculator), the values of the objective function presented in Table 2 are obtained.
Non-basic variables are zero. The value of the basic variables can be taken into account. The function F contains only non-basic variables. Thus, the value of the function F for the basis can be kept in mind.
S 1 = 0 , S 2 = 0 , S 3 = 500 , S 4 = 1000 x 1 = 2000 , x 2 = 2000 F 6000000 = 0 F = 6000000
From a geometric point of view, both solutions are points of space, i.e., they form a line segment. Any point (any solution) on this segment will also be a solution.
Finally, the result of linear programming is presented by the following function:
x 1 = 2500 t + 2000 ( 1 t ) x 2 = 1000 t + 2000 ( 1 t ) 0 t 1 F max = 6000000
The optimal solution is represented by the extreme points, A and B, of the convex set C (Figure 2). In other words, the optimal solution is represented by every point of the line segment AB.
The maximum value of the objective function is as follows:
F max ( X ) = F ( A ) = 2000 ( 2000 ) + 1000 ( 2000 ) = 6000000
i.e.,
F max ( X ) = F ( B ) = 2000 ( 2500 ) + 1000 ( 1000 ) = 6000000
We will obtain the same value of the criterion function if we include the coordinates of any point on the line segment AB, which represents a large number of admissible optimal solutions. The optimal value of the objective function is unique because it is six million regardless of how many products A and how many products B will be produced, but there are infinitely many admissible solutions that provide this function value.
The multiplicity of the optimal solution gives greater opportunities to the decision-maker, and this is where MCDM methods come into focus. The decision-maker chooses an optimal solution that best fits these conditions, e.g., if the extreme point B is taken as the optimal solution, then the capacity of the second group of machines will be fully utilized, while the first group of machines will have an unused capacity of 750 working hours. The market will be fully supplied by the first article, and it will lack 2000 units of article B. Which solution the decision-maker will choose depends on whether it is more important for him to use the capacity of the machines or to supply the market with articles individually. It is through the further development of the multi-criteria Fuzzy-Rough MCDM model that the optimal solution is reached depending on the current needs of the company, and it is defined through determining the significance of the criteria and extensive sensitivity analysis.

3.2. Formation of MCDM Model—Defining Criteria and Alternatives

Since we have as a solution many points that represent the optimal solution, a multi-criteria model can be applied further. Four criteria are defined: capacity utilization of the first group of machines M (C1), capacity utilization of the second group of machines N (C2), supplying the market with articles of group A (C3) and supplying the market with articles of group B (C4). A total of 21 alternatives were defined in accordance with the value t = 0–1 with an interval of 0.05 (Table 3):

3.3. Determining the Significance of Criteria Using the IMF SWARA Method

In this section of the paper, the calculation of weight values of the criteria was carried out using the IMF SWARA method. It is not group decision-making because previously it was a specific optimization problem using linear programming, after which the decision-maker takes into account the current state and needs of the company and evaluates the criteria as shown in Table 4.
First, the linguistic values of the criterion comparison j ¯ are defined, and the values j ¯ for j > 1 are obtained by adding TFN (1,1,1). The values j ¯ for j > 1 are obtained as follows: 1 ¯ = 1 1 ¯ 1 ¯ = 4 ¯ 1 ¯ = ( 1 1.286 , 1 1.250 , 1 1.222 ) = ( 0.778 , 0.800 , 0.818 ) . The final fuzzy criterion weights are obtained as follows:
w 1 ¯ = 1 ¯ j = 1 n j ¯ = ( 1 3.488 , 1 3.440 , 1 3.383 ) = ( 0.287 , 0.291 , 0.296 ) .
The results show that the decision-maker views the significance of the criteria quite equally, so did not decide to assign much greater significance to any one criterion than to the others. Even the highest-ranked criterion is not much more significant compared to the least significant one.

3.4. Determining the Optimal Solution Using the Rough CRADIS Approach

In Section 3.2, alternative values of t and values of x1 and x2 are defined, by including the values of x1 and x2 in each constraint, as e.g., for A1:
3 2 × 2500 + 3 2 × 1000 = 5250 6000 ( 750 )       2500 + 1 2 × 1000 = 3000 ( 0 ) 2500 = 2500 ( 0 ) 1000 3000 ( 2000 )
For A2:
3 2 × 2475 + 3 2 × 1050 = 5287.5 6000 ( 712.5 )       2475 + 1 2 × 1050 = 3000 ( 0 ) 2475 2500 ( 25 ) 1050 3000 ( 1950 )
For A21:
3 2 × 2000 + 3 2 × 2000 = 6000 ( 0 )     2000 + 1 2 × 2000 = 3000 ( 0 ) 2000 2500 ( 500 ) 2000 3000 ( 1000 )
In this way, the complete initial decision matrix is obtained as shown in Table 5, which essentially represents the difference in satisfying the set restrictions for each criterion separately. For example, for A1, the values of x1 and x2 are 2500 and 1000, respectively, and the utilization of the first group of machines M is not complete, but there are still 750 h available, while the second group of machines N is fully utilized. When it comes to meeting the requirements of the market with product A, it is fully satisfied, while there is a shortage of 2000 products B according to these parameters. The initial decision-making matrix represents equal low and upper numbers because it represents the difference in satisfying the set restrictions for each criterion separately which is essential for managers while having no influence on the application of the rough model.
By Equations (8)–(21), the final results of applying the new Rough CRADIS approach are obtained, which is shown in Table 6.
After applying the Rough CRADIS model, the optimal solution is alternative A1, while the alternative A21 represents the second best solution under the given conditions and parameters of the model.

4. Verification of the Developed Model and Discussion

4.1. Sensitivity Analysis

The sensitivity analysis has been performed through 40 scenarios that represent new simulated criterion values using Equation (22) [22,23,24]:
W ¯ n β = ( 1 W ¯ n α ) W ¯ n = β ( 1 W ¯ n )
The new criterion values are shown in Figure 3, where the values are simulated as follows: starting from the first scenario in which the value of criterion C1 is reduced by 5% to S10 in which the value of that criterion is reduced by 95%, while the values of criteria C2, C3 and C4 increase proportionally. In the same way, new criterion values are determined in S11–S20 for C2, in S21–S30 for C3 and in S31–S40 for C4.
Figure 4 shows the ranks of all alternatives depending on the simulated weight values of the criteria.
The change in the ranks of alternatives is obvious so, e.g., A1 in S4–S10 is in the first position, which is a consequence of the reduction of C1 criterion by 35–95% while, e.g., in S27–S30, its rank decreases drastically to the 7th, 11th, 14th and 19th position, respectively, which is a consequence of the decrease in the significance of criterion C3 by 65–95%. This kind of sensitivity analysis in this special case of searching for an optimal solution represents an exceptional contribution for the decision-maker because he can fully perceive the real needs and determine which solution is optimal in real time.
Since there are certain changes in the ranks of potential solutions when reducing the weight values of the criteria, a statistical correlation of the ranks is performed (Figure 5), SCC and WS, which has already been mentioned in the methodology section of the paper. The lowest SCC = 0.521 is in set 40 since A3 drastically increases its position with a total decrease in the value of criterion C4. The lowest WS = 0.762 is in set 31, when the value of the third criterion decreases. Considering the large set of potential solutions (21), the average values of both coefficients show a very high level of rank correlation: SCC = 0.970, WS = 0.944.

4.2. Comparative Analysis

Comparative analysis represents one of the ways of testing the results obtained [25]. As already mentioned, the comparative analysis involves a comparison with four other Rough MCDM methods: Rough WASPAS [11], Rough SAW [12], Rough ARAS [13] and Rough MARCOS [14]. The results of the comparative analysis are shown in Figure 6.
Compared to three Rough MCDM methods, there is no change in the best alternative, while when applying the Rough ARAS method, the two best alternatives exchange places. Moreover, the greatest changes are in the comparison of Rough CRADIS with the Rough ARAS method, where A3 changes its position by several places.
Figure 7 shows the values of SCC and WS coefficients in the comparative analysis.
Rough CRADIS has a full correlation of ranks with Rough SAW, while it has 0.999 with Rough MARCOS, 0.874 with Rough ARAS, and 0.961 with Rough WASPAS when it comes to the SCC coefficient. The situation is similar with the WS coefficient, which has a full correlation with R-WASPAS, R-SAW, R-MARCOS.

4.3. Limitations and Managerial Implications

Limitations of the solution proposed by this research can be traced to linear programming and MCDM methodology. From the aspect of linear programming, the number of variables can represent limitations, while the developed Rough CRADIS model can be applied for group decision-making only, which is the second limitation.
It is very important to note that sensitivity analysis enables real-time decision-making since the decision-maker has at his disposal a set of 40 formed scenarios in which the criteria change their original values. This means that in real time, if it is more important for the company to satisfy the requirements of the market with product B rather than A, it can precisely determine from the obtained results what quantity of product it needs. The same is the case with the use of own resources, which are presented in the paper as criteria.

5. Conclusions

The paper proposes a novel integrated model consisting of a combination of linear programming, fuzzy set theory, rough set theory and a MCDM method. A novel Rough CRADIS approach has been created that can be used for any problem involving multiple criteria. The problem of optimizing production with the aim of maximizing profit, which represents a specific problem of linear programming, has been presented and solved. Based on a numerical example using linear programming, a set of optimal solutions that lie on the line segment AB has been obtained. After that, a total of 21 different solutions were formed, which can represent the optimum depending on the real requirements of the market and the orientation of the company. The IMF SWARA method was applied in order to determine the significance of the criteria considered when making a decision. At the end, a calculation was made using the novel Rough CRADIS approach in order to determine the optimal solution under the given conditions, i.e., how many products should be produced. The results show that the integration of linear programming and a Fuzzy-Rough MCDM model can be an exceptional solution for solving specific optimization problems.
Future research can be reflected through defining a larger number of influential factors and defining an objective function that will minimize costs while maximizing the profit. In addition, the application of fuzzy linear programming with the development of new approaches is highlighted as a useful option.

Author Contributions

Conceptualization, M.D., R.T. and Ž.S.; methodology, M.D. and Ž.S.; validation, D.K.D., S.T. and M.J.; writing—original draft preparation, M.D., R.T. and Ž.S.; writing—review and editing, D.K.D. and S.V.; project administration, S.T. and M.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bayá, G.; Canale, E.; Nesmachnow, S.; Robledo, F.; Sartor, P. Production Optimization in a Grain Facility through Mixed-Integer Linear Programming. Appl. Sci. 2022, 12, 8212. [Google Scholar] [CrossRef]
  2. Chandrawat, R.K.; Kumar, R.; Garg, B.P.; Dhiman, G.; Kumar, S. An analysis of modeling and optimization production cost through fuzzy linear programming problem with symmetric and right angle triangular fuzzy number. In Proceedings of the Sixth International Conference on Soft Computing for Problem Solving; Springer: Singapore, 2017; pp. 197–211. [Google Scholar]
  3. Pop, O.B.; Uifălean, D.C. The Impact of Innovation on the Management of Misapplication of Change in the Production System in the Automotive Industry. Econ. Innov. Econ. Res. J. 2022, 10, 199–210. [Google Scholar]
  4. Lazzerini, B.; Pistolesi, F. A linear programming-driven MCDM approach for multi-objective economic dispatch in smart grids. In Proceedings of the 2015 SAI Intelligent Systems Conference (IntelliSys), London, UK, 10–11 November 2015; pp. 475–484. [Google Scholar]
  5. Khalilzadeh, M.; Ghasemi, P.; Afrasiabi, A.; Shakeri, H. Hybrid fuzzy MCDM and FMEA integrating with linear programming approach for the health and safety executive risks: A case study. J. Model. Manag. 2021, 16, 1025–1053. [Google Scholar] [CrossRef]
  6. Malekpoor, H.; Chalvatzis, K.; Mishra, N.; Mehlawat, M.K.; Zafirakis, D.; Song, M. Integrated grey relational analysis and multi objective grey linear programming for sustainable electricity generation planning. Ann. Oper. Res. 2018, 269, 475–503. [Google Scholar] [CrossRef]
  7. Razak, N.H.; Hashim, H.; Yunus, N.A.; Klemeš, J. Integrated linear programming and analytical hierarchy process method for diesel/biodiesel/butanol in reducing diesel emissions. J. Clean. Prod. 2022, 337, 130297. [Google Scholar] [CrossRef]
  8. Cheng, S.; Chan, C.W.; Huang, G.H. An integrated multi-criteria decision analysis and inexact mixed integer linear programming approach for solid waste management. Eng. Appl. Artif. Intell. 2003, 16, 543–554. [Google Scholar] [CrossRef]
  9. Kannan, D.; Khodaverdi, R.; Olfat, L.; Jafarian, A.; Diabat, A. Integrated fuzzy multi criteria decision making method and multi-objective programming approach for supplier selection and order allocation in a green supply chain. J. Clean. Prod. 2013, 47, 355–367. [Google Scholar] [CrossRef]
  10. Sanayei, A.; Mousavi, S.F.; Abdi, M.R.; Mohaghar, A. An integrated group decision-making process for supplier selection and order allocation using multi-attribute utility theory and linear programming. J. Frankl. Inst. 2008, 345, 731–747. [Google Scholar] [CrossRef]
  11. Stojić, G.; Stević, Ž.; Antuchevičienė, J.; Pamučar, D.; Vasiljević, M. A novel rough WASPAS approach for supplier selection in a company manufacturing PVC carpentry products. Information 2018, 9, 121. [Google Scholar] [CrossRef]
  12. Durmić, E.; Stević, Ž.; Chatterjee, P.; Vasiljević, M.; Tomašević, M. Sustainable supplier selection using combined FUCOM–Rough SAW model. Rep. Mech. Eng. 2020, 1, 34–43. [Google Scholar] [CrossRef]
  13. Zavadskas, E.K.; Nunić, Z.; Stjepanović, Ž.; Prentkovskis, O. A novel rough range of value method (R-ROV) for selecting automatically guided vehicles (AGVs). Stud. Inform. Control 2018, 27, 385–394. [Google Scholar] [CrossRef]
  14. Vojinović, N.; Sremac, S.; Zlatanović, D. A Novel Integrated Fuzzy-Rough MCDM Model for Evaluation of Companies for Transport of Dangerous Goods. Complexity 2021, 2021, 1–16. [Google Scholar] [CrossRef]
  15. Kumar, S.; Maity, S.R.; Patnaik, L. Optimization of Wear Parameters for Duplex-TiAlN Coated MDC-K Tool Steel Using Fuzzy MCDM Techniques. Oper. Res. Eng. Sci. Theory Appl. 2022. [Google Scholar] [CrossRef]
  16. Sałabun, W.; Urbaniak, K. A new coefficient of rankings similarity in decision-making problems. In Proceedings of the International Conference on Computational Science, Amsterdam, The Netherlands, 3–5 June 2020; Springer: Cham, Switzerland, 2020; pp. 632–645. [Google Scholar]
  17. Vrtagić, S.; Softić, E.; Subotić, M.; Stević, Ž.; Dordevic, M.; Ponjavic, M. Ranking road sections based on MCDM model: New improved fuzzy SWARA (IMF SWARA). Axioms 2021, 10, 92. [Google Scholar] [CrossRef]
  18. Zolfani, S.H.; Görçün, Ö.F.; Küçükönder, H. Evaluating logistics villages in Turkey using hybrid improved fuzzy SWARA (IMF SWARA) and fuzzy MABAC techniques. Technol. Econ. Dev. Econ. 2021, 27, 1582–1612. [Google Scholar] [CrossRef]
  19. Puška, A.; Božanić, D.; Nedeljković, M.; Janošević, M. Green Supplier Selection in an Uncertain Environment in Agriculture Using a Hybrid MCDM Model: Z-Numbers–Fuzzy LMAW–Fuzzy CRADIS Model. Axioms 2022, 11, 427. [Google Scholar] [CrossRef]
  20. Puška, A.; Stević, Ž.; Pamučar, D. Evaluation and selection of healthcare waste incinerators using extended sustainability criteria and multi-criteria analysis methods. Environ. Dev. Sustain. 2022, 24, 11195–11225. [Google Scholar] [CrossRef]
  21. Puška, A.; Nedeljković, M.; Prodanović, R.; Vladisavljević, R.; Suzić, R. Market assessment of pear varieties in Serbia using fuzzy CRADIS and CRITIC methods. Agriculture 2022, 12, 139. [Google Scholar] [CrossRef]
  22. Božanić, D.; Milić, A.; Tešić, D.; Salabun, W.; Pamučar, D. D numbers–FUCOM–fuzzy RAFSI model for selecting the group of construction machines for enabling mobility. Facta Univ. Ser. Mech. Eng. 2021, 19, 447–471. [Google Scholar] [CrossRef]
  23. Mešić, A.; Miškić, S.; Stević, Ž.; Mastilo, Z. Hybrid MCDM solutions for evaluation of the logistics performance index of the Western Balkan countries. Econ.-Innov. Econ. Res. J. 2022, 10, 13–34. [Google Scholar] [CrossRef]
  24. Torkayesh, A.E.; Zolfani, S.H.; Kahvand, M.; Khazaelpour, P. Landfill location selection for healthcare waste of urban areas using hybrid BWM-grey MARCOS model based on GIS. Sustain. Cities Soc. 2021, 67, 102712. [Google Scholar] [CrossRef]
  25. Gergin, R.E.; Peker, İ.; Kısa, A.C.G. Supplier selection by integrated IFDEMATEL-IFTOPSIS Method: A case study of automotive supply industry. Decis. Mak. Appl. Manag. Eng. 2022, 5, 169–193. [Google Scholar] [CrossRef]
Figure 1. Research methodology.
Figure 1. Research methodology.
Axioms 11 00510 g001
Figure 2. The function F that has a maximum value at any point on the line segment AB.
Figure 2. The function F that has a maximum value at any point on the line segment AB.
Axioms 11 00510 g002
Figure 3. Simulation of the weight values of criteria through 40 newly formed scenarios.
Figure 3. Simulation of the weight values of criteria through 40 newly formed scenarios.
Axioms 11 00510 g003
Figure 4. Ranks of alternatives when changing the values of the criteria.
Figure 4. Ranks of alternatives when changing the values of the criteria.
Axioms 11 00510 g004
Figure 5. Statistical SCC and WS correlation coefficients for changes in criterion significance.
Figure 5. Statistical SCC and WS correlation coefficients for changes in criterion significance.
Axioms 11 00510 g005
Figure 6. Comparison results of R-CRADIS with four Rough MCDM methods.
Figure 6. Comparison results of R-CRADIS with four Rough MCDM methods.
Axioms 11 00510 g006
Figure 7. Correlation of ranks in comparative analysis.
Figure 7. Correlation of ranks in comparative analysis.
Axioms 11 00510 g007
Table 1. Linguistics and the TFN scale.
Table 1. Linguistics and the TFN scale.
Linguistic VariableAbbreviationTFN Scale
Absolutely less significantALS111
Dominantly less significantDLS1/22/31
Much less significantMLS2/51/22/3
Really less significantRLS1/32/51/2
Less significantLS2/71/32/5
Moderately less significantMDLS1/42/71/3
Weakly less significantWLS2/91/42/7
Equally significantES000
Table 2. Value of function F.
Table 2. Value of function F.
x1X2S1S2S3S4Const.
001−33/20750
0102−201000
1000102500
000−2212000
000−200000F − 6000000
002/3−210500
0102−201000
1000102500
000−2212000
000−200000F − 6000000
002/3−210500
014/3−2002000
10−2/32002000
00−4/32011000
000−200000F − 6000000
Table 3. Defining alternatives depending on the value t.
Table 3. Defining alternatives depending on the value t.
tx1x2
A1125001000
A20.09524751050
A30.09024501100
A40.08524251150
A50.08024001200
A60.07523751250
A70.07023501300
A80.06523251350
A90.06023001400
A100.05522751450
A110.05022501500
A120.04522251550
A130.04022001600
A140.03521751650
A150.03021501700
A160.02521251750
A170.02021001800
A180.01520751850
A190.01020501900
A200.0520251950
A21020002000
Table 4. Weight values of the criteria calculated with the IMF SWARA method.
Table 4. Weight values of the criteria calculated with the IMF SWARA method.
j ¯ j ¯ j ¯ w j ¯
C3 1.0001.0001.0001.0001.0001.0000.2870.2910.296
C4ES1.0001.0001.0001.0001.0001.0000.2870.2910.296
C1WLS1.2221.2501.2860.7780.8000.8180.2230.2330.242
C2WLS1.2221.2501.2860.6050.6400.6690.1730.1860.198
SUM3.3833.4403.488
Table 5. Initial decision matrix.
Table 5. Initial decision matrix.
C1C2C3C4
A1750.00750.00111120002000
A2712.50712.5011252519501950
A3675.00675.0011505019001900
A4637.50637.5011757518501850
A5600.00600.001110010018001800
A6562.50562.501112512517501750
A7525.00525.001115015017001700
A8487.50487.501117517516501650
A9450.00450.001120020016001600
A10412.50412.501122522515501550
A11375.00375.001125025015001500
A12337.50337.501127527514501450
A13300.00300.001130030014001400
A14262.50262.501132532513501350
A15225.00225.001135035013001300
A16187.50187.501137537512501250
A17150.00150.001140040012001200
A18112.50112.501142542511501150
A1975.0075.001145045011001100
A2037.5037.501147547510501050
A211.001.001150050010001000
Table 6. Results of applying the Rough CRADIS approach.
Table 6. Results of applying the Rough CRADIS approach.
i + i o + o i + i χ i Rank
A10.94650.66850.61650.99850.65140.6695420.6601
A21.33430.28070.46210.2811460.3723
A31.34220.27290.45940.2732710.3667
A41.34460.27040.45850.2708080.36511
A51.34570.26930.45810.2697090.36414
A61.34630.26880.45800.2691670.36417
A71.34650.26850.45790.2689130.36319
A81.34660.26840.45780.2688340.36321
A91.34660.26850.45790.2688760.36320
A101.34640.26860.45790.2690110.36318
A111.34620.26880.45800.2692220.36416
A121.34590.26910.45810.2695060.36415
A131.34560.26950.45820.2698620.36413
A141.34510.26990.45830.2702970.36412
A151.34460.27040.45850.2708290.36510
A161.34390.27110.45870.2714870.3659
A171.34310.27190.45900.2723370.3668
A181.34190.27310.45940.2735240.3666
A191.34000.27510.46010.2754850.3685
A201.33510.28000.46180.2803830.3714
A211.01990.59510.60450.5960120.6002
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dordevic, M.; Tešić, R.; Todorović, S.; Jokić, M.; Das, D.K.; Stević, Ž.; Vrtagic, S. Development of Integrated Linear Programming Fuzzy-Rough MCDM Model for Production Optimization. Axioms 2022, 11, 510. https://doi.org/10.3390/axioms11100510

AMA Style

Dordevic M, Tešić R, Todorović S, Jokić M, Das DK, Stević Ž, Vrtagic S. Development of Integrated Linear Programming Fuzzy-Rough MCDM Model for Production Optimization. Axioms. 2022; 11(10):510. https://doi.org/10.3390/axioms11100510

Chicago/Turabian Style

Dordevic, Milan, Rade Tešić, Srdjan Todorović, Miloš Jokić, Dillip Kumar Das, Željko Stević, and Sabahudin Vrtagic. 2022. "Development of Integrated Linear Programming Fuzzy-Rough MCDM Model for Production Optimization" Axioms 11, no. 10: 510. https://doi.org/10.3390/axioms11100510

APA Style

Dordevic, M., Tešić, R., Todorović, S., Jokić, M., Das, D. K., Stević, Ž., & Vrtagic, S. (2022). Development of Integrated Linear Programming Fuzzy-Rough MCDM Model for Production Optimization. Axioms, 11(10), 510. https://doi.org/10.3390/axioms11100510

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