end yhat2
plot(1:15,y,'*')
a15=2*yhat1(end)-yhat2(end)
b15=2*(yhat1(end)-yhat2(end))/(n-1) y2005=a15-5*b15 y2006=a15-4*b15 y2007=a15-3*b15 y2008=a15-2*b15 y2009=a15-b15 y2010=a15
y2011=a15+b15 y2012=a15+2*b15 y2013=a15+3*b15 y2014=a15+4*b15 y2015=a15+5*b15
（5）关联分析程序：
clc,clear
x=[1309 1614 1620 2094 2537 2900 3270 3391.... 4089 5058 6000 6944 8100 9399.7 10815
91 92 110 124 136 142 140 140 145 147 186 190 200 215 243; 65 70 212 263 414 421 666 1040 1425 1559.... 1761 2206 2316 2433 3088; 106.6 99.6 98.8 96.8 98.5 98.4...
100.2 100.1 103 100.9 101.2 104 106.1 99.1 102.1; 4838.9 5160.3 5425.1 5854.02 6280 6859.6...
7702.8 8472.2 9421.6 10493 11759.5 13785.8 15780.76 17174.65 19109; 48.8 46.6 44.7 42.1 39.4 38.2 37.7 37.1 37.7 36.7 35.8 36.3 37.9 36.5 35.7; ];
for i=1:5
x(i,:)=x(i,:)/x(i,1); %标准化数据 end for i=6
x(i,:)=x(i,1)./x(i,:); %标准化数据 end data=x;
n=size(data,2); %求矩阵的列数，即观测时刻的个数 ck=data(1,:); %提出参考数列 bj=data(2:end,:); %提出比较数列 m2=size(bj,1); %求比较数列的个数 for j=1:m2
t(j,:)=bj(j,:)-ck;
25
end
mn=min(min(abs(t'))); %求最小差 mx=max(max(abs(t'))); %求最大差 rho=0.5; %分辨系数设置
ksi=(mn+rho*mx)./(abs(t)+rho*mx); %求关联系数 r=sum(ksi')/n %求关联度
[rs,rind]=sort(r,'descend') %对关联度进行排序
26