Public Transportation Hub Location with Stochastic Demand : An Improved Approach Based on Multiple Attribute Group Decision-Making

Urban public transportation hubs are the key nodes of the public transportation system.The location of such hubs is a combinatorial problem. Many factors can affect the decision-making of location, including both quantitative and qualitative factors; however, most current research focuses solely on either the quantitative or the qualitative factors. Little has been done to combine these two approaches. To fulfill this gap in the research, this paper proposes a novel approach to the public transportation hub location problem, which takes both quantitative and qualitative factors into account. In this paper, an improved multiple attribute group decision-making (MAGDM) method based on TOPSIS (Technique for Order Preference by Similarity to Ideal Solution) and deviation is proposed to convert the qualitative factors of each hub into quantitative evaluation values. A location model with stochastic passenger flows is then established based on the above evaluation values. Finally, stochastic programming theory is applied to solve the model and to determine the location result. A numerical study shows that this approach is applicable and effective.


Introduction
With the development of economy and the expansion of city, the role of the public transportation system has become increasingly important.The location and planning of public transportation hubs is an important step to improving transportation efficiency, optimizing resources allocation, and reducing energy consumption [1,2].
The hub location problems (HLPs) involve locating consolidating facilities when flow (travelers, airline passengers, telecom messages, cargo, farm products, mail, etc.) must be sent from origin to destination nodes (e.g., cities).Trips between origin-destination (OD) nodes may represent direct links between OD pairs or indirect paths passing through a hub [1].Researchers have used various methods in an attempt to solve the hub location problem [3][4][5].These methods can be categorized into different subject areas like single-HLP, p-HLP, multiple allocation p-HLP, p-HLP with fixed link cost, minimum-value flow on links model, p-HLP with limited capacity, continuous p-HLP, multiobjective p-HLP, p-Hub center location problem, p-Hub covering location problem, and hub set covering location problem [6,7].
However, current studies have some limitations; the existing researches on transportation hub location usually focus only on either the quantitative (e.g., cost, speed) or qualitative aspects (e.g., commercial potential).With regard to the quantitative factors, scholars usually apply p-HLP, gravity method, nonlinear programming and stochastic programming, and so forth to obtain optimal solutions [1,2,5,8], However, such models cannot incorporate all factors, especially those factors that cannot be measured precisely.While MAGDM, fuzzy evaluation methods [9][10][11][12][13] can make the qualitative factors quantified, the comparability of the quantitative elements is weakened in such methods.Neglecting the influence of quantitative or qualitative factors will eventually affect the rationality of the hub location result.
We suggest that the public transportation hub location problem should consider both quantitative and qualitative factors simultaneously given that they directly influence the optimal location result.Therefore, in order to solve the above problem, this paper proposes a novel line of thought that considers both quantitative and qualitative factors in the public transportation hub location problem.More specifically, we combine an improved multiple attribute group decisionmaking (MAGDM) method, which converts the qualitative factors into quantitative values, with the stochastic location model.
Our study contributes to extant research in three ways.First, unlike previous research that focuses solely on either quantitative or qualitative factors, our model takes both factors into account.Quantitative and qualitative factors should be considered simultaneously in the model because both have influential impacts on the location of public transportation hubs.Second, we propose an improved MAGDM that can objectively weight attributes and DMs (decisionmakers) simultaneously.That is, the decision-making process is conducted under the condition that the attributes and DMs weights are unknown in advance.Weighting attributes and DMs objectively can reduce the subjective elements in the decision-making process.Finally, unlike most of the existing research on public transportation hub location that assumes that passenger flow is known, we treat passenger flow as a stochastic variable in the model.Our method is more practical.Because passenger flow is influenced by many factors that are difficult to estimate and predict, the stochastic programming in our approach should lead to more reasonable implications, compared with research that assumes passenger flow is known.
The rest of the paper is structured as follows.Section 2 reviews related literature.Section 3 introduces the stochastic transportation location model.In Section 4, an improved MAGDM method with unknown attributes and DM weights is presented and discussed.Section 5 provides the general steps of the proposed method; and Section 6 demonstrates a numerical example.Finally, Section 7 presents our conclusions.

Literature Review
The goal of transportation hub location is to choose hubs in a set of alternatives to minimize costs or maximize benefits.In essence, the transportation hub location problem is one of several types of hub location problems.Many researchers have in fact studied hub location problems [1,3,4,14,15].Most of them address the p-hub median problem and its variants.As an extension of the hub location problem, Campbell et al. [8,16] propose a new hub location approach called the hub arc location model, and they analyze the characteristics of the optimal solution with special attention to spatial patterns and relationships.The authors provide integer programming formulations and optimal algorithms for these hub arc problems.With this method, one can not only locate discrete hub facilities but also reduce the unit flow cost.
Robinson and Bookbinder [17] propose a mixed integer programming model to solve the location problem of finishing plants and distribution centers between Canada, Mexico, and the United States; then, through a real-world example, overall supply chain costs (total system costs of inventory, transportation, and facilities) can be minimized under those circumstances of North American Free Trade Agreement (NAFTA).Alumur and Kara [18] review over 100 papers regarding the network hub location problem, including the p-hub median problem, the hub location problem with fixed costs, the p-hub center problem, and the hub covering problem.Based on these papers, the authors analyze the trends and provide a synthesis of the literature.
In order to make models that in a closer way reflect the real world, some scholars have studied the hub location problem with stochastic variables.Contreras et al. [19] study the hub location problem in which demands or transportation costs are stochastic, and they prove that stochastic problems with uncertain demands or dependent transportation costs are equivalent to a deterministic expected value problem (EVP).The authors propose a Monte-Carlo simulationbased algorithm, combined with a Benders decomposition algorithm, to solve the model.Louveaux and Peeters [20] study how to convert the uncapacitated facility location problem into a dual stage stochastic problem.The authors propose a solution to the dual-based procedure.Ball and Lin [21] study the emergency services facilities point location problem from the perspective of system reliability.They convert the stochastic problem of the model into a 0-1 integer programming problem and solve it by the branchand-bound method.Albareda-Sambola et al. [22] convert the stochastic facility location and routing problem into a bilevel programming problem and propose a two-phase heuristics solution.Zhou and Liu [23] propose a hybrid intelligence algorithm which combines the stochastic simulation and genetic algorithm.This new algorithm can effectively solve the chance-constrained model and expected value model of facility location and distribution problem under the circumstances of stochastic demands.Gelareh and Pisinger [24] propose a novel mathematical model to solve the hub location problem of liner shipping companies.They assume the demand is elastic and then propose a mixed integer linear programming formulation for the simultaneous design of network and fleet deployment.
Expanding the current research on facility location that only considers the demands as either stochastic or fuzzy, Wen and Kang [25] view the demands as being stochastic and fuzzy simultaneously.The authors propose the expected cost minimization model, (, )-cost minimization model, and chance maximization model and develop a hybrid genetic algorithm to solve these models.Gabor and Van Ommeren [26] propose a 2-approximation algorithm to solve the facility location problem when the demands and inventory are stochastic.The numerical example proves that their method can effectively obtain a near optimal solution.Wang and He [27] propose a robust optimization model using the formation of regret model on the basis of the stochastic optimization model and provide an enumeration method and genetic algorithm to solve the model.The numerical example demonstrates that their method is better than the stochastic optimization model.Yang et al. [28,29] propose a method which combines the multiobjective fuzzy optimization model with a genetic algorithm to solve the fire station location problem.The model fully incorporates the fuzzy nature of decision-makers and the different fire demands of the different regions.In their model, the authors consider that construction costs, circulation costs, and customer demands are all fuzzy variables.The authors propose a logistics distribution center location model under such a fuzzy environment, establish a fuzzy chanceconstrained model, and combine the Tabu search algorithm with the genetic algorithm to solve the model.
Introducing stochastic variables into the location model renders the decision results more accurate.However, there are still several influential factors, such as policy environments, hydrological conditions, and commercial potential, that cannot be defined as a variable, or a stochastic variable.Therefore, the location model with stochastic programming still does not take many factors into account.Furthermore, many factors cannot be quantified.On the other hand, modeling real-life situations will introduce many quantitative and qualitative factors, and the models will become more complicated [30].In order to improve the quality of location decision, we should combine both qualitative and quantitative factors in a model.
Expert assessment methods are the most widely used methods to measure qualitative factors.Among them, the group decision-making method (GDM) is one of the most effective and efficient ways.In order to make decision results more accurate, researchers began to use fuzzy set theory to study MAGDM problems.As an important component of MAGDM, fuzzy multiattribute group decision-making (FMAGDM) is a difficult and hot research area [11,31].
Herrera et al. [32][33][34] have conducted a number of studies on multiattribute group decision-making problems based on language preference relationship.They use fuzzy sets to represent language variables and convert the variables into triangular or trapezoidal fuzzy numbers.Herrera et al. [32] and Herrera and Martínez [35] propose a method called linguistic 2-tuple representation model to represent the language variables based on the above fuzzy model.
Karacapilidis [36] find that the decision-makers may establish a common belief in the decision process by following a series of well-defined communicative actions.Then, they defined these actions and take them into consideration in the group decision process, which make the decision results more reliable.Lahdelma and Salminen [37] propose a stochastic multiattribute method which effectively provides decision-making support for decision-makers in a group decision.The paper proposes a SMAA (Stochastic Multicriteria Acceptability Analysis), an approach which can obtain a better compromise solution compared with that of traditional methods.Fenton and Wang [38] use linguistic variables and triangular fuzzy numbers to describe the opinions and attitudes of decision-makers and propose a fuzzy group decision-making approach based on risk and credible analysis.The numerical example shows that their method can effectively solve the inaccuracy and subjectivity problem in the decision-making system.Hochbaum and Levin [39] improve the existing methods of group decisionmaking.They present a new paradigm using an optimization framework that provides a specific performance measure for the quality of the aggregate ranking as per its deviations from the individual decision-makers' rankings.The new model is based on rankings provided with intensity; that is, the degree of preference is quantified.The model allows for flexibility in decision protocols and can take into consideration imprecise beliefs, less than full confidence in some of the rankings, and differentiating between the expertise of the reviewers.This method is simpler and more efficient than other existing methods and can be effectively applied in large-scale multiattribute decision-making problems.
The above studies optimize and improve the MAGDM method through different ways.Most of them assume that the attribute weights or DM weights are already known.However, because every project has its own specific influential factors and experts from different fields have their individual specific background knowledge, skills, and experience, subjective assignment of values to attribute weights and DM weights will increase the uncertainty of the decision-making process.
For the determination of experts' weights, French Jr. [40] measures the decision-makers' relative importance through the influence degree among the DMs.Bodily [41] establishes another decision-making team to give weights to the initial decision members and work out the weights of the initial decision members by measuring the additional preference value deviation.Yue et al. [42] decide the weights of attributes by converting attribute values into intuitive fuzzy numbers.Parreiras et al. [43] propose a flexible consensus scheme to establish the orders.Yue [44,45] proposes a calibration value method to calculate weights.Xu [46] improves Bodily's method and proposes a more direct method to calculate weights.Mirkin and Fishburn [47] give two eigenvector methods to measure the relative importance of decisionmakers.Van den Honert [48] uses the REMBRANDT (ratio estimation in magnitudes of deci-bells to rate alternatives which are nondominated) method, which is the combination of AHP and SMART (the sample multiattribute rating technique) to quantify the decision-makers' ability.Ramanathan and Ganesh [49] present a simple method which uses decision-makers' own subjective opinions to calculate the decision-makers' weights.
As shown above, current studies usually obtain attribute weights and experts' weights in one of three ways: (1) both attribute weights and experts' weights are assumed to be known; (2) attribute weights are assumed to be known, but experts' weights are calculated; (3) experts' weights are assumed to be known, but attribute weights are calculated.Little research has considered the situation that attribute weights and experts' weights are both unknown.
In order to make the calculation of qualitative factors more objective, this paper proposes an improved MAGDM method which can overcome the problem that commonly exists in previous studies (i.e., the subjective assignment of weights for attributes or experts in advance).Our method can effectively reduce the uncertainty caused by subjective weights for attributes or experts in the decision-making process and makes the calculation of qualitative factors more precise.
In sum, we solve the location problem of public transportation hubs in the following way.First, we obtain the values of qualitative factors at each alternative hub through our improved MAGDM method.Second, the quantitative factors, such as land cost and transfer distance, are considered in the model.Third, a stochastic programming is built to model the public transportation hub location problem where the passenger flow is a stochastic variable.Finally, we use the stochastic programming theory to solve the model and get the optimal location result.

The Stochastic Location Model
This section discusses the combination of quantitative and qualitative factors and establishes a public transportation hub location model under the condition that passenger flow is a stochastic variable.(Passenger flow denotes the number of passengers from one point to another.For instance, the passenger flow between points  and  is the number of passengers who choose the route  to .) Suppose there are  passenger generation points, where the coordinate of point  is (  ,   ) and  alternative places for public transportation hubs, where the coordinate of alternative hub  is (  ,   ).  denotes the maximum service capacity at alternative hub .  denotes the transportation cost of passengers from passenger generation point  to the alternative transportation hub .  is the passenger flow from passenger generation point  to the alternative transportation hub .  is the generalized distance from passenger generation point  to the alternative transportation hub .  is the fixed cost for building an alternative transportation hub , including setup cost and land cost.  denotes the number of passengers generated at point .  denotes the evaluation value of alternative hub , which reflects the influence of qualitative factors and is determined by a proposed MAGDM approach, as discussed in the next section: We construct the mathematical model for the problem as follows: Subject to: = 1 or 0,  = 1, 2, 3, . . ., .
Objective function ( 2) is to minimize the total costs.Function ( 3) is to maximize the total benefits of building hubs.Constraint (5) ensures that the designed passenger flow between  and  should be greater than or equal to the real passenger flow; namely, even in rush hour when passenger flow is greater on some routes, passengers can still choose another hub that is not necessarily the nearest.And constraint (6) limits the capacity of each alternative transportation hub.Because passenger flow is generalized from  and transfers to , and  is the alternative, constraints (5) and ( 6) are the key differences between distribution center location problems and classic facility location problems.
Traditional models assume that the number of passengers for each passenger generation point is given; however, passenger flow is usually unpredictable.Ignoring this uncertainty, transportation planning and construction would involve significant risk.Therefore, in this study, we consider that passenger flow is unknown, which means that  = { 1 ,  2 , . . .,   } is a stochastic vector in the model.As a result, our model becomes a stochastic programming problem.
The procedure to solve this model is as follows.First, we transform the above stochastic programming model into a chance-constrained programming model [23,58,59].Theorem 1.Let  be a decision vector,  a stochastic vector, (, ) the objective function, and   (, ) ( = 1, 2, . . ., ) the stochastic constraint function; the chance constraint form will be as follows: Suppose the stochastic vector  = {  [23,60].
Proof.According to the assumptions, we can rewrite Pr{(, ) ≤ 0} ≥  as follows: Obviously, for each given confidence level  (0 ≤  ≤ 1), there must exist a number   to meet the following equation: If a smaller number to replace   is used, the probability Pr{  ≤ } at the left will become higher.Therefore, Pr{ℎ() ≤ } ≥  is true if and only if ℎ() ≤   .Notice that the equation Pr{  ≤ } = 1 − Φ(  ) is always true.So we get   = Φ −1 (1 − ), where Φ −1 is the inverse function of Φ. Sometimes, the function Φ −1 may have multiple solutions.In this case, we will choose the largest one: Thus, we can get the equivalence of (10), which is as follows: We use Theorem 1 as the basis for solving the model.In constraint (5) of model ( 2), if   is a stochastic vector, constraint condition ∑  =1   ≥   is not meaningful.Therefore, if the probability of the constraint conditions being true is not smaller than the confidence level , the stochastic programming model can be transformed into a chanceconstrained model.The constraint condition ∑  =1   ≥   is converted to the following chance constraint: Suppose   follows normal distribution, and its probability distribution function is   ; for Constraint (5), according to Theorem 1, for a given confidence level , there is Pr{  (, ) ≥ 0,  = 1, 2, . . ., } ≥ .Let (, ) = ℎ() − , (, ) ≤ 0, and the chance constraint can be rewritten as follows: where ℎ() is the linear or nonlinear function of decision variable ,  is the stochastic variable, and probability distribution function is ().According to (11), for each given confidence level  (0 ≤  ≤ 1), there must exist a number   to satisfy the following equation: If using a larger number to replace   , the probability at the left will become higher.Therefore, Pr{ℎ() ≥ } ≥  becomes true if and only if ℎ() ≥   .Notice that the equation Pr{  ≥ } = () is always true.Thus, we get   =  −1 (), where  −1 is the inverse function of .According to (12), we will choose the largest one: Thus, we can get the equivalence of (15), which is as follows: Because Pr{  ≥ } = , then   =  −1 ().When ∑  =1   ≥   ,  = 1, 2, . . ., , then Pr{∑  =1   ≥   ,  = 1, 2, . . ., } ≥ .Therefore, constraint (5) ∑  =1   ≤   ,  = 1, 2, . . ., , can be converted to Thus, the original model can be transformed to If we know each alternative hub 's evaluation value   , the optimal location result can be obtained by solving the above model.Because we consider both qualitative and quantitative factors in our model, the decision results should be more accurate and practical.In the next section, we introduce an improved MAGDM approach that can objectively evaluate the attribute and expert weights.The proposed method can avoid the subjective weight assignment and reduce decision uncertainty.

Determining Each Alternative's Evaluation Value Using an Improved MAGDM Approach
This section introduces the proposed MAGDM method which can objectively assign weights to attributes and experts when these values are completely unknown in advance. . . .
be the decision matrix of th DM.First, we normalize   as   : By the given weight vector of attributes  = ( 1 ,  2 , . . .,   )  , we can establish the weighted normalized decision matrix as follows: ,  ∈ . ( Then we can get a group decision matrix by the following equation: where  = ( 1 ,  2 , . . .,   )  is the weight vector of DMs,   = ∑  =1      .Denote Using (25), we can get the evaluation value of each alternative.Most research in MAGDM problems assumes  = ( 1 ,  2 , . . .,   ) or  = ( 1 ,  2 , . . .,   ) is known in advance.However, the attributes of the problem are often complicated.Thus, it is difficult to directly assign weights to the attributes.And experts from different fields have different background knowledge, skills, and experience.Therefore, subjectively determining the weights of attributes and experts will increase the uncertainty of the decision-making process.Currently, little research has objectively assigned weights to both attributes and experts simultaneously.This paper proposes an improved group decision-making method that combines standard and mean deviation with improved TOPSIS.The method can objectively assign weights to the attributes and experts simultaneously.In our method, first, the standard and mean deviation maximization method is used to determine the weight of each attribute.And then an improved TOPSIS method is used to calculate the expert weights.Finally, based on the attribute weights, expert weights, and the grades which the experts give to each alternative, the final decision can be obtained.The whole decision-making process avoids subjective weighing for experts and attributes.This method reduces the uncertainty of the whole decisionmaking process and increases decision accuracy in general.
In order to minimize the loss of linguistic information, we can extend the original language assessment scale to be continuous  = {  |  ∈ [−, ]} [54].If   ∈ , we call   an original linguistic term.Otherwise, we call   a virtual linguistic term [52,54].In general, the decision-maker uses the original linguistic terms to evaluate the alternatives, and the virtual linguistic terms can only appear in operation [54].
For expert   and attribute   , the standard deviation between alternative   and other alternatives is [62] where represents the mean value of the attribute   given by expert   .( ()   ,  ()  ) represents the deviation of the mean value  ()   to the attribute value  ()  of the alternative   for the attribute   of the expert   .So  ()  represents the standard deviation for the attribute   of expert   .
( ()  ,  ()  ) represents the deviation of the mean value  ()   to the attribute value  ()  of the alternative   for the attribute   of the expert   .So  ()  represents the mean deviation for attribute   of the expert   .
The choice of weighed vector  should make the total standard deviation and mean deviation of all evaluation indexes maximized.Therefore, the objective function can be constructed as follows: Subject to: where  ()  and  ()  denote the standard deviation and mean deviation for the attribute   of the expert   , respectively, and  and V denote the preference of the decision-makers;  = 0 means that the DMs only consider the mean deviation; V = 0 denotes that the DMs only consider the standard deviation; and  ̸ = 0, V ̸ = 0 denote that the DMs consider the standard deviation and mean deviation simultaneously.Thus, we have the following optimization model: Subject to: The solution of ( 38) can be obtained: Normalizing  ()  we obtain

Determining the DMs' Weights by an Improved TOPSIS
Method.The traditional TOPSIS method first calculates the closeness coefficient based on the distance between the positive ideal solution denoted as PIS  + and negative ideal solution denoted as NIS  − and then ranks the alternatives by closeness coefficient.Here, we describe an improved TOPSIS method to calculate the expert weights [44,45].
The decision matrix of th expert is Normally, there are benefit and cost attributes in a MAGDM problem; we next use (43) to normalize the attribute value   of decision matrix  = (  ) × to corresponding element   and then calculate the normalized matrix [63,64]: For benefit attributes  ()  ,  ∈ ,  ∈ , we can get For cost attributes  ()  ,  ∈ ,  ∈ , we can get Next, we can calculate the normalized decision matrix as follows based on the weight vectors of attribute  = ( 1 ,  2 , . . .,   ) calculated above: . . .
At average, the best decision result should be the average matrix of the group decision matrix that is PIS [65]: where By the utilitarian distance, the separation  +  of each individual decision matrix   from  * is given as follows: ,  ∈ .(48) Thus, the smaller the distance  +  , the better the decision   of th expert.Here, we divide the NISs into two parts L-NIS  −  and R-NIS  −  : where ≥  *  }.Actually,  −  and  −  are the minimum and maximum matrixes of the group decision, respectively.Similarly,  −  and  −  are given as follows: Obviously, the larger the value of  −  and  −  , the better the decision   of th expert.The closeness coefficient of the th expert with respect to  * is defined as follows: where  −  ≥ 0,  −  ≥ 0, and  +  ≥ 0 ( ∈ ).Obviously,   ∈ [0, 1] ( ∈ ).
By the closeness coefficient, we can rank the decisionmakers' order and get the weight vectors of DMs by transforming the closeness coefficient with the following equation: Therefore, 4.5.Determining the Evaluation Value.After the DMs and attributes' weights vectors are calculated based on the above analysis, the evaluation value of each alternative can be calculated by the following equation:

General Steps in Decision Analysis
Based on the above analysis, we propose a new transportation hub location method that considers both qualitative and quantitative factors.The steps of decision-making process are as follows.
Step 4. Use the LWAA operator to integrate the group decision-making information, and then apply (55) to calculate the evaluation value   of each alternative transportation hub.
Step 5. Use the evaluation value   of each candidate transportation hub and the given data to establish the stochastic location model (2).
Step 7. Solve model (20) and obtain the solution.

A Numerical Study
The evaluation values of four alternative hubs given by expert  2 .
environmental impact, and  5 , hydrological conditions.The language variable set is as follows: Three experts give their evaluation values in Tables 1, 2, and 3.
Finally, using (55) the evaluation values of the four alternatives can be got:    4 and 5.
The generalized transportation costs   of passengers from passenger generation point  to the alternative transportation hub  are shown in Table 6.Table 7 denotes the expected value , standard deviation , and the confidence level   of stochastic generated {  } of the passenger generated point.Here,   = 350 and   = [20,19,18,22] (unit: ten thousand Yuan, i.e., Chinese currency).
A mathematical model can be constructed using (20).The model can be solved in two ways.By assigning different values to  1 and  2 , we can obtain noninferior solutions. 1 ()'s dimension is "Yuan." 2 ()'s dimension can be considered as "satisfaction," which means the higher the evaluation value is, the more satisfaction the decision-maker has.If the decision-maker thinks that to increase  2 in satisfaction is equivalent to decreasing  1 Yuan in cost, the double objective problem can be converted into a single objective problem.This is because the satisfaction can be converted into cost [66].
Then, we can calculate the result as in Table 8.The decision-maker selects one project from Table 8 according to its own preferences.For example, if the decisionmaker chooses the 6th solution,  1 ,  2 ,  3 solution is selected.

When Decision-Maker's Preferences Are Unknown.
When there is no information about the objective function, finding the neutral compromise solution is a good way to get the final decision [67].Accordingly, we calculate the positive and negative ideal solution  + = (4675, 1.4395),  − = (11190, 1.1879) [68].Then, we can get the compromise programming model: where  mid  = ( −  +  +  )/2,  = 1, 2, represents the first and second objective,  represents the feasible region of (68), and  is a very small positive number.

Conclusions
This paper proposes a novel approach to solving the public transportation hub location problem, by combining stochastic programming with an improved multiattribute group decision-making method and fully incorporating the relevant quantitative and qualitative factors.Firstly, we define the passenger flow as a stochastic variable and establish the stochastic location model to figure out the influence of quantitative factors.Secondly, we propose an improved MAGDM approach to figure out the influence of qualitative factors, in which standard and mean deviation method is used to determine the attributes' weights, and TOPSIS method is used to determine the DMs' weights.Thirdly, combining quantitative calculation with the qualitative MAGDM method, the optimal alternative hub is computed.Moreover, a numerical example is given to verify the developed approach and to demonstrate its practicality and effectiveness.It can be inferred that further research may focus on the application of the method in other areas and introduce more stochastic variables in the model.For instance, prospective research may treat construction and transportation costs as stochastic variables.

6 . 2 .
Location Calculation.After calculating the evaluation values of four alternative hubs, in this section we compute the final result of the location problem.First, suppose there are five passenger generation points in the four alternative

Table 3 :
The evaluation values of four alternative hubs given by expert  3 .

Table 4 :
The coordinates of the five passenger generated points.

Table 5 :
The coordinates of the four alternative hubs.

Table 7 :
The stochastic distribution of passenger generated points. 1 ,  2 ,  3 ,  4 , and  5 .The problem is how to select the hubs in the four alternative hubs ( 1 ,  2 ,  3 ,  4 ) to minimize the total cost.The coordinates of the alternative hubs and the passenger generated points are shown in Tables

Table 8 :
The noninferior solutions in different combinations of weights.