| 查看: 807 | 回复: 6 | ||
| 【悬赏金币】回答本帖问题，作者100121114将赠送您 80 个金币 | ||
100121114金虫 (小有名气)
|
[求助]
求助大神帮忙做下数据拟合已有1人参与
|
|
|
我用matlab做拟合： xdata=[0.1 0.5 1 2 3 4 5 6 7 8 9 10 12 15 18 20]; ydata=[0.423,0.585,0.682,0.748,0.786,0.813,0.837,0.859,0.878,0.897,0.914,0.929,0.963,1.01,1.048,1.053]； 拟合公式为y=0.0257*ln(x/a)+b*x;其中啊a，b是需要拟合的参数；@月只蓝@beefly |
回复此楼» 猜你喜欢
氧化硅薄膜玻璃 求购
已经有3人回复
-
vasp极化强度
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有250人回复
-
vasp计算的一些问题
已经有0人回复
-
国科温州研究院 项目聘用有意“二战”的优秀本科毕业生
已经有31人回复
-
SCAPS-1D学习交流
已经有1人回复
-
深圳大学跟哈工大联合博士
已经有6人回复
★ ★ ★ ★ ★
独孤神宇: 金币+5, 鼓励交流 2019-08-07 12:07:00
独孤神宇: 金币+5, 鼓励交流 2019-08-07 12:07:00
|
xdata=[0.1 0.5 1 2 3 4 5 6 7 8 9 10 12 15 18 20]; ydata=[0.423,0.585,0.682,0.748,0.786,0.813,0.837,0.859,0.878,0.897,0.914,0.929,0.963,1.01,1.048,1.053]; f=inline('0.0275*(log(x)-a(1))+a(2)*x','a','x'); [q,r]=lsqcurvefit(f,[10,5],xdata,ydata); a = exp(q(1)); b = q(2); x = 0 : 0.1 : 20; y = 0.0275*log(x/a)+b*x; plot(xdata,ydata,'.'); hold on plot(x,y,'r'); 发自小木虫Android客户端 |
» 本帖已获得的红花（最新10朵）
1001211142楼2019-08-07 11:11:46
100121114
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1491.6
- 散金: 30
- 帖子: 65
- 在线: 21小时
- 虫号: 3202917
- 注册: 2014-05-13
- 专业: 半导体科学与信息器件
3楼2019-08-07 11:32:38
hzlhm
至尊木虫 (著名写手)
- 应助: 373 (硕士)
- 金币: 16171.8
- 红花: 52
- 帖子: 2569
- 在线: 572.1小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
4楼2019-08-07 11:47:53
hzlhm
至尊木虫 (著名写手)
- 应助: 373 (硕士)
- 金币: 16171.8
- 红花: 52
- 帖子: 2569
- 在线: 572.1小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
【答案】应助回帖
★ ★ ★ ★ ★
独孤神宇: 金币+5, 鼓励交流 2019-08-07 12:07:17
独孤神宇: 金币+5, 鼓励交流 2019-08-07 12:07:17
|
xdata=[0.1 0.5 1 2 3 4 5 6 7 8 9 10 12 15 18 20]; ydata=[0.423,0.585,0.682,0.748,0.786,0.813,0.837,0.859,0.878,0.897,0.914,0.929,0.963,1.01,1.048,1.053]; fun=@(a,x)a(1)*log(x/a(2))+a(3)*x; x0=[0.098 0.001 0.005]; [a,r,J] = nlinfit(xdata,ydata,fun,x0); a |
» 本帖已获得的红花（最新10朵）
5楼2019-08-07 11:49:06
100121114
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1491.6
- 散金: 30
- 帖子: 65
- 在线: 21小时
- 虫号: 3202917
- 注册: 2014-05-13
- 专业: 半导体科学与信息器件
6楼2019-08-07 21:54:16
hzlhm
至尊木虫 (著名写手)
- 应助: 373 (硕士)
- 金币: 16171.8
- 红花: 52
- 帖子: 2569
- 在线: 572.1小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
100121114(独孤神宇代发): 金币+20, 鼓励交流 2019-08-09 08:23:36
100121114(独孤神宇代发): 金币+20, 鼓励交流 2019-08-09 08:23:36
|
如那个等于0.0257系数是物理条件设定，那么根据你数据其拟合精度就要差很多。 拟合结果： a = 7.89169947611115e-12 b= 0.0192094598725463 R2 = 0.936354194940639 代码： xdata=[。。。。。。]; ydata=[。。。。。。]; a0=[1e-12 0.019]; fun=@(a,x)0.0257*log(x/a(1))+a(2)*x; [a,r] = nlinfit(xdata,ydata,fun,a0) |
7楼2019-08-08 22:45:17
5