| 查看: 318 | 回复: 1 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 wangzhining156 的 4 个金币 | |||
wangzhining156铁杆木虫 (小有名气)
|
[交流]
求助用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 ] |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有13人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有18人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有18人回复
-
2027博士申请-临床检验诊断学、病原生物学、微生物方向
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有10人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有22人回复
-
越来越担心硕士毕业后没有竞争力了
已经有14人回复
-
论文的挂名怎么评判贡献？基金申请中，排名第一的挂名文章
已经有3人回复
-
申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请
已经有3人回复
-
闲聊
已经有6人回复
woshilsh
荣誉版主 (职业作家)
- 应助: 43 (小学生)
- 贵宾: 10.33
- 金币: 21115.1
- 散金: 737
- 红花: 22
- 帖子: 3777
- 在线: 639.7小时
- 虫号: 293038
- 注册: 2006-11-04
- 专业: 计算机应用技术
- 管辖: 仿真模拟
2楼2008-12-02 12:15:28