24小时热门版块排行榜

返回列表

zfym

铜虫 (正式写手)

应助: 1 (幼儿园)
金币: 884.4
帖子: 556
在线: 283.1小时
虫号: 854647

[交流] 三个变量之间的非线性回归问题求助

三个变量之间的非线性回归问题求助

已知三个变量之间的函数关系为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论文润色/翻译怎么收费? 已经有182人回复
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）已经有6人回复
【修正版】单环用磁约束低速电子实验方案（简化版）已经有0人回复
桂林理工大学物理学专业招收调剂，还有三个名额！！！已经有22人回复
考博自荐已经有4人回复
山东大学第二批博士研究生招生已经有0人回复
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学已经有0人回复
《电磁学》教材推荐已经有1人回复
【急招】合肥工大核聚变材料计算方向2026级工程博士生已经有4人回复

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

【公式拟合问题】如何拟合下面这样的一个公式，多谢大牛。已经有6人回复
非线性回归参数..急求帮助！！谢谢已经有8人回复
【求助】注册化工工程师的报考条件，报考时间以及考试科目都是什么？已经有26人回复
【求助】关于SPSS多元非线性回归的问题已经有3人回复
【求助】求MAtlab求解程序！！！已经有20人回复
【求助】大家谈谈对注册环保工程师的了解已经有15人回复
【求助】什么样的数据适合主成分分析已经有8人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2011-03-16 20:57:49

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

change0618

铁杆木虫 (著名写手)

★
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

赞一下(1人)

回复此楼

2楼2011-03-17 08:50:49

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

zfym

铜虫 (正式写手)

应助: 1 (幼儿园)
金币: 884.4
帖子: 556
在线: 283.1小时
虫号: 854647

引用回帖:

Originally posted by change0618 at 2011-03-17 08:50:49:
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= ...

??? function ccccc
|
Error: Function definitions are not permitted in this context.

该怎么办

赞一下

回复此楼

3楼2011-03-17 10:25:06

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