24小时热门版块排行榜    

查看: 1311  |  回复: 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的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

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的回帖

1254

金虫 (正式写手)

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

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

请教老师,代码每次运行的结果都一样么?谢谢您啦
健康快乐
6楼2016-06-13 17:29:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1254

金虫 (正式写手)

引用回帖:
9楼: Originally posted by dingd at 2016-06-13 22:42:22
为什么还要用Matlab来画图验证?1stOpt自身就足够了,1stOpt是随机赋初值的,你的这道题非线性程度高,求解难度大,10次随机运行有8次能得到相同的最优解已经是很好了。你自己试试Matlab,随机赋初值,你计算100次也 ...

老师,您好,我还不太会用1stopt,希望多多向您学习。最近摆线拟合困扰我很久了,您给出的结果非常好。

但是,我忘了说前提条件了,拟合的参数方程,b的范围要求为-pi,pi,您能再算下结果么?非常感谢老师!

发自小木虫Android客户端
健康快乐
10楼2016-06-13 23:18:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

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的回帖
祝福
5楼2016-06-13 12:42:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

引用回帖:
6楼: Originally posted by 1254 at 2016-06-13 17:29:41
请教老师,代码每次运行的结果都一样么?谢谢您啦...

10次有8次一样吧。
7楼2016-06-13 18:12:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1254

金虫 (正式写手)

引用回帖:
7楼: Originally posted by dingd at 2016-06-13 18:12:37
10次有8次一样吧。...

总有不一样的,只能自己验算下,通过matlab画图,比对下数据和拟合的曲线。然后选一个拟合结果。您是这么做的么?

发自小木虫Android客户端
健康快乐
8楼2016-06-13 22:11:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

为什么还要用Matlab来画图验证?1stOpt自身就足够了,1stOpt是随机赋初值的,你的这道题非线性程度高,求解难度大,10次随机运行有8次能得到相同的最优解已经是很好了。你自己试试Matlab,随机赋初值,你计算100次也得不到一次1stOpt的结果,而且每次结果都会不同。

[ 发自手机版 http://muchong.com/3g ]
9楼2016-06-13 22:42:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 1254 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 335分 | 材料与化工专硕 | GPA 4.07 | 有科研经历 +4 cccchenso 2026-03-23 4/200 2026-03-23 23:00 by 徐ckkk
[考研] 一志愿武理材料工程348求调剂 +6  ̄^ ̄゜汗 2026-03-19 9/450 2026-03-23 19:53 by pswait
[考研] 接收2026硕士调剂(学硕+专硕) +4 allen-yin 2026-03-23 6/300 2026-03-23 15:04 by 汪!?!
[考研] 306求调剂 +9 chuanzhu川烛 2026-03-18 9/450 2026-03-23 13:17 by luoyongfeng
[考研] 307求调剂 +3 余意卿 2026-03-21 3/150 2026-03-23 10:32 by Iveryant
[考研] 求调剂材料学硕080500,总分289分 5+3 @taotao 2026-03-19 21/1050 2026-03-23 10:17 by 冠c哥
[考研] 070300,一志愿北航320求调剂 +3 Jerry0216 2026-03-22 5/250 2026-03-23 09:16 by 。。堂堂
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +11 枫桥ZL 2026-03-18 13/650 2026-03-22 20:26 by edmund7
[考研] 311求调剂 +6 冬十三 2026-03-18 6/300 2026-03-22 20:18 by edmund7
[考研] 307求调剂 +11 冷笙123 2026-03-17 11/550 2026-03-22 20:16 by edmund7
[考研] 269专硕求调剂 +6 金恩贝 2026-03-21 6/300 2026-03-22 14:31 by ColorlessPI
[考研] 286求调剂 +10 Faune 2026-03-21 10/500 2026-03-21 23:34 by 314126402
[考研] 考研调剂 +3 呼呼?~+123456 2026-03-21 3/150 2026-03-21 20:04 by 无际的草原
[考研] 294求调剂材料与化工专硕 +15 陌の森林 2026-03-18 15/750 2026-03-20 23:28 by JourneyLucky
[考研] 一志愿中海洋材料工程专硕330分求调剂 +8 小材化本科 2026-03-18 8/400 2026-03-20 23:16 by JourneyLucky
[考研] 一志愿苏州大学材料求调剂,总分315(英一) +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[考研] 一志愿南理工085701环境302求调剂院校 +3 葵梓卫队 2026-03-20 3/150 2026-03-20 19:28 by zhukairuo
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
信息提示
请填处理意见