24小时热门版块排行榜

>论坛更新日志 (2776)
>导师招生 (856)
>虫友互识 (538)
>休闲灌水 (408)
>招聘信息布告栏 (145)
>论文投稿 (133)
>文献求助 (112)
>博后之家 (92)
>学术会议 (51)
>硕博家园 (47)
>考博 (40)
>教师之家 (39)
>考研 (38)
>公派出国 (28)
>基金申请 (21)
>找工作 (17)

返回列表

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

回复此楼

» 猜你喜欢

散金币祈福已经有187人回复
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布已经有1人回复
物理学I论文润色/翻译怎么收费? 已经有74人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
深圳大学应用技术学院招聘凝聚态物理博士后已经有117人回复
深圳大学2027级光电信息工程/物理学研究生推免已经有70人回复
求助中医院临床中医中标的青年标书学习，谢谢已经有4人回复
浙大物院表面物理课题组长期招收普博及直博学生已经有257人回复
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry 已经有47人回复

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

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