| 查看: 917 | 回复: 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 |
回复此楼» 猜你喜欢
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有32人回复
-
物理学I论文润色/翻译怎么收费?
已经有258人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有78人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
-
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Matlab拟合动力学参数问题，求指导！
已经有9人回复
- 如何数值求解由麦克斯韦方程组谐波展开后得到的椭圆方程
已经有0人回复
- 求如何用有限差分法解matlab中偏微分方程与微分方程联立的方程组
已经有2人回复
- 采用runge-kutta法求解常微分方程组时，时间步长怎么设定，积分时间最长为多少
已经有0人回复
- 微分方程ode45求解，最小二乘法优化微分方程参数，程序运行求助
已经有5人回复
- 不知道matlab程序为什么错了（错误应该存在非线性方程组求解中）
已经有0人回复
- 牛顿到庞卡莱，ODE = NP？——常微分方程发展史和评述 转载
已经有12人回复
- 【资源】推荐一些经典的Matlab书籍、教程和工具箱免费下载【已搜索无重复】
已经有48人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2016-09-25 11:54:44
【答案】应助回帖
感谢参与，应助指数 +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
1135725495
铁杆木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 6868.7
- 散金: 18
- 红花: 1
- 帖子: 1492
- 在线: 209.6小时
- 虫号: 2756771
- 注册: 2013-10-27
- 性别: MM
- 专业: 有机化工
|
您好！不好意思，第一组数据确实有问题，您能帮忙计算一下吗？我这的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
