| 查看: 1085 | 回复: 10 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】求助:最小二乘法曲线拟合参数已有4人参与
|
|||
|
请高手帮忙编个程序，求出参数k，非常感谢！！！ 方程：y=y0+(ylimt-y0)/(2*c0)*(x+1/k+c0-((x+c0+1/k)^2-4*c0*x)^0.5) 其中y0=0.723; ylimt=0.330; c0=0.00002 xdata=[0.00005 0.00015 0.0002 0.00025 0.0003 0.00035 0.00045 0.00055 0.0006 0.00075 0.00085 0.00095 0.001 0.00115 0.00125 0.00145 0.00155 0.00175 0.00205 0.00245] ydata=[0.706 0.665 0.652 0.624 0.599 0.574 0.545 0.526 0.492 0.472 0.450 0.434 0.426 0.402 0.387 0.369 0.364 0.345 0.338 0.330] [ Last edited by 天边5328 on 2010-4-3 at 17:35 ] |
回复此楼» 猜你喜欢
论文投稿，期刊推荐
已经有6人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有3人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有6人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
sxlion811
金虫 (正式写手)
- 程序强帖: 3
- 应助: 12 (小学生)
- 贵宾: 0.128
- 金币: 1078.7
- 散金: 2188
- 红花: 3
- 帖子: 868
- 在线: 444.3小时
- 虫号: 479151
- 注册: 2007-12-14
- 专业: 食品科学
2楼2010-04-09 18:38:05
3楼2010-04-09 20:35:39
he866837
至尊木虫 (著名写手)
- 应助: 4 (幼儿园)
- 金币: 14170.4
- 散金: 184
- 红花: 1
- 帖子: 1040
- 在线: 441.5小时
- 虫号: 672062
- 注册: 2008-12-10
- 专业: 纳米医学
4楼2010-04-09 20:45:54
5楼2010-04-10 08:22:19
sgyin
铁杆木虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 5519.4
- 散金: 100
- 红花: 8
- 帖子: 680
- 在线: 279.4小时
- 虫号: 405447
- 注册: 2007-06-17
- 性别: GG
- 专业: 材料物理化学
★ ★ ★ ★
jjdg(金币+2):辛苦了 2010-04-11 03:21
wangen994(金币+2):活动期间额外奖励 2010-04-11 20:19
jjdg(金币+2):辛苦了 2010-04-11 03:21
wangen994(金币+2):活动期间额外奖励 2010-04-11 20:19
|
ff=inline('0.723+(0.330-0.723)./(2.*0.00002).*(x+1./k(1)+0.00002-((x+0.00002+1./k(1)).^2-4.*0.00002.*x).^0.5)','k','x') [xx,res]=lsqcurvefit(ff,[1e3],xdata,ydata);  这是固定y0=0.723; ylimt=0.330; c0=0.00002 的结果！拟合只能这样了，楼主把这些参数放开试试 [ Last edited by sgyin on 2010-4-10 at 23:49 ] |
6楼2010-04-10 23:45:48
sgyin
铁杆木虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 5519.4
- 散金: 100
- 红花: 8
- 帖子: 680
- 在线: 279.4小时
- 虫号: 405447
- 注册: 2007-06-17
- 性别: GG
- 专业: 材料物理化学
|
ff=inline('k(2)+(k(3)-k(2))./(2.*0.00002).*(x+1./k(1)+0.00002-((x+0.00002+1./k(1)).^2-4.*0.00002.*x).^0.5)','k','x') [xx,res]=lsqcurvefit(ff,[1e3,1,1],xdata,ydata); fplot(@(x)ff(xx,x),[0,0.003],'r'),hold on,plot(xdata,ydata,'o')  这是拟合结果k= [1295.69223047089 0.756467370036168 0.171649191261954] |
7楼2010-04-10 23:54:51
sgyin
铁杆木虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 5519.4
- 散金: 100
- 红花: 8
- 帖子: 680
- 在线: 279.4小时
- 虫号: 405447
- 注册: 2007-06-17
- 性别: GG
- 专业: 材料物理化学
★ ★ ★ ★ ★
天边5328(金币+50):非常感谢您 2010-04-11 07:42
wangen994(金币+5):活动期间额外奖励 2010-04-11 10:52
天边5328(金币+50):非常感谢您 2010-04-11 07:42
wangen994(金币+5):活动期间额外奖励 2010-04-11 10:52
|
ff=inline('k(2)+(k(3)-k(2))./(2.*k(4)).*(x+1./k(1)+k(4)-((x+k(4)+1./k(1)).^2-4.*k(4).*x).^0.5)','k','x') [xx,res]=lsqcurvefit(ff,[1e3,1,1,0.0001],xdata,ydata); fplot(@(x)ff(xx,x),[0,0.003],'r'),hold on,plot(xdata,ydata,'o')  k= [ 3599.30404451002 0.73350054683102 0.263275701253758 0.000684697129278193] |
8楼2010-04-10 23:59:33
9楼2010-04-11 07:58:56
sgyin
铁杆木虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 5519.4
- 散金: 100
- 红花: 8
- 帖子: 680
- 在线: 279.4小时
- 虫号: 405447
- 注册: 2007-06-17
- 性别: GG
- 专业: 材料物理化学
10楼2010-04-11 12:13:16