24小时热门版块排行榜    

查看: 1050  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-28 13:44 by winsaint
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +10 Kittylucky 2026-08-27 11/550 2026-08-28 10:45 by Kittylucky
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +7 苏东坡二世 2026-08-23 7/350 2026-08-27 18:14 by 瞬息宇宙
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 基金不中,共勉 +10 eulota 2026-08-26 10/500 2026-08-27 17:23 by lqllinqiaoli
[文学芳草园] 梦想 +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 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
请填处理意见