24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2214  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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

» 猜你喜欢

事繁勿慌,事闲勿荒,有言必信,无欲则刚,和若春风,肃若秋霜,取象于钱,外圆内方
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

引用回帖:
6楼: Originally posted by flyeagle00 at 2015-08-29 09:33:25
再请问一下,分段函数分三段以上,怎么表示,中间也是用逗号?

Function y=if(x<168,0.097/248+0.097/p1*(1-exp(-p1*x/p2)),168<=x<336,0.002+(0.193-0.097)/229+(0.193-0.097)/p1*(1-exp(-p1*(x-168) ...

Function y=if(x<168,0.097/248+0.097/p1*(1-exp(-p1*x/p2)),if((x>=168)and(x<336),0.002+(0.193-0.097)/229+(0.193-0.097)/p1*(1-exp(-p1*(x-168)/p2)),0.004+(0.259-0.193)/205+(0.259-0.193)/p1*(1-exp(-p1*(x-336)/p2))+(0.259-0.226)*(x-336)/20));
Data;
0.00 0.000195
0.17 0.00042
0.33 0.00056
0.67 0.000695
1.00 0.000775
.....
7楼2015-08-29 14:02:40
已阅   回复此楼   关注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
还是用最小二乘法!
就是系数矩阵规模可能较大!
善恶到头终有报,人间正道是沧桑.
3楼2015-08-28 08:48:06
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 283求调剂 +21 那个噜子 2026-04-09 21/1050 2026-04-10 21:18 by zhouxiaoyu
[考研] 293求调剂 +4 勇远库爱314 2026-04-06 4/200 2026-04-10 11:00 by 高维春
[考研] 化工调剂求导师收留!一志愿失利,踏实肯干,有植物提取科研经历 +17 yzyzx 2026-04-09 17/850 2026-04-10 10:29 by wp06
[考研] 调剂 +24 不逢春 2026-04-07 25/1250 2026-04-10 09:39 by 690616278
[考研] 070300化学 求调剂 +13 73372112 2026-04-08 13/650 2026-04-09 20:22 by maddjdld
[考研] 材料工程调剂 +12 小刘同学吖吖 2026-04-06 13/650 2026-04-09 17:07 by luoyongfeng
[考研] 材料考研求调剂总分280 +30 mkjlz1 2026-04-06 35/1750 2026-04-08 21:25 by cyh—315
[考研] 0703化学调剂 348分 +14 唉我超真没招了 2026-04-06 15/750 2026-04-08 19:16 by 我减肥1
[考研] 281求调剂 +10 椰子蘑菇 2026-04-06 10/500 2026-04-08 11:43 by zzucheup
[考研] 326分,一志愿沪9,求生物学调剂 +4 刘墨墨 2026-04-05 4/200 2026-04-08 06:22 by lijunpoly
[考研] 机械调剂 +3 zzzbcb 2026-04-07 3/150 2026-04-07 22:19 by hemengdong
[考研] 338求调剂 +5 小猪红色 678 2026-04-06 6/300 2026-04-07 21:18 by 乔哒哒哒
[考研] 301求调剂 +18 121. 2026-04-04 18/900 2026-04-07 17:49 by 蓝云思雨
[考研] 一志愿西电085401求调剂 +4 sunw1306 2026-04-07 4/200 2026-04-07 16:40 by 啵啵啵0119
[考研] 328求调剂 +4 ghhh88888 2026-04-06 5/250 2026-04-07 14:45 by ghhh88888
[考研] 一志愿华中农业大学0710(A)初试329分 求调剂 +5 一名26考研生 2026-04-04 5/250 2026-04-07 08:54 by 18828373951
[考研] 材料调剂 +5 小刘同学吖吖 2026-04-06 5/250 2026-04-06 18:34 by sherry_1901
[考研] 307求调剂 +3 所念及所望 2026-04-06 3/150 2026-04-06 17:30 by 土木硕士招生
[考研] 296求调剂 +3 汪!?! 2026-04-05 4/200 2026-04-05 20:13 by 啵啵啵0119
[考研] 考研调剂 +5 四川王涛 2026-04-04 5/250 2026-04-04 22:18 by 啵啵啵0119
信息提示
请填处理意见