CN113345596A - 一种基于seir模型的传染病系统的事件触发滤波方法 - Google Patents
一种基于seir模型的传染病系统的事件触发滤波方法 Download PDFInfo
- 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
Links
- 208000035473 Communicable disease Diseases 0.000 title claims abstract description 79
- 208000015181 infectious disease Diseases 0.000 title claims abstract description 69
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000001914 filtration Methods 0.000 title claims abstract description 23
- 230000008859 change Effects 0.000 claims abstract description 23
- 238000011161 development Methods 0.000 claims abstract description 7
- 238000013461 design Methods 0.000 claims description 28
- 230000001960 triggered effect Effects 0.000 claims description 20
- 230000008569 process Effects 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 11
- 201000010099 disease Diseases 0.000 claims description 7
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 claims description 7
- 238000009795 derivation Methods 0.000 claims description 6
- 238000002955 isolation Methods 0.000 claims description 6
- 238000012546 transfer Methods 0.000 claims description 6
- 238000012795 verification Methods 0.000 claims description 6
- 239000000654 additive Substances 0.000 claims description 4
- 230000000996 additive effect Effects 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 3
- 230000001629 suppression Effects 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- 208000025721 COVID-19 Diseases 0.000 description 5
- 238000010276 construction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 241000711573 Coronaviridae Species 0.000 description 2
- 206010035664 Pneumonia Diseases 0.000 description 2
- 238000011217 control strategy Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 241001678559 COVID-19 virus Species 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- 206010051379 Systemic Inflammatory Response Syndrome Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/80—ICT 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation 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型传染病系统中各类人群数量变化的模型,尤其涉及一种传染病的事件触发滤波方法。
背景技术
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型传染病的人群数量的状态空间模型,具体方法是:
其中,表示t时刻传染病的各类人群数量的变化率;表示在t时刻传染病的各类人群数量,n表示人群的类别数;代表传感器测量输出数量;表示期望输出数量;代表影响各类人群数量变化的外来人群等外在因素的干扰,m表示外在因素数量;函数rt为半马尔科夫跳变过程也即传染病预测系统的跳变信号,在有限集内取值,为系统权重矩阵,反应各类人群之间的影响,表示系统在t时刻的不确定摄动项,α(t)表示是服从伯努利分布的随机信号。 均为实数矩阵。记rt=i,,
步骤1.2设计半马尔科夫跳变信号rt,其转移速率满足以下条件:
步骤1.3设计步骤1.1中建立模型的加性不确定摄动ΔAi(t),具体形式如下:
步骤2、建立SEIR型传染病中各类人群良好发展的事件触发条件,其构建形式如下:
‖ey(t)‖1>β‖y(t)‖1 (4)
步骤3、设计估计SEIR型传染病中各类人群数量变化的事件触发滤波器,其具体构建如下:
步骤3.1滤波器系统设计:
步骤3.2给出事件触发滤波器如下:
其中,Gi表示传感器的测量输出。
步骤3.3基于步骤2.1和步骤3.1,设计滤波器误差系统如下:
步骤3.4设计预测SEIR型传染病滤波估计系统的条件,具体如下:
其中Φ=I+β1m×m,Ψ=I-β1m×m;
由此可见,在事件触发滤波器增益矩阵下,传染病滤波估计系统是正的且随机稳定的。
步骤4、SEIR型传染病系统正性验证过程如下:
根据步骤3.4的前五个条件可得:
步骤5、SEIR型传染病系统随机稳定性的验证过程如下:
步骤3.3中误差系统的上界:
构造一个随机余正Lyapunov函数如下:
对上述的Lyapunov函数取弱无穷小算子,具体如下:
下面分两个步骤讨论无/有外扰输入影响的传染病滤波估计系统的随机稳定性。
根据步骤3.4中的条件可得:
进而,可得:
从而可得:
步骤5.2当ω(t)≠0时,结合步骤3.4可得:
进一步,根据上述推导可得:
其中,E{·}表示数学期望,dτ表示时间变量算子、γ是扰动增益抑制水平,所述的传染病滤波估计系统是随机稳定的且具有l1增益性能γ。
本发明的有益效果如下:
本发明方法针对当前所述的SEIR型传染病系统中各类人群数量变化的情况,建立了SEIR型传染病系统的状态空间模型,借助矩阵论知识和构造的随机余正Lyapunov函数,设计了一类事件触发滤波器,并分析了其正性及稳定性,预测了传染病系统中各类人群数量变化的情况。
本发明提供了一种利用现代控制理论技术建立了传染病系统事件触发滤波估计的状态空间模型,设计了一类事件触发滤波器,分析了其正性及稳定性,降低了系统设计成本,预测了SEIR型传染病中各类人群数量的变化趋势。
附图说明
图1是本发明所述的SEIR传染病模型示意图。
具体实施方式
下面结合具体实施例来对本发明进行进一步说明。
本实例提供了一种预测SEIR型传染病系统中各类人群数量变化的事件触发滤波估计方法,其具体步骤如下:
步骤1、通过采集SEIR型传染病中各类人群数量的数据,利用采集到的数据结合隔离条件,建立SEIR型传染病的人群数量的状态空间模型:
具体构造如下步骤:
步骤1.1建立半马尔科夫跳变过程的传染病系统:
其中,表示t时刻传染病的各类人群数量的变化率;表示在t时刻传染病的各类人群数量,n表示人群的类别数;代表传感器测量输出数量;表示期望输出数量;代表影响各类人群数量变化的外来人群等外在因素的干扰,m表示外在因素数量;函数rt为半马尔科夫跳变过程也即传染病预测系统的跳变信号,在有限集内取值,为系统权重矩阵,反应各类人群之间的影响,表示系统在t时刻的不确定摄动项,α(t)表示是服从伯努利分布的随机信号。 均为实数矩阵。记rt=i,,
步骤1.2设计半马尔科夫跳变信号rt,其转移速率满足以下条件:
步骤1.3设计步骤1.1中建立模型的加性不确定摄动ΔAi(t),具体形式如下:
步骤2、建立SEIR型传染病中各类人群良好发展的事件触发控制条件,其构建形式如下:
步骤2、建立SEIR型传染病中各类人群良好发展的事件触发条件,其构建形式如下:
‖ey(t)‖1>β‖y(t)‖1 (4)
步骤3、设计估计SEIR型传染病中各类人群数量变化的事件触发滤波器,其具体构建如下:
步骤3.1滤波器系统设计:
步骤3.2给出事件触发滤波器如下:
其中,Gi表示传感器的测量输出。
步骤3.3基于步骤2.1和步骤3.1,设计滤波器误差系统如下:
步骤3.4设计预测SEIR型传染病滤波估计系统的条件,具体如下:
其中Φ=I+β1m×m,Ψ=I-β1m×m;
由此可见,在事件触发滤波器增益矩阵下,传染病滤波估计系统是正的且随机稳定的。
步骤4、SEIR型传染病系统正性验证过程如下:
根据步骤3.4的前五个条件可得:
步骤5、SEIR型传染病系统随机稳定性的验证过程如下:
步骤3.3中误差系统的上界:
构造一个随机余正Lyapunov函数如下:
对上述的Lyapunov函数取弱无穷小算子,具体如下:
下面分两个步骤讨论无/有外扰输入影响的传染病滤波估计系统的随机稳定性。
根据步骤3.4中的条件可得:
进而,可得:
从而可得:
步骤5.2当ω(t)≠0时,结合步骤3.4可得:
进一步,根据上述推导可得:
其中,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型传染病的人群数量的状态空间模型,具体方法是:
其中,表示t时刻传染病的各类人群数量的变化率;表示在t时刻传染病的各类人群数量,n表示人群的类别数;代表传感器测量输出数量;表示期望输出数量;代表影响各类人群数量变化的外来人群等外在因素的干扰,m表示外在因素数量;函数rt为半马尔科夫跳变过程也即传染病预测系统的跳变信号,在有限集内取值,为系统权重矩阵,反应各类人群之间的影响,表示系统在t时刻的不确定摄动项,α(t)表示是服从伯努利分布的随机信号; 均为实数矩阵;记rt=i,,
步骤1.2设计半马尔科夫跳变信号rt,其转移速率满足以下条件:
步骤1.3设计步骤1.1中建立模型的加性不确定摄动ΔAi(t),具体形式如下:
4.根据权利要求3所述的一种基于SEIR模型的传染病系统的事件触发滤波方法,其特征在于步骤3所述的设计估计SEIR型传染病中各类人群数量变化的事件触发滤波器,其具体构建如下:
步骤3.1滤波器系统设计:
步骤3.2给出事件触发滤波器如下:
其中,Gi表示传感器的测量输出;
步骤3.3基于步骤2.1和步骤3.1,设计滤波器误差系统如下:
步骤3.4设计预测SEIR型传染病滤波估计系统的条件,具体如下:
其中Φ=I+β1m×m,Ψ=I-β1m×m;
由此可见,在事件触发滤波器增益矩阵下,传染病滤波估计系统是正的且随机稳定的。
6.根据权利要求5所述的一种基于SEIR模型的传染病系统的事件触发滤波方法,其特征在于步骤5所述的SEIR型传染病系统随机稳定性的验证过程如下:
步骤3.3中误差系统的上界:
构造一个随机余正Lyapunov函数如下:
对上述的Lyapunov函数取弱无穷小算子,具体如下:
下面分两个步骤讨论无/有外扰输入影响的传染病滤波估计系统的随机稳定性:
根据步骤3.4中的条件可得:
进而,可得:
从而可得:
步骤5.2当ω(t)≠0时,结合步骤3.4可得:
进一步,根据上述推导可得:
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)
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 | 云境商务智能研究院南京有限公司 | 一种受欺骗攻击的非线性网络化系统状态估计方法 |
-
2021
- 2021-06-24 CN CN202110703319.0A patent/CN113345596B/zh active Active
Patent Citations (3)
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)
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 |