CN115358469B - 一种获得城市交通集成传感设施的铺设最优方案的方法 - Google Patents

一种获得城市交通集成传感设施的铺设最优方案的方法 Download PDF

Info

Publication number
CN115358469B
CN115358469B CN202211012278.1A CN202211012278A CN115358469B CN 115358469 B CN115358469 B CN 115358469B CN 202211012278 A CN202211012278 A CN 202211012278A CN 115358469 B CN115358469 B CN 115358469B
Authority
CN
China
Prior art keywords
integrated sensing
facilities
laying
integrated
sensing
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
CN202211012278.1A
Other languages
English (en)
Other versions
CN115358469A (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.)
Beijing University of Chemical Technology
Original Assignee
Beijing University of Chemical Technology
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 Beijing University of Chemical Technology filed Critical Beijing University of Chemical Technology
Priority to CN202211012278.1A priority Critical patent/CN115358469B/zh
Publication of CN115358469A publication Critical patent/CN115358469A/zh
Application granted granted Critical
Publication of CN115358469B publication Critical patent/CN115358469B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G06Q10/043Optimisation of two dimensional placement, e.g. cutting of clothes or wood
    • 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/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • 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/20Administration of product repair or maintenance
    • 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
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/10Internal combustion engine [ICE] based vehicles
    • Y02T10/40Engine management systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Human Resources & Organizations (AREA)
  • Theoretical Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Mathematical Physics (AREA)
  • Operations Research (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Quality & Reliability (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Development Economics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Algebra (AREA)
  • Game Theory and Decision Science (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Educational Administration (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Traffic Control Systems (AREA)
  • Train Traffic Observation, Control, And Security (AREA)

Abstract

本发明公开了一种获得城市交通集成传感设施的铺设最优方案的方法,构建两个共用集成传感设施的含稀疏的连续型冗余系统模型;以两个冗余系统模型为前提,确定城市交通集成传感设施铺设方案的预设约束条件及优化目标;根据预设约束条件和优化目标建立整数线性规划模型;基于目标区域的长度以及集成传感设施的感应覆盖范围确定至少一个可行铺设方案中集成传感设施的总数量;利用有限马尔可夫链嵌入法对整数线性规划模型进行求解,获得至少一个可行铺设方案的实际目标值;将符合预设约束且实际目标值最大的铺设方案确定为目标区域的城市交通集成传感设施的最优铺设方案。

Description

一种获得城市交通集成传感设施的铺设最优方案的方法
技术领域
本发明涉城市交通智能感知设备的技术领域,具体涉及一种获得城市交通集成传感设施的铺设最优方案的方法。
背景技术
近年来,伴随着我国公共交通工程建设的快速发展,信息化水平低、运维成本高、电能浪费严重等问题日渐凸显,如何对道路基础设施进行综合部署,降低运维成本并减少电能消耗是智慧城市演进过程中必须面对并坚决落实的难点问题。
在上述背景下,集成传感设施应运而生,它是依托综合管理平台对集成传感设施进行远程监控调试,并支持故障报警、故障检修和故障处理情况追踪功能,减少人工巡检工作,减少运维成本,提高了信息化水平。采用集成传感设施取代传统市政路灯、公安监控等独立杆建设模式,不仅提升了市容市貌还提高了城市的数字化管理水平。然而,此类集成传感设施的采购和维修成本高昂,因此,在交通道路上铺设适当数量的集成传感设施使得整个传感系统的可靠性最大对智慧城市的建设具有重要意义。
发明内容
针对现有技术的不足,本发明旨在提供一种获得城市交通集成传感设施的铺设最优方案的方法。
为了实现上述目的,本发明采用如下技术方案:
一种获得城市交通集成传感设施的铺设最优方案的方法,所述方法包括以下步骤:
通过将待铺设集成传感设施的目标区域抽象成一条直线,并在所述直线上以间距l铺设所述集成传感设施的方式,构建两个共用集成传感设施的线型n中取含稀疏d的连续k的F冗余系统模型;其中,n为所述直线上铺设集成传感设施的总数量;d为两个相邻失效的集成传感设施之间出现的连续正常工作集成传感设施的数目;k为含稀疏d的传感设施连续失效的数目;
以所述两个共用集成传感设施的线型n中取含稀疏d的连续k的F冗余系统模型为前提,确定所述目标区域城市交通集成传感设施的铺设方案的预设约束条件以及优化目标;
根据所述预设约束条件以及所述优化目标建立整数线性规划模型;所述整数线性规划模型为两个共用集成传感设施的线型n中取含稀疏d的连续k的F冗余系统模型的实际目标值最大化的线性表示;
基于所述目标区域的长度以及城市交通集成传感设施的感应覆盖范围确定至少一个可行铺设方案中集成传感设施的总数量;
利用有限马尔可夫链嵌入法对整数线性规划模型进行求解,获得至少一个可行铺设方案的实际目标值;将符合预设约束且实际目标值最大的铺设方案确定为目标区域的城市交通集成传感设施的最优铺设方案。
需要说明的是,整数规划模型中的决策变量为n;优化目标为最大化所述目标区域中集成传感系统正常工作的可靠度。
需要说明的是,整数线性规划模型如下表示:
Figure BDA0003811052900000031
其中:R为含稀疏d的F冗余系统模型的实际约束值;n为决策变量,表示铺设集成传感设施的总数量;A为有限马尔可夫链的转移概率矩阵;π为初始状态概率向量,表示目标区域的传感系统中所有集成传感设施都正常工作;u为向量,作用是将目标区域的传感系统中所有工作状态的概率进行求和;c为铺设一个集成传感设施的实际约束值;C为在目标区域的传感系统中铺设集成传感设施所要满足的预设约束值;L为目标区域的长度;l为铺设集成传感设施之间的间距;
Figure BDA0003811052900000032
表示整数取大。
需要说明的是,基于目标区域的长度以及所述集成传感设施的感应范围确定至少一个可行铺设方案中所述集成传感设施的总数量,包括:
基于所述目标区域的长度以及所述集成传感设备的感应范围确定多个备选铺设方案;
基于铺设一个所述集成传感设备的实际约束值确定每个备选铺设方案对应的实际约束值;
将满足预设约束条件的实际目标值对应地方案确定为可行铺设方案。
需要说明的是,利用有限马尔可夫链嵌入法对所述整数线性规划模型进行求解,获得所述至少一个可行铺设方案的实际目标值,包括:
定义状态空间S=SW∪SF上的有限马尔可夫链{Y(t),t=1,…,n};其中,工作状态空间SW和失效状态空间SF分别为:
Figure BDA0003811052900000041
SF={sN},
其中,Y(t)=(i1,j1;i2,j2)∈SW是前t个集成传感设施在t1≤t2情形下的一个工作状态,表示前t个集成传感设施组成的子系统在tr(r=1,2)时刻存在jr个集成传感设施t-jr+1,…,t正常工作,集成传感设施t-jr失效,且集成传感设施1,…,t-jr中,至多有含稀疏dr的连续ir个集成传感设施失效;失效状态空间SW={sN}中状态sN表示前t个集成传感设施组成的子系统只要在tr时刻集成传感设施1,...,t中,至少有含稀疏dr的连续kr个集成传感设施失效,即集成传感系统需要派维修人员进行维修;
对所述工作状态空间和所述失效状态空间中的任一状态进行编号;
将所有状态按编号由小到大的顺序重新排列后,利用状态转移的规确定所述有限马尔可夫链的转移概率矩阵;
基于有限马尔可夫链嵌入法对所述整数线性规划模型进行求解,获得所述可行铺设方案的实际目标值;
将符合预设约束且实际目标值最大的铺设方案确定为目标区域的城市交通集成传感设施的最优铺设方案。
本发明有益效果在于:
1、本发明构建了两个共用部件的n中取含稀疏d的连续k的F冗余系统,该系统的失效模式是出现含稀疏d的连续k个集成传感设施失效,即出现连续k个集成传感设施失效或者k个失效的集成传感设施之间最多有d个正常工作的集成传感设施时,整个集成传感系统才会失效,这样就使得在出现某一集成传感设施失效时不需要立即派人维修,即减少维修次数,节省人力、物力成本;
2、本发明提出了一种基于有限马尔可夫链嵌入法计算两个共用部件的n中取含稀疏d的连续k的F冗余系统可靠度的方法,相较于已有的递归算法在计算效率方面有明显优势。
附图说明
图1示出根据本公开一实施方式的城市交通集成传感设施的铺设方案确定方法的流程图;
图2示出根据本公开一实施方式的目标区域的集成传感系统的工作状态示意图。
具体实施方式
以下将对本发明作进一步的描述,需要说明的是,以下实施例以本技术方案为前提,给出了详细的实施方式和具体的操作过程,但本发明的保护范围并不限于本实施例。
本发明为一种获得城市交通集成传感设施的铺设最优方案的方法,所述方法包括以下步骤:
通过将待铺设集成传感设施的目标区域抽象成一条直线,并在所述直线上以间距l铺设所述集成传感设施的方式,构建两个共用集成传感设施的线型n中取含稀疏d的连续k的F冗余系统模型;其中,n为所述直线上铺设集成传感设施的总数量;d为两个相邻失效的集成传感设施之间出现的连续正常工作集成传感设施的数目;k为含稀疏d的传感设施连续失效的数目;
以所述两个共用集成传感设施的线型n中取含稀疏d的连续k的F冗余系统模型为前提,确定所述目标区域城市交通集成传感设施的铺设方案的预设约束条件以及优化目标;
根据所述预设约束条件以及所述优化目标建立整数线性规划模型;所述整数线性规划模型为两个共用集成传感设施的线型n中取含稀疏d的连续k的F冗余系统模型的实际目标值最大化的线性表示;
基于所述目标区域的长度以及城市交通集成传感设施的感应覆盖范围确定至少一个可行铺设方案中集成传感设施的总数量;
利用有限马尔可夫链嵌入法对整数线性规划模型进行求解,获得至少一个可行铺设方案的实际目标值;将符合预设约束且实际目标值最大的铺设方案确定为目标区域的城市交通集成传感设施的最优铺设方案。
需要说明的是,整数规划模型中的决策变量为n;优化目标为最大化所述目标区域中集成传感系统正常工作的可靠度。
需要说明的是,整数线性规划模型如下表示:
Figure BDA0003811052900000071
其中:R为含稀疏d的F冗余系统模型的实际约束值;n为决策变量,表示铺设集成传感设施的总数量;A为有限马尔可夫链的转移概率矩阵;π为初始状态概率向量,表示目标区域的传感系统中所有集成传感设施都正常工作;u为向量,作用是将目标区域的传感系统中所有工作状态的概率进行求和;c为铺设一个集成传感设施的实际约束值;C为在目标区域的传感系统中铺设集成传感设施所要满足的预设约束值;L为目标区域的长度;l为铺设集成传感设施之间的间距;
Figure BDA0003811052900000072
表示整数取大。
需要说明的是,基于目标区域的长度以及所述集成传感设施的感应范围确定至少一个可行铺设方案中所述集成传感设施的总数量,包括:
基于所述目标区域的长度以及所述集成传感设备的感应范围确定多个备选铺设方案;
基于铺设一个所述集成传感设备的实际约束值确定每个备选铺设方案对应的实际约束值;
将满足预设约束条件的实际目标值对应地方案确定为可行铺设方案。
需要说明的是,利用有限马尔可夫链嵌入法对所述整数线性规划模型进行求解,获得所述至少一个可行铺设方案的实际目标值,包括:
定义状态空间S=SW∪SF上的有限马尔可夫链{Y(t),t=1,…,n};其中,工作状态空间SW和失效状态空间SF分别为:
Figure BDA0003811052900000081
SF={sN},
其中,Y(t)=(i1,j1;i2,j2)∈SW是前t个集成传感设施在t1≤t2情形下的一个工作状态,表示前t个集成传感设施组成的子系统在tr(r=1,2)时刻存在jr个集成传感设施t-jr+1,…,t正常工作,集成传感设施t-jr失效,且集成传感设施1,…,t-jr中,至多有含稀疏dr的连续ir个集成传感设施失效;失效状态空间SW={sN}中状态sN表示前t个集成传感设施组成的子系统只要在tr时刻集成传感设施1,...,t中,至少有含稀疏dr的连续kr个集成传感设施失效,即集成传感系统需要派维修人员进行维修;
对所述工作状态空间和所述失效状态空间中的任一状态进行编号;
将所有状态按编号由小到大的顺序重新排列后,利用状态转移的规确定所述有限马尔可夫链的转移概率矩阵;
基于有限马尔可夫链嵌入法对所述整数线性规划模型进行求解,获得所述可行铺设方案的实际目标值;
将符合预设约束且实际目标值最大的铺设方案确定为目标区域。
需要进一步指出的是,基于上述,还包括:
状态编号
对任一工作状态重新编号(i1,j1;i2,j2)∈SW重新编号,其编号函数规则为:
Figure BDA0003811052900000091
失效状态sN∈SF的编号定义为N:=e(k1-1,d1;k2-1,d2)+1。
转移概率矩阵
假设所有集成传感设施是同一型号,即假设其满足独立同分布,且任一集成传感设施的寿命满足指数分布F(t)=1-e-λt,因此,令
Figure BDA0003811052900000092
表示集成传感设施在t1之前失效的概率;
Figure BDA0003811052900000093
表示集成传感设施在t1和t2之间失效的概率;
Figure BDA0003811052900000094
表示集成传感设施在t2之后失效的概率。将所有状态按编号从小到大的顺序排列后,该马尔可夫链的转移概率矩阵为
Figure BDA0003811052900000095
这里,A0表示条件
Figure BDA0003811052900000096
对应“集成传感设施t失效于t1时刻之前”的事件。若前t-1个集成传感设施构成子系统的工作状态为(i1,j1;i2,j2),且满足{i1=k1-1或i2=k2-1},在加入t1时刻之前已失效的集成传感设施t后,则由前t个集成传感设施构成子系统状态为sN;否则,其状态为(i1+1,0;i2+1,0).
A1表示条件
Figure BDA0003811052900000097
对应“集成传感设施t失效于t1和t2时刻之间”的事件。若前t-1个集成传感设施构成子系统的工作状态为(i1,j1;i2,j2),且满足{i2=k2-1},在加入t1和t2时刻之间失效的集成传感设施t后,则由前t个集成传感设施构成子系统状态为sN;否则,其状态为
Figure BDA0003811052900000101
A2表示条件
Figure BDA0003811052900000102
对应“集成传感设施t失效于t2时刻之后”的事件。若前t-1个集成传感设施构成子系统的工作状态为(i1,j1;i2,j2),再加入t2时刻之后失效的集成传感设施t后,则由前t个集成传感设施构成子系统不会失效,因此其状态为
Figure BDA0003811052900000103
系统可靠度
由有限马尔可夫链嵌入法可求得两个共用部件的n中取含稀疏d的连续k的F冗余系统的联合可靠度为:
R=πAnu,
其中,初始状态概率向量
Figure BDA0003811052900000104
表示所有集成传感是都正常工作;向量
Figure BDA0003811052900000105
的作用是所有工作状态的概率进行求和,即该传感系统正常工作的可靠度。
进一步地,求解整数线性规划模型;在满足预设约束的前提下,决策变量的取值有限,因此可以利用枚举法对整数线性规划模型进行求解,即给出M个铺设方案,分别计算出每一种方案的总约束值,将该总约束值与所述预设约束条件的约束阈值进行比较后,筛选出满足预设约束条件的M′种铺设方案,再分别利用联合可靠度的求解方法,计算出M′种可行方案的可靠度,选出其中可靠度最大的铺设方案,即得到整数线性规划模型的最优解。
实施例
下面将参考附图并结合具体实施例对本公开内容作详细说明,以使本领域的技术人员对其较易理解并应用。
1、将若干照明范围为24.5米,传感范围为32.5米的智慧路灯铺设在某一长为300米的路段上,为其提供照明和传感监控服务。当出现含稀疏d1的连续k1个智慧路灯失效则整个智慧路灯系统的照明功能失效;当出现含稀疏d2的连续k2个智慧路灯失效则整个智慧路灯系统的传感功能失效。因此,首先构建了两个共用部件的n中取含稀疏d的连续k的F冗余系统模型;
2、以两个共用部件的n中取含稀疏d的连续k的F冗余系统模型为前提,建立一个整数线性规划模型,这里决策变量为n,预设约束条件是铺设成本不超过投资预算C,假设每个集成传感设施的成本c=1万元,投资预算C=20万元,优化目标是传感系统的可靠度R达到最大;
3、根据优化目标和约束条件,建立如下整数线性规划模型:
Figure BDA0003811052900000111
其中:R为智慧路灯系统的可靠度;
n为决策变量,表示铺设智慧路灯的总数量;
A为有限马尔可夫链的转移概率矩阵;
π为初始状态概率向量,表示所有智慧路灯正常工作;
u为向量,作用是将智慧路灯系统中所有工作状态的概率进行求和;
c为铺设一个智慧路灯所需要的成本;
C为铺设智慧路灯的投资预算成本;
L为某一路段的长度;
l为铺设智慧路灯的间距;
Figure BDA0003811052900000121
表示整数取大,如,
Figure BDA0003811052900000122
4、求解上述整数线性规划模型,由路段的长度和智慧路灯的服务范围和不同的铺设距离(一般为15~20米)可以提供不同的备选铺设方案。这里可以给出4种备选铺设方案:
方案一:间距l=15米,铺设n=21个智慧路灯,则智慧路灯系统的照明功能失效当且仅当至少含稀疏为1的连续2个智慧路灯失效,传感功能失效当且仅当至少含稀疏为2的连续3个智慧路灯失效,即k1=2,d1=1,k2=3,d2=2;
方案二:间距l=16米,铺设n=19个智慧路灯,则智慧路灯系统的照明功能失效当且仅当至少含稀疏为1的连续2个智慧路灯失效,传感功能失效当且仅当至少含稀疏为2的连续3个智慧路灯失效,即k1=2,d1=1,k2=3,d2=2;
方案三:间距l=17米,铺设n=18个智慧路灯,则智慧路灯系统的照明功能失效当且仅当至少含稀疏为1的连续1个智慧路灯失效,传感功能失效当且仅当至少含稀疏为2的连续2个智慧路灯失效,即k1=1,d1=1,k2=2,d2=2;
方案四:间距l=18米,铺设n=17个智慧路灯,则智慧路灯系统的照明功能失效当且仅当至少含稀疏为1的连续1个智慧路灯失效,传感功能失效当且仅当至少含稀疏为2的连续2个智慧路灯失效,即k1=1,d1=1,k2=2,d2=2。分别计算出每一种方案的费用成本,与投资预算C=20万元进行比较后,筛选出满足预算的3种可行铺设方案,方案二、方案三和方案四。
5、分别利用有限马尔可夫链嵌入法求得这三种可行铺设方案的可靠度,选出其中可靠度最大的铺设方案,即整数线性规划模型的最优解,具体求解过程如下:
5.1方案二的可靠度
当铺设间距l=16米,铺设数量n=19个智慧路灯时,若连续出现含稀疏为1的2个智慧路灯失效,则智慧路灯系统的照明功能失效;
若连续出现含稀疏为2的3个智慧路灯失效,则智慧路灯系统的传感功能失效。其各个量的确定方法只需要符合间距与智慧路灯各功能的覆盖范围即可。具体求解过程如下:
5.1.1定义有限马尔可夫链
在状态空间S=SW∪SF上定义一个有限马尔可夫链{Y(t),t=1,…,19},其中工作状态空间
SW={(0,0;0,0),(0,0;1,0),(0,0;1,1),(0,0;1,2),(0,0;2,0),(0,0;2,1),(0,0;2,2),(1,0;0,0),(1,0;1,0),(1,0;1,1),(1,0;1,2),(1,0;2,0),(1,0;2,1),(1,0;2,2),(1,1;0,0),(1,1;1,0),(1,1;1,1),(1,1;1,2),(1,1;2,0),(1,1;2,1),(1,1;2,2)},
失效状态空间SF={sN},即sN为所有不在工作状态空间SW中的状态。
5.1.2状态编号
工作状态Y(t)=(i1,j1;i2,j2)∈SW的编号函数为
Figure BDA0003811052900000141
失效状态sN∈SF的编号定义为N:=e(1,1;2,2)+1。
5.1.3转移概率矩阵
所有智慧路灯是同一型号,即假设其满足独立同分布,因平均寿命为6000小时,则设任一智慧路灯的寿命满足指数分布
Figure BDA0003811052900000142
因此,令
Figure BDA0003811052900000143
表示智慧路灯在时刻t1之前失效的概率;
Figure BDA0003811052900000144
表示智慧路灯在时刻t1和t2之间失效的概率;
Figure BDA0003811052900000145
表示智慧路灯在时刻t2之后失效的概率。将所有状态按编号从小到大的顺序排列后,该马尔可夫链的转移概率矩阵为
Figure BDA0003811052900000146
这里,A0表示条件
Figure BDA0003811052900000147
对应“智慧路灯t失效于时刻t1之前”的事件。若前t-1个智慧路灯构成的子系统的工作状态为(i1,j1;i2,j2),且满足{i1=1或i2=2},在加入时刻t1之前已失效的智慧路灯t后,则由前t个智慧路灯构成的子系统的状态为sN;否则,其状态为(i1+1,0;i2+1,0).
A1表示条件
Figure BDA0003811052900000148
对应“智慧路灯t失效于时刻t1和t2之间”的事件。若前t-1个智慧路灯构成的子系统的工作状态为(i1,j1;i2,j2),且满足{i2=2},在加入时刻t1和t2之间失效的智慧路灯t后,则由前t个智慧路灯构成的子系统的状态为sN;否则,其状态为
Figure BDA0003811052900000151
A2表示条件
Figure BDA0003811052900000152
对应于“智慧路灯t失效于时刻t2之后”的事件。若前t-1个智慧路灯构成的子系统的工作状态为(i1,j1;i2,j2),在加入时刻t2之后失效的智慧路灯t后,则由前t个智慧路灯构成的子系统无法进入失效状态,只有工作状态
Figure BDA0003811052900000153
综上,可得出转移概率矩阵为
Figure BDA0003811052900000154
5.1.4系统可靠度
由有限马尔可夫链嵌入法可求得两个共用部件的n中取含稀疏d的连续k的F冗余系统的联合可靠度为:
R=πA19u,
综上可得,当t1=240(小时),t2=264(小时)时,智慧路灯系统正常工作的可靠度R=0.9500。
5.2方案三的可靠度
当铺设间距l=17米,铺设数量n=18个智慧路灯时,若连续出现含稀疏为1的1个智慧路灯失效,则智慧路灯系统的照明功能失效;若连续出现含稀疏为2的2个智慧路灯失效,则智慧路灯系统的传感功能失效。其各个量的确定方法只需要符合间距与智慧路灯各功能的覆盖范围即可。具体求解过程如下:
5.2.1定义有限马尔可夫链
在状态空间S=SW∪SF上定义一个有限马尔可夫链{Y(t),t=1,…,19},其中工作状态空间
SW={(0,0;0,0),(0,0;1,0),(0,0;1,1),(0,0;1,2)},
失效状态空间SF={sN},即sN为所有不在工作状态空间SW中的状态。
5.2.2状态编号
工作状态Y(t)=(i1,j1;i2,j2)∈SW的编号函数为
Figure BDA0003811052900000161
失效状态sN∈SF的编号定义为N:=e(0,0;1,2)+1。
5.2.3转移概率矩阵
所有智慧路灯是同一型号,即假设其满足独立同分布,因平均寿命为6000小时,则设任一智慧路灯的寿命满足指数分布
Figure BDA0003811052900000162
因此,令
Figure BDA0003811052900000163
表示智慧路灯在时刻t1之前失效的概率;
Figure BDA0003811052900000164
表示智慧路灯在时刻t1和t2之间失效的概率;
Figure BDA0003811052900000171
表示智慧路灯在时刻t2之后失效的概率。将所有状态按编号从小到大的顺序排列后,该马尔可夫链的转移概率矩阵为
Figure BDA0003811052900000172
这里,A0表示条件
Figure BDA0003811052900000173
对应“智慧路灯t失效于时刻t1之前”的事件。若前t-1个智慧路灯构成的子系统的工作状态为(i1,j1;i2,j2),且满足{i1=0或i2=1},在加入时刻t1之前已失效的智慧路灯t后,则由前t个智慧路灯构成的子系统的状态为sN;否则,其状态为(i1+1,0;i2+1,0).
A1表示条件
Figure BDA0003811052900000174
对应“智慧路灯t失效于时刻t1和t2之间”的事件。若前t-1个智慧路灯构成的子系统的工作状态为(i1,j1;i2,j2),且满足{i2=1},在加入时刻t1和t2之间失效的智慧路灯t后,则由前t个智慧路灯构成的子系统的状态为sN;否则,其状态为
Figure BDA0003811052900000175
A2表示条件
Figure BDA0003811052900000176
对应于“智慧路灯t失效于时刻t2之后”的事件。若前t-1个智慧路灯构成的子系统的工作状态为(i1,j1;i2,j2),在加入时刻t2之后失效的智慧路灯t后,则由前t个智慧路灯构成的子系统无法进入失效状态,只有工作状态
Figure BDA0003811052900000177
综上,可得出转移概率矩阵为:
Figure BDA0003811052900000181
5.2.4系统可靠度
由有限马尔可夫链嵌入法可求得两个共用部件的n中取含稀疏d的连续k的F冗余系统的联合可靠度为:
R=πA18u,
综上可得,当t1=240(小时),t2=264(小时)时,智慧路灯系统正常工作的可靠度R=0.4864。
5.3方案四的可靠度
当铺设间距l=18米,铺设数量n=17个智慧路灯时,若连续出现含稀疏为1的1个智慧路灯失效,则智慧路灯系统的照明功能失效;若连续出现含稀疏为2的2个智慧路灯失效,则智慧路灯系统的传感功能失效。其各个量的确定方法只需要符合间距与智慧路灯各功能的覆盖范围即可。具体求解过程如下:
5.3.1定义有限马尔可夫链
在状态空间S=SW∪SF上定义一个有限马尔可夫链{Y(t),t=1,…,19},其中工作状态空间
SW={(0,0;0,0),(0,0;1,0),(0,0;1,1),(0,0;1,2)},
失效状态空间SF={sN},即sN为所有不在工作状态空间SW中的状态。
5.3.2状态编号
工作状态Y(t)=(i1,j1;i2,j2)∈SW的编号函数为
Figure BDA0003811052900000191
失效状态sN∈SF的编号定义为N:=e(0,0;1,2)+1。
5.3.3转移概率矩阵
所有智慧路灯是同一型号,即假设其满足独立同分布,因平均寿命为6000小时,则设任一智慧路灯的寿命满足指数分布
Figure BDA0003811052900000192
因此,令
Figure BDA0003811052900000193
表示智慧路灯在时刻t1之前失效的概率;
Figure BDA0003811052900000194
表示智慧路灯在时刻t1和t2之间失效的概率;
Figure BDA0003811052900000195
表示智慧路灯在时刻t2之后失效的概率。将所有状态按编号从小到大的顺序排列后,该马尔可夫链的转移概率矩阵为
Figure BDA0003811052900000196
这里,A0表示条件
Figure BDA0003811052900000197
对应“智慧路灯t失效于时刻t1之前”的事件。若前t-1个智慧路灯构成的子系统的工作状态为(i1,j1;i2,j2),且满足{i1=0或i2=1},在加入时刻t1之前已失效的智慧路灯t后,则由前t个智慧路灯构成的子系统的状态为sN;否则,其状态为(i1+1,0;i2+1,0).
A1表示条件
Figure BDA0003811052900000198
对应“智慧路灯t失效于时刻t1和t2之间”的事件。若前t-1个智慧路灯构成的子系统的工作状态为(i1,j1;i2,j2),且满足{i2=1},在加入时刻t1和t2之间失效的智慧路灯t后,则由前t个智慧路灯构成的子系统的状态为sN;否则,其状态为
Figure BDA0003811052900000201
A2表示条件
Figure BDA0003811052900000202
对应于“智慧路灯t失效于时刻t2之后”的事件。若前t-1个智慧路灯构成的子系统的工作状态为(i1,j1;i2,j2),在加入时刻t2之后失效的智慧路灯t后,则由前t个智慧路灯构成的子系统无法进入失效状态,只有工作状态
Figure BDA0003811052900000203
综上,可得出转移概率矩阵为
Figure BDA0003811052900000204
5.3.4系统可靠度
由有限马尔可夫链嵌入法可求得两个共用部件的n中取含稀疏d的连续k的F冗余系统的联合可靠度为:
R=πA17u,
综上可得,当t1=240(小时),t2=264(小时)时,智慧路灯系统正常工作的可靠度R=0.5063。
通过上述计算结果可知,当智慧路灯的使用时长为t1=240(小时),t2=264(小时)时,智慧路灯系统正常工作的可靠度在第二种铺设方案下最大,因此选择第二种铺设方案作为最优方案。
对于本领域的技术人员来说,可以根据以上的技术方案和构思,给出各种相应的改变和变形,而所有的这些改变和变形,都应该包括在本发明权利要求的保护范围之内。

Claims (5)

1.一种获得城市交通集成传感设施的铺设最优方案的方法,其特征在于,所述方法包括以下步骤:
通过将待铺设集成传感设施的目标区域抽象成一条直线,并在所述直线上以间距l铺设所述集成传感设施的方式,构建两个共用集成传感设施的线型n中取含稀疏d的连续k的F冗余系统模型;其中,n为所述直线上铺设集成传感设施的总数量;d为两个相邻失效的集成传感设施之间出现的连续正常工作集成传感设施的数目;k为含稀疏d的传感设施连续失效的数目;
以所述两个共用集成传感设施的线型n中取含稀疏d的连续k的F冗余系统模型为前提,确定所述目标区域城市交通集成传感设施的铺设方案的预设约束条件以及优化目标;
根据所述预设约束条件以及所述优化目标建立整数线性规划模型;所述整数线性规划模型为两个共用集成传感设施的线型n中取含稀疏d的连续k的F冗余系统模型的实际目标值最大化的线性表示;
基于所述目标区域的长度以及城市交通集成传感设施的感应覆盖范围确定至少一个可行铺设方案中集成传感设施的总数量;
利用有限马尔可夫链嵌入法对整数线性规划模型进行求解,获得至少一个可行铺设方案的实际目标值;将符合预设约束且实际目标值最大的铺设方案确定为目标区域的城市交通集成传感设施的最优铺设方案。
2.根据权利要求1所述的获得城市交通集成传感设施的铺设最优方案的方法,其特征在于,整数规划模型中的决策变量为n;优化目标为最大化所述目标区域中集成传感系统正常工作的可靠度。
3.根据权利要求1所述的获得城市交通集成传感设施的铺设最优方案的方法,其特征在于,整数线性规划模型如下表示:
其中:R为含稀疏d的F冗余系统模型的实际约束值;n为决策变量,表示铺设集成传感设施的总数量;A为有限马尔可夫链的转移概率矩阵;π为初始状态概率向量,表示目标区域的传感系统中所有集成传感设施都正常工作;u为向量,作用是将目标区域的传感系统中所有工作状态的概率进行求和;c为铺设一个集成传感设施的实际约束值;C为在目标区域的传感系统中铺设集成传感设施所要满足的预设约束值;L为目标区域的长度;l为铺设集成传感设施之间的间距;表示整数取大。
4.根据权利要求1所述的获得城市交通集成传感设施的铺设最优方案的方法,其特征在于,基于目标区域的长度以及所述集成传感设施的感应范围确定至少一个可行铺设方案中所述集成传感设施的总数量,包括:
基于所述目标区域的长度以及所述集成传感设施的感应范围确定多个备选铺设方案;
基于铺设一个所述集成传感设施的实际约束值确定每个备选铺设方案对应的实际约束值;
将满足预设约束条件的实际目标值对应地方案确定为可行铺设方案。
5.根据权利要求1所述的获得城市交通集成传感设施的铺设最优方案的方法,其特征在于,利用有限马尔可夫链嵌入法对所述整数线性规划模型进行求解,获得所述至少一个可行铺设方案的实际目标值,包括:
定义状态空间S=SW∪SF上的有限马尔可夫链{Y(t),t=1,…,n};其中,工作状态空间SW和失效状态空间SF分别为:
SF={sN},
其中,Y(t)=(i1,j1;i2,j2)∈SW是前t个集成传感设施在t1≤t2情形下的一个工作状态,表示前t个集成传感设施组成的子系统在tr(r=1,2)时刻存在jr个集成传感设施t-jr+1,…,t正常工作,集成传感设施t-jr失效,且集成传感设施1,…,t-jr中,至多有含稀疏dr的连续ir个集成传感设施失效;失效状态空间SW={sN}中状态sN表示前t个集成传感设施组成的子系统只要在tr时刻集成传感设施1,...,t中,至少有含稀疏dr的连续kr个集成传感设施失效,即集成传感系统需要派维修人员进行维修;
对所述工作状态空间和所述失效状态空间中的任一状态进行编号;
将所有状态按编号由小到大的顺序重新排列后,利用状态转移的规确定所述有限马尔可夫链的转移概率矩阵;
基于有限马尔可夫链嵌入法对所述整数线性规划模型进行求解,获得所述可行铺设方案的实际目标值;
将符合预设约束且实际目标值最大的铺设方案确定为目标区域的城市交通集成传感设施的最优铺设方案。
CN202211012278.1A 2022-08-23 2022-08-23 一种获得城市交通集成传感设施的铺设最优方案的方法 Active CN115358469B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211012278.1A CN115358469B (zh) 2022-08-23 2022-08-23 一种获得城市交通集成传感设施的铺设最优方案的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211012278.1A CN115358469B (zh) 2022-08-23 2022-08-23 一种获得城市交通集成传感设施的铺设最优方案的方法

Publications (2)

Publication Number Publication Date
CN115358469A CN115358469A (zh) 2022-11-18
CN115358469B true CN115358469B (zh) 2023-04-07

Family

ID=84001977

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211012278.1A Active CN115358469B (zh) 2022-08-23 2022-08-23 一种获得城市交通集成传感设施的铺设最优方案的方法

Country Status (1)

Country Link
CN (1) CN115358469B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113158547A (zh) * 2021-03-11 2021-07-23 上海电力大学 计及经济性和可靠性的区域综合能源系统优化配置方法
CN113343426A (zh) * 2021-05-11 2021-09-03 山东大学 考虑现实需求导向的城市配电网线路规方法及系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2493985A (en) * 2011-08-26 2013-02-27 Intelligent Sustainable Energy Ltd Utility consumption identification

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113158547A (zh) * 2021-03-11 2021-07-23 上海电力大学 计及经济性和可靠性的区域综合能源系统优化配置方法
CN113343426A (zh) * 2021-05-11 2021-09-03 山东大学 考虑现实需求导向的城市配电网线路规方法及系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Juan Yin and Lirong Cui.Reliability for consecutive-k-out-of-n: F systems with shared components between adjacent subsystems.《Reliability Engineering &amp System Safety》.2021,第210卷全文. *
Kirtee K. Kamalja and Kalpesh P. Amrutkar.Reliability and Reliability Importance of Weighted-r -Within-Consecutive-k -out-of-n:F System.《IEEE Transactions on Reliability 》.2018,第67卷(第3期),第951-969页. *
林聪.复杂冗余系统可靠性建模与评估研究.《中国博士学位论文全文数据库 (基础科学辑)》.2019,(第09期),A003-3. *

Also Published As

Publication number Publication date
CN115358469A (zh) 2022-11-18

Similar Documents

Publication Publication Date Title
CN111737916B (zh) 一种基于大数据的道路桥梁病害分析与养护决策方法
CN110516820B (zh) 一种基于bim的钢结构桥梁信息化运维系统及处理方法
CN105512990A (zh) 一种基于灯联网传感器群的城市管理分析系统
CN111076096B (zh) 燃气管网泄漏识别方法和装置
CN107516167A (zh) 城市排水一体化运行智能调度管理系统
CN103152945A (zh) 一种消防应急照明和疏散指示系统
CN107450439A (zh) 一种路灯智能故障诊断方法
CN115147098B (zh) 一种设备检修方案自动规划方法及系统
CN115358469B (zh) 一种获得城市交通集成传感设施的铺设最优方案的方法
CN115392689A (zh) 一种古建筑安全风险动态检测方法及装置
CN115860715A (zh) 一种智慧城市路灯管控系统及方法
CN109508897B (zh) 基于双层多因素的电网与综合管廊协同规划评价方法
CN105674178A (zh) 一种地面风雨监控智能灯柱系统及监控方法
Chen et al. Analytic hierarchy process–simulation framework for lighting maintenance decision-making based on the clustered network
CN103150605A (zh) 电网规划辅助系统
CN108536980B (zh) 一种考虑可靠性因素的气体探测器离散选址优化方法
CN106781533A (zh) 一种用于道路交叉口的交通智能检测系统
CN109523422A (zh) 一种配电网故障影响因素的挖掘方法
CN103020290A (zh) 电网信息校验方法及系统
CN116886738A (zh) 一种智慧城市道路照明工程管理系统
CN116482567A (zh) 一种基于物联网的太阳能灯电量监测系统及方法
CN113139584B (zh) 用于供水管网污染物入侵点识别的传感器优化布置方法
KR102498003B1 (ko) 스마트 교차로에서 교통안전 및 사고예방을 위한 교통제어시스템
JP5352440B2 (ja) 保全管理システム、保全管理方法および保全管理プログラム
Tang et al. An Intelligent Fault Diagnosis Method for Street Lamps

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