| 查看: 3939 | 回复: 8 | |||
[求助]
用matlab里的nlinfit如何拟合分段函数
|
|
Sample Text t=[0 0.0833 0.25 0.5 1 2]; c=[0 0.0833 1.2491 1.7025 0.2623 0.0757]; 函数形式为c=A*(1-exp(-K1*t)) 0<=t<=T c=A*(1-exp(-k1*T))*exp(-k2*(t-T)) t>=T 如何用脚本式或是函数式M文件拟合，求出A，k1,k2,(如果可以请划一下图形) |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有261人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 对两个自变量，一个因变量进行非线性拟合回归，求出目标函数的两个参数。
已经有5人回复
- 原创：用MATLAB做的GUI做数据拟合，通过全局优化算法 （含原文件，code，教程）
已经有55人回复
- matlab非线性拟合计算求助
已经有10人回复
- matlab中关于非线性拟合，用nlinfit.程序总出错。谢谢各位大侠
已经有3人回复
- 用MATLAB或者ORIGIN怎么拟合复数数据
已经有10人回复
- matlab拟合求助
已经有11人回复
- matlab多元非线性拟合
已经有9人回复
- 如何用软件画出这个方程的图形
已经有13人回复
- MATLAB用于三元回归参数
已经有3人回复
- 诚恳求助如何做IC50拟合曲线
已经有7人回复
- 用matlab拟合双因变量实验结果求助
已经有7人回复
- matlab非线性参数拟合问题
已经有7人回复
- 采用matlab中nlinfit拟合，出现问题了，不知道是哪错了，该怎么解决，求高手指点一下
已经有9人回复
- matlab的简单函数拟合
已经有6人回复
- 求助matlab隐函数参数拟合
已经有16人回复
- 【求助】使用MATLAB拟合出按照公式的系数
已经有19人回复
- 【求助】求助matlab数值拟合
已经有8人回复
- 如何用MATLAB拟合指数函数
已经有3人回复
- 【求助】matlab函数拟合，拟合结果并保存，方便调用
已经有4人回复
- 【求助】如何使得matlab拟合的曲线强制经过一个点
已经有8人回复
- 【求助】使用MATLAB怎么实现拟合力场参数的程序？【已解决】
已经有15人回复
- Matlab 的自定义函数拟合
已经有3人回复
- 【其他】自定义的二元函数拟合
已经有6人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★
感谢参与，应助指数 +1
fegg7502: 金币+1, 鼓励交流 2013-04-02 09:24:49
感谢参与，应助指数 +1
fegg7502: 金币+1, 鼓励交流 2013-04-02 09:24:49
|
1stOpt更简单： 均方差(RMSE): 0.0282872879644411 残差平方和(SSE): 0.00480102396229929 相关系数(R): 0.999311354215736 相关系数之平方(R^2): 0.998623182664489 决定系数(DC): 0.998184285936899 卡方系数(Chi-Square): 0.0580621164226486 F统计(F-Statistic): 367.165006987895 参数 最佳估算 ---------- ------------- t1 0.317859946002189 a -0.0366738585259589 k1 -14.2282073650289 k2 3.70145325521581  jv1.jpg |
2楼2013-03-29 08:50:59
lgycjpcqu
金虫 (正式写手)
- 应助: 35 (小学生)
- 金币: 566.8
- 散金: 411
- 帖子: 371
- 在线: 176.9小时
- 虫号: 1327328
- 注册: 2011-06-20
- 性别: GG
- 专业: 电化学
3楼2013-03-29 10:38:11
lgycjpcqu
金虫 (正式写手)
- 应助: 35 (小学生)
- 金币: 566.8
- 散金: 411
- 帖子: 371
- 在线: 176.9小时
- 虫号: 1327328
- 注册: 2011-06-20
- 性别: GG
- 专业: 电化学
【答案】应助回帖
★ ★ ★ ★
cheng1378653: 回帖置顶 2013-03-29 18:09:29
cheng1378653: 金币+2, ★有帮助 2013-03-30 22:02:11
fegg7502: 金币+2, 应助指数+1, 鼓励交流 2013-04-02 09:24:40
cheng1378653: 回帖置顶 2013-03-29 18:09:29
cheng1378653: 金币+2, ★有帮助 2013-03-30 22:02:11
fegg7502: 金币+2, 应助指数+1, 鼓励交流 2013-04-02 09:24:40
|
还是2楼的结果好，用刚才我用1sopt运算的结果作为初始值，用MATLAB进行拟合得到结果和1楼一致，看来非线性最优化还是1sopt强大啊 A = -0.0367 k1 = -14.2282 k2 = 3.7015 T = 0.3179 你需要MATLAB 代码的话如下 clear tic fun1=inline('beta(1).*(1-exp(-beta(2).*t))','beta','t'); fun2=inline('beta(1).*(1-exp(-beta(2).*beta(4)))*exp(-beta(3)*(t-beta(4)))','beta','t'); fun=@(beta,t)((t<0.5).*fun1(beta,t)+(t>=0.5).*fun2(beta,t)); t=[0 0.0833 0.25 0.5 1 2]; c=[0 0.0833 1.2491 1.7025 0.2623 0.0757]; beta=nlinfit(t,c,fun,[-0.05 -6.8 3.7 0.5]); A=beta(1); k1=beta(2); k2=beta(3); T=beta(4); warning off all c1=beta(1).*(1-exp(-beta(2).*t(1:4))); c2=beta(1).*(1-exp(-beta(2).*beta(4))).*exp(-beta(3).*(t(4:6)-beta(4))); cfit(1:4)=c1; cfit(4:6)=c2; plot(t,c,'*') hold on; plot(t,cfit,'r-') toc |
4楼2013-03-29 10:47:34
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
5楼2013-03-29 11:15:58
lgycjpcqu
金虫 (正式写手)
- 应助: 35 (小学生)
- 金币: 566.8
- 散金: 411
- 帖子: 371
- 在线: 176.9小时
- 虫号: 1327328
- 注册: 2011-06-20
- 性别: GG
- 专业: 电化学
6楼2013-03-29 13:53:09
7楼2013-03-30 22:01:05
lgycjpcqu
金虫 (正式写手)
- 应助: 35 (小学生)
- 金币: 566.8
- 散金: 411
- 帖子: 371
- 在线: 176.9小时
- 虫号: 1327328
- 注册: 2011-06-20
- 性别: GG
- 专业: 电化学
8楼2013-03-31 10:34:12
9楼2013-03-31 12:25:54
20