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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 1stOpt 圆摆线拟合，共享参数模式
151/2返回列表12下一页
查看: 1489  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 18:22 by 7huP8NPKZCp9
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 HFw0lei2R37i 2026-08-14 6/300 2026-08-15 17:34 by 7huP8NPKZCp9
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +27 医学老男孩 2026-08-13 58/2900 2026-08-15 15:57 by GZUZM
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 7lpolszZVXgi 2026-08-14 3/150 2026-08-15 14:53 by QXtWNz7PI7wZ
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 HFw0lei2R37i 2026-08-14 8/400 2026-08-15 14:41 by QXtWNz7PI7wZ
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 k0dTPqJtl0jt 2026-08-14 3/150 2026-08-15 12:33 by KxMI1BYBxWX1
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 07:09 by 4wMiSEwB6436
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 4/200 2026-08-15 04:17 by 4wMiSEwB6436
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[硕博家园] 读博的好处 +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 要乐观耀哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
关闭
请填处理意见
关闭