作业帮 > 综合 > 作业

怎么用matlab求解Logistic模型中的两个参数Logistic模型的方程变为y=b/(1+((b-a)/a)*e

来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/05/13 06:20:01
怎么用matlab求解Logistic模型中的两个参数Logistic模型的方程变为y=b/(1+((b-a)/a)*exp(-kt)),求b与k.
a=68.65.t=1:10.y=[73.22 80.22 87.69 93.97 99.16 103.38 109.46 114.6 119.85 124.92]
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))          %相关系数 0.9989 挺高的
--------------------------------------
结果:
beta =
  162.4100    0.1498
R2 =
    1.0000    0.9989
    0.9989    1.0000
>>