24小时热门版块排行榜

返回列表

当前主题已经存档。

【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 wangzhining156 的 4 个金币

wangzhining156

铁杆木虫 (小有名气)

应助: 1 (幼儿园)
金币: 6433.6
散金: 999
帖子: 217
在线: 50.7小时
虫号: 352764
注册: 2007-04-22
性别: GG
专业: 知识管理

[交流] 求助用matlab实现灰色预测GM(1,1)

大论文收官部分需要用到灰色预测GM(1,1)
我尝试用matlab实现
但是，网上荡的程序都无法调试通过
所以请高手指点迷津

1、能发段可用的程序最好

2、帮俺看看下面的程序问题出在哪里也感激不尽

function [C, P, a, b, Xm0]=graymat(x, m)
n=size(x);
sum0=0;
x1=[];
for i=1:n
  sum0=sum0+x(i);
  x1=[x1 sum0];
end
for i=1:n-1
  b(i)=-(x1(i)+x1(i+1))/2;
  Yn(i)=x(i+1);
end
B=[b' ones(size(b'))];
Coe=pinv(B)*Yn';
a=Coe(1);
b=Coe(2);
for i=1:m+1
  Xm1(i)=(x(1)-b/a)*exp(-a*(i-1))+b/a;
end
for i=1:m
  Xm0(i)=Xm1(i+1)-Xm1(i);
end
Err=Xm0(1:n-1)-x(2:n,1)';
C=std(Err,n-1)/std(x);
RelErr=abs((Err)-mean(Err)*ones(size(Err)));
[Nrow,Ncol]=size(RelErr);
count=0;
for j=1:Ncol
  if RelErr(1,j)<0.6745*std(x)
count=count+1;
  else count=count;
  end
end
P=count/Ncol;

[ Last edited by sunxiao on 2009-3-9 at 09:00 ]

回复此楼