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代码
谢谢！

回复此楼

» 猜你喜欢

26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！已经有1人回复
如何从铁电相到顺电相。已经有1人回复
物理学I论文润色/翻译怎么收费? 已经有199人回复
各位大佬，求一份最新Wien2k版本作为学习用途已经有0人回复
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后已经有41人回复
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向已经有24人回复
PRB投稿，Acknowledgment sent to author状态十天了已经有3人回复
PVA溶解已经有1人回复
基金申请已经有44人回复
CSC与新西兰维多利亚大学PhD奖学金项目已经有0人回复
新西兰Robinson研究所招聘CSC公派访问人员已经有0人回复

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

【公式拟合问题】如何拟合下面这样的一个公式，多谢大牛。已经有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的回帖