原创力文档- 1、本文档共3页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
GM(1,1)法度模范
GM( 1, 1) 的matlab程序
gm.m:
function f=gm(x0,m) %定义为函数gm(x)
n=length(x0);
x1=zeros(1,n);
x1(1)=x0(1);
for i=2:n %计算累加序列x1
x1(i)=x1(i-1)+x0(i);
end
i=2:n; %对原始数列平行移位并负值给y
y(i-1)=x0(i);
y=y;
i=1:n-1;
c(i)=-0.5*(x1(i)+x1(i+1));
B=[c ones(n-1,1)];
au=inv(B*B)*B*y; %计算参数a,u矩阵
i=1:m; %计算预测累加数列的值
ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1);
yc(1)=ago(1);
i=1:m-1; %还原数列的值
yc(i+1)=ago(i+1)-ago(i);
i=2:n;
error(i)=yc(i)-x0(i); %计算残差值
yc(1)=ago(1);
i=1:m-1; %修正还原数列的值
yc(i+1)=ago(i+1)-ago(i);
c=std(error)/std(x0); %计算后验差比
relerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);
p=0;
for i=2:ncol
if relerror(1,i)0.6745*std(x0)
p=p+1;
end
end
p=p/(n-1);
w1=min(abs(error));
w2=max(abs(error));
i=1:n; %计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
w=sum(w)/(n-1);
au %输出参数a,u的值
ago; %输出累加数列ago的值
x0; %输出原始序列值
f=yc; %输出预测的值
error; %输出残差的值
c; %输出后验差比的值
p; %输出小误差概率的值
w %输出关联度
end
文档评论（0）