clc
clear all
x=[2002 2003 2004 2005 2006 2007 2008 2009 2010]';
y=[1573.52 1813.38 2129.33 2462.22 2644.45 2812.12 2938.74 3037 3217.06]';
X=[ones(length(x),1),x];                %构造自变量观察值矩阵
[b,bint,r,rint,stats]=regress(y,X);     %线性回归建模与评价
b,stats
%作图
figure;
plot(x,y,'-+');                          %数据散点图
legend('第二产业能源消费量');
figure(2)
rcoplot(r,rint);                        %作残差分析图
hold on;
figure(3)
x0=[x' 2011 2012 2013 2014 2015]';
z=b(1)+b(2)*x0;
plot(x,y,'*',x0,z,'-r^');
legend('第二产业能源消费量预测图','第二产业能源消费量回归线图');             %预测及作回归线图
%检验
SSE=sum((y-(b(1)+b(2)*x)).^2)           %残差平方和
STD=sqrt(SSE/(length(x)-2))             %标准误差
DELTA=2*STD                             %0.05显著水平下的边际误差
ci=[z-DELTA,z+DELTA]                    %0.95置信区间