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（‘預(yù)測值‘‘實際值‘）
xlabel（‘年份‘）;ylabel（‘二氧化碳排放量‘）;
title（‘二氧化碳預(yù)測值和實際值曲線圖（Nair法）‘）
set（gca‘XTick‘[1965:2:2017]）
grid?on
format?short;
forecast=YY（end-4:end）%CO2排放量的預(yù)測結(jié)果
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