^{1}

^{*}

^{1}

^{2}

Due to the problem complexity, simultaneous solution methods are limited. A hybrid algorithm is emphatically proposed for LRP. First, the customers are classified by clustering analysis with preference-fitting rules. Second, a chaos search (CS) algorithm for the optimal routes of LRP scheduling is presented in this paper. For the ergodicity and randomness of chaotic sequence, this CS architecture makes it possible to search the solution space easily, thus producing optimal solutions without local optimization. A case study using computer simulation showed that the CS system is simple and effective, which achieves significant improvement compared to a recent LRP with nonlinear constrained optimization solution. Lastly the pratical anlysis is presented relationship with regional logistics and its development in Fujian province.

All companies that aim to be competitive on the market have to pay attention to their organization related to the entire supply chain. In particular, companies have to increase the efficiency of their logistics operation. It is essential that reasonable optimal systems of logistics distribution for enterprises with the development of e-commerce and integrated logistics. Location-routing problem (LRP) is an important branch of routing optimization in integrated logistics systems, which has been solving for every logistics distribution corporations.

The conceptual foundation of LRP studies dates back to Von Boventer (1961), Maranzana (1965), Webb (1968), Lawrence and Pengilly (1969), Christofids and Eilon (1969) and Higgins (1972) [

Chaotic system is a common phenomenon exiting in nonlinear dynamical system with popular characteristic of dynamics. Due to the ergodicity and randomness of chaotic sequence, the chaos search (CS) architecture can be solved the traditional optimization [

The main contribution of the paper is statement of a new hybrid algorithm combined clustering analysis and chaos search, together with a proof of their correctness [

Section 2 introduces the meaning and descriptions of LRP and describes the mathematical model of LRP. Section 3 emphasizes on the hybrid algorithm based on clustering analysis and chaos search for LRP. Section 4 shows the computational results and analysis. Finally, the last section gives conclusions and future work.

Location-routing problem (LRP) is one of the problems in integrated logistics optimizations. The LRP can be defined as follows. A feasible set of potential facility sites and location; and expected demands of each customer are given. Each customer is to be assigned to a facility, which will supply its demand. The shipments of customer demand are carried out by vehicles, which are dispatched from the facilities and operated on routes that include multiple customers. There is a fixed cost associated with opening a facility at each potential site, and a distribution cost associated with any routing of vehicles that includes the cost of acquiring the vehicles used in the routing, and the cost of delivery operations. The cost of delivery operations is linear in the total distance traveled by the vehicles [

A constraint-based model is presented in this paper for the location routing problem. The hypotheses are as follows:

The transportation is just in time.

The facility is both starting point and destination of circular vehicle routing; each facility serves more than two customers.

Nature of demand/supply: deterministic.1 There are multiple facilities.

Size of vehicle fleets: single vehicle (A facility is served by one vehicle; meanwhile, it is satisfied that the total demand of the route is no less than one vehicle’s capacity).

Vehicle capacities are determined. The total quantity (amount) of goods is limited in every route by each vehicle’s capability.

Facility capacities are undetermined; not all facilities have been chosen in every decision.

Each customer is served by one and only one vehicle. Considering the complexity of supply/demand markets, it is assumed that they are retail markets.

Each facility is considered a separate entity, not linked to the other facilities.

Objective: minimize total costs.

is the set of m feasible sites of candidate facility,

is the set of N customers to be served,

is the set of all feasible sites and customers (it is also referred to nodes),

is the set of K vehicles available for routing from the facilities, average annual cost per distance traveling from

annual cost of establishing and operating a facility at site, average number of units demands by customer capacity of vehicle distance from node i to node j, traveling by vehicle k from node r to customer i or to customer j, respectively.

The objective function of LRP is

Subject to

In this model, the objective function minimizes the total cost of routing, and establishing and operating the facilities. Constraint (2) ensures that each customer is served by one and only one vehicle. Constraint (3) ensures that the vehicle capacity constraints are not exceeded for any of the vehicles used in routing, while (4) is the route continuity constraint, which implies that the same vehicle should leave every point that is entered by a vehicle. Constraint (5) guarantees that each vehicle is routed from one depot. Constraint (6) guarantees that there are no links between any two depots. Constraints (7) and (8) require that a vehicle can only from a depot if that depot is opened. The last two of constraints are the integer constraints.

Clustering analysis is one of the major techniques in pattern recognition. It is a technique for grouping data and finding where similar data are assigned to the same cluster whereas dissimilar data should belong to different clusters. The conventional clustering methods generate partitions, in a partition; each pattern belongs to one and only one cluster.

The customer aggregate C is divided to cluster by the euclidean measure. The steps of optimal facility location are performed in the following way, which have extensively analyzed and used in location theory to approximate distances between two spatial coordinates. For the sake of simplicity, it is assumed that each facility has its own set m of their coordinates of candidate sites for facility location, and no two facilities share the same location, which are denoted by

Every two coordinates of candidates sites for facility locations are denoted by and, and which satisfied with

where R_{ij} is one half of the distance between and.

Rules for classification are given as follows.

Rule1 (12)

Ｃ_{p} can only be divided into PF_{i} if formula (12) is true.

Rule 2 The coordinates of customers on the circumferences are satisfied that the distance to the nearest two customers is less than that to the sub nearest two points. sub-nearest meaning that they are not included a cluster by reason of distance).

Rule 3 If (13)

where r_{jq} is the distance from any customer C_{q} to any facility PF_{j}, maxR_{ij} is half of the most distance between any two facilities. Either these customers will not be served, or else the candidate’s facility location will be abolished. Last two rules only relate to abnormal customers points.

where Ｎ is the numbers of customers,

are expressed the points and their coordinates, respectively.

r_{ip} is the distance between any candidate facility PF_{i} and any customer p. The tangent cycles O_{1}, O_{2}, …, O_{i}, O_{j}, …, O_{m} are drawn with the radius R_{ij} , where R_{1}, R_{2}, …, R_{i}, R_{j}, …, R_{m} are their radiuses, respectively.

The fundamental idea is that chaotic sequences belonging to the ergodic points in [0,1] are constructed by the chaos system [

where is the so called driving parameter(a parameter capturing such factors as fertility rate) initial living area, etc. The initial value, or seed x is in., There if From about x is restricted to and to.

The characteristic of Logistic map is given as follow:

1) When λ to, there are stable points in Logistic map.

2) When λ to, two sequences with seeds only slightly apart diverge very rapidly, that is to say, there exists in the period-doubling bifurcation of the solutions.

3) When λ to, chaos start to come out, and Logistic map isn’t in stable cycling orbit. It is difficult to decide which points with Lebesgueis more than zero in chaotic sequences.

4) When λ to 4, there is a stable chaotic sequence.

In the stable chaotic sequence, any orbit f is dense, that is to say, the points of one orbit f is included in open ball, when and. On the contrary, all points that belong to are pressed on towards any orbit f as precision as possible. In this paper when, the population sequence is chaotic and ergodic and gives search optimal solutions (or near optimal solutions) solutions. A chaotic system is defined as one where two different, but close, seeds result in widely different population sequences. An ergodic system is one where the sequences will approach every possible point arbitrarily, little by little (more formally, for every pair of element of, there is an xk between them). Thus, an ergodic system is one that is “trying to converge to everything”. The Logistic model can be solved not only for the continuous variable optimization, but also discrete variable optimization.

The hybrid algorithm includes two steps. First, preference-fitting rules to solve location-allocation cluster the customers. Second, the optimal routes is obtained by chaos search to be minimized the function. The Logistic chaotic model is transformed for discrete optimization of LRP [

where

denoted the distance between customer and customer.

Steps with the hybrid algorithm based on clustering analysis and chaos search are described as follow.

Step 1 Customers are clustered by preference-fitting rules.

1) Initial calculations for R_{ij} and r_{ip }are performed.

2) If r_{ip} < R_{ij}, the customer p is divided into cluster O_{i}.

3) Otherwise, return step (1.1).

Step 2 is given by initial random value. The best optimal value of the function, given with the control parameter k gifted with a larger value.

Step 3 Put the results of order arranged into. It is obtainable to a current searching route, and to a current legal route by hybrid algorithm, compared with.

Step 4 Calculate the function values f of the current legal route.

Step 5 Remember the minimum function values and and also remembering the present legal route.

Step 6 Decide whether the cycle is stopped or not. If step = k, stop.

Step 7 A set of is generated by Logistic map, step = step + 1.

Step 8 Return steps 3, and cycling.

In the following example three potential facilities, and ten customs were analyzed. Each was served. With the same type of distribution vehicle that had a carrying capacity equal to 30 ton. The annual cost establishing and operating a facility F_{r} is equal to 160 yuan/each, The transportation cost C_{ij} is equal to 2 yuan per ton, per km. The hybrid algorithm has been discussed to search the optimal route based on the data in

The simulation results are shown in

The empirical analysis is given in relationship between with regional logistics and regional economy. The model

with analysis is proposed for the relationship with modern logistics and GDP. The factors are chosen include GDP, logistics production value (WLCZ), freight turnover (HZL), logistics mileage (WLC), and which the data is chosen from 1978 to 2010. The unchangeable price in1978 is given with each variable. The quantitative software is accepted with EVIEWS 6.0, and which the trend of each variable is give in

Test type (C, T, K), The letters of C, T, K means intercept, Trend, lagging numbers. The first four are contains intercept and trend; The after four are not contains trend.

From the

The Johansen co-integration test is accepted for GDP, WLCZ, HZL, WLC. This test results is presented in

From

Note: The letters of r means numbers of co-integration equation, ^{*}means reject null hypothesis.

This demonstrates that the variable is at least one cointegration equation between the relationships, that there is a long-term equilibrium relationship between variables.

(2.22) (0.88) (0.24) (17)

The above co-integration equation reflects the longrun equilibrium relationship between variables, but cannot fully reflect the short-term changes relationship between variables. Therefore, it is necessary to establish the error correction model (ECM), to estimate short-term changes between variables. Estimated ECM is give as following:

(3.01) (6.84) (1.90) (–2.58) (18)

From the equation, T statistics test of all variables are significant in level 10%, the error correction of ECMt-1 coefficient is negative, accord with reverse revision mechanism, thus further shows that the error correction model is good, with residual sequence is the white noise.

By comparing the co-integration equation and ECM, We can conclusion that t compared on short-term, freight turnover (HZL) have negative influence on the GDP in long term, but its elasticity coefficient is greater than the short-term elastic coefficient; Output value of logistics (WLCZ) have short-term elastic significantly greater than its long-term flexibility to GDP, and both have different symbols, changes of direction that does not accordance, it indicates that in recent years the logistics infrastructure investment lead to the logistics industry close to a saturation stage in Fujian province, although output value of logistics (WLCZ) promote the growth of GDP in short-term, but due to the expansion of the regular, the GDP growth plays a role of obstacles in the long term.

In this present paper, it has been proposed that the effective method for LRP using a hybrid algorithm based on clustering analysis and chaos search. Through computer simulations the method has been shown to solve largescale LRP. Using the clustering technique, the computational cost is dramatically reduced. Due to the feature of chaotic global ergodicity, the optimal routes are randomly seeked in one cluster, which ensures search efficiently and produces optimal solutions (or near optimal solutions) without local optimization. Thus, the emphasis is directed towards developing a methodology, which has been improved the quality of solution.

The application to other combinatorial optimization problems should be investigated in integrated logistics. The relationship is given with regional logistics and its development. There exists a long-term equilibrium relation between economics growth and logistics development in Fujian province.

This paper was supported in part by National Natural Science Foundation of China under grant number 71040009, in part by Program for New Century Excellent Talents in University under grant number NCET-10- 0118, in part by Youth Foundation of Chinese Education by HuoYing-dong under grant number 104009.