巨幼细胞危象

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN110097322A(43)申请公布日2019.08.06(21)申请号201910395954.XG06F17/16(2006.01)(22)申请日2019.05.14(71)申请人吉林大学地址130012吉林省长春市前进大街2699号(72)发明人于新　姜盛乾　王坤　杨兆军　陈传海　赵安然　(74)专利代理机构长春吉大专利代理有限责任公司22201代理人崔斌(51)Int.Cl.G06Q10/08(2012.01)G06Q10/06(2012.01)G06K9/62(2006.01)G06F17/15(2006.01)权利要求6页说明书12页附图2页(54)发明名称一种虚拟场景中零件库存节点选址方法(57)摘要本发明属于库存节点选址技术领域，具体的说是一种虚拟场景中零件库存节点选址方法。本发明提供一种将工厂内的加工单元零件需求量和未来加工单元零件需求量作为评价指标，构建模糊评价矩阵后，利用TOPSIS方法建立出GE矩阵计算出工厂内零件物流的差异化指标，再加上库存节点的距离、库存节点成本，通过基于聚类的重心算法求解出虚拟场景中零件库存节点位置。本发明能通过减少工厂内物流量，有效降低厂内的物流部门的工作人员和叉车数量，从而降低物流成本。工厂库存选址领域暂未有相关技术研究，本方法提供一种技术填补当前空白。CN110097322ACN110097322A权　利　要　求　书1/6页1.一种虚拟场景中零件库存节点选址方法，其特征在于，该方法包括以下步骤：步骤一、加工单元运输量评价与选择；为了对库存节点进行选址，将工厂内的加工单元零件需求量和未来加工单元零件需求量作为评价指标；加工单元零件需求量影响因素有加工单元与库存区的总物流量A1、加工零件外协比例A2、加工单元物流水平A3、零件运输速率对加工的影响A4、加工单元的设备利用率A5、加工单元的日生产单元体积A6、加工单元零件装卸率A7、加工单元排程情况A8；未来加工单元零件需求量影响因素有加工单元当前设备利用率B1、当前加工单元零件当前销售情况B2、加工单元与库存区的最短距离B3、加工单元当前所占面积B4、工厂的仓库成本B5、加工单元的运输反应时间B6、加工单元的加工停滞时间B7、加工单元的运输率B8；需要对加工单元零件需求量和未来加工单元零件需求量的影响因素进行Min-max标准化处理，得到指标矩阵A和指标矩阵B；其中，a11代表第一行第一列的数据，是加工单元零件需求量影响因素有加工单元与库存区的总物流量的第一个Min-max标准化数据；b11代表第一行第一列的数据，是未来加工单元零件需求量影响因素有加工单元当前设备利用率的第一个Min-max标准化数据；首先建立加工单元零件需求量和未来加工单元零件影响关系的网络分析法即ANP结构；建立主要因子的有向图；再经过定位分析得到加工单元零件需求量和未来加工单元零件的主因子P(P1,P2,…,Pm)，次因子P(P11,…,P1n,P21,P22,…,P2l,Pm1,…,Pmk)，其中，Pm为建工单元需求量和未来加工单元零件元素组；Pmk是第m个主因子的第k个次因子，初步构建ANP结构和准则之间的相互影响有向图；构建关系矩阵后，再建立P1,P2y…,Pm之间的直接影响矩阵Z，并将矩阵Z变换成标准矩阵Y；Z＝[zij]m×mi,j＝1,2,…,n  (1.3)其中，zij表示在主因子Pi对主因子Pj的影响程度Y＝k·Z,k＞0其中，2CN110097322A权　利　要　求　书2/6页计算最终关系矩阵，构建ANP主要因子结构图；建立总影响矩阵T，经过公式1.4转化得到Ti,i＝0,1,2,...,m；根据影响因素实际关系，选取临界值p，经过计算得到最后的关系矩阵G，接着再构建ANP准则间结构；T＝Y·(I-Y)-1  (1.4)T＝[tij]m×m  (1.5)G＝[gij]m×m  (1.6)其中Y为标准矩阵，tij由式1.4换算所得；根据ANP主因子间新结构，构建Pk主因子下Pk1,Pk2…,Pkm，k∈{1,2,...,m}之间的直接关系矩阵，得到完整ANP结构；建立ANP判断矩阵；构建决策t(t＝1,2,3,...,n)的判断矩阵；主因子P1,P2…,Pm对主因子t(jl)Pj判断矩阵为A和主因子P1的元素Pi1,Pi2…,Pm对主因子Pj元素Pjl影响的判断矩阵为计算加权矩阵，对判断矩阵Ai(jl)求归化特征向量进而得到矩阵同理求得At(jl)的权重矩阵At；再求超矩阵Wt和加权矩阵t截取Pm与P11,P12…,P1n,P21,P22…,P2l…,P(m-1)1,P(m-1)2…,P(m-1)k之间的加权矩阵E，规范化后得到Rt；3CN110097322A权　利　要　求　书3/6页t其中eij为加权矩阵E中的第i行第j列数值；再构建模糊评价矩阵；将全部评分人的规范化评价矩阵进行模糊处理，从而获得模糊评价矩阵R；其中，确定理想解和负理想解；根据TOPSIS方法得到理想解向量和负理想解向量；理想解：负理想解：其中J1为效益型指标集，J2为成本型指标集；计算贴近度C；C＝[c1]1×m；其中进而，根据公式1.11-1.14，求解第i部分需求量的贴近度和得到加工单元零件需求量和未来加工单元零件需求量的评价值，再建立“零件需求量和未来零件需求量”的GE矩阵，通过加和得到加工单元的评价值；步骤二、建立库存节点选址方法；将从厂内物流成本和各单元的实际生产情况出发，选择库存节点的距离，库存节点成本和加工单元的评价值作为工厂库存节点的评价标准；库存节点之间的连接关系分为单连通结构和直通式结构；因此构建两种模型如下：单连通模型：4CN110097322A权　利　要　求　书4/6页s.t.2≤n≤nc2≤nj≤nc式中直通式模型：s.t.2≤n≤nc2≤nj≤nc式中，12f(x)为效益型函数，f(x)为成本型函数，将w1、w2、w3设置成权重，库存节点选址模型的函数就转化为：单连通模型：5CN110097322A权　利　要　求　书5/6页直通式模型：其中：n：库存建设节点数目；Cbc：库存物流单位长度费用；CbT：库存建设数额；为第k次算得的库存节点l与库存节点m间的长度，Nj：加工单元运输量集，Nj＝1,2,...,nc且i,j∈Nj；dlj：库存节点l与加工单元j的距离，dlj＝|xl-xj|+|yl-yj|；θj：加工单元j的综合评价值；Il：库存节点l的加工单元运输量集合；L：库存节点集，L＝1,2,...,n且l,k∈L；nc：加工单元数目；dmax：库存节点最大服务范围；dmin：加工单元库存最小距离；Mmin：极小值，可取值为标准化的dmin；库存建设投资额；nl：库存节点l输送的加工单元运送量；xl：库存节点l的横坐标；yl：库存节点l的横坐标。dlk：库存节点l与库存节点k的欧式距离；库存节点成本：nCbT，库存物流费用：不同加工单元综合评价值与物流距离之间的反比函数来反映网络差异化水平：库存节点选址模型中的n,nj均不确定，并且dij,Nj,dij,等变量应随n,nj的变化而变化，为了得到库存节点数目、位置和库存范围等指标，将通过聚类分析的方式通过迭代搜索寻优，找到最优解；步骤三、算法设计；1：将n个库存节点坐标进行标准化处理，使其坐标落在[0～1，0～1]坐标系中，它的位置矩阵为Xn×2；建立综合评价矩阵Y1×n，设定w1,w2,w3初始值；(0)2：将库存节点按照距离进行聚类，建立距离矩阵D＝[dij]m×n(i,j＝1,2,...,n)；如果最小元素为和两类之间的距离，则将和合并成由此构建新类别：进而，转化得到D(b+1)；有4种方法可计算得到D(b+1)，通过比较聚类效果，采用簇平均距离算法进行聚类；簇平均距离算法：式中，da,b表示A类中的样本xa和B类中的样本xb之间的距离；na,nb分别表示A类中的样本xa和B类中的样本xb的个数；DA,B表示A类样本与B类样本之间不同聚类算法的距离；3：令b＝b+1，跳转至2，重复计算及合并，可一直将全部区域运输量最终聚为一类；4：令初始簇数为Ncluster，截取簇，得到Ncluster个库存节点的位置矩阵库存节点和Ncluster个簇；5：当各簇内区域运输量数量lt不大于2或者时，令Ncluster＝Ncluster+1，返回4，否则，计算idu＝1到lt的库存节点idu到簇idu各个区6CN110097322A权　利　要　求　书6/6页域运输量的距离dtij，当dtij＜0.01时客户j的综合评价值qj与dtij的商pij变成pij＝qj×M，通过反复分析计算dtij的取值范围，M取100；最终得到簇idu内pij的平均值Cd(i)＝sum(pij)/It和Ncluster个簇的平均值Cd1(Ncluster)＝sum(C')/Ncluster；6：求Ncluster个库存节点的最短长度Cd2(Ncluster)；构建矩阵任意向量g3×1＝[i,j,dij]，其中i,j为第i,j个库存节点，dij为库存节点i到库存节点j距离；将dij从大到小排列，排序后得到G'；7：从G'的第3行中找到距离最大，且删除之后不影响路径连通的边dij，否则选取另外一边；8：重复步骤7，直到G'只有n-1条边；求得最小路径长度Cd2(Ncluster)，令Cd3(Ncluster)＝Ncluster；9：令搜索向量O最小值，如果最小值为0，需判断是向量O中的空值还是计算值，空值舍去，并在向量O的正元素中寻找最小值；10：输出库存节点的位置坐标。7CN110097322A说　明　书1/12页一种虚拟场景中零件库存节点选址方法技术领域[0001]本发明属于库存节点选址技术领域，具体的说是一种虚拟场景中零件库存节点选址方法。背景技术[0002]设施规划与设计是工业工程学科的重要研究领域之一，它是对系统的各类设施、人员、系统等进行规划和设计，用于优化人流、物流和信息流，从而有效、经济、安全地实现系统的预期目标。[0003]常规的设施选址可以分为单一设施选址、多个设施选址等。通过考虑最低成本、专业化分工协作、分散与集中相结合、重点开发等原则对设置进行选址。常规的选址方法包括优缺点比较法、重心法、线性规划、德尔菲法、费用效果分析法、关联矩阵法、层次分析法、基于遗传算法的选址方法等。[0004]生产库存选址在工厂内更多的是采用SLP等传统方法进行选择，大部分在厂内仅存在1-3个大型零件库存供应区，其它的零件会根据需求，通过物流小车运向指定加工单元。发明内容[0005]本发明提供一种将工厂内的加工单元零件需求量和未来加工单元零件需求量作为评价指标，构建模糊评价矩阵后，利用TOPSIS方法建立出GE矩阵计算出工厂内零件物流的差异化指标，再加上库存节点的距离、库存节点成本，通过基于聚类的重心算法求解出虚拟场景中零件库存节点位置。[0006]本发明技术方案结合附图说明如下：[0007]一种虚拟场景中零件库存节点选址方法，该方法包括以下步骤：[0008]步骤一、加工单元运输量评价与选择；[0009]为了对库存节点进行选址，将工厂内的加工单元零件需求量和未来加工单元零件需求量作为评价指标；加工单元零件需求量影响因素有加工单元与库存区的总物流量A1、加工零件外协比例A2、加工单元物流水平A3、零件运输速率对加工的影响A4、加工单元的设备利用率A5、加工单元的日生产单元体积A6、加工单元零件装卸率A7、加工单元排程情况A8；未来加工单元零件需求量影响因素有加工单元当前设备利用率B1、当前加工单元零件当前销售情况B2、加工单元与库存区的最短距离B3、加工单元当前所占面积B4、工厂的仓库成本B5、加工单元的运输反应时间B6、加工单元的加工停滞时间B7、加工单元的运输率B8；[0010]需要对加工单元零件需求量和未来加工单元零件需求量的影响因素进行Min-max标准化处理，得到指标矩阵A和指标矩阵B；8CN110097322A说　明　书2/12页[0011][0012]a11代表第一行第一列的数据，是加工单元零件需求量影响因素有加工单元与库存区的总物流量的第一个Min-max标准化数据。[0013][0014]b11代表第一行第一列的数据，是未来加工单元零件需求量影响因素有加工单元当前设备利用率的第一个Min-max标准化数据。[0015]首先建立加工单元零件需求量和未来加工单元零件影响关系的网络分析法即ANP结构；建立主要因子的有向图；再经过定位分析得到加工单元零件需求量和未来加工单元零件的主因子P(P1,P2,…,Pm)，次因子P(P11,…,P1n,P21,P22,…,P2l,Pm1,…,Pmk)(Pmk是第m个主因子的第k个次因子)，其中Pm为建工单元需求量和未来加工单元零件元素组；初步构建ANP结构和准则之间的相互影响有向图；构建关系矩阵后，再建立P1,P2y…,Pm之间的直接影响矩阵Z，并将矩阵Z变换成标准矩阵Y；[0016]Z＝[zij]m×mi,j＝1,2,…,n  (1.3)[0017]其中，zij表示在主因子Pi对主因子Pj的影响程度[0018]Y＝k·Z,k＞0[0019]其中，[0020]计算最终关系矩阵，构建ANP主要因子结构图；建立总影响矩阵T，经过公式1.4转化得到Ti,i＝0,1,2,...,m；根据影响因素实际关系，选取临界值p，经过计算得到最后的关系矩阵G，接着再构建ANP准则间结构；[0021]T＝Y·(I-Y)-1  (1.4)[0022]T＝[tij]m×m  (1.5)[0023]G＝[gij]m×m  (1.6)[0024]其中(Y为标准矩阵，tij由式1.4换算所得)[0025]根据ANP主因子间新结构，构建Pk主因子下Pk1,Pk2…,Pkm，k∈{1,2,...,m}之间的直接关系矩阵，得到完整ANP结构；9CN110097322A说　明　书3/12页[0026]建立ANP判断矩阵；构建决策t(t＝1,2,3,...,n)的判断矩阵；主因子P1,P2…,Pm对t(jl)主因子Pj判断矩阵为A和主因子P1的元素Pi1,Pi2…,Pm对主因子Pj元素Pjl影响的判断矩阵为[0027]计算加权矩阵，对判断矩阵Ai(jl)求归化特征向量进而得到矩阵[0028]同理求得At(jl)的权重矩阵At；再求超矩阵Wt和加权矩阵[0029][0030][0031]t[0032]截取Pm与P11,P12…,P1n,P21,P22…,P2l…,P(m-1)1,P(m-1)2…,P(m-1)k之间的加权矩阵E，规范化后得到Rt；[0033][0034][0035]其中为加权矩阵Et中的第i行第j列数值；[0036]再构建模糊评价矩阵；将全部评分人的规范化评价矩阵进行模糊处理，从而获得模糊评价矩阵R；[0037]10CN110097322A说　明　书4/12页[0038]其中，[0039]确定理想解和负理想解；根据TOPSIS方法得到理想解向量和负理想解向量；[0040]理想解：[0041][0042]负理想解：[0043][0044]其中J1为效益型指标集，J2为成本型指标集；[0045]计算贴近度C；[0046]C＝[c1]1×m；[0047][0048]其中[0049]进而，根据公式1.11-1.14，求解第i部分需求量的贴近度和得到加工单元零件需求量和未来加工单元零件需求量的评价值，再建立“零件需求量和未来零件需求量”的GE矩阵，通过加和得到加工单元的评价值；[0050]步骤二、建立库存节点选址方法；[0051]将从厂内物流成本和各单元的实际生产情况出发，选择库存节点的距离，库存节点成本和加工单元的评价值作为工厂库存节点的评价标准；[0052]参阅图2、图3，库存节点之间的连接关系分为单连通结构和直通式结构；因此构建两种模型如下：[0053]单连通模型：[0054][0055][0056]s.t.[0057]2≤n≤nc[0058]2≤nj≤nc[0059]11CN110097322A说　明　书5/12页[0060][0061]式中[0062]直通式模型：[0063][0064][0065]s.t.[0066]2≤n≤nc[0067]2≤nj≤nc[0068][0069][0070]式中，12[0071]f(x)为效益型函数，f(x)为成本型函数，将w1、w2、w3设置成权重，库存节点选址模型的函数就转化为：[0072]单连通模型：[0073]直通式模型：[0074]其中：n：库存建设节点数目；Cbc：库存物流单位长度费用；CbT：库存建设数额；为第k次算得的库存节点l与库存节点m间的长度，Nj：加工单元运输量集，Nj＝1,2,...,nc且i,j∈Nj；dlj：库存节点l与加工单元j的距离，dlj＝|xl-xj|+|yl-yj|；θj：加工单元j的综合评价值；Il：库存节点l的加工单元运输量集合；L：库存节点集，L＝1,2,...,n且l,k∈L；nc：加工单元数目；dmax：库存节点最大服务范围；dmin：加工单元库存最小距离；Mmin：极小值，可取值为标准化的dmin；库存建设投资额；nl：库存节点l输送的加工单元运送12CN110097322A说　明　书6/12页量；xl：库存节点l的横坐标；yl：库存节点l的横坐标。dlk：库存节点l与库存节点k的欧式距离；[0075]库存节点成本：nCbT，库存物流费用：不同加工单元综合评价值与物流距离之间的反比函数来反映网络差异化水平：[0076]库存节点选址模型中的n,nj均不确定，并且等变量应随n,nj的变化而变化，为了得到库存节点数目、位置和库存范围等指标，将通过聚类分析的方式通过迭代搜索寻优，找到最优解；[0077]步骤三、算法设计；[0078]1：将n个库存节点坐标进行标准化处理，使其坐标落在[0～1，0～1]坐标系中，它的位置矩阵为Xn×2；建立综合评价矩阵Y1×n，设定w1,w2,w3初始值。(0)[0079]2：将库存节点按照距离进行聚类，建立距离矩阵D＝[dij]m×n(i,j＝1,2,...,n)；如果最小元素为和两类之间的距离，则将和合并成由此构建新类别：进而，转化得到D(b+1)；有4种方法可计算得到D(b+1)，通过比较聚类效果，采用簇平均距离算法进行聚类；簇平均距离算法：[0080][0081]式中，da,b表示A类中的样本xa和B类中的样本xb之间的距离；na,nb分别表示A类中的样本xa和B类中的样本xb的个数；DA,B表示A类样本与B类样本之间不同聚类算法的距离；[0082]3：令b＝b+1，跳转至2，重复计算及合并，可一直将全部区域运输量最终聚为一类；[0083]4：令初始簇数为Ncluster，截取簇，得到Ncluster个库存节点的位置矩阵库存节点和Ncluster个簇；[0084]5：当各簇内区域运输量数量lt不大于2或者时，令Ncluster＝Ncluster+1，返回4，否则，计算idu＝1到lt的库存节点idu到簇idu各个区域运输量的距离dtij，当dtij＜0.01时客户j的综合评价值qj与dtij的商pij变成pij＝qj×M，通过反复分析计算dtij的取值范围，M取100；最终得到簇idu内pij的平均值Cd(i)＝sum(pij)/It和Ncluster个簇的平均值Cd1(Ncluster)＝sum(C')/Ncluster；[0085]6：求Ncluster个库存节点的最短长度Cd2(Ncluster)；构建矩阵任意向量g3×1＝[i,j,dij]，其中i,j为第i,j个库存节点，dij为库存节点i到库存节点j距离；将dij从大到小排列，排序后得到G'；[0086]7：从G'的第3行中找到距离最大，且删除之后不影响路径连通的边dij，否则选取另外一边；[0087]8：重复步骤7，直到G'只有n-1条边；求得最小路径长度Cd2(Ncluster)，令Cd313CN110097322A说　明　书7/12页(Ncluster)＝Ncluster；[0088]9：令搜索向量O最小值，如果最小值为0，需判断是向量O中的空值还是计算值，空值舍去，并在向量O的正元素中寻找最小值；[0089]10：输出库存节点的位置坐标。[0090]本发明的有益效果为：[0091]本发明能对于工厂内库存多级节点选择提供数学求解方法，能通过减少工厂内物流量，有效降低厂内的物流部门的工作人员和叉车数量，从而降低物流成本。工厂库存选址领域暂未有相关技术研究，本方法提供一种技术方案填补当前空白。附图说明[0092]图1为本发明的布局结构示意图；[0093]图2为本发明的单连通库存之间转换结构图；[0094]图3为本发明的库存之间物流网络结构图。具体实施方式[0095]1、一种虚拟场景中零件库存节点选址方法，其特征在于，该方法包括以下步骤：[0096]步骤一、加工单元运输量评价与选择；[0097]为了对库存节点进行选址，将工厂内的加工单元零件需求量和未来加工单元零件需求量作为评价指标；加工单元零件需求量影响因素有加工单元与库存区的总物流量A1、加工零件外协比例A2、加工单元物流水平A3、零件运输速率对加工的影响A4、加工单元的设备利用率A5、加工单元的日生产单元体积A6、加工单元零件装卸率A7、加工单元排程情况A8；未来加工单元零件需求量影响因素有加工单元当前设备利用率B1、当前加工单元零件当前销售情况B2、加工单元与库存区的最短距离B3、加工单元当前所占面积B4、工厂的仓库成本B5、加工单元的运输反应时间B6、加工单元的加工停滞时间B7、加工单元的运输率B8；[0098]参阅图1，需要对加工单元零件需求量和未来加工单元零件需求量的影响因素进行处理，得到指标矩阵A和指标矩阵B；[0099][0100][0101]其中，a11代表第一行第一列的数据，是加工单元零件需求量影响因素有加工单元14CN110097322A说　明　书8/12页与库存区的总物流量的第一个Min-max标准化数据；[0102]b11代表第一行第一列的数据，是未来加工单元零件需求量影响因素有加工单元当前设备利用率的第一个Min-max标准化数据；[0103]首先建立加工单元零件需求量和未来加工单元零件影响关系的网络分析法即ANP结构；建立主要因子的有向图；再经过定位分析得到加工单元零件需求量和未来加工单元零件的主因子P(P1,P2,…,Pm)，次因子P(P11,…,P1n,P21,P22,…,P2l,Pm1,…,Pmk)，其中Pm为建工单元需求量和未来加工单元零件元素组；初步构建ANP结构和准则之间的相互影响有向图；构建关系矩阵后，再建立P1,P2y…,Pm之间的直接影响矩阵Z，并将矩阵Z变换成标准矩阵Y；[0104]Z＝[zij]m×mi,j＝1,2,…,n  (1.3)[0105]其中，zij表示在主因子Pi对主因子Pj的影响程度[0106]Y＝k·Z,k＞0[0107]其中，[0108]计算最终关系矩阵，构建ANP主要因子结构图；建立总影响矩阵T，经过公式1.4转化得到Ti,i＝0,1,2,...,m；根据影响因素实际关系，选取临界值p，经过计算得到最后的关系矩阵G，接着再构建ANP准则间结构；[0109]T＝Y·(I-Y)-1  (1.4)[0110]T＝[tij]m×m  (1.5)[0111]G＝[gij]m×m  (1.6)[0112]其中[0113]根据ANP主因子间新结构，构建Pk主因子下Pk1,Pk2…,Pkm，k∈{1,2,...,m}之间的直接关系矩阵，得到完整ANP结构；[0114]建立ANP判断矩阵；构建决策t(t＝1,2,3,...,n)的判断矩阵；主因子P1,P2…,Pm对t(jl)主因子Pj判断矩阵为A和主因子P1的元素Pi1,Pi2…,Pm对主因子Pj元素Pjl影响的判断矩阵为[0115]计算加权矩阵，对判断矩阵Ai(jl)求归化特征向量进而得到矩阵[0116]同理求得At(jl)的权重矩阵At；再求超矩阵Wt和加权矩阵15CN110097322A说　明　书9/12页[0117][0118][0119]t[0120]截取Pm与P11,P12…,P1n,P21,P22…,P2l…,P(m-1)1,P(m-1)2…,P(m-1)k之间的加权矩阵E，规范化后得到Rt；[0121][0122]t[0123]其中eij为加权矩阵E中的第i行第j列数值；[0124]再构建模糊评价矩阵；将全部评分人的规范化评价矩阵进行模糊处理，从而获得模糊评价矩阵R；[0125][0126]其中，[0127]确定理想解和负理想解；根据TOPSIS方法得到理想解向量和负理想解向量；TOPSIS(Technique for Order Preference by Similarity to an Ideal Solution)法是C.L.Hwang和K.Yoon于1981年首次提出，TOPSIS法根据有限个评价对象与理想化目标的接近程度进行排序的方法，是在现有的对象中进行相对优劣的评价。[0128]理想解：16CN110097322A说　明　书10/12页[0129][0130]负理想解：[0131][0132]其中J1为效益型指标集，J2为成本型指标集；[0133]计算贴近度C；[0134]C＝[c1]1×m；[0135][0136]其中[0137]进而，根据公式1.11-1.14，求解第i部分需求量的贴近度和得到加工单元零件需求量和未来加工单元零件需求量的评价值，再建立“零件需求量和未来零件需求量”的GE矩阵，通过加和得到加工单元的评价值；[0138]步骤二、建立库存节点选址方法；[0139]将从厂内物流成本和各单元的实际生产情况出发，选择库存节点的距离，库存节点成本和加工单元的评价值作为工厂库存节点的评价标准；[0140]库存节点之间的连接关系分为单连通结构和直通式结构；因此构建两种模型如下：[0141]单连通模型：[0142][0143][0144]s.t.[0145]2≤n≤nc[0146]2≤nj≤nc[0147][0148][0149]式中17CN110097322A说　明　书11/12页[0150]直通式模型：[0151][0152][0153]s.t.[0154]2≤n≤nc[0155]2≤nj≤nc[0156][0157][0158]式中，12[0159]f(x)为效益型函数，f(x)为成本型函数，将w1、w2、w3设置成权重，库存节点选址模型的函数就转化为：[0160]单连通模型：[0161]直通式模型：[0162]其中：n：库存建设节点数目；Cbc：库存物流单位长度费用；CbT：库存建设数额；为第k次算得的库存节点l与库存节点m间的长度，Nj：加工单元运输量集，Nj＝1,2,...,nc且i,j∈Nj；dlj：库存节点l与加工单元j的距离，dlj＝|xl-xj|+|yl-yj|；θj：加工单元j的综合评价值；Il：库存节点l的加工单元运输量集合；L：库存节点集，L＝1,2,...,n且l,k∈L；nc：加工单元数目；dmax：库存节点最大服务范围；dmin：加工单元库存最小距离；Mmin：极小值，可取值为标准化的dmin；库存建设投资额；nl：库存节点l输送的加工单元运送量；xl：库存节点l的横坐标；yl：库存节点l的横坐标。dlk：库存节点l与库存节点k的欧式距离；[0163]库存节点成本：nCbT，库存物流费用：不同加工单元综合评价值与物流距离之间的反比函数来反映网络差异化水平：18CN110097322A说　明　书12/12页[0164]库存节点选址模型中的n,nj均不确定，并且等变量应随n,nj的变化而变化，为了得到库存节点数目、位置和库存范围等指标，将通过聚类分析的方式通过迭代搜索寻优，找到最优解；[0165]步骤三、算法设计；[0166]1：将n个库存节点坐标进行标准化处理，使其坐标落在[0～1，0～1]坐标系中，它的位置矩阵为Xn×2；建立综合评价矩阵Y1×n，设定w1,w2,w3初始值。(0)[0167]2：将库存节点按照距离进行聚类，建立距离矩阵D＝[dij]m×n(i,j＝1,2,...,n)；如果最小元素为和两类之间的距离，则将和合并成由此构建新类别：进而，转化得到D(b+1)；有4种方法可计算得到D(b+1)，通过比较聚类效果，采用簇平均距离算法进行聚类；簇平均距离算法：[0168][0169]式中，da,b表示A类中的样本xa和B类中的样本xb之间的距离；na,nb分别表示A类中的样本xa和B类中的样本xb的个数；DA,B表示A类样本与B类样本之间不同聚类算法的距离；[0170]3：令b＝b+1，跳转至2，重复计算及合并，可一直将全部区域运输量最终聚为一类；[0171]4：令初始簇数为Ncluster，截取簇，得到Ncluster个库存节点的位置矩阵库存节点和Ncluster个簇；[0172]5：当各簇内区域运输量数量lt不大于2或者时，令Ncluster＝Ncluster+1，返回4，否则，计算idu＝1到lt的库存节点idu到簇idu各个区域运输量的距离dtij，当dtij＜0.01时客户j的综合评价值qj与dtij的商pij变成pij＝qj×M，通过反复分析计算dtij的取值范围，M取100；最终得到簇idu内pij的平均值Cd(i)＝sum(pij)/It和Ncluster个簇的平均值Cd1(Ncluster)＝sum(C')/Ncluster；[0173]6：求Ncluster个库存节点的最短长度Cd2(Ncluster)；构建矩阵任意向量g3×1＝[i,j,dij]，其中i,j为第i,j个库存节点，dij为库存节点i到库存节点j距离；将dij从大到小排列，排序后得到G'；[0174]7：从G'的第3行中找到距离最大，且删除之后不影响路径连通的边dij，否则选取另外一边；[0175]8：重复步骤7，直到G'只有n-1条边；求得最小路径长度Cd2(Ncluster)，令Cd3(Ncluster)＝Ncluster；[0176]9：令搜索向量O最小值，如果最小值为0，需判断是向量O中的空值还是计算值，空值舍去，并在向量O的正元素中寻找最小值；[0177]10：输出库存节点的位置坐标。19CN110097322A说　明　书　附　图1/2页图1图220CN110097322A说　明　书　附　图2/2页图321