疫情隔离阶段果然是可以心无旁骛安心学习的时光。我在《用matlab对曲线拟合做最小二乘法单目标优化》https://blog.csdn.net/weixin_41855010/article/details/104225901 中查阅了2003年的SARS数据,根据当年香港的病例数据建立并且优化一波指数模型,爬了Z市卫计委的公开数据,用非典参数估计了2月H是的新增和累积,结论是按照现在趋势,2月底H市的累积确诊预计在250见顶。所以还是好好呆在家里吧,疫情应该很快就会过去的。

下面是代码:

clear all;clf;clc;
T=3; %平移量
T0=0.5;%半衰期
Peak=20;%峰值
t=1:5;
y1=1./(1+exp(-(t-T)./T0))*Peak;% plot(t,y1)
t=1:10
ry2=[19    18  16  13  12  8   14  9   10  5];
plot(t,ry2);
hold on;for t=1:30y2(t)=max(0,ceil(19*exp(-t/10.5)+4*(rand-0.5)));
end
figure(1)
t=1:30
plot(t,y2)
title('H市新冠肺炎新增确诊断数量预测图')
legend('真实值','预测值')total = [6  12  26  28  32  51  69  85  98  110 118 132 141 151 156]
added = [0  6  14  2   4   19  18  16  13  12  8   14  9   10  5]for i = 11:30added(5+i) = y2(i)total(5+i) = total(4+i) + y2(i)
endTime = length(added);
t=1:Time;date = datetime(2020,1,23) + caldays(1:Time);
grid on
% figure(1)
% plot(t,added)
figure(2)
plot(date,total)
hold on
xx1 = [5,12,19,26,33]; % x = [6.321,11.15,21.15]
tt1 = date(xx1);
yy1 = total( xx1 ); % 与x对应的y值
tDay = [737819,737826,737833,737840,737847]
DATE = ['Jan.29th:','Feb.5th','Feb.12th','Feb.19th','Feb.26th']
% 画点 标注
for i = 1:length(xx1)
plot(tt1(i),yy1(i),'rs','MarkerEdgeColor','k','MarkerFaceColor','g','MarkerSize',6)
text(tDay(i),yy1(i),[num2str(yy1(i))],'EdgeColor','red','BackgroundColor',[.7 .9 .7],'VerticalAlignment','bottom');
endtitle('H市新冠肺炎累计确诊断数量预测图')
xlabel('日期')
ylabel('累计确证人数预测量')
grid on

得到的结果如下:


这个预测准不准,2月底来看看验证一下。

用matlab对2003年香港SARS数据建模预估新冠病毒在H市的疫情走势相关推荐

  1. [数据科学]001.新冠病毒全球疫情-数据可视化-plotly-完整详尽教程

    -1 关于plotly 的小问题 因为plotly用着还不熟练,我把我遇到的一些plotly调用问题总结了一下: [打补丁]plotly常见问题 0. 前言: 疫情期间,工作'过于'繁忙...教程早已 ...

  2. 扎克伯格:Facebook不会把位置数据交给政府,哪怕是用于追踪新冠病毒!

    数据猿报道 最新招聘公告: ①阿里数据中台品牌团队杭州四大职位虚位以待 ②蚂蚁金服大数据部2020实习生招募全面启动 ③TalkingData北京招资深银行行业BD... ④数数科技上海招大数据运维支 ...

  3. 香港科技大学数据建模硕士学位项目(MSc DDM)武汉站宣讲会通知(Mix Mode 线上+线下同步进行)

    宣讲时间:2021年12月4日(星期六)15:00-16:30 宣讲地点:武汉君宜王朝大饭店三楼绿岛厅(地铁2号线街道口站,银泰创意城旁) 嘉宾教授:刘军伟教授 Prof.Junwei LIU htt ...

  4. 香港科技大学数据建模硕士学位项目(MSc DDM)宣讲会通知(2022年9月入学申请)

    拉近理论.技术和实践的距离,成就更有竞争力的数据建模人才  香港科技大学(2021年度QS全球顶尖大学排名27位) 2019年泰晤士高等教育 ·全球大学就业能力排名第10位(连续7年大中华区院校之首) ...

  5. 基于改进SEIR模型的病毒传播动力学建模与疫情预测分析(以COVID-19新冠病毒为例,超详细,带matlab源码)

    文章目录 前言 一.数学基础知识 符号定义 二.传统SEIR模型的建立与求解 1.经典的SEIR传播动力学模型建立 2.根据经典的SEIR模型列出微分方程并求解 三.SEIR模型第一次修正 1.模型建 ...

  6. 港科夜闻|香港科大新研究显示人体T细胞免疫反应可有效应对新冠病毒变异株Omicron...

    关注并星标 每周阅读港科夜闻 建立新视野 开启新思维 1.香港科大新研究显示人体T细胞免疫反应可有效应对新冠病毒变异株Omicron.透过接种新冠疫苗或经感染新冠肺炎后所产生的T细胞,能有效清除被感染 ...

  7. 哈佛大学通过百度搜索的数据写了篇论文:新冠病毒起源于去年8月的武汉...

    "我想上百度去找找数据写篇论文." "你哈佛大学的?" 从全球最好的大学之一,到民间口碑大翻车的"野鸡团队",哈佛大学只花了一晚上的时间. ...

  8. 基于跨国人口流动数据的新冠病毒传播途径分析

    来源:运筹OR帷幄 本文约5500字,建议阅读5分钟. 本文通过建立跨国航空网络数据基础上的空间-时间动态模型,分析过去三个月新冠病毒在国家人口流动与国家感染增长方面的关联. [ 导读 ]本文通过建立 ...

  9. [转载] python+selenium定时爬取丁香园的新冠病毒每天的数据,并制作出类似的地图(部署到云服务器)

    参考链接: Python vars() python+selenium定时爬取丁香园的新冠病毒每天的数据,并制作出类似的地图(部署到云服务器) 声明:仅供技术交流,请勿用于非法用途,如有其它非法用途造 ...

最新文章

  1. The Interface name: -- index:21 you checked seems not up.
  2. 关于Ubuntu 安装tftp服务器的问题解决
  3. html中盒子模型立体结构图
  4. android 开发环境以及hello world,Android Studio开发环境搭建和HelloWorld
  5. 【深度学习】深度学习实验记录--自编码+分类器
  6. python百度百科-Python语言的发展简史
  7. 极客大学产品经理训练营 产品经理新人如何落地 第0课总结
  8. 不能注册 DLL/OCX:Regsvr32失败的解决办法
  9. MATLAB与STK互联23:卫星对象操作(14)—仿真分析案例1,建立8个地面站,分析对卫星的可见性
  10. VS2015安装教程(带图解+下载地址+超详细)
  11. iso硬盘安装 凤凰os_在Linux中安装凤凰系统(Phoenix OS)的方法
  12. LeetCode 第35题 搜索插入位置 做题记录
  13. 解锁教爸妈使用智能手机的新方法!
  14. 多周期时间序列分解算法——MSTL原理
  15. ORACLE数据库数据类型
  16. Vue.js 技术揭秘
  17. 小程序wxs使用教程
  18. DSPE-PEG-GSH
  19. Bias/variance tradeoff
  20. git push -u参数是什么意思?

热门文章

  1. 高速PCB基础——电阻,电容,电感(第三讲)
  2. java解析webp格式图片宽高;java解析webp图片转png格式
  3. sklearn 中F1-score的计算
  4. unity 获取多个坐标点的中心点
  5. spm + host
  6. 五个无刷马达驱动电路分享!
  7. 关于币圈炒币 你必须要知道的风险
  8. 运维工作中常见的一些定律
  9. 说话技巧:成为谈话高手必不可少的十个特点
  10. 强化学习中的backups