版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

flyeagle00

木虫 (正式写手)

小木虫

应助: 2 (幼儿园)
金币: 729.7
散金: 1360
红花: 1
帖子: 417
在线: 1108.1小时
虫号: 694828
注册: 2009-01-31
性别: GG
专业: 力学中的基本问题和方法

[求助] 请教如何回归参数，拟合公式是分段函数？已有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

» 猜你喜欢

想换工作。大多数高校都是评职称时认可5年内在原单位取得的成果吗？已经有7人回复
需要合成515-64-0，50g，能接单的留言已经有4人回复
自荐读博已经有4人回复
上海工程技术大学【激光智能制造】课题组招收硕士已经有4人回复
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有6人回复
带资进组求博导收留已经有10人回复
最近几年招的学生写论文不引自己组发的文章已经有11人回复
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）已经有3人回复
临港实验室与上科大联培博士招生1名已经有8人回复

事繁勿慌，事闲勿荒，有言必信，无欲则刚，和若春风，肃若秋霜，取象于钱，外圆内方

1楼 2015-08-27 22:37:19

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

matlab编程

禁虫 (小有名气)

本帖内容被屏蔽

2楼2015-08-28 00:05:59

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

wurongjun

专家顾问 (职业作家)

专家经验: +831
数学EPI: 9
应助: 791 (博后)
贵宾: 0.308
金币: 24609
散金: 310
红花: 75
帖子: 3004
在线: 881.2小时
虫号: 1368482
注册: 2011-08-14
性别: GG
专业: 计算数学与科学工程计算
管辖: 数学

【答案】应助回帖

★
感谢参与，应助指数 +1
flyeagle00: 金币+1, ★有帮助 2015-08-28 14:24:43

还是用最小二乘法!
就是系数矩阵规模可能较大!

赞一下(1人)

回复此楼

善恶到头终有报,人间正道是沧桑.

3楼2015-08-28 08:48:06

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.5小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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的回帖

flyeagle00

木虫 (正式写手)

小木虫

应助: 2 (幼儿园)
金币: 729.7
散金: 1360
红花: 1
帖子: 417
在线: 1108.1小时
虫号: 694828
注册: 2009-01-31
性别: GG
专业: 力学中的基本问题和方法

引用回帖:

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

能帮忙看一下程序如何改吗？

赞一下

回复此楼

事繁勿慌，事闲勿荒，有言必信，无欲则刚，和若春风，肃若秋霜，取象于钱，外圆内方

5楼2015-08-28 19:09:26

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

flyeagle00

木虫 (正式写手)

小木虫

应助: 2 (幼儿园)
金币: 729.7
散金: 1360
红花: 1
帖子: 417
在线: 1108.1小时
虫号: 694828
注册: 2009-01-31
性别: GG
专业: 力学中的基本问题和方法

引用回帖:

4楼: Originally posted by dingd at 2015-08-28 09:31:03
1stOpt拟合：

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 ...

再请问一下，分段函数分三段以上，怎么表示，中间也是用逗号？

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)/p2)),x>=336,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
2.00 0.00097
4.00 0.001155
8.00 0.001315
12.00 0.0015
24.00 0.0015
60.00 0.00178
72.00 0.00179
84.00 0.001825
96.00 0.001835
108.00 0.001845
120.00 0.0019
132.00 0.001965
144.00 0.001985
156.00 0.002075
167.00 0.00212
168.00 0.00212
168.01 0.002345
168.17 0.002745
168.33 0.002935
168.67 0.003115
169.00 0.003285
170.00 0.003595
172.00 0.003595
176.00 0.003735
180.00 0.003865
192.00 0.003985
204.00 0.004095
216.00 0.004195
228.00 0.00428
240.00 0.00431
252.00 0.00449
264.00 0.004515
276.00 0.00454
288.00 0.004555
300.00 0.004565
312.00 0.004575
336.00 0.00458
336.01 0.004875
336.17 0.00526
336.50 0.00614
337.17 0.006775
338.17 0.00795

赞一下

回复此楼

事繁勿慌，事闲勿荒，有言必信，无欲则刚，和若春风，肃若秋霜，取象于钱，外圆内方

6楼2015-08-29 09:33:25

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.5小时
虫号: 291104
注册: 2006-10-28

引用回帖:

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
.....

赞一下(1人)

回复此楼

7楼2015-08-29 14:02:40

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖