24小时热门版块排行榜    

查看: 2361  |  回复: 17

努力就有希望

金虫 (正式写手)

[求助] 非线性回归,超过3个参数的时候,应该如何估计模型参数,求助 已有3人参与

最近在拟合一个非线性模型,自变量只有一个,但是参数比较多,一共有5个,看了一下用nls 进行拟合,都涉及到初值的设定问题,
但是 由于参数较多,不知如何下手,求助各位大侠,是否存在比较稳定的方法,用来解决这种多参数模型,
回复此楼

» 猜你喜欢

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

路就在前方
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

难写啊

新虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
努力就有希望: 金币+5, ★★★很有帮助 2015-07-07 17:59:13
拟合的时候参数就变成了未知数,未知数就是变成了常数(因为被样本替代了)
你首先构造误差函数,这个函数是以参数为未知数的函数,
然后选用一种优化方法,比如随机逼近法,参数初始值可以随机选取,然后基于该初始值求解出参数的收敛值就是结果了
email:myronsaga1@sohu.com.qq:89260998
4楼2015-07-07 15:46:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
努力就有希望: 金币+5, ★★★很有帮助 2015-07-07 17:59:22
feixiaolin: 金币+10, 2015-07应助金币奖励 2015-08-02 15:33:13
试试1stOpt,不需要初值,应该是当今最好用的拟合工具
5楼2015-07-07 17:13:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
努力就有希望: 金币+10, ★★★★★最佳答案 2015-07-07 18:43:08
Function y=p1 + p2*tanh(p3*x + p4);
p1  = 1.88833808781911
p2  = 1.84560725208222
p3  = 1.59338527918929
p4  = -2.72381041627442
8楼2015-07-07 18:36:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

数据放上来看看。

[ 发自小木虫客户端 ]
2楼2015-07-07 13:20:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

努力就有希望

金虫 (正式写手)

引用回帖:
2楼: Originally posted by feixiaolin at 2015-07-07 13:20:48
数据放上来看看。

数据 不太方便,
路就在前方
3楼2015-07-07 13:36:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

努力就有希望

金虫 (正式写手)

引用回帖:
4楼: Originally posted by 难写啊 at 2015-07-07 15:46:13
拟合的时候参数就变成了未知数,未知数就是变成了常数(因为被样本替代了)
你首先构造误差函数,这个函数是以参数为未知数的函数,
然后选用一种优化方法,比如随机逼近法,参数初始值可以随机选取,然后基于该初 ...

你说的是用数值优化吗,我用R里面的optim 试过,用构造了RSS误差函数,但是一样的需要设定初值,这样才能进行搜索,虽然不要求很精确,但是对那个范围要求比较严格,好像行不通
路就在前方
6楼2015-07-07 17:58:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

努力就有希望

金虫 (正式写手)

引用回帖:
5楼: Originally posted by dingd at 2015-07-07 17:13:15
试试1stOpt,不需要初值,应该是当今最好用的拟合工具

这么好,太感谢了,我试试
路就在前方
7楼2015-07-07 17:58:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

努力就有希望

金虫 (正式写手)

引用回帖:
8楼: Originally posted by feixiaolin at 2015-07-07 18:36:07
Function y=p1 + p2*tanh(p3*x + p4);
p1  = 1.88833808781911
p2  = 1.84560725208222
p3  = 1.59338527918929
p4  = -2.72381041627442

厉害
路就在前方
9楼2015-07-07 19:36:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

粥易翻

新虫 (初入文坛)

引用回帖:
8楼: Originally posted by feixiaolin at 2015-07-07 18:36:07
Function y=p1 + p2*tanh(p3*x + p4);
p1  = 1.88833808781911
p2  = 1.84560725208222
p3  = 1.59338527918929
p4  = -2.72381041627442

没看懂啊,怎么冒出来这个的
10楼2015-07-14 09:17:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 努力就有希望 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-02 03:37 by xzcv0Tehd8lQ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-02 03:18 by xzcv0Tehd8lQ
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-02 02:47 by xzcv0Tehd8lQ
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-02 02:45 by xzcv0Tehd8lQ
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +7 X5bHxOULxFqb 2026-10-30 24/1200 2026-11-02 02:39 by xzcv0Tehd8lQ
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 22:27 by xzcv0Tehd8lQ
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 22:24 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 xvJS1BYDxvBS 2026-10-31 10/500 2026-11-01 21:52 by xzcv0Tehd8lQ
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 8/400 2026-11-01 21:35 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-01 21:26 by xzcv0Tehd8lQ
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-01 21:19 by xzcv0Tehd8lQ
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 20:03 by xzcv0Tehd8lQ
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 16:58 by BCdvBAEpt3zm
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 16:52 by BCdvBAEpt3zm
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 16:48 by BCdvBAEpt3zm
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 X5bHxOULxFqb 2026-10-30 22/1100 2026-11-01 16:35 by BCdvBAEpt3zm
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +9 zuesli 2026-10-27 15/750 2026-11-01 12:59 by Equinoxhua
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 3/150 2026-11-01 12:53 by BCdvBAEpt3zm
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 12:40 by BCdvBAEpt3zm
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 03:59 by YXM1FDCSE8Bn
信息提示
请填处理意见