| 查看: 2294 | 回复: 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
» 猜你喜欢
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有179人回复
-
课题组招1-2名硕士
已经有2人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
日内瓦国际发明展申报完整流程（科研人员通用指南）
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【分享】《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人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2012-06-19 08:44:43
ymh900423
木虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 4331.5
- 红花: 1
- 帖子: 733
- 在线: 293.2小时
- 虫号: 1359669
- 注册: 2011-08-03
- 性别: GG
- 专业: 有机合成
3楼2012-06-18 22:12:57
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
