| 查看: 766 | 回复: 8 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 hk1234 的 1 个金币 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
hk1234金虫 (初入文坛)
|
[交流]
【求助】求一个MATLAB 程序
|
||
|
准备做一个燃气锅炉腐蚀模型，采用GM(1,1)灰色模型，前面写着一个MATLAB的程序，由于不是科班出身，在调试过程中老是出错，后面机器出毛病后程序也丢掉了，老板催着要出篇文章，求助各位高手指教，能不能给个已调试好的程序，谢谢！ 小弟新虫，手上金币也不多，尽量发放一点金币给提供帮助的虫友，谢谢！ |
回复此楼» 猜你喜欢
奇怪，两个人的filecode固定段从头到尾一模一样
已经有3人回复
-
国自然结果
已经有5人回复
-
8月时间戳变的，举个手。玩一下，释放压力
已经有12人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有8人回复
-
据悉今年马上要出结果了
已经有8人回复
-
关于代码变化问题，想知道的进来
已经有23人回复
-
FileCode能看出啥？
已经有6人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有9人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有20人回复
-
2026国自然放榜时间
已经有9人回复
hk1234
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1017.8
- 散金: 22
- 帖子: 46
- 在线: 18小时
- 虫号: 745621
- 注册: 2009-04-11
- 性别: GG
- 专业: 数学/环境工程
4楼2009-08-18 11:02:56
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
试一下 这是个改进无偏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
20