发布网友 发布时间:2022-04-23 05:11
共1个回答
热心网友 时间:2023-09-26 01:33
参考代码: t = 2004:2012;x = [714.33 720.22 730.51 759.67 780.37 798.62 833.89 843.23 857.97];% 模型一: 指数增长模型。y = log(x);a = polyfit(t,y,1);r = a(1);x0 = exp(a(2));x1 = x0 * exp(r*t);% 模型二:阻滞增长模型f = @(a,t) a(1)./(1+(a(1)/x(1)-1)*exp(-a(2)*(t-t(1)))); a = lsqcurvefit(f,[880 1],t,x);plot(t,x,'o',t,x1,'r:.')x2 = f(a,t);plot(t,x,'o',t,x1,'r:.',t,x2,'g*--')legend('原始数据','指数增长模型','阻滞增长模型',2)xlabel 年份; ylabel 人口(万人)