SEIR模型案例_SOR模型

SEIR模型案例_SOR模型SEIR模型案例原理说明https://zhuanlan.zhihu.com/p/142117573%N是群体总样本数E=0;%潜伏期人数I=1;%患者人数S=N-I;%易感者人数

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE使用 1年只要46元 售后保障 童叟无欺

SEIR模型案例

原理说明

https://zhuanlan.zhihu.com/p/142117573

                               %N是群体总样本数
E=0;                           %潜伏期人数
I=1;                           %患者人数
S=N-I;                         %易感者人数
R=0;                           %康复者人数
r=20;                          %每个患者平均每天能接触的人数
B=0.03;                        %患者每天接触的人中被感染的概率
a=0.1;                         %潜伏期装变为患者的概率
y=0.1;                         %患者痊愈概率
r2=20;                         %处于潜伏期的人每天接触的人数

B2=0.03;                       %潜伏期的人每天接触的人中被感染的概率
T=1:140;                       %模拟从第一天到第140天的状况

for idx =1:length(T)-1         %循环遍历
    S(idx+1) = S(idx)- r*B*S(idx)*I(idx)/N - r2*B2*S(idx)*E(idx)/N;
    E(idx+1) = E(idx) + r*B*S(idx)*I(idx)/N - a*E(idx) + r2*B2*S(idx)*E(idx)/N;
    I(idx+1) = I(idx) + a*E(idx) - y*I(idx);
    R(idx+1) = R(idx) + y*I(idx);
end 

%{迭代计算每天的感染情况 idx代表天数
	S(idx+1) = S(idx)- r*B*S(idx)*I(idx)/N -r2*B2*S(idx)*E(idx)/N;
	这里S(idx+1)代表第idx+1天的易感者人数 
	其数值等于第idx天的易感者人数减去由于 潜伏者接触易感染者所造成的感染人数 和 患者接触易感染者所造成的感染人数
    r*I(idx)表示所有患者每天接触的人的总数  S(idx)/N表示所有人中易感染者占比
    r*B*S(idx)*I(idx)/N 即所有患者一天因为易感染者所造成的感染数
    同理r2*B2*S(idx)*E(idx)/N 即所有潜伏者一天因为易感染者所造成的感染数
    
    E(idx+1) = E(idx) + r*B*S(idx)*I(idx)/N - a*E(idx) + r2*B2*S(idx)*E(idx)/N;
    第idx+1天的潜伏者人数等于 
    第idx天的潜伏者人数 加上 潜伏者接触易感染者所造成的感染人数 和 患者接触易感染者所造成的感染人数
    减去由潜伏者转化为患者的人数
    
    y*I(idx)第idx天的患者痊愈人数
    
    以下分析同理
%}

plot(T,S,T,E,T,I,T,R);grid on;               %以下是绘图部分
xlabel('天');ylabel('人数')
legend('易感者','潜伏者','传染者','康复者')

在这里插入图片描述

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/166783.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)


相关推荐

  • 听说B站视频的完成率也是一项挺重要的指标,下面来提高它

    听说B站视频的完成率也是一项挺重要的指标,下面来提高它

    2021年11月11日
  • f1 score是什么_F1值

    f1 score是什么_F1值F1score是一个平均数;对精确率与召回率进行平均的一个结果;平均算法有四个,如图所示:调和平均数:Hn=n/(1/a1+1/a2+…+1/an)几何平均数:Gn=(a1a2…an)^(1/n)算术平均数:An=(a1+a2+…+an)/n平方平均数:Qn=√[(a1^2+a2^2+…+an^2)/n]这四种平均数满足Hn≤Gn≤An≤QnF…

    2022年10月14日
  • 为什么程序员做外包会被瞧不起?

    为什么程序员做外包会被瞧不起?二哥,有个事想询问下您的意见,您觉得应届生值得去外包吗?公司虽然挺大的,中xx,但待遇感觉挺低,马上要报到,挺纠结的。

  • 短视频创作的技巧是什么_短图文创作特点

    短视频创作的技巧是什么_短图文创作特点现在短视频越来越受到大众的喜爱,大概现在每个人坐车休假吃饭都在拿着手机刷着短视频,可见现在短视频对于现在的人来说还是挺普遍的,那么很多人都想从事短视频行业应该如何去进行创作呢,下面就和大家分享平时我会用到的一些小技巧。构思框架在做短视频的时候一定不要想着能够一夜爆火,当然如果你的作品足够优质,那也不排除这样的可能,首先需要你先考虑的是各种因素,主题、定位和内容连贯性,还有视觉效果。在确定主题后,要做好计划,如拍摄方向、表达形式。时间一定要把握住短视频的时长,因为现在短视频推送都是讲究一个完播

  • airplay影像_播放ftp服务器的视频

    airplay影像_播放ftp服务器的视频http://bbs.weiphone.com/read-htm-tid-1785042.htmliOS客户端软件:AirPlayer[支持RMVB,MKV,MP4,AVI等等格式],另外有个PulgPlayer不支持RMVBMKV等等格式。PC端媒体服务器:WIN7自带的媒体中心(论坛有设置方法的教程,这不介绍了)…

  • 如何查看linux服务器cuda_[882]如何查看windows和linux的CUDA版本

    如何查看linux服务器cuda_[882]如何查看windows和linux的CUDA版本linux查看cuda版本cuda一般安装在/usr/local/cuda/路径下,该路径下有一个version.txt文档,里面记录了cuda的版本信息cat/usr/local/cuda/version.txt如何查看windows的CUDA版本CUDA是什么?CUDA(ComputeUnifiedDeviceArchitecture),是显卡厂商NVIDIA推出的运算平台。CU…

    2022年10月21日

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号