clcclear
y=[15.19	15.91	16.88	18	16.85	17.12];?%數(shù)據(jù)1行6列的矩陣
n=length（y）;?%?計(jì)算出數(shù)據(jù)y長(zhǎng)度
disp（‘級(jí)比檢驗(yàn)‘）
lambda=y（1:end-1）./y（2:end）;
range=minmax（lambda）
range1=[exp（-2/（n+1））exp（2/（n+1））]
yy=ones（n1）;??%?n行1列的1矩陣
for?i=2:n????%for循環(huán)?累加數(shù)據(jù)
yy（i）=yy（i-1）+y（i）;
end???%原始數(shù)據(jù)累加
B=ones（n-12）;??%?n-1行2列的1矩陣
for?i=1:（n-1）??%for循環(huán)?計(jì)算B矩陣?詳見灰色預(yù)測(cè)模型ppt
B（i1）=-（yy（i）+yy（i+1））/2;
B（i2）=1;
end
BT=B‘;??%?‘表示轉(zhuǎn)置?行列互換
for?j=1:（n-1）???%相當(dāng)于取y矩陣的后面九個(gè)數(shù)，其實(shí)可以寫的更簡(jiǎn)單，YN=y（:2:10）;
YN（j）=y（j+1）;
end
YN=YN‘;?%轉(zhuǎn)置
A=inv（BT*B）*BT*YN;?%ppt中的公式
a=A（1）;??%?上面計(jì)算出來的a是個(gè)矩陣，取A矩陣的第一個(gè)數(shù)
u=A（2）;%?取第二個(gè)
t=u/a;??
t_test=input（‘輸入預(yù)測(cè)個(gè)數(shù)‘）;%預(yù)測(cè)后續(xù)數(shù)據(jù)
i=1:t_test+n;??%?相當(dāng)于i是1到（t_test+n）的數(shù)列，間隔為1
yys（i+1）=（y（1）-t）.*exp（-a.*i）+t;?%數(shù)學(xué)公式詳見ppt??yys為預(yù)測(cè)數(shù)據(jù)
yys（1）=y（1）;?%將y矩陣的第一個(gè)數(shù)字賦給yys矩陣的第一個(gè)位置
for?j=n+t_test:-1:2??%for循環(huán)?將預(yù)測(cè)數(shù)據(jù)yys還原（預(yù)測(cè)出來的是累加數(shù)據(jù)）
ys（j）=yys（j）-yys（j-1）;
end
x=1:n;?%x為1到n的數(shù)列
xs=2:n+t_test;?%?xs為2到n+t_text的數(shù)列
yn=ys（2:n+t_test）;?%將ys矩陣第2到n+t_text的數(shù)賦值給yn
xs=[1xs];