24小时热门版块排行榜

返回列表

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

sujing928

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 98.9
帖子: 65
在线: 124.1小时
虫号: 445718

[交流] 【求助】微分方程拟合求参数

请各位高手帮帮忙，我数据处理中有一个微分方程需要得出参数z的值：
dx/dt=1/(2*96485*3*(1+z/x))

已经有的数据
x：
0.005
0.003218421
0.002230776
0.001452469
0.000973591
0.000605047
0.00038134
0.000300164
0.000303939

t:
0
43200
86400
129600
172800
216000
259200
302400
345600

能否帮我编个代码，文献中看到人家求出的解，卡在那儿了，弄不明白

回复此楼

» 猜你喜欢

自荐读博已经有5人回复
求个博导看看已经有16人回复
上海工程技术大学张培磊教授团队招收博士生已经有4人回复
上海工程技术大学【激光智能制造】课题组招收硕士已经有5人回复
求助院士们，这个如何合成呀已经有4人回复
临港实验室与上科大联培博士招生1名已经有9人回复
需要合成515-64-0，50g，能接单的留言已经有4人回复
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有6人回复
带资进组求博导收留已经有10人回复
最近几年招的学生写论文不引自己组发的文章已经有11人回复

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2010-12-23 21:55:22

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

dingd

铁杆木虫 (职业作家)

应助: 1641 (讲师)
金币: 15037.3
帖子: 3410
在线: 1223.5小时
虫号: 291104

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
xiegangmai(金币+1): 谢谢参与！ 2011-12-04 23:56:21

用1stOpt直接拟合：

CODE:

Variable t,x;

ODEFunction x'=1/(2*96485*3*(1+z/x));

Data;

//t:        x：

0        0.005

43200        0.003218421

86400        0.002230776

129600        0.001452469

172800        0.000973591

216000        0.000605047

259200        0.00038134

302400        0.000300164

345600        0.000303939

均方差(RMSE): 5.66052147217277E-5
残差平方和(SSE): 2.56332026695432E-8
相关系数(R): 0.998700248018427

参数最佳估算
-------------------- -------------
z -0.183470471493407