版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

>论坛更新日志 (5049)
>导师招生 (571)
>文献求助 (449)
>论文投稿 (236)
>学术会议 (222)
>留学DIY (207)
>虫友互识 (203)
>硕博家园 (193)
>基金申请 (184)
>博后之家 (161)
>考博 (158)
>招聘信息布告栏 (140)
>休闲灌水 (124)
>育儿交流 (71)
>论文道贺祈福 (71)
>教师之家 (47)

返回列表

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表达式。

回复此楼

» 猜你喜欢

基金申请已经有44人回复
CSC与新西兰维多利亚大学PhD奖学金项目已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有213人回复
新西兰Robinson研究所招聘CSC公派访问人员已经有0人回复
帮我的英语口语老师找学生已经有0人回复
什么时候开奖？已经有13人回复
散金币祈福已经有94人回复
青基已中已经有1人回复
散金币祈福已经有185人回复
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布已经有1人回复

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

跪求： matlab三元非线性回归，请指教。。。已经有3人回复
跪求matlab 非线性回归程序已经有5人回复

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

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

回帖支持 ( 显示支持度最高的前 50 名 )

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.7小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

感谢参与，应助指数 +1

"NLLS"是否就是指"NonLinear Least Squares"？即非线性最小二乘法，试试1stOpt吧，不需要初值，简单易懂。

赞一下(1人)

回复此楼

2楼2015-03-06 09:42:43

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

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68171.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3292.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

感谢参与，应助指数 +1

http://muchong.com/bbs/viewthread.php?tid=6572963

赞一下(1人)

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

3楼2015-03-06 10:54:10

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

普通回帖

spinach45

捐助贵宾 (初入文坛)

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

好，谢谢楼上两位！我明天趁周末先捣弄一下。

赞一下

回复此楼

4楼2015-03-07 23:05:34

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

spinach45

捐助贵宾 (初入文坛)

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

引用回帖:

2楼: Originally posted by dingd at 2015-03-06 09:42:43
"NLLS"是否就是指"NonLinear Least Squares"？即非线性最小二乘法，试试1stOpt吧，不需要初值，简单易懂。

我已经用NLLS试过了，但算不出和看的论文一样的结果。你有最新版本的1stOpt吗？能不能帮我运行一下？

赞一下

回复此楼

5楼2015-03-10 11:24:46

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

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68171.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3292.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

引用回帖:

5楼: Originally posted by spinach45 at 2015-03-10 11:24:46
我已经用NLLS试过了，但算不出和看的论文一样的结果。你有最新版本的1stOpt吗？能不能帮我运行一下？

Parameters b1, b2, b3;
Variable a, x, y, z;
Constant R=1.073533812;
Function R*a=(1+exp(b1*(b2-b3 ...

比较一下你算出结果的相关系数和论文给出的，论文给出的计算结果，不一定是最优的，对文献不能太相信。

赞一下

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

6楼2015-03-10 11:27:07

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

spinach45

捐助贵宾 (初入文坛)

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

引用回帖:

2楼: Originally posted by dingd at 2015-03-06 09:42:43
"NLLS"是否就是指"NonLinear Least Squares"？即非线性最小二乘法，试试1stOpt吧，不需要初值，简单易懂。

之前那个回帖不知道为什么显示不了。
我把我的编码再发一次，并附上数据。
另外，其实a, x, y, z是xt, xt-1, xt-2, xt-3的时间序列，是不是需要因为这样而做另外的处理？如果是的话，怎么处理？
$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; $

赞一下

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : 发问excel.xls

2015-03-10 11:31:19, 27.5 K

7楼2015-03-10 11:32:42

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

spinach45

捐助贵宾 (初入文坛)

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

引用回帖:

6楼: Originally posted by 月只蓝 at 2015-03-10 11:27:07
比较一下你算出结果的相关系数和论文给出的，论文给出的计算结果，不一定是最优的，对文献不能太相信。...

你能帮我看看么？我1stopt里做出的和论文不大一样，也做不出论文里面其他的处理（论文里的几个变量其实是xt, xt-1, xt-2, xt-3的时间序列，是不是要因此做特殊处理）。（数据省略，在xls文件里）
$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;$
我根据你3楼的那个帖子写了个Matlab代码，但最后被说是“输入参数的数目不足”，不知道你是否能帮我看看呢？（数据省略，在xls文件里）
$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))); $

最后，这是论文本身的式子和结果截图。
$\"初学新手求教。求用Matlab求解一个非线性回归，或请指点门路。-1\"$
$\"初学新手求教。求用Matlab求解一个非线性回归，或请指点门路。-2\"$
$\"初学新手求教。求用Matlab求解一个非线性回归，或请指点门路。-3\"$ $\"初学新手求教。求用Matlab求解一个非线性回归，或请指点门路。-4\"$

赞一下

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : 发问excel.xls

2015-03-10 11:41:02, 27.5 K

8楼2015-03-10 11:41:51

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

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.7小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

1stOpt的代码没显示出来，附件形式再发一次。

赞一下

回复此楼

9楼2015-03-10 11:47:17

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

spinach45

捐助贵宾 (初入文坛)

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

引用回帖:

9楼: Originally posted by dingd at 2015-03-10 11:47:17
1stOpt的代码没显示出来，附件形式再发一次。

好的，已经发了txt了。请查收。

赞一下

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : 发问.txt

2015-03-10 11:51:53, 6.69 K

10楼2015-03-10 11:52:15

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

相关版块跳转我要订阅楼主 spinach45 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3	ero8OE6tv9cu 2026-09-02	5/250	2026-09-02 22:50 by rM1TE0WVDIIY
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3	ero8OE6tv9cu 2026-09-02	3/150	2026-09-02 22:50 by rM1TE0WVDIIY
[基金申请] 科研人应该花精力去思考如何解决问题，而不是去凝练问题 +7	瞬息宇宙 2026-09-01	14/700	2026-09-02 17:33 by ma0526
[基金申请] 要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够 +9	瞬息宇宙 2026-08-31	17/850	2026-09-02 14:04 by anjeeshine
[基金申请] 学科评审组评审是指会评吗？ +5	瞬息宇宙 2026-08-31	5/250	2026-09-02 10:13 by 雨冰共舞
[基金申请] 面上合作单位盖章 +6	ssyjh 2026-08-27	9/450	2026-09-01 20:01 by huagongfeihu
[基金申请] 面上函评意见出来了，像什么等级？ 20+4	Tsingking1 2026-08-27	17/850	2026-09-01 19:51 by 超级无敌华子
[论文投稿] 小白求助投论文要求的highlights应该如何写 5+3	l1963982152 2026-08-29	4/200	2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +7	gdd2018 2026-08-28	12/600	2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 国社科又开始会评了，不知道这次命运如何 +7	雨打竹帘 2026-08-30	11/550	2026-08-31 23:16 by hittle2008
[基金申请] 基金未中，这种答复是模板吗？ +6	zhaosm1982 2026-08-27	7/350	2026-08-31 21:18 by qdxxmc
[基金申请] 能否申诉？ +7	echo8914667 2026-08-30	8/400	2026-08-31 17:00 by yihongxu
[基金申请] 中青基了要发朋友圈吗？ +7	349506619 2026-08-28	7/350	2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4	笨笨唐 2026-08-28	4/200	2026-08-31 09:30 by huixian257
[基金申请] 为什么到现在没收到通知？ +5	tannykie 2026-08-29	5/250	2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7	德尚中行 2026-08-27	8/400	2026-08-30 20:52 by purplejack
[考博] 找导师 +6	yuanjiabao 2026-08-29	7/350	2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +6	gouxfjh 2026-08-28	11/550	2026-08-30 07:57 by gouxfjh
[基金申请] 国自然评审意见 +13	wangmingqi 2026-08-28	19/950	2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4	bud_bud 2026-08-27	7/350	2026-08-29 07:23 by foolishmani

24小时热门版块排行榜

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

» 猜你喜欢

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

【答案】应助回帖

【答案】应助回帖

【答案】应助回帖

» 本帖附件资源列表

» 本帖附件资源列表

【答案】应助回帖

» 本帖附件资源列表

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