Logistic曲线的三种参数估计方法及matlab算法实现

主要是关于logistic的三种求参数的算法：Yule算法，Rhodes算法，Nair算法matlab算法的实现。

clear;clc;
%Nair算法
%author:朱伟杰
%date:2018-1-24
X=[480.9522468.8469.5573.8737.8869.8933.7977.2...
    997.71120.31176.11284.81422.11462.11499.7...
    1473.11539.2163717711886.51994.62145.72292...
    2396.823872484.42580.82750.22915.73163.83231.9...
    3319.53319.63484.3550.63613.93833.14471.25283...
    5803.26415.56797.97033.57636.38209.88979.1]
n=length（X）-1
for t=1:n
    Z（t）=1/X（t）-1/X（t+1）
    S（t）=1/X（t）+1/X（t+1）
end
X1=[ones（461） S（1:n）‘]
Y=Z‘
[BBintrrintstats]=regress（YX1）%最小二乘（OLS）
gamma=B（11）
beta=B（21）
b=log（（1-beta）/（1+beta））
c=gamma*（1+exp（b））/（2*（exp（b）-1））
a=exp（（sum（log（1./X（1:n）-c））-n*（n+1）*b/2）/n）
XX=1965:2016
YY=1./（c+a*exp（b*（[XX-1965]）））
plot（XXYY‘r-o‘）
hold on
plot（XX（1:length（X））X‘g-^‘）
legend（‘预测值‘‘实际值‘）
xlabel（‘年份‘）;ylabel（‘二氧化碳排放量‘）;
title（‘二氧化碳预测值和实际值曲线图（Nair法）‘）
set（gca‘XTick‘[1965:2:2017]）
grid on
format short;
forecast=YY（end-4:end）%CO2排放量的预测结果
MAPE=sum（abs（YY（1:n+1）-X）./X）/length（X）%平均相对差值
abc

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-05-29 19:33  logistic模型及其matlab算法\logistic模型及其matlab算法\
     目录           0  2018-01-25 00:13  logistic模型及其matlab算法\logistic模型及其matlab算法\logistic的matlab代码\
     文件        1095  2018-01-25 00:04  logistic模型及其matlab算法\logistic模型及其matlab算法\logistic的matlab代码\Nair.m
     文件        1062  2018-01-25 00:03  logistic模型及其matlab算法\logistic模型及其matlab算法\logistic的matlab代码\Rhodes.m
     文件        1047  2018-01-25 00:03  logistic模型及其matlab算法\logistic模型及其matlab算法\logistic的matlab代码\Yule.m
     文件      149399  2018-05-29 19:33  logistic模型及其matlab算法\logistic模型及其matlab算法\logistic模型及其matlab算法.docx