发明名称 一种外贸物流路径的优化方法
摘要 本发明涉及一种外贸物流路径的优化方法,其解决了目前计算方法事实操作可行性不强,面对大规模的物流网络无法精确且成本高,无自适用学能力。其通过建立无向图模型推到出期望值模型进行计算得出最优路径。其可广泛应用物流运输领域。
申请公布号 CN103399868B 申请公布日期 2016.08.10
申请号 CN201310284966.8 申请日期 2013.07.08
申请人 哈尔滨工业大学(威海) 发明人 初佃辉;叶允明;李春山;周学权;王德泉
分类号 G06F17/30(2006.01)I;G06Q10/08(2012.01)I 主分类号 G06F17/30(2006.01)I
代理机构 北京怡丰知识产权代理有限公司 11293 代理人 于振强
主权项 一种外贸物流路径的优化方法,其特征在于包括以下步骤:(1)识别外贸物流网络中的实体,所述实体包括:起点,代理,第一运输公司,仓储站,第二运输公司和目的地;(2)构建外贸物流网络;(3)将所述外贸物流网络抽象为层次无向图;(4)针对所述层次无向图,面对不同场景,计算最优路径;所述层次无向图中,不同层次的节点抽象为无向图的节点集合V,不同层次节点之间的连线抽象为无向图边的集合E,形成G=&lt;V,E&gt;;节点集合V可以分为k个不相交的子集:V=S<sub>1</sub>∪S<sub>2</sub>∪…∪S<sub>k</sub>,用V<sub>i</sub>来表示一个节点,则有V<sub>i</sub>∈S<sub>i</sub>,边集E={<V<sub>i</sub>,V<sub>i+1</sub>|V<sub>i</sub>∈S<sub>i</sub>,V<sub>i+1</sub>∈S<sub>i+1</sub>>},i∈{1,...,k‑1},每条边e=&lt;v<sub>i</sub>,v<sub>j</sub>&gt;的权为ω<sub>ij</sub>表示该边在物流路径上的概率;所述步骤(4)中的计算方法如下:所述层次无向图中,每对&lt;S<sub>i</sub>,S<sub>i+1</sub>&gt;生成一个转移概率矩阵M;所述外贸物流网络存在6个实体,每一对实体之间有一个转移概率矩阵,因此,存在5个转移概率矩阵:M<sub>1</sub>,M<sub>2</sub>,M<sub>3</sub>,M<sub>4</sub>和M<sub>5</sub>,矩阵M<sub>i</sub>为&lt;S<sub>i</sub>,S<sub>i+1</sub>&gt;上的转移概率矩阵;定义如下的优化路径迭代方程:<maths num="0001"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msup><mi>v</mi><mn>1</mn></msup><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>c</mi><mo>)</mo></mrow><msub><mi>M</mi><mn>1</mn></msub><msup><mi>v</mi><mn>2</mn></msup><mo>+</mo><mi>c</mi><mi>p</mi></mrow></mtd></mtr><mtr><mtd><mrow><msup><mi>v</mi><mn>2</mn></msup><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><msub><mi>M</mi><mn>1</mn></msub><mi>T</mi></msup><msup><mi>v</mi><mn>1</mn></msup><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>M</mi><mn>2</mn></msub><msup><mi>v</mi><mn>3</mn></msup></mrow></mtd></mtr><mtr><mtd><mrow><msup><mi>v</mi><mn>3</mn></msup><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><msub><mi>M</mi><mn>2</mn></msub><mi>T</mi></msup><msup><mi>v</mi><mn>2</mn></msup><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>M</mi><mn>3</mn></msub><msup><mi>v</mi><mn>4</mn></msup></mrow></mtd></mtr><mtr><mtd><mrow><msup><mi>v</mi><mn>4</mn></msup><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><msub><mi>M</mi><mn>3</mn></msub><mi>T</mi></msup><msup><mi>v</mi><mn>3</mn></msup><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>M</mi><mn>4</mn></msub><msup><mi>v</mi><mn>5</mn></msup></mrow></mtd></mtr><mtr><mtd><mrow><msup><mi>v</mi><mn>5</mn></msup><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><msub><mi>M</mi><mn>5</mn></msub><mi>T</mi></msup><msup><mi>v</mi><mn>3</mn></msup><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>M</mi><mn>5</mn></msub><msup><mi>v</mi><mn>5</mn></msup></mrow></mtd></mtr><mtr><mtd><mrow><msup><mi>v</mi><mn>6</mn></msup><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>c</mi><mo>)</mo></mrow><msub><mi>M</mi><mn>5</mn></msub><msup><mi>v</mi><mn>2</mn></msup><mo>+</mo><mi>c</mi><mi>q</mi></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000985660050000011.GIF" wi="1541" he="437" /></maths>向量p和q具有初始值,c是一个常数,根据经验设定为0.5,v<sup>i</sup>是一个向量,v<sup>i</sup>中的每一个维度分量v<sub>h</sub><sup>i</sup>表示S<sub>i</sub>中第h个节点被选中成为外贸物流路径一个结点的可能性;转移概率定义为:<maths num="0002"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>f</mi><mi>&psi;</mi></msub><mrow><mo>(</mo><msub><mi>v</mi><mi>i</mi></msub><mo>,</mo><msub><mi>v</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>F</mi><mo>(</mo><mrow><msub><mi>v</mi><mi>i</mi></msub><mo>,</mo><msub><mi>v</mi><mi>j</mi></msub><mo>,</mo><mi>&psi;</mi></mrow><mo>)</mo><mo>)</mo></mrow></mrow></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>&psi;</mi><mi>k</mi></msub><msup><mrow><mo>(</mo><mrow><msubsup><mi>v</mi><mi>i</mi><mi>k</mi></msubsup><mo>-</mo><msubsup><mi>v</mi><mi>j</mi><mi>k</mi></msubsup></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow></mrow></mfrac><mo>.</mo></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000985660050000021.GIF" wi="1623" he="300" /></maths>参数ψ可以采用最大似然估计,对数似然函数是:<img file="FDA0000985660050000022.GIF" wi="751" he="123" />其中m是边的数量;最大化步骤可表示为:<maths num="0003"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mfrac><mrow><mo>&part;</mo><mi>F</mi><mrow><mo>(</mo><msub><mi>v</mi><mi>i</mi></msub><mo>,</mo><msub><mi>v</mi><mi>j</mi></msub><mo>,</mo><mi>&psi;</mi><mo>)</mo></mrow></mrow><mrow><mo>&part;</mo><msub><mi>&psi;</mi><mi>k</mi></msub></mrow></mfrac><mo>=</mo><mi>s</mi><mi>i</mi><mi>m</mi><mrow><mo>(</mo><msubsup><mi>v</mi><mi>i</mi><mi>k</mi></msubsup><mo>,</mo><msubsup><mi>v</mi><mi>j</mi><mi>k</mi></msubsup><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><mrow><mo>&part;</mo><mi>L</mi></mrow><mrow><mo>&part;</mo><msub><mi>&psi;</mi><mi>k</mi></msub></mrow></mfrac><mo>=</mo><mrow><mo>(</mo><mo>-</mo><mn>1</mn><mo>)</mo></mrow><munder><mo>&Sigma;</mo><mi>m</mi></munder><mfrac><mrow><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>F</mi><mo>(</mo><mrow><msub><mi>v</mi><mi>i</mi></msub><mo>,</mo><msub><mi>v</mi><mi>j</mi></msub><mo>,</mo><mi>&psi;</mi></mrow><mo>)</mo><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>+</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>F</mi><mo>(</mo><mrow><msub><mi>v</mi><mi>i</mi></msub><mo>,</mo><msub><mi>v</mi><mi>j</mi></msub><mo>,</mo><mi>&psi;</mi></mrow><mo>)</mo><mo>)</mo></mrow></mrow></mfrac><mfrac><mrow><mo>&part;</mo><mi>F</mi><mrow><mo>(</mo><msub><mi>v</mi><mi>i</mi></msub><mo>,</mo><msub><mi>v</mi><mi>j</mi></msub><mo>,</mo><mi>&psi;</mi><mo>)</mo></mrow></mrow><mrow><mo>&part;</mo><msub><mi>&psi;</mi><mi>k</mi></msub></mrow></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mrow><mo>(</mo><mo>-</mo><mn>1</mn><mo>)</mo></mrow><munder><mo>&Sigma;</mo><mi>m</mi></munder><mfrac><mrow><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>F</mi><mo>(</mo><mrow><msub><mi>v</mi><mi>i</mi></msub><mo>,</mo><msub><mi>v</mi><mi>j</mi></msub><mo>,</mo><mi>&psi;</mi></mrow><mo>)</mo><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>+</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>F</mi><mo>(</mo><mrow><msub><mi>v</mi><mi>i</mi></msub><mo>,</mo><msub><mi>v</mi><mi>j</mi></msub><mo>,</mo><mi>&psi;</mi></mrow><mo>)</mo><mo>)</mo></mrow></mrow></mfrac><mi>s</mi><mi>i</mi><mi>m</mi><mrow><mo>(</mo><msubsup><mi>v</mi><mi>i</mi><mi>k</mi></msubsup><mo>,</mo><msubsup><mi>v</mi><mi>j</mi><mi>k</mi></msubsup><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&psi;</mi><mi>t</mi></msub><mo>=</mo><msub><mi>&psi;</mi><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><mi>&eta;</mi><mfrac><mrow><mo>&part;</mo><mi>L</mi></mrow><mrow><mo>&part;</mo><mi>&psi;</mi></mrow></mfrac><mo>,</mo></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000985660050000023.GIF" wi="1722" he="475" /></maths>其中η是迭代参数,当|ψ<sub>t</sub>‑ψ<sub>t‑1</sub>|小于预定义的ε值时,迭代结束;A:基于随机游走模型的运输路线算法当没有特殊要求或约束时,可以通过RWTR算法获得最佳路径;该算法直接采用随机游走模型优化路径,算法如下所示:算法的输入项为:出口商品属性P,业务数据集D;算法的输出项为:图G中的路径和转移概率矩阵;算法的运行步骤为:1)生成转移概率矩阵M<sub>1</sub>,M<sub>2</sub>,M<sub>3</sub>,M<sub>4</sub>,M<sub>5</sub>;2)根据出口商品属性P生成向量p和q;3)初始化向量v<sup>1</sup>,v<sup>2</sup>,v<sup>3</sup>,v<sup>4</sup>,v<sup>5</sup>,v<sup>6</sup>;4)判断v<sup>1</sup>,v<sup>2</sup>,v<sup>3</sup>,v<sup>4</sup>,v<sup>5</sup>,v<sup>6</sup>的值是否收敛,如果收敛则跳转到步骤6);5)根据所述优化路径迭代方程来更新v<sup>1</sup>,v<sup>2</sup>,v<sup>3</sup>,v<sup>4</sup>,v<sup>5</sup>,v<sup>6</sup>的值;6)在向量v<sup>i</sup>中,找到值最大的维度f,那么节点<img file="FDA0000985660050000024.GIF" wi="59" he="75" />被选中作为优化路径中的一个节点;7)根据步骤6),在v<sup>2</sup>,v<sup>3</sup>,v<sup>4</sup>,v<sup>5</sup>,v<sup>6</sup>上获取整个优化路径的节点;或者B:面向约束的运输路径算法,如下所示的CTR算法算法的输入项为:出口商品属性P,概率矩阵:M<sub>1</sub>,M<sub>2</sub>,M<sub>3</sub>,M<sub>4</sub>,M<sub>5</sub>,约束L={l<sub>i</sub>(s<sub>j</sub>),s<sub>j</sub>∈V‑S<sup>1</sup>},其中S<sup>1</sup>表示外贸物流中的所有出发点组成的集合;算法的输出项为:图G中的路径;算法的运行步骤为:Ⅰ、根据出口商品属性P生成向量p和q;Ⅱ、初始化向量v<sup>1</sup>,v<sup>2</sup>,v<sup>3</sup>,v<sup>4</sup>,v<sup>5</sup>,v<sup>6</sup>;Ⅲ、判断v<sup>1</sup>,v<sup>2</sup>,v<sup>3</sup>,v<sup>4</sup>,v<sup>5</sup>,v<sup>6</sup>的值是否收敛,如果收敛则跳转到步骤Ⅴ;Ⅳ、根据所述优化路径迭代方程来更新v<sup>1</sup>,v<sup>2</sup>,v<sup>3</sup>,v<sup>4</sup>,v<sup>5</sup>,v<sup>6</sup>的值;Ⅴ、将向量v<sup>i</sup>中的每一个维度按照值的大小进行降序排列,顺序得找到值最大且满足约束l<sub>i</sub>的维度h,那么节点<img file="FDA0000985660050000032.GIF" wi="57" he="73" />被选中作为优化路径中的一个节点;Ⅵ、根据步骤Ⅴ,在v<sup>2</sup>,v<sup>3</sup>,v<sup>4</sup>,v<sup>5</sup>,v<sup>6</sup>上获取整个优化路径的节点;或者C:增量式运输路径算法一旦接收到增量数据集D<sub>t</sub>,转换概率矩阵M1,M2,M3,M4,M5被公式(4)更新;<maths num="0004"><math><![CDATA[<mrow><msup><msub><mi>M</mi><mi>i</mi></msub><mo>&prime;</mo></msup><mo>=</mo><mfrac><mrow><mo>|</mo><mi>D</mi><mo>|</mo></mrow><mrow><mo>|</mo><mi>D</mi><mo>|</mo><mo>+</mo><mo>|</mo><msub><mi>D</mi><mi>t</mi></msub><mo>|</mo></mrow></mfrac><msub><mi>M</mi><mi>i</mi></msub><mo>+</mo><mfrac><mrow><mo>|</mo><msub><mi>D</mi><mi>t</mi></msub><mo>|</mo></mrow><mrow><mo>|</mo><mi>D</mi><mo>|</mo><mo>+</mo><mo>|</mo><msub><mi>D</mi><mi>t</mi></msub><mo>|</mo></mrow></mfrac><msub><mi>M</mi><mrow><mi>i</mi><mi>t</mi></mrow></msub><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000985660050000031.GIF" wi="1518" he="143" /></maths>其中,M<sub>it</sub>是由D<sub>t</sub>计算出来的转换概率矩阵;如下所示的ICTR算法:算法的输入项为:出口商品属性P,业务数据集D,新业务数据D<sub>t</sub>,约束L={l<sub>i</sub>(s<sub>j</sub>),s<sub>j</sub>∈V‑S<sup>1</sup>},其中S<sup>1</sup>表示外贸物流中的所有出发点组成的集合;算法的输出项为:图G中的路径;算法的运行步骤为:步骤1,根据新业务数据生成转移概率矩阵M<sub>1</sub>’,M<sub>2</sub>’,M<sub>3</sub>’,M<sub>4</sub>’,M<sub>5</sub>’;步骤2,根据出口商品属性P生成向量p和q;步骤3,初始化向量v<sup>1</sup>,v<sup>2</sup>,v<sup>3</sup>,v<sup>4</sup>,v<sup>5</sup>,v<sup>6</sup>;步骤4,判断v<sup>1</sup>,v<sup>2</sup>,v<sup>3</sup>,v<sup>4</sup>,v<sup>5</sup>,v<sup>6</sup>的值是否收敛,如果收敛则跳转到步骤6;步骤5,根据所述优化路径迭代方程来更新v<sup>1</sup>,v<sup>2</sup>,v<sup>3</sup>,v<sup>4</sup>,v<sup>5</sup>,v<sup>6</sup>的值;步骤6,将向量v<sup>i</sup>中的每一个维度按照值的大小进行降序排列,顺序得找到值最大且满足约束l<sub>i</sub>的维度z,那么节点<img file="FDA0000985660050000041.GIF" wi="60" he="82" />被选中作为优化路径中的一个节点;步骤7,根据步骤6,在v<sup>2</sup>,v<sup>3</sup>,v<sup>4</sup>,v<sup>5</sup>,v<sup>6</sup>上获取整个优化路径的节点。
地址 264209 山东省威海市文化西路2号
您可能感兴趣的专利