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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 1stopt参数拟合
81/1返回列表
查看: 1856  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

小破孩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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2016-01-16 18:23:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

月只蓝

主管区长 (职业作家)

【答案】应助回帖

引用回帖:
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

新虫 (初入文坛)
截图上显示的就是在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的回帖
相关版块跳转 我要订阅楼主 小破孩rkj 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 07:59 by gy1nBQXYQJqL
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 07:55 by popo0625
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 07:46 by gy1nBQXYQJqL
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +16 爱看书的可乐 2026-08-26 18/900 2026-08-29 07:32 by winsaint
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +4 gouxfjh 2026-08-28 6/300 2026-08-29 07:15 by gouxfjh
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 基金未中，这种答复是模板吗？ +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[职场人生] 学生撞见辅导员送外卖，第二天全班都沉默了 +3 鲸鱼融金_浙江_� 2026-08-22 3/150 2026-08-26 09:01 by zzuzxg
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+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
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
关闭
请填处理意见
关闭