matlab病毒扩散,[MATLAB数学相关] 现有Allen的病毒扩散模型SIRS的马尔科夫过程程序,求指导...

本文主要是介绍matlab病毒扩散,[MATLAB数学相关] 现有Allen的病毒扩散模型SIRS的马尔科夫过程程序,求指导...,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

icon1.gif 现有Allen的病毒扩散模型SIRS的马尔科夫过程程序,求指导

由Allen编写的SIRS模型的马尔可夫过程如下:

for k=1:4 % 解马尔科夫链的样本轨迹

clear t s i r

t(1)=0; i(1)=i0; s(1)=N-i0;r(1)=0;

j=1;

while i(j)>0 && t(j)

nn=s(j)+i(j)+r(j);

u1=rand;u2=rand;

den=((beta/nn)*i(j)*s(j)+gam*i(j)+alpha*i(j)+delta*r(j));%所有patch的概率之和

t(j+1)=-log(u1)/den+t(j); % Time to next event

e1=(beta/nn)*s(j)*i(j)/den;%以下是每个patch概率与总概率之比

e2=e1+gam*i(j)/den;

e3=e2+alpha*i(j)/den;

e4=e3+delta*r(j)/den;

if (u2<=e1) %每个patch情况下的个体变化情况

s(j+1)=s(j)-1;

i(j+1)=i(j)+1;

r(j+1)=r(j);

elseif (u2>e1 && u2<=e2)

s(j+1)=s(j);

i(j+1)=i(j)-1;

r(j+1)=r(j)+1;

elseif (u2>e2 && u2<=e3)

s(j+1)=s(j);

i(j+1)=i(j)-1;

r(j+1)=r(j);

else

s(j+1)=s(j)+1;

i(j+1)=i(j);

r(j+1)=r(j)-1;

end

j=j+1;

end

(省略出图部分)

end

本人菜鸟,对上述标黄部分的分支概率不太理解,变量“den”是各个转换概率之和,e1-e4分别是四种转换概率与“den”的比值,请问各位大神,e1-e4的顺序有区别么?另附本程序的转换图,见“SIR模型”图,若想在本程序的基础上,改写成“SIRP模型”图所应用的程序,应该怎么写呢?

目前存在两个问题:1.S→I之间有两种转换概率的情况该怎么解决呢?

2. P作为独立变量,对S→I有影响,但其增加或减少的变动情况与S,I,R相对独立,P的每一步数量应该单独统计么?

感谢各位!论文卡在这里很久了,请大家帮忙。

SIR模型.jpg

SIRP模型.jpg

这篇关于matlab病毒扩散,[MATLAB数学相关] 现有Allen的病毒扩散模型SIRS的马尔科夫过程程序,求指导...的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/397446

相关文章

Python的time模块一些常用功能(各种与时间相关的函数)

《Python的time模块一些常用功能(各种与时间相关的函数)》Python的time模块提供了各种与时间相关的函数,包括获取当前时间、处理时间间隔、执行时间测量等,:本文主要介绍Python的... 目录1. 获取当前时间2. 时间格式化3. 延时执行4. 时间戳运算5. 计算代码执行时间6. 转换为指

将Java程序打包成EXE文件的实现方式

《将Java程序打包成EXE文件的实现方式》:本文主要介绍将Java程序打包成EXE文件的实现方式,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录如何将Java程序编程打包成EXE文件1.准备Java程序2.生成JAR包3.选择并安装打包工具4.配置Launch4

Java程序进程起来了但是不打印日志的原因分析

《Java程序进程起来了但是不打印日志的原因分析》:本文主要介绍Java程序进程起来了但是不打印日志的原因分析,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录Java程序进程起来了但是不打印日志的原因1、日志配置问题2、日志文件权限问题3、日志文件路径问题4、程序

PyInstaller打包selenium-wire过程中常见问题和解决指南

《PyInstaller打包selenium-wire过程中常见问题和解决指南》常用的打包工具PyInstaller能将Python项目打包成单个可执行文件,但也会因为兼容性问题和路径管理而出现各种运... 目录前言1. 背景2. 可能遇到的问题概述3. PyInstaller 打包步骤及参数配置4. 依赖

SpringBoot实现微信小程序支付功能

《SpringBoot实现微信小程序支付功能》小程序支付功能已成为众多应用的核心需求之一,本文主要介绍了SpringBoot实现微信小程序支付功能,文中通过示例代码介绍的非常详细,对大家的学习或者工作... 目录一、引言二、准备工作(一)微信支付商户平台配置(二)Spring Boot项目搭建(三)配置文件

Spring Security基于数据库的ABAC属性权限模型实战开发教程

《SpringSecurity基于数据库的ABAC属性权限模型实战开发教程》:本文主要介绍SpringSecurity基于数据库的ABAC属性权限模型实战开发教程,本文给大家介绍的非常详细,对大... 目录1. 前言2. 权限决策依据RBACABAC综合对比3. 数据库表结构说明4. 实战开始5. MyBA

将Mybatis升级为Mybatis-Plus的详细过程

《将Mybatis升级为Mybatis-Plus的详细过程》本文详细介绍了在若依管理系统(v3.8.8)中将MyBatis升级为MyBatis-Plus的过程,旨在提升开发效率,通过本文,开发者可实现... 目录说明流程增加依赖修改配置文件注释掉MyBATisConfig里面的Bean代码生成使用IDEA生

C# WinForms存储过程操作数据库的实例讲解

《C#WinForms存储过程操作数据库的实例讲解》:本文主要介绍C#WinForms存储过程操作数据库的实例,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录一、存储过程基础二、C# 调用流程1. 数据库连接配置2. 执行存储过程(增删改)3. 查询数据三、事务处

JSON Web Token在登陆中的使用过程

《JSONWebToken在登陆中的使用过程》:本文主要介绍JSONWebToken在登陆中的使用过程,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录JWT 介绍微服务架构中的 JWT 使用结合微服务网关的 JWT 验证1. 用户登录,生成 JWT2. 自定义过滤

java中使用POI生成Excel并导出过程

《java中使用POI生成Excel并导出过程》:本文主要介绍java中使用POI生成Excel并导出过程,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录需求说明及实现方式需求完成通用代码版本1版本2结果展示type参数为atype参数为b总结注:本文章中代码均为