| 查看: 1930 | 回复: 7 | ||
[求助]
matlab Rung-kutta算法
|
|
由于做论文，做动力学这一块，有动力学模型，有几组浓度随时间变化的曲线数据，来求动力学模型里的参数，进而得到动力学方程。 数据是：(c,t)=（0,0.75）（2,0.7）（3,0.65）（4,0.64）（5,0.62） 动力学模型是：-dc/dt=k(c-c‘)，c'指平衡浓度也就是t=5时的值0.62，采用rung-kutta算法在matlab中求出k进而得到此方程。希望得到高手的指点，不胜感激！我的邮箱是yan.feng21@163.com |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : matlab.docx
2012-06-18 21:28:46, 134.98 K
» 猜你喜欢
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组招收2026硕士研究生1名
已经有0人回复
-
山西大同大学物理学硕士研究生招收调剂生
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有298人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
［调剂信息］211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
一志愿物理所凝聚态物理352求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有2人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
-
070200求调剂，一志愿某211，288分
已经有18人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【分享】《MATLAB 7.0基础教程》【转载】
已经有109人回复
- 【请教】 2D Gaussian Fitting based on Matlab
已经有4人回复
- matlab论文分享（包括神经网络，遗传算法，粒子群算法）
已经有442人回复
- 【分享】《matlab解密》教程PDF版【转载】
已经有705人回复
- 求粒子群优化算法进行变量优选的matlab程序
已经有20人回复
- MATLAB经典算法程序
已经有531人回复
- matlab的函数调用问题
已经有3人回复
- 求回归方程参数估计的最优化算法matlab代码
已经有9人回复
- 数学建模matlab算法大全
已经有720人回复
- 关于MATLAB和OpenCV对相同算法的运算速度对比问题
已经有20人回复
- 【求助】matlab 算法求教
已经有4人回复
- 【求助】多目标pso算法的MATLAB程序
已经有3人回复
- 【求助】matlab 遗传算法工具箱绘图命令是什么？
已经有6人回复
 |
2楼2012-06-18 21:57:24
ymh900423
木虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 4331.5
- 红花: 1
- 帖子: 733
- 在线: 293.2小时
- 虫号: 1359669
- 注册: 2011-08-03
- 性别: GG
- 专业: 有机合成
3楼2012-06-18 22:12:57
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2012-06-19 08:44:43
5楼2012-06-19 09:34:18
【答案】应助回帖
★ ★
感谢参与，应助指数 +1
yupeifeng: 金币+2, ★有帮助 2012-06-19 11:16:23
感谢参与，应助指数 +1
yupeifeng: 金币+2, ★有帮助 2012-06-19 11:16:23
|
Sample TextSample Text 我是这么算的，先把公式化成-1*d(c-c')/dt=k(c-c');d(c-c')/(c-c')=-1*k*dt; ln(c-c')=-k*t;c=exp(-k*t)+c';然后利用matlab非线性最小二乘法进行拟合，方法如下先建立一个名为jack的函数，内容 function c=jack(k,t) c=exp(-k*t)+0.62; 然后在matlab主窗口中输入t=[0,2,3,4,5]; c=[0.75,0.7,0.65,0.64,0.62]; [xx,res]=lsqcurvefit(@jack,[0],t,c) 结果为xx=1.2359 res=0.7571 我把拟合的函数与原数据进行对比发现t=0的点不太符合，其他的符合的非常好，误差虽然是0.7571但是主要是t=0点处的误差。下图是对比图其中红色是拟合出来的图。  |
6楼2012-06-19 10:00:56
7楼2012-06-19 11:50:09
flingbaby
铁杆木虫 (著名写手)
小虫子
- 应助: 4 (幼儿园)
- 金币: 6595.2
- 散金: 10
- 红花: 3
- 帖子: 1209
- 在线: 404.1小时
- 虫号: 904212
- 注册: 2009-11-16
- 性别: GG
- 专业: 聚合物共混与复合材料
8楼2014-09-28 16:37:20
