本发明涉及物流技术领域，尤其涉及一种物流配送中心选址方法、装置及计算机可读存储介质。
背景技术：
物流业被人们看作在本世纪经济发展的新领域，是最重要的新经济增长点之一，也被称之为货物运输企业在经营方式发展上的高级阶段。
目前关于配送中心的选址方法普遍是采用偏微分方程法、重心法和微分法进行选址，上述方法对物流网络的规模的限制要求较高，难以对较大规模的物流网络配送中心选址问题进行求解。
上述内容仅用于辅助理解本发明的技术方案，并不代表承认上述内容是现有技术。
技术实现要素：
本发明的主要目的在于提供一种物流配送中心选址方法、装置及计算机可读存储介质，旨在解决现有技术中，难以对大规模物流网络配送中心选址问题进行求解的问题。
为实现上述目的，本发明提供一种物流配送中心选址方法，所述方法包括以下步骤：
获取待选址交通网络中的物流配送点，所述物流配送点包括n个预设配送中心点，所述n为大于等于1，且小于所述物流配送点个数的整数；
获取各预设配送中心点对应的无向图，并根据所述无向图获取各预设配送中心点对应的可行路线集，所述可行路线集包括若干条可行配送路线；
对各预设配送中心点对应的可行路线集进行筛选，将满足预设条件的可行配送路线作为目标路线，并根据所述目标路线获取对应的物流配送中心点地址。
优选地，所述获取各预设配送中心点对应的无向图，并根据所述无向图获取各预设配送中心点对应的可行路线集的步骤，具体包括：
获取各预设配送中心点对应的无向图，根据所述无向图利用膜粒子群算法获取各预设配送中心点对应的可行路线集。
优选地，所述膜粒子群算法为基于膜计算的粒子群算法；
相应地，所述根据所述无向图利用膜粒子群算法获取各预设配送中心点对应的可行路线集的步骤，具体包括：
对各预设配送中心点进行遍历，获取包含遍历到的目标配送中心点的目标无向图；
根据所述目标无向图，利用所述膜粒子群算法获取满足预设约束条件的s条放射状配送路线，所述s为不小于1的整数；
从所述s条放射状配送路线中选取一条配送路线，获取包含被选取的配送路线的环状路线；
根据所述环状路线，获取所述环状路线中可删除路线的集合；
从所述集合中选取一条可删除路线，并将其在所述环状路线中对应的路线删除，获得被选取的配送路线的可行路线；
从剩余放射状配送路线中选取配送路线并获取对应的可行路线，直至获取到s条放射状配送路线对应的可行路线；
根据获取到的s条放射状配送路线对应的可行路线，得到所述配送中心点对应的可行路线集；
选取剩余的配送中心点，并获取对应的可行路线集，直至获取各预设配送中心点对应的可行路线集。
优选地，所述根据所述目标无向图，利用所述膜粒子群算法获取满足预设约束条件的s条放射状配送路线的步骤，具体包括：
根据所述目标无向图，利用所述膜粒子群算法，以物流配送中心选址为研究对象，建立对应的数学模型；
根据所述数学模型，获取满足预设约束条件的s条放射状配送路线。
优选地，所述根据所述数学模型，获取满足预设约束条件的s条放射状配送路线的步骤，具体包括：
根据所述数学模型，以所述目标无向图对应的预设配送中心点为起点，搜索满足预设约束条件的待配送客户点，并将所述待配送客户点加入到集合t中；
从所述集合t中选取一个待配送客户点加入到集合q中，将选取到的待配送客户点从所述集合t中删除；
以所述选取到的待配送客户点为起点，搜索满足预设约束条件的相邻待配送客户点，并将所述相邻待配送客户点加入到所述集合t中；
返回从所述集合t中选取一个待配送客户点加入到集合q中的步骤，直至所述目标无向图中所有满足预设约束条件的待配送客户点均加入到所述集合q中，所述集合t为空集；
根据最终获得的集合q，获取满足预设约束条件的s条放射状配送路线。
优选地，所述从所述s条放射状配送路线中选取一条配送路线，获取包含被选取的配送路线的环状路线的步骤，具体包括：
从所述s条放射状配送路线中选取一条配送路线，并获取被选取的配送路线对应的待加入路线集，从所述待加入路线集中选取一条路线加入所述被选取的配送路线，获取包含所述被选取的配送路线的环状路线。
优选地，所述根据所述环状路线，获取所述环状路线中可删除路线的集合的步骤，具体包括：
根据所述环状路线，从所述环状路线中预设配送中心点所在位置出发，使用深度优先搜索获取所述环状路线所在的位置，根据所述位置获取从所述待加入路线集中选取的路线的节点；并获取所述节点在所述s条放射状配送路线中对应的位置信息；
根据所述位置信息，采用回溯法获取所述环状路线中可删除路线的集合。
优选地，所述对各预设配送中心点对应的可行路线集进行筛选，将满足预设条件的可行配送路线作为目标路线的步骤，具体包括：
根据预设条件对各预设配送中心点对应的可行路线集进行初次筛选，获取每个可行路线集中满足所述预设条件的可行路线；
获取所述初次筛选后的可行路径，从所述可行路径中再次筛选出满足预设条件的可行配送路线，并将所述可行配送路线作为目标路线。
此外，为实现上述目的，本发明还提出一种物流配送中心选址装置，所述装置包括：存储器、处理器及存储在所述存储器上并可在所述处理器上运行的物流配送中心选址程序，所述物流配送中心选址程序配置为实现如上文所述的物流配送中心选址方法的步骤。
此外，为实现上述目的，本发明还提出一种计算机可读存储介质，所述计算机可读存储介质上存储有物流配送中心选址程序，所述物流配送中心选址程序被处理器执行时实现如上文所述的物流配送中心选址方法的步骤。
本发明通过获取待选址交通网络中大量的物流配送点，根据所述物流配送点中的预设配送中心点，获取各预设配送中心点对应的无向图，根据所述无向图获取各预设配送中心点对应的可行路线集，并对各预设配送中心点对应的可行路线集进行筛选，将满足预设条件的可行配送路线作为目标路线，并根据所述目标路线获取对应的物流配送中心点地址，由于通过获取每个预设配送中心点对应的可行路线集，进而获取可行路线集中满足预设条件的可行配送路线对应的配送中心点的地址，从而解决了难以对大规模物流网络中选取配送中心地址的问题进行求解的问题，提高了物流配送中心选址的效率。
附图说明
图1为本发明实施例方案涉及的硬件运行环境的物流配送中心选址装置结构示意图；
图2为本发明物流配送中心选址方法第一实施例的流程示意图；
图3为本发明物流配送中心选址方法第二实施例的流程示意图；
图4为本发明物流配送中心选址方法第二实施例中膜结构的结构框图；
图5为本发明物流配送中心选址方法第二实施例实际应用流程图；
图6为本发明物流配送中心选址方法第二实施例实际应用中放射状配送路线示意图。
本发明目的的实现、功能特点及优点将结合实施例，参照附图做进一步说明。
具体实施方式
应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
参照图1，图1为本发明实施例方案涉及的硬件运行环境的物流配送中心选址装置结构示意图。
如图1所示，该物流配送中心选址装置可以包括：处理器1001，例如cpu，通信总线1002、用户接口1003，网络接口1004，存储器1005。其中，通信总线1002用于实现这些组件之间的连接通信。用户接口1003可以包括显示屏(display)、输入单元比如键盘(keyboard)，可选用户接口1003还可以包括标准的有线接口、无线接口。网络接口1004可选的可以包括标准的有线接口、无线接口(如wi-fi接口)。存储器1005可以是高速ram存储器，也可以是稳定的存储器(non-volatilememory)，例如磁盘存储器。存储器1005可选的还可以是独立于前述处理器1001的存储装置。
本领域技术人员可以理解，图1中示出的物流配送中心选址装置结构并不构成对物流配送中心选址装置的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件布置。
如图1所示，作为一种存储介质的存储器1005中可以包括操作系统、数据存储模块、网络通信模块、用户接口模块以及物流配送中心选址程序。
所述物流配送中心选址装置可以是能够进行物流配送中心选址计算及程序运行的装置，例如，电脑等，本实施例对此不加以限制。
在图1所示的物流配送中心选址装置中，网络接口1004主要用于与后台服务器进行数据通信；用户接口1003主要用于与用户进行数据交互；本发明物流配送中心选址装置中的处理器1001、存储器1005可以设置在物流配送中心选址装置中，所述物流配送中心选址通过处理器1001调用存储器1005中存储的物流配送中心选址程序，并执行以下操作：
进一步地，处理器1001可以调用存储器1005中存储的物流配送中心选址程序，还执行以下操作：
本实施例的有益效果是：通过获取待选址交通网络中大量的物流配送点，根据所述物流配送点中的预设配送中心点，获取各预设配送中心点对应的无向图，根据所述无向图获取各预设配送中心点对应的可行路线集，并对各预设配送中心点对应的可行路线集进行筛选，将满足预设条件的可行配送路线作为目标路线，并根据所述目标路线获取对应的物流配送中心点地址，由于通过获取每个预设配送中心点对应的可行路线集，进而获取可行路线集中满足预设条件的可行配送路线对应的配送中心点的地址，从而解决了难以从大规模物流网络中选取配送中心地址的问题，提高了物流配送中心选址的效率。
基于上述硬件结构，提出本发明物流配送中心选址方法实施例。
参照图2，图2为本发明一种物流配送中心选址方法第一实施例的流程示意图。
本实施例中，所述方法包括以下步骤：
步骤s10：获取待选址交通网络中的物流配送点，所述物流配送点包括n个预设配送中心点，所述n为大于等于1，且小于所述物流配送点个数的整数；
步骤s20：获取各预设配送中心点对应的无向图，并根据所述无向图获取各预设配送中心点对应的可行路线集，所述可行路线集包括若干条可行配送路线；
需要说明的是，所述物流配送点除预先设定的配送中心点外，还包括众多的待配送客户点，在获取到预设配送中心点后，可从预设配送中心点中选取一个配送中心点，然后将该被选取的配送中心点与待配送客户点连接，获取包含该被选取的配送中心点对应的无向图，按照同样的方法分别获取各个预设中心点对应的无向图，然后根据所述无向图获取各预设配送中心点对应的可行路线集，所述可行路线集包括有若干条可行配送路线，所述可行配送路线可以是满足一定配送条件的物流配送线路。
步骤s30：对各预设配送中心点对应的可行路线集进行筛选，将满足预设条件的可行配送路线作为目标路线，并根据所述目标路线获取对应的物流配送中心点地址。
在具体实现中，可根据预设条件对各预设配送中心点对应的可行路线集进行初次筛选，获取每个可行路线集中满足所述预设条件的可行路线；获取所述初次筛选后的可行路径，并从所述可行路径中再次筛选出满足预设条件的可行配送路线，将所述可行配送路线作为目标路线，此处，所述预设条件可以是路线距离最短，也可以是运输成本最小，本实施例对此不加以限制，当筛选出满足预设条件的可行配送路线后，所述可行配送路线就是我们寻求的目标路线，即最佳配送路线，然后再根据最佳配送路线对应的配送中心点所在的位置，就可以确定出物流配送中的最佳地址。
本实施例通过获取待选址交通网络中大量的物流配送点，根据所述物流配送点中的预设配送中心点，获取各预设配送中心点对应的无向图，根据所述无向图获取各预设配送中心点对应的可行路线集，并对各预设配送中心点对应的可行路线集进行筛选，将满足预设条件的可行配送路线作为目标路线，并根据所述目标路线获取对应的物流配送中心点地址，从而解决了难以对大规模物流网络中选取配送中心地址的问题进行求解的问题，提高了物流配送中心选址的效率。
参照图3，图3为本发明物流配送中心选址方法第二实施例的流程示意图。
基于上述第一实施例，为进一步提高物流配送中心选址问题的求解效率与精度，在获取各预设配送中心点对应的无向图之后，可根据所述无向图利用膜粒子群算法获取各预设配送中心点对应的可行路线集。
需要说明的是，膜粒子群算法是采用膜计算和粒子群算法相结合的一种算法，所述膜计算是从细胞的结构和功能中抽象出来的一种新的计算模型，它的运行机理是运用细胞膜的特性将整个环境划分为多个不同的区室，每个小的区域(类似细胞)根据细胞的功能进行不同的生化反应，细胞之间通过转运规则可以快速的交换区域内的物质，当所有细胞的操作都被执行之后，则完成系统的一个格局，随之膜系统进入下一格局的转换，反复重复上述操作，最后输出膜系统计算的最终的格局。
所述粒子群算法(particleswarmoptimization,pso)是一种模拟鸟群捕食行为的进化算法，位置和方向是个体的重要指标，依据全局最优和局部最优指引个体飞行，算法具有简单易行、操作方便和精度较高的特点。粒子的进化主要分为以下几个步骤：
(1)产生初始的粒子；(2)根据当前的速度，计算公式获取粒子在下一时刻的速度；(3)获取接收下一时刻粒子的速度，根据公式计算出粒子新飞行的位置；不断的执行上述步骤，直至所有粒子都更新完，结束。在传统的计算装置中，大量粒子的飞行位置的更新是以单个个体为单位的串行的方式更新粒子在下一时刻的位置，这使得粒子的进化速度较慢；因此，受到膜计算的高效性以及飞鸟集群活动的规律性启发，进而利用两种技术的结合建立一个新的算法-膜粒子群算法(又称膜算法)。它是一种基于膜结构的粒子群算法，即利用不同的细胞去指导粒子的进化。在对粒子活动行为观察的基础上，利用细胞膜中粒子对信息的共享使整个群体的运动在问题求解空间中产生以不断优化的次序进行演化，从而获得最优解。具体执行过程如下：
将不同的进化动作放在不同的细胞中进行，如一个细胞用于一次性产生所有的粒子(对象)，一个细胞用于更新当前所有粒子的速度，一个细胞用于更新下一时刻粒子飞行的位置，在更新时，用另一个细胞用于获取粒子飞行的可能位置的集合，这样每次执行以整体作为研究对象，从而证所有的粒子在同一时刻同时飞向新的位置。不断更新系统的格局，直至达到终止状态，结束。输出最优的粒子位置。
本实施例中，基于上述第一实施例，通过利用膜粒子群算法，并根据所述无向图获取各预设配送中心点对应的可行路线集的步骤s20具体包括：
步骤s201：对各预设配送中心点进行遍历，获取包含遍历到的目标配送中心点的目标无向图；
需要说明的是，为获取每个预设配送中心点对应的可行路线集，需要对各预设配送中心点进行遍历，并将每次遍历到的配送中心点作为目标配送中心点，获取该目标配送中心点对应的目标无向图，例如当第一次遍历到的配送中心点为v1，则将所述配送中心点为v1作为本次求解的目标配送中心点，获取目标配送中心点v1对应的目标无向图。
步骤s202：根据所述目标无向图，利用所述膜粒子群算法获取满足预设约束条件的s条放射状配送路线，所述s为不小于1的整数；
需要说明的是，在获取到包含遍历到的目标配送中心点的目标无向图之后，需根据所述目标无向图，利用所述膜粒子群算法，以物流配送中心选址为研究对象，建立对应的数学模型：
∏＝(o,σ1,σ2,σ3,σ4,syn,i0)
其中，o＝{xkj}是一个由0和1组成的输入字母表，0 参考图4，图4为本发明物流配送中心选址方法第二实施例中膜结构的结构框图。其中，细胞膜m1用于根据所述目标无向图，通过普利姆随机搜索规则，获取满足预设约束条件的s条放射状配送路线。 本实施例将所述预设约束条件定义为具有放射状以及长度限制的约束表达式： 其中，f(x)是目标函数，xij表示节点vi和vj之间是否连通，wij表示节点vi和vj的权值。 为使获得预设约束条件的放射状配送路线，对于表达式(1)获取的函数值还需要满足公式(2)和(3)： 上述公式(2)中，z表示目标无向图中所有物流配送点的个数，v表示初始时获取的交通网络中所有的物流配送点个数，t(xij)＝{xij|xij＝1}是满足放射状和长度约束的网络图中的边集。 上述公式(3)中，wij表示节点vi和vj的权值，p0表示辐射状分支上的最大路经值，存储从节点v1开始的一条路径中的所有节点。 本实施例中，根据所述数学模型，以所述目标无向图中预先选取的配送中心点v1为起点，搜索满足预设约束条件的待配送客户点，并将所述待配送客户点加入到集合t中；从所述集合t中选取一个待配送客户点加入到集合q中，将选取到的待配送客户点从所述集合t中删除；以所述选取到的待配送客户点为起点，搜索满足预设约束条件的相邻待配送客户点，并将所述相邻待配送客户点加入到所述集合t中；返回从所述集合t中选取一个待配送客户点加入到集合q中的步骤，直至所述无向图中所有满足预设约束条件的待配送客户点均加入到所述集合q中，所述集合t为空集；根据最终获得的集合q，获取满足预设约束条件的s条放射状配送路线。 步骤s203：从所述s条放射状配送路线中选取一条配送路线，获取包含被选取的配送路线的环状路线； 在确定被选取的配送路线后，需根据当前配送路线包含的各物流配送点位置信息，获取该配送路线对应的待加入路线集。 结合图4，在具体实现中，细胞膜m2用于从所述s条放射状配送路线中选取一条配送路线,获取被选取路线对应的待加入路线集。根据当前优化路径的特点和分析，细胞膜m2用于获取待加入路线集。膜中的规则的进化机制r1如公式(4)所示：其中，vi(t)是第i个粒子在第t时刻的速度分量。 r1≡[vi(t)]2→[vi(t+1)]2(4) 需要说明的是，所述公式(4)表示：在细胞膜m2中，第t时刻的所有粒子的速度分量vi(t)根据公式(11)中第一个速度更新公式计算出粒子在第t+1时刻的所有粒子的速度分量vi(t+1),其中下标2的含义表示对象的更新是在细胞膜m2的区间内。 进一步地，图4中的细胞膜m4用于从所述环状路线中预设配送中心点所在位置出发，使用深度优先搜索获取所述环状路线所在的位置，并根据所述位置获取从所述待加入路线集中选取的路线的节点；并获取所述节点在所述s条放射状配送路线中对应的位置信息，采用回溯法获取所述环状路线中可删除路线的集合。 在具体实现中，可以从所述s条放射状配送路线中随机选取一条配送路线，从被选取的配送路线对应的待加入路线集中随机选取一条路线加入所述被选取的配送路线，获取包含所述被选取的配送路线的环状路线。 步骤s204：根据所述环状路线，获取所述环状路线中可删除路线的集合； 本实施例根据所述环状路线，从所述环状路线中预设配送中心点所在位置出发，使用深度优先搜索获取所述环状路线所在的位置，根据所述位置获取从所述待加入路线集中选取的路线的节点；并获取所述节点在所述s条放射状配送路线中对应的位置信息；根据所述位置信息，采用回溯法获取所述环状路线中可删除路线的集合； 步骤s205：从所述集合中选取一条可删除路线，并将其在所述环状路线中对应的路线删除，获得被选取的配送路线的可行路线； 结合图4，在具体实现中，细胞膜m3用于接收从细胞膜m1和细胞膜m2传递过来的信息并更新当前的路径，根据从细胞膜m2传过来的待加入路线集，细胞膜m3中的规则开始执行相应的进化操作，其具体的实施方法如下： 首先，细胞膜m3随机选取从细胞膜m2传递进来的待加入路线集中的边加入到当前放射状网络中，此时形成环状线路，然后细胞膜m3中的转运规则此线路发送到细胞膜m4中，通过细胞膜m4中规则的执行再反馈到细胞膜m3，最后细胞膜m3依次执行切割操作，最终得到可行的线路，即所求得的最佳线路，当前线路中的配送中心点就是最佳中心位置。根据上述叙述，细胞膜m3的进化规则机制r2如公式(5)所示：其中xi表示当前线路，集合u表示边eik的可删除路线集。 需要说明的是，所述公式(5)表示：对象xi在接收细胞膜m4传递过来的可删除边集u和xik(xik＝1)之后，将边eij加入到当前放射状网络中，即xij从xij＝1变为xij’＝0；然后细胞膜m3从集合u中随机选取一条边eij(xij＝0)删除，此时xij＝0变为xij’＝1(xij’∈u)表示此边eij从放射状配送线路中删除，其中，在图论中，xij＝0表示此边eij在放射状配送路线中,xij＝1表示边eij不在当前线路中，xi’表示新产生的对象。 此外，细胞膜m0，用于储存函数值最小的配送线路，且细胞膜m0与细胞膜m1之间存在一定的膜信息交换，细胞膜m0接收细胞膜m1转运过来的局部因子，并更新细胞膜m0中的对象。 步骤s206：从剩余放射状配送路线中选取配送路线并获取对应的可行路线，直至获取到s条放射状配送路线对应的可行路线； 步骤s207：根据获取到的s条放射状配送路线对应的可行路线，得到所述配送中心点对应的可行路线集； 步骤s208：选取剩余的配送中心点，并获取对应的可行路线集，直至获取各预设配送中心点对应的可行路线集。 需要说明的是，本实施例从所述可删除路线集合中，选择可删除路线时，可通过轮盘赌的方式选取一条边进行删除，从而到一条新的配送路线，依照此种方式，直至获取被选取的配送路线的所有可行路线。在获取到获取被选取的配送路线的所有可行路线后，将所有的可行路线加入集合，即获得了先前被选取到的配送中心点对应的可行路线集。为了获取所有预设配送中心点对应的可行路线集，需要对预设配送中心点进行遍历，获取所有预设配送中心点对应的可行路线集。 参考图5，图5为本发明物流配送中心选址方法第二实施例实际应用流程图。基于上述各实施例，本实施例结合实际情况，对上述各实施例的物流配送中心选址过程进行详细描述。 根据实际的物流交通网络布局的特点，本实施例对该网络布局中所采用到的待配送客户点(以下简称客户点)的信息以及对网络布局问题的约束进行了设置与分析； 步骤s400：获取实际问题的参数集； 参照图6，图6为本发明物流配送中心选址方法第二实施例实际应用中放射状配送路线示意图； 在以下步骤的执行过程中，物流配送的中心选址都是基于以参数集作为研究对象，建立相应的优化机理模型，从而为后面的优化提供强大的数据和理论支持。 步骤s401：以参数集作为研究对象，获取带辐射状编码成位串。 获取配送中心点的配送路线的表示的形式，并将其配送线路的表达形式构造简述为： 基于配送中心的交通网络布局问题应用二进制编码表达，对于实际问题中的多重设计方案，其群体的表达形式为：p＝{x0,x1,...,xs-1}，其中每个个体xi的表示形式为：xi＝{xi0,xi1,...,xi(n-1)}(xij∈{0,1}且|xi|＝|n-m+1|)，若xij＝1，则表示顶点为vi和vj构成的边eij在当前子树中。 满足约束条件的个体x2＝{1,0,1,0,1,1,1,0,1,1,0,0,0,0,1,1,0,1,1,1}，其中，在放射状配送路线中的边eij(xij＝0)，未在放射状配送路线中的边eij(xij＝1)。对于xij＝0而言，e(2,1)(v1,v4),e(2,3)(v1,v5),e(2,7)(v2,v7),e(2,10)(v3,v4),e(2,11)(v3,v6),e(2,12)(v4,v9),e(2,14)(v5,v8)以及e(2,16)(v8,v7)构成的线路与所有客户点构成的网络即为一个带约束的放射状网络图，其中，e(2,14)(v5,v8)表示第二个对象的第14个分量的边的两个顶点分别为v5和v8。 需要说明的是，本实施例中的其他形式的编码与其编码方式在约束条件下的构造类似。 步骤s402：基于并行机制的膜结构获取在第i次迭代后的网络布局的形状的实际值。 参照图4，根据对膜结构的设计与分析可知，膜系统一共包含5个执行不同操作的膜(细胞膜m0、细胞膜m2、细胞膜m2、细胞膜m3和细胞膜m4)，其中每个膜都包含对象和规则，对象代表一个可行解，规则用于制约对象的进化范围。用于接收从其它膜传递过来的信息。对于进化膜中的对象，计算出所有对应的配送线路以及配送中心的编号；然后对获取对象，记录对用的函数实际值f。根据优化领域的先验知识，基于配送中心的配送路线的优劣变化随其函数值的变化呈正先关趋势，若f越小，则对应的配送线路就越优以及配送中心位置的选取就越佳，且这种变化趋势开始大后来小，在数学上表示为下降的曲线。 根据对膜的执行操作与分析可知，基于生物细胞机制原理的设计的膜与膜之间存在一定的通信，如图4中的箭头指向部分。对于膜规则的转运过程，基于上述模型假设，其转运机制设计如下： 在本实施例中，基于膜计算装置系统的设置，细胞膜m1和细胞膜m4之间存在双向通信机制，即细胞膜m1首先向细胞膜m4发送多重对象集ω(t)，细胞膜m4开始进化通过接收细胞膜m3和细胞膜m2传过来的信息，其规则r3的进化过程如公式(6)所示： 公式(4)中，“≡”表示隶属关系，细胞膜m1中存储的物质为对象集ω(t)，“go”表示细胞膜m1中的对象集发送到的指定的细胞中，基底细胞膜m4表示细胞膜m4接收细胞膜m1传递过来的对象集。 在本实施例中，基于粒子机制的进化机理，细胞膜m4中对象集的更新需要接收从细胞膜m3传递过来的粒子的速度集vi’(t)，基于此，其转运规则r4的设计如公式(7)所示： 其中，细胞膜m3括号中的对象为速度集vi’(t)，“go”表示将细胞膜m3中的vi’(t)发送到的指定的膜中，基底细胞膜m4表示细胞膜m4接收细胞膜m3传递过来的vi’(t)。 在本实施例中，基于优化领域求得的是最优解，细胞膜m0中存储的是至当代函数值最小的配送线路，基于此，细胞膜m0与细胞膜m1之间存在一定的膜信息交换，细胞膜m0接收细胞膜m1转运过来的局部因子pbest，然后并更新细胞膜m0中的对象pg，其相应的转运规则r5如公式(8)所示： 公式(8)中，pbest(t)表示第t代的全局最优路径，基底细胞膜m0表示接收细胞膜m1转运的对象膜。 在本实施例中，基于粒子速度集进化的数学模型，细胞膜m3的进化需要接收细胞膜m0中的对象，基于此，细胞膜m0和细胞膜m3之间存在通信机制，其转运规则r6如公式(9)所示： 公式(9)为所述带通信机制的计算规则，此规则不断向细胞膜m3中传输全局因子从而保证这个系统的顺利执行。 在本实施例中，基于换树机制的模型机理，细胞膜m3中对象的更新需要细胞膜m2中传递过来的边集才能保证细胞膜m3中的对象既满足放射状约束的要求又满足分支长度约束的要求，细胞膜m2规则的机理是用于获取配送路线中出现的环的可删除路线集，根据上述分析，其细胞膜m2和细胞膜m4的转运规则r7如公式(10)： 公式(10)中，u(xi)表示环路中的可去线路的集合，基底(xi,xij＝1)表示将细胞膜m4中的可行对象和一条余边发送到细胞膜m2中，细胞膜m2通过执行换树算子，计算出集合u，然后返回细胞膜m4，则细胞膜m4即可得到满足约束条件的物流配送线路集从而也获取配送中心的最佳位置。 所述公式(10)具体实现过程中，细胞膜m4接收从细胞膜m3传递过来的一个放射状配送路线xi和一个不在xi中的边eij(xij＝1)，然后将边eij加入到放射状路线中，此时xij从1变为0，则配送路线xi产生一条环路，然后通过执行换树机制，产生可删除边集u，然后通过细胞膜m2和细胞膜m4的通信机制将细胞膜m2的信息发送到细胞膜m4中，从而保证细胞膜m4能够产生新的对象。 步骤s403：利用膜粒子群算法迭代计算使得所述函数值最小，获得最佳位置。 基于步骤s200-步骤s220，利用膜粒子群算法迭代计算： 获得带约束的多重对象集； 其中，vi(t+1)为下一格局的函数实际值，wt(t)为影响下一格局变化的惯性权系数，c1和c2为对象向全局因子和局部因子学习的程度，r1和r2为干扰程度，△pbest和△pg为对象之间的偏差程度，带参数的函数s用于计算对象的变化分量。 基于数据库中的实验数据，结合图论中二叉树的思想，通过建立速度模型与位移模型之间的函数关系，由于对象编码是离散的，无法直接构建模型之间的函数关系，采用s函数将两模型关联建立带树状(放射状)的配送线路，其中二叉树要带有长度约束的条件：树中每个分支的权值之和不能大于p0，然后不断的迭代，计算出最优线路的函数值。 具体地，首先将公式(11)中的第二个公式带入第一个公式计算出粒子在下一时刻的速度分量集vi(t+1)，然后将vi(t+1)带入第公式(11)的第三个公式中计算出待改变的粒子位置的分量集：若s(vij(t+1))≤rij，则xij从1改为0，所有的分量计算方式依此。上述操作都是利用基于膜粒子群算法，通过速度和粒子位置的更新公式，迭代求解上式的最优解。解出的最优解，就确定了最佳的配送中心的位置。 但是，此时的数据很少，得到的配送次序的变化和函数实际值的变化很少。为了进一步得到精确解，可以多做几次实验，补充实验数据，从而得到相对精确的最优解。 例如，第一次实验产生一组数据(p(1)g,f1),第二到三十次实验产生29组实验(p(2)g,f2)，(p(3)g,f3)...(p(30)g,f30)。f1，f2...f30是有实际的数据，p(1)g,p(2)g...p(30)g只表示配送路线的形状的二进制编码没有具体的数值。此时，通过公式(11)不断的更新配送路线的编码。基于以上的方法，再次计算最优解，不断修改参数以及更新编码，得到最优解。 在本实施例中，从微观建模，基于配送中心的交通网络布局问题的数学模型，以基于配送中心的配送路线优化为研究对象，建立评估对象优劣的函数表达式；设计实验方案，通过带约束的放射状网络的布局配送模型参数，从而为选取较优的配送中心位置提高强大的数据和理论支持。 步骤s404：基于参数函数及函数的实际值评价对象的优劣。 获取最优配送线路中的配送中心的位置。 对于配送中心的选址过程，基于上述模型假设，根据步骤s200-步骤s230求得的最佳配送路线集，选取评价函数计算出路线集中的最佳配送方案，通过放射状结构的“中序遍历”方式计算出配送方案中的配送中心的结点编号，记为v1，则v1就作为系统最终格局的输出。 本实施例通过引入具有分布性和并行性计算模式的膜结构获取最佳带辐射状和长度约束要求的路径，从而能够在处理大规模算例时，高效地计算出使物流运输的总成本达到最小的最佳的配送中心位置。 此外，本发明实施例还提出一种计算机可读存储介质，所述计算机可读存储介质上存储有物流配送中心选址程序，所述物流配送中心选址程序被处理器执行时实现如下操作： 进一步地，所述物流配送中心选址程序被处理器执行时还实现如下操作： 本实施例的有益效果是：通过获取待选址交通网络中大量的物流配送点，根据所述物流配送点中的预设配送中心点，获取各预设配送中心点对应的无向图，根据所述无向图获取各预设配送中心点对应的可行路线集，并对各预设配送中心点对应的可行路线集进行筛选，将满足预设条件的可行配送路线作为目标路线，并根据所述目标路线获取对应的物流配送中心点地址，从而解决了难以对大规模物流网络中选取配送中心地址的问题进行求解的问题，提高了物流配送中心选址的效率 需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者系统不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者系统所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括该要素的过程、方法、物品或者系统中还存在另外的相同要素。 上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。 通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到上述实施例方法可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质(如rom/ram、磁碟、光盘)中，包括若干指令用以使得一台终端设备(可以是手机，计算机，服务器，空调器，或者网络设备等)执行本发明各个实施例所述的方法。