Svoboda | Graniru | BBC Russia | Golosameriki | Facebook
Next Article in Journal
STM-Suite, an Online Platform for the Assessment of Memory Functions Discriminates among Subgroups of Children with Different Types of Specific Learning Disorders
Previous Article in Journal
Serum Tumor Necrosis Factor-Related Apoptosis-Inducing Ligand and the Cardiovascular Disease Continuum: Insights from Hypertensive Urgencies and Acute Heart Failure Events
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solving the Fuzzy Transportation Problem by a Novel Particle Swarm Optimization Approach

by
Chrysanthi Aroniadi
and
Grigorios N. Beligiannis
*
Department of Food Science and Technology, University of Patras, Agrinio Campus, G. Seferi 2, 30100 Agrinio, Greece
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(13), 5885; https://doi.org/10.3390/app14135885
Submission received: 5 June 2024 / Revised: 2 July 2024 / Accepted: 3 July 2024 / Published: 5 July 2024

Abstract

:
The fuzzy transportation problem (FTP) represents a significant extension of the classical transportation problem (TP) by introducing uncertainly and imprecision into the parameters involved. Various algorithms have been proposed to solve the FTP, including fuzzy linear programming, metaheuristic algorithms and fuzzy mathematical programming techniques combined with artificial neural networks. This paper presents the application of trigonometric acceleration coefficients-PSO (TrigAC-PSO) to solve the FTP. TrigAC-PSO is a variation of the classical particle swarm optimization algorithm, which has already been applied to solve the TP showing remarkable success. This fact constitutes the main reason that drives the utilization of TrigAC-PSO in current contribution to further investigate its performance in solving the FTP. TrigAC-PSO’s adaptability to handle fuzzy data by solving the FTP via instances with classic fuzzy numbers and generalized fuzzy numbers is explored through a comprehensive comparison between TrigAC-PSO and established methods applied to solve the FTP. The comparative analysis, with recent state-of-the-art algorithms, demonstrates the efficiency and robustness of the proposed method in solving the FTP across various scenarios. Through experimental results and performance metrics, the superiority of the proposed method is presented by achieving optimal solutions. The innovation of current research contributes to advancing the field of fuzzy optimization while providing variable insights into the application of TrigAC-PSO in real-world scenarios.

1. Introduction

In the realm of logistics and supply chain management, the transportation problem (TP) stands as a fundamental challenge, aiming to satisfy the distribution of goods from suppliers to consumers while minimizing overall transportation costs. Over the years, extensive research has been conducted to develop efficient algorithms and methods for solving various transportation problems with precise and deterministic values of supply and demand units such as transportation costs.
The northwest corner method (NWC) is one of the methods that discovers a basic feasible solution to various transportation problems [1]. The method is named after the fact that it stands from the northwest corner cell of the transportation table and proceeds in a systematic manner. Despite the fact that is an easy and applicable method, it does not always guarantee an optimal solution [1]. The least-cost method (LCM) is an alternative method, which focuses on selecting the cell with the minimum cost during the allocation process [2]. Vogel’s approximation method (VAM) is known for yielding values relatively close to the optimum, often aligning precisely with it [3]. Although VAM is more complex than the LCM, since it imposes penalties on cells characterized by high transportation costs, it tends to achieve more favorable initial solutions than LCM [3].
However, in real-world scenarios, it is difficult to predict supply and demand quantities “a priori” as they are inextricably linked by the season and by the trends, prevailing by consumers and the financial circumstances of each corporation. Equivalent situations unfold with the prices of transportation which usually depend on factors such as fluctuating fuel prices, the weather, dynamic market conditions and varying transportation routes. To handle with these uncertainties, researchers have turned to fuzzy logic, a mathematical framework that allows for the representation of vagueness and ambiguity in decision-making processes. Going into further detail, fuzzy set theory provides a robust foundation for modeling and analyzing problems with imprecise information, emerging as a promising avenue to optimize transportation problems.
A fuzzy transportation problem (FTP) is a variation of the traditional transportation problem in which the transportation cost, supply and demand are fuzzy quantities [4]. There is some research that has relied on fuzzy environments, proving that the fuzzy decision-making method is becoming paramount. The inspirers of the idea of decision-making in fuzzy environments were Bellman and Zadeh in 1970 [5]. Following their study, new findings emerged, with the most recent studies being noteworthy of mention. Chanas et al., in 1996, introduced the concept of an optimal solution for the TP with fuzzy coefficients represented as fuzzy numbers and developed an algorithm to attain an optimal solution [6]. Ahmed, Khan and Uddin presented a new algorithm for finding an initial basic feasible solution of the TP when the transportation matrix contains both fuzzy and crisp numbers [7]. Another intriguing respective was offered by Chakraborty and Chakraborty in 2009 [8]. Their method was based on the minimization of transportation cost as well as time of transportation when the demand, supply and shipping cost per unit of the quantities are fuzzy. Basirzadeh, two years later, introduced an approach for solving a wide range of transportation problems by transforming fuzzy quantities into crisp values using their new method [4]. During the same year, Nagoor Gani et al. suggested their individual version for solving the fuzzy transportation problem, utilizing the simplex algorithm [9]. Shanmugasundaram also worked on fuzzy exact algorithms. He successfully developed the fuzzy version of Vogel’s and MODI methods for obtaining a fuzzy initial basic solution and fuzzy optimal solution, respectively [10]. Balasubramanian and Subramanian formulated an alternative method for dealing with a special type of FTP, in which the cost is depicted as triangular fuzzy number. In their research, they utilized ranking methods for numbers to evaluate the fuzzy objective values of the objective function and determine the optimal alternative [11]. Pandian and Natarajan introduced a novel algorithm known as the fuzzy zero point method for finding a fuzzy optimal solution for the FTP, where the transportation cost, supply and demand are represented by trapezoidal fuzzy numbers [12]. Gani and Razak formulated a two-stage-cost-minimizing FTP, where supply and demand are trapezoidal fuzzy numbers, while the costs remain crisp [13]. Malini and Kennedy introduced a method for solving the FTP by utilizing octagonal fuzzy numbers [14]. Ekanayake and Ekanayake based their research on Yager’s robust ranking method [15]. They proposed an alternative algorithm for finding an initial basic solution to the fuzzy triangular and trapezoidal transportation problem [15].
In the literature, there are many other surveys consisted of problems which are solved by various methods via generalized trapezoidal fuzzy numbers. Specifically, generalized fuzzy numbers extend the concept of classical fuzzy numbers by introducing additional parameters to provide a more versatile representation of uncertainty. Kaur and Kuman proposed a new approach for solving the FTP using generalized fuzzy numbers, which had not been applied till that point [16]. Building upon the concept of Kaur and Kuman, Ali Ebrahimnejad strived to diminish the computational complexity of the existing method. The technique applied in his study was simpler, computationally efficient and novel compared to the proposed method by Kaur and Kuman [17]. Chen contributed to the previously mentioned research by proposing the concept of generalized fuzzy numbers when there are no prescribed limits on the membership function to the normal form [18].
FTP approximation is a crucial aspect in the realm of computational intelligence algorithms. Various optimization techniques including metaheuristic algorithms and mathematical programming approaches have been adapted or developed to handle fuzzy transportation problems effectively. Lin solved the transportation problem with fuzzy coefficients using genetic algorithms (GA) [19]. Halder and Jana emphasized on four-dimensional multi-item TP through GA and particle swarm optimization (PSO) [20]. New approaches in metaheuristics to solve the fixed TP in a fuzzy environment were developed by Sadeghi-Moghaddamm et al. [21]. Gupta et al. design a hybrid GA-PSO algorithm to solve the traveling salesman problem (TSP) [22]. Last but not least, Singh and Singh, in 2021, suggested an extension of PSO for solving the FTP [23]. Their research, though, focused on solutions presented in a crisp form from fuzzy sets. All the abovementioned techniques, achieved satisfactory results in optimizing some instances of the fuzzy transportation problem. However, they were applied to a rather small set of instances with more less-similar characteristics. So, it is obvious that a more generic optimization algorithm should be designed and applied to optimize all the above fuzzy transportation instances, achieving results comparable or even better to the ones found so far.
Motivated by the aforementioned applications of metaheuristic algorithms, this study focuses on the effective application of PSO to solve the FTP. Twenty-eight distinct instances are tacked using ten well-established methods found in the existing literature. The spotlight will be on the trigonometric acceleration coefficients-PSO (TrigAC-PSO) algorithm, a variation of the classical PSO algorithm, introduced by Aroniadi and Beligiannis, in 2023 [24]. Extensive experimentation was carried out to assess the stability and the performance of the proposed algorithm in navigating transportation problems with fuzzy costs.
Furthermore, various particle configurations were examined to evaluate TrigAC-PSO’s capacity to achieve optimal results, demonstrating its superior performance.
The contribution of this paper is as follows:
  • According to our knowledge, PSO has already been applied for solving the FTP. However, the conducted experiments lack depth in their data analysis, as they predominantly focus on individual instances of the problem rather than a diverse range of scenarios. Introducing a novel approach, TrigAC-PSO is applied for the first time for solving the FTP, highlighting exceptional performance across a comprehensive set of instances;
  • Nevertheless, the selected problem instances encompass various types of fuzzy numbers such as triangular and trapezoidal fuzzy sets. Moreover, the fuzzy numbers extend beyond conventional representations, encompassing both classic fuzzy numbers and fuzzy generalized numbers, thereby adding depth and complexity to this study. It is the first time that fuzzy generalized numbers are applied to a particle swarm optimization variation;
  • The evaluation of results, from each method, is not based entirely on their attainment of the optimal solution. Instead, this study estimates the method’s performance according to the degree of membership in which their solutions approach optimality as determined by membership functions. Once more, the TrigAC-PSO method demonstrates remarkable completeness against other methods.
The remainder of this paper is organized as follows: Section 2 presents the mathematical formulation of the TP. Definitions for fuzzy logic and fuzzy numbers are briefly described in Section 3. Section 4 presents the FTP. PSO as well as TrigAC-PSO are represented in detail in Section 5. An extensive and intriguing study case is conducted in Section 6, solving the FTP by using the TrigAC-PSO algorithm. Subsequently, a comparative analysis is presented, comparing results obtained by the application of TrigAC-PSO against well-established methods documented in the respective literature. Lastly, conclusive remarks and future recommendations are presented in Section 7.

2. Transportation Problem (TP)

The transportation problem (TP) is a classic optimization problem in the field of operation research and logistics. It deals with the optimal allocation of goods or resources from multiple suppliers to multiple consumers, in an effort to minimize the overall cost while satisfying supply and demand constraints. The TP is formulated as a linear programming problem.
Let us denote the following:
  • m = n u m b e r   o f   s u p p l i e r s ;
  • n = n u m b e r   o f   c o n s u m e r s ;
  • s i = t h e   c a p a c i t y   o f   s u p p l i e r   i ,   i = 1,2 , m ;
  • d j = t h e   c a p a c i t y   o f   c o n s u m e r   j ,   j = 1,2 n ;
  • x i j = n u m b e r   o f   u n i t s   t o   b e   t r a n s p o r t e d   f r o m   s u p p l i e r   i   t o   c o n s u m e r   j ;
  • c i j = c o s t   o f   t r a n s p o r t i n g .
The mathematical model of the TP can be formulated as follows:
m i n Z = i = 1 m j = 1 n c i j · x i j
i = 1 m x i j d j     f o r   j = 1 ,   2 ,   ,   n
j = 1 n x i j s i     f o r   i = 1 ,   2 , ,   m
x i j 0     f o r   i = 1 ,   2 ,   , m   ,   j = 1 ,   2 , , n
Equation (1) outlines the objective function that needs to be minimized. Equation (2) encompasses the supply constraints, ensuring that the available quantity at the source points is greater than or equal to the demanded quantity at the destination points. Similarly, Equation (3) ensures that the sum of the quantities transferred from source s i to destination d j does not exceed the available quantity. Equation (4) introduces a necessary condition, specifying that the units x i j must be positive integers. For the sake of simplicity, we assume a balanced condition model, where supply and demand are equal in this paper.
As already mentioned, various algorithms and methods, such as the northwest corner method, the least-cost method and Vogel’s approximation method have been developed for finding a basic feasible solution. In addition, extensions and variations of the TP have emerged to enhance the model’s relevance across diverse scenarios.

3. Fuzzy Logic Definitions

In this section, fundamental definitions will be provided about fuzzy numbers, fuzzy sets and a variation of fuzzy numbers known as generalized fuzzy numbers. This variation will be explored in greater detail using specific examples in Section 6.

3.1. Fuzzy Sets and Membership Functions

Definition 1.
A fuzzy set A ~ in the universe of discourse X can be defined as a set of ordered pairs.
A ~ = { ( x ,   μ A ~ ( x )   ) | x X }
where μ A ~ x is called membership function of A ~ . The membership function μ A ~ x defines all the information contained in a fuzzy set.
Definition 2.
The fuzzy number A is an extension of a regular number in the sense that it does not refer to one single value but rather to a connected set of possible values, where each possible value has each own degree of membership between 0 and 1. The membership function  μ A ~ ( x ) satisfies the following conditions [5]:
(i)
μ A ~ ( x )  is piecewise continuous;
(ii)
μ A ~ ( x )  is a convex fuzzy subset;
(iii)
μ A ~ x : x 0,1 , where
  • μ A ~ ( x ) = 1   i f   x   i s   t o t a l l y   i n   A ;
  • μ A ~ ( x ) = 0   i f   x   i s   n o t   i n   A ;
  • 0 < μ A ~ x < 1   i f   x   i s   p a r t l y   i n   A .
Definition 3.
A fuzzy number A ~ = ( a , b , c ) , is said to be a triangular fuzzy number if its membership function is given by the following equation and that presented in Figure 1:
μ A ~ ( x ) = 0 ,   i f   x a x a b a ,     i f   a x b c x c b ,   i f   b x c 0   ,   i f   x c
Definition 4.
A fuzzy number A ~ = ( a , b , c , d ) is said to be a trapezoidal fuzzy number if its membership function is given by the following equation and that presented in Figure 2:
μ A ~ x = x a b a ,   a x < b 1 ,   b x c x d c d ,   c < x d 0 ,       o t h e r w i s e

3.2. Generalized Fuzzy Numbers

In our comprehensive investigation of the TP involving a broad spectrum of fuzzy numbers, the exploration of generalized fuzzy numbers comprises an indispensable aspect of our research. Kaur and Kuman [16] had the idea of implying the generalized fuzzy numbers in the TP using exact algorithms. In our study, in some of the presented problems the transportation cost will be illustrated by generalized fuzzy numbers. Our basic aim is to assess their effectiveness following the implementation of TrigAC-PSO algorithm.
But how are generalized fuzzy numbers defined? Generalized fuzzy numbers extend the concept of classical fuzzy numbers by providing more specific details through additional parameters.
In [18], Chen represented a generalized trapezoidal fuzzy number A ~ as A ~ = a , b , c , d ; w , where 0 < w 1 and a ,   b ,   c ,   d are real numbers.
According to Chen, the membership function of a fuzzy set A ~ has the following properties.
  • μ A ~ : R 0 ,   w is continuous;
  • μ A ~ x = 0 for all x ,   a u d ,   ;
  • μ A ~ ( x ) is strictly increasing on a , b and strictly decreasing on [c, d];
  • μ A ~ x = w , for all x b , c where 0 < w 1 .
Definition 5.
A fuzzy number A ~ = a ,   b ,   c ,   d ; w is said to be a generalized trapezoidal fuzzy number if its membership function is provided by
μ A ~ x = w · x a b a ,   a x < b w ,             b x c w · x d c d ,   c < x d 0 ,       o t h e r w i s e

3.3. Arithmetic Operations

In this subsection, the method for defining the two basic mathematical operations, addition and multiplication, is provided. The following table (Table 1) defines these operations for triangular fuzzy numbers and trapezoidal fuzzy numbers, as well as generalized fuzzy numbers. The operator ⊕ denotes addition in fuzzy sets, while the operator ⊗ denotes multiplication [16].
These arithmetic operations offer a systematic approach to combining and manipulating different types of fuzzy numbers, preserving their essential attributes of fuzziness and uncertainty that are fundamental to fuzzy logic. By performing addition and multiplication on fuzzy numbers, we generate new fuzzy numbers that encapsulate the collective or altered characteristics of the original fuzzy sets, as defined by their membership functions. This method empowers us to flexibly handle and analyze uncertain or imprecise data, enhancing our ability to model complex systems and make informed decisions across various domains. The numbers A ~ 1 and A ~ 2 are elements of the set of fuzzy numbers, while λ belongs to the set of real numbers.

4. Fuzzy Transportation Problem

The fuzzy transportation problem (FTP) is an extension of the classical TP which introduces fuzzy sets and fuzzy numbers for representing transportation cost, supply and demand [25]. Solving the FTP involves developing mathematical models and algorithms which can affectively deal with the fuzzy nature of the input data.
In our study, we treat the supply and demand values as crisp numbers, whereas the transportation costs are represented as fuzzy sets. The solution of the objective function z ~ is also depicted as a fuzzy number.
Let c ~ i j be the fuzzy cost of transporting one unit of product from warehouse i to consumer j .
Under changing circumstances, the mathematical model of the FTP could be structured as follows.
m i n Z ~ = i = 1 m j = 1 n c ~ i j x i j
ι = 1 m x i j d j     f o r   j = 1 ,   2 ,   ,   n
j = 1 n x i j s i     f o r   i = 1 ,   2 ,   ,   m
x i j 0     f o r   i = 1 ,   , m   ,   j = 1 , , n

5. Particle Swarm Optimization (PSO)

The particle swarm optimization (PSO) algorithm stands out as a contemporary and innovative heuristic approach, gaining widespread popularity over the years due to its straightforward implementation and ability to yield satisfactory solutions [26]. The PSO algorithm is inspired by the collective behavior of animals from studying their interactions to function, resulting in a strong approach that handles optimization challenges in various applications [27].
In the PSO framework, every potential solution is represented as a particle, and the complete set of particles constitutes the algorithm’s population, known as swarm. The optimization process relies heavily on collaboration and exchange of information, with the particles improving their effectiveness through active cooperation within the swarm. Sharing information within the swarm, particles are empowered to consistently enhance their performance, striving towards optimal efficiency. In finer detail, each particle modifies its movement according to its own experience and its neighboring particle’s experience.
A particle is defined by three essential parameters:
  • Position, which indicates its location in the search space;
  • Velocity, which dictates the direction and extent of particle movement;
  • Its previous best position, which operates as a memory mechanism for the positions that the particle has already “visited”.
Consequently, in an n-dimensional search space, each particle of the swarm is represented by x i j = x i 1 ,   x i 2 ,   ,   x i j , and the equation of its position is as follows:
x i j t + 1 = x i j t + v i j t + 1 ,   i = 1 ,   2 ,   n   and   j = 1 ,   2 , n
where
  • x i j t + 1 is its current position;
  • x i j t is its previous position;
  • v i j t + 1 is the velocity in the current iteration t + 1 .
In sequence, the particle’s velocity is represented as v i j , and its expression is determined by the following equation:
v i j t + 1 = w v i j t + c 1 r 1 p b e s t i j   ( t ) x i j   t + c 2 r 2 g b e s t i j   ( t ) x i j   t i = 1 ,   2 , n     and   j = 1 ,   2 , n
  • where v i j t + 1 represents the velocity of the particle in the current iteration, while v t is the velocity in the previous iteration;
  • w denotes the inertia weight, balancing global and local exploitation by incorporating the memory of the previous particles’ direction to prevent major changes in suggested directions;
  • r 1 and r 2 are variables, randomly generated from uniform distribution, in the range 0 ,   1 ;
  • c 1 and c 2 are defined as acceleration coefficients, impacting the efficiency of the PSO method. More precisely, c 1 signifies the particle’s confidence in itself, while c 2 expresses the particle’s confidence in the swarm;
  • p b e s t i j   ( t ) denotes the best position of the particle up to iteration t , while g b e s t i j   ( t ) is the finest position of the entire swarm up to the same iteration;
  • The term c 1 r 1 p b e s t i j ( t ) x i j   t is known as the cognitive component, acting as a form of memory storing the particle’s best previous positions. The cognitive component reflects the tendency of the particles to return to their best positions;
  • The term c 2 r 2   g b e s t i j ( t ) x i j   t is known as the social component. In this context, particles follow the guidance of the swarm’s best position, incorporating knowledge acquired from the swarm.
The acceleration coefficients c 1 and c 2 , along with random variables r 1 and r 2 significantly influence the evolution of the cognitive and social components, determining the velocity value, which is crucial for the search direction of the particles and the convergence of the algorithm. Furthermore, inertia weight plays a key role in the process of providing balance between exploration and exploitation. Periodically, several variations have been developed and introduced, successfully attaining a dynamic approach to the optimal solution by striking a balance between algorithm’s exploration and exploitation searches.
In our previous research [24], we introduced two new variations. They succeeded in finding the optimal solution under extensive experimental analysis and testing of different combinations of values for variables c 1 , c 2 and w .

Trigonometric Acceleration Coefficient–Particle Swarm Optimization (TrigAC-PSO)

Trigonometric acceleration coefficient–particle swarm optimization (TrigAC-PSO), has been proven to be the variant yielding the most favorable outcomes in solving the transportation problem [24]. In fact, the TrigAC-PSO variation underwent testing across a set of 32 instances, and the outcomes were compared with widely recognized exact methods and some preceding variations of the PSO algorithm. In comparison to all other methods, the TrigAC-PSO variation achieved the best results [24].
In TrigAC-PSO, initially, each particle is guided by the collective knowledge and experience acquired by the swarm, were the value of c 2 significantly exceeds the value of c 1 . Subsequently, leveraging the learning mechanism, each particle formulates its unique strategy and accumulates individual experience, with the value of c 2 decreasing while the value of c 1 increases (refer to Equations (11) and (12)). The adjustment of c 1 and c 2 continues until both parameters are equalized to 2 in the algorithm’s last generation.
The subsequent equations are employed to compute the cognitive and social acceleration parameters [24]:
c 1 = c 1 f 2 + sin 2 c 1 i · t t m a x   · π 2  
c 2 = c 2 i + c o s c 2 f · π · t 2 · t m a x 1 2
  • In the initial iteration, the personal acceleration value c 1 i , is set to 0.5, while the social acceleration value c 2 i , is set to 3.5;
  • In the last iteration of the algorithm, both the personal c 1 f and social c 2 f acceleration values are adjusted to 2;
  • The inertia weight (w) is dynamic and depends on the current iteration t and the maximum number of iterations t m a x and is defined by the following equation:
    w = t m a x t t m a x

6. Case Studies and Experimental Results

In this section, the TrigAC-PSO algorithm will be utilized to solve the fuzzy transportation problem. A total of twenty-eight numerical examples will be resolved with results derived from four other established methods. To ensure a thorough investigation, the numerical examples encompass both triangular and trapezoidal fuzzy numbers. Moreover, the inclusion of generalized fuzzy numbers in our domain adds a distinctive dimension, making the application of our variation particularly intriguing in this context.
Initially, the problems will be evaluated in comparison to methods whose reliability has remained consistent over the years. In this context, the TrigAC-PSO algorithm will be assessed alongside the northwest method, the least-cost method, Vogel’s approximation method and the maximum supply with minimum cost method (MOMC) which was introduced in 2015 by Giancarlo de Franca Aguiar et al. [28]. Their classification technique showed a spectacular computation advantage characterized by a higher processing speed and less use of memory. Both MOMC and the previously mentioned methods have been employed in attempting to solve the fuzzy transportation problem. This poses a challenge for us, as we strive to optimize the values of the given problems using TrigAC-PSO.
At this stage, it would be beneficial to make a reference to the initialization process. In the initial steps of Algorithm 1, two vectors, namely S u p p l y and D e m a n d , are created as input (line 1 and 2). Subsequently, the variables m and n are calculated, representing the values of the parameters S u p p l y and D e m a n d , respectively. Then, a matrix is generated containing random Monte Carlo real numbers (line 7). Line 10 rounds the elements of the candidate solution matrix to the nearest integer, as the commodity quantities must be non-negative integers. The subsequent lines of the algorithm commence a process of realignment and redistribution of matrix L to align with the specified supply and demand amounts. In lines 11 and 12, the sum of all elements in each row of matrix L is stored in vector S u m r o w , and similarly, the sum of all elements in each column is stored in vector S u m c o l . Then, two new vectors, s and d , are created by subtracting S u m r o w from S u p p l y and S u m c o l from D e m a n d , respectively. Following this, for each cell in the final matrix, any discrepancies are identified and adjusted accordingly, ensuring that the excess amounts in vectors s and d are zeroed out.
The output of Algorithm 1 is a matrix comprising the initial solution, known as an initial basic feasible solution (IBFS). Subsequently, the PSO algorithm is employed to find the particle achieving the optimal position and its corresponding optimal transportation cost. The supply and demand values are crisp values, whereas transportation cost is depicted using fuzzy sets. However, the expanded positions of the particles, while adhering to demand and supply constraints, exhibited negative or fractional values. Such values are incompatible with the nature of the solution, which necessitates quantities represented solely by positive integers. Consequently, necessary adjustments were implemented to ensure the integrity of the particle’s positions. These sub-algorithms were presented in 2023 by Aroniadi and Beligiannis [24]. To derive the final solution, fuzzy operations outlined in Table 1 will be applied, culminating in the representation of the final iteration cost as a fuzzy number.
Algorithm 1: Initialization algorithm
Applsci 14 05885 i001
The entire algorithmic approach is implemented utilizing MATLAB R2021b. TrigAC-PSO was tested on a set of different dimensional problems, both triangular and trapezoidal fuzzy sets as classical and generalized fuzzy numbers. Table 2 presents the 28 test instances that are examined in our research [15,17,25,29,30,31,32,33,34,35,36,37,38,39].
In Table 3 and Figure 3, provided below, the performance of both exact methods and the Trig-PSO approach across 20 Monte Carlo runs is showcased. In every cell, the best value of each method is depicted using bold for the cases where the algorithms achieved the optimal solution. The last column provides the optimal solution for each numerical instance. The values highlighted in bold correspond to the optimal ones.
As illustrated, the NW method attains optimal solutions in 3 out of 28 test instances (10.71%). LCM achieved the optimal solution in 5 out of 28 instances (17.86%). VAM surpasses the previous methods, securing optimal solution in 15 out of 28 instances (53.57%). Despite being introduced later than the preceding methods, MOMC fell short of exceeding Vogel’s method, as it succeeded in identifying the optimal solution in six instances (21.43%). Notably, the TrigAC-PSO method emerges as the most efficient, achieving optimal results in 24 out of 28 instances (81.71%).
The deviation measurement quantifies the extent of variation between observed values and expected or target values. It provides a numerical representation of how much individual data values differ from the average or desired outcome and comprises a metric of each algorithm’s stability. Deviation is given by the following formula:
D e v = x i j o p t i m a l o p t i m a l      
where x i j is the current solution.
Based on the data presented in Table 4, it is apparent that VAM, MOMC and TrigAC-PSO appear to be more efficient than the other methods. Therefore, the solutions achieved by NWC differ from the optimal solution by 27.3%, the results of LCM by 10.12% and MOMC by 7.5%. Vogel’s approximation method presented higher levels of efficiency since the average of the export solutions differ from the optimal solutions by 2.7%. The TrigAC-PSO method is now highlighted as the most efficient method, surpassing all other methods, with a percentage close to zero. This indicates its remarkable stability and accuracy in predicting the optimal value with near-perfect precision.
In Boolean (classical) logic, transportation costs are represented as precise, fixed values. In fuzzy logic, transportation costs, as supply and demand quantities, are represented as approximate values, allowing for gradual transitions between costs and shipments.
The pivotal aspect of our study was the visual representation of the final solution using fuzzy numbers. In classical logic, a solution is deemed optimal only if its value aligns precisely with that of the optimal solution. However, the scientific quandary arises: How equitable is it to discard other solutions even when their deviation from the optimal one is infinitesimal?
This quandary sparked our original concept, i.e., to represent the costs incurred by exact methods using fuzzy sets and to ascertain the degree of truth inherent in each algorithm’s solution belonging to this set. Consequently, every solution is now endowed with its unique value, irrespective of its optimality status. The binary classification of classical logic, where 0 signifies non-optimality and 1 denotes optimality, is supplanted by a continuum with the interval 0 ,   1 , wherein each number signifies the solution’s degree of belongingness to the fuzzy set. In this study, the Gaussian membership function will be employed to ascertain the truthfulness of each solution.
The Gaussian membership function is defined by a bell-shaped curve characterized by its mean μ and standard deviation σ . The function assigns a degree of membership to each element in a fuzzy set, based on its proximity to the mean. Mathematically, the Gaussian membership function is given by the formula:
μ x = e x μ 2 2 σ 2  
where
  • x is the input value;
  • μ is the mean of the fuzzy numbers;
  • σ is the standard deviation of the fuzzy numbers.
The following table (Table 5) shows the degree of truth for each solution of the 28 problems.
We notice that among the methods that were unable to precisely identify the optimal solution, they achieved values that were remarkably close to the optimal cost. The degree of truth studies how closely the solution obtained by each method approaches our optimal solution offering a clearer perspective. When the degree of truth is 1, the NWM solutions exhibited optimality with a degree of truth of 0.6061 (equivalent to 60.61% optimality). In contrast, both LCM and MOMC attained optimality with a substantially higher degree, scoring 0.8261 (82.61% optimality) and 0.8546 (85.46% optimality), respectively, surpassing the NWM. While VAM succeeded in discovering the optimal solution in half of the instances, its remaining solutions were notably proximal, as VAM achieved the optimal with degree of truth as high as 0.9638 (96.38%). Yet again, TrigAC-PSO stands out significantly, having achieved the optimal solution in 24 out of 28 problems. However, in instances where TrigAC-PSO did not achieve the desired results, its attained solution diverged slightly from the optimal one, as indicated by the Gaussian membership function. TrigAC-PSO’s precision is evident, reaching the optimal solution with a remarkable truth score of 0.9991 (99.91% optimality).
The selection of these methods was deliberate, guided by two principal factors. Firstly, they were chosen because they have been established for years, verifying their reliability. Secondly, our intention was to demonstrate a compelling narrative: Despite their perceived obsolescence and occasional failure to attain optimality, these methods possess the potential for resurgence and redemption through fuzzy sets, transcending the boundaries of classical logic to reclaim their erstwhile prominence.
Upon culmination of the exhaustive comparative analysis involving the aforementioned methods, it became apparent that TrigAC-PSO appeared as the best choice, distinguishing itself as the quintessential option. In a quest to examine further its rustiness and precision, TrigAC-PSO underwent additional tests against contemporaneous methods introduced concurrently. Firstly, the proposed algorithm is tested on twelve instances, from Table 2, against the algorithm introduced by Ekanayake and Ekanayake, in 2023 [15]. This algorithm was based on Yager’s robust ranking method [16].
But how is a ranking function defined?
An efficient approach for comparing fuzzy numbers is via the use of the racking function R : F ( R ) R where F ( R ) is a set of fuzzy numbers defined on a set of real numbers, which convents a fuzzy number into a crisp value, under specific circumstances [16].
The following table (Table 6) gives the formulae of the ranking functions, which have been applied in this manuscript.
The subsequent table (Table 7) presents the outcomes of the solutions derived from both TrigAC-PSO and the research conducted by Ekanayake and Ekanayake [15].
Once more, TrigAC-PSO algorithm stands out by attaining the optimal solution with a membership degree of 1, demonstrating its superiority. Similarly, the method of Ekanayake and Ekanayake yields impressive outcomes, with an aggregated membership degree reaching 0.9995.
At this juncture, it is intriguing to delve into the realm of generalized fuzzy numbers, where the application of a metaheuristic approach to the FTP marks a pioneering endeavor. Five problems, from Table 2, will be subjected to comparison against TrigAC-PSO, shedding light on the algorithm’s adaptability and efficacy in undertaking challenges modeled by fuzzy generalized numbers [16,17,30,39].
The ensuing table (Table 8) illustrates the outcomes generated by four distinct methods applied to five of the aforementioned problem sets. These solutions are depicted as fuzzy numbers, subsequently subjected to defuzzification through the ranking function (Table 6). This process facilitates an in-depth analysis of the deviation from the optimal values, as shown in Table 9, thus enhancing comprehension for the reader.
Upon reviewing the table above, the following conclusions can be drawn:
  • Ebrahimnejad’s approach as well as Thota and Raja’s approach exhibited zero deviation from the optimal solution, establishing them as the preferred methods for solving the TP involving fuzzy generalized numbers;
  • TrigAC-PSO’s performance in this scenario was exceptional, yet again. This method almost reached an ultimate solution with a deviation rate of 1.71%;
  • Kaur and Kuman’s method demonstrated commendable efficiency with a deviation from the optimum standing at 2.89%;
  • Conversely, the outcomes derived from Mathur’s method exhibited a substantial deviation from the optimal solution, amounting to 20.54%.
Based on these assessments, it is apparent that our method, the TrigAC-PSO algorithm, ensures maximum accuracy, completeness and stability, adeptly handling both the classical TP and the more intricate FTP with remarkable proficiency.
Metaheuristic algorithms consistently emerge as optimal solutions for optimization problems with the PSO algorithm standing out as particularly advantageous. In our innovative variation, TrigAC-PSO, we have methodically designed the cognitive and social component to align with our philosophy, prioritizing the particle’s assimilation of swarm behavior and collective experience to follow its own decision-making process. Additionally, the inertia weight w intricately intertwines with the number of iterations outlined in our algorithm. The remaining consideration is to assess whether the number of particles, comprising our swarm, influences the evolution and quality of our solutions. The solutions for 28 problems were found using 20, 35 and 50 particles, generated by Algorithm 1. These solutions are summarized in the table below (Table 10). The highlighted indications show the cases where the method failed to identify the optimal solution but achieved a value very close to it.
A brief examination reveals minimal deviation among solutions, often rendering them almost identical. However, a more thorough examination was deemed necessary. Each of the 28 problems was subjected to runs with varying numbers of particles and 20 tests conducted for each instance.
Table 11 expands the accuracy rate of each algorithm for 20, 35 and 50 particles. The accuracy rate is commonly defined based on the following relationship:
A c c u r a c y = T O R T R
where T O R is the total number of runs where an optimal solution was found and T R is the number of runs in general. In our scenario, the experiments were carried out for 20 runs.
The data represented in Table 11 indicate that the algorithm achieves an accuracy rate of 66.07% for 20 particles. For 35 particles, the algorithm’s accuracy rate improves to 74.03%. Furthermore, solutions obtained using 50 particles exhibit even higher accuracy, reaching 75.54%. This pattern indicates that as the number of particles increases, there is enhanced interaction and information exchange within the swarm, resulting in higher solution accuracy rates. Figure 4 also reveals the high levels of accuracy.
In the following table (Table 12), the most important statistical measures for 20 independent runs are represented for 20, 35 and 50 particles, respectively.
The experimental findings highlight the remarkable stability of our proposed algorithm in solving the FTP. Across all scenarios, the mean values closely approximate the optimal solution, regardless of whether 20, 35 or 50 particles are utilized. This emphasizes not only the efficiency but also the robustness of the TrigAC-PSO method. Notably, the CV value, a crucial metric for measuring an algorithm’s stability, remains consistently low. Specifically, for 20 particles, the mean CV is 1.82%, reducing to 1.53% for 35 particles and further decreasing to 0.91% for 50 particles. This tendency suggests that increasing the number of particles improves algorithm’s stability. The consistently low CV values affirm the reliability of our algorithm, establishing it as a robust competitor among established methods.

7. Conclusions

In summary, the transportation problem stands as one of the pivotal challenges within the realm of operational research. While numerous methods have successfully solved it through advanced algorithms, the question remains: can these solutions be effectively implemented under non-ideal conditions in real-world scenarios? This question finds its answer in the realm of the fuzzy transportation problem (FTP), where classical logic is supplanted by fuzzy logic.
In this publication, we solved the FTP using TrigAC-PSO algorithm, a variant of particle swarm optimization (PSO), which effectively determined the classical transportation problem. Our method underwent rigorous testing across a spectrum of problem instances and was subsequently compared against well-established techniques documented in the respective literature. Our findings unequivocally demonstrate TrigAC-PSO’s superiority in terms of accuracy, negligible deviation from optimal solutions, and the quality of its solutions. Moreover, beyond traditional fuzzy numbers, TrigAC-PSO was tested under conditions where numbers were classical fuzzy and generalized fuzzy ones, respectively. In this context, TrigAC-PSO showed excellent performance. Based on comprehensive experimentation and meticulous research, it becomes noticeable that TrigAC-PSO emerges as the ideal option for resolving the complexities in the FTP.
As we look to the future, the application of other PSO models, such as ant colony optimization, to solve both the TP and the FTP, will be part of our future work. Moreover, embracing and improving TrigAC-PSO presents hopeful opportunities for improving transportation logistics and management, thus enabling smoother and more productive operations in real-world scenarios, an aspect that presents another challenge for us to explore in our future research.

Author Contributions

Conceptualization, C.A. and G.N.B.; methodology, C.A. and G.N.B.; software, C.A.; validation, C.A. and G.N.B.; formal analysis, C.A. and G.N.B.; investigation, C.A. and G.N.B.; resources, C.A.; data curation, C.A.; writing—original draft preparation, C.A.; writing—review and editing, C.A. and G.N.B.; visualization, C.A.; supervision, G.N.B.; project administration, G.N.B. 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

The majority of the data presented in this research article are available within the manuscript. The data and the programming code used in this paper can be sent, upon request, to the interested reader. Please contact [email protected].

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Deshpande, V.A. An optimal method for obtaining initial basic feasible solution of the transportation problem. In Proceedings of the National Conference on Emerging Trends in Mechanical Engineering Patel College of Engineering & Technology (GCET), Vallabh Vidyanagar (ETME-2009), Vallabh Vidyanagar, India, 13 October 2010; Volume 20, p. 21. [Google Scholar]
  2. Taha, H.A. Operations Research: An Introduction, 8th ed.; Pearson Prentice Hall: Hoboken, NJ, USA, 2007. [Google Scholar]
  3. Amaliah, B.; Fatichah, C.; Suryani, E. A new heuristic method of finding the initial basic feasible solution to solve the transportation problem. J. King Saud Univ.-Comput. Inf. Sci. 2022, 34, 2298–2307. [Google Scholar] [CrossRef]
  4. Basirzadeh, H. An approach for solving fuzzy transportation problem. Appl. Math. Sci. 2011, 5, 1549–1566. [Google Scholar]
  5. Bellman, R.E.; Zadeh, L.A. Decision-making in fuzzy environment. Manag. Sci. 1970, 17, B141–B164. [Google Scholar] [CrossRef]
  6. Chanas, S.; Kuchta, D. A concept of optimal solution of the transportation with Fuzzy cost coefficient. Fuzzy Sets Syst. 1996, 82, 299–305. [Google Scholar] [CrossRef]
  7. Ahmed, N.; Khan, A.; Uddin, M. Solution of Mixed Type Transportation Problem: A Fuzzy Approach. Bul. Institutului Politeh. Din Iaşi 2015, 61, 19–31. [Google Scholar]
  8. Chakraborty, A.; Chakraborty, M. Cost-time minimization in a transportation problem with fuzzy parameters: A case study. J. Transp. Syst. Eng. Inf. Technol. 2010, 10, 53–63. [Google Scholar] [CrossRef]
  9. Gani, A.N. Simplex Type Algorithm for Solving Fuzzy Transportation Problem. Tamsui Oxf. J. Inf. Math. Sci. 2011, 27, 89–98. [Google Scholar]
  10. Shanmugasundari, M.; Ganesan, K. A novel approach for the fuzzy optimal solution of fuzzy transportation problem. Int. J. Eng. Res. Appl. 2013, 3, 1416–1421. [Google Scholar]
  11. Balasubramanian, K.; Subramanian, S. Optimal Solution of Fuzzy Transportation Problems Using Ranking Function. Int. J. Mech. Prod. 2018, 8, 551–558. [Google Scholar] [CrossRef]
  12. Pandian, P.; Natarajan, G. A new algorithm for finding a fuzzy optimal solution for fuzzy transportation problems. Appl. Math. Sci. 2010, 4, 79–90. [Google Scholar]
  13. Gani, A.; Razak, K.A. Two stage fuzzy transportation problem. J. Phys. Sci. 2006, 10, 63–69. [Google Scholar]
  14. Malini, S.; Kennedy, F. An approach for Solving Fuzzy Transportation Problem using Octagonal Fuzzy Numbers. Appl. Math. Sci. 2013, 7, 2661–2673. [Google Scholar] [CrossRef]
  15. Ekanayake, D.; Ekanayake, U. An Average Based Method for Finding the Basic Feasible Solution for the Fuzzy Transportation Problems. Am. J. Appl. Sci. Res. 2023, 9, 1–13. [Google Scholar]
  16. Kaur, A.; Kumar, A. A new approach for solving fuzzy transportation problems using generalized trapezoidal fuzzy numbers. Appl. Soft Comput. 2012, 12, 1201–1213. [Google Scholar] [CrossRef]
  17. Ebrahimnejad, A. A simplified new approach for solving fuzzy transportation problems with generalized trapezoidal fuzzy numbers. Appl. Soft Comput. 2014, 19, 171–176. [Google Scholar] [CrossRef]
  18. Chen, S.H. Operations of fuzzy numbers with function principal. Tamkang J. Manag. Sci. 1985, 6, 13–25. [Google Scholar]
  19. Lin, F.T. Solving the transportation problem with fuzzy coefficients using genetic algorithms. In Proceedings of the IEEE International Conference on Fuzzy Systems, Jeju, Republic of Korea, 20–24 August 2009; pp. 1468–1473. [Google Scholar]
  20. Sharmistha, H.J.; Biswapati, J. Application of fuzzy programming techniques to solve solid transportation problem with additional constraints. Oper. Res. Decis. 2020, 30, 67–84. [Google Scholar]
  21. Sadeghi-Moghaddam, S.; Hajiaghaei-Keshteli, M.; Mahmoodjanloo, M. New approaches in metaheuristics to solve the fixed charge transportation problem in a fuzzy environment. Neural Comput. Appl. 2019, 31, 477–497. [Google Scholar] [CrossRef]
  22. Rajshri, G.; Onkar, C.; Dhawade, N. Optimizing Fuzzy Transportation Problem of Trapezoidal Numbers. Int. J. Fuzzy Math. Syst. 2017, 7, 15–23. [Google Scholar]
  23. Singh, G.; Singh, A. Extension of particle swarm optimization algorithm for solving transportation problem in fuzzy environment. Appl. Soft Comput. 2021, 110, 107619. [Google Scholar] [CrossRef]
  24. Aroniadi, C.; Beligiannis, G.N. Applying Particle Swarm Optimization Variations to Solve the Transportation Problem Effectively. Algorithms 2023, 16, 372. [Google Scholar] [CrossRef]
  25. Kumar, P.S. PSK Method for Solving Intuitionistic Fuzzy Solid Transportation Problems. Int. J. Fuzzy Syst. Appl. 2018, 7, 62–99. [Google Scholar] [CrossRef]
  26. Rosendo, M.; Pozo, A. A hybrid particle swarm optimization algorithm for combinatorial optimization problems. In Proceedings of the IEEE Congress on Evolutionary Computation (CEC), Barcelona, Spain, 18–23 July 2010; pp. 1–8. [Google Scholar]
  27. Wang, D.; Tan, D.; Liu, L. Particle swarm optimization algorithm: An overview. Soft Comput. 2017, 22, 387–408. [Google Scholar] [CrossRef]
  28. de Franca Aguiar, G. Barbara de Cassia Xavier Cassins Aguiar, Volmir Eugenio Wilhelm. The MOMC method: A new methodology to find initial solution. Appl. Math. Sci. 2015, 9, 901–914. [Google Scholar]
  29. Farikhin, M.S.; Bayu, S.; Bambang, I. Solving of Fuzzy Transportation Problem Using Fuzzy Analytical Hierarchy Process. In Proceedings of the 2nd International Seminar on Science and Technology, (ISSTEC 2019), South Tangerang, Indonesia, 25 November 2019; Atlantis Press: Dordrecht, The Netherlands, 2020; pp. 10–15. [Google Scholar]
  30. Mathur, N.; Srivastava, P. An Inventive Approach to Optimize Fuzzy Transportation Problem. Int. J. Math. Eng. Manag. Sci. 2020, 5, 985–994. [Google Scholar] [CrossRef]
  31. Bisht, D.; Srivastava, P. One Point Conventional Model to Optimize Trapezoidal Fuzzy Transportation Problem. Int. J. Math. Eng. Manag. Sci. 2019, 4, 1251–1263. [Google Scholar] [CrossRef]
  32. Muhammad, S.; Farikhin; Hariyanto, S.; Surarso, B. Optimal solution of full fuzzy transportation problems using total integral ranking. J. Phys. Conf. Ser. 2018, 983, 012075. [Google Scholar]
  33. Pandian, P.; Natarajan, G. An Appropriate Method for Real Life Fuzzy Transportation Problems. Int. J. Inf. Sci. Appl. 2011, 3, 127–134. [Google Scholar]
  34. Mathur, N.; Srivastava, P.K.; Paul, A. Trapezoidal fuzzy model to optimize transportation problem. Int. J. Model. Simul. Sci. Comput. 2016, 7, 1650028. [Google Scholar] [CrossRef]
  35. Singh, R.; Saxena, V. A new ranking based fuzzy approach for fuzzy transportation problem. Comput. Model. New Technol. 2017, 21, 16–21. [Google Scholar]
  36. Srinivasan, R.; Karthikeyan, N.; Renganathan, K.; Vijayan, D.V. Method for solving fully fuzzy transportation problem to transform the materials. Mater. Today Proc. 2020, 37, 431–433. [Google Scholar] [CrossRef]
  37. Hussain, R.; Jayaraman, P. Fuzzy Transportation Problem Using Improved Fuzzy Russells Method. Int. J. Math. Trends Technol. 2014, 5, 50–59. [Google Scholar] [CrossRef]
  38. Kumar, P.S. A Simple Method for Solving Type-2 and Type-4 Fuzzy Transportation Problems. Int. J. Fuzzy Log. Intell. Syst. 2016, 16, 225–237. [Google Scholar] [CrossRef]
  39. Thota, S.; Raja, P. New Method for Finding an Optimal Solution of Generalized Fuzzy Transportation Problems. Asian J. Math. Sci. 2020, 4, 19–24. [Google Scholar]
Figure 1. Triangular fuzzy number.
Figure 1. Triangular fuzzy number.
Applsci 14 05885 g001
Figure 2. Trapezoidal fuzzy number.
Figure 2. Trapezoidal fuzzy number.
Applsci 14 05885 g002
Figure 3. Number of optimal solutions achieved by each method.
Figure 3. Number of optimal solutions achieved by each method.
Applsci 14 05885 g003
Figure 4. The accuracy of TrigAC-PSO for 20, 35 and 50 particles.
Figure 4. The accuracy of TrigAC-PSO for 20, 35 and 50 particles.
Applsci 14 05885 g004
Table 1. Arithmetic operations.
Table 1. Arithmetic operations.
Triangular Fuzzy NumbersTrapezoidal Fuzzy NumbersGeneralized Trapezoidal Fuzzy Numbers
Let A ~ 1 = ( a 1 , b 1 , c 1 ) and A ~ 2 = ( a 2 , b 2 , c 2 ) be two triangular fuzzy numbers; then,
  • A ~ 1 A ~ 2 = ( a 1 + a 2 , b 1 + b 2 , c 1 + c 2 )
  • λ A ~ 1 = λ a 1 , λ b 1 , λ c 1
Let   A ~ 1 = ( a 1   , b 1 , c 1 , d 1 ) and A ~ 2 = ( a 2 , b 2 , c 2 , d 2 ) be two trapezoidal fuzzy numbers; then,
  • A ~ 1 A ~ 2 = ( a 1   + a 2 ,   b 1   + b 2 , c 1   + c 2 ,   d 1   + d 2 )
  • λ A ~ 1 = ( λ a 1   , λ   b 1   , λ c 1   ,   λ d 1   )
Let A ~ 1 = ( a 1   , b 1 , c 1 , d 1 ;   w 1 ) and A ~ 2 = ( a 2 , b 2 , c 2 , d 2 ;   w 2 ) be two generalized trapezoidal fuzzy numbers; then,
  • A ~ 1 A ~ 2 = ( a 1   + a 2     ,   b 1   + b 2     , c 1   + c 2     ,   d 1   + d 2     ; m i n ( w 1 ,   w 2 ) )
  • λ A ~ 1 = ( λ a 1   , λ   b 1   , λ c 1   ,   λ d 1   ; w 1 )
Table 2. Details of the 28 test instances of the FTP.
Table 2. Details of the 28 test instances of the FTP.
NoFrom JournalNameProblem SizeTypeOptimal Solution
1Ebrahimnejad (2014)Pr. 013 × 3Generalized Trapezoidal Fuzzy Number64.35
2Kumar and Subramanian (2018)Pr. 024 × 4Classic Triangular Fuzzy Number853.35
3Farikhin et al. (2020)Pr. 033 × 4Classic Triangular Fuzzy Number817.17
4Mathur and Srivastava (2020)Pr. 043 × 4Generalized Trapezoidal Fuzzy Number196.00
5Srivastava and Bisht (2018)Pr. 053 × 3Classic Triangular Fuzzy Number166.00
6Srivastava and Bisht (2018)Pr. 063 × 4Classic Triangular Fuzzy Number101.00
7Sam’an et al. (2018)Pr. 073 × 3Classic Trapezoidal Fuzzy Number1770.00
8Pandian and Natarajan (2011)Pr. 083 × 4Classic Trapezoidal Fuzzy Number199.00
9Mathur et al. (2016)Pr. 093 × 3Classic Trapezoidal Fuzzy Number155.25
10Singh and Saxena (2017)Pr. 103 × 4Classic Trapezoidal Fuzzy Number126.00
11Ekanayake and Ekanayake (2023)Pr. 114 × 4Classic Triangular Fuzzy Number294.00
12Ekanayake and Ekanayake (2023)Pr. 123 × 4Classic Triangular Fuzzy Number65.80
13Ekanayake and Ekanayake (2023)Pr. 132 × 3Classic Triangular Fuzzy Number7735.50
14Ekanayake and Ekanayake (2023)Pr. 144 × 4Classic Triangular Fuzzy Number130.68
15Srinivasan et al. (2020)Pr. 156 × 6Classic Triangular Fuzzy Number2170.00
16Ekanayake and Ekanayake (2023)Pr. 163 × 3Classic Trapezoidal Fuzzy Number951.25
17Ekanayake and Ekanayake (2023)Pr. 174 × 3Classic Trapezoidal Fuzzy Number821.25
18Ekanayake and Ekanayake (2023)Pr. 183 × 4Classic Triangular Fuzzy Number149.00
19Ekanayake and Ekanayake (2023)Pr. 193 × 4Classic Trapezoidal Fuzzy Number67.25
20Hussain and Jayaraman (2014)Pr. 203 × 3Classic Triangular Fuzzy Number3640.56
21Hussain and Jayaraman (2014)Pr. 214 × 4Classic Trapezoidal Fuzzy Number3844.00
22Ekanayake and Ekanayake (2023)Pr. 223 × 3Classic Triangular Fuzzy Number295.90
23Ekanayake and Ekanayake (2023)Pr. 233 × 3Classic Triangular Fuzzy Number551.03
24Ebrahimnejad (2014)Pr. 244 × 6Generalized Trapezoidal Fuzzy Number4300.20
25Kumar (2016)Pr. 253 × 4Classic Trapezoidal Fuzzy Number68.00
26Kumar (2016)Pr. 263 × 4Classic Trapezoidal Fuzzy Number141.00
27Thota and Raja (2020)Pr. 273 × 3Generalized Trapezoidal Fuzzy Number91.45
28Thota and Raja (2020)Pr. 283 × 4Generalized Trapezoidal Fuzzy Number75.60
Table 3. The optimal solution of each method for the 28 test instances.
Table 3. The optimal solution of each method for the 28 test instances.
PrNWCLCMVAMMOMCTrigAC-PSOOptimal
0164.3567.6067.6067.6065.1064.35
021046.67870.05853.35855.00853.35853.35
03861.53894.66817.171000.67817.17817.17
04233.00266.50268.00266.50196.00196.00
05166.00190.00172.00172.00166.00166.00
06125.00120.5.0101.00105.00101.00101.00
072025.001790.001770.001800.001770.001770.00
08233.00223.00203.00199.00199.00199.00
09155.25178.25159.25164.50155.25155.25
10144.25140.00130.00126.00126.00126.00
11376.00294.00294.00375.00294.00294.00
12110.6765.8065.8066.0065.8065.80
137736.677735.507735.507736.677735.507735.50
14196.00130.68130.68130.68130.68130.68
154285.002455.002310.002620.002261.002170.00
161035.00971.25951.25951.25951.25951.25
17967.50887.50821.25826.25821.25821.25
18176.00152.00149.00150.00149.00149.00
1993.0067.2567.2577.0067.2567.25
205070.333740.583644.583944.343640.563640.56
214172.004172.004091.003932.003844.003844.00
22486.67339.92295.90340.00295.90295.90
23592.67557.71557.71581.00551.03551.03
246549.907567.804414.954452.904386.454300.20
2593.0073.0070.0068.0068.0068.00
26169.00148.00141.00141.00141.00141.00
27108.8097.5097.4597.5091.4591.45
28134.1895.0075.6095.0082.5075.60
Table 4. The deviation (dev) of the methods for the 28 test instances.
Table 4. The deviation (dev) of the methods for the 28 test instances.
NWMLCMVAMMOMCTrigAC-PSO
Pr. 010.0000000000.0505050500.0505050500.0505050500.011655011
Pr. 020.2265424500.0195699300.0000000000.0019335560.000000000
Pr. 030.0542849100.0948272690.0000000000.2245554780.000000000
Pr. 040.1887755100.3596938770.3673469400.3596938770.000000000
Pr. 050.0000000000.0144578310.0361445800.0361445780.000000000
Pr. 060.2376237620.1930693060.0000000000.0396039600.000000000
Pr. 070.1440677960.0112994350.0000000000.0169491520.000000000
Pr. 080.1708542710.1206030150.0201005000.0000000000.000000000
Pr. 090.0000000000.1481481480.0257649000.0595813200.000000000
Pr. 100.1448412690.1111111110.0317460300.0000000000.000000000
Pr. 110.2789115640.0000000000.0000000000.2755102040.000000000
Pr. 120.6819148930.0000000000.0000000000.0030395140.000000000
Pr. 130.0001512510.0000000000.0000000000.0001512510.000000000
Pr. 140.4998469540.0000000000.0000000000.0000000000.000000000
Pr. 150.9746543770.1313364050.0606060600.2073732710.041935483
Pr. 160.0880420490.0210194420.0000000000.0000000000.000000000
Pr. 170.1780821920.0806697110.0000000000.0060882800.000000000
Pr. 180.1812080540.0201342280.0000000000.0067114090.000000000
Pr. 190.3828996280.0000000000.0000000000.1449814130.000000000
Pr. 200.3927335350.0275232380.0011042300.0834432070.000000000
Pr. 210.0853277840.0853277840.0642559800.0228928200.000000000
Pr. 220.6447110510.1487664750.0000000000.1490368370.000000000
Pr. 230.0755675730.0121227520.0121227500.0543890530.000000000
Pr. 240.5018138690.7385237900.0053369600.0355099760.020057207
Pr. 250.3676470590.0735294120.0294117600.0000000000.000000000
Pr. 260.1773049650.0496453900.0000000000.0000000000.000000000
Pr. 270.1897211590.0661563700.0656096200.0661563700.000000000
Pr. 280.7748015870.2566137570.0000000000.2566137570.091269841
Average0.2729403400.1012376330.0275019800.0750308690.005889912
Table 5. The Gaussian membership function for the 28 test instances.
Table 5. The Gaussian membership function for the 28 test instances.
NoNWMLCMVAMMOMCTrigAC-PSO
Pr. 011.00000.97670.97670.97670.9950
Pr. 020.71150.99851.00000.99991.0000
Pr. 030.60220.34881.00000.03041.0000
Pr. 040.78440.46990.45530.46991.0000
Pr. 051.00000.91370.99440.99441.0000
Pr. 060.61400.71831.00000.98621.0000
Pr. 070.75610.99831.00000.99611.0000
Pr. 080.83190.83190.99491.00001.0000
Pr. 091.00000.99810.96970.98681.0000
Pr. 100.74830.84450.98621.00001.0000
Pr. 110.60451.00001.00000.61391.0000
Pr. 120.00001.00001.00001.00001.0000
Pr. 130.99801.00001.00000.99801.0000
Pr. 140.22711.00001.00001.00001.0000
Pr. 150.00000.13110.61470.00600.9870
Pr. 160.99470.99511.00001.00001.0000
Pr. 170.99740.97091.00000.99951.0000
Pr. 180.48680.99111.00000.99991.0000
Pr. 190.75591.00001.00000.96921.0000
Pr. 200.47620.99240.99990.96651.0000
Pr. 210.99990.99990.99990.99991.0000
Pr. 220.97180.99211.00000.99201.0000
Pr. 230.90930.99670.99670.95141.0000
Pr. 240.00000.00000.99950.99300.9999
Pr. 250.22350.98931.00001.00001.0000
Pr. 260.25900.97501.00001.00001.0000
Pr. 270.01900.99980.99980.99971.0000
Pr. 280.00000.00001.00000.00000.9930
Average0.60610.82620.96390.85460.9991
Table 6. The formulae of ranking functions.
Table 6. The formulae of ranking functions.
Ranking Function for Two Classical Fuzzy
Numbers [7]
Ranking Function for Generalized Trapezoidal Fuzzy Numbers [17]
  • Let A ~ = a , b , c a triangular fuzzy number;
    then,
       R A ~ = a + b + c 3
  • Let A = a , b , c , d a trapezoidal fuzzy number;
         then,
       R A ~ = a + b + c + d 4
  • Let
        A ~ 1 = a 1 , b 1 , c 1 , d 1 ; w 1
    and
        A ~ 2 = a 2 , b 2 , c 2 , d 2 ; w 2
    be two generalized trapezoidal fuzzy numbers and w = min w 1 , w 2 .
    Then,
    R A ~ 1 = w · a 1 + b 1 + c 1 + d 1 4
    and R A ~ 2 = w · a 2 + b 2 + c 2 + d 2 4
Table 7. Comparison between Ekanayake and Ekanayake and the TrigAC-PSO method.
Table 7. Comparison between Ekanayake and Ekanayake and the TrigAC-PSO method.
Pr.Ekanayake Optimal SolutionTrigAC-PSO Optimal SolutionEkanayake Membership ValueTrigAC-PSO Membership Value
Pr. 02(400, 845, 1315)(400, 845, 1315)1.00001.0000
Pr. 11(144, 285, 453)(144, 285, 453)1.00001.0000
Pr. 12(20, 89, 89)(20, 89, 89)1.00001.0000
Pr. 13(5960, 7620, 9630)(5960, 7620, 9630)1.00001.0000
Pr. 14(64, 124, 206)(64, 124, 206)1.00001.0000
Pr. 16(370, 735, 1145, 1595)(370, 715, 1085, 1635)0.99511.0000
Pr. 17(400, 640, 875, 1370)(400, 640, 875, 1370)1.00001.0000
Pr. 18(105, 150, 195)(104, 149, 194)0.99991.0000
Pr. 22(148, 322, 418)(148, 322, 418)1.00001.0000
Pr. 23(347, 566, 760)(347, 554, 752)0.99501.0000
Pr. 25(12, 55, 88, 117)(12, 55, 88, 117)1.00001.0000
Pr. 26(52, 106, 176, 230)(52, 106, 176, 230)1.00001.0000
Table 8. Comparison TrigAC-PSO among four established methods.
Table 8. Comparison TrigAC-PSO among four established methods.
Pr. Ebrahimnejad [17]Thota and Raja [39]Kaur and Kuman [16]Mathur et al. [30]TrigAC-PSO
Pr. 01(117, 205, 352, 613; 0.2)(117, 205, 35, 613; 0.2)(147, 220, 382, 603; 0.2)(197, 240, 382, 643; 0.2)(147, 220, 382, 553; 0.2)
Pr. 04(315, 810, 1220, 1575; 0.2)(318, 813, 1220, 1582; 0.2)(315, 810, 1220, 1575; 0.2)(415, 970, 1460, 1865; 0.2)(315, 810, 1220, 1575; 0.2)
Pr. 24(5148, 6475, 7802, 9244; 0.6)(5148, 6475, 7802, 9244; 0.6)(5307, 6794, 7.872, 9460; 0.6)(5307, 6794, 7.872, 9460; 0.6)(3170.4, 4052.4, 4717.2, 5628; 0.6)
Pr. 27(376, 436, 474, 543; 0.2)(376, 436, 474, 543; 0.2)(413, 459, 506, 572; 0.2)(411, 455, 509, 563; 0.2)(404, 460, 502, 559; 0.2)
Pr. 28(294, 348, 408, 462; 0.2)(294, 348, 408, 462; 0.2)(294, 348, 408, 462; 0.2)(294, 348, 408, 462; 0.2)(294, 348, 408, 462; 0.2)
Table 9. The deviation (dev) of the methods for 5 numerical examples.
Table 9. The deviation (dev) of the methods for 5 numerical examples.
Pr.Ebrahimnejad [17]Thota and Raja [39]Kaur and Kuman [16]Mathur et al. [30]TrigAC-PSO
Pr. 010.00000.00000.05050.13600.0117
Pr. 040.00000.00280.00000.20150.0000
Pr. 240.00000.00000.02780.02780.0213
Pr. 270.00000.00000.06620.66170.0524
Pr. 280.00000.00000.00000.00000.0000
Average0.00000.00060.02890.20540.0171
Table 10. Solutions for 20, 35 and 50 particles.
Table 10. Solutions for 20, 35 and 50 particles.
Pr.20 Particles35 Particles50 ParticlesOptimal
0165.1065.10065.1064.35
02853.35853.35853.35853.35
03817.17817.17817.17817.17
04200.50196.00196.00196.00
05166.00166.00166.00166.00
06101.00101.00101.00101.00
071770.001770.001770.001770.00
08199.00199.00199.00199.00
09155.25155.25155.25155.25
10126.00126.00126.00126.00
11294.00294.00294.00294.00
1265.8065.8065.8065.80
137735.507735.507735.507735.50
14130.68130.68130.68130.68
152327.002330.002261.002170.00
16951.25951.25951.25951.25
17821.25821.25821.25821.25
18149.00149.00149.00149.00
1967.2567.2567.2567.25
203640.563640.563640.563640.56
213844.003844.003844.003844.00
22295.90295.90295.90295.90
23551.03551.03551.03551.03
244389.604386.454399.354392.00
2568.0068.0068.0068.00
26141.00141.00141.00141.00
2796.2596.2596.2596.25
2882.5082.5082.5075.60
Table 11. The accuracy of the 28 problems for 20, 35 and 50 particles.
Table 11. The accuracy of the 28 problems for 20, 35 and 50 particles.
Pr.20 Particles35 Particles50 Particles
010.951.001.00
021.001.001.00
030.951.000.95
040.050.000.05
051.001.001.00
060.150.750.90
070.400.600.70
081.000.750.90
091.001.001.00
100.700.750.90
110.800.950.95
120.450.600.80
131.001.001.00
140.850.900.95
150.000.000.00
161.001.001.00
170.350.730.55
180.901.000.95
190.750.951.00
200.650.650.70
210.200.200.20
221.001.000.95
231.001.001.00
240.000.000.00
250.550.850.90
260.751.000.80
270.050.050.00
281.001.001.00
Average0.660.740.76
Table 12. Statistical measures for 20, 35 and 50 particles.
Table 12. Statistical measures for 20, 35 and 50 particles.
20 Particles35 Particles50 Particles
Pr. 01Mean65.22565.10065.100
St. Dev0.5590.0000.000
Min65.10065.10065.100
Max67.60065.10065.100
CV%0.8570.0000.000
Pr. 02Mean853.350853.350853.350
St. Dev0.0000.0000.000
Min853.350853.350853.350
Max853.350853.350853.350
CV%0.0000.0000.000
Pr. 03Mean817.659817.170817.490
St. Dev2.1870.0001.431
Min817.170817.170817.170
Max826.950817.170823.570
CV%0.2670.0000.175
Pr. 04Mean204.285203.225200.765
St. Dev3.3063.2471.579
Min200.500196.000196.000
Max210.200210.100203.400
CV%1.6181.5980.786
Pr. 05Mean166.000166.000166.000
St. Dev0.0000.0000.000
Min166.000166.000166.000
Max166.000166.000166.000
CV%0.0000.0000.000
Pr. 06Mean102.575101.450101.200
St. Dev0.7120.8090.616
Min101.000101.000101.000
Max103.000103.000103.000
CV%0.6940.7980.608
Pr. 07Mean1818.4131809.38751782.35
St. Dev76.91669.0840675426.25888321
Min1770.00017701770
Max2020.00020201870
CV%4.2303.8180913451.473273106
Pr. 08Mean199.000200.500199.600
St. Dev0.0002.6661.847
Min199.000199.000199.000
Max199.000205.000205.000
CV%0.0001.3290.925
Pr. 09Mean155.250155.250155.250
St. Dev0.0000.0000.000
Min155.250155.250155.250
Max155.250155.250155.250
CV%0.0000.0000.000
Pr. 10Mean126.300126.250126.100
St. Dev0.4700.4440.308
Min126.000126.000126.000
Max127.000127.000127.000
CV%0.3720.3520.244
Pr. 11Mean294.550295.350295.400
St. Dev1.2346.0376.261
Min294.000294.000294.000
Max298.000321.000322.000
CV%0.4192.0442.119
Pr. 12Mean69.53567.64066.720
St. Dev6.3092.3121.888
Min65.80065.80065.800
Max94.50070.40070.400
CV%9.0733.4182.829
Pr. 13Mean7735.5007735.5007735.500
St. Dev0.0000.0000.000
Min7735.5007735.5007735.500
Max7735.5007735.5007735.500
CV%0.0000.0000.000
Pr. 14Mean132.418132.317131.483
St. Dev4.4835.0663.589
Min130.680130.680130.680
Max147.390148.710146.730
CV%3.3863.8292.730
Pr. 15Mean2703.4502738.7002445.350
St. Dev251.641400.107138.771
Min2327.0002330.0002261.000
Max3275.0003585.0002795.000
CV%9.30814.6105.675
Pr. 16Mean951.250951.250951.250
St. Dev0.0000.0000.000
Min951.250951.250951.250
Max951.250951.250951.250
CV%0.0000.0000.000
Pr. 17Mean829.900826.775823.813
St. Dev15.60111.1203.498
Min821.250821.250821.250
Max878.000868.000834.750
CV%1.8801.3450.425
Pr. 18Mean149.450149.000149.150
St. Dev1.4680.0000.671
Min149.000149.000149.000
Max155.000149.000152.000
CV%0.9820.0000.450
Pr. 19Mean68.10067.40067.250
St. Dev2.6850.6710.000
Min67.25067.25067.250
Max79.00070.25067.250
CV%3.9420.9950.000
Pr. 20Mean3641.0963641.2973641.029
St. Dev0.8021.2660.787
Min3640.563640.5603640.560
Max3643.243644.5803643.240
CV%0.0220.0350.022
Pr. 21Mean3896.4133896.4133896.413
St. Dev53.85253.85253.852
Min3844.0003844.0003844.000
Max4020.0004020.0004020.000
CV%1.3821.3821.382
Pr. 22Mean295.900298.667295.900
St. Dev0.00012.3720.000
Min295.900295.900295.900
Max295.900351.230295.900
CV%0.0004.1420.000
Pr. 23Mean551.030551.030551.030
St. Dev0.0000.0000.000
Min551.030551.030551.030
Max551.030551.030551.030
CV%0.0000.0000.000
Pr. 24Mean4518.1504465.5154446.878
St. Dev106.24959.47743.525
Min4389.6004386.4504399.350
Max4857.3004629.0004580.550
CV%2.3521.3320.979
Pr. 25Mean70.40068.15068.200
Var4.1600.3660.696
Min68.00068.00068.000
Max82.00069.00071.000
CV%5.9090.5381.020
Pr. 26Mean143.500141.000143.200
St. Dev4.4430.0004.584
Min141.000141.000141.000
Max151.000141.000155.000
CV%3.0960.0003.201
Pr. 27Mean96.19096.13096.550
St. Dev1.1981.3430.533
Min91.45091.45096.250
Max97.45097.45097.450
CV%1.2461.3970.552
Pr. 28Mean82.50082.50082.500
St. Dev0.0000.0000.000
Min82.50082.50082.500
Max82.50082.50082.500
CV%0.0000.0000.000
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Aroniadi, C.; Beligiannis, G.N. Solving the Fuzzy Transportation Problem by a Novel Particle Swarm Optimization Approach. Appl. Sci. 2024, 14, 5885. https://doi.org/10.3390/app14135885

AMA Style

Aroniadi C, Beligiannis GN. Solving the Fuzzy Transportation Problem by a Novel Particle Swarm Optimization Approach. Applied Sciences. 2024; 14(13):5885. https://doi.org/10.3390/app14135885

Chicago/Turabian Style

Aroniadi, Chrysanthi, and Grigorios N. Beligiannis. 2024. "Solving the Fuzzy Transportation Problem by a Novel Particle Swarm Optimization Approach" Applied Sciences 14, no. 13: 5885. https://doi.org/10.3390/app14135885

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