CN113657827A - 一种物流系统切换建模与滤波方法 - Google Patents

一种物流系统切换建模与滤波方法 Download PDF

Info

Publication number
CN113657827A
CN113657827A CN202110914386.7A CN202110914386A CN113657827A CN 113657827 A CN113657827 A CN 113657827A CN 202110914386 A CN202110914386 A CN 202110914386A CN 113657827 A CN113657827 A CN 113657827A
Authority
CN
China
Prior art keywords
switching
logistics system
matrix
filter
time
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202110914386.7A
Other languages
English (en)
Inventor
李临敏
张俊锋
李强
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hangzhou Dianzi University
Original Assignee
Hangzhou Dianzi University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hangzhou Dianzi University filed Critical Hangzhou Dianzi University
Priority to CN202110914386.7A priority Critical patent/CN113657827A/zh
Publication of CN113657827A publication Critical patent/CN113657827A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/08Logistics, e.g. warehousing, loading or distribution; Inventory or stock management
    • G06Q10/083Shipping
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Business, Economics & Management (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Economics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Quality & Reliability (AREA)
  • Algebra (AREA)
  • Development Economics (AREA)
  • Computing Systems (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Operations Research (AREA)
  • Software Systems (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明属于物流工程技术领域,公开了一种物流系统切换建模与滤波方法。本发明基于正切换系统建模方法、事件触发策略和异步滤波器技术,针对物流系统中庞大的快递量进行数据采集,提出了一种物流系统中快递高效传输的方法,该方法考虑了网购狂欢节带来的庞大的快递量问题,根据快递量和各个转运中心的压力,利用切换技术为各个转运中心合理的分配任务量,从而提高了处理快递的效率,也考虑了恶劣天气等其它影响物流系统正常运转的因素,从而大大的提高了物流系统的运转性能,及时的派遣快递到对应的收货人手中。

Description

一种物流系统切换建模与滤波方法
技术领域
本发明属于物流工程技术领域,尤其涉及一种物流系统切换建模与滤波方法。
背景技术
现代物流系统在城市基础建设中具有重要的作用。物流系统一般要求两个或两个以上的物流功能单元,这些物流功能单元有机结合在一起,进而完成物流服务的目的。物流系统主要由采购、运输、储存、流通加工、装卸、搬运、包装、销售、物流信息处理等物流环节。其中所需的劳务、设备、材料、资源等要素,这需要由外部环境向系统提供。物流系统就是要求在一定的时间和空间里,把需输送的物料比如快递以及相关的设备、或者是输送工具、人员以及通信联系等若干相互制约的动态要素构成的具有特定功能的部件构成一个有机整体。物流系统的成功要素是使物流系统整体优化以及合理化,并服从或改善社会大系统的环境。
随着我国物流系统建设不断取得进步,快递行业蓬勃发展,越来越多的人选择了网上购物,极大的方便了人们的日常生活,面对着大量的网上购物,物流系统也面临着极大的压力,为了迎合广大人民的需求,物流的建设也在不断增加,不管是城市还是乡村,出现了越来越多的快递站点,但是快递站点的增加并没有解决快递无法及时到站的情况,快递往往会因为各种原因而不能及时到达最终的客户手中,比如客户下单之后,快递从商家出货到快递仓库中,再经过分拣、打包、分类、运输等基本流程,这其中的每一个环节都会有各种可能性发生,比如遭遇恶劣的天气、交通阻塞、人员问题等导致快递无法及时到达到收货人手中。还有一些其它的情况,比如双十一、双十二这种网购狂欢日,一天成交的订单量极其庞大,如果处理不得当,甚至会导致物流系统瘫痪,在这种情况下,如何高效的运输快递到收货人显得尤为重要。网络购物无法离开物流系统,拥有一个好的物流系统就意味着拥有了市场,客户会优先选择服务更好的、收发货速度更快的物流公司。
发明内容
本发明目的在于提供一种物流系统切换建模与滤波方法,以解决快递爆仓导致的物流系统瘫痪的技术问题。
为了能够设计出令客户满意的物流系统,提出一种事件触发机制,当满足事件触发条件时来调动各个物流站点来分担不同的物流任务。考虑到快递出货时的数量是非负的,而且需要经过多个处理流程,因此使用正切换系统进行建模,为了防止类似双十一这种网购狂欢日带来的庞大的快递数量问题,采取饱和控制。为了更好的表明物流系统的高效特点,设计了一种事件触发异步滤波器来估计快递及时到达的数量。本发明的一种物流系统切换建模与滤波方法的具体技术方案如下:
一种物流系统切换建模与滤波方法,包括如下步骤:
步骤1:构造现代物流系统的状态空间模型;
步骤2:建立切换信号满足的切换条件;
步骤3:建立现代物流系统的事件触发条件;
步骤4:建立事件触发异步滤波器模型;
步骤5:引入饱和函数模型;
步骤6:构造现代物流系统的增广系统模型;
步骤7:设计针对现代物流系统的事件触发异步滤波器。
进一步地,所述步骤1的现代物流系统的状态空间模型如下:
Figure BDA0003204900340000031
y(t)=Cσ(t)x(t)+Dσ(t)w(t)
z(t)=Eσ(t)x(t)+Fσ(t)w(t)
其中,x(t)=[x1(t),x2(t),...,xn(t)]T∈Rn为时刻t运输到始发转运中心的快递数量,xj(t)为时刻t第j个始发转运中心上收到的快递数,其中1≤j≤n,n代表转运中心的个数;y(t)∈Rm为时刻t达到收货人的快递数量,m代表y(t)的维数,z(t)∈Rs是对y(t)的估计,s代表z(t)的维数,
Figure BDA0003204900340000032
是不可测的造成快递延误的外部扰动因素;σ(t)表示一个切换序列,在一个有限集
Figure BDA0003204900340000038
内取值,对于每一个切换序列0≤t0≤t1≤...它是右连续的;
Figure BDA0003204900340000037
表示对向量x(t)求导数运算;饱和函数sat(u)=(sat(u1),sat(u2),...sat(um))T,其中sat(ui)=sgn(ui)min{|ui,1|},i=1,2,...m,其中ui表示向量u的第i个分量,sgn是符号函数,min是最小值函数;Aσ(t),Bσ(t),Cσ(t),Dσ(t),Eσ(t),Fσ(t)为已知的系统矩阵;对于连续系统,假设矩阵Aσ(t)是Metzler矩阵,
Figure BDA0003204900340000033
Figure BDA0003204900340000034
Rn,
Figure BDA0003204900340000035
N+,Rn×n分别表示n维向量、n维非负向量、正整数和n×n维欧氏矩阵空间。
进一步地,所述步骤2切换信号σ(t)满足的切换条件如下:
Figure BDA0003204900340000036
其中,Nσ(t1,t2)表示时刻t1与时刻t2之间的切换次数,τ表示驻留时间,N0表示抖动界。
进一步地,所述步骤3的事件触发条件构造形式如下:
||ey(t)‖1>β||y(t)‖1
其中,常量0<β<1,ey(t)是测量的快递数量误差,
Figure BDA0003204900340000041
其中
Figure BDA00032049003400000414
Figure BDA00032049003400000413
表示自然数,它表示时刻
Figure BDA0003204900340000044
到达收货人的测量的快递数量,y(t)为时刻t到达收货人的实际快递数量。
进一步地,所述步骤4的事件触发异步滤波器模型的结构形式如下:
Figure BDA0003204900340000045
Figure BDA0003204900340000046
其中,xf(t)表示滤波器的状态信号,zf(t)表示对滤波器对输出信号y(t)的估计,Afi,Bfi,Efi,Ffi是要设计的滤波器矩阵。
进一步地,所述步骤5的饱和函数模型如下:
Figure BDA0003204900340000047
其中向量w∈Rn,k∈Rn,||k||≤1,Dl∈Rn×n是元素为0或1的对角矩阵,
Figure BDA0003204900340000048
令H∈Rn×n满足||H||≤1则有
Figure BDA0003204900340000049
其中矩阵A∈Rn×n
进一步地,所述步骤6的增广系统模型如下:
Figure BDA00032049003400000410
Figure BDA00032049003400000411
其中,
Figure BDA00032049003400000412
xe(t)=xf(t)-x(t),e(t)=zf(t)-z(t),
Figure BDA0003204900340000051
Figure BDA0003204900340000052
Figure BDA0003204900340000053
l是用来累加的中间变量,
Figure BDA0003204900340000054
是一个常数。
进一步地,所述步骤7包括如下具体步骤:
步骤7.1:设计的事件触发异步滤波器系统矩阵如下:
Figure BDA0003204900340000055
Figure BDA0003204900340000056
其中,α>0,Rn向量
Figure BDA0003204900340000057
ι,s是设计滤波器的中间变量,1n表示元素全为1的n维向量,
Figure BDA0003204900340000058
表示第μ个元素为1,其余元素为0的n维向量,I表示单位矩阵;
步骤7.2:设计常数α>0,β>0,γ>0,Rn向量
Figure BDA0003204900340000059
Figure BDA00032049003400000510
使得以下不等式:
Figure BDA0003204900340000061
Figure BDA0003204900340000062
Figure BDA0003204900340000063
Figure BDA0003204900340000064
Figure BDA0003204900340000065
Figure BDA0003204900340000066
Figure BDA0003204900340000067
Figure BDA0003204900340000068
Figure BDA0003204900340000069
Figure BDA00032049003400000610
Figure BDA00032049003400000611
Figure BDA00032049003400000612
Figure BDA00032049003400000613
Figure BDA00032049003400000614
Figure BDA00032049003400000615
Figure BDA00032049003400000616
对于每一个μ1>0,μ2>0,以及任意的l=1,2,...2n,在步骤7.1设计的滤波器下成立,其中α是证明故障检测系统为正使用的中间变量,β为事件触发条件中的系数,γ是L1增益性能指标,p,q是物流系统或滤波器的模态,(p,q)表示此时第p个子系统和第q个滤波器模态运行,Mql表示Mq的第l个分量,Φ=I-β1m×m,Ψ=I+β1m×m,1m×m是m行m列的全1矩阵,I为单位矩阵;
步骤7.3:设计正切换系统的异步切换律如下:
Figure BDA00032049003400000617
Figure BDA00032049003400000618
其中,κ-(t0,t)表示切换系统同步的总时间,κ+(t0,t)表示切换系统异步的总时间,τa,
Figure BDA0003204900340000071
表示切换子系统的驻留时间,Δm表示tm时刻异步的持续时间;
步骤7.4:根据步骤3、步骤5、步骤7.1和步骤7.2的前八个条件,得到保证故障检测系统为正的条件:
当t∈[tl,tll)
Figure BDA0003204900340000072
Figure BDA0003204900340000073
其中,
Figure BDA0003204900340000074
当t∈[tll,tl+1)
Figure BDA0003204900340000075
Figure BDA0003204900340000076
其中,
Figure BDA0003204900340000077
步骤7.5:考虑外部各种不确定因素对现代物流系统快速运输效率的影响,考虑如下约束性能:
Figure BDA0003204900340000078
其中,
Figure BDA0003204900340000079
ρ是正常数;
步骤7.6:根据步骤3、步骤7.1得到保证故障检测系统随机稳定的条件:
当t∈[tl,tll)
Figure BDA00032049003400000710
Figure BDA00032049003400000711
其中,
Figure BDA0003204900340000081
当t∈[tll,tl+1)
Figure BDA0003204900340000082
Figure BDA0003204900340000083
其中,
Figure BDA0003204900340000084
步骤7.7:设计多重李亚普诺夫函数:
Figure BDA0003204900340000085
其中,
Figure BDA0003204900340000086
定义Γ(t)=||e(t)||1-γ||w(t)||1
结合步骤7.1、步骤7.2的后八个条件、步骤7.6得到:
Figure BDA0003204900340000087
根据步骤7.2及步骤7.3中的条件,可以得到:
Figure BDA0003204900340000088
根据步骤7.7得到现代物流系统在设计的事件触发滤波器下是L1稳定的。
本发明的一种物流系统切换建模与滤波方法具有以下优点:本发明提出了一种现代物流系统中物流高效传输方法。本发明基于正切换系统模型、事件触发策略和基于异步滤波器的检测方法,针对物流系统中传输的快递量进行数据采集,提出了现代物流系统中物流高效传输方法,即便物流系统中有干扰因素或快递量过大,该方法都可以有效的提升物流系统的性能,提高快递传输的效率。
附图说明
图1为本发明的物流系统示意图。
具体实施方式
为了更好地了解本发明的目的、结构及功能,下面结合附图,对本发明一种物流系统切换建模与滤波方法做进一步详细的描述。
本发明通过对现代物流系统的输入输出数据进行采集描述实际的运作流程:
考虑一个物流系统的流程图,一个现代物流系统是主要由揽件端网点、始发转运中心、目的地转运中心、派件端网点组成,图1展示了快递全流程中各个揽件端网点、始发转运中心、目的地转运中心、派件端网点之间的联系。图1物流系统示意图中各个站点能够在快递运输过程中起到调节当前运输快递量的作用,对防止某一站点积攒大量快递而导致物流系统瘫痪有着至关重要的作用。我们知道在快递运输的过程中,往往会受到不同因素的干扰,比如恶劣的天气、交通拥挤等因素,导致快递无法及时到到达某个站点,或者受到网购狂欢日的影响而导致某一快递站点堆积大量的快递,造成了物流系统的局部瘫痪。比如图1中,当因为网购狂欢日到来时,始发转运中心就会堆积大量的快递,如果运输掉的快递数量一直小于来自揽件端网点的快递量,那么转运中心将会堆积越来越多的快递。对于其它快递站点也一样,如果运输前天气极其恶劣,就会造成快递的堆积,影响物流系统的正常运作。为了能够让物流系统更好的运转、为了快递能及时到达客户,必须要合理的安排转运中心的任务量,人口密度大的转运中心往往有着极大的转运压力,通过将转运分配给人口密度小或任务量较少的转运中心,这样当某个转运中心任务量过大时,任务量较小的转运中心将分担一些任务,从而避免了快递在站点堆积的问题,提高了快递的运输效率。采用一种事件触发策略来控制快递的运输,对时间要求严格的快递比如海鲜、水果等及时的派遣,按照这种策略依次对快递进行派遣,在一定程度上减缓了物流系统的压力。考虑到快递出货时的数量是非负的,而且需要经过多个处理流程,因此使用正切换系统进行建模,为了防止类似双十一这种网购狂欢日带来的庞大的快递数量问题,采取饱和控制。为了更好的表明物流系统的高效特点,设计了一种事件触发异步滤波器来估计快递及时到达的数量,进一步完善物流系统。
如图1所示,一种物流系统切换建模与滤波方法,以现代物流系统始发转运中心接收的快递数量和到收货人手中的快递数量为研究对象,以进入转运中心的快递数量为控制输入,以实际到达客户手中的快递数量为输出,来建立现代物流系统的动态模型,包括如下步骤:
步骤1:构造现代物流系统的状态空间模型:
Figure BDA0003204900340000101
y(t)=Cσ(t)x(t)+Dσ(t)w(t)
z(t)=Eσ(t)x(t)+Fσ(t)w(t)
其中,x(t)=[x1(t),x2(t),...,xn(t)]T∈Rn为时刻t运输到始发转运中心的快递数量,xj(t)为时刻t第j个始发转运中心上收到的快递数,其中1≤j≤n,n代表转运中心的个数。y(t)∈Rm为时刻t达到收货人的快递数量,m代表y(t)的维数,z(t)∈Rs是对y(t)的估计,s代表z(t)的维数,
Figure BDA0003204900340000102
是不可测的造成快递延误的外部扰动因素(例如高峰期车流量的突然增加而造成交通阻塞,恶劣的天气等)。σ(t)表示一个切换序列,在一个有限集S={1,2,...,J},J∈N+内取值,对于每一个切换序列0≤t0≤t1≤...它是右连续的。
Figure BDA0003204900340000103
表示对向量x(t)求导数运算。饱和函数sat(u)=(sat(u1),sat(u2),...sat(um))T,其中
Figure BDA00032049003400001111
其中ui表示向量u的第i个分量,sgn是符号函数,min是最小值函数。Aσ(t),Bσ(t),Cσ(t),Dσ(t),Eσ(t),Fσ(t)为已知的系统矩阵。对于连续系统,假设矩阵Aσ(t)是Metzler矩阵,
Figure BDA0003204900340000111
Figure BDA0003204900340000112
Rn,
Figure BDA0003204900340000113
N+,Rn×n分别表示n维向量、n维非负向量、正整数和n×n维欧氏矩阵空间。
步骤2:建立切换信号σ(t)满足的切换条件:
Figure BDA0003204900340000114
其中,Nσ(t1,t2)表示时刻t1与时刻t2之间的切换次数,τ表示驻留时间,N0表示抖动界。
步骤3:建立现代物流系统的事件触发条件,其构造形式如下:
||ey(t)‖1>β||y(t)‖1
其中,常量0<β<1,ey(t)是测量的快递数量误差,
Figure BDA0003204900340000115
其中
Figure BDA00032049003400001113
Figure BDA00032049003400001112
表示自然数,它表示时刻
Figure BDA0003204900340000118
到达收货人的测量的快递数量,y(t)为时刻t到达收货人的实际快递数量。
步骤4:建立事件触发异步滤波器模型,其结构形式如下:
Figure BDA0003204900340000119
Figure BDA00032049003400001110
其中,xf(t)表示滤波器的状态信号,zf(t)表示对滤波器对输出信号y(t)的估计,Afi,Bfi,Efi,Ffi是要设计的滤波器矩阵。
步骤5:引入饱和函数模型
Figure BDA0003204900340000121
其中向量w∈Rn,k∈Rn,||k||≤1,Dl∈Rn×n是元素为0或1的对角矩阵,
Figure BDA0003204900340000122
令H∈Rn×n满足||H||≤1则有
Figure BDA0003204900340000123
其中矩阵A∈Rn×n
步骤6:构造现代物流系统的增广系统模型:
Figure BDA0003204900340000124
Figure BDA0003204900340000125
其中,
Figure BDA0003204900340000126
xe(t)=xf(t)-x(t),e(t)=zf(t)-z(t),
Figure BDA0003204900340000127
Figure BDA0003204900340000128
Figure BDA0003204900340000129
l是用来累加的中间变量,
Figure BDA00032049003400001210
是一个常数。
步骤7:设计针对现代物流系统的事件触发异步滤波器:
7.1设计的事件触发异步滤波器系统矩阵如下:
Figure BDA00032049003400001211
Figure BDA00032049003400001212
其中,α>0,Rn向量
Figure BDA00032049003400001213
ι,s是设计滤波器的中间变量,1n表示元素全为1的n维向量,
Figure BDA0003204900340000131
表示第μ个元素为1,其余元素为0的n维向量,I表示单位矩阵。
7.2设计常数α>0,β>0,γ>0,Rn向量
Figure BDA0003204900340000132
Figure BDA0003204900340000133
使得以下不等式:
Figure BDA0003204900340000134
Figure BDA0003204900340000135
Figure BDA0003204900340000136
Figure BDA0003204900340000137
Figure BDA0003204900340000138
Figure BDA0003204900340000139
Figure BDA00032049003400001310
Figure BDA00032049003400001311
Figure BDA00032049003400001312
Figure BDA00032049003400001313
Figure BDA00032049003400001314
Figure BDA00032049003400001315
Figure BDA00032049003400001316
Figure BDA00032049003400001317
Figure BDA00032049003400001318
Figure BDA00032049003400001319
对于每一个μ1>0,μ2>0,以及任意的l=1,2,...2n,在步骤7.1设计的滤波器下成立,其中α是证明故障检测系统为正使用的中间变量,β为事件触发条件中的系数,γ是L1增益性能指标,p,q是物流系统或滤波器的模态,(p,q)表示此时第p个子系统和第q个滤波器模态运行,Mql表示Mq的第l个分量,Φ=I-β1m×m,Ψ=I+β1m×m,1m×m是m行m列的全1矩阵,I为单位矩阵。
7.3设计正切换系统的异步切换律如下:
Figure BDA0003204900340000141
Figure BDA0003204900340000142
其中,κ-(t0,t)表示切换系统同步的总时间,κ+(t0,t)表示切换系统异步的总时间,τa,
Figure BDA0003204900340000143
表示切换子系统的驻留时间,Δm表示tm时刻异步的持续时间。
7.4根据步骤3、步骤5、步骤7.1和步骤7.2的前八个条件,得到保证故障检测系统为正的条件:
当t∈[tl,tll)
Figure BDA0003204900340000144
Figure BDA0003204900340000145
其中,
Figure BDA0003204900340000146
当t∈[tll,tl+1)
Figure BDA0003204900340000147
Figure BDA0003204900340000148
其中,
Figure BDA0003204900340000149
7.5考虑外部各种不确定因素对现代物流系统快速运输效率的影响,考虑如下约束性能:
Figure BDA0003204900340000151
其中,
Figure BDA0003204900340000152
ρ是正常数。
7.6根据步骤3、步骤7.1得到保证故障检测系统随机稳定的条件:
当t∈[tl,tll)
Figure BDA0003204900340000153
Figure BDA0003204900340000154
其中,
Figure BDA0003204900340000155
当t∈[tll,tl+1)
Figure BDA0003204900340000156
Figure BDA0003204900340000157
其中,
Figure BDA0003204900340000158
7.7设计多重李亚普诺夫函数:
Figure BDA0003204900340000159
其中,
Figure BDA00032049003400001510
定义Γ(t)=||e(t)||1-γ||w(t)||1
结合步骤7.1、步骤7.2的后八个条件、步骤7.6得到:
Figure BDA0003204900340000161
根据步骤7.2及步骤7.3中的条件,可以得到:
Figure BDA0003204900340000162
根据步骤7.7说明现代物流系统在设计的事件触发滤波器下是L1稳定的。
可以理解,本发明是通过一些实施例进行描述的,本领域技术人员知悉的,在不脱离本发明的精神和范围的情况下,可以对这些特征和实施例进行各种改变或等效替换。另外,在本发明的教导下,可以对这些特征和实施例进行修改以适应具体的情况及材料而不会脱离本发明的精神和范围。因此,本发明不受此处所公开的具体实施例的限制,所有落入本申请的权利要求范围内的实施例都属于本发明所保护的范围内。

Claims (8)

1.一种物流系统切换建模与滤波方法,其特征在于,包括如下步骤:
步骤1:构造现代物流系统的状态空间模型;
步骤2:建立切换信号满足的切换条件;
步骤3:建立现代物流系统的事件触发条件;
步骤4:建立事件触发异步滤波器模型;
步骤5:引入饱和函数模型;
步骤6:构造现代物流系统的增广系统模型;
步骤7:设计针对现代物流系统的事件触发异步滤波器。
2.根据权利要求1所述的物流系统切换建模与滤波方法,其特征在于,所述步骤1的现代物流系统的状态空间模型如下:
Figure FDA0003204900330000013
y(t)=Cσ(t)x(t)+Dσ(t)w(t)
z(t)=Eσ(t)x(t)+Fσ(t)w(t)
其中,x(t)=[x1(t),x2(t),...,xn(t)]T∈Rn为时刻t运输到始发转运中心的快递数量,xj(t)为时刻t第j个始发转运中心上收到的快递数,其中1≤j≤n,n代表转运中心的个数;y(t)∈Rm为时刻t达到收货人的快递数量,m代表y(t)的维数,z(t)∈Rs是对y(t)的估计,s代表z(t)的维数,
Figure FDA0003204900330000011
是不可测的造成快递延误的外部扰动因素;σ(t)表示一个切换序列,在一个有限集
Figure FDA0003204900330000014
内取值,对于每一个切换序列0≤t0≤t1≤...它是右连续的;
Figure FDA0003204900330000012
表示对向量x(t)求导数运算;饱和函数sat(u)=(sat(u1),sat(u2),...sat(um))T,其中sat(ui)=sgn(ui)min{|ui,1|},i=1,2,...m,其中ui表示向量u的第i个分量,sgn是符号函数,min是最小值函数;Aσ(t),Bσ(t),Cσ(t),Dσ(t),Eσ(t),Fσ(t)为已知的系统矩阵;对于连续系统,假设矩阵Aσ(t)是Metzler矩阵,
Figure FDA0003204900330000015
Figure FDA0003204900330000029
Rn,
Figure FDA0003204900330000021
N+,Rn×n分别表示n维向量、n维非负向量、正整数和n×n维欧氏矩阵空间。
3.根据权利要求2所述的物流系统切换建模与滤波方法,其特征在于,所述步骤2切换信号σ(t)满足的切换条件如下:
Figure FDA0003204900330000022
其中,Nσ(t1,t2)表示时刻t1与时刻t2之间的切换次数,τ表示驻留时间,N0表示抖动界。
4.根据权利要求3所述的物流系统切换建模与滤波方法,其特征在于,所述步骤3的事件触发条件构造形式如下:
||ey(t)‖1>β||y(t)‖1
其中,常量0<β<1,ey(t)是测量的快递数量误差,
Figure FDA0003204900330000023
其中
Figure FDA0003204900330000024
Figure FDA0003204900330000025
表示自然数,它表示时刻
Figure FDA0003204900330000026
到达收货人的测量的快递数量,y(t)为时刻t到达收货人的实际快递数量。
5.根据权利要求4所述的物流系统切换建模与滤波方法,其特征在于,所述步骤4的事件触发异步滤波器模型的结构形式如下:
Figure FDA0003204900330000027
Figure FDA0003204900330000028
其中,xf(t)表示滤波器的状态信号,zf(t)表示对滤波器对输出信号y(t)的估计,Afi,Bfi,Efi,Ffi是要设计的滤波器矩阵。
6.根据权利要求5所述的物流系统切换建模与滤波方法,其特征在于,所述步骤5的饱和函数模型如下:
Figure FDA0003204900330000031
其中向量w∈Rn,k∈Rn,||k||≤1,Dl∈Rn×n是元素为0或1的对角矩阵,
Figure FDA0003204900330000032
令H∈Rn×n满足||H||≤1则有
Figure FDA0003204900330000033
其中矩阵A∈Rn×n
7.根据权利要求6所述的物流系统切换建模与滤波方法,其特征在于,所述步骤6的增广系统模型如下:
Figure FDA0003204900330000034
Figure FDA0003204900330000035
其中,
Figure FDA0003204900330000036
xe(t)=xf(t)-x(t),e(t)=zf(t)-z(t),
Figure FDA0003204900330000037
Figure FDA0003204900330000038
Figure FDA0003204900330000039
l是用来累加的中间变量,
Figure FDA00032049003300000310
是一个常数。
8.根据权利要求7所述的物流系统切换建模与滤波方法,其特征在于,所述步骤7包括如下具体步骤:
步骤7.1:设计的事件触发异步滤波器系统矩阵如下:
Figure FDA00032049003300000311
Figure FDA00032049003300000312
其中,α>0,Rn向量
Figure FDA00032049003300000418
ι,s是设计滤波器的中间变量,1n表示元素全为1的n维向量,
Figure FDA0003204900330000041
表示第μ个元素为1,其余元素为0的n维向量,I表示单位矩阵;
步骤7.2:设计常数α>0,β>0,γ>0,Rn向量
Figure FDA00032049003300000419
Figure FDA00032049003300000420
使得以下不等式:
Figure FDA0003204900330000042
Figure FDA0003204900330000043
Figure FDA0003204900330000044
Figure FDA0003204900330000045
Figure FDA0003204900330000046
Figure FDA0003204900330000047
Figure FDA0003204900330000048
Figure FDA0003204900330000049
Figure FDA00032049003300000410
Figure FDA00032049003300000411
Figure FDA00032049003300000412
Figure FDA00032049003300000413
Figure FDA00032049003300000414
Figure FDA00032049003300000415
Figure FDA00032049003300000416
Figure FDA00032049003300000417
对于每一个μ1>0,μ2>0,以及任意的l=1,2,...2n,在步骤7.1设计的滤波器下成立,其中α是证明故障检测系统为正使用的中间变量,β为事件触发条件中的系数,γ是L1增益性能指标,p,q是物流系统或滤波器的模态,(p,q)表示此时第p个子系统和第q个滤波器模态运行,Mql表示Mq的第l个分量,Φ=I-β1m×m,Ψ=I+β1m×m,1m×m是m行m列的全1矩阵,I为单位矩阵;
步骤7.3:设计正切换系统的异步切换律如下:
Figure FDA0003204900330000051
Figure FDA0003204900330000052
其中,κ-(t0,t)表示切换系统同步的总时间,κ+(t0,t)表示切换系统异步的总时间,τa,
Figure FDA0003204900330000053
表示切换子系统的驻留时间,Δm表示tm时刻异步的持续时间;
步骤7.4:根据步骤3、步骤5、步骤7.1和步骤7.2的前八个条件,得到保证故障检测系统为正的条件:
当t∈[tl,tll)
Figure FDA0003204900330000054
Figure FDA0003204900330000055
其中,
Figure FDA0003204900330000056
当t∈[tll,tl+1)
Figure FDA0003204900330000057
Figure FDA0003204900330000058
其中,
Figure FDA0003204900330000059
步骤7.5:考虑外部各种不确定因素对现代物流系统快速运输效率的影响,考虑如下约束性能:
Figure FDA0003204900330000061
其中,
Figure FDA0003204900330000062
ρ是正常数;
步骤7.6:根据步骤3、步骤7.1得到保证故障检测系统随机稳定的条件:
当t∈[tl,tll)
Figure FDA0003204900330000063
Figure FDA0003204900330000064
其中,
Figure FDA0003204900330000065
当t∈[tll,tl+1)
Figure FDA0003204900330000066
Figure FDA0003204900330000067
其中,
Figure FDA0003204900330000068
步骤7.7:设计多重李亚普诺夫函数:
Figure FDA0003204900330000069
其中,
Figure FDA00032049003300000610
定义Γ(t)=||e(t)||1-γ||w(t)||1
结合步骤7.1、步骤7.2的后八个条件、步骤7.6得到:
Figure FDA0003204900330000071
根据步骤7.2及步骤7.3中的条件,可以得到:
Figure FDA0003204900330000072
根据步骤7.7得到现代物流系统在设计的事件触发滤波器下是L1稳定的。
CN202110914386.7A 2021-08-10 2021-08-10 一种物流系统切换建模与滤波方法 Pending CN113657827A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110914386.7A CN113657827A (zh) 2021-08-10 2021-08-10 一种物流系统切换建模与滤波方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110914386.7A CN113657827A (zh) 2021-08-10 2021-08-10 一种物流系统切换建模与滤波方法

Publications (1)

Publication Number Publication Date
CN113657827A true CN113657827A (zh) 2021-11-16

Family

ID=78479305

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110914386.7A Pending CN113657827A (zh) 2021-08-10 2021-08-10 一种物流系统切换建模与滤波方法

Country Status (1)

Country Link
CN (1) CN113657827A (zh)

Similar Documents

Publication Publication Date Title
Zhu et al. A flexsim-based optimization for the operation process of cold-chain logistics distribution centre
Zhang et al. Multi-size container transportation by truck: modeling and optimization
CN104008475A (zh) 一种商品出库方法
CA2601866A1 (en) Next generation visibility package tracking
CN110140139B (zh) 用于在自动分配系统中对载货进行排序的方法
Walha et al. A rail-road PI-hub allocation problems: model and heuristic
CN110182530A (zh) 一种仓库管理方法和装置
US7084365B2 (en) Methods and apparatus for a non-discrete sortation process
Walha et al. The cross docking under uncertainty: State of the art
CN111552910B (zh) 在收集收集器上负载期间可减少无序的对自动化分配系统中负载进行排序的方法
Johnson et al. Performance analysis of split-case sorting systems
CN111047102B (zh) 基于精英驱动粒子群算法的快递配送路线优化方法
CN114926118A (zh) 车货装载信息的处理方法、装置、设备及存储介质
CN115215035A (zh) 一种接力拣选方法及系统
Moutaoukil et al. Modeling a logistics pooling strategy for Agri-Food SMEs
CN111784269A (zh) 物流转运管理系统
CN113657827A (zh) 一种物流系统切换建模与滤波方法
Thevenin et al. Dynamic scheduling of a picking robot with limited buffer and rejection: an industrial case study
CN110443445A (zh) 一种基于互联网的订单管理方法
CN117730311A (zh) 物流网络模拟方法和系统
CN111311019B (zh) 一种基于人口分布的物流配送系统
Liu Research on scheduling policy of automated warehouse system
CN113971539A (zh) 快递包裹拦截管理方法、装置、设备及存储介质
CN113962619A (zh) 运单的运输班次确定方法、装置、电子设备及存储介质
CN112044786A (zh) 阵列式无人物流分拣传送方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination