| 查看: 1185 | 回复: 6 | |||
iceicy99金虫 (小有名气)
|
[求助]
求大神用1stopt拟合数据 紧急
|
|
各位大神帮我用1stopt拟合下三组数据吧，各种头大中，都不知道往哪个方向解了所以直接给原方程组吧，麻烦各位了 1.A=a*b*C(t)^n 2.dB/dt=c*(b*d*C(t)^n-B) 3.A+B=1200*(e-C(t))-f 4.a+c=1-g 说明：大写的A B C 代表是时间t 的函数 小写的 b n e f g是已知的参数。 a c d是需要求参数 实验数据是与时间t对应的函数值C(t) 自变量是t 应变量是A B C(t)，要求的是参数a c d 数据如下 第一组： b=6476.73 n=0.413 e=14.23 f=11247.4 g=0.6042 t C(t) 0 0 0.5 1.34 1 1.68 2 2.06 3 2.41 5 2.66 7 2.74 9 2.89 11 2.91 第二组 b=17499.19 n=0.294 e=21.93 f=25779.61 g=0.7938 t C(t) 0 0 0.5 0.97 1 1.16 2 1.37 3 1.45 5 1.85 7 2.11 9 2.13 11 2.14 第三组 b=19555.81 n=0.261 e=24.23 f=26467.84 g=0.8573 t C(t) 0 0 0.5 0.73 1 0.89 2 1.16 3 1.35 5 1.68 7 1.87 9 1.95 11 2.00 |
回复此楼» 收录本帖的淘帖专辑推荐
 实验方法技巧 |
» 猜你喜欢
帮我的英语口语老师找学生
已经有4人回复
-
计算三阶力常数用thirdorder扩胞时参数如何选择？临近原子数怎么取值
已经有7人回复
-
物理学I论文润色/翻译怎么收费?
已经有238人回复
-
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有5人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 1stopt能进行隐函数非线性拟合吗？谁能帮拟合一下新型液体的粘度混合计算公式
已经有5人回复
- 求大神用1stOpt帮忙模拟一个复数拟合。
已经有3人回复
- 采用1stopt拟合
已经有8人回复
- 谁用1stopt这个软件啊，离散数据拟合方程，方程未知
已经有5人回复
- 求大神帮忙用1stOpt 拟合复数数据
已经有13人回复
- 求助，matlab/ 1stOpt 自定义微分方程拟合实验数据来求方程中的参数
已经有4人回复
- 求高版本1stopt拟合，
已经有11人回复
- 1stopt中使用sharedmodel拟合曲线出现问题，求指示！！
已经有6人回复
- 用1stopt拟合非线性方程结果与其他软件拟合结果差异大
已经有5人回复
- 求高版本1stopt，解决一个问题，微分方程组参数拟合
已经有3人回复
- 大家在使用1stOpt进行数据拟合的时候有没有遇到过这种情况？
已经有8人回复
- 小妹菜鸟一枚，请教大神帮忙看我的数据可以拟合成什么样的方程，多谢
已经有6人回复
- 有1stOpt4.0版本的朋友帮忙拟合一下复数吧 感激不尽呐
已经有4人回复
- 1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？
已经有13人回复
- 哪位大哥能帮我用1stopt3.0以上版本跑一个拟合程序
已经有11人回复
- 用origin或1stopt拟合一条隐函数方程曲线
已经有4人回复
- 用1stopt进行数据拟合时，无法运行sharedmodel ,求助高手
已经有7人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
iceicy99: 金币+30, ★★★★★最佳答案, 非常谢谢 2013-10-14 10:10:54
csgt0: 金币+2, 谢谢 2013-10-15 13:26:34
感谢参与，应助指数 +1
iceicy99: 金币+30, ★★★★★最佳答案, 非常谢谢 2013-10-14 10:10:54
csgt0: 金币+2, 谢谢 2013-10-15 13:26:34
|
你的问题似乎不是那么简单地用软件可以直接算的，得先整理模型公式： 1：由公式“A=a*b*C(t)^n”和“A+B=1200*(e-C(t))-f”可得： B=1200*(e-C(t))-f-a*b*C(t)^n; (1) 2：将上面B代入“dB/dt=c*(b*d*C(t)^n-B)”可得： dB/dt=c*(b*d*C(t)^n-(1200*(e-C(t))-f-a*b*C(t)^n)) (2) 3：由上面（1）式求dB/dt，得： dB/dt=-1200*C'-a*b*n*C^(n-1)*C' (3) 其中C'=dC/dt 4：（2）（3）式右边相等，得： -1200*C'-a*b*n*C^(n-1)*C' =c*(b*d*C(t)^n-(1200*(e-C(t))-f-a*b*C(t)^n)) 也可写成： C'=(c*(b*d*C^n-(1200*(e-C)-f-a*b*C^n)))/(-1200-a*b*n*C^(n-1)) (4) 5：（4）式可作为最终的微分方程拟合公式，用1stOpt求解，代码如下，只是数据第一列t=0时C也为0，会导致公式（4）计算出错（“C^(n-1)”项），故起始C改为0.1 均方差(RMSE):0.125389408179871 残差平方和(SSE):0.125780029469586 相关系数(R): 0.976606061220058 相关系数之平方(R^2): 0.953759398811755 确定系数(DC): 0.947778509408695 F统计(F-Statistic): 114.895226698086 参数 最佳估算 -------------------- ------------- a -0.086158084535358 d 0.338137688474233 传递参数(PassParameter): ((1-0.6042-a)): 0.481958084535358  c5.jpg [ Last edited by dingd on 2013-10-13 at 22:41 ] |
2楼2013-10-13 22:40:20
iceicy99
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1012.1
- 散金: 48
- 帖子: 69
- 在线: 42.8小时
- 虫号: 1151752
- 注册: 2010-11-20
- 性别: GG
- 专业: 环境污染化学
|
非常感谢，但我先直接复制你的代码，跑一下，显示：没有参数parameters被定义 最后一个方程我对应错了，  第四个方程应该为a+d=1-g （a,d,g为三种组分的比例，和为1） 整理得 C'=(c*(b*d*C^n-(1200*(e-C)-f-a*b*C^n)))/(-1200-a*b*n*C^(n-1)) 增加限制条件 a d g为比例因子，你第五点中C初始值设置为0.1太大，设成0.01 代码如下： parameters a[0,0.3958],c1[0,],d[0,] Constant b=6476.73, n=0.413, e=14.23, f=11247.4, g=0.6042; ConstStr d=1-g-a; Variable t,C; ODEFunction C'=(c1*(b*d*C^n-(1200*(e-C)-f-a*b*C^n)))/(-1200-a*b*n*C^(n-1)); Data; 0 0.01 0.5 1.34 1 1.68 2 2.06 3 2.41 5 2.66 7 2.74 9 2.89 11 2.91 跑不动，跑出结果为空白页，我用的是win7 ，专业破解版1.5，运行时是以管理员身份运行的，帮我看看怎么回事，麻烦了 |
3楼2013-10-14 10:41:37
iceicy99
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1012.1
- 散金: 48
- 帖子: 69
- 在线: 42.8小时
- 虫号: 1151752
- 注册: 2010-11-20
- 性别: GG
- 专业: 环境污染化学
|
非常感谢，但我先直接复制你的代码，跑一下，显示：没有参数parameters被定义 最后一个方程我对应错了， 第四个方程应该为a+d=1-g （a,d,g为三种组分的比例，和为1） 整理得 C'=(c*(b*d*C^n-(1200*(e-C)-f-a*b*C^n)))/(-1200-a*b*n*C^(n-1)) 增加限制条件 a d g为比例因子，你第五点中C初始值设置为0.1太大，设成0.01 代码如下： parameters a[0,0.3958],c1[0,],d[0,] Constant b=6476.73, n=0.413, e=14.23, f=11247.4, g=0.6042; ConstStr d=1-g-a; Variable t,C; ODEFunction C'=(c1*(b*d*C^n-(1200*(e-C)-f-a*b*C^n)))/(-1200-a*b*n*C^(n-1)); Data; 0 0.01 0.5 1.34 1 1.68 2 2.06 3 2.41 5 2.66 7 2.74 9 2.89 11 2.91 跑不动，跑出结果为空白页，我用的是win7 ，专业破解版1.5，运行时是以管理员身份运行的，帮我看看怎么回事，麻烦了 |
4楼2013-10-14 10:43:28
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
|
老版本或破解版无法进行微分方程运算的。 均方差(RMSE):0.584281116229558 残差平方和(SSE):2.73107538225967 相关系数(R): 0.854504288964097 相关系数之平方(R^2): 0.730177579858037 确定系数(DC): -0.133890872662782 F统计(F-Statistic): 5.29151450519206 参数 最佳估算 -------------------- ------------- c1 1.32523654455355 a 0.0329638406282273 传递参数(PassParameter): d: 0.362836159371773  c3.jpg |
» 本帖已获得的红花（最新10朵）
iceicy995楼2013-10-14 11:41:02
iceicy99
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1012.1
- 散金: 48
- 帖子: 69
- 在线: 42.8小时
- 虫号: 1151752
- 注册: 2010-11-20
- 性别: GG
- 专业: 环境污染化学
|
这样啊，那这结果唯一不？因为相关系数太差，不能用啊，再拟合的话相关系数是不是基本上变化不大？ 再麻烦你帮我把下面两组也拟合下吧，我是按你上面的改的，改了常量的值跟实验数据，//PassParameter d1=d+3+6这式子不知什么意思就没改 第二组： Constant b=17499.1, n=0.294,e=21.93,f=25779.61,g=0.7938; ConstStr d=1-g-a; PassParameter d; //PassParameter d1=d+3+6; Variable t,C; ODEFunction C'=(c1*(b*d*C^n-(1200*(e-C)-f-a*b*C^n)))/(-1200-a*b*n*C^(n-1)); Data; 0 0.01 0.5 0.97 1 1.16 2 1.37 3 1.45 5 1.85 7 2.11 9 2.13 11 2.14 第三组 Constant b=19555.81, n=0.261,e=24.23,f=26467.84,g=0.8573; ConstStr d=1-g-a; PassParameter d; //PassParameter d1=d+3+6; Variable t,C; ODEFunction C'=(c1*(b*d*C^n-(1200*(e-C)-f-a*b*C^n)))/(-1200-a*b*n*C^(n-1)); Data; 0 0 0.5 0.73 1 0.89 2 1.16 3 1.35 5 1.68 7 1.87 9 1.95 11 2.00 |
6楼2013-10-14 15:16:12
iceicy99
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1012.1
- 散金: 48
- 帖子: 69
- 在线: 42.8小时
- 虫号: 1151752
- 注册: 2010-11-20
- 性别: GG
- 专业: 环境污染化学
7楼2013-10-14 15:18:15
