1、1农产品物流配送中心 农产品物流配送中心选址数学建模作业题目：农产品物流配送中心选址队员：姓名：贺海龙姓名：邹高永姓名：姜昌海2011 年 08 月 09 日农产品物流配送中心选址摘要应用层次分析法对农产品配送中心选址问题进行研究。2结合农产品具有的区域性，季节性和分散性等特点选择主要应考虑因素，建立多要素的层次结构模型，构建各层次判断矩阵，经过层次的一致性检验确定判断矩阵的一致性是否可接受，计算出各因素重要程度的权数，而后利用熵权法计算各指标权重，为多指标综合评价提供依据。结果表明，此方法将定性信息定量化，将各要素层次化、明朗化，具有较高的应用价值。关键词：农产品配送中心选址；层次分析法；一
2、致性检验；指标权重；熵权法；一 问题的提出 农产品配送中心选址是以提高物流系统的经济效益和社会效益为目标，根据农产品供货状况、需求分布、运输条件、自然环境等因素，用系统工程的方法，对配送中心的地理位置进行决策的过程。不仅如此，物流运输费用的高低直接的影响了产品的销售价格及销售区域。因此如何降低物流运输费用，选择合理的配送地址，成为农产品经营者的首要解决的问题。配送中心选址合理与否会直接影响到配送中心各项活动的成本、作业效率、服务水平和经济效益。因此，配送中心选址与布局必须在充分调查分析的基础上综合自身经营的特点、交通状况等因素，在详细分析现状及预测的基础上进行配送中心选址。 特别对于农产品来说
3、，它不同于一3般的商品，因为农产品不仅仅要达到运费低，而且必须保证没有发生腐烂，因此对于农产品物流配送中心的选址就更加重要。不仅要减低农产品的运输成本，而且要保鲜保质，是目前需要解决的重要问题。二 评价指标构造2.1 初选指标根据主要的分为 V（交通状况，经济效益，自然环境，其他因素） 交通状况：交通条件是影响物流效率和配送成本的重要因素，特别是大宗物资的配送。物流配送中心选址应接近交通运输枢纽，尽量使配送 交通条件是影响物流效率和配送成本的重要因素，特别是大宗物资的配送。物流配送中心选址应接近交通运输枢纽，尽量使配送中心形成物流过程中的一个恰当的结点。在有条件的情况下，配送中心应尽可能靠近交
4、通要道，如高速公路、铁路货运站、港口、空港等。（1）交通运输时间（物流的效率）：农产品运输时间是物流配送中最关键的一点，能不能及时将产品运输到配送中心，或是将产品从配送中心运输到销售处是衡量物流的一个重要指标。路况越好，运输时间越短，物流反应速度也快，物流效率也就越高。这对农产品的运输来说至观重要。4（2）交通可能度（交通的拥挤程度）：运输可能度反映了交通拥挤程度，交通越便利，单位时间内运送的货物越多，物流的运输量也就越多。（3）公共设施的好坏：配送中心的所在地，要求城市的道路、通信等公共设施齐备，有充足的供电、水、热、燃气的能力，且场区周围要有污水、固体废弃物等的处理能力。对于水电等资源，如
5、果它是设施大量需要的，在选址过程中，不仅要考虑其供应量是否充足，而且还要考虑其费用是否可以承受。（4）交通的便利性（靠近港口，码头不同的运输方式）：配送中心必须具备方便的交通运输条件。最好靠近交通枢纽进行布局，如紧临港口、交通主干道枢纽、铁路编组站或机场，有两种以上运输方式相连接，以备紧急预案。经济效益：（1）地价的选择（不同的地段会有不同的地价）：（2）发展潜力（当地经济的未来发展状况）：主要是指当地农产品发展状况，依托环境区域发展的预测还有当地经济水平及广大消费人群对农产品的需求度的发展预测，从而来得到当地发展潜力，发展潜力越大，就越应该在当地建立农产品配送中心。（3）劳动力条件：随着物流
6、配送中心的建设，许多大规5模的配送中心聚集在一起，现代化农产品的运作需要机械化处理设备，拥有数量充足、素质较高的劳动力已成为影响配送中心位置选择的重要条件。（4）物流费用：这是配送中心配送的对象，如货物来源和去向的分布情况、历史和现在以及将来的预测和发展等。配送中心应该尽可能地与生产地和配送区域形成短距离优化。货物数量是随配送规模的增长而不断增长的。货物增长率越高，越是要求配送中心选址的合理性，从而减少输送过程中不必要的浪费。环境因素：（1）气象条件：配送中心选址过程中，主要考虑的气象条件有温度，风力，降水量，年平均蒸发量等指标。选址时温度的高低，一般情况下过高的温度会加快农产品的腐蚀过程。（
7、2）地形条件：配送中心应选择地势较高、地形平坦，且应具有适当的面积与外形。若选在完全平坦的地形上是最理想的，其次选择稍有坡度或起伏的地方；对于山区陡坡地区则应该完全避开。在外形上可选择长方形，不宜选择狭长或不规则形状。（3）服务水平：服务水平也是配送中心选址所需考虑的因素。由于现代农产品物流过程中能否实现准时运送是衡量配送中心服务水平的重要指标，因此，在急需农产品配6送中心选址时，应保证客户可在任何时候向配送中心提出物流需求，都能获得快速满意的服务。（4）经营环境：配送中心所在地区的优惠物流产业政策对物流企业的经济效益将产生重要影响，数量充足和素质较高的劳动力条件也是配送中心选址考虑因素之一。
8、经营不同类型商品的配送中心最好能分别布局在不同地域。如生产型配送中心的选址应与产业结构、产品结构httP:/、工业布局等紧密结合进行考虑。物流费用是配送中心选址的重要考虑因素之一。大多数配送中心选择接近物流服务需求地，例如接近大型工业、商业区，以便缩短运距、降低运费等物流费用。其他因素：（1）国土资源利用：配送中心的规划应贯彻节约用地、充分利用国土资源的原则。配送中心一般占地面积较大，周围还需留有足够的发展空间，为此地价的高低对布局规划有重要影响。此外，配送中心的布局还要兼顾区域与城市规划用地的其他要素。（2）环境保护：配送中心的选址需要考虑保护自然环境与人文环境等因素，尽可能降低对城市生活的
9、干扰。对于大型转运枢纽，应适当设置在远离市中心区的地方，使得大城市交通环境状况能够得到改善，城市的生态建设得以维持和增进。7（3）人流量:人流量的多少直接决定了经济效益的好坏。（4）治安管理：好的治安管理能给配送中心的管理带来很多方便。2.2 指标筛选一级指标有四个，二级指标有比较多，现在对二级指标进行筛选，首先构造评价矩阵为 X(xij)2016，然后对评价矩阵求协方差，记为:S(sij)1616将初选指标集分成两部分，分别记为 C(1),C(2)，则相应的方差可以写成SS11S21S12 则 C(1)的协方差阵为 S11，C(2)的协方差阵为S22， S12，S21 为 S22来自 C(1
10、),C(2)的元素的协方差阵。给定 C(1)后，C(2)对C(1)的条件协方差阵为：1S(C(2)|C(1)S22S21S11S12标记 D(i)S(C(2)|C(1)(i1,16)，然后把 D(i)从大到小排序，去掉最小的达到筛选指标的目的。2.3 最终指标一级指标有：交通状况，经济效益，自然环境，其他因素；二级指标有，交通运输时间，公共设施的好坏，交通的便利性，地价的选择，发展潜力，物流费用，气象条件，地形条件，经营环境，国土资源利用环境保护人流量。三 指标权重83.1 层次分析法首先将所要分析的问题层次化，根据问题的性质和所要达到的总目标，将问题分解为不同的组成因素，即层次每一层次里都包
11、含若干组成要素，并按因素间的相互关联影响以及隶属关系将因素按不同层次聚集组合，对每一层次中的各个要素进行两两比较，确定其相互之间的重要性，最后，综合计算各层相对的重要性权值，确定各要素相当于总目标的权重。以此作为选择方案的依据。最终 AHP 法分析的系统归结为最低层次相对于最高层次相对重要程序的权值或相对优劣次序的排序问题。根据判断矩阵标度，构造互反对称矩阵检验 A，然后求 A 得最大特征值max ，接着对 A3.1.1 一致性检验判断矩阵中的数值多半是根据专家意见而给出的，而衡量判断矩阵的合格标准是它是否具有一致性，一般只需近似的满足一致性即可。在给出判断矩阵后应进行一致性检验，方法如下：a
12、.计算一致性指标 CI。公式为 CImaxnn1 公式中，n 为判断矩阵的阶数，max 为判断矩阵的最大特征值。b. 查询平均随机一致性指标 RI(见下表)9c. 计算一致性比例 CR.公式为将初选指标集分成两部分，分别记为C(1),C(2)，则相应的方差可以写成SS11S21S12则 C(1)的协方差阵为 S11，C(2) 的协方差阵为S22， S12，S21 为 S22来自 C(1),C(2)的元素的协方差阵。给定 C(1)后，C(2)对C(1)的条件协方差阵为：1S(C(2)|C(1)S22S21S11S12标记 D(i)S(C(2)|C(1)(i1,16)，然后把 D(i)从大到小排序
13、，去掉最小的达到筛选指标的目的。2.3 最终指标一级指标有：交通状况，经济效益，自然环境，其他因素；二级指标有，交通运输时间，公共设施的好坏，交通的便利性，地价的选择，发展潜力，物流费用，气象条件，地形条件，经营环境，国土资源利用环境保护人流量。三 指标权重103.1 层次分析法首先将所要分析的问题层次化，根据问题的性质和所要达到的总目标，将问题分解为不同的组成因素，即层次每一层次里都包含若干组成要素，并按因素间的相互关联影响以及隶属关系将因素按不同层次聚集组合，对每一层次中的各个要素进行两两比较，确定其相互之间的重要性，最后，综合计算各层相对的重要性权值，确定各要素相当于总目标的权重。以此作
14、为选择方案的依据。最终 AHP 法分析的系统归结为最低层次相对于最高层次相对重要程序的权值或相对优劣次序的排序问题。根据判断矩阵标度，构造互反对称矩阵检验A，然后求 A 得最大特征值max ，接着对 A3.1.1 一致性检验判断矩阵中的数值多半是根据专家意见而给出的，而衡量判断矩阵的合格标准是它是否具有一致性，一般只需近似的满足一致性即可。在给出判断矩阵后应进行一致性检验，方法如下：a.计算一致性指标 CI。公式为 CImaxnn1公式中，n 为判断矩阵的阶数，max 为判断矩阵的最大特征值。11b. 查询平均随机一致性指标 RI(见下表)c. 计算一致性比例 CR.公式为CRCIRI理论说明
15、：当 CR 小于 0.1 时，一般认为判断矩阵的一致性可以接受；否则就需要重新进行成对比较，对判断矩阵加以调整，使之具有满意的一致性。根据公式计算 CR，因CR 小于 0.1.满足一致性指标，故判断矩阵是可行的。用 Matlab 计算得到 max，然后计算 max 对应的特征向量为 V,归一化后特征向量即相对权重为WT，利用层次分析法进行配送中心选址决策的流程图如下：3.1.2 AHP 计算权重123.2 熵权法采集物流配送地址的评价值设为 xij，i 表示第 i 个备选中心，指标在该组的序号。则物流选址配送中心评价矩阵为的具体计算步骤：步骤一：由于各指标的优劣存在差异，必须对原始数据做 xi
16、j 归一化处理，计算比重j 表示第 j 个评价X(xij)m6 。下面给出熵权法pij。pijxijxi1m0iijm ，步骤二：计算第j 项指标的熵值 ej：ejkpijlnpij j1,2,3,4,5,6i1m13其中，k0 。若 xij 对于给定的 i 全部相等，那么，则 pij1。此时，ej 取 m1lnm，则得极大值，即ejki1m11lnmm（j1,2,3,4,5,6 ） 。若设 k0ej1 。步骤三：计算第 i 项指标的偏差度 gi对于给定的j，指标在评价体系中的作用越小，越接近完全无序状态，ej 越大，指标值的偏差度应越小，因此定义 gj 步骤四：计算熵权值1ej ，j1,2,
17、3,4,5,6wj14gjgj16jj1,2,3,4,5,63.3 组合法用层次分析得到的权重带有很大的主观性，用熵权法求得权重比较客观，给主观和客观的到的权重在加权，的到综合权重WWWT1(Wj)3.4 理想解法假设理想点（系统）为*(y1,y2,ym)，对于一个被评价对象(yi1,yi2,yim) ，定义二者之间的加权距离为yiwjf(yijy*j),i1,2,n ，j1m其中 wj 为权系数，15*为与 yf(yij,y*)yi 案。然而笔者必须指出，本文仅将一般农产品配送中心选择问题应考虑的主要因素考虑在内，而有可能没有将某些其他特殊配送中心的特殊影响因素包含在内。对不同类型的配送中心
18、进行选址时对于某些影响因素的重要程度赋值也要根据实际情况做出调整，然而，由于客观事物的复杂性，特殊性和人的认识的有限性和差异性，将会对评价结果产生一定程度的不同影响。如在对农产品配送中心选址时天气因素的权重赋值应适当增大。在各种其他具体应用中，操作者只需根据配送中心的性质及作用等实际情况增删影响因素并调节其权重即可。五 分析通过层次分析法，熵权法，得到了个指标的权值，然后用理想解法求得最大正理想解和最大负理想解，通过求距离的到了理想解的贴近度 C，对贴近度从大到小排序；贴近度越大则对应得方案越好，综上本文就解决了农产品物流配送中心选址的最有方案。参考文献：1 李涛，李博，农产品物流运作模式与对
19、策 期刊论文科技信息，2009（35） 2 贾会棉，路剑，我国农产品物流发展中的问题及建议J，中国经贸导刊，2006（2）：403 蒋长兵，现代物流学导论 M，北京：中国物资出版16社，2006 4 郑畅，物流中心选址方法研究D，武汉：武汉理工大学，20045 冯耕中，物流配送中心规划与设计M，西安：西安交通大学出版社，2006六 附件6.1 算法程序clear all clcA=rand(20,16);%随机产生 20*16 的矩阵，二级指标有 16 个 b(:,:,16)=zeros(20,16); % 交换列 for j=1:16 B=A(:,17-j); A(:,17-j)=A(:,1
20、6); A(:,16)=B; b(:,:,j)=A; end b;D(:,:,16)=zeros(1); for i=1:16s=cov(b(:,:,i); %计算协方差 s11=s(1:15,1:15); s22=s(16,16); s12=s(1:15,16); s21=s(16,1:15); ss=s11 s12;s21 s22;D(:,:,i)=s22-s21*inv(s11)*s12; end D;%miu=0.012;%C=Dmiu; %筛选指标%每个二级指标都有 4 个指标，通过排序去掉最小的达到筛选目的17AA=D(:,:,1),D(:,:,2),D(:,:,3),D(:,:,
21、4),D(:,:,5),D(:,:,6),D(:,:,7),D(:,:,8),D(:,:,9),D(:,:,10),D(:,:,11),D(:,:,12),D(:,:,13),D(:,:,14),D(:,:,15),D(:,:,16);dd1=AA(1,1:4);dd2=AA(1,5:8);dd3=AA(1,9:12);dd4=AA(1,13:16);b1,index=sort(dd1); %对每一个排序，去掉最小的达到筛选目的b2,index=sort(dd2);b3,index=sort(dd3);b4,index=sort(dd4);%AHP%判断矩阵 AA=1 3 5;1/3 1 2;
22、1/5 1/2 1;V,D=eig(A);%lmd=max(D)%lmd2=max(lmd)18lamd=max(eig(A); %最大特征值%RI=0 0 0.58 0.9 1.12 1.24 1.32 1.41 1.45 1.49;% 一致性检验CR=(lamd-8)/7)/0.58;w1=V(:,1); % 最大特征值对于的特征向量 （AHP 权重）for i=1:3w1=w1./sum(w1); %归一化end%熵权法%m=input(最终确定的指标个数 m=) 4%n=input(物流备选地址的个数 n=)r=rand(6,12);% 随机产生一个评价矩阵 rand(n,m)for
23、i=1:12 %mR(:,i)=r(:,i)./sum(r(:,i);end % 对矩阵进行归一化处理for j=1:12%me(j)=-sum(R(:,j).*log(R(:,j)./(log(6);%log(n)19end %计算熵值R;e;g=1-e ;% 计算偏差度w2=g./sum(g); % 计算熵权值 每个二级指标分支的权重w21=w2(1,1:3)w22=w2(1,4:6)w23=w2(1,7:9)w24=w2(1,10:12)%组合权重% a=input(请输入组合权重 a=)a=0.4;% ww(:,:,4)=zeros(1,3);ww1=(a*w1+(1-a)*w21);
24、 %主观（AHP） ，客观（熵权法）综合权重ww2=(a*w1+(1-a)*w22);ww3=(a*w1+(1-a)*w23);ww4=(a*w1+(1-a)*w24);20%二级指标下的方案选择rr=r;rr1=rr(1:3,1:6);rr2=rr(4:6,1:6);rr3=rr(7:9,1:6);rr4=rr(10:12,1:6);%在指标下各方案的权qz1=ww1*rr1;qz2=ww2*rr2;qz3=ww3*rr3;qz4=ww4*rr4;QZ=qz1;qz2;qz3;qz4%理想解法% ww=input(m 个评价指标的综合权重向量为 w=) %此处的权重要求是综合主观、客观判断矩
25、阵后的权重ww=0.2 0.4 0.3 0.1; %一级指标的权重，主观给出for i=1:4RR(i,:) =QZ(i,:).*ww(i)end%计算加权评价矩阵21RR=RRMAX_V=max(RR)%计算加权评价矩阵各列最大值，即正理想解MIN_V=min(RR)%计算加权评价矩阵各列最小值，即负理想解for i=1:6S_MAX(i)=sum(RR(i,:)-MAX_V). );%计算各备选地址到正理想解的距离 S_MIN(i)=sum(RR(i,:)-MIN_V). );%计算各备选地址到负理想解的距离endfor i=1:6C(i)=S_MIN(i)/(S_MAX(i)+S_MIN(i);%计算各备选地址到理想解的贴近度 endC %贴近度CC index=sort(C) %贴近度 C 从小到大排序CCC=fliplr(index)百度搜索“就爱阅读”,专业资料,生活学习,尽在就爱阅读网,您的在线图书馆 22
