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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 1stOpt 圆摆线拟合，共享参数模式
151/2返回列表12下一页
查看: 1486  |  回复: 14
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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

» 猜你喜欢
健康快乐
1楼 2016-06-12 16:37:33
已阅   回复此楼   关注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客户端
赞一下(1人) 回复此楼
健康快乐
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
...

请教老师，代码每次运行的结果都一样么？谢谢您啦
赞一下(1人) 回复此楼
健康快乐
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客户端
赞一下(1人) 回复此楼
健康快乐
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的回帖

vjebejw570
祝福
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 的主题更新
151/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳又变了8-15 +10 archvillain 2026-08-15 18/900 2026-08-15 17:21 by 医学老男孩
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 P4JhWrbJGjSe 2026-08-14 5/250 2026-08-15 15:46 by QXtWNz7PI7wZ
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +7 majunge000 2026-08-11 9/450 2026-08-15 13:34 by kiz2011
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-15 11:41 by KxMI1BYBxWX1
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 HFw0lei2R37i 2026-08-14 9/450 2026-08-15 11:33 by KxMI1BYBxWX1
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 04:45 by 4wMiSEwB6436
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[基金申请] 是这周出结果还是下周出结果？ +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-14 22:32 by 4wMiSEwB6436
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥？ +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
信息提示
关闭
请填处理意见
关闭