【物流选址】基于matlab麻雀搜索算法求解物流选址问题【含Matlab源码 H003期】

💥💥💥💥💥💥💞💞💞💞💞💞💞💞欢迎来到海神之光博客之家💞💞💞💞💞💞💞💞💥💥💥💥💥💥

✅博主简介：热爱科研的Matlab仿真开发者，修心和技术同步精进；
🍎个人主页：海神之光
🏆代码获取方式：
海神之光Matlab王者学习之路—代码获取方式

⛳️座右铭：行百里者，半于九十。
更多Matlab路径规划仿真内容点击👇
①Matlab路径规划（进阶版）
②付费专栏Matlab路径规划（初级版）

⛳️关注优快云海神之光，更多资源等你来！！

⛄一、物流选址简介

1 LRP问题本质
LRP问题的本质是为了满足区域类客户的配送需求，需要建立一个多级配送网络，一般由配送枢纽、配送中心、配送末端网点三级构成（其他类型可类比）。要求满足客户的需求量，不同设施有其吞吐容量和服务能力约束，不同级设施之间配送单位费用均不同，客户也有配送的时间约束。在满足多种约束下，求整个配送网络总成本最小的选址和需求分配方案。

2 LRP重难点
LRP问题求解的重难点一般是两个：选址和需求分配。选址的难度是在不同级设施备选点中选择合适数量和位置的设施组合；需求分配的难度是在不同分配路径方案中选择不同路径方案的组合，其中，因为每个客户的需求量从源头流经多级配送设施的流量是未知的，而且这个未知的变量是MGH个的。（M,G,H分别是1，2，3级配送设施的备选点数量），流量取值是连续的变量。因此，通过转化为多变量的目标函数求最值的线性或者非线性函数求解是不可能做到的。

3 LRP思路
因此，如何克服以上难点成为了我们的关键问题。首先，选址的找到一个最有的排列组合，可以通过万能钥匙的遗传算法求解。然后需求分配则需要按照最小费用最大流问题的思路进行分配。最小费用最大流问题指的是，在一个由一个源和汇和中间节点构成的配送网络中，每个路径都有其成本权重和流量约束，在满足流量约束的条件下找到最小配送成本的配送路径，这种方法就称为最小费用最大流方法。
借鉴了这种方法之后，再根据模型的约束进行改进，得到一个新的最小费用最大流方法。

⛄二、粒子群算法简介

粒子群优化(PSO)是一种基于群体智能的数值优化算法，由社会心理学家James Kennedy和电气工程师Russell Eberhart于1995年提出。自PSO诞生以来，它在许多方面都得到了改进，这一部分将介绍基本的粒子群优化算法原理和过程。
1.1 粒子群优化
粒子群优化(PSO)是一种群智能算法，其灵感来自于鸟类的群集或鱼群学习，用于解决许多科学和工程领域中出现的非线性、非凸性或组合优化问题。

1.1.1 算法思想
许多鸟类都是群居性的，并由各种原因形成不同的鸟群。鸟群可能大小不同，出现在不同的季节，甚至可能由群体中可以很好合作的不同物种组成。更多的眼睛和耳朵意味着有更多的及时发现食物和捕食者的机会。鸟群在许多方面对其成员的生存总是有益的：
觅食：社会生物学家E.O.Wilson说，至少在理论上，群体中的个体成员可以从其他成员在寻找食物过程中的发现和先前的经验中获益[1]。如果一群鸟的食物来源是相同的，那么某些种类的鸟就会以一种非竞争的方式聚集在一起。这样，更多的鸟类就能利用其他鸟类对食物位置的发现。
抵御捕食者：鸟群在保护自己免受捕食者侵害方面有很多优势。
更多的耳朵和眼睛意味着更多的机会发现捕食者或任何其他潜在的危险；
一群鸟可能会通过围攻或敏捷的飞行来迷惑或压制捕食者；
在群体中，互相间的警告可以减少任何一只鸟的危险。
空气动力学：当鸟类成群飞行时，它们经常把自己排成特定的形状或队形。鸟群中鸟的数量不同，每只鸟煽动翅膀时产生不同的气流，这都会导致变化的风型，这些队形会充分利用不同的分型，从而使得飞行中的鸟类能够以最节能的方式利用周围的空气。
粒子群算法的发展需要模拟鸟群的一些优点，然而，为了了解群体智能和粒子群优化的一个重要性质，值得提一下是鸟群的一些缺点。当鸟类成群结队时，也会给它们带来一些风险。更多的耳朵和眼睛意味着更多的翅膀和嘴，这导致更多的噪音和运动。在这种情况下，更多的捕食者可以定位鸟群，对鸟类造成持续的威胁。一个更大的群体也会需要更多的食物，这导致更多食物竞争，从而可能淘汰群体中一些较弱的鸟类。这里需要指出的是，PSO并没有模拟鸟类群体行为的缺点，因此，在搜索过程中不允许杀死任何个体，而在遗传算法中，一些较弱的个体会消亡。在PSO中，所有的个体都将存活，并在整个搜索过程中努力让自己变得更强大。在粒子群算法中，潜在解的改进是合作的结果，而在进化算法中则是因为竞争。这个概念使得群体智能不同于进化算法。简而言之，在进化算法中，每一次迭代都有一个新的种群进化，而在群智能算法中，每一代都有个体使自己变得更好。个体的身份不会随着迭代而改变。Mataric[2]给出了以下鸟群规则：

安全漫游：鸟类飞行时，不存在相互间或与障碍物间的碰撞；
分散：每只鸟都会与其他鸟保持一个最小的距离；
聚合：每只鸟也会与其他鸟保持一个最大的距离；
归巢：所有的鸟类都有可能找到食物来源或巢穴。
在设计粒子群算法时，并没有采用这四种规则来模拟鸟类的群体行为。在Kennedy和Eberhart开发的基本粒子群优化模型中，对agent的运动不遵循安全漫游和分散规则。换句话说，在基本粒子群优化算法的运动过程中，允许粒子群优化算法中的代理尽可能地靠近彼此。而聚合和归巢在粒子群优化模型中是有效的。在粒子群算法中，代理必须在特定的区域内飞行，以便与任何其他代理保持最大距离。这就相当于在整个过程中，搜索始终停留在搜索空间的边界内或边界处。第四个规则，归巢意味着组中的任何代理都可以达到全局最优。

在PSO模型的发展过程中，Kennedy和Eberhart提出了五个判断一组代理是否是群体的基本原则：

就近原则：代理群体应该能够进行简单的空间和时间计算；
质量原则：代理群体能够对环境中的质量因素作出反应；
多响应原则：代理群体不应在过于狭窄的通道从事活动；
稳定性原则：代理群体不能每次环境变化时就改变其行为模式；
适应性原则：计算代价不大时，代理群体可以改变其行为模式。

1.1.2 粒子群优化过程
考虑到这五个原则，Kennedy和Eberhart开发了一个用于函数优化的PSO模型。在粒子群算法中，采用随机搜索的方法，利用群体智能进行求解。换句话说，粒子群算法是一种群智能搜索算法。这个搜索是由一组随机生成的可能解来完成的。这种可能解的集合称为群，每个可能解都称为粒子。
在粒子群优化算法中，粒子的搜索受到两种学习方式的影响。每一个粒子都在向其他粒子学习，同时也在运动过程中学习自己的经验。向他人学习可以称为社会学习，而从自身经验中学习可以称为认知学习。由于社会学习的结果，粒子在它的记忆中存储了群中所有粒子访问的最佳解，我们称之为gbest。通过认知学习，粒子在它的记忆中储存了迄今为止它自己访问过的最佳解，称为pbest。

任何粒子的方向和大小的变化都是由一个叫做速度的因素决定的，速度是位置相对于时间的变化率。对于PSO，迭代的是时间。这样，对于粒子群算法，速度可以定义为位置相对于迭代的变化率。由于迭代计数器单位增加，速度v的维数与位置x相同。

对于D维搜索空间，在时间步t下群体中的第ith个粒子由D维向量x i t = ( x i 1 t , ⋯   , x i D t ) T x_i^t = {(x_{i1}^t, \cdots ,x_{iD}^t)T}xit​=(xi1t​,⋯,xiDt​)T来表示，其速度由另一个D维向量v i t = ( v i 1 t , ⋯   , v i D t ) T v_i^t = {(v_{i1}^t, \cdots ,v_{iD}^t)T}vit​=(vi1t​,⋯,viDt​)T表示。第ith个粒子访问过的最优解位置用p i t = ( p i 1 t , ⋯   , p i D t ) T p_i^t = {\left( {p_{i1}^t, \cdots ,p_{iD}^t} \right)^T}pit​=(pi1t​,⋯,piDt​)T表示，群体中最优粒子的索引为“g”。第ith个粒子的速度和位置分别由下式进行更新：
v i d t + 1 = v i d t + c 1 r 1 ( p i d t − x i d t ) + c 2 r 2 ( p g d t − x i d t ) (1) v_{id}^{t + 1} = v_{id}^t + {c_1}{r_1}\left( {p_{id}^t - x_{id}^t} \right) + {c_2}{r_2}\left( {p_{gd}^t - x_{id}^t} \right)\tag 1vidt+1​=vidt​+c1​r1​(pidt​−xidt​)+c2​r2​(pgdt​−xidt​)(1)

x i d t + 1 = x i d t + v i d t + 1 (2) x_{id}^{t + 1} = x_{id}^t + v_{id}^{t + 1}\tag 2xidt+1​=xidt​+vidt+1​(2)

其中d=1,2,…,D为维度，i=1,2,…,S为粒子索引，S是群体大小。c1和c2为常数，分别称为认知和社交缩放参数，或简单地称为加速系数。r1和r2是满足均匀分布[0,1]之间的随机数。上面两个式子均是对每个粒子的每个维度进行单独更新，问题空间中不同维度之间唯一的联系是通过目标函数引入的，也就是当前所找到的最好位置gbest和pbest[3]。PSO的算法流程如下：
1.1.3 解读更新等式
速度更新等式（1）的右侧包括三部分3：
前一时间的速度v，可以认为是一动量项，用于存储之前的运动方向，其目的是防止粒子剧烈地改变方向。
第二项是认知或自我部分，通过这一项，粒子的当前位置会向其自己的最好位置移动，这样在整个搜索过程中，粒子会记住自己的最佳位置，从而避免自己四处游荡。这里需要注意的是，pidt-xidt是一个方向从xidt到pidt的向量，从而将当前位置向粒子的最佳位置吸引，两者的顺序不能改变，否则当前位置会远离最佳位置。
第三项是社交部分，负责通过群体共享信息。通过该项，粒子向群体中最优的个体移动，即每个个体向群体中的其他个体学习。同样两者应该是pgbt-xidt。
可以看出，认知尺度参数c1调节的是粒子在其最佳位置方向上的最大步长，而社交尺度参数c2调节的是全局最优粒子方向上的最大步长。图2给出了粒子在二维空间中运动的典型几何图形。

图2 粒子群优化过程中粒子移动的几何说明

从更新方程可以看出，Kennedy和Eberhart的PSO设计遵循了PSO的五个基本原则。在粒子群优化过程中，在d维空间中对一系列时间步进行计算。在任何时间步，种群都遵循gbest和pbest的指导方向，即种群对质量因素作出反应，从而遵循质量原则。由于速度更新方程中有均布随机数r1和r2，在pbest和gbest之间的当前位置随机分配，这证明了响应原理的多样性。在粒子群优化过程中，只有当粒子群从gbest中接收到较好的信息时，才会发生随机运动，从而证明了粒子群优化过程的稳定性原则。种群在gbest变化时发生变化，因此遵循适应性原则。
1.2 粒子群优化中的参数
任何基于种群的算法的收敛速度和寻优能力都受其参数选择的影响。通常，由于这些算法的参数高度依赖于问题参数，因此不可能对这些算法的参数设置给出一般性的建议。但是，已有的理论和/或实验研究，给出了参数值的一般范围。与其他基于种群的搜索算法类似，由于在搜索过程中存在随机因素r1和r2，因此通用PSO的参数调整一直是一项具有挑战性的任务。PSO的基础版本只需要很少的参数。本章只讨论了[4]中介绍的PSO基础版本的参数。

一个基本的参数是群体规模，它通常是根据问题中决策变量的数量和问题的复杂性经验地设置的。一般建议20-50个粒子。

另一个参数是缩放因子c1和c2。如前所述，这些参数决定了下一个迭代中粒子的步长。也就是说，c1和c2决定了粒子的速度。在PSO的基础版本中，选择c1=c2=2。在这种情况下，粒子s速度的增加是不受控制的，这有利于更快的收敛速度，但不利于更好地利用搜索空间。如果我们令c1=c2>0，那么粒子会吸引到pbest和gbest的平均值。c1>c2设置有利于多模态问题，而c2>c1有利于单模态问题。在搜索过程中，c1和c2的值越小，粒子轨迹越平滑，而c1和c2的值越大，粒子运动越剧烈，加速度越大。研究人员也提出了自适应加速系数[5]。
停止准则不仅是粒子群算法的参数，也是任何基于种群的元启发式算法的参数。常用的停止准则通常基于函数评估或迭代的最大次数，该次数与算法所花费的时间成正比。一个更有效的停止准则是基于算法的搜索能力，如果一个算法在一定的迭代次数内没有显著地改进解，那么应该停止搜索。

⛄三、部分源代码

%==========================================================================
clear all;
%这里是参数的预备部分，调试时就调整这个里面的参数即可preparation
itmax=400;%进化代数，就是预设的迭代次数。
W(1)=0.729;% 粒子先前速度保持。惯性权重
a(1)=0.316;% 用于计算W。
c1=2; %1.45449; %认知部分 加速系数
c2=2; %1.45449;% 社会部分 加速系数
m=3; % 待选择DC的个数（一定要知道，这是基数，不是序数，即不是需求点的标号）
M=[13 13 13];%各DC的容量均为13个单位，M是m维向量，此处就是3维了。%要×1000的。
R=1000000000; %罚系数，要求充分大即可
h=[ 0 1 6 7 4 3 4 6 6 9 8 9;
1 0 5 6 5 4 5 7 7 10 9 10;
6 5 0 3 6 9 10 12 12 15 14 15;
7 6 3 0 3 10 11 13 13 16 15 12;
4 5 6 3 0 7 8 10 10 13 12 9;
3 4 9 10 7 0 6 4 9 10 6 6;
4 5 10 11 8 6 0 2 9 5 4 9;
6 7 12 13 10 4 2 0 10 6 2 7;
6 7 12 13 10 9 9 10 0 4 8 13;
9 10 15 16 13 10 5 6 4 0 4 9;
8 9 14 15 12 6 4 2 8 4 0 5;
9 10 15 12 9 6 9 7 13 9 5 0]; % h必须是对称矩阵。
%h矩阵表示单位运输量的成本，第（i,j）个元素指从第i个需求点到第j个需求点的单位运输量的成本。
dem=[5 4 2 3 2 4 3 5 4 3 2 2];% dem 原来是demand ，需求量
F=[11 16 14 14 15 13 18 12 11 14 16 11]; %固定费用
num_particle=100; % number of particle, should be more than
D=length(h); % 需求点个数
%**********************************************************
%下面生成随机初始值x,假设初始的服从[1,3]均匀分布,即[1,m].
%初始v 在[-2,2]之间的整数，即[-(m-1),m-1].
xmax=m;
xmin=1;
IndexIteration=1;
tic;
particle=zeros(2num_particle,D,itmax);
particle(:,:,1)=xmin+(xmax-xmin)rand(2num_particle,D);
V(:,:,1)=round((xmin-xmax)+2(xmax-xmin)rand(2num_particle,D));
%**********************************************************
fit=zeros(num_particle,itmax);% 用于存储粒子的适应值
pbest=zeros(2num_particle,D,itmax); % 用于存储粒子的位置
num=numbering(m,D,2num_particle,particle(:,:,1));% 粒子的编号及对应关系
fit(:,1)=fitness(h,F,m,D,2num_particle,dem,num,M,R);
%********************************************************
pbest(:,:,1)=particle(:,:,1);
pbest_value(:,1)=fit(:,1); %个体最优值
[C,I]=min(pbest_value(:,1));
gbest_value(1)=C; % 群最优值
for i=1:num_particle
gbest(2i-1:2i,:,1)=particle(2I-1:2I,:,1); % 群最优粒子位置
end
%第一次迭代中最优就是本身的结果，因为就一个嘛，没有选择
%******************************************************
for IndexIteration=2:itmax
%V(:,:,IndexIteration)=W(IndexIteration-1)V(:,:,IndexIteration-1)+c1rand*(pbest(:,:,IndexIteration-1)…
% -particle(:,:,IndexIteration-1))+…
%c2rand(gbest(:,:,IndexIteration-1)-particle(:,:,IndexIteration-1));

V(:,:,IndexIteration)=0.729V(:,:,IndexIteration-1)+c1rand*(pbest(:,:,IndexIteration-1)…
-particle(:,:,IndexIteration-1))+…
c2rand(gbest(:,:,IndexIteration-1)-particle(:,:,IndexIteration-1));

V(:,:,IndexIteration)=min(V(:,:,IndexIteration),xmax-xmin);%
V(:,:,IndexIteration)=max(V(:,:,IndexIteration),xmin-xmax);% 所有粒子的速度均在[-(m-1),m-1]内
particle(:,:,IndexIteration)=particle(:,:,IndexIteration-1)+V(:,:,IndexIteration);
particle(:,:,IndexIteration)=min(particle(:,:,IndexIteration),xmax);%
particle(:,:,IndexIteration)=max(particle(:,:,IndexIteration),xmin);% 所有粒子均在[1,m]内

%for i=1:num_particle
%V(2i,:,IndexIteration)=min(V(2i,:,IndexIteration),xmax-xmin);%
%V(2i,:,IndexIteration)=max(V(2i,:,IndexIteration),xmin-xmax);% 所有粒子的速度均在[-(m-1),m-1]内
%end
%particle(:,:,IndexIteration)=particle(:,:,IndexIteration-1)+V(:,:,IndexIteration);
%for i=1:num_particle
%particle(2i,:,IndexIteration)=min(particle(2i,:,IndexIteration),xmax);%
%particle(2i,:,IndexIteration)=max(particle(2i,:,IndexIteration),xmin);%
%所有粒子均在[1,m]内
num=numbering(m,D,2num_particle,particle(:,:,IndexIteration));% 粒子的编号及对应关系
fit(:,IndexIteration)=fitness(h,F,m,D,2num_particle,dem,num,M,R);
%下面更新 pbest and pbest_value
pbest_value(:,IndexIteration)=min(pbest_value(:,IndexIteration-1),fit(:,IndexIteration));
for i=1:num_particle
if pbest_value(i,IndexIteration)==fit(i,IndexIteration)
pbest(2i-1:2i,:,IndexIteration)=particle(2i-1:2i,:,IndexIteration);
else
pbest(2i-1:2i,:,IndexIteration)=pbest(2i-1:2i,:,IndexIteration-1);
end
end
function y=numbering(m,D,n,x);
[xs,xid]=sort(x’);
xz=floor(x’)+1;
xx=zeros(D,n);
for i=1:D
for j=1:2:n
if i<=D-m
xx(xid(i,j),j)=0;
elseif i>D-m
xx(xid(i,j),j)=i-(D-m);
end % x是按列看的，最后要转置的。
end
for j=2:2:n
if xz(i,j)>m
xx(i,j)=m;
elseif xz(i,j)<1
xx(i,j)=1;
else
xx(i,j)=xz(i,j);
end
end
end

⛄四、运行结果

⛄五、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1]阮信波,刘丽华,陈丽瑾.麻雀搜索算法在物流配送中心选址的应用[J].物流技术. 2021,40(12)

3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除

🍅 仿真咨询
1 各类智能优化算法改进及应用
生产调度、经济调度、装配线调度、充电优化、车间调度、发车优化、水库调度、三维装箱、物流选址、货位优化、公交排班优化、充电桩布局优化、车间布局优化、集装箱船配载优化、水泵组合优化、解医疗资源分配优化、设施布局优化、可视域基站和无人机选址优化

2 机器学习和深度学习方面
卷积神经网络（CNN）、LSTM、支持向量机（SVM）、最小二乘支持向量机（LSSVM）、极限学习机（ELM）、核极限学习机（KELM）、BP、RBF、宽度学习、DBN、RF、RBF、DELM、XGBOOST、TCN实现风电预测、光伏预测、电池寿命预测、辐射源识别、交通流预测、负荷预测、股价预测、PM2.5浓度预测、电池健康状态预测、水体光学参数反演、NLOS信号识别、地铁停车精准预测、变压器故障诊断

3 图像处理方面
图像识别、图像分割、图像检测、图像隐藏、图像配准、图像拼接、图像融合、图像增强、图像压缩感知

4 路径规划方面
旅行商问题（TSP）、车辆路径问题（VRP、MVRP、CVRP、VRPTW等）、无人机三维路径规划、无人机协同、无人机编队、机器人路径规划、栅格地图路径规划、多式联运运输问题、车辆协同无人机路径规划、天线线性阵列分布优化、车间布局优化

5 无人机应用方面
无人机路径规划、无人机控制、无人机编队、无人机协同、无人机任务分配

6 无线传感器定位及布局方面
传感器部署优化、通信协议优化、路由优化、目标定位优化、Dv-Hop定位优化、Leach协议优化、WSN覆盖优化、组播优化、RSSI定位优化

7 信号处理方面
信号识别、信号加密、信号去噪、信号增强、雷达信号处理、信号水印嵌入提取、肌电信号、脑电信号、信号配时优化

8 电力系统方面
微电网优化、无功优化、配电网重构、储能配置

9 元胞自动机方面
交通流 人群疏散 病毒扩散 晶体生长

10 雷达方面
卡尔曼滤波跟踪、航迹关联、航迹融合