24小时热门版块排行榜    

查看: 870  |  回复: 3

1135725495

铁杆木虫 (著名写手)

[求助] 求解下面微分方程组的参数 已有2人参与

求解下面微分方程组的参数:
已知方程为:
CB'=-K1*CB^a*(K1/K2*CB^a)^(b/(r-b));
CH'=K2*(K1/K2*CB^a)^(r/(r-b));
其中K1=m*exp(-E1/(8.314*T);
       K2=n*exp(-E2/(8.314*T);
求解m, n, E1,E2,r, a, b.
数据如下:
70℃        CB          CH                              80℃           CB                 CH                       90℃             CB            CH                       100℃  CB            CH
60        11.09        0                         60        10.9237          0                        60        10.86266           0                         60 10.55546          0
90        11.09        0                         90        10.7462           0                        90        10.61313           0                         90        9.87675        0.02843
120        11.09        0.0237                120        10.4357        0.03778                120        10.07194         0.05744                120        8.93854        0.11961
135        11.09        0.0321                135        10.1695        0.0432                135        9.65052         0.12009                135        8.36509        0.30762
150        11.09        0.0876                150        9.7481        0.1203                150        9.07752         0.22735                150        7.68274        0.532
165        11.09        0.1765                165        9.2346        0.21242                165        8.37552         0.38549                165        6.90094        0.81513
180        11.09        0.2954                180        8.6613        0.34579                180        7.56567         0.54642                180        6.07504        1.1321
195        11.09        0.3456                195        8.0058        0.56225                195        6.66582         0.71494                195        5.21599        1.46915
210        11.09        0.5438                210        7.2423        0.83487                210        5.73957         0.93929                210        4.34059        1.8712
225        11.09        0.8765                225        6.4188        1.12793                225        4.78197         1.12055                225        3.51154        2.39291
240        11.09        1.0987                240        5.5353        1.38079                240        3.88077         1.40866                240        2.7776        2.9503
255        11.09        1.423                255        4.5768        1.869                255        3.04407         1.78465                255        2.274        3.36023
270        11.09        1.998                270        4.0146        2.5                        270        2.36982         2.26452                270        1.83049        3.753
285        11.09        2.565                285        3.5703        3.01                        285        1.82637         2.73157                285        1.4453        4.343
300        11.09        3.124                300        3.1158        3.54452                300        1.47042         3.21773                300        1.11904        4.85487
330        11.09        3.923                330        2.4438        4.312                330        0.95172         4.07042                330        0.71974        5.36852
360        11.09        4.367                360        1.9878        4.70402                360        0.68772         4.80723                360        0.4935        5.73301
390        11.09        4.498                390        1.6668        4.8548                390        0.49692         5.39523                390        0.3273        5.9604
420        11.09        4.502                                                                420        0.4328         5.63261                420        0.23855        6.13173 @beefly
回复此楼

» 猜你喜欢

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

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
1135725495: 金币+5, 有帮助 2016-10-12 13:43:33
第一组的CB值都一样为11.09,没问题吗?

微分方程拟合问题自己用1stOpt试试,很容易解决
2楼2016-09-25 11:54:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lu_yu_lan

新虫 (初入文坛)

【答案】应助回帖

感谢参与,应助指数 +1
70℃时,公式中的T是否是 T=70+273.15 ?
可用OpenLu尝试求解,参考:http://muchong.com/bbs/viewthread.php?tid=10232397&fpage=1
3楼2016-09-26 05:32:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1135725495

铁杆木虫 (著名写手)

引用回帖:
2楼: Originally posted by dingd at 2016-09-25 11:54:44
第一组的CB值都一样为11.09,没问题吗?

微分方程拟合问题自己用1stOpt试试,很容易解决

您好!不好意思,第一组数据确实有问题,您能帮忙计算一下吗?我这的1stOpt软件只能计算四个参数。
其中要将温度换算成单位为K,T=273.15+
T=70+273.15
343.15K               
60        11.02346        0
90        10.94583        0
120        10.79057        0.0237
135        10.52441        0.0321
150        10.16953        0.0876
165        9.795797        0.1765
180        9.379922        0.2954
195        8.9829        0.3456
210        8.26205        0.5438
225        7.78518        0.8765
240        7.10869        1.0987
255        6.49874        1.423
270        6.01078        1.998
285        5.51173        2.565
300        5.03486        3.124
330        4.472597        3.923
360        4.033433        4.367
390        3.79278        4.498
4楼2016-09-28 08:39:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 1135725495 的主题更新
信息提示
请填处理意见