| 查看: 778 | 回复: 8 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 hk1234 的 1 个金币 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
hk1234金虫 (初入文坛)
|
[交流]
【求助】求一个MATLAB 程序
|
||
|
准备做一个燃气锅炉腐蚀模型，采用GM(1,1)灰色模型，前面写着一个MATLAB的程序，由于不是科班出身，在调试过程中老是出错，后面机器出毛病后程序也丢掉了，老板催着要出篇文章，求助各位高手指教，能不能给个已调试好的程序，谢谢！ 小弟新虫，手上金币也不多，尽量发放一点金币给提供帮助的虫友，谢谢！ |
回复此楼» 猜你喜欢
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有52人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有15人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有29人回复
-
关于Filecode分析方法
已经有12人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有11人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有10人回复
-
阿姨
已经有4人回复
-
一作与独作在应聘高校教师时区别大吗
已经有4人回复
-
静等基金结果
已经有21人回复
-
FileCode能看出啥？
已经有28人回复
试一下 这是个改进无偏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
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
