24小时热门版块排行榜    

查看: 1495  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

姜玉旺

木虫 (著名写手)

[求助] 关于matlab微分、及非线性拟合问题

最近一直苦苦思索化学反应速率的问题,想用matlab来求微分拟合非线性曲线求出常数,可我确实不懂matlab,望各位大侠不吝赐教,帮我编写出程序,最好能有图形结合,小弟感激涕零:
我的方程式是(dx)/(dt)=a*c*[(K-1)*x^2-32*K*x+31*K]
t,x的取值为t=0,x=0;t=1,x=0.0355;t=2,x=0.1014;t=3,x=0.1394;t=4,x=0.2147;t=5,x=0.2823;t=6,x=0.3274;t=7,x=0.3717;t=8,x=0.4095.
其中a,K未知,c=1.213.
小弟已经计算好几天了,其在无能为力,拜托了,帮我编出程序。。。感谢了。。。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
姜玉旺: 金币+10, ★★★很有帮助 2013-11-27 18:37:37
用1stOpt直接微分方程拟合:
CODE:
Constant c=1.213;
Variable t,x;
ODEFunction x'=a*c*((K-1)*x^2-32*K*x+31*K);
Data;
0,0
1,0.0355
2,0.1014
3,0.1394
4,0.2147
5,0.2823
6,0.3274
7,0.3717
8,0.4095

均方差(RMSE):0.0144298065185047
残差平方和(SSE):0.00166575452929184
相关系数(R): 0.995283478470763
相关系数之平方(R^2): 0.990589202516862
确定系数(DC): 0.986895420543766
F统计(F-Statistic): 457.855211610462

参数                  最佳估算
--------------------        -------------
a        -0.101555842272903
k        -0.0151477051447373

====== 结果输出 ======

文件: 数据文件-1
No        t        目标 x        计算 x
1        1        0.0355        0.0562847266411599
2        2        0.1014        0.110052998758528
3        3        0.1394        0.162114481657899
4        4        0.2147        0.213185538231459
5        5        0.2823        0.263930448271775
6        6        0.3274        0.314996671293445
7        7        0.3717        0.367048236670994
8        8        0.4095        0.420800978410109
6楼2013-11-22 12:46:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 姜玉旺 的主题更新
信息提示
请填处理意见