| 查看: 1176 | 回复: 5 | |||
[交流]
三个变量之间的非线性回归问题求助
|
|
三个变量之间的非线性回归问题求助 已知三个变量之间的函数关系为z=(p1*p2*p3*x*y^0.5)/[(1+p2*x)(1+p3*y^0.5)],且知道n组x,y,z的值,能否利用matlab等数学工具非线性回归原函数,求出其中的参数p1,p2,p3呢? 数据如下 x y z 9.9621 89.9837 0.00585 9.0558 90.8959 0.00574 8.3014 91.6560 0.00552 7.6633 92.2989 0.00531 7.1158 92.8498 0.00520 请告诉我matlab代码 谢谢! |
» 猜你喜欢
帮我的英语口语老师找学生
已经有4人回复
计算三阶力常数用thirdorder扩胞时参数如何选择?临近原子数怎么取值
已经有7人回复
物理学I论文润色/翻译怎么收费?
已经有186人回复
2026年循环经济功能材料国际会议(ICFMCE 2026)
已经有0人回复
2026年第五届电气、电子与信息工程国际会议(ISEEIE 2026)
已经有0人回复
哈尔滨理工大学物理系招收物理学考研调剂
已经有5人回复
0702一志愿吉大B区求调剂
已经有5人回复
求调剂
已经有0人回复
0702一志愿吉大B区求调剂有论文
已经有0人回复
请问还有没有用Latex写文章的小伙伴们?
已经有0人回复
» 本主题相关价值贴推荐,对您同样有帮助:
【公式拟合问题】如何拟合下面这样的一个公式,多谢大牛。
已经有6人回复
非线性回归参数..急求帮助!!谢谢
已经有8人回复
【求助】注册化工工程师的报考条件,报考时间以及考试科目都是什么?
已经有26人回复
【求助】关于SPSS多元非线性回归的问题
已经有3人回复
【求助】 求MAtlab求解程序!!!
已经有20人回复
【求助】大家谈谈对注册环保工程师的了解
已经有15人回复
【求助】什么样的数据适合主成分分析
已经有8人回复
» 抢金币啦!回帖就可以得到:
南京医科大学生殖医学与子代健康国重实验室-董飞宏教授课题组-硕士招生
+1/107
郑州大学付永柱课题组诚招2026年入学博士生
+1/80
厦门大学电池工况表征博士后招聘(固态电池 / 层状氧化物正极 / 原位-工况表征方向)
+1/78
天津工业大学航空航天博士招生
+1/35
上海交通大学电催化方向博后招聘-有锂电、高分子、燃料电池背景者优先
+1/35
南开大学罗景山教授课题组招聘博士后(光/电催化、钙钛矿光伏方向)
+2/34
中科院生态环境研究中心国重实验室招聘客座研究生1-2名
+1/31
中山大学院士团队微纳器件、脑机接口方向博士招生、博后招聘
+1/27
国家级领军人才团队高分子材料方向2026年博士研究生招生
+1/16
2026年西南科技大学功能涂层课题组简介
+1/12
2026年上海大学王亮课题组招收电催化或光催化研究背景的博士生
+2/10
西京学院土木水利 2026 级研究生招生相关说明
+1/8
教育部长江学者和创新团队发展计划”入选团队招收 材料,化学与化工博士研究生
+1/7
青岛科技大学可持续高分子团队 考研招生
+1/7
上海第二工业大学-朱大海课题组招生(过线就能调剂)
+1/5
南方科技大学生命科学学院基础免疫与微生物学系招收博士后
+1/5
中国石油大学(北京)招收机械设计(形状记忆合金驱动机构)方向博士生
+1/3
杭州北航国新院徐梁教授招聘博士(控制/机器视觉/机器人导航操作方向)
+1/1
温州大学招收2026年入学博士研究生(化学、材料、环境)
+1/1
澳洲维多利亚大学计算机类全额奖学金博士招生 (邮箱+微信可联系)
+1/1
★
zfym(金币+1):谢谢参与
zfym(金币+9): 谢谢了,太好了,麻烦你了! 2011-03-17 10:13:44
zfym(金币+1):谢谢参与
zfym(金币+9): 谢谢了,太好了,麻烦你了! 2011-03-17 10:13:44
|
function ccccc clear clc format long data=[9.9621 89.9837 0.00585 9.0558 90.8959 0.00574 8.3014 91.6560 0.00552 7.6633 92.2989 0.00531 7.1158 92.8498 0.00520]; x=data(:,1:2); z=data(:,3); beta0 = [0.0088 0.2012 45.7356]; % 参数初值 options=optimset('TolFun',1e-8,'display','off'); [beta,resnorm,residual,EXITFLAG] = lsqnonlin(@Objfun,beta0,[],[],options,x,z); zcal=myfun(x,beta); % 模型计算值 % 结果输出 fprintf('\tp1=%10.6f\n\tp2=%10.6f\n\tp3=%10.6f\n',beta) fprintf('\tresidual=\n') fprintf('\t%10.6f\n',residual) fprintf('\tresnorm=%10.6f\n',resnorm) fprintf('\n\tEXITFLAG=%d\n',EXITFLAG) figure(1) plot(z,zcal,'b.') % 实验值与模型计算值比较 xlabel('z 实验值') ylabel('z 计算值') hold on refline([1 0]) % 对角线 figure(2) plot(z,residual./z*100,'ro') % 实验值与误差值绘图 ylim([-10 10]) xlabel('z') ylabel('相对误差 %') hold on refline([0 0]) % 误差为0的线 function f=Objfun(beta,lamda,z) % 目标函数,sum((T(i)-Tcal(i))^2)最小化 zcal=myfun(lamda,beta); f=z-zcal; function z=myfun(x,beta) % 模型 p1=beta(1); p2=beta(2); p3=beta(3); z=(p1.*p2.*p3.*x(:,1).*x(:,2).^0.5)./((1+p2.*x(:,1)).*(1+p3*x(:,2).^0.5)); p1= 0.008833 p2= 0.201290 p3= 45.735600 residual= -0.000030 0.000049 0.000007 -0.000037 0.000010 resnorm= 0.000000 EXITFLAG=1 |
2楼2011-03-17 08:50:49
3楼2011-03-17 10:25:06
4楼2011-03-17 10:39:07
5楼2011-03-17 11:35:30
zfym(金币+3): 谢谢。 2011-03-19 23:47:16
|
本帖内容被屏蔽 |
6楼2011-03-17 15:07:38













回复此楼