| 查看: 1293 | 回复: 14 | |||
heisekaishui银虫 (小有名气)
|
[求助]
关于1STOPT拟合的小问题 已有2人参与
|
|
各位前辈，我最近用1stopt拟合了一组数据，但是发现拟合的数据使用不同阶的龙格库塔方程结果不同，但是精度都差不多，不知道是什么原因，希望大家不吝赐教。 代码如下： Parameters k1,k2,k3; Variable t,a,b,c; ODEFunction a'=-k1*a-k2*a; b'=k1*a; c'=k2*a; Data; t=0,2,4,6,8,10; a=1.161575096,1.098850041,1.022186084,0.952491579,0.885120223,0.825879893; b=0,0.060153328,0.132977117,0.197793007,0.254614938,0.306825415; c=0,0.002571727,0.006411895,0.01129051,0.021839935,0.028869787; 同时，我在设定参数的范围时，并没有起到效果，而且1stopt的结果纵坐标范围是0-4而不是0-10不知道是什么原因 |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有242人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 利用Excel数据前处理，结合1stopt数据自定义拟合初值，美美origin数据作图分析！
已经有72人回复
- 关于用1stOpt进行参数拟合问题
已经有4人回复
- 求助！求助！小女子快崩溃了，关于1STOPT的拟合问题，谢谢各位大神啦，求解答啊！
已经有4人回复
- 个人发现软件1stopt的隐函数拟合功能存在严重问题
已经有8人回复
- 求助一个多元函数拟合的问题。matlab或1stopt
已经有2人回复
- [求助] 1stopt 微分方程代码求助
已经有5人回复
- 求助大神，关于1stopt带积分的函数拟合
已经有7人回复
- 1stopt 非线性约束拟合
已经有6人回复
- 关于用debye公式拟合热熔数据 求解debye温度的问题
已经有6人回复
- 关于1stopt拟合隐函数方程时如何设置目标函数的问题
已经有2人回复
- 关于1stOpt使用的问题
已经有1人回复
- 1stopt软件进行非线性拟合的问题
已经有1人回复
- 学1stopt要多久？
已经有12人回复
- 求助关于1stopt拟合~··~~
已经有4人回复
ybkooo
至尊木虫 (著名写手)
^^
- 应助: 32 (小学生)
- 金币: 11056.5
- 散金: 178
- 红花: 15
- 帖子: 2312
- 在线: 545.4小时
- 虫号: 2710519
- 注册: 2013-10-09
- 性别: GG
- 专业: 控制理论与方法
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
heisekaishui: 金币+10, ★★★很有帮助, 非常感谢，我还想问一下如果我要设置参数的取值范围应该怎么设置，比如我想让所有的 2015-07-18 13:32:29
heisekaishui: 金币+10, ★★★很有帮助, 非常感谢，我还想问一下如果我要设置参数的取值范围应该怎么设置，比如我想让所有的 2015-07-18 13:32:48
感谢参与，应助指数 +1
heisekaishui: 金币+10, ★★★很有帮助, 非常感谢，我还想问一下如果我要设置参数的取值范围应该怎么设置，比如我想让所有的 2015-07-18 13:32:29
heisekaishui: 金币+10, ★★★很有帮助, 非常感谢，我还想问一下如果我要设置参数的取值范围应该怎么设置，比如我想让所有的 2015-07-18 13:32:48
|
看看这样能不能解决你的问题 Parameters k1,k2,k3; Variable t,a,b,c; SharedModel; ODEFunction a'=-k1*a-k2*a; b'=k1*a; c'=k2*a; RowData; 0,2,4,6,8,10; 1.161575096,1.098850041,1.022186084,0.952491579,0.885120223,0.825879893; 0,0.060153328,0.132977117,0.197793007,0.254614938,0.306825415; 0,0.002571727,0.006411895,0.01129051,0.021839935,0.028869787; |
2楼2015-07-17 21:15:36
ybkooo
至尊木虫 (著名写手)
^^
- 应助: 32 (小学生)
- 金币: 11056.5
- 散金: 178
- 红花: 15
- 帖子: 2312
- 在线: 545.4小时
- 虫号: 2710519
- 注册: 2013-10-09
- 性别: GG
- 专业: 控制理论与方法
|
Parameters k1,k2,k3,k4; 修改方程，这个可以得到稳定结果 Variable t,a,b,c; ODEFunction a'=-k1*a-k2*a; b'=k3*a; c'=k4*a; RowData; 0,2,4,6,8,10; 1.161575096,1.098850041,1.022186084,0.952491579,0.885120223,0.825879893; 0,0.060153328,0.132977117,0.197793007,0.254614938,0.306825415; 0,0.002571727,0.006411895,0.01129051,0.021839935,0.028869787; |
3楼2015-07-17 21:41:22
ybkooo
至尊木虫 (著名写手)
^^
- 应助: 32 (小学生)
- 金币: 11056.5
- 散金: 178
- 红花: 15
- 帖子: 2312
- 在线: 545.4小时
- 虫号: 2710519
- 注册: 2013-10-09
- 性别: GG
- 专业: 控制理论与方法
|
坐标显示0-4可能是他的作图功能不好，或不是正版软件，只看他的计算值就行了（计算值不是有，每时刻的计算值吗） 2 1.098850041 .060153328 .002571727 1.08633992640023 6.97307002743706E-02 5.60446897778746E-03 For a: 4 1.022186084 .132977117 .006411895 1.01597773545576 .134944946131689 1.08459367441951E-02 DC: 0.99487234139559 6 .952491579 .197793007 .01129051 .950172900633618 .195935268051574 1.57479148657442E-02 R: 0.999959508427279 8 .885120223 .254614938 .021839935 .888630242171105 .25297524940825 2.03323920723293E-02 RMSE: 0.00692414160103415 10 .825879893 .306825415 .028869787 .831073699086233 .306320753592402 2.46199328847424E-02 SSE: 0.000239718684555859 |
4楼2015-07-17 21:48:42
ybkooo
至尊木虫 (著名写手)
^^
- 应助: 32 (小学生)
- 金币: 11056.5
- 散金: 178
- 红花: 15
- 帖子: 2312
- 在线: 545.4小时
- 虫号: 2710519
- 注册: 2013-10-09
- 性别: GG
- 专业: 控制理论与方法
5楼2015-07-17 21:51:27
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
感谢参与，应助指数 +1
|
稳定唯一： 均方差(RMSE):0.00523400858904149 残差平方和(SSE):0.000410922688652401 相关系数(R): 0.99229618847388 相关系数之平方(R^2): 0.98465172565979 确定系数(DC): 0.949184328781669 F统计(F-Statistic): 585.146344863956 参数 最佳估算 -------------------- ------------- k1 0.0310189175554483 k2 0.0024794201096689 |
» 本帖已获得的红花（最新10朵）
heisekaishui6楼2015-07-17 23:23:35
heisekaishui
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 215
- 散金: 30
- 帖子: 52
- 在线: 59.2小时
- 虫号: 2445260
- 注册: 2013-05-02
- 性别: GG
- 专业: 催化化学
7楼2015-07-18 13:33:32
heisekaishui
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 215
- 散金: 30
- 帖子: 52
- 在线: 59.2小时
- 虫号: 2445260
- 注册: 2013-05-02
- 性别: GG
- 专业: 催化化学
8楼2015-07-18 13:37:22
heisekaishui
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 215
- 散金: 30
- 帖子: 52
- 在线: 59.2小时
- 虫号: 2445260
- 注册: 2013-05-02
- 性别: GG
- 专业: 催化化学
9楼2015-07-18 13:39:07
heisekaishui
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 215
- 散金: 30
- 帖子: 52
- 在线: 59.2小时
- 虫号: 2445260
- 注册: 2013-05-02
- 性别: GG
- 专业: 催化化学
10楼2015-07-18 13:41:26
5
投票: