24小时热门版块排行榜    

查看: 1857  |  回复: 7

小破孩rkj

新虫 (初入文坛)

[求助] 1stopt参数拟合 已有4人参与

Title "pici;
Parameters a,b,c;
Variable x,y;
Function (1+b)*ln(100*y)-b*ln((a-0.01+y)/a)=c*x;
Data;
0        0.01
2        0.012
4        0.037
6        0.057
8        0.15
10        0.2
12.5        0.25
14.5        0.3
16.5        0.38
18.5        0.45
19.5        0.58
20.5        0.73
21.5        0.75
22.5        0.78
23.5        0.88
24.5        1.03
25.5        1.25
26.5        1.58
27.5        2.3
28.5        3.31
29.5        4.44
30.5        5.2
31.5        6.07
32.5        7.2
33.5        8.25
34.5        8.64
35.5        9.163
36.5        9.4
37.5        9.856
38.5        10.15
39.5        10.6
40.5        11.405
41.5        11.104
42.5        11.654
43.5        11.571
44.5        11.609
45.5        11.541
47.5        11.82
49.5        11.458
51.5        11.578
53.5        11.352
这是我在1stopt上进行的一个参数拟合,在1.5版本可以拟合,但是每次运行参数都不一样,在5.0版本上却不能拟合,提示不行,哪位大神看一下是什么情况?拜托了

1stopt参数拟合
图片1.png
回复此楼

» 猜你喜欢

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

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

小破孩rkj

新虫 (初入文坛)

截图上显示的就是在5.0的版本上运行的时候弹出的
2楼2016-01-16 18:24:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
提示很清楚啊: 试用版求解参数不能超过4个

[ 发自手机版 http://muchong.com/3g ]
3楼2016-01-16 22:29:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lipenggg

专家顾问 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
左边1+b,b,和右边c,是精简的。出现了过拟合现象。所以拟合出来的a,b,c结果不唯一。
一团和气
4楼2016-01-18 12:19:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lipenggg

专家顾问 (正式写手)

引用回帖:
4楼: Originally posted by lipenggg at 2016-01-18 12:19:39
左边1+b,b,和右边c,是精简的。出现了过拟合现象。所以拟合出来的a,b,c结果不唯一。

左边1+b,b,和右边c,不是精简的。出现了过拟合现象。所以拟合出来的a,b,c结果不唯一。
一团和气
5楼2016-01-18 12:26:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wyj210

禁虫 (小有名气)

感谢参与,应助指数 +1
本帖内容被屏蔽

6楼2016-01-18 12:48:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
相关系数(R): 0.981519660589979
相关系数之平方(R^2): 0.963380844124668
决定系数(DC): 0.963201238900815
卡方系数(Chi-Square): 6.15775970768337
F统计(F-Statistic): 1026.01635736153

参数        最佳估算
----------        -------------
a                  -12.0480910721035
b                  8.9494617557344
c                  2.16619084043786
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
7楼2016-01-18 14:12:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

引用回帖:
6楼: Originally posted by wyj210 at 2016-01-18 12:48:41
"pici


迭代数: 3547
计算用时(时:分:秒:微秒): 00:01:04:697
计算结束原因: 达到收敛判断标准
优化算法: 快速简面体爬山法 + 通用全局优化算法(SM2)
目标函数值(最小): 1.05029006736345E-11
均方差 ...

个人觉得,这个结果应该是错误的。
从楼主给出的原始实验数据来看,最后几个点震荡严重,不可能达到确定系数(DC): 0.9999这样的高拟合度。这应该是直接复制楼主的代码运行得出的结果。原拟合代码等式右边的c应当移项到左边,就是说原拟合代码是需要修改的。
@wyj210
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
8楼2016-01-18 14:17:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小破孩rkj 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然评审意见 +12 wangmingqi 2026-08-28 17/850 2026-08-29 08:43 by appleren
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 8/400 2026-08-29 08:32 by foolishmani
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +13 晴天加油 2026-08-26 14/700 2026-08-28 18:36 by huagongfeihu
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 有没有仍没收到信息的 +3 德尚中行 2026-08-27 3/150 2026-08-28 15:24 by rongshuxia
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +8 苏东坡二世 2026-08-23 8/400 2026-08-28 14:41 by lkforward
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
信息提示
请填处理意见