| 查看: 606 | 回复: 1 | ||
[求助]
微分方程求助1sopt高版本运行程序求参数及相关系数
|
|
第一个程序： Constant c0=6.79,cp=5.79; Variable t,c; ODEFunction c'=-k*(c0-c)^(1/3)*(c-cp) Data; t,c 0 6.79 1 6.39 2 6.25 4 6.1 12 5.9 22 5.8 48 5.79 第二个程序： Constant c0=6.79,cp=5.79; Variable t,c; ODEFunction c'=-k*(c0-c)^(2/3)*(c-cp)^p Data; t,c 0 6.79 1 6.39 2 6.25 4 6.1 12 5.9 22 5.8 48 5.79 第三个程序： Constant c0=6.79,cp=5.79; Variable t,c; ODEFunction c'=-k*(c0-c)^(4/3)*(c-cp)^p Data; t,c 0 6.79 1 6.39 2 6.25 4 6.1 12 5.9 22 5.8 48 5.79 第四个程序： Constant c0=6.79,cp=5.79; Variable t,c; ODEFunction c'=-k*(c-cp)^n Data; t,c 0 6.79 1 6.39 2 6.25 4 6.1 12 5.9 22 5.8 48 5.79 |
回复此楼» 猜你喜欢
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
物理学I论文润色/翻译怎么收费?
已经有249人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 一阶微分方程求解，并优化参数值
已经有14人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】怎么计算二阶的偏微分方程
已经有7人回复
- 【求助】微分方程拟合求参数
已经有7人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】求助一个一阶常微分方程的初值问题
已经有16人回复
- 【求助】急请微分方程高手看下这个方程组！
已经有14人回复
- 【求助】Python 中 求解微分方程 odeint 的用法
已经有11人回复
- 【求助】偏微分方程解法
已经有19人回复
hoopman
金虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 1927.7
- 红花: 1
- 帖子: 91
- 在线: 163.1小时
- 虫号: 584906
- 注册: 2008-08-01
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
yupeifeng: 金币+2, ★有帮助 2012-08-14 08:21:59
csgt0: 金币+2, 多谢应助 2012-08-14 13:13:47
yupeifeng: 金币+2, ★有帮助 2012-08-28 10:22:48
感谢参与，应助指数 +1
yupeifeng: 金币+2, ★有帮助 2012-08-14 08:21:59
csgt0: 金币+2, 多谢应助 2012-08-14 13:13:47
yupeifeng: 金币+2, ★有帮助 2012-08-28 10:22:48
|
模型公式: c' = -k*(6.79-c)^0.3333333333*(c-5.79) c' = -0.425868408288807*(6.79-c)^0.3333333333*(c-5.79) 迭代数: 17 常微分方程算法: 龙格-库塔-费尔博格法(Runge-Kutta-Fehlberg Method) 优化算法: 通用全局优化法(UGO1) 计算结束原因: 达到收敛判定标准 计算用时(时:分:秒:微秒): 00:00:00:327 均方差(RMSE): 6.04260843234664 残差平方和(SSE): 219.0787 相关系数(R): 0 相关系数之平方(R^2): 0 决定系数(DC): -705.971548432206 参数 最佳估算 -------------------- ------------- k 0.425868408288807 ====== 输出结果 ===== 文件：数据文件 - 1 No t 目标c 计算c 1 1 6.39 0 2 2 6.25 0 3 4 6.1 0 4 12 5.9 0 5 22 5.8 0 6 48 5.79 0 第二个程序： 模型公式: c' = -k*(6.79-c)^0.6666666667*(c-5.79)^p c' = -(-0.162244447506964)*(6.79-c)^0.6666666667*(c-5.79)^0.682513076346368 迭代数: 17 常微分方程算法: 龙格-库塔-费尔博格法(Runge-Kutta-Fehlberg Method) 优化算法: 通用全局优化法(UGO1) 计算结束原因: 达到收敛判定标准 计算用时(时:分:秒:微秒): 00:00:00:319 均方差(RMSE): 6.04260843234664 残差平方和(SSE): 219.0787 相关系数(R): 0 相关系数之平方(R^2): 0 决定系数(DC): -705.971548432206 F统计(F-Statistic): 0.00565793631847063 参数 最佳估算 -------------------- ------------- k -0.162244447506964 p 0.682513076346368 ====== 输出结果 ===== 文件：数据文件 - 1 No t 目标c 计算c 1 1 6.39 0 2 2 6.25 0 3 4 6.1 0 4 12 5.9 0 5 22 5.8 0 6 48 5.79 0 第三个程序： 迭代数: 17 常微分方程算法: 龙格-库塔-费尔博格法(Runge-Kutta-Fehlberg Method) 优化算法: 通用全局优化法(UGO1) 计算结束原因: 达到收敛判定标准 计算用时(时:分:秒:微秒): 00:00:00:352 均方差(RMSE): 6.04260843234664 残差平方和(SSE): 219.0787 相关系数(R): 0 相关系数之平方(R^2): 0 决定系数(DC): -705.971548432206 F统计(F-Statistic): 0.00565793631847063 参数 最佳估算 -------------------- ------------- k 0.415244598407298 p 0.344202137552202 ====== 输出结果 ===== 文件：数据文件 - 1 No t 目标c 计算c 1 1 6.39 0 2 2 6.25 0 3 4 6.1 0 4 12 5.9 0 5 22 5.8 0 6 48 5.79 0 第四个程序： 模型公式: c' = -k*(c-5.79)^n c' = -0.554391360086748*(c-5.79)^1.79392199053452 迭代数: 17 常微分方程算法: 龙格-库塔-费尔博格法(Runge-Kutta-Fehlberg Method) 优化算法: 通用全局优化法(UGO1) 计算结束原因: 达到收敛判定标准 计算用时(时:分:秒:微秒): 00:00:03:204 均方差(RMSE): 0.0266573801378389 残差平方和(SSE): 0.0042636954948795 相关系数(R): 0.993570472945939 相关系数之平方(R^2): 0.987182284710018 决定系数(DC): 0.986240965433616 F统计(F-Statistic): 290.718071874963 参数 最佳估算 -------------------- ------------- k 0.554391360086748 n 1.79392199053452 ====== 输出结果 ===== 文件：数据文件 - 1 No t 目标c 计算c 1 1 6.39 6.42165149290096 2 2 6.25 6.24143630782393 3 4 6.1 6.06831214499588 4 12 5.9 5.88880308718633 5 22 5.8 5.84061663939316 6 48 5.79 5.81023213437628 |
2楼2012-08-13 19:47:20
50