24小时热门版块排行榜    

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

www42884297

铜虫 (著名写手)

[求助] 急求,常微分方程的参数计算问题,有1stopt的帮忙解答一下

(一)实验测得的数据表(物种A3的浓度不计)
t(时间)        c1(物种A1浓度)        c2(物种A2浓度)
15        8.446       
30        7.939        1.458
50        7.687        1.535
80        7.289        1.602
140        6.658        1.717
170        6.531        1.722
200        6.218        1.671
230        5.979        1.620
260        5.591        1.550
290        5.414        1.488
320        4.968        1.433
360        4.692        1.350
400        4.438        1.319
460        4.144        1.294
520        4.041        1.294
580        4.052        1.287
600        4.030        1.292
(二)根据数据建立的动力学模型

-dc1/dt=k1c1+k4c1-k2c2
-dc2/dt=-k1c1+k2c2+k3c2
即:已知c1、c2分别与t的关系,求k1, k2, k3, k4.
回复此楼

» 猜你喜欢

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

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
www42884297: 金币+50, ★★★★★最佳答案 2013-10-02 23:23:17
一楼公式:
-dc1/dt=k1c1+k4c1-k2c2
-dc2/dt=-k1c1+k2c2+k3c2

完全可写成:
-dc1/dt=k4c1-k2c2
-dc2/dt=-k1c1+k3c2

即使这样,似乎存在过拟合现象,参数值不唯一,下面1stOpt结果供参考:

Root of Mean Square Error (RMSE):0.173934237971959
Sum of Square Error:0.665568621055493
Correlation Coef. (R): 0.999292870935962
R-Square: 0.998586241903436
Determination Coef. (DC): 0.998586241903436
F-Statistic: 223.945234354797

Parameter                  Best Estimate
--------------------        -------------
k4        -14.47007732379
k2        -59.6882900889183
k1        621.631020330754
k3        2552.37338260508

====== Output Results ======

File: Data File-1
No        t        Target c1        Calculated c1        Target c2        Calculated c2
1        30        7.939        7.92229452593219        1.458        1.92952644595724
2        50        7.687        7.54395623882783        1.535        1.83738130134023
3        80        7.289        7.18368576432058        1.602        1.74963344345992
4        140        6.658        6.84062043098884        1.717        1.66607915963243
5        170        6.531        6.51393845062598        1.722        1.58651209757037
6        200        6.218        6.20285756048272        1.671        1.51074774553163
7        230        5.979        5.90663258346555        1.62        1.43859877504974
8        260        5.591        5.62455419298996        1.55        1.36989809120694
9        290        5.414        5.3559466805486        1.488        1.30447561574853
10        320        4.968        5.10016685579276        1.433        1.24218004678904
11        360        4.692        4.85660200215829        1.35        1.18285694802669
4楼2013-09-18 23:18:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 www42884297 的主题更新
信息提示
请填处理意见