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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 工科数学 » 请教如何回归参数，拟合公式是分段函数？
51/1返回列表
查看: 2206  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

flyeagle00

木虫 (正式写手)

小木虫

[求助] 请教如何回归参数，拟合公式是分段函数？ 已有2人参与

问题：已知曲线f~t, 和拟合公式，拟合公式是分段函数， 求拟合公式里待定参数？
比如matlab里inline或lsqcurvefit函数可以实现曲线回归，但是对于拟合公式是分段函数的形式该如何编写程序？谢谢解答
不是分段函数可以如下来实现。
clear all
mu=inline('0.097/248+0.097/beta(1)*(1-exp(-beta(1)/beta(2)*t))','beta','t');
t=[0        0.17        0.33        0.67        1        2        4        8        12        24        36        48        60        72        84        96        108        120        132        144        156        168        168        168.17        168.33        168.67        169        170        172        176        180        192        204        216        228        240        252        264        276        288        300        312        324        336
];
y=[0.000195        0.00042        0.00056        0.000695        0.000775        0.00097        0.001155        0.001315        0.00145        0.0015        0.001655        0.00169        0.00178        0.00179        0.001825        0.001835        0.001845        0.0019        0.001965        0.001985        0.002075        0.00212        0.002345        0.002745        0.002935        0.003115        0.003285        0.003445        0.003595        0.003735        0.003865        0.003985        0.004095        0.004195        0.00428        0.00431        0.00449        0.004515        0.00454        0.004555        0.004565        0.004575        0.00458        0.00458
];
beta0=[ 100  1000 ];
[beta,r]=nlinfit(t,y,mu,beta0  );
beta
plot(t,y,'rp');
y1=0.097/248+0.097/beta(1)*(1-exp(-beta(1)/beta(2)*t));
hold on
plot(t,y1,'b+-')
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : 问题.doc
    • 2015-08-27 22:27:51, 51 K

» 猜你喜欢
事繁勿慌，事闲勿荒，有言必信，无欲则刚，和若春风，肃若秋霜，取象于钱，外圆内方
1楼 2015-08-27 22:37:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
flyeagle00: 金币+50, ★★★很有帮助, 谢谢，请问软件能分享一下吗 2015-08-28 14:24:24
1stOpt拟合：
CODE:
Function y=if(x<168,0.097/248+0.097/p1*(1-exp(-p1*x/p2)),0.193/229+0.193/p1*(1-exp(-p1*x/p2)));
Data;
0        0.000195
0.17        0.00042
0.33        0.00056
0.67        0.000695
1        0.000775
2        0.00097
4        0.001155
8        0.001315
12        0.00145
。。。。

Root of Mean Square Error (RMSE):0.00053117761646034
Sum of Square Error:1.24145850500535E-5
Correlation Coef. (R): 0.952603849622092
R-Square: 0.907454094314828
Determination Coef. (DC): 0.854457197552128
Chi-Square: 0.00413271275566929
F-Statistic: 288.574902321555

Parameter        Best Estimate
----------        -------------
p1        43.4005751213438
p2        7756.71840625649
请教如何回归参数，拟合公式是分段函数？
c1.jpg
赞一下 回复此楼
4楼2015-08-28 09:31:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

matlab编程

禁虫 (小有名气)
本帖内容被屏蔽
2楼2015-08-28 00:05:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wurongjun

专家顾问 (职业作家)

【答案】应助回帖

★
感谢参与，应助指数 +1
flyeagle00: 金币+1, ★有帮助 2015-08-28 14:24:43
还是用最小二乘法!
就是系数矩阵规模可能较大!
赞一下(1人) 回复此楼
善恶到头终有报,人间正道是沧桑.
3楼2015-08-28 08:48:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

flyeagle00

木虫 (正式写手)

小木虫
引用回帖:
3楼: Originally posted by wurongjun at 2015-08-28 08:48:06
还是用最小二乘法!
就是系数矩阵规模可能较大!

能帮忙看一下程序如何改吗？
赞一下 回复此楼
事繁勿慌，事闲勿荒，有言必信，无欲则刚，和若春风，肃若秋霜，取象于钱，外圆内方
5楼2015-08-28 19:09:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学工程与技术324调剂 +23 孙常华 2026-04-09 25/1250 2026-04-11 00:07 by 骑牛渡寒江
[考研] 材料085601调剂 +25 何润采123 2026-04-10 27/1350 2026-04-10 23:17 by Ftglcn90
[考研] 中科院总分315求调剂 +8 lallalh 2026-04-09 8/400 2026-04-10 19:30 by dick_runner
[考研] 284求调剂 +9 让我上岸吧阿西 2026-04-09 11/550 2026-04-10 19:18 by 靖jing
[考研] 362求调剂 +10 我要考大 2026-04-06 14/700 2026-04-10 17:00 by luoyongfeng
[考研] 材料调剂 +11 一样YWY 2026-04-05 11/550 2026-04-10 09:32 by 钟洲2011
[考研] 一志愿双非085400电子信息344 求调剂，对材料和化学方向也感兴趣 +8 无情的小羊 2026-04-09 9/450 2026-04-10 09:30 by 松花缸1201
[考研] 293求调剂 +5 勇远库爱314 2026-04-08 5/250 2026-04-10 08:46 by vgtyfty
[考研] 085400电子信息类（川大控制工程）求调剂可跨专业 求老师联系 +3 626776879 2026-04-08 3/150 2026-04-09 16:05 by 猪会飞
[考研] 262求调剂 +10 天下第一文 2026-04-04 13/650 2026-04-09 15:16 by 探123
[考研] 一志愿西南大学生物学学硕344 求生物学相关调剂/生物与医药 +7 超人不会飞@ 2026-04-08 7/350 2026-04-09 09:35 by gong120082
[考研] 一志愿华南师范大学0702物理学305调剂 +4 念常安 2026-04-07 6/300 2026-04-08 22:53 by bljnqdcc
[考研] 机械专硕273请求调剂 +6 庚申壬申 2026-04-07 6/300 2026-04-08 22:41 by bljnqdcc
[考研] 22408 318分求调剂 +4 勤奋的小笼包 2026-04-06 6/300 2026-04-07 15:05 by 纸鹤555
[论文投稿] Decision: Revise for Editor还会送审吗 100+3 CccccccccFD 2026-04-04 5/250 2026-04-07 10:58 by 北京莱茵润色
[考研] 一志愿安徽某211 0703化学总分339求调剂 +7 晚风不晚 2026-04-04 7/350 2026-04-06 14:06 by houyaoxu
[考研] 322求调剂 +3 嗯哼哼恒 2026-04-05 3/150 2026-04-05 19:52 by nepu_uu
[考研] 298求调剂 +7 manman511 2026-04-05 7/350 2026-04-05 10:29 by 唐沐儿
[考研] 313求调剂 +3 海日海日 2026-04-04 3/150 2026-04-05 07:48 by 544594351
[考研] 可跨专业调剂 +3 周的得地 2026-04-04 6/300 2026-04-04 22:21 by barlinike
信息提示
关闭
请填处理意见
关闭