24小时热门版块排行榜    

查看: 1516  |  回复: 14
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

1254

金虫 (正式写手)

[求助] 1stOpt 圆摆线拟合,共享参数模式 已有1人参与

尊敬的老师们:
大家好!
我有一个1stOpt 圆摆线的拟合问题,请教大家:
离散的数据点(x, y):  (3.47, 5.48), (3.86, 5.52), (4.18, 5.81), (4.32, 6.28), (4.13, 6.72), (3.82, 6.94), (3.50, 6.97), (3.18, 6.80), (2.98, 6.59), (2.92, 6.31), (3.03, 5.95), (3.36, 5.76), (3.56, 5.75), (3.68, 5.83), (3.74, 5.88)
  自定义拟合参数拟合方程:
x = a*(cos(b)+2/3*cos(3/2*b+w))+c;
y = a*(sin(b)+2/3*sin(3/2*b+w))+d;
其中,x, y 是关于b的参数方程,w, c, d 可通过以上散点得出拟合曲线得出么?
此致
敬礼!
便于您运行代码,附上data:
3.47,5.48
3.86,5.52
4.18,5.81
4.32,6.28
4.13,6.72
3.82,6.94
3.50, 6.97
3.18,6.80
2.98,6.59
2.92,6.31
3.03,5.95
3.36,5.76
3.56,5.75
3.68,5.83
3.74,5.88
回复此楼

» 猜你喜欢

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

1254

金虫 (正式写手)

引用回帖:
11楼: Originally posted by dingd at 2016-06-14 10:40:36
稳定唯一解,效果比b无范围限制要稍差点:

Root of Mean Square Error (RMSE): 0.0731555590881265
Sum of Square Error: 0.160552074764891
Correlation Coef. (R): 0.989302363729821
R-Square: 0.97871916 ...

老师,您好!不好意思,我上次说的约束条件不清楚,请允许我再说一次,这次说全。考虑到问题的实际意义,拟合的参数a, 0~1; b, —pi~pi; w, 0~pi。除此之外,还有一个非常重要的条件,待拟合的参数方程是一个摆线,周期为4pi,但考虑到实际意义,只用—pi~pi之间的半个摆线来拟合。问题有点复杂,期待您的帮助,非常感谢老师!

发自小木虫Android客户端
健康快乐
12楼2016-06-14 11:38:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

1254

金虫 (正式写手)

困扰了好久的问题,盼望虫友解答~~
健康快乐
2楼2016-06-12 16:44:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
1254: 金币+50, ★★★★★最佳答案, 谢谢老师!!! 2016-06-13 17:03:39
1stOpt使用手册里就有类似的例子啊,几乎一样,就是数据和公式略有差别而已。
下面1stOpt计算的结果参考下,也许还有更好的结果:

均方差(RMSE): 0.00826982514393626
残差平方和(SSE): 0.00205170023733841
相关系数(R): 0.999876605059272
相关系数之平方(R^2): 0.999753225344856
修正R平方(Adj. R^2): 0.999685925431167
确定系数(DC): 0.999693595599629
F统计(F-Statistic): -941.585521196443

参数                  最佳估算
--------------------        -------------
a        0.466596275593374
w        -0.0530999992746435
c        3.55158293129851
d        6.21522916474733
3楼2016-06-12 22:07:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1254

金虫 (正式写手)

引用回帖:
3楼: Originally posted by dingd at 2016-06-12 22:07:49
1stOpt使用手册里就有类似的例子啊,几乎一样,就是数据和公式略有差别而已。
下面1stOpt计算的结果参考下,也许还有更好的结果:

均方差(RMSE): 0.00826982514393626
残差平方和(SSE): 0.00205170023733841
...

谢谢老师啦,我在手册里也看到啦,同样的例子,高版本1stopt更加方便,低版本1.5只能用手册的第一种方法,因为b消除不了。。。不知我的理解对不?

发自小木虫Android客户端
健康快乐
4楼2016-06-13 08:56:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:50 by jCd0dEvKHShX
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:40 by jCd0dEvKHShX
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-29 11:39 by jCd0dEvKHShX
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 9/450 2026-08-29 09:59 by jklily
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 09:26 by G6APbkg8SA6w
[基金申请] 有没有仍没收到信息的 +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
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
请填处理意见