24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

一直在纠结

铁杆木虫 (正式写手)

应助: 87 (初中生)
金币: 5516
散金: 190
红花: 20
帖子: 681
在线: 466.9小时
虫号: 2320219
注册: 2013-03-05
专业: 分离过程

[求助] 用MATLAB回归方程的参数，是个隐函数，请附上程序已有3人参与

数据和方程如下，
x=
53.428 54.076 55.046 57.571 78.753 79.154 79.485 80.022 80.488 80.938 81.162 81.387 81.613 82.068 82.296 82.524 83.209 83.436 84.116 84.344 84.809 85.048 85.544 85.989 86.577 87.222 87.987 88.951 89.272 89.832 91.809 101.409 138.539
y=
1.29E+01 1.30E+01 1.32E+01 1.37E+01 1.93E+01 1.93E+01 1.93E+01 1.93E+01 1.93E+01 1.93E+01 1.93E+01 1.92E+01 1.92E+01 1.92E+01 1.92E+01 1.92E+01 1.92E+01 1.91E+01 1.91E+01 1.91E+01 1.91E+01 1.91E+01 1.89E+01 1.85E+01 1.84E+01 1.83E+01 1.83E+01 1.83E+01 1.83E+01 1.84E+01 1.88E+01 2.14E+01 3.36E+01
方程如下：
logy=a1+a2/(a3+x)+a4*x+a5*x^2+a6*logx+a7*y/x^2
其中a1~a7是方程要回归的7个参数
能帮忙回归出来就可以，不限于用MATLAB

但是请把解决的程序附上，多谢！！！

回复此楼

» 猜你喜欢

求国际会议网站已经有1人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有217人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有19人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
26申博推荐：南京航空航天大学国际前沿科学研究院光学方向招收博士生！已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼2014-06-28 14:41:08

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

一直在纠结

铁杆木虫 (正式写手)

应助: 87 (初中生)
金币: 5516
散金: 190
红花: 20
帖子: 681
在线: 466.9小时
虫号: 2320219
注册: 2013-03-05
专业: 分离过程

引用回帖:

5楼: Originally posted by Holly2012 at 2014-06-28 20:09:04
参考 http://muchong.com/bbs/viewthread.php?tid=3240154 一楼有完整代码，只要修改一下函数和数据就可以

呃  仿照他的输入了一下有个错误但是不知道是哪里不行  能帮忙看一下吗
p=[1.29E+01       1.30E+01       1.32E+01       1.37E+01       1.93E+01       1.93E+01       1.93E+01       1.93E+01       1.93E+01       1.93E+01       1.93E+01       1.92E+01       1.92E+01       1.92E+01       1.92E+01       1.92E+01       1.92E+01       1.91E+01       1.91E+01       1.91E+01       1.91E+01       1.91E+01       1.89E+01       1.85E+01       1.84E+01       1.83E+01       1.83E+01       1.83E+01       1.83E+01       1.84E+01       1.88E+01       2.14E+01       3.36E+01];
c=[53.428       54.076       55.046       57.571       78.753       79.154       79.485       80.022       80.488       80.938       81.162       81.387       81.613       82.068       82.296       82.524       83.209       83.436       84.116       84.344       84.809       85.048       85.544       85.989       86.577       87.222       87.987       88.951       89.272       89.832       91.809       101.409       138.539];
cc=c';
pp=p';
F=@(q,x)(q(1)+q(2)/(q(3)+x(:,1))+q(4)*x(:,1)+q(5)*(x(:,1)^2)+q(6)*log(x(:,1))+q(7)*x(:,2)/x(:,1)^2)-log(x(:,2));
x=[cc pp];
q0=[1 1    1 1 1 1 1];
warning off
q=nlinfit(x,zeros(size(x,1),1),F,q0);
disp(num2str(q));
显示错误，应该是方程输入错误，但是不知道错误出在了哪里