查看: 18423|回复: 2

[问题求助] 急需求助灰色模型

字体大小: 正常放大

1 主题	5 听众	75 积分

升级 73.68%

TA的每日心情

	开心 2015-9-15 14:48

签到天数: 17 天

[LV.4]偶尔看看III

电梯直达

1^#

发表于 2012-5-6 21:41 |只看该作者 |正序浏览

|招呼Ta 关注Ta

我要预测10年之后的人口，用灰色模型
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));                c(i)=-0.5*(x1(i)+x1(i+1))
B=[c' ones(n-1,1)];
au=inv(B'*B)*B'*y;             %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度

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)=abs(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);                   %计算后验差比 *2
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
x0=[31.26
32.09
33.39
35.45
40.52
43.52
47.86
51.45
55.6
60.14
64.82
68.65
73.22
80.22
87.69
93.97
99.16
103.38
109.46
114.6
119.85
124.92
132.04
139.45
150.93
165.13
181.93
196.83
212.38
228.07
241.45
251.03
]
我在窗口输入：gm(x0) 为啥会出现
Error in ==> gm at 22
error(i)=abs(yc(i)-x0(i));                   %计算残差值
那位高手帮忙谢谢

zan

计算, function, 渐开线, error