x=[3.9,5.3,7.2,9.6,12.9,17.1,23.2,31.4,38.6,50.2,62.9,76.0,92.0,...
106.5,123.2,131.7,150.7,179.3,204.0,226.5,251.4,281.4,309.35]';
n=length(x);
y=x(1:n);%1790 年到 2010 年資料
t=(1:n)';
beta0=[5.3,0.22,400,]; %[x0,r,xm]
[beta,R,J]=nlinfit(t,y,'logisfun',beta0);
%R 為殘差,beta 為待求引數
py=beta(3)./(1+(beta(3)/beta(1)-1)*exp(-beta(2)*t));%預測各年人口
p24=beta(3)./(1+(beta(3)/beta(1)-1)*exp(-beta(2)*24));%預測 2020 年人口
rmse=sqrt(sum(R.^2)/n); %均方誤差根
plot(1:n,y,'*',1:n,py); %作對比圖
%logisfun.m
function yhat=logisfun(beta,x)
yhat=beta(3)./(1+(beta(3)./beta(1)-1).*exp(-beta(2)*x));
完全復制此步驟出現錯誤：The model function 'logisfun'was not found
轉載請註明出處，本文鏈接：https://www.uj5u.com/qita/122818.html
下一篇：HTML與css布局問題