| 查看: 672 | 回复: 8 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 hk1234 的 1 个金币 | |||
hk1234金虫 (初入文坛)
|
[交流]
【求助】求一个MATLAB 程序
|
||
|
准备做一个燃气锅炉腐蚀模型，采用GM(1,1)灰色模型，前面写着一个MATLAB的程序，由于不是科班出身，在调试过程中老是出错，后面机器出毛病后程序也丢掉了，老板催着要出篇文章，求助各位高手指教，能不能给个已调试好的程序，谢谢！ 小弟新虫，手上金币也不多，尽量发放一点金币给提供帮助的虫友，谢谢！ |
回复此楼» 猜你喜欢
26调剂 086003
已经有3人回复
-
265求调剂
已经有20人回复
-
271分求调剂学校
已经有11人回复
-
320分人工智能调剂
已经有6人回复
-
求调剂，一志愿南京航空航天大学 ，080500材料科学与工程学硕
已经有8人回复
-
280求调剂
已经有19人回复
-
求材料调剂，一志愿郑州大学289分
已经有15人回复
-
26考研调剂0710 0860
已经有6人回复
-
282求调剂
已经有22人回复
-
283求调剂
已经有5人回复
kuhailangyu
荣誉版主 (著名写手)
小木虫蛀虫防疫站站长
- 仿真EPI: 1
- 应助: 0 (幼儿园)
- 贵宾: 5.006
- 金币: 11615.2
- 红花: 11
- 帖子: 1716
- 在线: 14.7小时
- 虫号: 214037
- 注册: 2006-03-09
- 性别: GG
- 专业: 能源化工
- 管辖: 仿真模拟
2楼2009-08-17 08:31:46
3楼2009-08-17 11:25:48
hk1234
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1017.8
- 散金: 22
- 帖子: 46
- 在线: 18小时
- 虫号: 745621
- 注册: 2009-04-11
- 性别: GG
- 专业: 数学/环境工程
4楼2009-08-18 11:02:56
试一下 这是个改进无偏GM（1.1）
★ ★ ★ ★ ★
hk1234(金币+3,VIP+0):谢谢指教 9-20 09:38
nono2009(金币+2,VIP+0):Thanks for the hand. 9-20 14:17
hk1234(金币+3,VIP+0):谢谢指教 9-20 09:38
nono2009(金币+2,VIP+0):Thanks for the hand. 9-20 14:17
|
clear; clc; x=[];%需要自己输入原始数据 [m,n]=size(x); sum=0; x1=[]; for k=1:n; sum=sum+x(k); x1=[x1 sum];%累加生成 end for k=1:n-1; b(k)=-(x1(k)+x1(k+1))/2; Yn(k)=x(k+1); end B=[b' ones(size(b'))]; Coe=pinv(B)*Yn'; a=Coe(1); u=Coe(2);%求微分方程的系数a和u a1=2-a;%改进的无偏参数 a2=2+a; a3=log(a1/a2);%改进的无偏参数a A=2*u/a2; x0(1)=x(1); for k=2:6 x0(k)=A*exp(a3*(k-1)); end |
5楼2009-09-19 16:32:34
liang_2004
木虫 (正式写手)
始作俑者
- 应助: 16 (小学生)
- 金币: 1803.3
- 散金: 56
- 红花: 5
- 帖子: 377
- 在线: 124.4小时
- 虫号: 848455
- 注册: 2009-09-15
- 专业: 无机非金属类光电信息与功
6楼2009-09-19 19:14:16
★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
nono2009(金币+5,VIP+0):Thanks for sharing your code. 9-21 06:39
hk1234(金币+2,VIP+0):谢谢指教 9-21 09:05
小木虫(金币+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; |
7楼2009-09-21 05:16:31
nono2009(金币+0,VIP+0):欢迎参与讨论。但请不要在此灌水！下次就要扣分了。 9-21 09:56
 |
8楼2009-09-21 09:32:08
nono2009(金币+0,VIP+0):欢迎参与讨论。但请不要在此灌水！下次就要扣分了。 9-21 10:11
  |
9楼2009-09-21 09:53:49
