| 查看: 2500 | 回复: 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
» 猜你喜欢
时间戳又变了
已经有19人回复
-
建议基金发布提前给出明确的时间点
已经有4人回复
-
filecode，4个jtjc了
已经有14人回复
-
看来今天不会放榜了？
已经有9人回复
-
让我中一个面上吧！
已经有13人回复
-
人气不行了
已经有7人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有18人回复
-
感觉是下周放榜了
已经有12人回复
-
我面上完蛋了
已经有8人回复
-
只有每年这种时候来逛逛小木虫
已经有25人回复
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
|
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
|
本帖内容被屏蔽 |
2楼2015-08-28 00:05:59
wurongjun
专家顾问 (职业作家)
-
专家经验: +831 - 数学EPI: 9
- 应助: 791 (博后)
- 贵宾: 0.308
- 金币: 24609
- 散金: 310
- 红花: 75
- 帖子: 3004
- 在线: 881.4小时
- 虫号: 1368482
- 注册: 2011-08-14
- 性别: GG
- 专业: 计算数学与科学工程计算
- 管辖: 数学
3楼2015-08-28 08:48:06
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
flyeagle00: 金币+50, ★★★很有帮助, 谢谢，请问软件能分享一下吗 2015-08-28 14:24:24
感谢参与，应助指数 +1
flyeagle00: 金币+50, ★★★很有帮助, 谢谢，请问软件能分享一下吗 2015-08-28 14:24:24
|
1stOpt拟合： 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
40