Abstract
FUZZY INFORMATION AND ENGINEERING 2021, VOL. 13, NO. 3, 323–334 https://doi.org/10.1080/16168658.2021.1938868 A New Method to Solve Multi-Objective Linear Fractional Problems Mojtaba Borza and Azmin Sham Rambely Department of Mathematical Sciences, Faculty of Science & Technology, UKM Bangi, Selangor, Malaysia ABSTRACT ARTICLE HISTORY Received 23 December 2020 Background: In the literature, there exists several approaches to Revised 9 May 2021 address the multi-objective linear fractional programming problem Accepted 1 June 2021 (MOLFPP). However, there is a drawback to these methods. Aim: This paper presents an efficient method treating the MOLFPP. KEYWORDS Methodology: To construct our approach,the membership func- Eﬃcient solution; tions of the objectives, suitable non-linear variable transformations, membership function; and max-min technique are used. max–min technique; linear programming; fractional Results: In our proposed method, the MOLFPP is finally changed into programming a linear programming problem (LPP). It is proven that the optimal solution of the LPP is an efficient solution for the MOLFPP. Conclusion: Numerical examples are solved, and the results demon- strate that our method with less computational expenses and cost reach the efficient solutions. 1. Introduction The concept of a multi-objective programming problem (MOPP) arises when a decision maker is going to consider more than one objective over a common set of restrictions. A number of real-world problems in transportation, finance, engineering, commercials, house planning, energy systems, etc. have appropriately been modelled as MOPPs [1]. For this class of problems, the concept of efficient solution is considered instead of exact opti- mal solution. A solution is efficient if moving to another solution does not improve all the objectives. In MOPP, if the objectives are linear fractional functions and the constraints are affine, then this model represents the multi-objective linear fractional programming prob- lem (MOLFPP). In Ref. [2], applications of a linear fractional programming problem (LFPP) in economy, business, engineering, management, etc. were demonstrated. Radhakrishnan and Anukokila [3] addressed a solid transportation problem with interval cost by the use of a fractional goal programming method. Wang et al. [4] developed a framework of bi- level MOLFPP to optimise a water consumption structure. Ahmad et al. [5] investigated the fractional-order tumour-immune-vitamin model trough fixed point results. Das et al. [6] presented an application of the LFPP with fuzzy nature in industry sector. CONTACT Azmin Sham Rambely asr@ukm.edu.my This article has been republished with minor changes. These changes do not impact the academic content of the article. © 2021 The Author(s). Published by Taylor & Francis Group on behalf of the Fuzzy Information and Engineering Branch of the Operations Research Society, Guangdong Province Operations Research Society of China. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons. org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 324 M. BORZA AND A. S. RAMBELY Because of the importance of the LFPP and also MOLFPP, many studies have been accomplished to come out with efficient methods and techniques for these optimisation problems. Chakraborty and Gupta [7] developed a method to address MOLFPP. In their method, the multi-objective problem is transformed into a multi-objective linear program- ming problem (MOLPP). Subsequently, the membership functions are specified after identi- fying the fuzzy aspiration levels of the linear objectives. Finally, the MOLPP is changed into a linear programming problem (LPP) using a max–min technique. Motivated by Chakraborty and Gupta’s methodology, Veeramani and Sumathi [8] and De and Deb [9] introduced approaches to deal with LFPP with fuzzy coefficients and MOLFPP, respectively. Follow- ing the methodology of Dinkelbach [10], Güzel [11] and Nayak and Ojha [12] developed approaches to MOLFPP. In fact, Nayak and Ojha attempted to improve the results of Guzel’s approach by employing the ε-constrain technique. However, applying the ε-constrain tech- nique encompasses some difficulties in practice when the decision maker is trying to specify the value of ε for each constraint. Pal et al. [13] transformed the MOLFPP into a LPP using a fuzzy goal programming approach in addition to suitable variable transformations. Tok- sari [14] introduced an approach to tackle the MOLFPP where the membership functions of the objectives are defined and then linearised using the first-order Taylor series about the individual optimal solutions. For some examples, Borza et al. [15] reported that the results of using the first-order Taylor series proposed by Toksari are to some extent more accurate than the results of the fuzzy goal programming used by Pal et al. Nayak and Ojha [16]intro- duced a method dealing with the MOLFPP with fuzzy coefficients where the fuzzy problem is altered into interval valued LFPP using the concept of α-cuts. In their method, the fuzzy problem is reduced into the MOLFPP. Afterwards, they reach a MOLPP employing the first- order Taylor series. Finally, weighted sum technique is utilised to transform the MOLPP into a LPP. Borza and Rambely [17] designed a non-iterative method to obtain the global opti- mal solution of the sum of the linear fractional programming problem (S-LFPP) by the use of variable transformation. Liu et al. [18] constructed an iterative algorithm for the large-scale S-LFPP using a branch and bound technique. In the literature, many researchers have tried to transform the MOLFPP into a LPP using different methodology and techniques such as the first-order Taylor series method, Dinkel- bach’s methodology, and Chakraborty and Guptas’ approach. However, there are draw- backs regarding these methods and methodologies. Using the first-order Taylor expansion reduces the accuracy of the method automatically. The method of Chakraborty and Gupta was designed in such a way that it has not been possible to prove their methodology results in efficient solutions. Following the methodology of Dinkelbach, a fractional pro- gramme is changed into a parametric non-fractional programming problem. However, in the existing methods, a non-parametric model of Dinkelbach’s methodology has been used, which reduces the accuracy. In this paper, we aim to present a new efficient and straightforward method with less computational expenses and appropriate accuracy to transform the MOLFPP into a LPP. In addition, we use the membership functions of the objectives to construct our approach in order to a wide range of problems be covered. Therefore, the membership functions of the objectives are specified after identifying the maxima and minima of the objectives and then a new MOLFPP is designed. This new prob- lem is changed into a MOLPP by the use of non-linear variable transformations. Finally, the max–min approach is used to tackle the MOLPP. It is proven that the solution resulted is efficient for the MOLFPP. Numerical examples are given to illustrate the method in addition to make comparison to some existing methods. FUZZY INFORMATION AND ENGINEERING 325 This article is organised in four sections. Following the introduction, in Section 2, some preliminaries are given. In Section 3, the main result and outcome of this survey are released. In Section 3, numerical examples are solved to illustrate the method and make comparison. Finally, Section 4 concludes the study. 2. Preliminaries 2.1. Linear Fractional Programming Consider the general form of the LFPP as follows: C X + α Maximize D X + β (1) s.t. AX ≤ b, X ≥ 0, D X + β> 0. According to the method introduced by Charnes and Cooper [19], Equation (1) is changed into the following linear problem by the use of variable transformations t = (1/(D X + β)), Y = tX. Maximize C Y + αt (2) s.t. AY − bt ≤ 0, D Y + βt = 1, Y, t ≥ 0. ∗ ∗ ∗ ∗ ∗ Theorem 1: (Ref. [19]). Let (Y , t ) be the optimal solution of (2), then X = (Y /t ) is optimum for (1). 2.2. Multi-Objective Programming Let us consider the general form of the MOPP as follows: Maximize {F (X), ... , F (X)} s.t. X ∈ S.(3) 1 k Definition 1: (Ref. [20]). For (3), a solution X ∈ S is called efficient if and only if X ∈ S such ∗ ∗ that F (X ) ≤ F (X),j = 1, ... , k, and ∃l ∈{1, ... , k} such that F (X )< F (X). j j j j 3. Main Results In this section, an approach is introduced in order to change the MOLFPP into a LPP such that the optimal solution of the LPP becomes an efficient solution for the MOLFPP. Consider the general type of the MOLFPP as follows: N X + m Maximize Z (X) = for i = 1, ... , k (4) P X + q s.t. S ={AX ≤ b, X ≥ 0}, where S is a regular set (non-empty and bounded set). Furthermore, for X = (X , ... , X ) ∈ 1 n S, it is assumed P X + q > 0for i = 1, ... , k. Our aim is to design a method so as to come out with efficient solution for (4). To do this, we need N X + m ≥ 0, ∀X ∈ S, i = 1, ... , k. But, these conditions are restrictive. To i 326 M. BORZA AND A. S. RAMBELY overcome this difficulty, an equivalent problem to (4) is constructed in which numerators are non-negative. Therefore, the membership functions of the objectives are specified and then are utilised instead of the objectives. max min T T Let max Z = z and min Z = z , i = 1, ... , k.Thus, μ (X) = (C X + d )/(P X + X∈S i i i i i i i i X∈S max min q ) is the membership function for objective Z , where X ∈ S, C = (1/(z − z ))N − i i i i i i min max min min T z P,and d = (m /z − z ) − z q . Accordingly, C X + d ≥ 0because μ (X) ∈ i i i i i i i i i i i [0, 1], and P X + q > 0for i = 1, ... , k. The equivalent of (4) in terms of the membership functions is C X + d Maximize for i = 1, ... , k P X + q (5) s.t. X ∈ S ={AX ≤ b, X ≥ 0}. Let us define new variables λ and Y as the functions of variable X as follows: λ = min{λ = (1/P X + q ), i = 1, ... , k} and λX = Y.(6) i i Thus, (5) is transformed into Maximize{C Y + λd for i = 1, ... , k} s.t. F ={AY − λb ≤ 0, Y, λ ≥ 0, P Y + λq ≤1for i = 1, ... , k}.(7) Lemma 1: In (7), variable λ = 0, ∀(Y, λ) ∈ F. ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ Proof: If (Y,0) ∈ F, thenAY ≤ 0. Now, if X ∈ S, then A(X + βY) = AX + β(AY) ≤ AX ≤ b ˆ ˆ for all β ≥ 0; this means X + βY ∈ S, ∀β ≥ 0. This results that the feasible region S is an unbounded set, which is a contradiction to the regularity of S. ¯ ¯ ¯ ¯ Lemma 2: If (Y, λ) ∈ F, then (Y/λ) ∈ S. ¯ ¯ ¯ ¯ ¯ ¯ ¯ ¯ ¯ Proof: Since (Y, λ) ∈ F, then Y ≥ 0, λ> 0, and AY − λb ≤ 0. Therefore, (1/λ)(AY − λb) = ¯ ¯ A(Y/λ) − b ≤ 0. By setting β = min{C Y + λd , i = 1, ... , k}, (7) is altered into X∈F Maximize β T T s.t. φ ={Y, λ, β ≥ 0, AY − λb ≤ 0, β ≤ C Y + λd , P Y + λq ≤1for i = 1, ... , k}.(8) i i i i Theorem 2: The optimal solution of (8) is unique. ∗ ∗ ∗ Proof: Let (Y , λ , β ) be the optimal solution and is not unique; this means constraint β ≥ 0 is active at the optimum, i.e. β = 0. In other words, if (Y, λ, β) ∈ φ, then β = 0.Therefore, T T ∃j ∈{1, ... , k} such that C Y + λd = 0 for all (Y, λ,0) ∈ φ.Since λ> 0, then C X + d = 0 j j j j for all X ∈ S; this means μ (X) = 0 for all X ∈ S. As the consequence, (5) is reduced into (k − 1) objective LFPP. This is a contradiction. ∗ ∗ ∗ ∗ ∗ ∗ Theorem 3: If (Y , λ , β ) is optimal for (8), then X = (Y /λ ) is an efficient solution for (5). FUZZY INFORMATION AND ENGINEERING 327 ∗ ∗ ∗ Proof: Let X = (Y /λ ) not be an efficient solution for (5). Therefore, ∃X ∈ S such that T ∗ T C (X ) + d C X + d i i i i ≤ for i = 1, ... , k,and T T P (X ) + q P X + q i i i i T ∗ T C (X ) + d C X + d j j j j ∃j ∈{1, ... , k} such that < .(9) T T P (X ) + q P X + q j j j j Consider ∗ ∗ ∗ ∗ ∗ ∗ T ∗ ∗ (Y , λ , β ) ∈ φ ⇒ λ ≤ λ = and 0 ≤ β ≤ C Y + λ q , i = 1, ... , k. i i P (X ) + q (10) ¯ ¯ ¯ ¯ ¯ Let us define θ = max{λ = (1/(P X + q )), i = 1, ... , k} and λ = θ − , where i i T ∗ C X + d ∗ i ¯ ¯ ¯ θ − λ ≤ ≤ θ − λ , i = 1, ... , k. (11) C X + d We need to show that (11) is well defined. In other words, there must exist satisfying (11). To do this, two below conditions must hold true. (I) C X + d = 0, i = 1, ... , k. T ∗ C X +d ¯ ¯ ¯ (II) θ − λ ≤ θ − λ , i = 1, ... , k. C X+d T T T T Since μ (X) = ((C X + d )/(P X + q )) ∈ [0, 1], P X + q > 0, then C X + d ≥ 0, ∀X ∈ S. i i i i i i i i i T T ∗ Now, let ∃j ∈{1, ... , k} such that C X + d = 0. Due to (9), it is possible that ((C (X ) + j j T ∗ T T ∗ ¯ ¯ d )/(P (X ) + q )) < ((C X + d )/(P X + q )); this means μ (X )< 0. This contradicts the j j j j j j j j non-negativity of membership functions. Therefore, (I) is verified. It follows directly from (9) and (10) that T T ¯ ¯ C X + d C X + d i i ∗ T ∗ ∗ T ∗ i i T ¯ ¯ λ (C (X ) + d ) ≤ λ (C X + d ) = ≤ = λ (C X + d ), (12) i i i i i i i i T T ¯ ¯ P X + q P X + q i i i i T ∗ C ( X ) + d ∗ i λ ≤ λ . (13) C X + d ∗ T ∗ T ¯ ¯ ¯ ¯ (13) ⇒ θ − λ ≤ θ − λ (C (X ) + d /C X + d ), i = 1, ... , k. Therefore, (II) is demon- i i i i i strated. It is time to show: ¯ ¯ (III) λ(P X + q ) ≤ 1, i = 1, ... , k. ∗ T ∗ T ¯ ¯ (IV) λ (C (X ) + d ) ≤ λ(C X + d ), i = 1, ... , k. i i i i To do this: ¯ ¯ ¯ ¯ (11) implies θ − ≤ λ . Furthermore, according to the definitions θ = max{λ,for i = i i T T ¯ ¯ ¯ ¯ ¯ ¯ 1, ... , k}, λ = θ − ,and λ = (1/P X + q ), i = 1, ... , k,itisconcludedthat λ(P X + q ) = i i i i i T T ¯ ¯ ¯ ¯ (θ − )(P X + q ) ≤ λ (P X + q ) = 1, i = 1, ... , k.Thus, (III) is demonstrated. i i i i i 328 M. BORZA AND A. S. RAMBELY ∗ T ∗ T ∗ T ∗ T ∗ T ∗ ¯ ¯ ¯ ¯ Since ≤ θ − λ (C X + d /C X + d ), then λ (C X + d /C X + d ) ≤ θ − ⇒ λ (C X i i i i i i i i i T T ¯ ¯ ¯ ¯ + d ) ≤ (θ − )(C X + d ) = λ(C X + d ), i = 1, ... , k.Thus, (IV) is verified. i i i i i ¯ ¯ ¯ ¯ ¯ Now, let us define Y = λX.Toshow (Y, λ) ∈ F, the followings must be true: (a) λ ≥ 0. T T ∗ ∗ ¯ ¯ ¯ Due to (11), let us set max = max{θ − λ (C X + d /C X + d ), i = 1, ... , k}= θ − i i i i ∗ T ∗ T ∗ T ∗ T ∗ ¯ ¯ ¯ ¯ ¯ ¯ λ (C X + d /C X + d ).Thus, λ ≥ θ − max = θ − (θ − λ (C X + d /C X + d )) = λ l l l l l l l l T ∗ T (C X + d /C X + d ) ≥ 0. l l l l (b) Y ≥ 0. ¯ ¯ ¯ ¯ ¯ Since X ∈ S, then X ≥ 0. Consequently, Y = λX ≥ 0. ¯ ¯ (c) (P Y + λq ) ≤1for i = 1, ... , k. ¯ ¯ ¯ Considering Y = λX and (III) proves c. ¯ ¯ (d) AY − λb ≤ 0. ¯ ¯ ¯ ¯ ¯ ¯ X ∈ S ⇒ AX − b ≤ 0. Therefore, AY − λb = λ(AX − b) ≤ 0. ¯ ¯ ¯ ¯ ¯ In what follows, we create β such that β ≥ β and (Y, λ, β) ∈ φ. (IV) ⇒ T ∗ ∗ ∗ T ∗ T T ¯ ¯ ¯ ¯ C Y + λ d = λ (C (X ) + d ) ≤ λ(C X + d ) = C Y + λd , i = 1, ... , k. (14) i i i i i i i i (10) and (14) ⇒ ∗ T ¯ ¯ 0 ≤ β ≤ C Y + λd ,i = 1, ... , k. (15) Let us set T ∗ ∗ ¯ ¯ ¯ γ = min{C Y + λd − β , i = 1, ... , k} and β = β + γ . (16) (15) and (16) ⇒ γ ≥ 0, and subsequently β ≤ β. (17) ∗ T ¯ ¯ It follows directly from (16) that 0 ≤ γ + β ≤ C Y + λd , i = 1, ... , k.Thus, ¯ ¯ ¯ 0 ≤ β ≤ C Y + λd , i = 1, ... , k. (18) ¯ ¯ ¯ ¯ ¯ Equation (18) in addition to (Y, λ) ∈ F results (Y, λ, β) ∈ φ. ¯ ¯ ¯ ¯ In brief, we found (Y, λ, β) ∈ φ such that β ≤ β. This contradicts the unique optimality ∗ ∗ ∗ of (Y , λ , β ) for (8). The proof is then complete. 4. Numerical Example In this section, four examples are considered taken from different references in order to illus- trate and evaluate this method. The third and fourth examples are mathematical models of the real-world organisations problems. FUZZY INFORMATION AND ENGINEERING 329 4.1. Example 1 (Ref. [14]) 12X + 13X 12X + 13X 1 2 1 2 Maximize Z (X) = , Z (X) = 1 2 40X + 55X + 500 1.5X + 1.6X 1 2 3 4 s.t. S ={2X + X ≤ 250, 5X + 4X ≤ 500, 45X + 30X ≤ 1500, 1 2 1 2 1 2 0.1X + 0.1X − X − X ≤ 0, 0.1X − X ≤ 0, 0.05X − X ≤ 0, 1 2 3 4 1 3 2 4 − X + X ≤ 0, − X + X ≤ 0, X , X , X , X ≥ 0}. (19) 1 3 2 4 1 2 3 4 max min First, the values of z and z for i = 1, 2 are individually determined by the i i max min use of Ref. [19] so as to define the membership functions: z = 0.2182, z = 1 1 max min 0, z = 83.6735, z = 8 . Thus, μ (X) = (54.9954X + 59.579X /40X + 55X + 500) Z 1 2 1 2 2 2 1 and μ (X) = (0.1586X + 0.1718X − 0.1586X − 0.1691X /1.5X + 1.6X ). Z 1 2 3 4 3 4 Equation (8) is formulated for (19) as follows: Maximize β s.t. {2Y + Y − 250λ ≤ 0, 5Y + 4Y − 500λ ≤ 0, 1 2 1 2 45Y + 30Y − 1500λ ≤ 0, 0.1Y + 0.1Y − Y − Y ≤ 0, 1 2 1 2 3 4 0.1Y − Y ≤ 0, 0.05Y − Y ≤ 0, 1 3 2 4 − Y + Y ≤ 0, − Y + Y ≤ 0, 1 3 2 4 40Y + 55Y + 500λ ≤ 1, 1.5Y + 1.6Y ≤ 1, 1 2 3 4 β ≤ 54.9954Y + 59.579Y , β ≤ 0.1586Y + 0.1718Y − 0.1586Y − 0.1691Y , 1 2 1 2 3 4 Y , Y , Y , Y , λ, β ≥ 0}.. (20) 1 2 3 4 ∗ ∗ ∗ Equation (20) is solved and the unique solution obtained is (Y , λ , β ) = (0.0008, 0.0146, ∗ ∗ ∗ 0.0008, 0.0007, 0.0003, 0.0029). Furthermore, the solution for (19) is X = (Y /λ ) = (2.5642, 46.1537, 2.5641, 2.3077). At the solution X , Z (X) = 0.2008, Z (X) = 83.6735, μ (X) = 0.9203, and μ (X) = 1. 1 2 Z Z 1 2 The average of μ (X) and μ (X) is 0.9602. Z Z 1 2 4.1.1. Comparison The solution resulted by Toksari is X = (0, 50, 0, 5). At the solution X, Z (X) = 0.2, Z (X) = 81.25, μ (X) = 0.9166, and μ (X) = 0.9669. 1 2 Z Z 1 2 The average of μ (X) and μ (X) is 0.9417. Z Z 1 2 As we see, the solution of Toksari is dominated by our proposed solution, i.e. ∗ ∗ ˆ ˆ Z (X)< Z (X ), Z (X)< Z (X ). 1 1 2 2 4.2. Example 2 (Ref. [7]) −3X + 2X 7X + X 1 2 1 2 Maximize Z (X) = , Z (X) = 1 2 X + X + 3 5X + 2X + 1 1 2 1 2 s.t. S ={−X + X ≤−1, 2X + 3X ≤ 15, −X ≤−3, X , X ≥ 0}. (21) 1 2 1 2 1 1 2 330 M. BORZA AND A. S. RAMBELY max min max min For (21), Z =−0.6087, Z =−2.1429, and Z = 1.3636, Z = 1.148. Accord- 1 1 2 2 ingly, μ (X) = (−0.5587X + 2.7004X + 4.1903/X + X + 3) and μ (X) = (5.8473X − Z 1 2 1 2 Z 1 1 2 6.041X − 5.3482/5X + 2X + 1). 2 1 2 Equation (8) is formulated for (21) as follows: Maximize β s.t. φ ={−Y + Y + λ ≤ 0, 2Y + 3Y − 15λ ≤ 0, −Y + 3λ ≤ 0, 1 2 1 2 1 Y + Y + 3λ ≤ 1, 5Y + 2Y + λ ≤ 1, 1 2 1 2 β ≤− 0.5587Y + 2.7004Y + 4.1903λ, β ≤ 5.8473Y − 6.041Y − 5.3482λ, 1 2 1 2 Y , Y , λ, β ≥ 0}. (22) 1 2 Equation (22) is solved and the unique optimal solution obtained is (0.1647, 0.0608, 0.0549, 0.3022). Thus, the solution proposed for problem (21) is X = = (3, 1.1073). At the solution X , Z (X) =−0.9547, Z (X) = 1.2137, μ (X) = 0.7746, and μ (X) = 0.3022. 1 2 Z Z 1 2 The average of μ (X) and μ (X) is 0.5384. Z Z 1 2 4.2.1. Comparison The solution of Chakraborty and Gupta is X = (3, 2). At the solution X, Z (X) =−0.625, Z (X) = 1.15, μ (X) = 0.9894, and μ (X) = 0.0056. 1 2 Z Z 1 2 The average of μ (X) and μ (X) is 0.4975. Z Z 1 2 As we observe, the solution of Chakrabory and Gupta does not dominate our proposed solution and vice versa. However, the average of membership functions shows that our proposed method has a better efficiency and function. 4.3. Example 3 (Ref. [21]) Maximize {Z (X), Z (X)} 1 2 59890X + 23390X + 30750X + 59750X + 40700X + 59435X 1 2 3 4 5 6 = , 35345X + 13420X + 18455X + 39455X + 23840X + 24070X + 500000 1 2 3 4 5 6 59890X + 23390X + 30750X + 59750X + 40700X + 59435X 1 2 3 4 5 6 96X + 120X + 144X + 144X + 84X + 120X + 480 1 2 3 4 5 6 s.t. S ={0.3X + 0.4X + 0.4X + 0.98X + 0.97X + 0.98X ≤ 600, 1 2 3 4 5 6 2280000X + 9200X + 16000X + 22500X + 20000X + 20000X ≤ 20000000, 1 2 3 4 5 6 650X + 630X + 320X + 660X + 360X + 640X ≤ 500000, 1 2 3 4 5 6 20X + 22X + 20X + 18X + 20X + 17X ≤ 15000, 1 2 3 4 5 6 FUZZY INFORMATION AND ENGINEERING 331 11400X + 3220X + 1800X + 12750X + 3250X + 3000X ≤ 6000000, 1 2 3 4 5 6 148X + 238X + 135X ≤ 50000, 1 4 6 180X + 220X + 200X + 150X + 100X + 160X ≤ 120000, 1 2 3 4 5 6 60X + 40X + 35X + 50X + 30X + 45X ≤ 30000, 1 2 3 4 5 6 30X + 32X + 28X + 35X + 26X + 20X ≤ 200000, 1 2 3 4 5 6 15X + 18X + 16X + 14X + 17X + 18X ≤ 10000, 1 2 3 4 5 6 42X + 38X + 36X + 40X + 37X + 35X ≤ 25000, 1 2 3 4 5 6 X ≥ 0, i = 1, ... ,6}. (23) min max min max For (23), Z = 0, Z = 2.3381, Z = 0, Z = 491.5151, 1 1 2 2 25615X + 10004X + 13152X + 25555X + 17407X + 25420X 1 2 3 4 5 6 μ (X) = , 35345X + 13420X + 18455X + 39455X + 23840X + 24070X + 500000 1 2 3 4 5 6 121.8477X + 47.5876X + 62.5617X + 121.5629X + 82.8052X + 120.992X 1 2 3 4 5 6 μ (X) = . 96X + 120X + 144X + 144X + 84X + 120X + 480 1 2 3 4 5 6 Equation (8) is formed for the above problem as follows: Maximize β s.t. φ ={0.3Y + 0.4Y + 0.4Y + 0.98Y + 0.97Y + 0.98Y − 600λ ≤ 0, 1 2 3 4 5 6 2280000Y + 9200Y + 16000Y + 22500Y + 20000Y + 20000Y − 20000000λ ≤ 0, 1 2 3 4 5 6 650Y + 630Y + 320Y + 660Y + 360Y + 640Y − 500000λ ≤ 0, 1 2 3 4 5 6 20Y + 22Y + 20Y + 18Y + 20Y + 17Y − 15000λ ≤ 0, 1 2 3 4 5 6 11400Y + 3220Y + 1800Y + 12750Y + 3250Y + 3000Y − 6000000λ ≤ 0, 1 2 3 4 5 6 148Y + 238Y + 135Y − 50000λ ≤ 0, 1 4 6 180Y + 220Y + 200Y + 150Y + 100Y + 160Y − 120000λ ≤ 0, 1 2 3 4 5 6 60Y + 40Y + 35Y + 50Y + 30Y + 45Y − 30000λ ≤ 0, 1 2 3 4 5 6 30Y + 32Y + 28Y + 35Y + 26Y + 20Y − 200000λ ≤ 0, 1 2 3 4 5 6 15Y + 18Y + 16Y + 14Y + 17Y + 18Y − 10000λ ≤ 0, 1 2 3 4 5 6 42Y + 38Y + 36Y + 40Y + 37Y + 35Y − 25000λ ≤ 0, 1 2 3 4 5 6 35345Y + 13420Y + 18455Y + 39455Y + 23840Y + 24070Y + 500000λ ≤ 1, 1 2 3 4 5 6 96Y + 120Y + 144Y + 144Y + 84Y + 120Y + 480λ ≤ 1, 1 2 3 4 5 6 β ≤ 25615Y + 10004Y + 13152Y + 25555Y + 17407Y + 25420Y , 1 2 3 4 5 6 β ≤ 121.8477Y + 47.5876Y + 62.5617Y + 121.5629Y + 82.8052Y + 120.0492Y , 1 2 3 4 5 6 Y ≥ 0, i = 1, ... ,6, λ, β ≥ 0}. (24) ∗ ∗ ∗ Equation (24) is solved and the solution X = (Y /λ ) = (0, 0, 0, 0, 0, 370) is obtained as an efficient solution for (23). 332 M. BORZA AND A. S. RAMBELY At the solution X , Z (X) = 2.3380, Z (X) = 489.9944, μ = 0.9999, and μ = 0.9897. 1 2 Z Z 1 2 The average of μ (X) and μ (X) is 0.9948. Z Z 1 2 4.3.1. Comparison The solution proposed by Pramy and Islam is X = (0, 0, 0, 0, 196.078, 370.37). At the solution X, Z (X) = 2.1288, Z (X) = 488.531, μ = 0.9105, and μ = 0.9887. 1 2 Z Z 1 2 The average of μ (X) and μ (X) is 0.9496. Z Z 1 2 The results show that our solution X dominates the solution X proposed by Pramy and Islam due to the fact that ∗ ∗ ˆ ˆ Z (X)< Z (X ), Z (X)< Z (X ). 1 1 2 2 4.4. Example 4 (Ref. [22]) In this section, a real life production planning in Taiwan is considered. The original prob- lem is modelled as a LFPP with fuzzy coefficients and fuzzy decision variables. In order to be able to solve the problem with the method provided, we change the fuzzy numbers into the intervals using the concept of α-cuts Moreover, the decision variables are set to be non-fuzzy. Therefore, we transformed the problem into the MOLFPP by the use of interval operations as follows: Maximize{Z (X), Z (X)} 1 2 ⎪ 9.2X + 21.4X + 9.2X + 19.5X + 14.6X + 19.3X + 11.2X 1 2 3 4 5 6 7 +7.2X + 19.4X + 11X + 9.1X + 14.6X 8 9 10 11 12 = , 2.2X + 5.4X + 2.2X + 4.4X + 3.4X + 3.4X + 3.4X ⎪ 1 2 3 4 5 6 7 +2.2X + 4.4X + 3.4X + 2.5X + 3.4X 8 9 10 11 12 9.7 + 22.8X + 9.76X + 20.8X + 15.4X + 20.8X ⎪ 2 3 4 5 6 +12.4X + 8.32X + 20.4X + 12.4X + 10.32X + 15.4X 7 8 9 10 11 12 1.8X + 4.6X + 1.72X + 3.6X + 2.8X + 2.6X 1 2 3 4 5 6 ⎪ +2.7X + 1.8X + 3.6X + 2.6X + 1.8X + 2.6X 7 8 9 10 11 12 s.t. S ={X + X + X + X ≤ 8.32, X + X + X + X ≤ 14.8, 1 2 3 4 5 6 7 8 X + X + X + X ≤ 12.72, 9 10 11 12 X + X + X ≥ 7.32, X + X + X ≥ 10.44, X + X + X ≥ 8.6, X + X + X ≥ 9.48, 1 5 9 2 6 10 3 7 11 4 8 12 X ≥ 0, i = 1, ... ,12}. (25) Equation (25) is solved by the proposed method and the solution obtained is X = (0.92, 0, 7.4, 1.56, 13.24, 0, 0, 4.2, 0, 0, 8.52). FUZZY INFORMATION AND ENGINEERING 333 At the solution X , Z (X) = 4.8271, Z (X) = 6.6052, 1 2 μ (X) = 1, μ (X) = 0.9660. Z Z 1 2 The average of μ (X) and μ (X) is 0.983. Z Z 1 2 As we observe, our proposed method addressed (25) in an excellent way since the average of the membership functions is very close to one. It is noticeable that the genetic algorithm of the global optimisation toolbox of MATLAB R2016 failed to reach a solution for this example. 5. Conclusion In this paper, a new method was presented to solve the MOLFPP. In the approach, the MOLFPP was changed finally into a LPP using suitable non-linear variable transformations. It was proven that the optimal solution of the LPP is unique and is efficient for the MOLFPP. We need to mention that the proposed method is easy and straightforward with less computa- tional complexities compared to the other existing methods. Moreover, this approach can be applied to address the LFPP with fuzzy coefficients if the fuzzy coefficients are changed into intervals using the concept of α-cuts. In this case, the fuzzy problem is further changed into a bi-objective LFPP. Four examples were solved to illustrate the approach in addition to make comparisons. For numerical examples, our proposed solutions gave better outcomes compared to Tok- sari, Chakraborty and Gupta, and Pramy and Islam. Furthermore, the results demonstrate that the method of Chakraborty and Gupta is reliable, but we cannot consider the methods of Toksari and Pramy and Islam as the effective approaches since their solutions proposed for Examples 1 and 3 were completely dominated by our proposed solutions. As a future research, one can employ the results of this study to cope with multi-level MOLFPP. Disclosure statement No potential conflict of interest was reported by the author(s). Funding This work was supported by Universiti Kebangsaan Malaysia [Grant Number ST-2019-016]. Notes on contributors Mojtaba Borza is Ph.D. of applied mathematics doing research in optimization and operations research under the supervision of Dr. Azmin Sham Rambely. Dr. Azmin Sham Rambely is associate professor of department of mathematics, University Kebangsaan Malaysia. Her research interest is optimization and bio-mechanic. References [1] Tanino T, Tanaka T, Inuiguchi M. Multi-objective programming and goal programming: theory and applications. Vol. 21. Berlin: Springer Science & Business Media; 2013. 334 M. BORZA AND A. S. RAMBELY [2] Stancu-Minasian IM. Fractional programming: theory, methods and applications. Vol. 409. Dor- drecht: Kluwer Academic Publishers; 1997. [3] Radhakrishnan B, Anukokila P. Fractional goal programming for fuzzy solid transportation prob- lem with interval cost. Fuzzy Inf Eng. 2014;6(3):359–377. [4] Wang Y, Liu L, Guo S, et al. A bi-level multi-objective linear fractional programming for water con- sumption structure optimization based on water shortage risk. J Clean Prod. 2019;237:117829. [5] Ahmad S, Ullah A, Akgül A, et al. Analysis of the fractional tumour-immune-vitamins model with Mittag-Leffler kernel. Results Phys. 2020;19:103559. [6] Das SK, Edalatpanah SA, Mandal T. Application of linear fractional programming problem with fuzzy nature in industry sector. Filomat. 2020;34(15):5073–5084. [7] Chakraborty M, Gupta S. Fuzzy mathematical programming for multi objective linear fractional programming problem. Fuzzy Sets Syst. 2002;125:335–342. [8] Veeramani C, Sumathi M. Fuzzy mathematical programming approach for solving fuzzy linear fractional programming problem. RAIRO – Oper Res. 2014;48(1):109–122. [9] De PK, Deb M. Solution of multi objective linear fractional programming problem by Tay- lor series approach. 2015 international conference on man and machine interfacing (MAMI); Bhubaneswar; 2015.p.1–5. [10] Dinkelbach W. On nonlinear fractional programming. Manage Sci. 1967;13(7):492–498. [11] Güzel N. A proposal to the solution of multiobjective linear fractional programming problem. Abstr Appl Anal. 2013;2013:1–4. [12] Nayak S, Ojha AK. Solution approach to multi-objective linear fractional programming problem using parametric functions. Opsearch. 2019;56(1):174–190. [13] Pal BB, Moitra BN, Maulik U. A goal programming procedure for fuzzy multiobjective linear fractional programming problem. Fuzzy Sets Syst. 2003;139(2):395–405. [14] Toksari MD. Taylor series approach to fuzzy multiobjective linear fractional programming. Inf Sci. 2008;178:1189–1204. [15] Borza M, Rambely AS, Saraj M. Fuzzy approaches to the multi objectives linear fractional pro- gramming problems with interval coefficients. Asian J Math Comput Res. 2015;4:83–94. [16] Nayak S, Ojha AK. Multi-objective linear fractional programming problem with fuzzy parameters. In: Bansal JC, Das KN, Nagar AK, Deep K, Ojha AK, editors. Soft computing for problem solving. Singapore: Springer; 2019. p. 79–90. [17] Borza M, Rambely AS. A linearization to the sum of linear ratios programming problem. Mathe- matics. 2021;9(9):1004. [18] Liu X, Gao YL, Zhang B, et al. A new global optimization algorithm for a class of linear fractional programming. Mathematics. 2019;7(9):867. [19] Charnes A, Cooper WW. Programming with linear fractional functionals. Nav Res Logist Q. 1962;9(3–4):181–186. [20] Antunes CH, Alves MJ, Clímaco J. Multiobjective linear and integer programming. New York (NY): Springer; 2016. [21] Pramy FA, Islam MA. Determining efficient solutions of multi-objective linear fractional program- ming problems and application. Open J Optim. 2017;6(4):164–175. [22] Das SK, Mandal T, Edalatpanah SA. A new approach for solving fully fuzzy linear fractional programming problems using the multi-objective linear programming. RAIRO – Oper Res. 2017;51(1):285–297.
Journal
Fuzzy Information and Engineering
– Taylor & Francis
Published: Jul 3, 2021
Keywords: Efficient solution; membership function; max–min technique; linear programming; fractional programming