24小时热门版块排行榜    

查看: 1828  |  回复: 4

似然以

新虫 (初入文坛)

[求助] 未知分段点情况下的数据分段拟合 已有1人参与

1stopt进行分段数据拟合求助:
要求拟合为如下形式:
当x<x0时:a+b*x+c*x^2+d*x^3+e*x^4
当 x>x0时:a1+b1*x+c1*x^2+d1*x^3+e1*x^4
如果以上形式不好也可以使用别的形式,但要求必须分段。
自己在拟合的过程中遇到如下问题:
       如果自己设置分段点则会导致在分段点处曲线不光滑
       如果1stopt自动拟合,然后把拟合后的公式导入MATLAB则结果完全不符合
谢谢大家如果有人可以提供 完美满意的 结果答案,可以悬赏自己目前全部金币(11个)进行答谢,谢谢大家。
数据如下:
    X           Y
0.00          0.0276806
0.467692 0.0280567
0.960316 0.0284907
1.45294 0.0289826
1.97149 0.0293588
2.41226 0.0295613
2.41226 0.0295903
2.93081 0.029706
3.47529 0.0297928
3.94199 0.0298507
4.46054 0.0298796
4.92724 0.0297928
5.44579 0.0296481
5.91249 0.0293588
6.45696 0.0286933
6.97552 0.0275937
7.44221 0.0264942
7.93484 0.0253657
8.42746 0.0242662
8.97194 0.0232535
9.49049 0.0223854
9.98312 0.0215463
10.4757 0.020794
10.9684 0.0200127
11.461 0.0193472
11.9536 0.0187396
12.4722 0.0181898
12.9389 0.017669
13.4833 0.0170613
13.9241 0.0166563
14.4686 0.0161354
14.9353 0.0157303
15.4279 0.0152963
15.9724 0.0148912
16.4909 0.0144861
16.9576 0.01411
17.4503 0.0137627
17.9429 0.0133866
18.98 0.0127789
19.4985 0.0124606
19.9912 0.0122002
20.4579 0.0119398
20.9764 0.0116794
21.469 0.011419
21.9876 0.0112164
22.5061 0.010985
23.4914 0.0106088
24.4507 0.0102326
24.9693 0.0100012
25.436 0.00976968
25.9286 0.00968287
26.9398 0.00939352
27.4583 0.00924884
27.9509 0.0090463
28.4436 0.00893056
28.9621 0.00878588
29.9214 0.00852546
30.9585 0.00829398
31.4771 0.00812037
31.9956 0.00800463
32.4623 0.00791782
33.4217 0.00771528
34.4847 0.0074838
35.444 0.00725231
36.4552 0.0070787
37.4664 0.00687616
38.4775 0.00676042
39.4628 0.00661574
40.448 0.00641319
41.4592 0.00632639
42.4704 0.00615278
43.4297 0.00606597
44.4668 0.0059213
45.478 0.00580556
46.4373 0.00568981
47.4744 0.00557407
47.967 0.0055162
48.9782 0.00537153
49.4708 0.00534259
49.9894 0.00528472
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
似然以: 金币+5, ★★★★★最佳答案, 结果非常好 2017-12-18 20:02:25
7.0正版运行效果不错啊,分断点出也可保证连续:
均方差(RMSE): 8.61124930656204E-5
残差平方和(SSE): 5.93228916958122E-7
相关系数(R): 0.999949541890363
相关系数之平方(R^2): 0.999899086326747
修正R平方(Adj. R^2): 0.999896465192377
确定系数(DC): 0.999899086326716
卡方系数(Chi-Square): 2.30579989208035E-5
F统计(F-Statistic): 68368.3971469104
约束条件: a+b*x0+c*x0^2+d*x0^3+e*x0^4-(a1+b1*x0+c1*x0^2+d1*x0^3+e1*x0^4) = 4.16333634234434E-17

参数        最佳估算
----------        -------------
x0        8.17152204425655
a        0.0276200264921194
b        0.00119470890039911
c        -0.000232173030346018
d        3.41452460858937E-5
e        -3.60501562535228E-6
a1        0.0446379145385405
b1        -0.00335111891003501
c1        0.00012475448514077
d1        -2.22732292462842E-6
e1        1.51797444155127E-8
未知分段点情况下的数据分段拟合
c298.jpg

2楼2017-12-17 22:29:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

似然以

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by dingd at 2017-12-17 22:29:41
7.0正版运行效果不错啊,分断点出也可保证连续:
均方差(RMSE): 8.61124930656204E-5
残差平方和(SSE): 5.93228916958122E-7
相关系数(R): 0.999949541890363
相关系数之平方(R^2): 0.999899086326747
修正R平 ...

谢谢谢谢你~请问这个7.0的软件是怎么获得的呢?再次感谢你!
3楼2017-12-18 20:04:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

正版的除了购买没其它途径吧。
4楼2017-12-18 20:46:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

宝恋枫

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by dingd at 2017-12-17 22:29:41
7.0正版运行效果不错啊,分断点出也可保证连续:
均方差(RMSE): 8.61124930656204E-5
残差平方和(SSE): 5.93228916958122E-7
相关系数(R): 0.999949541890363
相关系数之平方(R^2): 0.999899086326747
修正R平 ...

您好,请问您的拟合代码可以提供吗,因为我最近也需要进行这样的拟合,谢谢!
5楼2018-02-06 22:30:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 似然以 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 重要消息,中午系统在维护 +6 yuleib84 2026-08-18 7/350 2026-08-18 15:02 by sytp
[教师之家] 两姐妹,一个辛辛苦苦30岁工作当大学教师,另一个23岁就工作当开车教练 +3 瞬息宇宙 2026-08-11 5/250 2026-08-18 14:55 by LNP@mRNA
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 8/400 2026-08-18 11:47 by gltch
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +4 Tide man 2026-08-14 5/250 2026-08-18 11:35 by Tide man
[基金申请] 时间戳变了,能看出什么问题? +5 基诺咪客 2026-08-17 6/300 2026-08-18 11:06 by hunter无悔
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
信息提示
请填处理意见