- 博客(51)
- 收藏
- 关注
原创 手搓传染病模型(SEIRD)
SEIRD模型是传染病研究中的重要工具,通过MATLAB实现其参数拟合与动态预测,可以更深入地理解疾病传播规律。模型通过微分方程组描述易感者、暴露者、感染者、康复者和死亡者五类人群的动态变化,并考虑自然出生率和死亡率,使模拟更贴近现实。参数拟合是关键步骤,通过调整传染率等参数,使模型预测与实际数据最佳匹配。模型中的关键参数如传染率、潜伏期转化率、康复率和病死率直接影响预测结果,通过调整这些参数可以模拟不同防控措施的效果。该模型具有重要的实际应用价值,如疾病预测、政策评估和资源分配。通过MATLAB的数值计算
2025-05-10 00:45:00
231
原创 MATLAB制作柱状图与条图:数据可视化的基础利器
柱状图和条图都是用来表示的常见图表形式,它们的核心目的是通过矩形的长度来比较各类别的数值大小。其实就是,它们的本质是一样的:用矩形的长度表示数值大小,不同的只是:也叫,是一种竖直排列的图表。每个类别对应一个竖直的矩形柱,柱子的高度表示该类别的数值大小。柱状图常用于比较不同类别的数据,或展示某一指标在不同时间点的变化。:也称为,是一种横向排列的图表。每个类别对应一个水平的矩形条,条的长度表示该类别的数值大小。条图在类别名称较长或类别较多时更易阅读,适合显示排名或比较数据。
2025-05-10 00:15:00
1399
原创 手搓传染病模型(SEIS)
代码中首先定义了一个嵌套函数SEISS = y(1);E = y(2);I = y(3);dydt = [dS;dE;dI];enddS:易感者的变化率。表示易感者被感染成为暴露者的速率,+ mu * I表示感染者恢复后再次成为易感者(SEIS 模型的关键)。dE:暴露者的变化率。是新增的暴露者,是暴露者转为感染者的速率。dI:感染者的变化率。omega * E是暴露者转为感染者,表示感染者恢复后成为易感者(占比1 - alpha),表示感染者永久免疫或死亡(占比alpha。
2025-05-09 01:00:00
741
原创 手搓传染病模型(SEI)
β = 0.8:传染率,反映感染者将疾病传播给易感者的能力。数值越高,传播越快。ω = 0.1:潜伏期转化率,即暴露者转为感染者的速率。
2025-05-08 02:30:00
1094
原创 MATLAB小试牛刀系列(5)
市场上有 n 种资产 (如股票、债券等)\(S_{i}(i = 1,2,\cdots,n)\)供投资者选择,某公司有数额为 M 的一笔相当大的资金可用作一个时期的投资。公司财务分析人员对这 n 种资产进行了评估,估算出在这一时期内购买资产\(S_{i}\)的平均收益率为\(r_{i}\),并预测出购买\(S_{i}\)的风险损失率为\(q_{i}\)。考虑到投资越分散,总的风险越小,公司确定,当用这笔资金购买若干种资产时,总体风险可用所投资的\(S_{i}\)中最大的一个风险来度量。
2025-05-08 01:30:00
727
原创 MATLAB制作直方图
统计每个区间内数据的个数理解数据的实际数量分布适合数据量不太大时用比较不同样本规模的数据分布强调数据比例而非具体数量避免样本量影响解读估计概率密度函数分布拟合或对比分析用于统计分析前的可视化检验对比不同组数据的分布展示组间的相似或差异可与透明度叠加更清晰呈现你想了解什么?用哪种直方图?数据在各区间的数量分布频数直方图各区间数据所占比例频率直方图数据的概率密度估计概率密度直方图对比两组或多组数据的分布分组直方图数据分布整体趋势,是否近似正态分布概率密度直方图或频率直方图。
2025-05-07 13:06:06
1242
原创 从 txt 到 mat:脑电数据批量提取的高效实战攻略
这套组合拳能帮你高效处理脑电数据,告别手动查找的麻烦!无论是科研还是项目实战,都能大幅提升效率。赶紧试试,让数据处理飞起来~ 有任何问题,欢迎评论区交流,一起解锁更多数据处理技巧!💻✨。
2025-05-07 03:15:00
228
原创 手搓SIP模型
SIP 模型主要涉及三个状态变量:S、I 和 P。在这个模型中,这些变量随时间的变化遵循一定的动力学规则。r = 2;alpha = 1;beta = 1;c = 1;d = 0.25;mu = 0.4;% [0.4,0.8] 之间尝试不同值r:代表系统的某个基础增长率。alpha和beta:用于描述不同状态变量之间的相互作用系数。c和d:与状态变量 P 的变化相关。mu:一个关键参数,可在[0.4, 0.8]范围内取值,它会影响系统的动态行为。
2025-05-06 16:41:53
294
原创 MATLAB制作折线图:从入门到精通(下)
展示跨年级别的数据趋势分析宏观变化(气候、经济、人口等)数据密集时更突出整体走势数学建模类比赛遇到超多数据时,也可以看情况选择此类折线图两组数据量纲差异大,无法放在同一Y轴强调变量间的同步性或关联性避免图形失真或误导数据包含“总量”和“比例”两种维度同时可视化数量和趋势让报表更直观、信息更丰富你手上有什么数据?用哪种折线图?跨度很长的时间序列(例:20年气温变化)超长时间序列折线图两种单位不同的数据(例:广告支出 vs 销售额)双坐标轴折线图。
2025-05-06 00:30:00
1117
原创 手搓传染病模型(SEIAVR)
随着时间推进,易感人群的数量逐渐减少并趋于平稳。这表明在传染病传播过程中,通过感染、隔离措施、疫苗接种以及人群自身免疫转化等方式,大多数人脱离了易感状态,进入感染、免疫等其他状态。当剩余的易感人群数量不再大幅下降时,意味着传染病在人群中的传播基本得到控制,传播态势趋于稳定,进入相对缓和的阶段。
2025-05-05 00:30:00
615
原创 MATLAB制作折线图:从入门到精通(上)
展示单变量的变化趋势分析时间序列数据直观显示某指标的上升或下降多组数据的趋势对比强调不同数据集在同一时间轴上的变化可视化多个实验或项目的数据变化既要看趋势又要看具体数据点数据点之间的差异很重要时想让数据展示更醒目你手上有什么数据?用哪种折线图?一组数据(例:某产品月销量)单线折线图多组数据(例:两家公司的股价走势)多线折线图一组数据,且需要清楚标出每个点(例:10次跑步成绩)带标记的折线图。
2025-05-05 00:15:00
1517
原创 MATLAB制作箱线图:箱线图类型与应用全解
比较多个组的数据分布:观察不同组之间是否存在显著的分布差异。分析多个组的中位数、四分位数差异:对比各组的集中趋势和离散程度。查看异常值差异:比较不同组之间的异常值数量和分布情况。比较不同类别的数据分布:查看不同类别下连续变量的分布差异。观察类别之间的集中趋势与分散程度:对比不同类别的中位数和四分位数。识别每个类别的异常值:通过箱线图的离群点标记,识别每个类别的异常数据。类目标签较长时的展示:当类别名称很长时,水平箱线图避免了标签重叠。横向比较不同组数据。
2025-05-04 07:00:00
1413
原创 MATLAB小试牛刀系列(4)
解 做变量变换 \(u_i=\frac{x_i + |x_i|}{2}\),\(v_i=\frac{|x_i| - x_i}{2}\),\(i = 1,2,3,4\),记 \(u = [u_1,u_2,u_3,u_4]^T\),\(v = [v_1,v_2,v_3,v_4]^T\),则可把模型变换为线性规划模型。求得最优解 \(x_1 = -2\),\(x_2 = x_3 = x_4 = 0\),最优值 \(z = 2\)。分析完毕,我们开始手搓。
2025-05-04 03:00:00
412
原创 手搓传染病模型(SEICR)
模型参数N = 100000;% 总人数I0 = 36;% 初始急性期感染人数C0 = 0;% 初始慢性期感染人数R0 = 0;% 初始康复人数E0 = 100;% 初始暴露人数S0 = 99864;% 初始易感人数% 暴露转化速率p = 0.7;% 感染者进入慢性期的比例% 康复速率f = 0.15;% 病死率% 人口自然出生率% 人口自然死亡率% 模拟天数% 传染率lb = 0.01;ub = 1;
2025-05-03 10:35:00
779
原创 手搓传染病模型(SVEIR)
SVEIR 模型扩展了经典的 SIR 模型,通过增加暴露者和疫苗接种者的状态,更贴合现实中添加疫苗措施的情况。模型能有效模拟传染病潜伏期和康复期对疫情的影响。
2025-05-03 00:30:00
1072
原创 MATLAB绘制局部放大图
这段代码展示了如何将技术性能比较以直观且专业的方式呈现,对于控制系统研究或数据可视化爱好者而言,具有很好的参考价值,同时也为MATLAB图表定制提供了实用范例。今天,我将分享一段 MATLAB 代码,该代码生成了一个主副图结合的可视化展示,用于比较不同控制系统性能表现。
2025-05-02 13:03:12
528
原创 MATLAB制作散点图:从基础到进阶的三种类型讲解
观察趋势和相关性检查是否有异常值做初步的数据探索分类变量(性别、地区、类型)+ 两个数值变量想看不同群体的差异或分布趋势非常适合做对比分析同时分析三个数值变量适合探索复杂的多变量关系更适合有交互式可视化支持的场景你手上有什么数据?用哪种散点图?两个数字(例:身高、体重)基本二维散点图两个数字 + 一个分类(例:性别)分组着色散点图三个数字(例:身高、体重、年龄)三维散点图。
2025-05-02 00:15:00
1753
原创 手搓传染病模型(SEIAR)
初始阶段,虽仅10人感染,但随着病毒迅速传播,暴露者(E)人数在第0天达到20人,成为疫情的“前哨”。随着疫情的持续,康复者(R)人数在第22天达到25人,预示着群体免疫屏障的初步形成。SEIAR 模型中,各个字母分别对应不同人群:S 代表易感者,E 代表暴露者,I 代表显性感染者,A 代表隐性感染者,R 代表康复者。在模型动态运行中,S 仓室会以单位时间 betaS(I + κA) 的速率向 E 仓室流入,而 E 仓室同样以单位时间 betaS(I + κA) 的速率向 E 仓室流入。
2025-05-01 11:42:59
297
原创 CodeArts for Cangjie环境配置(全网最全)
static,静态库,是一组预先编译好的目标文件的集合;dynamic,动态库,是一种在程序运行时才被加载到内存中的库文件,多个程序共享一个动态库副本,而不是像静态库那样每个程序都包含一份完整的副本。Attach,针对正在运行的程序,附加到已启动的程序。可以在顶部看到调试按钮,在左侧看到调试过程中的变量、监视、调用堆栈、断点信息。4.注册好之后回到第2步的页面,再次点击就会跳转到下载界面。7.这里的SDK位置,默认是没有的,需要下载,点击链接。按钮,首次会提示配置调试,点击。代码中设置断点,再次点击。
2025-05-01 00:15:00
570
原创 MATLAB绘制饼图(二维/三维)
在数据分析与展示领域,饼图是一种直观且高效的可视化工具,能够在瞬间传递各部分与整体的比例关系。今天,我将分享一段 MATLAB 绘制二维及三维饼图的代码,助你轻松将数据以饼图形式呈现于众人眼前。无论是二维饼图的简洁明了,还是三维饼图的立体生动,这段代码都能满足你的需求。
2025-04-30 15:33:30
388
原创 MATLAB小试牛刀系列(3)
表示产地\(A_i\)运到销地\(B_j\)的量,\(c_{ij}\)表示产地\(A_i\)到销地\(B_j\)的单位运价,\(d_j\)表示销地\(B_j\)的需求量,\(e_i\)表示产地\(A_i\)的产量。(1) 需求量约束,\(B_j\)销地的需求量等于所有产地运到\(B_j\)销地的运量和,即。(2) 产量约束,\(A_i\)产地运到所有销地的运量和少于等于该地的产量,即。单位商品运价如表所示。解 这是一个运输问题,总的产量大于总的需求量,是满足供应的运输问题。目标函数是使总的运费最小化,即。
2025-04-30 10:59:02
622
原创 手搓传染病模型(SEIAS)
在曲线基本平行于坐标轴的部分,表示对应人群数量变化极其缓慢,疫情进入稳定期或衰退期,此时新感染者的出现基本停止,或感染者数量趋近于零。话不多说,我们开始手搓。
2025-04-30 09:46:54
219
原创 MATLAB小试牛刀系列(2)
捷运公司在下一年度 1 - 4 月的 4 个月内拟租用仓库堆放物资。仓库租借费用随合同期而定,期限越长,折扣越大,具体数字见表 1.1。租借合同每月初都可办理,每份合同具体规定租用面积和期限。因此,该公司可根据需要,在任何一个月初办理租借合同。每次办理时可签一份合同,也可签若干份租用面积和租借期限不同的合同,试确定该公司签订租借合同的最优决策,目的是使所付租借费用最小。这个模型中的约束条件分别表示当月初签订的租借合同的面积加上该月前签订的未到期的合同的租借面积总和,应不少于该月所需的仓库面积。
2025-04-29 11:03:26
564
原创 手搓传染病模型(SEIR-拓展)
我们在今天的拓展SEIR模型中,通过添加干预措施(如戴口罩),模拟了戴口罩人群比例为0.00%, 20%, 40%, 60%, 80%的情况下,该传染病流行趋势的变化,从而更加直观的得到不同程度干预措施下的疾病负担变化情况和干预效果,为政策制定提供参考。在干预措施(戴口罩)力度不断加强的影响下,E和I人群的人数峰值会降低,也会延缓峰值的到来。但同时,传染病也会保持一段更长时间的低水平流行。话不多说,我们开始手搓。
2025-04-29 09:32:28
115
原创 MATLAB小试牛刀系列(1)
某机床厂生产甲、乙两种机床,每台机床销售后的利润分别为 4 千元与 3 千元。生产甲机床需用 A、B 机器加工,加工时间分别为每台 2h 和每台 1h;生产乙机床需用 A、B、C 三种机器加工,加工时间均为每台 1h。若每天可用于加工的机器时数分别为 A 机器 10h、B 机器 8h 和 C 机器 7h,问该厂应生产甲、乙机床各几台才能使总利润最大?假设该企业的产品不存在积压,即产量等于销量。题目来源:《数学建模算法与应用》——司守奎。分别表示甲、乙机床每天的产量。
2025-04-28 14:12:46
380
原创 手搓传染病模型(SEIR-拓展)
这里只做展示,小伙伴们可以根据自己的需要调整误差啥的。上述代码中的forecast函数见下面。今天我们将根据真实数据拟合出最优参数。我们在前文手搓了最基本的SEIR模型。话不多说,我们直接开始手搓。以传染率(β)为例。
2025-04-28 13:53:25
133
原创 手搓传染病模型(SEIR)
在本模型中,人群有四种自然史状态:易感者(S),暴露者(E),感染者(I)以及康复者(R) 2.模型假设人群分布是同质均匀的,未考虑人群出生、死亡、迁入迁出对疾病传播的影响 3.康复者永久免疫:康复者永久免疫该传染病,无再感染风险。
2025-04-27 15:26:03
207
原创 手搓传染病模型(SEIA-拓展)
2.暴露者转化为有症状感染者的概率( p )可能随时间变化。继上文手搓SEIA基础模型后,想要加两个创新点。1.传播过程中加入干预措施。我们话不多说,开始手搓。这里有一个难点,就是。
2025-04-27 14:50:39
203
原创 全网最全:osgeo的安装教程汇总
由于小编的python是基于conda环境的,所以在cmd中,和在Anacoda Prompt中所支持的版本是不一样的。但是还是不行,但是隐隐约约感觉到可能还是解释器的原因,就新建了一个系统解释器,链接到.venv文件下的python.exe。作者有提到:一定要根据python的版本去下载相对应的,但是小编并看不懂这些whl文件的名称,所以参考了。2.继续按照核心教程的方法,在后续过程中就可以成功的安装GDAL了。在查看可支持的版本的过程中,需要在cmd中键入。在复刻该教程的过程中,出现的一些问题。
2025-04-25 05:00:00
1004
原创 手搓种族相互依存模型
现在只考虑两个种群的情况,甲乙两种群的相互依存有三种形式:1) 甲可以独自生存,乙不能独自生存(或者反过来);2) 甲乙均可以独自生存;3) 甲乙均不能独自生存。我们讨论第一个情况:甲可以独自生存,乙不能独自生存。下面我们分别对这三种形式进行讨论。
2025-04-22 14:52:20
102
原创 全网最全:Tensorflow-gpu安装教程汇总
博主推荐install Spyder来测试,由于本人并不会使用Spyder,所以小伙伴们可以不必在Spyder中去输入相关代码,我的方法是简单粗暴的在。1.我安装的Anaconda所适配的python版本是3.10,第一遍安装的时候,没有把环境变量中python的path删除干净,导致在cmd中键入。在运行Anaconda Navigator的时候,会不断的出现prompt的弹窗,一定不要关闭。这一步,我愿称之为最关键的一步,关键点不仅仅是图中框起来的,还有。在复刻该教程的过程中,出现的一些问题。
2025-04-21 08:24:32
1432
湖北省乡界(shp)-可直接导入Arcgis
2025-04-29
辽宁省村界(shp)-可直接导入Arcgis
2025-04-29
Emilia Vynnycky, Richard White - An Introduction to Infectious Disease Modelling (2010, OUP Oxford)
2025-04-14
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人