| 查看: 3683 | 回复: 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,(如果可以请划一下图形) |
回复此楼» 猜你喜欢
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
物理学I论文润色/翻译怎么收费?
已经有108人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
桂林理工大学物理学专业招收调剂
已经有18人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有2人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 对两个自变量，一个因变量进行非线性拟合回归，求出目标函数的两个参数。
已经有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人回复
lgycjpcqu
金虫 (正式写手)
- 应助: 35 (小学生)
- 金币: 566.8
- 散金: 411
- 帖子: 371
- 在线: 176.9小时
- 虫号: 1327328
- 注册: 2011-06-20
- 性别: GG
- 专业: 电化学
8楼2013-03-31 10:34:12
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
20