| 查看: 1370 | 回复: 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 |
回复此楼» 收录本帖的淘帖专辑推荐
 实验方法技巧 |
» 猜你喜欢
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
物理学I论文润色/翻译怎么收费?
已经有155人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有1人回复
-
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生
已经有5人回复
-
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究
已经有12人回复
-
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有1人回复
-
如何从铁电相到顺电相。
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
