CN114664089B - 一种城市道路交通系统车流量的pi控制方法 - Google Patents

一种城市道路交通系统车流量的pi控制方法 Download PDF

Info

Publication number
CN114664089B
CN114664089B CN202210357833.8A CN202210357833A CN114664089B CN 114664089 B CN114664089 B CN 114664089B CN 202210357833 A CN202210357833 A CN 202210357833A CN 114664089 B CN114664089 B CN 114664089B
Authority
CN
China
Prior art keywords
urban road
road traffic
matrix
traffic
traffic flow
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.)
Active
Application number
CN202210357833.8A
Other languages
English (en)
Other versions
CN114664089A (zh
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 CN202210357833.8A priority Critical patent/CN114664089B/zh
Publication of CN114664089A publication Critical patent/CN114664089A/zh
Application granted granted Critical
Publication of CN114664089B publication Critical patent/CN114664089B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/0104Measuring and analyzing of parameters relative to traffic conditions
    • 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
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10Services
    • G06Q50/26Government or public services
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/0104Measuring and analyzing of parameters relative to traffic conditions
    • G08G1/0125Traffic data processing
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/0104Measuring and analyzing of parameters relative to traffic conditions
    • G08G1/0137Measuring and analyzing of parameters relative to traffic conditions for specific applications
    • G08G1/0145Measuring and analyzing of parameters relative to traffic conditions for specific applications for active traffic flow control
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/065Traffic control systems for road vehicles by counting the vehicles in a section of the road or in a parking area, i.e. comparing incoming count with outgoing count

Landscapes

  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Tourism & Hospitality (AREA)
  • Health & Medical Sciences (AREA)
  • Educational Administration (AREA)
  • Development Economics (AREA)
  • Economics (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • General Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Traffic Control Systems (AREA)

Abstract

本发明公开了一种城市道路交通系统车流量的PI控制方法,基于正马尔可夫跳变系统建立了城市道路交通系统的状态空间模型。借助多重线性共余正Lyapunov函数和线性规划方法设计了PI控制器。针对城市道路交通系统中的车流量进行数据采集,提出了一种对城市道路交通系统车流量进行控制的方法。此方法可以有效控制车流量大小,避免出行高峰期造成的道路堵塞以及道路安全隐患。本发明所建立的模型充分考虑了实际系统的正性、随机性、传感器造成的测量信息滞后性等特点,具有更高的应用价值。

Description

一种城市道路交通系统车流量的PI控制方法
技术领域
本发明属于自动化技术与现代控制领域,涉及一种城市道路交通系统车流量的建模以及一种PI控制方法。
背景技术
城市中的各个组成部分通过城市道路系统,将它们联结成一个相互协调、有机联系的整体。城市道路交通是城市社会、经济和物质结构的基本组成部分,是城市交通系统的主题。最近十几年,随着人们生活质量的提升,购买机动车辆的人群呈几何级数增长,显示出对城市道路交通需求的急剧增长,这也暴露出城市道路交通基础设施与交通需求之间的冲突和矛盾,出现了交通阻塞、事故增加、环境污染加剧等问题。因此,为提高城市道路交通系统的交通分流能力、避免出现交通堵塞的现象,必须对城市道路交通系统的流量进行科学调控,以保证交通安全和顺畅。城市道路同一般公路相比,主要特点是:道路交叉点多,区间段短,交通流速较低,通行能力较小,通行压力较大;道路上行人和公共交通车辆,机动车和非机动车等各种交通流相互交织,交通组织比较复杂;在交通安全和交通管理方面要求较高。由于城市道路交通系统具备以上几个特点,因此难以对实际系统建立精确模型,从而无法提出科学有效的控制车流量的策略。本发明主要是动态获取城市道路交通系统中车流量值,根据车流量实际情况来采取更为合理的管控方法,提出一种可靠的城市道路交通车流量控制方法,以对车流量的大小采取必要的控制措施,以达到城市道路交通安全运行的目的。
因为城市道路交通系统中车流量总是非负的,所以采用正系统建模更为精确。考虑到城市道路交通系统中同一路段的车流量并不是一成不变也不是按照既定规律变化的,而是随机发生变化的,采用正马尔可夫跳变系统建模与实际系统相符。在实际城市道路交通系统中,为避免某一时间段或某一路段车流量过大导致交通堵塞、交通事故等现象,需要车流量维持在一定的数值或者一定的范围内。传统的控制系统通常采用比例控制来达到预期目标,但由于检测设备故障、突发紧急事件等因素的影响,实际城市道路交通系统车流量的监测可能会存在一定的滞后性,使系统误差增大。因此,仅用比例控制来控制实际系统无法实现对车流量的精确控制。为了解决这一问题,本发明提出了基于PI控制器的控制策略。PI控制器(比例-积分控制器)是通过设定比例单元P、积分单元I,对整个控制系统进行偏差调节,从而使被控变量的实际值与给定值一致。PI控制器将比例(P)控制和积分(I)控制结合起来,既有比例控制作用的控制及时,又有积分控制作用消除偏差的能力,而且易于设计、实现方便。
综上,本发明采用现代控制理论技术建立城市道路交通系统的状态空间模型,基于正马尔可夫跳变系统设计了PI控制器以及反馈增益矩阵,并分析其正性和稳定性,最终使系统中的车流量值维持在给定值上,从而可以有效避免交通堵塞、交通事故等现象的发生。
发明内容
本发明针对城市道路交通系统车流量的控制问题,建立了城市道路交通系统的正马尔可夫跳变系统模型,并提出了一种城市道路交通系统车流量的PI控制方法。
本发明方法的具体步骤包括如下:
步骤1、结合城市道路交通系统,建立状态空间模型;
步骤2、构建城市道路交通控制系统的控制律;
步骤3、设计事件触发PI控制器的积分部分;
步骤4、构造基于PI控制律的闭环控制系统;
步骤5、设计针对城市道路交通系统的PI控制器。
步骤1具体步骤如下:
1.1对城市道路交通系统车流量的输入输出数据量进行采集;
通过道路监测系统获取进入某个城区、开出某个城区的车流量数据。考虑到实际城市道路交通系统中车流量的非负性以及随机性,可以将实际系统抽象为一个正马尔科夫跳变系统。在实际系统中,由于传感器测量得到的信息存在一定的滞后性,导致测量得到的车流量数据不精确,为了解决这一问题,提出了PI控制策略。
1.2采集城市道路交通系统车流量数据,建立城市道路交通系统的状态空间模型,形式如下:
δx(t)=Ar(t)x(t)+Br(t)u(t),
y(t)=Cr(t)x(t),
其中,x(t)=[x1(t)x2(t)…xn(t)]T∈Rn00000000由t时刻城市道路交通系统的车流量组成,其中,xn(t)表示系统中第n个传感器测量得到的车流量数据。δx(t)表示对向量x(t)求导数运算。u(t)∈Rr表示t时刻城市道路交通系统车流量的控制信号。y(t)∈Rs表示t时刻从城市道路交通系统离开的实际车流量。函数r(t)是马尔科夫跳变信号即城市道路交通系统的跳变信号,表示[0,∞]到有限集S={1,2,···,N+}的映射。令r(t)=i,i∈S,表示系统处于第i个子系统状态,则系统矩阵可被记作Ai,Bi,Ci,Ai为第i个子系统系统的状态矩阵,Bi为第i个子系统系统的输入矩阵,Ci为第i个子系统系统输出矩阵。矩阵Ai是Metzler矩阵,Bi≥0,Ci≥0。Rn,N+,Rn×n分别表示n维向量、正整数和n×n维欧氏矩阵空间。[x1(t)x2(t)…xn(t)]T表示向量[x1(t)x2(t)…xn(t)]的转置。
1.3设计马尔科夫跳变信号r(t),其转移速率满足以下条件:
Figure BDA0003582577720000041
其中,转移率λij≥0表示跳变信号从子系统i到子系统j,
Figure BDA0003582577720000042
λii表示跳变信号从子系统i跳变到自身。
步骤2构建城市道路交通控制系统的控制律,其结构形式如下:
u(t)=KPx(t)+FIe(t),
其中,KP和FI分别是需要设计的第i个子系统的比例增益矩阵和积分增益矩阵,e(t)是PI控制器的积分部分。
步骤3、设计事件触发PI控制器的积分部分,其构建形式如下:
δe(t)=y(t)-αe(t),
其中,α为调优参数且α>0,δe(t)表示对向量e(t)求导数运算。
步骤4、构造基于PI控制律的闭环控制系统,具体如下:
结合步骤1.2和步骤3,得到城市道路交通系统的闭环控制系统:
Figure BDA0003582577720000043
Figure BDA0003582577720000044
则进一步表示为:
Figure BDA0003582577720000045
其中,
Figure BDA0003582577720000046
表示对向量/>
Figure BDA0003582577720000047
求导数运算,/>
Figure BDA0003582577720000048
表示/>
Figure BDA0003582577720000049
的系数矩阵,Is表示s维单位向量。
步骤5、设计针对城市道路交通系统的PI控制器。
5.1设计的PI控制器增益矩阵如下:
Figure BDA0003582577720000051
其中,piiil分别为由不同实数组成的n维向量,θiik分别为由不同实数组成的s维向量,l,k是中间变量。1r表示元素全为1的r维向量,
Figure BDA0003582577720000052
表示第l个元素为1,其余元素为0的n维向量,/>
Figure BDA0003582577720000053
表示第k个元素为1,其余元素为0的n维向量。
5.2设计常数
Figure BDA0003582577720000054
μ>0,存在Rn向量pi>0,εi>0,εil>0,Rs向量/>
Figure BDA0003582577720000055
θil>0,θi>0,使得以下不等式,对于任意的i,j∈S成立:
Figure BDA0003582577720000056
Figure BDA0003582577720000057
Figure BDA0003582577720000058
0<εil≤εi,l=1,2,...,r,
0<θik≤θi,k=1,2,...,r,
则在步骤2设计的PI控制律满足步骤5.1所设计的比例增益和积分增益下,步骤4得到的城市道路交通闭环控制系统是正的,并且是随机稳定的。其中,I表示单位矩阵。
5.3城市道路交通闭环控制系统的正性验证过程如下:
根据步骤5.2中保证系统正性的条件,得到:
Figure BDA0003582577720000059
即矩阵Ai+BiKP为Metzler矩阵。由Is是单位矩阵得到矩阵-αIs为Metzler矩阵。
根据步骤1.2,步骤5.1和步骤5.2,得到:
Ci≥0,BiFI≥0.
Figure BDA0003582577720000061
为Metzler矩阵。因此城市道路交通闭环控制系统是正的。
5.4城市道路交通闭环控制系统的随机稳定性验证过程如下:
构造一个共余正Lyapunov函数:
Figure BDA0003582577720000062
其中,
Figure BDA0003582577720000063
表示向量pi的转置,/>
Figure BDA0003582577720000064
表示向量/>
Figure BDA0003582577720000065
的转置。
得到:
Figure BDA0003582577720000066
其中,Γ表示弱无穷小算子。
5.5结合步骤1.3和步骤5.4,得到;
Figure BDA0003582577720000067
5.6根据幂级数展开公式,
Figure BDA0003582577720000068
的一阶泰勒展开估计为:
Figure BDA0003582577720000069
5.7结合步骤4,步骤5.5以及步骤5.6,得到以下式子:
Figure BDA00035825777200000610
进一步,得到:
Figure BDA00035825777200000611
5.8根据Dynkin公式,得到:
Figure BDA0003582577720000071
5.9结合步骤5.7和步骤5.8,得到以下不等式:
Figure BDA0003582577720000072
进一步可以得到:
Figure BDA0003582577720000073
/>
5.10结合步骤4,步骤5.1,步骤5.2和步骤5.9得到:
Figure BDA0003582577720000074
其中,
Figure BDA0003582577720000075
表示数学期望,dτ表示时间变量算子,/>
Figure BDA0003582577720000076
表示初始状态,r(0)表示初始转移概率。
5.11根据步骤5.2以及步骤5.10,得到:
Figure BDA0003582577720000077
根据步骤5.11说明城市道路闭环控制系统在设计PI控制器下是随机稳定的。
本发明有益效果如下:
本发明基于正马尔可夫跳变系统建立了城市道路交通系统的状态空间模型。借助多重线性共余正Lyapunov函数和线性规划方法设计了PI控制器。针对城市道路交通系统中的车流量进行数据采集,提出了一种对城市道路交通系统车流量进行控制的方法。此方法可以有效控制车流量大小,避免出行高峰期造成的道路堵塞以及道路安全隐患。本发明所建立的模型充分考虑了实际系统的正性、随机性、传感器造成的测量信息滞后性等特点,具有更高的应用价值。
附图说明
图1为某城市道路交通系统示意图;
图2PI控制框架示意图。
具体实施方式
下面结合附图对本发明做进一步说明。
如图1所示,以城市道路交通系统为研究对象,以进入交叉路口的车流量为输入,以传感器实际测量到的离开交叉路口的车流量为输出,建立了该城市道路交通系统的动态模型。图2是PI控制框架示意图。
一种城市道路交通车流量的PI控制方法,包括如下步骤:
步骤1、结合城市道路交通系统,建立状态空间模型。
1.1对城市道路交通系统车流量的输入输出数据量进行采集以描述实际的系统:
考虑某城市的道路交通系统如图1所示(见说明书附图),其中五个城市建筑标志代表A、B、C、D、E五个城区,黑色线段表示城区之间的交通道路,白色汽车模型表示公共交通工具,黄色汽车模型表示私家车或者出租车,此示意图展示了该城市道路交通系统车流量的输入输出关系。在实际系统中,可以通过道路监测系统获取进入某个城区、开出某个城区的车流量数据,一般的检测汽车流量的方法是在地面下埋入一个线圈通电,当汽车通过时,电流会发生改变,测电流比较简单准确。在上下班或者节假日的出行高峰期,车流量大幅度增加,流入某一区域的车流量过大容易造成交通拥堵的现象。因此,根据城区之间的车流量数据构建适当的数学模型进而及时对流入流出的车流量进行控制,使某一城区的车流量维持在一定的大小,对城市道路交通系统的安全运行具有重要意义。考虑到实际城市道路交通系统中车流量的非负性以及随机性,可以将实际系统抽象为一个正马尔科夫跳变系统。在实际系统中,由于传感器测量得到的信息存在一定的滞后性,导致测量得到的车流量数据不精确。为了解决这一问题,提出了PI控制策略。图2(见说明书附图)是PI控制框架图。给定一个希望车流量维持的数值,当系统中车流量的状态x(t)和输出y(t)产生的误差较大时,PI控制器启用,用来调节输入车流量的大小,进而使系统的给定值保持在一定的数值上。平时可认为系统状态处于子系统1,当下班高峰期或者节假日出行时,系统切换到子系统2,PI控制器也会切换到相应的控制器。
1.2采集城市道路交通系统车流量数据,建立城市道路交通系统的状态空间模型,形式如下:
δx(t)=Ar(t)x(t)+Br(t)u(t),
y(t)=Cr(t)x(t),
其中,x(t)=[x1(t)x2(t)…xn(t)]T∈Rn0为n维的车流量状态向量,,由t时刻城市道路交通系统的车流量组成,其中,xn(t)表示系统中第n个传感器测量得到的车流量数据。δx(t)表示对向量x(t)求导数运算。u(t)∈Rr表示t时刻城市道路交通系统车流量的控制信号。y(t)∈Rs表示t时刻从城市道路交通系统离开的实际车流量。函数r(t)是马尔科夫跳变信号即城市道路交通系统的跳变信号,表示[0,∞]到有限集S={1,2,···,N+}的映射。令r(t)=i,i∈S,表示系统处于第i个子系统状态,则系统矩阵可被记作Ai,Bi,Ci,Ai为第i个子系统系统的状态矩阵,Bi为第i个子系统系统的输入矩阵,Ci为第i个子系统系统输出矩阵。矩阵Ai是Metzler矩阵,Bi≥0,Ci≥0。Rn,N+,Rn×n分别表示n维向量、正整数和n×n维欧氏矩阵空间。[x1(t)x2(t)…xn(t)]T表示向量[x1(t)x2(t)…xn(t)]的转置。
1.3设计马尔科夫跳变信号r(t),其转移速率满足以下条件:
Figure BDA0003582577720000091
其中,转移率λij≥0表示跳变信号从子系统i到子系统j,
Figure BDA0003582577720000092
λii表示跳变信号从子系统i跳变到自身。
步骤2、构建城市道路交通控制系统的控制律,其结构形式如下:
u(t)=KPx(t)+FIe(t),
其中,KP和FI分别是需要设计的第i个子系统的比例增益矩阵和积分增益矩阵,e(t)是PI控制器的积分部分。
步骤3、设计事件触发PI控制器的积分部分,其构建形式如下:
δe(t)=y(t)-αe(t),
其中,α为调优参数且α>0,δe(t)表示对向量e(t)求导数运算。
步骤4、构造基于PI控制律的闭环控制系统,具体如下:
结合步骤1.2和步骤3,得到城市道路交通系统的闭环控制系统:
Figure BDA0003582577720000101
Figure BDA0003582577720000102
则进一步表示为:/>
Figure BDA0003582577720000103
其中,
Figure BDA0003582577720000104
表示对向量/>
Figure BDA0003582577720000105
求导数运算,/>
Figure BDA0003582577720000106
表示/>
Figure BDA0003582577720000107
的系数矩阵,Is表示s维单位向量。
步骤5、设计针对城市道路交通系统的PI控制器。
5.1设计的PI控制器增益矩阵如下:
Figure BDA0003582577720000108
其中,piiil分别为由不同实数组成的n维向量,θiik分别为由不同实数组成的s维向量,l,k是中间变量。1r表示元素全为1的r维向量,
Figure BDA0003582577720000109
表示第l个元素为1,其余元素为0的n维向量,/>
Figure BDA00035825777200001010
表示第k个元素为1,其余元素为0的n维向量。
5.2设计常数
Figure BDA00035825777200001011
μ>0,存在Rn向量pi>0,εi>0,εil>0,Rs向量/>
Figure BDA00035825777200001012
θi>0,θi>0,使得以下不等式,对于任意的i,j∈S成立:
Figure BDA0003582577720000111
Figure BDA0003582577720000112
Figure BDA0003582577720000113
0<εil≤εi,l=1,2,...,r,
0<θik≤θi,k=1,2,...,r,
则在步骤2设计的PI控制律满足步骤5.1所设计的比例增益和积分增益下,步骤4得到的城市道路交通闭环控制系统是正的,并且是随机稳定的。其中,I表示单位矩阵。
5.3城市道路交通闭环控制系统的正性验证过程如下:
根据步骤5.2中保证系统正性的条件,得到:
Figure BDA0003582577720000114
即矩阵Ai+BiKP为Metzler矩阵。由Is是单位矩阵得到矩阵-αIs为Metzler矩阵。
根据步骤1.2,步骤5.1和步骤5.2,得到:
Ci≥0,BiFI≥0.
Figure BDA0003582577720000115
为Metzler矩阵。因此城市道路交通闭环控制系统是正的。
5.4城市道路交通闭环控制系统的随机稳定性验证过程如下:
构造一个共余正Lyapunov函数:
Figure BDA0003582577720000116
/>
其中,
Figure BDA0003582577720000117
表示向量pi的转置,/>
Figure BDA0003582577720000118
表示向量/>
Figure BDA0003582577720000119
的转置。得到:
Figure BDA0003582577720000121
其中,Γ表示弱无穷小算子。
5.5结合步骤1.3和步骤5.4,得到;
Figure BDA0003582577720000122
5.6根据幂级数展开公式,
Figure BDA0003582577720000123
的一阶泰勒展开估计为:
Figure BDA0003582577720000124
5.7结合步骤4,步骤5.5以及步骤5.6,得到以下式子:
Figure BDA0003582577720000125
进一步,得到:
Figure BDA0003582577720000126
5.8根据Dynkin公式,得到:
Figure BDA0003582577720000127
5.9结合步骤5.7和步骤5.8,得到以下不等式:
Figure BDA0003582577720000128
进一步可以得到:
Figure BDA0003582577720000129
5.10结合步骤4,步骤5.1,步骤5.2和步骤5.9得到:
Figure BDA0003582577720000131
其中,
Figure BDA0003582577720000132
表示数学期望,dτ表示时间变量算子,/>
Figure BDA0003582577720000133
表示初始状态,r(0)表示初始转移概率。
5.11根据步骤5.2以及步骤5.10,得到:
Figure BDA0003582577720000134
根据步骤5.11说明城市道路闭环控制系统在设计PI控制器下是随机稳定的。

Claims (1)

1.一种城市道路交通系统车流量的PI控制方法,其特征在于,步骤如下:
步骤1、结合城市道路交通系统,建立状态空间模型;
步骤2、构建城市道路交通控制系统的控制律;
步骤3、设计事件触发PI控制器的积分部分;
步骤4、构造基于PI控制律的闭环控制系统;
步骤5、设计针对城市道路交通系统的PI控制器;
步骤1具体步骤如下:
1.1对城市道路交通系统车流量的输入输出数据量进行采集;
通过道路监测系统获取进入某个城区、开出某个城区的车流量数据;考虑到实际城市道路交通系统中车流量的非负性以及随机性,可以将实际系统抽象为一个正马尔科夫跳变系统;在实际系统中,由于传感器测量得到的信息存在一定的滞后性,导致测量得到的车流量数据不精确,为了解决这一问题,提出了PI控制策略;
1.2采集城市道路交通系统车流量数据,建立城市道路交通系统的状态空间模型,形式如下:
δx(t)=Ar(t)x(t)+Br(t)u(t),
y(t)=Cr(t)x(t),
其中,x(t)=[x1(t) x2(t)…xn(t)]T∈Rn00000000由t时刻城市道路交通系统的车流量组成,其中,xn(t)表示系统中第n个传感器测量得到的车流量数据;δx(t)表示对向量x(t)求导数运算;u(t)∈Rr表示t时刻城市道路交通系统车流量的控制信号;y(t)∈Rs表示t时刻从城市道路交通系统离开的实际车流量;函数r(t)是马尔科夫跳变信号即城市道路交通系统的跳变信号,表示[0,∞]到有限集S={1,2,···,N+}的映射;令r(t)=i,i∈S,表示系统处于第i个子系统状态,则系统矩阵可被记作Ai,Bi,Ci,Ai为第i个子系统系统的状态矩阵,Bi为第i个子系统系统的输入矩阵,Ci为第i个子系统系统输出矩阵;矩阵Ai是Metzler矩阵,Bi≥0,Ci≥0;Rn,N+,Rn×n分别表示n维向量、正整数和n×n维欧氏矩阵空间;[x1(t) x2(t)…xn(t)]T表示向量[x1(t) x2(t)…xn(t)]的转置;
1.3设计马尔科夫跳变信号r(t),其转移速率满足以下条件:
Figure FDA0003990438300000021
其中,转移率λij≥0表示跳变信号从子系统i到子系统j,
Figure FDA0003990438300000022
λii表示跳变信号从子系统i跳变到自身;
步骤2构建城市道路交通控制系统的控制律,其结构形式如下:
u(t)=KPx(t)+FIe(t),
其中,KP和FI分别是需要设计的第i个子系统的比例增益矩阵和积分增益矩阵,e(t)是PI控制器的积分部分;
步骤3、设计事件触发PI控制器的积分部分,其构建形式如下:
δe(t)=y(t)-αe(t),
其中,α为调优参数且α>0,δe(t)表示对向量e(t)求导数运算;
步骤4、构造基于PI控制律的闭环控制系统,具体如下:
结合步骤1.2和步骤3,得到城市道路交通系统的闭环控制系统:
Figure FDA0003990438300000023
Figure FDA0003990438300000024
则进一步表示为:
Figure FDA0003990438300000025
其中,
Figure FDA0003990438300000026
表示对向量
Figure FDA0003990438300000027
求导数运算,
Figure FDA0003990438300000028
表示
Figure FDA0003990438300000029
的系数矩阵,Is表示s维单位向量;
步骤5、设计针对城市道路交通系统的PI控制器;
5.1设计的PI控制器增益矩阵如下:
Figure FDA0003990438300000031
其中,piiil分别为由不同实数组成的n维向量,θiik分别为由不同实数组成的s维向量,l,k是中间变量;1r表示元素全为1的r维向量,
Figure FDA0003990438300000032
表示第l个元素为1,其余元素为0的n维向量,
Figure FDA0003990438300000033
表示第k个元素为1,其余元素为0的n维向量;
5.2设计常数
Figure FDA0003990438300000034
存在Rn向量pi>0,εi>0,εil>0,Rs向量
Figure FDA0003990438300000035
使得以下不等式,对于任意的i,j∈S成立:
Figure FDA0003990438300000036
Figure FDA0003990438300000037
Figure FDA0003990438300000038
0<εil≤εi,l=1,2,...,r,
0<θik≤θi,k=1,2,...,r,
则在步骤2设计的PI控制律满足步骤5.1所设计的比例增益和积分增益下,步骤4得到的城市道路交通闭环控制系统是正的,并且是随机稳定的;其中,I表示单位矩阵;
5.3城市道路交通闭环控制系统的正性验证过程如下:
根据步骤5.2中保证系统正性的条件,得到:
Figure FDA0003990438300000039
即矩阵Ai+BiKP为Metzler矩阵;由Is是单位矩阵得到矩阵-αIs为Metzler矩阵;
根据步骤1.2,步骤5.1和步骤5.2,得到:
Ci≥0,BiFI≥0,
Figure FDA00039904383000000412
为Metzler矩阵;因此城市道路交通闭环控制系统是正的;
5.4城市道路交通闭环控制系统的随机稳定性验证过程如下:
构造一个共余正Lyapunov函数:
Figure FDA0003990438300000041
其中,
Figure FDA0003990438300000042
Figure FDA0003990438300000043
表示向量pi的转置,
Figure FDA0003990438300000044
表示向量
Figure FDA0003990438300000045
的转置;
得到:
Figure FDA0003990438300000046
其中,Γ表示弱无穷小算子;
5.5结合步骤1.3和步骤5.4,得到;
Figure FDA0003990438300000047
5.6根据幂级数展开公式,
Figure FDA0003990438300000048
的一阶泰勒展开估计为:
Figure FDA0003990438300000049
5.7结合步骤4,步骤5.5以及步骤5.6,得到以下式子:
Figure FDA00039904383000000410
进一步,得到:
Figure FDA00039904383000000411
5.8根据Dynkin公式,得到:
Figure FDA0003990438300000051
5.9结合步骤5.7和步骤5.8,得到以下不等式:
Figure FDA0003990438300000052
进一步可以得到:
Figure FDA0003990438300000053
5.10结合步骤4,步骤5.1,步骤5.2和步骤5.9得到:
Figure FDA0003990438300000054
其中,
Figure FDA0003990438300000055
表示数学期望,dτ表示时间变量算子,
Figure FDA0003990438300000056
表示初始状态,r(0)表示初始转移概率;
5.11根据步骤5.2以及步骤5.10,得到:
Figure FDA0003990438300000057
根据步骤5.11说明城市道路闭环控制系统在设计PI控制器下是随机稳定的。
CN202210357833.8A 2022-04-06 2022-04-06 一种城市道路交通系统车流量的pi控制方法 Active CN114664089B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210357833.8A CN114664089B (zh) 2022-04-06 2022-04-06 一种城市道路交通系统车流量的pi控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210357833.8A CN114664089B (zh) 2022-04-06 2022-04-06 一种城市道路交通系统车流量的pi控制方法

Publications (2)

Publication Number Publication Date
CN114664089A CN114664089A (zh) 2022-06-24
CN114664089B true CN114664089B (zh) 2023-03-24

Family

ID=82034919

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210357833.8A Active CN114664089B (zh) 2022-04-06 2022-04-06 一种城市道路交通系统车流量的pi控制方法

Country Status (1)

Country Link
CN (1) CN114664089B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116740922B (zh) * 2023-05-08 2024-04-05 海南大学 一种基于模糊观测协议的智慧交通系统的控制方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103700274A (zh) * 2014-01-17 2014-04-02 重庆大学 一种交通流检测诱导方法
CA3104842A1 (en) * 2020-01-03 2021-07-03 Battelle Memorial Institute Blockchain applicability framework and cybersecurity vulnerability mitigation framework

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102375442A (zh) * 2010-08-23 2012-03-14 同济大学 一种对混杂非线性系统的实时在线控制系统及控制方法
CN105159090A (zh) * 2015-09-24 2015-12-16 盐城工学院 基于马尔可夫模型的随机系统的pi跟踪控制器设计方法
CN106842947B (zh) * 2017-02-28 2019-11-05 杭州电子科技大学 一种城市排水管网的安全运行控制方法
CN109828493B (zh) * 2019-01-11 2020-04-21 杭州电子科技大学 一种城市水务管网系统供水故障的控制方法
CN109901395B (zh) * 2019-03-22 2022-03-11 杭州电子科技大学 异步系统的自适应容错控制方法
CN110020475B (zh) * 2019-04-03 2023-10-10 北京工业大学 一种交通流预测的马尔科夫粒子滤波方法
CN110187641B (zh) * 2019-07-12 2021-12-14 杭州电子科技大学 一种城市水务管网供水系统在外加扰动输入下的控制方法
CN110610019A (zh) * 2019-07-29 2019-12-24 南京理工大学 转移概率部分未知的马尔科夫跳变系统的动力学分析方法
CN112907417B (zh) * 2020-12-31 2023-12-08 杭州电子科技大学 基于正系统的城市交通建模和事件触发非脆弱控制方法
CN113033976B (zh) * 2021-03-10 2022-05-06 杭州电子科技大学 基于事件触发机制的城市道路系统的可靠滤波设计方法
CN113126500B (zh) * 2021-04-20 2022-06-07 河北工业大学 一种交通系统的非脆弱同步控制方法及装置
CN113238485A (zh) * 2021-05-13 2021-08-10 杭州电子科技大学 一种空中交通流量输入饱和的控制方法
CN113393691B (zh) * 2021-05-31 2022-05-03 杭州电子科技大学 一种城市道路系统交通信号灯故障检测方法
CN113378376B (zh) * 2021-06-08 2024-02-13 杭州电子科技大学 一种城市轨道交通车站客流量的混合触发控制方法
CN113641102B (zh) * 2021-08-10 2024-01-30 杭州电子科技大学 城市轨道交通车站正系统建模与混杂滤波器设计方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103700274A (zh) * 2014-01-17 2014-04-02 重庆大学 一种交通流检测诱导方法
CA3104842A1 (en) * 2020-01-03 2021-07-03 Battelle Memorial Institute Blockchain applicability framework and cybersecurity vulnerability mitigation framework

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
输入约束下高阶智能体均方一致性研究;输入约束下高阶智能体均方一致性研究;《电子科技大学学报》;20201130;全文 *

Also Published As

Publication number Publication date
CN114664089A (zh) 2022-06-24

Similar Documents

Publication Publication Date Title
Xu et al. Modeling relationship between truck fuel consumption and driving behavior using data from internet of vehicles
Gao et al. Reinforcement-learning-based cooperative adaptive cruise control of buses in the Lincoln tunnel corridor with time-varying topology
Yang et al. Eco-driving system for connected automated vehicles: Multi-objective trajectory optimization
Zhang et al. An improved k-nearest neighbor model for short-term traffic flow prediction
CN102629418B (zh) 基于模糊卡尔曼滤波的交通流预测方法
Bao et al. An improved deep belief network for traffic prediction considering weather factors
CN107085943B (zh) 一种道路旅行时间短期预测方法和系统
Jeong et al. Prediction model of bus arrival time for real-time applications
CN106503840A (zh) 停车场可用车位预测方法及系统
DE102015208380A1 (de) Fahrzeugenergie-Handhabungsvorrichtung
CN102136190A (zh) 城市公交客运事件应急响应调度管理系统及方法
CN113033976B (zh) 基于事件触发机制的城市道路系统的可靠滤波设计方法
CN109493566A (zh) 一种基于gps数据的疲劳驾驶预测方法
CN114664089B (zh) 一种城市道路交通系统车流量的pi控制方法
CN103942953A (zh) 一种基于浮动车数据的城市路网动态交通拥挤预测方法
CN106019938B (zh) 基于数据驱动的acc系统离散二阶滑模控制系统的控制方法
CN108427273A (zh) 一种减少交通拥堵现象的反馈控制设计方法
Zhi-Peng et al. An improved adaptive exponential smoothing model for short-term travel time forecasting of urban arterial street
Suh et al. Vehicle speed prediction for connected and autonomous vehicles using communication and perception
Bae et al. Multicontextual machine-learning approach to modeling traffic impact of urban highway work zones
CN114418606B (zh) 基于时空卷积网络的网约车订单需求预测方法
Yufang et al. Prediction of vehicle energy consumption on a planned route based on speed features forecasting
CN108364490A (zh) 城市公路交通系统车辆运行调控方法
CN111754776A (zh) 一种基于加权分数阶灰色模型的交通流预测方法及系统
DE102020202650A1 (de) Computerimplementiertes Verfahren zur Routenberechnung für ein autonom fahrendes Kraftfahrzeug, Verfahren zum Fahren eines autonom fahrenden Kraftfahrzeugs, Computerprogrammprodukt sowie Kraftfahrzeug

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
GR01 Patent grant
GR01 Patent grant