function?grey（）
clcclear
[y]=textread（‘input.txt‘‘%f‘）;
%y=[15.19	15.91	16.88	18	16.85	17.12];?%數(shù)據(jù)1行6列的矩陣
n=length（y）;?%?計算出數(shù)據(jù)y長度
disp（‘級比檢驗‘）
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)?計算B矩陣?詳見灰色預(yù)測模型ppt
B（i1）=-（yy（i）+yy（i+1））/2;
B（i2）=1;
end
BT=B‘;??%?‘表示轉(zhuǎn)置?行列互換
for?j=1:（n-1）???%相當于取y矩陣的后面九個數(shù)，其實可以寫的更簡單，YN=y（:2:10）;
YN（j）=y（j+1）;
end
YN=YN‘;?%轉(zhuǎn)置
A=inv（BT*B）*BT*YN;?%ppt中的公式
a=A（1）;??%?上面計算出來的a是個矩陣，取A矩陣的第一個數(shù)
u=A（2）;%?取第二個
t=u/a;??
t_test=50;
%t_test=input（‘輸入預(yù)測個數(shù)‘）;%預(yù)測后續(xù)數(shù)據(jù)
i=1:t_test+n;??%?相當于i是1到（t_test+n）的數(shù)列，間隔為1
yys（i+1）=（y（1）-t）.*exp（-a.*i）+t;?%數(shù)學(xué)公式詳見ppt??yys為預(yù)測數(shù)據(jù)
yys（1）=y（1）;?%將y矩陣的第一個數(shù)字賦給yys矩陣的第一個位置
for?j=n+t_test:-1:2??%for循環(huán)?將預(yù)測數(shù)據(jù)yys還原（預(yù)測出來的是累加數(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];
yn=[y（1）yn];
fid=fopen（‘output.txt‘‘wt‘）;
fprintf（fid‘%g\n‘yn）;
fclose（fid）;
plot（xy‘^r‘xsyn‘.-b‘）;?
det=0;%畫線
DD=[0];
for?i=2:n?
??DD=[DD;abs（yn（i）-y（i））];
??det=det+abs（yn（i）-y（i））;??
??end
det=det/（n-1）;??
det1=0;%畫線
EE=[];
for?i=2:n?
EE=[EE;abs（yn（i）-y（i））/yn（i）];
det1=det1+abs（yn（i）-y（i））/yn（i）;??
end
det1=det1/（n-1）;?
disp（[‘絕對誤差為：‘num2str（det）]）;
disp（[‘百分絕對誤差為：‘num2str（det1）‘%‘]）;
disp（[‘預(yù)測值為：‘num2str（ys（n+1:n+t_test））]）;
p=0.5;
max_err=max（abs（DD））;
r=0;
for?k=1:n
r=r+p*max_err/（abs（DD（k））+p*max_err）;
end
r=r/n;?%r表示關(guān)聯(lián)度
aver_e0=sum（DD）/n;%預(yù)測誤差均值
aver_x0=sum（y）/n;%原始數(shù)據(jù)均值
s1=sqrt（sum（（y-aver_x0）.^2）/n）;%原始數(shù)據(jù)標準差
s2=sqrt（sum（（DD-aver_e0）.^2）/n）;%預(yù)測數(shù)據(jù)標準差
c=s2/s1;?%方差比
s0=0.6745*s1;
m=0;
for?k=1:n
if?abs（DD（k）-aver_e0）m=m+1;
end
end
P=m/n;?%小誤差概率
if?P>0.95&c<0.35
????disp（‘預(yù)測精度好‘）;
else?if?P>0.8&c<0.5
????????disp（‘預(yù)測合格‘）;
????else?if?P>0.7&c<0.65
????????????disp（‘預(yù)測勉強合格‘）;
????????else
????????????disp（‘預(yù)測不合格‘）;
????????end
????end
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????1621327??2018-05-17?22:49??人口預(yù)測模型\2007nian.rar
?????文件????????2336??2018-05-19?14:47??人口預(yù)測模型\grey.m
?????文件????????3522??2018-05-22?21:15??人口預(yù)測模型\leslie.m
?????文件??????847069??2018-05-17?22:53??人口預(yù)測模型\二胎政策對中國人口的影響.pdf
?????文件?????1239552??2018-05-18?20:39??人口預(yù)測模型\人口預(yù)測模型.doc
?????文件????56824569??2018-05-12?22:50??人口預(yù)測模型\人口預(yù)測點評.rar
?????文件?????1338147??2018-05-18?23:02??人口預(yù)測模型\第五屆數(shù)學(xué)建模競賽論文.doc
?????目錄???????????0??2018-06-04?15:50??人口預(yù)測模型\