matlab 中logistic模型求解这是我写的源程序,直接从网上下载的建立m文件function y=zout(beta,t) % 保存为 zout.m 文件a=68.65;b=beta(1);k=beta(2);y=b./(1+((b-a)/a)*exp(-k*t));----------------------------------------------------
来源:学生作业帮助网 编辑:作业帮 时间:2024/06/21 06:51:26
![matlab 中logistic模型求解这是我写的源程序,直接从网上下载的建立m文件function y=zout(beta,t) % 保存为 zout.m 文件a=68.65;b=beta(1);k=beta(2);y=b./(1+((b-a)/a)*exp(-k*t));----------------------------------------------------](/uploads/image/z/5922822-30-2.jpg?t=matlab+%E4%B8%ADlogistic%E6%A8%A1%E5%9E%8B%E6%B1%82%E8%A7%A3%E8%BF%99%E6%98%AF%E6%88%91%E5%86%99%E7%9A%84%E6%BA%90%E7%A8%8B%E5%BA%8F%2C%E7%9B%B4%E6%8E%A5%E4%BB%8E%E7%BD%91%E4%B8%8A%E4%B8%8B%E8%BD%BD%E7%9A%84%E5%BB%BA%E7%AB%8Bm%E6%96%87%E4%BB%B6function+y%3Dzout%28beta%2Ct%29+%25+%E4%BF%9D%E5%AD%98%E4%B8%BA+zout.m+%E6%96%87%E4%BB%B6a%3D68.65%3Bb%3Dbeta%281%29%3Bk%3Dbeta%282%29%3By%3Db.%2F%281%2B%28%28b-a%29%2Fa%29%2Aexp%28-k%2At%29%29%3B----------------------------------------------------)
matlab 中logistic模型求解这是我写的源程序,直接从网上下载的建立m文件function y=zout(beta,t) % 保存为 zout.m 文件a=68.65;b=beta(1);k=beta(2);y=b./(1+((b-a)/a)*exp(-k*t));----------------------------------------------------
matlab 中logistic模型求解
这是我写的源程序,直接从网上下载的
建立m文件
function y=zout(beta,t) % 保存为 zout.m 文件
a=68.65;
b=beta(1);k=beta(2);
y=b./(1+((b-a)/a)*exp(-k*t));
----------------------------------------------------------
运行以下:
clear;clc
t=1:10;
y=[73.22 80.22 87.69 93.97 99.16 103.38 109.46 114.6 119.85 124.92];
beta0=[200 0.1];
beta=nlinfit(t,y,'zout',beta0)
plot(t,y,'*',t,zout(beta,t))
R2=corrcoef(y,zout(beta,t))
--------------------------------------
结果:
beta =
162.4100 0.1498
R2 =
1.0000 0.9989
0.9989 1.0000
beta0=[200 0.里面的数值200,0.1是怎么确定的?
还有如果结果显示
371.01 +1.1315e-005i -0.1308 -5.6501e-008i
matlab 中logistic模型求解这是我写的源程序,直接从网上下载的建立m文件function y=zout(beta,t) % 保存为 zout.m 文件a=68.65;b=beta(1);k=beta(2);y=b./(1+((b-a)/a)*exp(-k*t));----------------------------------------------------
beta0=[200 0.1];是b=beta(1);k=beta(2);的初始值.初始值确定是一个比较头疼的事,完全凭自己的感觉来初选,一般可用随机数来初步确定,看结果是否你要数据,然后再调整,直到满意为止.
如果结果显示
371.01 +1.1315e-005i -0.1308 -5.6501e-008i
这说明b=beta(1);k=beta(2);的初始值选择有问题.