原创力文档- 1、本文档共8页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
PAGE
PAGE9
（1）灰色预测程序：
y=input(请输入数据?);
?n=length(y);?
yy=ones(n,1);?
yy(1)=y(1);?
for?i=2:n??
yy(i)=yy(i-1)+y(i);?
end??
B=ones(n-1,2);?
for?i=1:(n-1)??????
B(i,1)=-(yy(i)+yy(i+1))/2;?????
B(i,2)=1;?
end?
BT=B;?
for?j=1:n-1?????
YN(j)=y(j+1);
?end
YN=YN;??
A=inv(BT*B)*BT*YN;?
a=A(1);?
u=A(2);
t=u/a;??
t_test=input(请输入需要预测个数：);?
i=1:t_test+n;??
yys(i+1)=(y(1)-t).*exp(-a.*i)+t;?
yys(1)=y(1);?
for?j=n+t_test:-1:2?????
ys(j)=yys(j)-yys(j-1);?
end?
x=1:n;?
xs=2:n+t_test;?
yn=ys(2:n+t_test);?
plot(x,y,^r,xs,yn,*-b);?
det=0;?for?i=2:n??????
det=det+abs(yn(i)-y(i));?
end??
det=det/(n-1);
disp([百分绝对误差为：,num2str(det),%]);?
disp([预测值为：?,num2str(ys(n+1:n+t_test))]);
（2）灰色预测程序：
x0=[255225792600270427192878244025492741276228032816295729123034314230163131325640013393338034703403357433933273341434813355348833963652316233423579357735553658360636453695364837413617362537393726389837853704383539943812387637473625?3573?3756?3456?3669?3527?3425];
?????n=length(x0);?
?????lamda=x0(1:n-1)./x0(2:n)?
?????range=minmax(lamda)?
?????x1=cumsum(x0)?
?????for?i=2:n?
?????????z(i)=0.5*(x1(i)+x1(i-1));?
?????end?
?????B=[-z(2:n),ones(n-1,1)];?
?????Y=x0(2:n);?
?????u=B\Y?
?????x=dsolve(Dx+a*x=b,x(0)=x0);?
?????x=subs(x,{a,b,x0},{u(1),u(2),x1(1)});?
?????yuce1=subs(x,t,[0:n-1]);?
?????digits(6),y=vpa(x)??%为提高预测精度，先计算预测值，再显示微分方程的解?
?????yuce=[x0(1),diff(yuce1)]?
?????epsilon=x0-yuce?????%计算残差?
?????delta=abs(epsilon./x0)?%计算相对误差?
?????rho=1-(1-0.5*u(1))/(1+0.5*u(1))*lamda??%计算级比偏差值?
（3）BP神经网络程序：
p?=[2552?2579?2600;?2579?2600?2704;?2600?2704?2719;?2704?2719?2878;?2719?2878?2440;?2878?2440?2549;?2440?2549?2741;?2549?2741?2762;?2741?2762?2803;?2762?2803?2816;?2803?2816?2957;?2816?2957?2912;?2957?2912?3034;?2912?3034?3142;?3034?3142?3016;?3142?3016?3131;?3016?3131?3256;?3131?3256?4001;?3256?4001?3393;?4001?3393?3380;?3393?3380?3470;?3380?3470?3403;?3470?3403?3574;?3403?3574?3393;?3574?3393?3273;?3393?3273?3414;?3273?3414?3418;?3414?3418?3355;?34
文档评论（0）