24小时热门版块排行榜    

查看: 1088  |  回复: 8

wodaifei

银虫 (小有名气)

[求助] 非线性回归参数..急求帮助!!谢谢

各位学哥学姐:小弟刚开始学matlab(菜鸟级别),遇到问题求帮助:
公式:lny=a+b/x+clnx                求解a,b,c
x        y       
283.15        5.37997
293.15        4.79126
303.15        4.23979        ,,,,,,,
313.15        3.80361       
323.15        3.42896       
333.15        3.10245       
343.15        2.81782       
353.15        0.03688       
按非线性方式回归最后运行出来出现:
Matrix dimensions must agree.
对于这个问题要想得到参数该怎么解决啊,
小弟很着急,各位学哥学姐帮帮忙啊
万分感谢!!!!
回复此楼

» 猜你喜欢

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

要想拥有一切,就要让自己变得足够优秀!!!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiuyouxu

铁杆木虫 (职业作家)

【答案】应助回帖


soliton923(金币+1): 谢谢参与~~ 2011-06-29 20:17:53
wodaifei(金币+1): 2011-06-30 10:14:53
Matrix dimensions must agree.

矩阵维数不一致, 检查矩阵运算是否有问题.
忘记自己,忘记一切烦恼(欢迎访问我的网站兆字节:http://www.mathbeta.com/)
2楼2011-06-29 16:44:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖


soliton923(金币+1): 谢谢参与~~ 2011-06-29 20:18:05
wodaifei(金币+1): 2011-06-30 08:39:53
你把方程“lny=a+b/x+clnx”改为:

y=exp(a+b/x+c*ln(x))

试试!
3楼2011-06-29 20:08:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pengyehui

木虫 (正式写手)

【答案】应助回帖

wodaifei(金币+2): 2011-06-30 10:14:38
检查数据的输入吧!
三楼的建议值得商榷,因为原来问题已经可以线性化了,按三楼的方法,把问题搞成非线性的,复杂了!商榷而已!
4楼2011-06-29 22:07:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wodaifei

银虫 (小有名气)

引用回帖:
Originally posted by pengyehui at 2011-06-29 22:07:59:
检查数据的输入吧!
三楼的建议值得商榷,因为原来问题已经可以线性化了,按三楼的方法,把问题搞成非线性的,复杂了!商榷而已!

这应该不是线性吧,否则应该能算出来,我就按照正确步骤一步一步来的,最后出现以上问题,我也知道是矩阵的问题,可小弟就是不知道该怎么改。。。
要想拥有一切,就要让自己变得足够优秀!!!
5楼2011-06-30 08:43:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pengyehui

木虫 (正式写手)

【答案】应助回帖


lovibond(金币+1): 鼓励交流 2011-07-01 08:32:14
令Y=lny; X1=1/x; X2=lnx
不就线性了吗
6楼2011-06-30 09:38:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wodaifei

银虫 (小有名气)

引用回帖:
Originally posted by pengyehui at 2011-06-30 09:38:33:
令Y=lny; X1=1/x; X2=lnx
不就线性了吗

非常感谢啊,呵呵,小弟愚昧。。。
要想拥有一切,就要让自己变得足够优秀!!!
7楼2011-06-30 10:14:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

wodaifei(金币+1): 2011-06-30 19:20:44
如果原问题是非线性的,不建议线性化后求解,因为误差有可能非常大,以此题为例,线性化后公式及数据:
y1=a+b*x1+c*x2;
数据:
x1        x2        y1
0.00353169698039908        5.64597679255981        1.68268279794997
0.00341122292341805        5.68068442341022        1.56679342484846
0.0032986970146792        5.71422773251774        1.44451373971541
0.00319335781574325        5.74668230897142        1.33595061576793
0.00309453814018258        5.7781166117089        1.23225700820423
0.0030016509079994        5.80859283900855        1.13219212193131
0.00291417747340813        5.83816766935464        1.03596353634767
0.00283165793572137        5.86689289585493        -3.30008588029374

拟合结果:
a        2392.53187520471
b        -107971.570616336
c        -356.014799219908

直接非线性拟合:
公式: y=exp(a+b/x+c*ln(x))
数据:
x        y
283.15        5.37997
293.15        4.79126
303.15        4.23979
313.15        3.80361
323.15        3.42896
333.15        3.10245
343.15        2.81782
353.15        0.03688

结果:
a        392.19036548745
b        -16708.6338358188
c        -58.7228208089647

非线性直接误差平方和:3.26128834414432
线性化后求得参数代入原方程误差平方和: 30.7310566203494

误差是惊人的。
8楼2011-06-30 15:48:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wodaifei

银虫 (小有名气)

引用回帖:
Originally posted by dingd at 2011-06-30 15:48:50:
如果原问题是非线性的,不建议线性化后求解,因为误差有可能非常大,以此题为例,线性化后公式及数据:
y1=a+b*x1+c*x2;
数据:
x1        x2        y1
0.00353169698039908        5.64597679255981        1.68268279794997
0.0034112 ...

哥,真的很感谢啊。。。你说的对,我做的时候也感觉有问题。。同学给的数据确实存在问题,不过方法是对的。。。再次感谢!!!
要想拥有一切,就要让自己变得足够优秀!!!
9楼2011-06-30 19:22:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wodaifei 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] 跑板能跑开,过柱过不纯怎么办 +7 小胡在努力 2024-06-18 9/450 2024-06-20 23:04 by if今年的吧
[论文投稿] Pattern Recognition期刊,二审审稿邀请发出一个月了,但有一个审稿人一直未接受审稿 5+3 PLVS_VLTRA 2024-06-19 5/250 2024-06-20 22:54 by holypower
[基金申请] 我标书代码变了 +63 学员NHuqdk 2024-06-16 86/4300 2024-06-20 22:03 by 841790061
[找工作] 药学硕士找不到工作,打算去做科研助理了 +11 pom戴墨镜 2024-06-14 24/1200 2024-06-20 21:59 by 海阔天空4924
[基金申请] 工材口青年基金大概什么样能上会? +12 今晚推荐22 2024-06-20 17/850 2024-06-20 21:50 by 今晚推荐22
[硕博家园] 关于硕博连读的一些疑问? +9 Lwenter 2024-06-14 11/550 2024-06-20 21:20 by rhhua11111
[硕博家园] 怎么带研究生? +10 豆豆小小2 2024-06-20 10/500 2024-06-20 20:50 by rhhua11111
[基金申请] 江南大学到瑞士招聘,称取消非升即走,改预聘+长聘 +21 babu2015 2024-06-18 23/1150 2024-06-20 18:25 by zpv2011
[论文投稿] 投稿求助 6+3 小短腿安琪拉 2024-06-19 4/200 2024-06-20 17:40 by 不一样烟火12345
[有机交流] 怎么萃取出锡盐内包裹的化合物 +3 硕六过 2024-06-19 4/200 2024-06-20 15:45 by 2461777831
[找工作] 高校两个offer选择 +15 cowox2021 2024-06-18 16/800 2024-06-20 15:24 by ase123456
[考博] 这个博士要读吗 +15 Sea Breeze 2024-06-16 29/1450 2024-06-20 14:30 by Sea Breeze
[硕博家园] 豫北虫友互识 +11 xuhongli903 2024-06-18 12/600 2024-06-20 13:32 by 1429524838
[基金申请] F03青年基金函评结果 +3 暨阳一只柴 2024-06-19 4/200 2024-06-20 11:15 by 暨阳一只柴
[考博] 关于读博感觉自己很抓马 +9 小九月 2024-06-19 12/600 2024-06-20 11:06 by SCITOPPP
[考博] 2025年博士申请,可先做科研助理 +4 limit888 2024-06-18 7/350 2024-06-20 10:23 by 安塔瓦拉多
[论文投稿] Scientific Reports二审90天了,催稿也没反应怎么办? +3 智晨君 2024-06-18 6/300 2024-06-19 20:08 by licheng78
[基金申请] F口401需要啥文章水平 +3 lhjr123 2024-06-16 7/350 2024-06-18 16:05 by hon920603
[硕博家园] 博士毕业高校和就业的相关问题 +7 SCITOPPP 2024-06-14 11/550 2024-06-18 07:51 by yinxing1995
[基金申请] 关于博后基金的bug问题 +6 lxr1991 2024-06-14 9/450 2024-06-15 21:17 by since—2010
信息提示
请填处理意见