| 查看: 1197 | 回复: 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代码 谢谢！ |
回复此楼» 猜你喜欢
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组招收2026硕士研究生1名
已经有0人回复
-
山西大同大学物理学硕士研究生招收调剂生
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有242人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
［调剂信息］211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
一志愿物理所凝聚态物理352求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有2人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
-
070200求调剂，一志愿某211，288分
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【公式拟合问题】如何拟合下面这样的一个公式，多谢大牛。
已经有6人回复
- 非线性回归参数..急求帮助！！谢谢
已经有8人回复
- 【求助】注册化工工程师的报考条件，报考时间以及考试科目都是什么？
已经有26人回复
- 【求助】关于SPSS多元非线性回归的问题
已经有3人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
- 【求助】大家谈谈对注册环保工程师的了解
已经有15人回复
- 【求助】什么样的数据适合主成分分析
已经有8人回复
» 抢金币啦！回帖就可以得到:
-
延安大学化学与化工学院接收化学专业调剂生（工科、理学均可调剂）
+5/160
-
欢迎调剂佛山大学资源与环境专业硕士
+2/142
-
上海交通大学化学化工学院膜分离-纳流控课题组诚聘博士后2-3名
+1/91
-
某211招收光学工程 仪器科学与技术 tiaoji 研究生 若干名！！！
+1/48
-
河南工业大学化学化工学院 2026 年硕士调剂正在进行！
+1/39
-
新疆碳纳米材料制备与应用重点实验室2026年招聘启事
+1/37
-
三峡大学国家青年人才课题组招2026年化学/材料/化工类硕士调剂生
+1/37
-
河南中部某本科高校 资源与环境（专硕）研究生 招生信息
+1/37
-
【博士招生】武汉科技大学招收材料、化工、环境、冶金类2026年“申请-考核”制博士生
+1/13
-
上海工程技术大学李聪玲课题组招收调剂生
+1/8
-
上海应用技术大学化工学院招收调剂
+1/5
-
丽水学院材料与能源工程学院 材料与化工、资源与环境专业2026年硕士研究生接受调剂
+1/5
-
西北农林科技大学水产养殖专业学术型硕士调剂名额1个
+1/4
-
军事医学研究员毒物药物研究所药理学生物技术药物代谢方向招收调剂生 4名
+1/4
-
湖北师范大学招调剂啦
+1/3
-
齐齐哈尔大学李莉课题组诚招2026级考研调剂生（学硕和专硕）
+1/2
-
大连工业大学轻工与化学工程学院招工科调剂生
+2/2
-
湖南工业大学生医学院接受有机化学、生物与医药等方向的硕士调剂生
+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
