基于Brusselator的物流配送中心系统稳定性分析.docxVIP

基于Brusselator的物流配送中心系统稳定性分析布鲁塞尔器机理简介在耗散结构理论中，普利高津等人曾从实际的化学反应的动力学过程中总结和概括出一个动力学模型，即著名的布鲁塞尔器(Brusselator)模型，其模型如下：式中：A、B——为反应物，在反应过程中不断被消耗并能充分补充；D、E——为生成物，在反应过程中不断被取走，以使反应能正常进行；X、Y——为中间产物，它们可以有随时间变化的浓度。——为反应中的催化剂，它的数量多少影响反映速度。这里最重要的是反应系统中的第三个自催化环节,他导致非线性,而这对系统的特殊行为至关重要。因为在第三环节上形成了3个X分子，所以这个模型也称之为三分子模型。这种反应路径可以用图1的环来表示。由图1可知布鲁塞尔器的反应路径，从A合成X，X进而分解为E，生成物X与一个交叉催化关系相联而生成Y。X在一个三分子步骤中由Y生成，但反过来，X和B的反应又合成Y。于是有如下方程组：式中字母X,Y,A,B均表示相应物质的数量或是浓度，，分别表示X,Y两种物质的扩散系数。为方便计算，可以略去扩散项，相当于均匀搅拌的反应。进行变量转换，取式(1)转换为式(2)：方程定态解为，表示系统存在均匀定态。利用线性稳定性分析方法解多维非线性系统稳定性，令，，为小量。在方程右端附近可展开多元函数的泰勒级数，得式(3):其特征方程为：根据的值确定系统的稳定性，方程的解为。物流配送中心服务能力形成机理分析2.1模型转义进行“布鲁塞尔器”模型转义即将模型中A,B,D,E,X,Y所代表的意义转变为物流配送中心的服务能力和服务能力形成过程中的因素。设A——为物流需求（可用指标“货物重量”来量化，单位：吨）； B——为引进物流技术（可用指标“投入资金”来量化，单位：万元）；D——收益增加（单位：万元）； E——物流服务能力（可用指标“顾客满意度”来量化）；X——车辆调度模式1（可用指标“总配送时间”来量化，单位：小时）；Y——车辆调度模式2（可用指标“总配送时间”来量化，单位：小时）于是配送中心的物流服务能力可以用布鲁塞尔器模型描述，其模型如下：2.2模型解释对于转义后的模型可以这样理解：某物流配送中心在有物流需求A的情况下，采用车辆调度模式1即X来进行配送，这样产生E，即形成一定的物流服务能力；后来由于引进了新的物流技术B（例如一种新的路径优化技术、信息系统管理技术等），在X的基础上产生了新的车辆调度模式2即Y,缩短了总配送时间，带来收益的增加；而配送中心通过加强对司机的管理（如GPRS跟踪定位），以及司机师傅自身对路况的不断摸索，X从而实现一个自催化的过程，使得总配送时间减少，增强物流服务能力。2.3赋值计算及稳定性分析 本文已将各个指标量化，以便进行实证研究，由于时间精力的有限，在此仅通过赋值计算来分析不同情况下系统的稳定性。(1)令，于是皆为正。对定态解的任意扰动，使，远离定态解，为不稳定节点，也就是说该配送中心的车辆调度模式是不稳定的，节点附近随时间变动的具体情况如图1所示。(2)令，于是具有正实部的共轭复根。随着时间的增长，，也变大，此时定态解为不稳定焦点，焦点附近随时间变化情况如图2所示。系统会出现耗散结构分支，进而形成新的有序结构。当系统失稳后，由于系统内部非线性因素的相互作用，当系统的涨落达到临界时，系统会从初始状态演化至新的状态，并不断向高级有序方向发展，如图3所示。说明该配送中心的车辆调度模式会不断的更新。(3)令，于是皆为正，对定态解的任意扰动，，仍会趋近定态解，称此时定态解为稳定的结点。得到与情况(1)一样的结果，只是演化的方向趋近定态解，结点附近随时间变动的具体情况如图4所示。(4)令，于是具有负实部的共轭复根。随着时间的增长，，变小且趋向于零，但具有震荡特点。此时定态解为稳定焦点，如图5所示。说明此时的两种车辆调度模式是稳定的，尽管会有外界小小的扰动，但随着时间的流逝，在发生震荡后最终会恢复当前状态。