CN113345596A - 一种基于seir模型的传染病系统的事件触发滤波方法 - Google Patents

一种基于seir模型的传染病系统的事件触发滤波方法 Download PDF

Info

Publication number
CN113345596A
CN113345596A CN202110703319.0A CN202110703319A CN113345596A CN 113345596 A CN113345596 A CN 113345596A CN 202110703319 A CN202110703319 A CN 202110703319A CN 113345596 A CN113345596 A CN 113345596A
Authority
CN
China
Prior art keywords
infectious disease
seir
filter
representing
event
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.)
Granted
Application number
CN202110703319.0A
Other languages
English (en)
Other versions
CN113345596B (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 CN202110703319.0A priority Critical patent/CN113345596B/zh
Publication of CN113345596A publication Critical patent/CN113345596A/zh
Application granted granted Critical
Publication of CN113345596B publication Critical patent/CN113345596B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/80ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for detecting, monitoring or modelling epidemics or pandemics, e.g. flu
    • 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/15Correlation function computation including computation of convolution operations

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Public Health (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Databases & Information Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Computing Systems (AREA)
  • Biomedical Technology (AREA)
  • Algebra (AREA)
  • Pathology (AREA)
  • Epidemiology (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Investigating Or Analysing Biological Materials (AREA)

Abstract

本发明公开了一种基于SEIR模型的传染病系统的事件触发滤波方法。本发明包括如下步骤:步骤1、建立SEIR型传染病中各类人群数量的状态空间模型;步骤2、建立SEIR型传染病中各类人群良好发展的事件触发条件;步骤3、设计估计SEIR型传染病中各类人群数量变化的事件触发滤波器;步骤4、SEIR型传染病系统正性验证;步骤5、SEIR型传染病系统随机稳定性。本发明不仅有效预测SEIR型传染病系统中各类人群数量变化,且在系统稳定运行的情况下,节省了一定的系统资源。

Description

一种基于SEIR模型的传染病系统的事件触发滤波方法
技术领域
本发明属于传染病预测建模技术领域,涉及一种预测SEIR型传染病系统中各类人群数量变化的模型,尤其涉及一种传染病的事件触发滤波方法。
背景技术
2019年12月以来,一场突如其来的新型冠状病毒肺炎(COVID-19)疫情蔓延开来,全球经历了由SARS-CoV-2病毒引起的新型冠状病毒肺炎大流行病的空前威胁。这给人们的正常生活和生命健康造成了重大影响。这一问题引起了全球广泛关注,研究者们综合利用数学和统计模型,基于多源信息对COVID-19流行病的发展规律进行更加精确地建模分析,为COVID-19流行病的防控提供了重要的理论依据和解决途径。但是现存的模型(SI、SIS、SIR、SIRS、SEIR)大多以确定性微分方程为基础,忽略了各类人群数量变化为非负的特点,也没有实施相应的控制措施,具有较大的保守性、少应用性。另外,在SEIR型传染病模型中,人群划分方法包括以下几类:易感人群S、潜伏人群E、已经感染人群I、康复/死亡的人群R。
事件触发机制是一种可以节省设计成本、节约通信资源的控制策略。在对COVID-19流行病采取措施的过程中,一些限制条件,比如可以用事件触发机制描述满足隔离人群的条,避免出现漏隔离或误隔离。因此,引入事件触发机制是必要的。另外,针对某一地域的人群采取控制措施的过程中,不可避免地会有该区域外的人员流入。将外来人员看作扰动输入,进而可用滤波器估计各类人群的数量变化,可以借助正半马尔科夫跳变系统的滤波器设计估计四类人群的数量。
针对此问题,虽有不少研究者建立SEIR数学模型预测COVID-19流行病的发展趋势。但采用一般系统的分析方法来研究SEIR模型预测流行病的各类人群数量,这忽略了该模型的本质特征,即S、E、I、R四类人群数量具有非负特性。而且,上述四类人群数量变化符合正半马尔科夫跳变系统的特性,且在实际生活中流行病爆发后需要采取一定的控制策略抑制流行病的传播、最终消灭流行病。因此,本发明利用现代控制理论技术建立出染病的事件触发滤波估计的状态空间模型,分析了其正性及稳定性,设计了滤波器增益矩阵,降低了系统设计成本,使传染病滤波估计系统保持稳定运行。综上所述,设计一种预测SEIR型传染病系统中各类人群数量变化的事件触发滤波估计方法对人类生活、社会发展具有重要的意义。
发明内容
针对现有技术存在的问题,本发明提供了一种利用现代控制理论技术建立了传染病系统事件触发滤波估计的状态空间模型,设计了一类事件触发滤波器,分析了其正性及稳定性,降低了系统设计成本,预测了SEIR型传染病中各类人群数量的变化趋势。
本发明方法的具体步骤包括如下:
步骤1、SEIR型传染病中各类人群数量的状态空间模型,具体步骤如下:
步骤1.1采集SEIR型传染病中各类人群数量的数据,利用采集到的数据结合隔离条件,建立SEIR型传染病的人群数量的状态空间模型,具体方法是:
Figure BDA0003131086150000021
其中,
Figure BDA0003131086150000022
表示t时刻传染病的各类人群数量的变化率;
Figure BDA0003131086150000023
表示在t时刻传染病的各类人群数量,n表示人群的类别数;
Figure BDA0003131086150000024
代表传感器测量输出数量;
Figure BDA0003131086150000025
表示期望输出数量;
Figure BDA0003131086150000026
代表影响各类人群数量变化的外来人群等外在因素的干扰,m表示外在因素数量;函数rt为半马尔科夫跳变过程也即传染病预测系统的跳变信号,在有限集
Figure BDA0003131086150000027
内取值,
Figure BDA0003131086150000028
为系统权重矩阵,反应各类人群之间的影响,
Figure BDA0003131086150000029
表示系统在t时刻的不确定摄动项,α(t)表示是服从伯努利分布的随机信号。
Figure BDA00031310861500000210
Figure BDA00031310861500000211
均为实数矩阵。记rt=i,,
Figure BDA00031310861500000212
Figure BDA00031310861500000213
步骤1.2设计半马尔科夫跳变信号rt,其转移速率满足以下条件:
Figure BDA00031310861500000214
其中,转移率λij(ζ)≥0表示跳变信号从模态i到模态j,
Figure BDA0003131086150000031
并且ζ>0。
步骤1.3设计步骤1.1中建立模型的加性不确定摄动ΔAi(t),具体形式如下:
Figure BDA0003131086150000032
其中,
Figure BDA0003131086150000033
是一个未知的时变向量函数,并且满足条件
Figure BDA0003131086150000034
Figure BDA0003131086150000035
步骤2、建立SEIR型传染病中各类人群良好发展的事件触发条件,其构建形式如下:
‖ey(t)‖1>β‖y(t)‖1 (4)
其中,常量0<β<1,ey(t)属于采样误差,
Figure BDA0003131086150000036
Figure BDA0003131086150000037
表示采样输出,y(t)表示估计输出数量,‖·‖1表示取1范数运算。
步骤3、设计估计SEIR型传染病中各类人群数量变化的事件触发滤波器,其具体构建如下:
步骤3.1滤波器系统设计:
Figure BDA0003131086150000038
其中,xf(t)是滤波器状态的变化率;
Figure BDA0003131086150000039
是滤波器状态;
Figure BDA00031310861500000310
是滤波器输出;Afi、Bfi、Efi以及Ffi是需要求解的滤波器增益矩阵。
步骤3.2给出事件触发滤波器如下:
Figure BDA00031310861500000311
其中,Gi表示传感器的测量输出。
步骤3.3基于步骤2.1和步骤3.1,设计滤波器误差系统如下:
Figure BDA0003131086150000041
其中,
Figure BDA0003131086150000042
以及
Figure BDA0003131086150000043
Figure BDA0003131086150000044
Figure BDA0003131086150000045
步骤3.4设计预测SEIR型传染病滤波估计系统的条件,具体如下:
设计常数
Figure BDA0003131086150000046
1,∈2,设计
Figure BDA0003131086150000047
向量
Figure BDA0003131086150000048
同时设计
Figure BDA0003131086150000049
向量
Figure BDA00031310861500000410
使得以下不等式对于
Figure BDA00031310861500000411
Figure BDA00031310861500000412
成立:
Figure BDA00031310861500000413
Figure BDA00031310861500000414
Figure BDA00031310861500000415
Figure BDA00031310861500000416
Figure BDA00031310861500000417
Figure BDA00031310861500000418
Figure BDA00031310861500000419
Figure BDA00031310861500000420
Figure BDA00031310861500000421
Figure BDA0003131086150000051
Figure BDA0003131086150000052
Figure BDA0003131086150000053
其中Φ=I+β1m×m,Ψ=I-β1m×m
Figure BDA0003131086150000054
由此可见,在事件触发滤波器增益矩阵下,传染病滤波估计系统是正的且随机稳定的。
步骤4、SEIR型传染病系统正性验证过程如下:
对于
Figure BDA0003131086150000055
Figure BDA0003131086150000056
不等式
Figure BDA0003131086150000057
成立。
进而,有
Figure BDA0003131086150000058
因此,能够得到步骤3.3中误差系统的下界:
Figure BDA0003131086150000059
其中,
Figure BDA00031310861500000510
Figure BDA00031310861500000511
Figure BDA00031310861500000512
根据步骤3.4的前五个条件可得:
Figure BDA0003131086150000061
进而可得:对于
Figure BDA0003131086150000062
Figure BDA0003131086150000063
通过递归推导,可得
Figure BDA0003131086150000064
Figure BDA0003131086150000065
故,所述的传染病滤波估计系统是正的。
步骤5、SEIR型传染病系统随机稳定性的验证过程如下:
步骤3.3中误差系统的上界:
Figure BDA0003131086150000066
其中,
Figure BDA0003131086150000067
Figure BDA0003131086150000068
Figure BDA0003131086150000069
构造一个随机余正Lyapunov函数如下:
Figure BDA00031310861500000610
其中,
Figure BDA00031310861500000611
并且取值满足
Figure BDA00031310861500000612
即向量中的每一个元素均为正数,且
Figure BDA00031310861500000613
对上述的Lyapunov函数取弱无穷小算子,具体如下:
Figure BDA00031310861500000614
下面分两个步骤讨论无/有外扰输入影响的传染病滤波估计系统的随机稳定性。
步骤5.1考虑ω(t)=0,结合
Figure BDA0003131086150000071
的事实,可得:
Figure BDA0003131086150000072
根据步骤3.4中的条件可得:
Figure BDA0003131086150000073
Figure BDA0003131086150000074
Figure BDA0003131086150000075
进而,可得:
Figure BDA0003131086150000076
从而可得:
Figure BDA0003131086150000077
步骤5.2当ω(t)≠0时,结合步骤3.4可得:
Figure BDA0003131086150000078
Figure BDA0003131086150000079
Figure BDA00031310861500000710
Figure BDA00031310861500000711
进一步,根据上述推导可得:
Figure BDA00031310861500000712
即当t→∞时,
Figure BDA00031310861500000713
成立。
其中,E{·}表示数学期望,dτ表示时间变量算子、γ是扰动增益抑制水平,所述的传染病滤波估计系统是随机稳定的且具有l1增益性能γ。
本发明的有益效果如下:
本发明方法针对当前所述的SEIR型传染病系统中各类人群数量变化的情况,建立了SEIR型传染病系统的状态空间模型,借助矩阵论知识和构造的随机余正Lyapunov函数,设计了一类事件触发滤波器,并分析了其正性及稳定性,预测了传染病系统中各类人群数量变化的情况。
本发明提供了一种利用现代控制理论技术建立了传染病系统事件触发滤波估计的状态空间模型,设计了一类事件触发滤波器,分析了其正性及稳定性,降低了系统设计成本,预测了SEIR型传染病中各类人群数量的变化趋势。
附图说明
图1是本发明所述的SEIR传染病模型示意图。
具体实施方式
下面结合具体实施例来对本发明进行进一步说明。
本实例提供了一种预测SEIR型传染病系统中各类人群数量变化的事件触发滤波估计方法,其具体步骤如下:
步骤1、通过采集SEIR型传染病中各类人群数量的数据,利用采集到的数据结合隔离条件,建立SEIR型传染病的人群数量的状态空间模型:
具体构造如下步骤:
步骤1.1建立半马尔科夫跳变过程的传染病系统:
Figure BDA0003131086150000081
其中,
Figure BDA0003131086150000082
表示t时刻传染病的各类人群数量的变化率;
Figure BDA0003131086150000083
表示在t时刻传染病的各类人群数量,n表示人群的类别数;
Figure BDA0003131086150000084
代表传感器测量输出数量;
Figure BDA0003131086150000085
表示期望输出数量;
Figure BDA0003131086150000086
代表影响各类人群数量变化的外来人群等外在因素的干扰,m表示外在因素数量;函数rt为半马尔科夫跳变过程也即传染病预测系统的跳变信号,在有限集
Figure BDA0003131086150000087
内取值,
Figure BDA0003131086150000091
为系统权重矩阵,反应各类人群之间的影响,
Figure BDA0003131086150000092
表示系统在t时刻的不确定摄动项,α(t)表示是服从伯努利分布的随机信号。
Figure BDA0003131086150000093
Figure BDA0003131086150000094
均为实数矩阵。记rt=i,,
Figure BDA0003131086150000095
Figure BDA0003131086150000096
步骤1.2设计半马尔科夫跳变信号rt,其转移速率满足以下条件:
Figure BDA0003131086150000097
其中,转移率λij(ζ)≥0表示跳变信号从模态i到模态j,
Figure BDA0003131086150000098
并且ζ>0。
步骤1.3设计步骤1.1中建立模型的加性不确定摄动ΔAi(t),具体形式如下:
Figure BDA0003131086150000099
其中,
Figure BDA00031310861500000910
是一个未知的时变向量函数,并且满足条件
Figure BDA00031310861500000911
Figure BDA00031310861500000912
步骤2、建立SEIR型传染病中各类人群良好发展的事件触发控制条件,其构建形式如下:
‖ey(t)‖1>β‖y(t)‖1,其中,常量0<β<1,ey(t)属于采样误差,
Figure BDA00031310861500000913
Figure BDA00031310861500000914
表示采样输出。
步骤2、建立SEIR型传染病中各类人群良好发展的事件触发条件,其构建形式如下:
‖ey(t)‖1>β‖y(t)‖1 (4)
其中,常量0<β<1,ey(t)属于采样误差,
Figure BDA00031310861500000915
t∈[tk,tk+1),
Figure BDA00031310861500000916
表示采样输出,y(t)表示估计输出数量,‖·‖1表示取1范数运算。
步骤3、设计估计SEIR型传染病中各类人群数量变化的事件触发滤波器,其具体构建如下:
步骤3.1滤波器系统设计:
Figure BDA0003131086150000101
其中,xf(t)是滤波器状态的变化率;
Figure BDA0003131086150000102
是滤波器状态;
Figure BDA0003131086150000103
是滤波器输出;Afi、Bfi、Efi以及Ffi是需要求解的滤波器增益矩阵。
步骤3.2给出事件触发滤波器如下:
Figure BDA0003131086150000104
其中,Gi表示传感器的测量输出。
步骤3.3基于步骤2.1和步骤3.1,设计滤波器误差系统如下:
Figure BDA0003131086150000105
其中,
Figure BDA0003131086150000106
以及
Figure BDA0003131086150000107
Ei=(Efi+FfiGiCi-Ei Efi),
Figure BDA0003131086150000108
步骤3.4设计预测SEIR型传染病滤波估计系统的条件,具体如下:
设计常数
Figure BDA0003131086150000109
Figure BDA00031310861500001010
1,∈2,设计
Figure BDA00031310861500001011
向量
Figure BDA00031310861500001012
同时设计
Figure BDA00031310861500001013
向量
Figure BDA00031310861500001014
使得以下不等式对于
Figure BDA00031310861500001015
Figure BDA00031310861500001016
成立:
Figure BDA00031310861500001017
Figure BDA00031310861500001018
Figure BDA0003131086150000111
Figure BDA0003131086150000112
Figure BDA0003131086150000113
Figure BDA0003131086150000114
Figure BDA0003131086150000115
Figure BDA0003131086150000116
Figure BDA0003131086150000117
Figure BDA0003131086150000118
Figure BDA0003131086150000119
Figure BDA00031310861500001110
其中Φ=I+β1m×m,Ψ=I-β1m×m
Figure BDA00031310861500001111
由此可见,在事件触发滤波器增益矩阵下,传染病滤波估计系统是正的且随机稳定的。
步骤4、SEIR型传染病系统正性验证过程如下:
对于
Figure BDA00031310861500001112
Figure BDA00031310861500001113
不等式
Figure BDA00031310861500001114
成立。
进而,有
Figure BDA00031310861500001115
因此,能够得到步骤3.3中误差系统的下界:
Figure BDA00031310861500001116
其中,
Figure BDA0003131086150000121
Figure BDA0003131086150000122
根据步骤3.4的前五个条件可得:
Figure BDA0003131086150000123
进而可得:对于
Figure BDA0003131086150000124
Figure BDA0003131086150000125
通过递归推导,可得
Figure BDA0003131086150000126
Figure BDA0003131086150000127
故,所述的传染病滤波估计系统是正的。
步骤5、SEIR型传染病系统随机稳定性的验证过程如下:
步骤3.3中误差系统的上界:
Figure BDA0003131086150000128
其中,
Figure BDA0003131086150000129
Figure BDA00031310861500001210
构造一个随机余正Lyapunov函数如下:
Figure BDA00031310861500001211
其中,
Figure BDA00031310861500001212
并且取值满足
Figure BDA00031310861500001213
即向量中的每一个元素均为正数,且
Figure BDA00031310861500001214
对上述的Lyapunov函数取弱无穷小算子,具体如下:
Figure BDA0003131086150000131
下面分两个步骤讨论无/有外扰输入影响的传染病滤波估计系统的随机稳定性。
步骤5.1考虑ω(t)=0,结合
Figure BDA0003131086150000132
的事实,可得:
Figure BDA0003131086150000133
根据步骤3.4中的条件可得:
Figure BDA0003131086150000134
Figure BDA0003131086150000135
Figure BDA0003131086150000136
进而,可得:
Figure BDA0003131086150000137
从而可得:
Figure BDA0003131086150000138
步骤5.2当ω(t)≠0时,结合步骤3.4可得:
Figure BDA0003131086150000139
Figure BDA00031310861500001310
Figure BDA00031310861500001311
Figure BDA00031310861500001312
进一步,根据上述推导可得:
Figure BDA0003131086150000141
即当t→∞时,
Figure BDA0003131086150000142
成立。
其中,E{·}表示数学期望,dτ表示时间变量算子、γ是扰动增益抑制水平,所述的传染病滤波估计系统是随机稳定的且具有l1增益性能γ。
本发明利用正半马尔科夫跳变系统对SEIR型传染病系统进行建模,设计了事件触发滤波器。首先,利用传感器设备采集系统自身固有信息,建SEIR型传染病系统的状态空间模型;然后,针对SEIR型传染病系统中各类人群数量变化情况设计一类事件触发滤波器,使得系统能够在加性和乘性模型不确定两种情况下正常运行。本发明不仅有效解决了模型不确定下系统的稳定运行,节省了一定的系统资源,预测了各类人群数量的变化趋势。

Claims (6)

1.一种基于SEIR模型的传染病系统的事件触发滤波方法,其特征在于包括如下步骤:
步骤1、建立SEIR型传染病中各类人群数量的状态空间模型;
步骤2、建立SEIR型传染病中各类人群良好发展的事件触发条件;
步骤3、设计估计SEIR型传染病中各类人群数量变化的事件触发滤波器;
步骤4、SEIR型传染病系统正性验证;
步骤5、SEIR型传染病系统随机稳定性。
2.根据权利要求1所述的一种基于SEIR模型的传染病系统的事件触发滤波方法,其特征在于步骤1具体步骤如下:
步骤1.1采集SEIR型传染病中各类人群数量的数据,利用采集到的数据结合隔离条件,建立SEIR型传染病的人群数量的状态空间模型,具体方法是:
Figure FDA0003131086140000011
其中,
Figure FDA0003131086140000012
表示t时刻传染病的各类人群数量的变化率;
Figure FDA0003131086140000013
表示在t时刻传染病的各类人群数量,n表示人群的类别数;
Figure FDA0003131086140000014
代表传感器测量输出数量;
Figure FDA0003131086140000015
表示期望输出数量;
Figure FDA0003131086140000016
代表影响各类人群数量变化的外来人群等外在因素的干扰,m表示外在因素数量;函数rt为半马尔科夫跳变过程也即传染病预测系统的跳变信号,在有限集
Figure FDA0003131086140000017
内取值,
Figure FDA0003131086140000018
为系统权重矩阵,反应各类人群之间的影响,
Figure FDA0003131086140000019
表示系统在t时刻的不确定摄动项,α(t)表示是服从伯努利分布的随机信号;
Figure FDA00031310861400000110
Figure FDA00031310861400000111
均为实数矩阵;记rt=i,,
Figure FDA00031310861400000113
Figure FDA00031310861400000112
步骤1.2设计半马尔科夫跳变信号rt,其转移速率满足以下条件:
Figure FDA0003131086140000021
其中,转移率λij(ζ)≥0表示跳变信号从模态i到模态j,
Figure FDA0003131086140000022
并且ζ>0;
步骤1.3设计步骤1.1中建立模型的加性不确定摄动ΔAi(t),具体形式如下:
Figure FDA0003131086140000023
其中,
Figure FDA0003131086140000024
是一个未知的时变向量函数,并且满足条件
Figure FDA0003131086140000025
Figure FDA0003131086140000026
3.根据权利要求2所述的一种基于SEIR模型的传染病系统的事件触发滤波方法,其特征在于步骤2具体构建形式如下:
‖ey(t)‖1>β‖y(t)‖1 (4)
其中,常量0<β<1,ey(t)属于采样误差,
Figure FDA0003131086140000027
Figure FDA0003131086140000028
表示采样输出,y(t)表示估计输出数量,‖·‖1表示取1范数运算。
4.根据权利要求3所述的一种基于SEIR模型的传染病系统的事件触发滤波方法,其特征在于步骤3所述的设计估计SEIR型传染病中各类人群数量变化的事件触发滤波器,其具体构建如下:
步骤3.1滤波器系统设计:
Figure FDA0003131086140000029
其中,xf(t)是滤波器状态的变化率;
Figure FDA00031310861400000210
是滤波器状态;
Figure FDA00031310861400000211
是滤波器输出;Afi、Bfi、Efi以及Ffi是需要求解的滤波器增益矩阵;
步骤3.2给出事件触发滤波器如下:
Figure FDA0003131086140000031
其中,Gi表示传感器的测量输出;
步骤3.3基于步骤2.1和步骤3.1,设计滤波器误差系统如下:
Figure FDA0003131086140000032
其中,
Figure FDA0003131086140000033
以及
Figure FDA0003131086140000034
Figure FDA0003131086140000035
Figure FDA0003131086140000036
步骤3.4设计预测SEIR型传染病滤波估计系统的条件,具体如下:
设计常数
Figure FDA00031310861400000321
1,∈2,设计
Figure FDA0003131086140000038
向量
Figure FDA0003131086140000039
同时设计
Figure FDA00031310861400000310
向量
Figure FDA00031310861400000311
使得以下不等式对于
Figure FDA00031310861400000312
Figure FDA00031310861400000313
成立:
Figure FDA00031310861400000314
Figure FDA00031310861400000315
Figure FDA00031310861400000316
Figure FDA00031310861400000317
Figure FDA00031310861400000318
Figure FDA00031310861400000319
Figure FDA00031310861400000320
Figure FDA0003131086140000041
Figure FDA0003131086140000042
Figure FDA0003131086140000043
Figure FDA0003131086140000044
Figure FDA0003131086140000045
其中Φ=I+β1m×m,Ψ=I-β1m×m
Figure FDA0003131086140000046
由此可见,在事件触发滤波器增益矩阵下,传染病滤波估计系统是正的且随机稳定的。
5.根据权利要求4所述的一种基于SEIR模型的传染病系统的事件触发滤波方法,其特征在于步骤4所述的SEIR型传染病系统正性验证过程如下:
对于
Figure FDA0003131086140000047
Figure FDA0003131086140000048
不等式
Figure FDA0003131086140000049
成立;
进而,有
Figure FDA00031310861400000410
因此,能够得到步骤3.3中误差系统的下界:
Figure FDA00031310861400000411
其中,
Figure FDA00031310861400000412
Figure FDA00031310861400000413
Figure FDA00031310861400000414
根据步骤3.4的前五个条件可得:
Figure FDA0003131086140000051
进而可得:对于
Figure FDA0003131086140000052
Figure FDA0003131086140000053
通过递归推导,可得
Figure FDA0003131086140000054
Figure FDA0003131086140000055
故,所述的传染病滤波估计系统是正的。
6.根据权利要求5所述的一种基于SEIR模型的传染病系统的事件触发滤波方法,其特征在于步骤5所述的SEIR型传染病系统随机稳定性的验证过程如下:
步骤3.3中误差系统的上界:
Figure FDA0003131086140000056
其中,
Figure FDA0003131086140000057
Figure FDA0003131086140000058
Figure FDA0003131086140000059
构造一个随机余正Lyapunov函数如下:
Figure FDA00031310861400000510
其中,
Figure FDA00031310861400000511
并且取值满足
Figure FDA00031310861400000512
即向量中的每一个元素均为正数,且
Figure FDA00031310861400000513
对上述的Lyapunov函数取弱无穷小算子,具体如下:
Figure FDA00031310861400000514
下面分两个步骤讨论无/有外扰输入影响的传染病滤波估计系统的随机稳定性:
步骤5.1考虑ω(t)=0,结合
Figure FDA0003131086140000061
的事实,可得:
Figure FDA0003131086140000062
根据步骤3.4中的条件可得:
Figure FDA0003131086140000063
Figure FDA0003131086140000064
Figure FDA0003131086140000065
进而,可得:
Figure FDA0003131086140000066
从而可得:
Figure FDA0003131086140000067
步骤5.2当ω(t)≠0时,结合步骤3.4可得:
Figure FDA0003131086140000068
Figure FDA0003131086140000069
Figure FDA00031310861400000610
Figure FDA00031310861400000611
进一步,根据上述推导可得:
Figure FDA00031310861400000612
即当t→∞时,
Figure FDA0003131086140000071
成立;
其中,E{·}表示数学期望,dτ表示时间变量算子、γ是扰动增益抑制水平,所述的传染病滤波估计系统是随机稳定的且具有
Figure FDA0003131086140000072
增益性能γ。
CN202110703319.0A 2021-06-24 2021-06-24 一种基于seir模型的传染病系统的事件触发滤波方法 Active CN113345596B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110703319.0A CN113345596B (zh) 2021-06-24 2021-06-24 一种基于seir模型的传染病系统的事件触发滤波方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110703319.0A CN113345596B (zh) 2021-06-24 2021-06-24 一种基于seir模型的传染病系统的事件触发滤波方法

Publications (2)

Publication Number Publication Date
CN113345596A true CN113345596A (zh) 2021-09-03
CN113345596B CN113345596B (zh) 2024-09-13

Family

ID=77478154

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110703319.0A Active CN113345596B (zh) 2021-06-24 2021-06-24 一种基于seir模型的传染病系统的事件触发滤波方法

Country Status (1)

Country Link
CN (1) CN113345596B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107065545A (zh) * 2017-04-01 2017-08-18 同济大学 基于马尔科夫跳变的分布式事件触发滤波系统及设计方法
CN112147901A (zh) * 2020-10-09 2020-12-29 杭州电子科技大学 一种状态饱和水务系统的事件触发滤波器设计方法
CN113009825A (zh) * 2021-02-08 2021-06-22 云境商务智能研究院南京有限公司 一种受欺骗攻击的非线性网络化系统状态估计方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107065545A (zh) * 2017-04-01 2017-08-18 同济大学 基于马尔科夫跳变的分布式事件触发滤波系统及设计方法
CN112147901A (zh) * 2020-10-09 2020-12-29 杭州电子科技大学 一种状态饱和水务系统的事件触发滤波器设计方法
CN113009825A (zh) * 2021-02-08 2021-06-22 云境商务智能研究院南京有限公司 一种受欺骗攻击的非线性网络化系统状态估计方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
HUIYING CHEN ET AL.: "Event-Triggered Dissipative Filter Design for Semi-Markovian Jump Systems with Time-Varying Delays", 《MATHEMATICAL PROBLEMS IN ENGINEERING》 *
YU SHAO ET AL.: "Event-Triggered Filter of Switched Positive Systems with State Saturation", 《2020 7TH INTERNATIONAL CONFERENCE ON INFORMATION, CYBERNETICS, AND COMPUTATIONAL SOCIAL SYSTEMS (ICCSS)》 *
李雪等: "一类随机SEIQR传染病模型的动力学行为分析", 《纺织高校基础科学学报》, no. 01 *

Also Published As

Publication number Publication date
CN113345596B (zh) 2024-09-13

Similar Documents

Publication Publication Date Title
Arino et al. A simple model for COVID-19
Luo et al. Time series prediction of COVID-19 transmission in America using LSTM and XGBoost algorithms
Fennell et al. Limitations of discrete-time approaches to continuous-time contagion dynamics
Kiss et al. Generalization of pairwise models to non-Markovian epidemics on networks
Gillijns et al. What is the ensemble Kalman filter and how well does it work?
Toth et al. The role of heterogeneity in contact timing and duration in network models of influenza spread in schools
CN102469103B (zh) 基于bp神经网络的木马事件预测方法
Chen et al. Event-triggered state estimation: Experimental performance assessment and comparative study
De la Sen et al. On an Sir Epidemic Model for the COVID‐19 Pandemic and the Logistic Equation
CN112259248B (zh) 一种covid-19境外输入风险评估预测方法
CN113223731A (zh) 深度时序n-seir传染病预测模型
Röst et al. Impact of non-Markovian recovery on network epidemics
CN113096821A (zh) 一种基于多仓室模型动态特征的流行病动力学预测方法
Kepler et al. Decomposing the sources of SARS-CoV-2 fitness variation in the United States
Xu et al. An extension sample classification‐based extreme learning machine ensemble method for process fault diagnosis
Zachreson et al. The effects of local homogeneity assumptions in metapopulation models of infectious disease
van der Laan Causal inference for networks
CN113345596A (zh) 一种基于seir模型的传染病系统的事件触发滤波方法
CN117153423A (zh) 基于贝叶斯推理的新发传染病暴发时间预测方法
Zhou et al. Application of hidden Markov models to analyze, group and visualize spatio-temporal COVID-19 data
Shahtori et al. Sequential monte carlo filtering estimation of ebola progression in west africa
An et al. A nonlinear model predictive control model aimed at the epidemic spread with quarantine strategy
Singham et al. Density estimation of simulation output using exponential epi-splines
CN116564542A (zh) 一种基于动力学模型的传染病传播规模预测方法及系统
McGowan et al. Quantifying uncertainty in infectious disease mechanistic models

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