| 查看: 2192 | 回复: 7 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
matlab非线性分段方程拟合 急~谢谢~ 已有2人参与
|
|||
初学者求帮忙！ 如图所示是一组非线性分段方程，其中待拟合参数有：T0，Ta，tm，δT，以及分段点ts。 我尝试使用lsqcurvefit进行拟合，但是代码编写有一些问题。。查阅相关资料也不太知道这种分段点就是未知参数的方程应该如何编程，希望大家给予帮助 xdata = [ 5.5000 6.5000 7.5000 8.5000 9.5000 10.5000 11.5000 12.5000 13.5000 14.5000 15.5000 16.5000 17.5000 18.5000 19.5000 20.5000 21.5000 22.5000 23.5000 24.5000 1.5000 2.5000 3.5000 4.5000]; ydata = [299.0545 300.8932 304.5193 308.6715 311.6907 315.0223 315.7924 317.1786 315.0809 313.1319 312.8464 309.4712 305.4651 304.1893 302.5152 301.9055 303.4350 299.9959 299.2613 299.7677 299.4591 298.5113 298.7621 298.6650]; x0 = [299.05 40 14 0.5 17.4]; function F = fun(x,xdata) if xdata > x(5) w = 4/3*(x(3) - 5.4); F = a + b*cos( 3.14/w.*( xdata -c ) ); elseif xdata <= x(5) w = 4/3*(c - 5.4); k = w/3.14 * ( atan( 3.14/w*(x(5)-x(3)) ) - d/b*asin( 3.14/w*( x(5)-x(3) ) ) ); F = a + d + ( b*cos( 3.14/w*(ee-c) ) - d ) * k/( k + xdata - ee ); end end [x,resnorm] = lsqcurvefit(@fun,x0,xdata,ydata); 这样运行之后会报错： Error using feval Output argument "F" (and maybe others) not assigned during call to "F:\DTC\pros\fun1012.m>fun1012". Error in lsqcurvefit (line 203) initVals.F = feval(funfcn_x_xdata{3},xCurrent,XDATA,varargin{:}); Caused by: Failure in initial user-supplied objective function evaluation. LSQCURVEFIT cannot continue. 最后应该就是没有F的返回值了。。。我这个语句写的应该是肯定有问题但是我不太知道是什么问题。。希望能够请高人指点一下 谢谢！！！！！ |
回复此楼» 猜你喜欢
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
物理学I论文润色/翻译怎么收费?
已经有210人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
桂林理工大学物理学专业招收调剂
已经有17人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
此贴已删除
已经有7人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何用Matlab拟合大量散点的边界曲线
已经有5人回复
- 在用nlinfit对非线性函数拟合时遇到的一些问题请教
已经有5人回复
- Matlab 二元线性函数拟合 求助
已经有3人回复
- 不可线性化的非线性曲线方程拟合
已经有5人回复
- 求助：以下非线性拟合方程可以用1stopt求解么？用什么方法
已经有16人回复
- 求助用matlab 多元非线性拟合
已经有4人回复
- 求Origin或Matlab高手非线性拟合公式
已经有5人回复
- matlab 拟合求助
已经有8人回复
- 关于matlab拟合
已经有4人回复
- matlab非线性拟合求助
已经有7人回复
- 用matlab里的nlinfit如何拟合分段函数
已经有8人回复
- 关于matlab微分、及非线性拟合问题
已经有5人回复
- matlab数据拟合 希望高手帮忙写一下matlab的程序 跪谢
已经有12人回复
- matlab数据拟合问题求助
已经有7人回复
- 关于拟合：同一个方程，非线性拟合与转换后线性拟合
已经有8人回复
- 等温吸附 线性 非线性拟合结果不一致，希望大家多多指教
已经有15人回复
- 1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？
已经有13人回复
- 用matlab拟合双因变量实验结果求助
已经有7人回复
- matlab非线性参数拟合问题
已经有7人回复
- 采用matlab中nlinfit拟合，出现问题了，不知道是哪错了，该怎么解决，求高手指点一下
已经有9人回复
6楼2014-10-15 10:09:39
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68131.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2014-10-13 09:56:31
3楼2014-10-13 12:02:24
戴钢盔的猪头
木虫 (知名作家)
- 应助: 200 (高中生)
- 金币: 5083.9
- 散金: 1310
- 红花: 18
- 帖子: 6001
- 在线: 319.9小时
- 虫号: 2854040
- 注册: 2013-12-07
- 性别: GG
- 专业: 波谱分析与成像分析
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
bloomi: 金币+10, ★★★★★最佳答案 2014-10-15 10:09:56
感谢参与，应助指数 +1
bloomi: 金币+10, ★★★★★最佳答案 2014-10-15 10:09:56
|
分段函数最小二乘拟合目标函数导数不连续，可采用两段函数的加权加和作为拟合模型。加权因子假设分别为a和b，可设计连续光滑函数作为加权因子。这个函数需要表现出阶越性，可采用某个三角函数。三角函数都是初等函数，连续可导，故可设计加权函数满足：当t小于ts时，a=1，b=0，否则a=0，b=1。具体哪个三角函数呢？楼主自己去研究下吧，授人以鱼不如授人以渔。这个方法分两段，分多段都可以用。 [ 发自手机版 http://muchong.com/3g ] |
4楼2014-10-14 15:48:59
50
