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人回复
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标已经有10人回复
物理学I论文润色/翻译怎么收费? 已经有120人回复
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名已经有1人回复
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生已经有5人回复
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究已经有12人回复
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！已经有1人回复
如何从铁电相到顺电相。已经有1人回复

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

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