| 查看: 1251 | 回复: 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代码 谢谢！ |
回复此楼» 猜你喜欢
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有180人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有6人回复
-
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有22人回复
-
考博自荐
已经有4人回复
-
山东大学第二批博士研究生招生
已经有0人回复
-
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【公式拟合问题】如何拟合下面这样的一个公式，多谢大牛。
已经有6人回复
- 非线性回归参数..急求帮助！！谢谢
已经有8人回复
- 【求助】注册化工工程师的报考条件，报考时间以及考试科目都是什么？
已经有26人回复
- 【求助】关于SPSS多元非线性回归的问题
已经有3人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
- 【求助】大家谈谈对注册环保工程师的了解
已经有15人回复
- 【求助】什么样的数据适合主成分分析
已经有8人回复
» 抢金币啦！回帖就可以得到:
-
湖北人在北京
+1/258
-
东南大学废弃碳资源利用-杰青团队-27级博士生(提前入学)-及博后、科研助理招募
+1/181
-
工科男，工作稳定，希望能遇到有趣的她
+1/174
-
东南大学废弃碳资源利用团队（杰青/海优团队）27级博士研究生和科研助理招生
+1/130
-
西澳大学 环境工程 全奖博士 2027 CSC博士招生（接收联培培养博士生）
+2/118
-
太行山徒步之凌水河（一）
+1/96
-
湖南师范大学化学化工学院李斌课题组招收2026级博士生
+1/89
-
中科大 冯伟 课题组招生-机器人界面控制、电学机械界面设计
+1/84
-
山东征女友，坐标济南
+1/73
-
电子科技大学（深圳）高等研究院陈明豪课题组招收推免硕士研究生
+1/62
-
华南理工大学环境与能源学院液流电池博士招生
+1/60
-
水处理技术员
+1/30
-
南京邮电大学李巍教授招收2026博士生（5月12日前有效）
+1/21
-
南京林业大学-国家级青年人才团队 招2026级博士 （合成化学）
+1/11
-
河南省医学科学院王宁利院士科研团队2026年博士、博后、硕士招聘
+1/9
-
大湾区大学-中山大学联培博士招生（电磁方向）
+1/2
-
天津医科大学基础医学院招收博士研究生
+1/2
-
南方科技大学-珞方生物（深圳）联合实验室 招聘高分子方向科研助理
+1/2
-
985院长团队！（大连理工，光电融合芯片方向） 招收2026秋入学博士生 5月初截止！
+1/1
-
中南林院士团队招招生物质能源与材料方向博士生
+1/1
5楼2011-03-17 11:35:30
★
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