24小时热门版块排行榜    

查看: 2866  |  回复: 19

人间四月

木虫 (著名写手)

[求助] 非线性回归解不唯一

要拟合的公式见图片。要确定的参数是Ym, ETm1, ETm2, ETm3, ETm4,和入1,入2,入3,入4,共9个。我们有18组实验数据,我在EXCEL里用规划求解每次给的初始值不同,得到的结果都不同。请教高人,解不唯一的原因,以及如何得到唯一解(如果有可能)。



数据为:

        ETa1        ETa2        ETa3        ETa4        Ya
1        90.9        174        163.2        126.5        11175
2        82.8        168.2        158.4        124.2        11191
3        71.3        150.8        148.8        121.9        10482
4        112.7        171.1        156        121.9        10185
5        112.7        150.8        148.8        119.6        9867
6        110.4        127.6        141.6        117.3        10168.5
7        115        182.7        139.2        124.2        11622
8        115        188.5        124.8        117.3        9732
9        112.7        176.9        110.4        115        9672
10        110.4        188.5        163.2        101.2        13249.5
11        110.4        188.5        163.2        92        12571.5
12        110.4        185.6        163.2        73.6        11494.5
13        80.5        124.7        134.4        108.1        8182.5
14        112.7        147.9        105.6        103.5        8869.5
15        115        188.5        122.4        69        10041
16        89.7        133.4        88.8        92        7176
17        112.7        171.1        115.2        59.8        8206.5
18        112.7        185.9        160.8        112.7        12910.5

我用EXCEL得到的一组解为:
Ym   ETm1        ETm2        ETm3        ETm4        入1        入2        入3        入4
13116.6         120.6         189.1         173.0         127.7         0.1074         0.4906         0.6184         0.0986
换另外一个初始值又可以收敛到另外一组解。

试过1STOPT,得到的解数量完全不对,不能用,也不知道怎么样设约束条件。

[ Last edited by 人间四月 on 2012-2-28 at 13:26 ]
回复此楼
有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

小雨萌萌(金币+1): 3Q 2012-02-29 09:41:57
人间四月(金币+20): ★★★很有帮助 有办法解决这个问题吗 2012-02-29 10:31:46
给你个简单示例:
拟合公式: y = a +b*c*x
a, b, c为待求参数,这种情况下b和c值就不是唯一的,你的情况就和这类似。
7楼2012-02-29 08:27:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
soliton923(金币+1): 谢谢参与讨论~~~ 2012-02-28 22:16:49
公式有问题,存在过拟合现象,解不唯一。如果有其它约束条件一并给出。
2楼2012-02-28 15:54:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiuyouxu

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
soliton923(金币+1): 谢谢参与讨论~~~ 2012-02-28 22:16:56
等式两边取对数后,即可化为一个多元线性回归问题
忘记自己,忘记一切烦恼(欢迎访问我的网站兆字节:http://www.mathbeta.com/)
3楼2012-02-28 15:58:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

人间四月

木虫 (著名写手)

引用回帖:
3楼: Originally posted by xiuyouxu at 2012-02-28 15:58:47:
等式两边取对数后,即可化为一个多元线性回归问题

ETm是分母啊,MS不能化成线性回归啊
有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有
4楼2012-02-28 18:14:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiuyouxu

铁杆木虫 (职业作家)

【答案】应助回帖

小雨萌萌(金币+1): 3Q 2012-02-29 09:41:48
人间四月(金币+10): 有帮助 2012-02-29 10:33:26
两边取对数后,采用最小二乘法求解
忘记自己,忘记一切烦恼(欢迎访问我的网站兆字节:http://www.mathbeta.com/)
5楼2012-02-28 18:24:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

人间四月

木虫 (著名写手)

引用回帖:
: Originally posted by dingd at 2012-02-28 15:54:27:
公式有问题,存在过拟合现象,解不唯一。如果有其它约束条件一并给出。

谢谢回复,过拟合意思就是需要确定的参数多,而我们的数据不够?
有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有
6楼2012-02-28 21:09:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

人间四月

木虫 (著名写手)

引用回帖:
7楼: Originally posted by dingd at 2012-02-29 08:27:06:
给你个简单示例:
拟合公式: y = a +b*c*x
a, b, c为待求参数,这种情况下b和c值就不是唯一的,你的情况就和这类似。

多谢了,是这个样子,有办法解决这个问题吗?
有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有
8楼2012-02-29 10:32:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

人间四月

木虫 (著名写手)

只要能解决问题,追加悬赏,直到散尽金币。
有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有木有
9楼2012-02-29 10:45:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaohsu2006

木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
人间四月(金币+10): ★★★很有帮助 2012-02-29 13:48:54
多元回归肯定不一,1stopt比较可信吧,它的初值是自己给定
10楼2012-02-29 10:58:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 人间四月 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天系统多次维护,明天很可能放榜! +5 zju2000 2026-08-16 5/250 2026-08-16 23:23 by Haru815
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +32 医学老男孩 2026-08-13 74/3700 2026-08-16 21:05 by 医学老男孩
[论文投稿] 岩土工程学报什么时候才能终审完呐 25+3 yeager111 2026-08-10 4/200 2026-08-16 20:42 by tfang
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +9 Ldrop2023 2026-08-13 9/450 2026-08-16 13:24 by wukun_007
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
信息提示
请填处理意见