24小时热门版块排行榜

返回列表

当前主题已经存档。

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

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

hk1234

金虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1017.8
散金: 22
帖子: 46
在线: 18小时
虫号: 745621
注册: 2009-04-11
性别: GG
专业: 数学/环境工程

[交流] 【求助】求一个MATLAB 程序

准备做一个燃气锅炉腐蚀模型，采用GM(1,1)灰色模型，前面写着一个MATLAB的程序，由于不是科班出身，在调试过程中老是出错，后面机器出毛病后程序也丢掉了，老板催着要出篇文章，求助各位高手指教，能不能给个已调试好的程序，谢谢！
小弟新虫，手上金币也不多，尽量发放一点金币给提供帮助的虫友，谢谢！

回复此楼

» 猜你喜欢

假如你的研究生提出不合理要求已经有12人回复
实验室接单子已经有7人回复
全日制（定向）博士已经有5人回复
萌生出自己或许不适合搞科研的想法，现在跑or等等看？已经有4人回复
Materials Today Chemistry审稿周期已经有4人回复
参与限项已经有3人回复
对氯苯硼酸纯化已经有3人回复
所感已经有4人回复
要不要辞职读博？已经有7人回复
北核录用已经有3人回复

高山之止，云在风端！

1楼 2009-08-17 00:31:45

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

hyde1022

铜虫 (小有名气)

应助: 12 (小学生)
金币: 494.4
散金: 18
红花: 6
帖子: 56
在线: 19.3小时
虫号: 587875
注册: 2008-08-22
专业: 理论和计算化学

★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
nono2009(金币+5,VIP+0):Thanks for sharing your code. 9-21 06:39
hk1234(金币+2,VIP+0):谢谢指教 9-21 09:05

这是我前两年编过的一个灰色预测函数，和大家分享一下吧

%alpha为均值生成数列的权重,一般取0.5即可
%gmax为生成预测数列的元素个数
function p=hsyc(x0,gmax,alpha)
x1=zeros(size(x0));
z1=zeros(size(x0));
x1yc=zeros(size(x0));
xyc=zeros(size(x0));
cancha=zeros(size(x0));
n=length(x0);
for i=1:n %做一次累加（AGO）生成数列x1
   if i==1
      x1(i)=x0(i);
   else
      x1(i)=x1(i-1)+x0(i);
   end
end
for i=2:n %生成均值数列z1
   z1(i)=alpha*x1(i)+(1-alpha)*x1(i-1);
end
c=sum(z1);
d=sum(x0)-x0(1);
z1x0=z1.*x0;
e=sum(z1x0);
z12=z1.^2;
f=sum(z12);
a=(c*d-(n-1)*e)/((n-1)*f-c^2);
b=(d*f-c*e)/((n-1)*f-c^2);
x1yc(1)=x1(1);
for k=1

gmax-1)
   x1yc(k+1)=(x0(1)-b/a)*exp(-a*k)+b/a;
end
for i=1:gmax
   if i==1
      xyc(1)=x1yc(1);
   else
      xyc(i)=x1yc(i)-x1yc(i-1);
   end
end
for i=1:n
   cancha(i)=(x0(i)-xyc(i))/x0(i);
end
a
b
cancha %如果残差r<0.2，则认为达到一般要求;如果r<0.1，则认为达到较高要求
x0
p=xyc;