24小时热门版块排行榜    

查看: 1052  |  回复: 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的回帖

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的回帖
查看全部 4 个回答

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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +11 Kittylucky 2026-08-27 12/600 2026-08-28 16:26 by lyxiaomc
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +8 苏东坡二世 2026-08-23 8/400 2026-08-28 14:41 by lkforward
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
请填处理意见