clear;clc
load data30.mat
global true_s true_i true_r
true_s = p30(:,1);
true_i = p30(:,2);
true_r = p30(:,3);
% plot(1:30,true_s,'r-',1:30,true_i,'b-',1:30,true_r,'k-')
% legend('S','I','R')
% hold on
% lb = [0.1 0.03]; ub = [0.6 0.2];
% 粒子群
% options = optimoptions('particleswarm','Display','iter','SwarmSize',200,'PlotFcn','pswplotbestf');
% [x, fval] = particleswarm(@sir_Obj_fun,2,lb,ub,options)
% fmincom函数
% x0 = [0.1 0.1]
% [x, fval] = fmincon(@sir_Obj_fun,x0,[],[],[],[],lb,ub)
T = 100;
beta = 0.3;
gamma = 0.1;
[~,p]=ode45(@(t,p) sir_fun(t,p,beta,gamma), [1:1:T],[true_s(1) true_i(1) true_r(1)]);
plot(1:T,p(:,1),'r*',1:T,p(:,2),'b*',1:T,p(:,3),'k*')
没有合适的资源?快使用搜索试试~ 我知道了~
资源详情
资源评论
资源推荐
收起资源包目录
新建文件夹.rar (3个子文件)
新建文件夹
main.m 710B
sir_fun.m 211B
sir_Obj_fun.m 323B
共 3 条
- 1
程籽籽
- 粉丝: 68
- 资源: 4722
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功
评论13