24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

spinach45

捐助贵宾 (初入文坛)

应助: 0 (幼儿园)
金币: 162
帖子: 9
在线: 3.1小时
虫号: 3713251
注册: 2015-03-05

[求助] 初学新手求教。求用Matlab求解一个非线性回归，或请指点门路。已有2人参与

第一次上论坛发帖，请多多指教啊。
想请教一下下面这个回归用怎样用Matlab编程？

论文里面写着用NLLS解决。
这个怎么用Matlab编程啊，初始值应该怎么确定啊。自己研究了三天，但发觉没有方向，还想请教一下各位高手，应该怎样系统地看什么书啊？

我研究了三天，把变量变为时间序列后，自己用了这样编程。

CODE:

X = [R1,Xt_1,Xt_2,Xt_3];

Y = y;

modelfun = @(1+exp(b(1)*(b(2)-b(3))*X(:,4)*(X(:,2)-X(:,1)*X(:,3)))).^(-1).*b(2)*X(:,2)+(1-(1+exp(b(1)*(b(2)-b(3))*X(:,4)*(X(:,2)-X(:,1)*X(:,3)))).^(-1)).*b(3)*X(:,2);

beta0=[-10, 0.5, 0.5];

mdl = NonLinearModel.fit(X,Y,modelfun,beta0);

modelfun那行要不就是括号不对称要不就是不应为Matlab表达式。

回复此楼

» 猜你喜欢

求2026年在台湾举行的物理和材料领域国际学术会议信息已经有0人回复
求国际会议网站已经有1人回复
物理学I论文润色/翻译怎么收费? 已经有164人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有19人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复

1楼2015-03-05 22:48:14

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

spinach45

捐助贵宾 (初入文坛)

应助: 0 (幼儿园)
金币: 162
帖子: 9
在线: 3.1小时
虫号: 3713251
注册: 2015-03-05

$\"初学新手求教。求用Matlab求解一个非线性回归，或请指点门路。-1\"$
$\"初学新手求教。求用Matlab求解一个非线性回归，或请指点门路。-2\"$
$\"初学新手求教。求用Matlab求解一个非线性回归，或请指点门路。-3\"$ $\"初学新手求教。求用Matlab求解一个非线性回归，或请指点门路。-4\"$

各位，实在抱歉，之前回的帖子图都挂了，式子也看不到。上图是原文的论述和式子。我根据原文提到的数据用1stopt进行了回归，得出的结果不大相同，而且也不知道怎么得到其他例如AIC之类的结果。请问能不能帮我看看应该如何处理，或者用更高版本的1stopt帮忙走走程序？
我写的1stopt程序（省略了data，参见附件）

CODE:

Parameters b1, b2[0,1], b3[1,2];

Variable a, x, y, z;

Constant R=1.073533812;

Function R*a=(1+exp(b1*(b2-b3)*z*(x-R*y)))^(-1)*b2*x+(1-(1+exp(b1*(b2-b3)*z*(x-R*y)))^(-1))*b3*x;

data;

另外，我在尝试用matlab的时候也出现了问题，说“输入参数的数目不足”，还有原因，“原因: Failure in initial user-supplied objective function evaluation. LSQNONLIN cannot continue”，昨晚查了一晚也不知道怎么回事，能麻烦帮忙看看吗？（省略了data，参见附件）

CODE:

function fit_nonl

clear all;clc

format long

data=[

];

R=1.073533812;

Xt_1=data(:,2);

Xt_2=data(:,3);

Xt_3=data(:,4);

Xt=data(:,1);

X=[Xt_1, Xt_2, Xt_3];

Y=R*Xt;

b0=[-10; 0.5; 1.5];

lb=[-100; 0; 1]*1e9;

ub=[0; 1; 2]*1e9;

%-------------------------------------------------------------------------

% 使用函数lsqnonlin()进行参数估计

OPTIONS=optimset('MaxFunEvals',1000);

[b,resnorm,residual,jacobian] = ...

    lsqnonlin(@ObjFunc,b0,lb,ub,OPTIONS,X,Y);

ci = nlparci(b,residual,jacobian);

%residual;

%-------------------------------------------------------------------------

function f = ObjFunc(b,X,Y)

f=KineticsEqs(X,b)-Y;

%------------------------------------------------------------------------

function xt = KineticsEqs(b,X)

b1=b(2); 

b2=b(3);

x1=X(:,1);

xt=(myfunc)^(-1)*b1*x1+(1-(myfunc)^(-1))*b2*x1;

function n=myfunc(b,X)

R=1.073533812;

B=b(1); 

b1=b(2); 

b2=b(3);

n=1+exp(B*(b1-b2)*X(:,3)*(X(:,2)-R*X(:,1)));