首页 > 专利 > 杭州电子科技大学 > 一种新型冠状病毒肺炎的预测控制模型专利详情

一种新型冠状病毒肺炎的预测控制模型   0    0

有效专利 查看PDF
专利申请流程有哪些步骤?
专利申请流程图
申请
申请号:指国家知识产权局受理一件专利申请时给予该专利申请的一个标示号码。唯一性原则。
申请日:提出专利申请之日。
2021-03-23
申请公布
申请公布指发明专利申请经初步审查合格后,自申请日(或优先权日)起18个月期满时的公布或根据申请人的请求提前进行的公布。
申请公布号:专利申请过程中,在尚未取得专利授权之前,国家专利局《专利公报》公开专利时的编号。
申请公布日:申请公开的日期,即在专利公报上予以公开的日期。
2021-08-24
授权
授权指对发明专利申请经实质审查没有发现驳回理由,授予发明专利权;或对实用新型或外观设计专利申请经初步审查没有发现驳回理由,授予实用新型专利权或外观设计专利权。
2022-07-12
预估到期
发明专利权的期限为二十年,实用新型专利权期限为十年,外观设计专利权期限为十五年,均自申请日起计算。专利届满后法律终止保护。
2041-03-23
基本信息
有效性 有效专利 专利类型 发明专利
申请号 CN202110309282.3 申请日 2021-03-23
公开/公告号 CN113223721B 公开/公告日 2022-07-12
授权日 2022-07-12 预估到期日 2041-03-23
申请年 2021年 公开/公告年 2022年
缴费截止日
分类号 G16H50/50G16H50/80G16H50/30 主分类号 G16H50/50
是否联合申请 独立申请 文献类型号 B
独权数量 1 从权数量 6
权利要求数量 7 非专利引证数量 1
引用专利数量 2 被引证专利数量 0
非专利引证 1、2019.01.08陈友春.状态空间模型及其在传染病发病率预测中的应用《.南京医科大学学报》.2015,;
引用专利 WO2016022650A、US10178120B 被引证专利
专利权维持 1 专利申请国编码 CN
专利事件 事务标签 公开、实质审查、授权
申请人信息
申请人 第一申请人
专利权人 杭州电子科技大学 当前专利权人 杭州电子科技大学
发明人 张俊锋、张素焕、焦晨旭 第一发明人 张俊锋
地址 浙江省杭州市下沙高教园区2号大街 邮编 310018
申请人数量 1 发明人数量 3
申请人所在省 浙江省 申请人所在市 浙江省杭州市
代理人信息
代理机构
专利代理机构是经省专利管理局审核,国家知识产权局批准设立,可以接受委托人的委托,在委托权限范围内以委托人的名义办理专利申请或其他专利事务的服务机构。
杭州君度专利代理事务所 代理人
专利代理师是代理他人进行专利申请和办理其他专利事务,取得一定资格的人。
朱月芬
摘要
本发明公开了一种新型冠状病毒肺炎的预测控制模型。本发明模型的建立实现如下:步骤1、建立新型冠状病毒肺炎模型的状态空间模型;步骤2、设计时变矩阵处理机制;步骤3、构造新型冠状病毒肺炎模型预测控制的事件触发条件;步骤4、设计模型预测控制的框架;步骤5、设计事件触发模型预测控制器;步骤6、验证新型冠状病毒肺炎模型预测控制的正性;步骤7、验证新型冠状病毒肺炎模型预测控制的鲁棒稳定性。本发明方法利用正系统对新型冠状病毒肺炎疫情进行SEIR型建模,模型预测控制器设计、性能指标最优分析等方法,提供了一种新型冠状病毒肺炎疫情发展趋势的预测估计以及采取相应控制策略的理论方法。
  • 摘要附图
    一种新型冠状病毒肺炎的预测控制模型
  • 说明书附图:图1
    一种新型冠状病毒肺炎的预测控制模型
法律状态
序号 法律状态公告日 法律状态 法律状态信息
1 2022-07-12 授权
2 2021-08-24 实质审查的生效 IPC(主分类): G16H 50/50 专利申请号: 202110309282.3 申请日: 2021.03.23
3 2021-08-06 公开
权利要求
权利要求书是申请文件最核心的部分,是申请人向国家申请保护他的发明创造及划定保护范围的文件。
1.一种新型冠状病毒肺炎的预测控制模型,其特征在于该模型的建立实现如下:
步骤1、建立新型冠状病毒肺炎模型的状态空间模型;
步骤2、设计时变矩阵处理机制;
步骤3、构造新型冠状病毒肺炎模型预测控制的事件触发条件;
步骤4、设计模型预测控制的框架;
步骤5、设计事件触发模型预测控制器;
步骤6、验证新型冠状病毒肺炎模型预测控制的正性;
步骤7、验证新型冠状病毒肺炎模型预测控制的鲁棒稳定性;
步骤1针对新型冠状病毒肺炎疫情建立SEIR模型,根据SEIR模型的各类人群的变化采集实时数据,建立新型冠状病毒肺炎SEIR的状态空间模型,形式如下:
x(k+1)=A(k)x(k)+B(k)(α(k)sat(u(k))+(1‑α(k))β(k)u(k)), (1)其中, 表示k时刻SEIR模型中各类人群数量,n表示人群数量的种类个数,是k时刻采取的隔离措施,α(k)表示采取的隔离措施u(k)是否饱和,β(k)表示是T
否采取隔离措施,sat(u(k))=(sat(u1(k)),...,sat(um(k))) 表示饱和向量函数,和 表示k时刻的加权矩阵,由k时刻传感器采集到的数据得到;考虑到新型冠状病毒肺炎SEIR模型要求各类人群数量都为非负值的实际情况,依据正系统模型来对新型冠状病毒肺炎进行建模。

2.根据权利要求1所述的一种新型冠状病毒肺炎的预测控制模型,其特征在于步骤2设计时变矩阵处理机制,其具体实现如下:
设传感器能够采集到一组值,设计时变矩阵A(k)和B(k)位于下列的区间不确定集合:
其中,A≥0和B≥0;其中,A和B表示下界, 和 表示上界。

3.根据权利要求2所述的一种新型冠状病毒肺炎的预测控制模型,其特征在于步骤3构造新型冠状病毒肺炎模型预测控制的事件触发条件,即各类人群满足隔离条件才采取隔离措施,具体构建形式如下:
其中,常量∈给定且满足 误差 满足:
其中, 是采样状态,ks和ks+1分别是第s个和第(s+1)个事件触发时刻,其中,k∈[ks,ks+1),

4.根据权利要求3所述的一种新型冠状病毒肺炎的预测控制模型,其特征在于步骤4设计模型预测控制的框架,具体设计步骤是:
步骤4.1给定系统的状态和系统输入的隔离措施约束:
其中,δ>0和η>0;
步骤4.2设计事件触发控制律:
其中, 是控制器增益;
步骤4.3通过所设计的事件触发控制律解决下面的优化问题:
这里的性能指标函数满足:
其中,x(k+i|k)表示SEIR模型中基于第k时刻各类人群数量的第i步预测的各类人群的数量, 表示将来k时刻预测
控制措施,此外, 和
步骤4.4构造一个线性余正Lyapunov函数V(i,k):
T
V(i,k)=x(k+i|k) v, (9)
其中 引入一个鲁棒稳定性条件:
然后对(10)求取期望和求和运算,可得:
进而,得到
再计算最小的γ(k),得到
T
V(0,k)=x(k|k) v≤γ(k);(13)。

5.根据权利要求4所述的一种新型冠状病毒肺炎的预测控制模型,其特征在于步骤5设计事件触发模型预测控制器,具体实现如下:
(ι) (ι)
如果存在常数μ1>1,0<μ2<1,μ3>0,γ(k)>0和 向量v>0,ξ(k),ξ (k),ζ(k),ζ(k), 使得不等式(14)‑(23)以及步骤4.4成立,则基于步骤4.2、反馈控制器增益和吸引域增益下,步骤1中建立的状态空间模型是正的、稳定的,并满足步骤4.3中的性能指标;对于任意初始状态集Φ,状态保持在集Ψ(Hi)中;所述的不等式包括如下:
minγ(k),(14)
(ι)
μ3v+ζ (k)>0, (23)
所述的反馈控制器增益如下:
所述的吸引域增益如下:
其中, ι=1,...,m,以及p=1,...,L,


6.根据权利要求5所述的一种新型冠状病毒肺炎的预测控制模型,其特征在于步骤6模型预测控制正性分析过程如下:
步骤6.1假设x(t)∈Ψ(Hi),然后有:
其中Dl和 为对角元素为0或1的对角矩阵并且 (I为单位矩阵),
然后,在区间k∈(kp,kp+1)内有:
给定初始条件x(k0|k0)≥0,可得:
对k∈(kp,kp+1)可得
在事件触发k0时刻有
再利用F(k|k)<0和H(k)<0,可得
因此,建立的模型在(kp,kp+1)内的事件触发时刻是正的;
步骤6.2考虑时变矩阵的区间不确定方法,可得:
进一步可得:
进而可得x(k0+1|k0)≥0,通过递归推导可得x(k0+s|k0)≥0, 因此,区间系统模型是正的。

7.根据权利要求6所述的一种新型冠状病毒肺炎的预测控制模型,其特征在于步骤7新型冠状病毒肺炎模型预测控制的鲁棒稳定性分析过程如下:
根据步骤6.1有
结合步骤4.4可得
公式(35)等价于:
进而等价于:
然后考虑Dl的三种取值情况:
情况1:Dl=0时有
进而,有:
结合步骤6可得:
因此可得步骤4.4中的鲁棒稳定性条件成立;
情况2:Dl=I时,类似情况1的方法得
进而:
结合步骤6可得:
因此可得步骤4.4中的鲁棒稳定性条件成立;
情况3:Dl≠0和Dl≠I时,根据权力要求6可得
进一步:
结合步骤6可得:
因此可得步骤4.4中的鲁棒稳定性条件成立。
说明书

技术领域

[0001] 本发明属于自动化技术领域,涉及一种基于新型冠状病毒肺炎疫情的模型预测控制、优化控制、状态反馈控制等方法,综合这些方法针对SEIR模型采取有效的控制策略,可用于新型冠状病毒肺炎疫情。

背景技术

[0002] 近年来,随着互联网技术和基因测序技术的发展,病毒基因数据等多源信息都被用来分析病毒的传播过程,传染病建模的方法进入了多元化发展的时代。研究者们综合利用数学和统计模型,基于多源信息对传染病的流行规律进行更加精确地建模分析。自2019年新型冠状病毒肺炎爆发以来,严重影响了全世界人们的日常生活和经济发展。随着疫情的发展与不断地扩散,如何做到有效地疫情预测和控制成为了焦点问题。这其中SEIR模型是普遍采用的传染病建模方法。此次新型冠状病毒肺炎疫情发生之后,国内外多个研究团队利用动力学方程模型对新型冠状病毒肺炎疫情的发展态势进行了研究和分析。
[0003] 在实际生活中,存在很多量是非负的,例如,城市交通流量、网络通信、物种数目、物质浓度等等。诸如此类非负量构成的系统可以用正系统建模分析。然而,现存的传染病模型均是用一般系统方法来分析病毒扩散过程、预测感染人数。这些模型忽视了系统的本质特性,即,S、E、I、R四类人群数量均是非负的。本申请考虑到了这一特性,建立正系统模型,采用正系统的方法进行分析,更具精确性,并采取相应的控制措施,更具应用性。

发明内容

[0004] 本发明的目的是针对新型冠状病毒肺炎疫情的发展趋势,提供了一种新型冠状病毒肺炎的预测控制模型。
[0005] 本发明模型的具体步骤包括:
[0006] 步骤1、建立新型冠状病毒肺炎模型的状态空间模型;
[0007] 步骤2、设计时变矩阵处理方法;
[0008] 步骤3、构造新型冠状病毒肺炎模型预测控制的事件触发条件;
[0009] 步骤4、设计模型预测控制的框架;
[0010] 步骤5、设计事件触发模型预测控制器;
[0011] 步骤6、分析新型冠状病毒肺炎模型预测控制的正性;
[0012] 步骤7、分析新型冠状病毒肺炎模型预测控制的鲁棒稳定性。
[0013] 步骤1、针对新型冠状病毒肺炎疫情建立SEIR模型,根据SEIR模型的各类人群的变化采集实时数据,建立新型冠状病毒肺炎SEIR的状态空间模型,形式如下:
[0014] x(k+1)=A(k)x(k)+B(k)(α(k)sat(u(k))+(1‑α(k))β(k)u(k)), (1)[0015] 其中, 表示k时刻SEIR模型中各类人群数量,n表示人群数量的种类个数,是k时刻采取的隔离措施,α(k)表示采取的隔离措施u(k)是否饱和,β(k)表示是T否采取隔离措施,sat(u(k))=(sat(u1(k)),...,sat(um(k))) 表示饱和向量函数,和 表示k时刻的加权矩阵,由k时刻传感器采集到的数据得到(所有时
刻的加权矩阵即为时变矩阵);考虑到新型冠状病毒肺炎SEIR模型要求各类人群数量都为非负值的实际情况,依据正系统模型来对新型冠状病毒肺炎进行建模。
[0016] 步骤2、设计时变矩阵处理机制,其具体实现如下:
[0017] 设传感器能够采集到一组值,设计时变矩阵A(k)和B(k)位于下列的区间不确定集合:
[0018]
[0019] 其中, 和 其中,A和B表示下界, 和 表示上界;
[0020] 步骤3、构造新型冠状病毒肺炎模型预测控制的事件触发条件,即各类人群满足隔离条件才采取隔离措施,其具体构建形式如下:
[0021]
[0022] 其中,常量∈给定且满足 误差满足
[0023]
[0024] 其中, 是采样状态,ks和ks+1分别是第s个和第(s+1)个事件触发时刻,其中,k∈[ks,ks+1),
[0025] 步骤4、设计模型预测控制的框架,其具体设计步骤是:
[0026] 步骤4.1给定系统的状态和系统输入的隔离措施约束:
[0027]
[0028] 其中,δ>0和η>0。
[0029] 步骤4.2设计事件触发控制律:
[0030]
[0031] 其中, 是控制器增益。
[0032] 步骤4.3通过所设计的事件触发控制律解决下面的优化问题:
[0033]
[0034] 这里的性能指标函数满足:
[0035]
[0036] 其中,x(k+i|k)表示SEIR模型中基于第k时刻各类人群数量的第i步预测的各类人群的数量,
[0037] 表示将来k时刻预测控制措施,此外, 和
[0038] 步骤4.4构造一个线性余正Lyapunov函数V(i,k):
[0039] V(i,k)=x(k+i|k)Tv, (9)
[0040] 其中 引入一个鲁棒稳定性条件:
[0041]
[0042] 然后对(10)求取期望和求和运算,可得:
[0043]
[0044] 进而,得到
[0045]
[0046] 再计算最小的γ(k),得到
[0047] V(0,k)=x(k|k)Tv≤γ(k)。 (13)
[0048] 步骤5、设计事件触发模型预测控制器,具体如下:
[0049] 如果存在常数μ1>1,0<μ2<1,μ3>0,γ(k)>0和 向量 ξ(k),ξ(ι)(k),ζ(ι)(k),ζ (k), 使得如下不等式和步骤4.4成立,则基于步骤4.2、反馈控制器增益和和吸引域增益下,步骤1中建立的状态空间模型是正的、稳定的,并满足步骤4.3中的性能指标。对于任意初始状态集Φ,状态保持在集Ψ(Hi)中。所述的不等式包括如下:
[0050] minγ(k),(14)
[0051]
[0052]
[0053]
[0054]
[0055]
[0056]
[0057]
[0058]
[0059]
[0060] 所述的反馈控制器增益如下:
[0061]
[0062] 所述的吸引域增益如下:
[0063]
[0064] 其 中 , ι= 1 , . . . , m ,以 及 p = 1 , . . . , L ,和
[0065] 步骤6、模型预测控制正性分析过程如下:
[0066] 步骤6.1假设x(t)∈Ψ(Hi),然后有:
[0067]
[0068] 其中Dl和 为对角元素为0或1的对角矩阵并且 (I为单位矩阵),然后,在区间k∈(kp,kp+1)内有:
[0069]
[0070] 给定初始条件 可得:
[0071]
[0072] 对k∈(kp,kp+1)可得
[0073]
[0074] 在事件触发k0时刻有
[0075]
[0076] 再利用 和 可得
[0077]
[0078] 因此,建立的模型在(kp,kp+1)内的事件触发时刻是正的。
[0079] 步骤6.2考虑时变矩阵的区间不确定方法,可得
[0080] 进一步可得:
[0081]
[0082] 进而可得 通过递归推导可得 因此,区间系统模型是正的。
[0083] 步骤7、新型冠状病毒肺炎模型预测控制的鲁棒稳定性分析过程如下:
[0084] 根据步骤6.1有
[0085]
[0086]
[0087] 结合步骤4.4可得
[0088]
[0089] 这等价于
[0090]
[0091] 进而,
[0092]
[0093]
[0094] 接着,考虑Dl的三种取值情况:
[0095] 情况1:Dl=0时有
[0096]
[0097] 进而,有
[0098]
[0099] 结合步骤6可得,
[0100]
[0101] 因此可得步骤4.4中的鲁棒稳定性条件成立。
[0102] 情况2:Dl=I时,类似情况1的方法得
[0103]
[0104] 进而,
[0105]
[0106] 结合步骤6可得,
[0107]
[0108] 因此可得步骤4.4中的鲁棒稳定性条件成立。
[0109] 情况3:Dl≠0和Dl≠I时,根据权力要求6可得
[0110]
[0111]
[0112]
[0113]
[0114]
[0115]
[0116] 进一步,
[0117]
[0118] 结合步骤6可得,
[0119]
[0120] 因此可得步骤4.4中的鲁棒稳定性条件成立。
[0121] 步骤4.3中的性能指标可以通过步骤5解决。
[0122] 本发明的有益效果如下:
[0123] 本发明通过模型构建、模型预测控制器设计、性能指标最优分析等方法,提供了一种新型冠状病毒肺炎疫情发展趋势的预测估计方法,该方法针对SEIR模型中各类人群数量的动态变化预测疫情发展趋势并采取相应的控制策略,进行了有效的控制与设计。该方法能够有效的预测新型冠状病毒肺炎疫情的发展趋势并进行相应的控制。

实施方案

[0125] 如图1所示,本实例提供了一种基于新型冠状病毒肺炎的模型预测控制方法,其具体步骤如下:
[0126] 步骤1、针对新型冠状病毒肺炎疫情建立SEIR模型,根据SEIR模型的各类人群的变化采集实时数据,建立新型冠状病毒肺炎SEIR的状态空间模型,形式如下:
[0127] x(k+1)=A(k)x(k)+B(k)(α(k)sat(u(k))+(1‑α(k))β(k)u(k)), (52)[0128] 其中, 表示k时刻SEIR模型中各类人群数量,n表示人群数量的种类个数,是k时刻采取的隔离措施,α(k)表示采取的隔离措施u(k)是否饱和,β(k)表示是T否采取隔离措施,sat(u(k))=(sat(u1(k)),...,sat(um(k))) 表示饱和向量函数,和 表示k时刻的加权矩阵,由k时刻传感器采集到的数据得到(所有
时刻的加权矩阵即为时变矩阵);考虑到新型冠状病毒肺炎SEIR模型要求各类人群数量都为非负值的实际情况,依据正系统模型来对新型冠状病毒肺炎进行建模。
[0129] 步骤2、设计时变矩阵处理机制,其具体实现如下:
[0130] 设传感器能够采集到一组值,设计时变矩阵A(k)和B(k)位于下列的区间不确定集合:
[0131]
[0132] 其中, 和 其中,A和B表示下界, 和 表示上界;
[0133] 步骤3、构造新型冠状病毒肺炎模型预测控制的事件触发条件,即各类人群满足隔离条件才采取隔离措施,其具体构建形式如下:
[0134]
[0135] 其中,常量∈给定且满足 误差是采样状态,ks和ks+1分别是第s个和第(s+1)个事件触发时刻,其中,k∈[ks,ks+1),[0136] 步骤4、设计模型预测控制的框架,其具体设计步骤是:
[0137] 步骤4.1给定系统的状态和系统输入的隔离措施约束:
[0138]
[0139] 其中,δ>0和η>0。
[0140] 步骤4.2设计事件触发控制律:
[0141]
[0142] 其中, 是控制器增益。
[0143] 步骤4.3通过所设计的事件触发控制律解决下面的优化问题:
[0144]
[0145] 这里的性能指标函数满足:
[0146]
[0147] 其中,x(k+i|k)表示SEIR模型中基于第k时刻各类人群数量的第i步预测的各类人群的数量,
[0148] 表示将来k时刻预测控制措施,此外, 和
[0149] 步骤4.4构造一个线性余正Lyapunov函数V(i,k):
[0150] V(i,k)=x(k+i|k)Tv, (60)
[0151] 其中 引入一个鲁棒稳定性条件:
[0152]
[0153] 然后对(10)求取期望和求和运算,可得:
[0154]
[0155] 进而,得到
[0156]
[0157] 再计算最小的γ(k),得到
[0158] V(0,k)=x(k|k)Tv≤γ(k)。 (64)
[0159] 步骤5、设计事件触发模型预测控制器,具体如下:
[0160] 如果存在常数μ1>1,0<μ2<1,μ3>0,γ(k)>0和 向量 ξ(k),ξ(ι)(k),ζ(ι)(k),ζ (k), 使得如下不等式和步骤4.4成立,则基于步骤4.2、反馈控制器增益和和吸引域增益下,步骤1中建立的状态空间模型是正的、稳定的,并满足步骤4.3中的性能指标。对于任意初始状态集Φ,状态保持在集Ψ(Hi)中。所述的不等式包括如下:
[0161] minγ(k), (65)
[0162]
[0163]
[0164]
[0165]
[0166]
[0167]
[0168]
[0169]
[0170]
[0171] 所述的反馈控制器增益如下:
[0172]
[0173] 所述的吸引域增益如下:
[0174]
[0175] 其 中 , ι= 1 , . . . , m ,以 及 p = 1 , . . . , L ,和
[0176] 步骤6、模型预测控制正性分析过程如下:
[0177] 步骤6.1假设x(t)∈Ψ(Hi),然后有:
[0178]
[0179] 其中Dl和 为对角元素为0或1的对角矩阵并且 (I为适当维度的单位矩阵), 然后,在区间k∈(kp,kp+1)内有:
[0180]
[0181] 给定初始条件 可得:
[0182]
[0183] 对k∈(kp,kp+1)可得
[0184]
[0185] 在事件触发k0时刻有
[0186]
[0187] 再利用 和 可得
[0188]
[0189] 因此,建立的模型在(kp,kp+1)内的事件触发时刻是正的。
[0190] 步骤6.2考虑时变矩阵的区间不确定方法,可得
[0191]
[0192] 进一步可得:
[0193]
[0194] 进而可得 通过递归推导可得 因此,区间系统模型是正的。
[0195] 步骤7、新型冠状病毒肺炎模型预测控制的鲁棒稳定性分析过程如下:
[0196] 根据步骤6.1有
[0197]
[0198] 结合步骤4.4可得
[0199]
[0200]
[0201] 这等价于
[0202]
[0203] 进而,
[0204]
[0205] 接着,考虑Dl的三种取值情况:
[0206] 情况1:Dl=0时有
[0207]
[0208] 进而,有
[0209]
[0210] 结合步骤6可得:
[0211]
[0212] 因此可得步骤4.4中的鲁棒稳定性条件成立。
[0213] 情况2:Dl=I时,类似情况1的方法得
[0214]
[0215] 进而,
[0216]
[0217]
[0218] 结合步骤6可得:
[0219]
[0220] 因此可得步骤4.4中的鲁棒稳定性条件成立。情况3:Dl≠0和Dl≠I时,根据权力要求6可得
[0221]
[0222]
[0223]
[0224]
[0225]
[0226] 进一步,
[0227]
[0228]
[0229] 结合步骤6可得:
[0230]
[0231] 因此可得步骤4.4中的鲁棒稳定性条件成立。
[0232] 步骤4.3中的性能指标可以通过优化权力要求6解决。

附图说明

[0124] 图1是本发明框架图。
版权所有:盲专网 ©2023 zlpt.xyz  蜀ICP备2023003576号