| 查看: 982 | 回复: 6 | |||
[交流]
matlab非线性回归 数据处理已有1人参与
|
|||
如图所示，想向各位请教一下，这样的数据如何实现非线性回归呢？式子中的d是表格中的因变量，变量为T温度以及t保温时间，d0是保温时间为0时候对应的数据。想要通过Matlab实现分线性拟合 求得A Q n三个数值。还望各位虫友不吝赐教啊。  d2888407a6f1569c6eff0c349589bc05.png  557ea5269e3cc8b5720b16b7c3d48147.png |
回复此楼» 猜你喜欢
with editor 两个月了，什么原因？
已经有12人回复
-
博士白读了
已经有42人回复
-
Chemical Engineering Journal投稿3周了，一直显示With editor状态。这是送审了吗？
已经有11人回复
-
"颜宁:基础研究应顶天立地"能做到基础研究同时顶天立地的才是牛人
已经有4人回复
-
普通院校药学硕士，做合成的，感觉找不到工作
已经有15人回复
-
申博求助
已经有10人回复
-
Nature一直在编辑手里，考虑好几天了，是悬了吗
已经有12人回复
-
刚博士毕业，入职普通一本高校
已经有19人回复
-
催稿后秒拒
已经有12人回复
-
取博导收留
已经有8人回复
小小企鹅2: 回帖置顶 2018-11-17 17:18:58
|
%%%第一段程序 function yhat=volum(beta,x) b0=beta(1); b1=beta(2); b2=beta(3); %b0为n b1为A b2为Q x0=x(:,1); x1=x(:,2); x2=x(:,3); %x0等于d0 x1等于保温时间 x2等于温度 y=(x0.^b0+b1*x1*exp(-1*b2./(8.314*x2)))^(1/b0); %%%第二段程序 x1=[0,0,0,0,0,0,300,300,300,300,300,300,600,600,600,600,600,600,900,900,900,900,900,900,1200,1200,1200,1200,1200,1200]; x2=[14.6,16,19.2,34,46.3,54.8,14.6,16,19.2,34,46.3,54.8,14.6,16,19.2,34,46.3,54.8,14.6,16,19.2,34,46.3,54.8,14.6,16,19.2,34,46.3,54.8]; x3=[850,950,1000,1050,1100,1500,850,950,1000,1050,1100,1500,850,950,1000,1050,1100,1500,850,950,1000,1050,1100,1500,850,950,1000,1050,1100,1500]; y=[14.6,16,19.2,34,46.3,54.8,14.1,16.9,20.9,36.1,48.6,61.2,14.8,17.7,21.8,37.5,50.7,65.5,14.2,18.3,23.6,38.1,54.2,68.8,15.4,18.6,24.3,40.2,55.7,72.7]; xx=[x1',x2',x3']; yy=[y']; beta0=[4,10000,50000]; [beta,r,j]=nlinfit(xx,yy,@volum,beta0); 但是总是报错，还请大家帮忙看一下哪里错了呢？还望赐教 |
7楼2018-11-15 15:07:23
hollya
至尊木虫 (文坛精英)
葛城?{美
- 应助: 744 (博后)
- 金币: 43440.6
- 散金: 1143
- 红花: 127
- 沙发: 779
- 帖子: 24904
- 在线: 834.8小时
- 虫号: 458790
- 注册: 2007-11-13
- 性别: MM
- 专业: 计算机软件
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
月只蓝: 金币+1, 鼓励交流 2018-11-15 08:50:07
小木虫: 金币+0.5, 给个红包，谢谢回帖
月只蓝: 金币+1, 鼓励交流 2018-11-15 08:50:07
|
%% data set of format 1 ddata=[14.6 14.1 14.8 14.2 15.4 16 16.9 17.7 18.3 18.6... 19.2 20.9 21.8 23.6 24.3 34 36.2 37.5 38.1 40.2... 46.3 48.6 50.7 54.2 55.7 54.8 61.2 65.5 68.6 72.7]' t=[0:300:1200 0:300:1200 0:300:1200 0:300:1200 0:300:1200 0:300:1200]' T=[850*ones(1,5) 950*ones(1,5) 1000*ones(1,5) 1050*ones(1,5)... 1100*ones(1,5) 1050*ones(1,5)]' nn=size(t) gscatter(t,ddata,T) xdata=[t T] hold on dfit=@(a,xdata)a(1)+a(2)*xdata(:,1).*exp(-a(3)./(8.314*xdata(:,2))) a=nlinfit(xdata,ddata,dfit,[1 1 1]) 能fit是能，但结果不太好，希望牛人出现帮一个忙。 |
2楼2018-11-15 06:13:54
。3楼2018-11-15 10:28:59
4楼2018-11-15 10:31:51
hollya
至尊木虫 (文坛精英)
葛城?{美
- 应助: 744 (博后)
- 金币: 43440.6
- 散金: 1143
- 红花: 127
- 沙发: 779
- 帖子: 24904
- 在线: 834.8小时
- 虫号: 458790
- 注册: 2007-11-13
- 性别: MM
- 专业: 计算机软件
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
jjdg: 金币+1, 感谢参与 2018-11-18 13:14:58
小木虫: 金币+0.5, 给个红包，谢谢回帖
jjdg: 金币+1, 感谢参与 2018-11-18 13:14:58
|
还有我觉得里面的T不对，不应该用摄氏度，应该用开尔文，也许会更接近你的数值 发自小木虫IOS客户端 |
5楼2018-11-15 11:13:49
6楼2018-11-15 13:42:10
10