引言
主要是大三的课设,结合的是之前做的SEIR同时相较之前加入了一些改进。
课程设计
学校的2020的校赛前两小问,大三菜鸡做这个大二大佬们的竞赛题,内心深觉羞愧,但还是厚着凉皮把这些东西做做整理。
建模过程
主要参考笔者的前一篇博文与知网的一些论文,做了一些结合和改进,最后作为自己的课程设计的结果。
建模的手写过程如下:
建模过程相较上次,做了一些小小的改进,(让他看起来像真的一样,其实还是假的y y出来的)
建模结果
根据上述的建模过程,在MatlabR2016b中进行实现:
- 首先是先看看自然条件的SEIR模型的模拟过程:
- 由模型不难看出,如果放着肺炎不管的话,按照模型的推算会有很多人患病和死亡。于是我们尝试加入隔离措施得到下图:
- 由于问题给出了附件数据,对附件尝试进行清洗,开始尝试对拟合,拟合方式主要为SEIR差分式的迭代和采取隔 离措施后的参数调整,主要调整的参数为接触率C与有效接触比例b作为(同时增加了治愈率)不同城市采取隔 离措施后的改变体现。(外部增援后续治疗跟进,采取隔 离后人与人间的接触减少同时又更多比例的人被抓去隔 离)
拟合效果如下: - 同时将模型推广,比如选取哈尔滨市做模型的推广结果进行尝试得到:
- 最后进行灵敏度的测试与分析,笔者这里主要对采取隔离的时间与强度进行相关的灵敏度分析,结果如下:
小结
因为菜,所以需要多敲代码,多做事,结论太过真实,导致自己活成消化。
一段
据说世界上最好的安慰,并不是告诉对方,一切都会好起来的。而是愁着脸说,哭个屁,你看,我比你还惨。