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人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有62人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复

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

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

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

查看全部散金贴

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

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

change0618

铁杆木虫 (著名写手)

你的数组太少，而且你总要给出个参数的范围吧，有物理意义吧

赞一下

回复此楼

5楼2011-03-17 11:35:30

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

查看全部 6 个回答

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的回帖