版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

doudounbamvp

银虫 (小有名气)

应助: 12 (小学生)
金币: 307.4
散金: 247
红花: 3
帖子: 130
在线: 45.2小时
虫号: 2636357
注册: 2013-09-05
性别: GG
专业: 污染控制化学

[求助] 1stOpt 非线性拟合问题已有1人参与

第一次用1stOpt对数据进行拟合，但是拟合完成以后没有R和R^2，以及没有计算值。
而且点了几次运算，每次计算出的常数都不一样，这是怎么回事儿，怎样得到针对自定义方程的最优拟合结果？
请给位前辈不吝赐教。

回复此楼

» 猜你喜欢

反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有122人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复

» 本主题相关价值贴推荐，对您同样有帮助:

matlab 非线性拟合已经有16人回复
matlab非线性分段方程拟合急~谢谢~ 已经有7人回复
求助 1stopt低版本方程的拟合（含若干未知参数）已经有10人回复
在用nlinfit对非线性函数拟合时遇到的一些问题请教已经有5人回复
1stopt拟合的动力学参数相关系数不高，有哪些原因已经有3人回复
求大神帮忙拟合一个非线性方程，求出模型参数已经有15人回复
用1stOpt进行复数曲线拟合时，结果出错。急求教！！！已经有5人回复
在进行多元非线性拟合的时候，用matlab拟合的值残差都大于0 已经有14人回复
1stopt中使用sharedmodel拟合曲线出现问题，求指示！！已经有6人回复
用1stopt拟合非线性方程结果与其他软件拟合结果差异大已经有5人回复
大家在使用1stOpt进行数据拟合的时候有没有遇到过这种情况？已经有8人回复
用1stopt或者origin非线性拟合磁导率已经有5人回复
求指点 langmuir线性拟合与非线性拟合不一致已经有12人回复
1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？已经有13人回复

1楼 2015-06-19 19:22:01

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

doudounbamvp

银虫 (小有名气)

应助: 12 (小学生)
金币: 307.4
散金: 247
红花: 3
帖子: 130
在线: 45.2小时
虫号: 2636357
注册: 2013-09-05
性别: GG
专业: 污染控制化学

Title "Type your title here";
Parameters a,b,D;
Variable t,y;
Function y=a+0.608*(b-a)*sum(N=1:9)(1/(N^2)*exp(-N^2*D*t));
Data;
20 0.000274584
40 0.003983045
60 0.005528237
80 0.004498109
100 0.004395096
120 0.007176442
160 0.011451473
180 0.011502979
200 0.014696376
220 0.014541857
240 0.016756632
260 0.017168683
280 0.017580734
300 0.018816888

赞一下

回复此楼

2楼2015-06-19 19:22:25

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

doudounbamvp

银虫 (小有名气)

应助: 12 (小学生)
金币: 307.4
散金: 247
红花: 3
帖子: 130
在线: 45.2小时
虫号: 2636357
注册: 2013-09-05
性别: GG
专业: 污染控制化学

"Type your title here"

迭代数: 2
计算用时(时:分:秒:微秒): 00:00:00:32
优化算法: 麦夸特法(Levenberg-Marquardt)
计算结束原因: 达到收敛判断标准
均方差(RMSE): 0.0121513148247567
残差平方和(SSE): 0.00206716232758492
相关系数(R): 0
相关系数之平方(R^2): 0
决定系数(DC): -3.17832909513088
卡方系数(Chi-Square): 0.148371135
F统计(F-Statistic): 0

参数初始值      最终估算
---------- -------------      --------------
a       0.550798209384084 0.550798209384084
b       0.468684116611257 0.468684116611257
D       0.856468993704766 0.856468993704766

====== 结果输出 =====

No 实测值y 计算值y
1 0.000274584 0.0000000
2 0.003983045 0.0000000
3 0.005528237 0.0000000
4 0.004498109 0.0000000
5 0.004395096 0.0000000
6 0.007176442 0.0000000
7 0.011451473 0.0000000
8 0.011502979 0.0000000
9 0.014696376 0.0000000
10 0.014541857 0.0000000
11 0.016756632 0.0000000
12 0.017168683 0.0000000
13 0.017580734 0.0000000
14 0.018816888 0.0000000

赞一下

回复此楼

3楼2015-06-19 19:22:39

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68123.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

感谢参与，应助指数 +1

引用回帖:

3楼: Originally posted by doudounbamvp at 2015-06-19 19:22:39
"Type your title here"

迭代数: 2
计算用时(时:分:秒:微秒): 00:00:00:32
优化算法: 麦夸特法(Levenberg-Marquardt)
计算结束原因: 达到收敛判断标准
均方差(RMSE): 0.0121513148247567
残差平方 ...

我这边用你的代码，1.5版本的1stopt可以得到稳定的结果：

计算结束原因: 达到收敛判断标准
均方差(RMSE): 0.00109117205004521
残差平方和(SSE): 1.66691901991982E-5
相关系数(R): 0.983009041900267
相关系数之平方(R^2): 0.966306776457681
决定系数(DC): 0.96630677645768
卡方系数(Chi-Square): 0.00182573685408263
F统计(F-Statistic): 344.154702292814

参数最佳估算
---------- -------------
a       0.489584631183171
b       -0.0338156395873558
D       -2.5817072976322E-5

====== 结果输出 =====

No 实测值y 计算值y
1 0.000274584 0.0010553
2 0.003983045 0.0024986
3 0.005528237 0.0039190
4 0.004498109 0.0053172
5 0.004395096 0.0066937
6 0.007176442 0.0080492
7 0.011451473 0.0106990
8 0.011502979 0.0119945
9 0.014696376 0.0132710
10 0.014541857 0.0145291
11 0.016756632 0.0157692
12 0.017168683 0.0169918
13 0.017580734 0.0181973
14 0.018816888 0.0193863

赞一下

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

4楼2015-06-19 19:46:44

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

doudounbamvp

银虫 (小有名气)

应助: 12 (小学生)
金币: 307.4
散金: 247
红花: 3
帖子: 130
在线: 45.2小时
虫号: 2636357
注册: 2013-09-05
性别: GG
专业: 污染控制化学

引用回帖:

4楼: Originally posted by 月只蓝 at 2015-06-19 19:46:44
我这边用你的代码，1.5版本的1stopt可以得到稳定的结果：

计算结束原因: 达到收敛判断标准
均方差(RMSE): 0.00109117205004521
残差平方和(SSE): 1.66691901991982E-5
相关系数(R): 0.983009041900267
相关 ...

如果我设置的参数b实际代表的是初始浓度，也就是说b必须大于零，我该怎么拟合？
我这种拟合不是一次拟合，是一系列的参数拟合，目的是求出参数a，我应该设置什么样的计算方法，以什么方法判断拟合过程相同，怎样判定拟合完成。

赞一下

回复此楼

5楼2015-06-19 21:22:10

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68123.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

引用回帖:

5楼: Originally posted by doudounbamvp at 2015-06-19 21:22:10
如果我设置的参数b实际代表的是初始浓度，也就是说b必须大于零，我该怎么拟合？
我这种拟合不是一次拟合，是一系列的参数拟合，目的是求出参数a，我应该设置什么样的计算方法，以什么方法判断拟合过程相同，怎样判 ...

要求b大于0：
Parameters a,b[0,],D;
Variable t,y;
Function y=a+0.608*(b-a)*sum(N=1:9)(1/(N^2)*exp(-N^2*D*t));
Data;
20       0.000274584
40       0.003983045
60       0.005528237
80       0.004498109
100       0.004395096
120       0.007176442
160       0.011451473
180       0.011502979
200       0.014696376
220       0.014541857
240       0.016756632
260       0.017168683
280       0.017580734
300       0.018816888

赞一下

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

6楼2015-06-19 21:24:46

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖