SIR模型实现(matlab)
时间:2020-04-12
本文章向大家介绍SIR模型实现(matlab),主要包括SIR模型实现(matlab)使用实例、应用技巧、基本知识点总结和需要注意事项,具有一定的参考价值,需要的朋友可以参考一下。
matlab代码
clc clear close all; A = 0.4; B = 0.1; I = 0.4; S = 0.5; %ode tspan = [0 50]; y0 = [I S]; [t, y] = ode45(@(t,y)odefun(t,y,A,B), tspan, y0); r = 1-y(:,1)-y(:,2); %euler n = size(r,1); h = 50 / (n-1); t_0 = [0:h:50]'; y_i = zeros(n,1); y_s = zeros(n,1); y_i(1) = I; y_s(1) = S; for i = 1:n-1 y_i(i+1) = h*[A*y_i(i)*y_s(i) - B*y_i(i)]+y_i(i); y_s(i+1) = h*[-A*y_i(i)*y_s(i)]+y_s(i); end r_0 = 1 - y_i(:,1) - y_s(:,1); %画图 subplot(2,2,1); plot(t,y(:,1),'-o',t,y(:,2),'-.',t,r,'g'); hold on; legend('生病人数:i(t)','健康人数:s(t)','移除人数:r(t)','Location','Best'); ylabel('占人口比例%'); xlabel('时间t'); str = ['接触数λ/μ:',num2str(A/B),' 初始生病人数:',num2str(I),',初始健康人数:',num2str(S)]; text(15,0.4,str,'FontSize',10); title('SIR模型(ode)'); subplot(2,2,2); plot(t_0,y_i,'-o',t_0,y_s,'-.',t_0,r_0,'g'); hold on; legend('生病人数:i(t)','健康人数:s(t)','移除人数:r(t)','Location','Best'); ylabel('占人口比例%'); xlabel('时间t'); str = ['接触数λ/μ:',num2str(A/B),' 初始生病人数:',num2str(I),',初始健康人数:',num2str(S)]; text(15,0.4,str,'FontSize',10); title('SIR模型(euler)'); subplot(2,2,3); plot(t_0,y_i,'r-',t,y(:,1),'-.'); diff = sum(abs(y_i - y(:,1))); str1 = ['生病人数对比图i(t), 误差:',num2str(diff)]; title(str1); legend('euler','ode','Location','Best'); ylabel('占人口比例%'); xlabel('时间t'); subplot(2,2,4); plot(t_0,y_s,'r-',t,y(:,2),'-.'); diff = sum(abs(y_s - y(:,2))); str1 = ['健康人数对比图s(t), 误差:',num2str(diff)]; title(str1); legend('euler','ode','Location','Best'); ylabel('占人口比例%'); xlabel('时间t'); function dydt = odefun(t,y,A,B) dydt = zeros(2,1); dydt(1) = A*y(1)*y(2) - B*y(1); dydt(2) = -A*y(1)*y(2); end
结果
原文地址:https://www.cnblogs.com/shish/p/12685865.html
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法
- Crawlab准备之python+scrapy环境搭建
- 一致性协议之 ZAB
- Elasticsearch Analyzer原理分析并实现中文分词
- php laravel5.5使用rabbitmq消息队列
- 还不知道怎么运行Python代码,快来撸代码
- Linux常用命令速查表
- 机器学习中必须要会用的模块
- 【MathorCup】2020年 A题 无车承运人平台线路定价问题,特征间的相关性分析
- 决策树预测成本价 2020年 MathorCup数学建模 A题 无车承运人平台线路定价问题
- 2020年 MathorCup数学建模 A题 无车承运人平台线路定价问题,定价评价
- RabbitMQ 消费端限流、TTL、死信队列
- 全栈必备 Redis基础
- JSP的异常处理
- (二)Mybatis-Plus的安装配置
- 如何在JSP里自定义标签